采用了LS-DYNA中13號材料*MAT_ISOTROPIC_ELASTIC_FAILURE簡單模擬了一下，冰塊撞擊碎裂的過程，相較于生死單元cohesive方法，計算更加方便、便捷，效果不差。模型如下：

model1.png

model2.png

材料參數(shù)供參考：（cm-g-us）

*MAT_ISOTROPIC_ELASTIC_FAILURE_TITLE

ICE

$#     mid        ro         g      sigy      etan      bulk    

         1       0.9     0.0222.12000E-5    0.0426    0.0526

$#     epf       prf       rem      trem    

      0.35     -4E-5       0.0       0.0

冰塊運動由以下關(guān)鍵字實現(xiàn)

*BOUNDARY_PRESCRIBED_MOTION_SET_ID

$#      id                                                               heading

         1translation

$#    nsid       dof       vad      lcid        sf       vid     death     birth

         2         1         0         1      -1.0         0   9000000       0.0

*DEFINE_CURVE_TITLE

zero

$#    lcid      sidr       sfa       sfo      offa      offo    dattyp     lcint

         1         0       1.0       1.0       0.0       0.0         0         0

$#                a1                  o1  

                 0.0    9.9999997474e-05

    9.9999997782e+21    9.9999997474e-05