SimWe 仿真论坛---(论坛实行邀请码注册)'s Archiver SimWe 仿真论坛---(论坛实行邀请码注册) ? A04：ABAQUS ? ABAQUS 子程序VUMAT 的坐标系的问题的讨论billowriver 发表于 2009-4-6 05:50

ABAQUS 子程序VUMAT 的坐标系的问题的讨论

在ABAQUS/Explicit 里面，Vumat 里的所有变量都是相对于一个corotational 坐标系 来说的，而在ABAQUS 的odb 里面显示的是cauchy stress, 是相对于当前坐标系来讲的。

因为变形梯度 F=RU ，意思是材料先变形U ，再做一个刚体旋转R ，而这个corotational 坐标系是在材料做刚体旋转前的坐标系，ABAQUS 上面说由于做刚体旋转不会影响应力在相应坐标系里的大小，因此刚体旋转前后应力大小应该一样，即VUMAT 里面的StressNew 和odb 里面显示的cauchy stress 大小应该一样。 我在做单轴拉伸模拟得到的结果二者一样，但是如果做剪切模拟，如果剪切力和正应力同时存在的话，StressNew 和odb 里面的应力不一样，即使把StressNew 设成常数，在odb 里得到的应力也不是常数。很奇怪。 这种影响在小变形下差别很小，但是大变形时候会产生很大误差。如何用VUMAT 计算出当前的cauchy stress, 哪位高手知道或者有兴趣的可以一起讨论下。

S_true=RSR'，但是ABAQUS 里面没有直接提供R ，而是提供了F_old ，F_new, U_old, U_new ，但是Vumat 里面应变也是基于corotational 坐标系的，而且是lnV,在这里=lnU

shawn2008 发表于 2009-4-6 06:05

[i=s] 本帖最后由 shawn2008 于 2009-4-6 06:07 编辑 [/i]

太多概念错误，我都不知道从哪开始说起了。

billowriver 发表于 2009-4-9 01:02

请多指教，这个都把我搞糊涂了

billowriver 发表于 2009-4-9 01:09

上面我讲的都是出自ABAQUS document ，而且从连续介质力学角度讲也应该是这样，只是ABAQUS 处理VUMAT 时多了个corotational 坐标系，得到的结果有些问题

不知道哪写概念有问题，还请shawn2008 说说

shawn2008 发表于 2009-4-9 01:40

呵呵，别误会，我以前也是像你写的这样理解的；document 的东西是有前提的。 ok, 我们先说第一句的前半句话，’在ABAQUS/Explicit 里面，Vumat 里的所有变量都是相对于一个corotational 坐标系 来说的‘。

这里有几个概念请你要先说明一下， ‘所有变量’指的是什么变量？ corotational 坐标系是怎么定义的？ 相对于是指的什么含义？

billowriver 发表于 2009-4-9 02:29

[i=s] 本帖最后由 billowriver 于 2009-4-9 02:43 编辑 [/i] 谢谢shawn2008的答复，所有变量就是应力，应变，应变率这些，当然ABAQUS 给出的变形梯度F 应该还是度量current 坐标系和

initial 坐标系的。 有限元分析核心增值经销商

北京盛维提供3D CAD/CAE/CAM/PDM/CA... 整

体解决方案。010-********/68253809