python turtle绘图中角度坐标系的绝对零度方向是_ABAQUS 子程序VUMAT 的坐标系的问题的讨论...-CSDN博客

本文链接：https://blog.csdn.net/weixin_29801567/article/details/112957960

在ABAQUS/Explicit中，VUMAT子程序的变量基于corotational坐标系，而ODB显示的Cauchy应力对应当前坐标系。在大变形情况下，VUMAT计算的StressNew与ODB中的应力可能存在差异，特别是在剪切力与正应力共存的场景。讨论涉及到如何利用VUMAT计算当前Cauchy应力，以及文档和实际结果之间的差异。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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 坐标系是怎么定义的？相对于是指的什么含义？