libaochao 发表于 2012-6-18 22:19:21

一个损伤演化的UMAT程序例子

小弟最近正在研究istropic损伤的一个模型,根据manual中的例子,自己编了一个针对自己模型的一个UMAT程序, 传上来跟大家分享下.如果程序中有什么需要改进的地方,还请各位大侠多多请教.

再附两张结果图吧,第一张是应力-应变曲线,第二张是损伤参数-应变曲线


hanmumu 发表于 2012-6-19 16:55:50

本帖最后由 hanmumu 于 2012-6-19 16:56 编辑

同学，你看的手册哪里啊，我怎么一直找不到啊，不会查手册，感觉找到的都不是自己需要的，请指导下吧

libaochao 发表于 2012-6-19 18:37:40

hanmumu 发表于 2012-6-19 16:55 static/image/common/back.gif
同学，你看的手册哪里啊，我怎么一直找不到啊，不会查手册，感觉找到的都不是自己需要的，请指导下吧 ...

就是ABAQUS的documentation， 然后搜索UMAT damage，然后在EXAMPLE里面找到的

vigorao 发表于 2012-6-19 19:29:42

你的变量
   XT=PROPS(3)
      XC=PROPS(4)
      XS=PROPS(5)
      GX=PROPS(6)   !Fracture energy in matrix
      ETA=0.001
代表什么意思啊，能说明下吗

libaochao 发表于 2012-6-19 20:00:43

vigorao 发表于 2012-6-19 19:29 static/image/common/back.gif
你的变量
   XT=PROPS(3)
      XC=PROPS(4)


XT,XC,XS分别是拉伸,压缩和剪切的强度,GX是损伤演化模型中的一个能量常数,而ETA则是为了提高收敛性而设置的一个粘性常数, ETA 并不是必需的

hanmumu 发表于 2012-6-20 08:35:12

libaochao 发表于 2012-6-19 18:37 static/image/common/back.gif
就是ABAQUS的documentation， 然后搜索UMAT damage，然后在EXAMPLE里面找到的

谢谢啦，可以留下qq么，小弟不才，以后向您老请教下，嘎嘎

vigorao 发表于 2012-6-20 12:39:06

能不能把你的，程序变量说明下啊，有点看不懂，谢谢

libaochao 发表于 2012-6-20 15:59:33

本帖最后由 libaochao 于 2012-6-20 16:03 编辑

vigorao 发表于 2012-6-20 12:39 static/image/common/back.gif
能不能把你的，程序变量说明下啊，有点看不懂，谢谢
矩阵C（i,j）为初始刚度矩阵，CD(i,j)为有损伤时的有效刚度矩阵，STRANT是步长末应变，DFMNDE之类的变量就是应变准则变量FMN对应变的导数，d(FMN)/d(STRANT)。DDMDE就是损伤变量DM对应变的导数d(DM)/d(STRANT)，之后的都类似，这些变量的目的就是求出最后UMAT要求的DDSDDE矩阵。DDSDDE的定义是根据手册里给出的公式来算的，程序里的变量对应相应的一些微分项http://pc-35630.interne.emd:2080/v6.11/books/exa/graphics/exa_eqn00373.gif，公式中有分成基体跟纤维两部分的，我的程序里就只有基体部分，没有纤维部分

libaochao 发表于 2012-6-20 16:01:34

hanmumu 发表于 2012-6-20 08:35 static/image/common/back.gif
谢谢啦，可以留下qq么，小弟不才，以后向您老请教下，嘎嘎

373469352，其实我研究也不多，只是看到了手册上的一个例子，拿过来模仿下，有什么经验可以多交流

ljhzyq 发表于 2012-7-12 20:34:46

lz的模型很有意思 我曾经尝试做损伤的vumat模型 后来被各种问题困扰 就没有继续下去
还没有仔细看你的代码
不过想问个初级的问题
你的这个损伤模型里好像没有mises的塑性变形 如果我想加入塑性变形和等效强化
然后比如说当应变达到一个特定值发生损伤
你觉得这个能够成立吗

刘晗 发表于 2012-7-16 19:06:13

libaochao 发表于 2012-6-20 16:01 static/image/common/back.gif
373469352，其实我研究也不多，只是看到了手册上的一个例子，拿过来模仿下，有什么经验可以多交流 ...

朋友，做过GTN的umat吗？

libaochao 发表于 2012-8-3 22:03:16

ljhzyq 发表于 2012-7-12 20:34 static/image/common/back.gif
lz的模型很有意思 我曾经尝试做损伤的vumat模型 后来被各种问题困扰 就没有继续下去
还没有仔细看你的代码
...

前些天放假了，刚回来。我没怎么研究过VUMAT，不过我觉得可以加入塑性，只要选择合适的塑性应该就可以了吧。我写的这个UMAT其实也是用的应变，不过是根据应力转化了一下，因为好多材料给出的都是应力极限

ljhzyq 发表于 2012-8-22 17:10:00

libaochao 发表于 2012-8-3 22:03 static/image/common/back.gif
前些天放假了，刚回来。我没怎么研究过VUMAT，不过我觉得可以加入塑性，只要选择合适的塑性应该就可以了 ...

多谢回复！
只是加了塑性和强化以后，Cep的求解就更复杂一些了，不想只有弹性的情况...

amanada 发表于 2016-11-20 14:20:07

感谢楼主分享，自己来研究研究，自己编的程序有点点问题，看看借鉴下能否解决

newwalker 发表于 2017-5-10 11:48:24

请问一下楼主的损伤变量D是基于什么理论？

romancezyl 发表于 2017-5-27 10:58:55

谢谢楼主谢谢楼主

linxinhao 发表于 2020-11-3 16:00:36

感谢楼主
页: [1]
查看完整版本: 一个损伤演化的UMAT程序例子