求助！把白金泽老师的例题材料换成空气K文件出错！

ansys平安

各位好，感谢你浏览我的帖子！
我把白金泽老师的书上第五章的共节点模拟爆炸例题中的被爆材料换成了空气，不知为什么运行出错！由于是新手，请各位批评指正！

以下是部分k文件：
$                             MATERIAL DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$
$
*MAT_NULL
         1 1.2929-03
*EOS_LINEAR_POLYNOMIAL
         1 -1.00E-06       0.0       0.0       0.0       0.4       0.4       0.0
2.500E-06       1.0
*MAT_HIGH_EXPLOSIVE_BURN
$      MID        RO         D       PCJ      BETA         K         G      SIGY
         2 1.6310E+3 0.6718E+4 0.185E+11 0.0000000 0.0000000
0.0000000 0.0000000
*EOS_JWL
$    EOSID         A         B        R1        R2      OMEG        E0        V0
         2 5.409E+11 0.094E+11 4.5000000 1.1000000 0.3500000
0.080E+11 0.0000000
$         
*INITIAL_DETONATION
$      PID         X         Y         Z        LT
         2       0.0       0.0       0.0       0.0
$
===================================================================
$                              PARTS DEFINITIONS                               $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$
$
$
*PART
Part          1 for Mat         1 and Elem Type         1
         1         1         1         1         0         0         0
$
*PART
Part          2 for Mat         2 and Elem Type         2
         2         2         2         2         0         0         0
==================================================================
$
*CONTROL_ENERGY
         2         2         2         2
*CONTROL_SHELL
  20.0             1        -1         1         2         2         1
*CONTROL_TIMESTEP
    0.0000    0.9000         0  0.00      0.00   
*CONTROL_TERMINATION
0.800E-03         0   0.00000   0.00000   0.00000
*DATABASE_HISTORY_NODE
         1
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$
$                                 TIME HISTORY                                 $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$
$
*DATABASE_MATSUM  
0.8000E-05
*DATABASE_GLSTAT  
0.8000E-05
*DATABASE_BINARY_D3PLOT
0.8000E-05
*DATABASE_BINARY_D3THDT
0.8000E-05


以下是message信息：
1 t 0.0000E+00 dt 4.82E-06 flush i/o buffers
       1 t 0.0000E+00 dt 4.82E-06 write d3plot file
       2 t 4.8183E-06 dt 4.82E-06 write d3plot file
       4 t 1.3981E-05 dt 4.78E-06 write d3plot file
       6 t 2.3717E-05 dt 1.76E-06 write d3plot file
*** Error negative volume in brick element #       6 cycle       7
*** Error negative volume in brick element #      29 cycle       7
*** Error negative volume in brick element #      63 cycle       7
*** Error negative volume in brick element #      75 cycle       7
*** Error negative volume in brick element #      97 cycle       7
*** Error negative volume in brick element #      98 cycle       7
*** Error negative volume in brick element #     109 cycle       7
*** Error negative volume in brick element #     394 cycle       7
*** Error negative volume in brick element #     395 cycle       7
*** Error negative volume in brick element #     433 cycle       7
*** Error negative volume in brick element #     435 cycle       7
*** Error negative volume in brick element #     457 cycle       7
*** Error negative volume in brick element #     721 cycle       7
*** Error negative volume in brick element #     722 cycle       7
*** Error negative volume in brick element #     745 cycle       7
******** minimum time step size reached *********
On a restart minimum timestep will be  0.000E+00
       7 t 2.6265E-05 dt-2.78E-01 write d3plot file

highbrow

"-1.00E-06"是对的吗?我不确定

ansys平安

状态方程没有错误

zpyc01

出现负体积了
*CONTROL_TIMESTEP中0.9改为0.67试试

ansys平安

您好！按您的提示问题已经解决，真的很感谢您！
由于我是初学者，还从没遇见过这种问题，在书上也没找到这方面资料，您可以简单的介绍一下负体积问题吗？谢谢！

zpyc01

5# ansys平安

负体积：
对于承受很大变形的材料，比如说泡沫，一个单元可能变得非常扭曲以至于单元的体积计算得到一个负
值。这可能发生在材料还没有达到失效标准前。对一个拉格朗日(Lagrangian)网格在没有采取网格光滑(mesh
smoothing)或者重划分(remeshing)时能适应多大变形有个内在的限制。LS-DYNA 中计算得到负体积(negative
volume)会导致计算终止，除非在*control_timestep 卡里面设置ERODE 选项为1，而且在*control_termination
里设置DTMIN 项为任何非零的值，在这种情况下，出现负体积的单元会被删掉而且计算继续进行(大多数情
况)。有时即使ERODE 和DTMIN 换上面说的设置了，负体积可能还是会导致因错误终止。

解决方法：
* 简单的把材料应力－应变曲线在大应变时硬化。这种方法会非常有效。
* 有时候修改初始网格来适应特定的变形场将阻止负体积的形成。此外，负体积通常只对非常严重的变
形情况是个问题，而且特别是仅发生在像泡沫这样的软的材料上面。
* 减小时间步缩放系数(timestep scale factor)。缺省的0.9 可能不足以防止数值不稳定。
* 避免用全积分的体单元(单元类型2 和3)，它们在包含大变形和扭曲的仿真中往往不是很稳定。全积分
单元在大变形的时候鲁棒性不如单点积分单元，因为单元的一个积分点可能出现负的Jacobian 而整个单元还
维持正的体积。在计算中用全积分单元因计算出现负的Jacobian 而终止会比单元积分单元来得快。
* 用缺省的单元方程(单点积分体单元)和类型4 或者5 的沙漏(hourglass)控制(将会刚化响应)。对泡沫材
料首先的沙漏方程是：如果低速冲击type 6，系数1.0; 高速冲击type 2 或者3。
* 对泡沫用四面体(tetrahedral)单元来建模，使用类型10 体单元。
* 增加DAMP 参数(foam model 57)到最大的推荐值0.5。
* 对包含泡沫的接触，用*contact 选项卡B 来关掉shooting node logic。
* 使用*contact_interior 卡
用part set 来定义需要用contact_interior 来处理的parts，在set_part 卡1 的第5 项DA4来定义contact_interior
类型。缺省类型是1，推荐用于单一的压缩。在版本970 里，类型1 的体单元可以设置type=2，这样可以处
理压缩和减切混合的模式。
* 如果用mat_126，尝试ELFORM=0
* 尝试用EFG 方程(*section_solid_EFG)。因为这个方程非常费时，所以只用在变形严重的地方，而且只
用于六面体单元。

另外关于负体积等计算不稳定的讨论，参见大牛的贴http://forum.simwe.com/thread-960209-1-1.html

ansys平安

由衷的谢谢您！在像您这样的前辈的不吝帮助下！我相信我们初学者会进步很快！

zpyc01

过奖了，称不上前辈，论坛里的大牛们才是真正的前辈。互相学习，才能共同进步
7# ansys平安