清华大学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，怀疑有误

: 还望大侠提醒！

--

你只有一个单元么？

【在saturn (好好工作) 的大作中提到: 】

: ~!~~~~~~~~~~~~~~~~~~~~~~~~~~~~这我以前还不知道

: 应该不是0，如果在每个增量步的最后你重新更新过他

: 输入文件*.inp中，*depvar不要忘了写statev()的维数哦！

--

【在air1 (无影如风※勤奋羊羊※课题在即) 的大作中提到: 】

: 可是这个写的命令有用么

: 对于输出而言

: 你只有一个单元么？

已经写了，定义了statev的维数，

*depvar