本帖最后由 feiyafei007 于 2014-6-11 09:56 编辑
使用VUMAT子程序计算SHPB撞击时出现单元变形过大，不承受应力的情况。 1、 SHPB实验得到的应力波（图1）：
图1 2、 使用UMAT对本构模型编制的程序进行计算得到的应力波（图2）：
图2 3、 使用VUMAT对本构模型编制的程序进行计算得到的应力波（图3）：
图3 4、对比两种程序所得到的实验结果可知，试件在UMAT计算过程中，虽然试件发生了较大的变形（图4），但单元变形均匀且可承受应力，得到的结果如图2所示：
图4 试件使用VUMAT计算的过程中，在试件的小变形阶段，单元可承受并传递应力（图5）；大变形阶段，单元不再承受应力且不能传递应力（图6）。 图5 图6 即使用UMAT编制的程序计算时所得到的结果与实验结果较为符合，但VUMAT计算结果与实验结果不同。引起上述现象出现的原因是在VUMAT计算过程中，单元不规则的大变形使其失去了承受应力的能力，从而不能传递应力波。本人试过增大试件（中间直径小的part为试件，试件两端的part为铝杆）网格密度，对试件的外表面进行绑定等方法，但单元在大变形阶段仍旧失效，找不到失效原因。请各位大神莅临指导~~ 补充：本构模型为ZWT模型。
| 
IP卡
狗仔卡
发表于 2014-6-11 09:50:23
提升卡
沉默卡
喧嚣卡
变色卡
抢沙发
千斤顶
显身卡
