15.利用命令进行计算时如何设置调用内存量
abaqus job=job-name memory=memory-size
16.清华大学bbs的abaqus精华
【 在 luxz (panda) 的大作中提到: 】
: 哪位有什么好办法？
: 感觉abaqus自己的前处理用得很不爽
: patran可以么？
我感觉abaqus/cae不错的呀,实体建模颇有pro/e的风格
感觉patran也好不到哪去,太麻烦,
也许patran功能更强一些(要不然也对不起它那硕大的身躯呀),
但对于不是很庞大的模型,根本体现不出来什么优势
: air1大侠，本人给你作广告，为何不可？
: 请问：弹塑性矩阵【D]与ddsdde有何联系，
: 你用过板壳单元吗？
stress=D*stran?
d(stress)=ddsdde*d(stran)
--
那应该就是一样的，因为全量理论，
Sij=DijklEkl（满足张量求和约定）
即Stress=D*Strain;
而在增量理论中，
△S=D*△E(在有限变形中，△其实应该为应力的客观率）
--
似乎不对吧
大变形下此D非彼D
你看过黄克智的固体本构关系这本书么
如果你从全量理论和增量理论的角度上讲
那似乎第一个Digkl就不对
你有第一个式子么
如果有，求导不久完了？
: 那应该就是一样的，因为全量理论，
: Sij=DijklEkl（满足张量求和约定）
: 即Stress=D*Strain;
: 而在增量理论中，
: △S=D*△E(在有限变形中，△其实应该为应力的客观率）
是啊，大变形下的[D]与普通意义下的[D]在构型上是不一样的，
毕竟[D]大是变形历史的函数，而[D]小则不是，
我推导一种新的本构关系，
△Sij=Dijkl△Ekl (其中△为Jaumann率）
假设材料一开始就屈服（即屈服面为0）
想用壳单元，
: 似乎不对吧
: 大变形下此D非彼D
: 你看过黄克智的固体本构关系这本书么
: 如果你从全量理论和增量理论的角度上讲
: 那似乎第一个Digkl就不对
: 你有第一个式子么
: 如果有，求导不久完了？
: 是啊，大变形下的[D]与普通意义下的[D]在构型上是不一样的，
: 毕竟[D]大是变形历史的函数，而[D]小则不是，
: 我推导一种新的本构关系，
: △Sij=Dijkl△Ekl (其中△为Jaumann率）
: 假设材料一开始就屈服（即屈服面为0）
: 想用壳单元，
唉，别提了，问题就出在，在abaqus中，明明写着可以考虑剪切效应，
可我打印出剪切力个数是，nshr=1,即只有S12,
那我的S13,S23就不知怎么计算，
（DDSDDE(5,5)无法计算，因为ntens=3,最多只能计算DDSDDE(3,3))
你编umat编进去不久行了
他让用NDI,NSHR,NTENS表示变量，你就用这些表示变量
这样他就可以任意的计算了呀，而不在意实际计算的变量数
: 唉，别提了，问题就出在，在abaqus中，明明写着可以考虑剪切效应，
: 可我打印出剪切力个数是，nshr=1,即只有S12,
: 那我的S13,S23就不知怎么计算，
: （DDSDDE(5,5)无法计算，因为ntens=3,最多只能计算DDSDDE(3,3))
因为我的UMAT从abaqus传来的变量（当我选用壳单元时）ntens=3,
ndi=2,nshr=1,这样的话，STRESS为3个，STRESS(3),不可能计算
STRESS(4),STRESS(5),同理，STRAIN,DDSDDE也存在同样的情况
在中厚度板壳元中(MINDLIN)中，DDSDDE为5×5
你的FOR文件收到，可惜研究的对象相差太大，看不懂！
不过你在文件中定义的变量有的根本就没有用到，
还有有的地方似乎是毫无用处的代码，比如：SMISES一段
我觉得你不应该太急，第一步应该保证代码的正确性。
另外看了你上面的讨论，我觉得你应该换一种单元试一试。
比如用梁单元。很有可能NDI和NSHR会不一样。
不知道你是不是一定要用到壳单元上？
如果必须用在壳单元上，可能S23，S31就必须自己在UAMT
中输出了。具体的做法我也不清楚。
还有对DDSDDE，我的理解也是: dS=DDSDDE*dE
我的意思是不管他是几乘几的矩阵
一般的D矩阵就是6*6
你编的umat需要适应于任何单元
所以当然6*6的矩阵应该给完全
当算体的时候
如果是壳单元，他只有ntens=3
那就等于三呗，他用哪个就是哪个
反正你的D6*6的矩阵放在哪了
如果是三维实体单元，ntens=6
所以你的D矩阵就全用上啦
记住你编的umat应该适应于所有的单元，所以不用考虑哪几个参数的值
而且这个stress(3)也许是人家只用到了你D6*6的第一行，第二行，和的四行哪
sigh
不知道我说明白了没有
: 因为我的UMAT从abaqus传来的变量（当我选用壳单元时）ntens=3,
: ndi=2,nshr=1,这样的话，STRESS为3个，STRESS(3),不可能计算
: STRESS(4),STRESS(5),同理，STRAIN,DDSDDE也存在同样的情况
: 在中厚度板壳元中(MINDLIN)中，DDSDDE为5×5
对于DDSDDE，其实就是给了一个应力应变的斜率，我感觉
真正计算的是利用更新的应力，DDSDDE我觉得就是为了确保收敛
真开心有这么多人讨论ABAQUS，特别是umat
我国外的师兄说Umaterial在国外应用得也很广
: 对于DDSDDE，其实就是给了一个应力应变的斜率，我感觉
: 真正计算的是利用更新的应力，DDSDDE我觉得就是为了确保收敛
: 真开心有这么多人讨论ABAQUS，特别是umat
: 我国外的师兄说Umaterial在国外应用得也很广

各位大侠，看到我提出的uamt,及statev有这么多人感兴趣，实在很
高兴，希望大家能把摸索到的东西介绍给大家，现在我的程序能通了，
（因没有屈服面的概念),虽然结果是错误的，但总可以静下心来，不像以前
一运行自己的umat程序就standard.exe error.(并且程序.dat.msg.中无
任何错误信息，log说error see msg,but msg中无）。
现在我明白了，在用户程序中，WRITE(*,*)则写入log,write(6,*)则写入
dat,等，
在监控程序时，发现，statev()，是每次迭代和增量运行它，他就为0，
还是为上一步的值，因为我每次都为0，怀疑有误
还望大侠提醒！
: dat,等，
: 在监控程序时，发现，statev()，是每次迭代和增量运行它，他就为0，
应该不是0，如果在每个增量步的最后你重新更新过他
输入文件*.inp中，*depvar不要忘了写statev()的维数哦！
: 还是为上一步的值，因为我每次都为0，怀疑有误

未完，请看最新Abaqus常用技巧大全（四）

相关标签搜索：Abaqus分析实例  Abaqus培训  静态应力/位移分析  广州有限元培训  solidworks培训  CAD培训  ansys培训  solidworks  proe培训  运动仿真  有限元FEA  

编辑