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

SimWe仿真论坛(论坛实行邀请码注册)s Archiver 有限元分析核心增值经销商 北京盛维提供3D CAD/CAE/CAM/PDM/CA... 整 体解决方案。01068253809 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坐标系的。 corotational坐标系是材料只做变形，在发生刚体转动前的坐标系，这样我们就有了三个坐标系，initial coordinates, corotational coordinates, current coordinates， “相对于corotational坐标系”, 位移u_corotational=x_corotational-X_initial, 定义应变要用u_corotational 和 x_corotational，定义应力时，微元面积和垂直方向也要用x_corotational，此时的变形梯度其实就是U shawn2008 发表于 2009-4-9 09:51 [i=s] 本帖最后由 shawn2008 于 2009-4-9 09:54 编辑 [/i] 好的，所有的变量可以归纳为材料相关的张量，比如应力，应变，应力率，应变率，变形梯度F，速度梯度L，变形率 张量D等等。 张