负体积问题

wu_lei

以下是我的模型

wu_lei

当我施加一个x方向的位移时，发生如下变化

wu_lei

以下是我的k文件，有兴趣的可以探讨一下，我的qq:286963858，请高手指导一下，谢谢！

jzbai

把模型简单介绍一下吧，做什么用的呢？

changTT

把*mat_rigid con1=5.07(7.07)改成5(7).

wu_lei

里面的两个小圆是刚体，左边的一个全部约束不动，右边一个只剩x方向平动自由度，施加沿x轴的平动位移，就产生了上序现象。
TO changTT:”“把*mat_rigid con1=5.07(7.07)改成5(7)”“不好意思，这个是笔误，谢谢指教！！

wu_lei

这是修改后的结果，计算到最后出现负体积，好像接触有点问题，刚体都跑到材料里面去了，请指教，谢谢！

wu_lei

附上k文件，以下是message提示：
*** Error negative volume in brick element #    1684 cycle 1522850

*** Error negative volume in brick element #    2557 cycle 1522850

******** minimum time step size reached *********
On a restart minimum timestep will be  0.000E+00
1522850 t 8.2924E-02 dt-1.36E-08 write d3plot file

changTT

指教，不敢当。只是对ls-dyna感兴趣，大家多交流学习嘛。

wu_lei

好，多多交流！
请问各位高手这个接触穿透怎么解决呢？

wu_lei

谢谢wenhongyu！
这个问题解决了，我把关键字奉上，希望对大家有用：
*RIGIDWALL_GEOMETRIC_CYLINDER_MOTION
$^CONTACT2
$     NSID    NSIDEX     BOXID
                   0         0
$       XT        YT        ZT        XH        YH        ZH      FRIC
  48.20522 -0.680765     -40.0  48.20522 -0.680765     -50.0       0.0
$   RADCYL    LENCYL
      18.0      60.0
$     LCID       OPT        VX        VY        VZ
         1         1       1.0                    
*RIGIDWALL_GEOMETRIC_CYLINDER
$^CONTACT3
$     NSID    NSIDEX     BOXID
                   0         0
$       XT        YT        ZT        XH        YH        ZH      FRIC
-47.79478 -0.680765     -10.0 -47.79478 -0.680765     -20.0       0.0
$   RADCYL    LENCYL
      18.0      60.0