第九章 LS-DYNA动力学分析

5．命令流文件

/PREP7     ！进入前处理器
ET,1,SOLID164    ！定义单元类型
MP,DENS,1,0.99821   ！定义材料密度
MP,EX,1,0     ！定义弹性模量
MP,NUXY,1,0    ！定义泊松比
MP,DENS,2,1.80
MP,EX,2,
MP,NUXY,2,
MP,DENS,3,0.001252
MP,EX,3,0
MP,NUXY,3,0
MP,DENS,4,2.65
MP,EX,4,40
MP,NUXY,4,0.3
CSYS,4     ！激活前面定义的数据
BLOCK,0,255,-250,140,0,1.5,  ！创建模型体
wpro,,-90.000000,    ！将工作平面绕X轴顺时针旋转90度
wpoff,0,0,-6    ！将工作平面向Z轴负向移动6cm
VSBW,1     ！用工作平面分离体
WPSTYLE,,,,,,,,0    ！重新设置工作平面
wpoff,0,0,12    ！将工作平面向Z轴正向移动6cm
FLST,2,2,6,ORDE,2
FITEM,2,2
FITEM,2,-3
VSBW,P51X    ！用工作平面分离体
wpoff,0,0,72    ！将工作平面向Z轴正向移动78cm
FLST,2,3,6,ORDE,3
FITEM,2,1
FITEM,2,-2
FITEM,2,4
VSBW,P51X    ！用工作平面分离体
WPSTYLE,,,,,,,,0    ！重新设置工作平面
wpoff,6,0,0     ！将工作平面沿WX方向平移6cm
wpro,,,-90.000000    ！将工作平面绕WY轴顺时针旋转90度
FLST,2,4,6,ORDE,2
FITEM,2,2
FITEM,2,-5
VSBW,P51X    ！用工作平面分离体
WPSTYLE,,,,,,,,0    ！重新设置工作平面
WPCSYS,-1,0    ！将当前坐标系X-Y平面定义为工作平面
BLOCK,0,250,78,90,0,1.5,  ！创建混凝土板
MAT,1     ！对炸药模型进行网格划分
REAL,
ESYS,0
SECNUM,
FLST,5,8,6,ORDE,3
FITEM,5,1
FITEM,5,6
FITEM,5,-12
VSEL,S,,,P51X
ASLV,S
LSLA,S
FLST,5,59,4,ORDE,14
FITEM,5,1
FITEM,5,3
FITEM,5,8
FITEM,5,-14
FITEM,5,16
FITEM,5,-17
FITEM,5,19
FITEM,5,-22
FITEM,5,24
FITEM,5,-25
FITEM,5,27
FITEM,5,-33
FITEM,5,35
FITEM,5,-69
CM,_Y,LINE
LSEL,,,,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,1.5,,,,,,,1    ！设置单元长度为1.5
MSHAPE,0,3D
MSHKEY,1
CM,_Y,VOLU
VSEL,,,,9
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
MAT,2      ！对土壤模型进行网格划分
REAL,
ESYS,0
SECNUM,
FLST,5,5,6,ORDE,4
FITEM,5,1
FITEM,5,6
FITEM,5,10
FITEM,5,-12
CM,_Y,VOLU
VSEL,,,,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
MAT,3      ！对空气模型进行网格划分
REAL,
ESYS,0
SECNUM,
FLST,5,2,6,ORDE,2
FITEM,5,7
FITEM,5,-8
CM,_Y,VOLU
VSEL,,,,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
MAT,4      ！对混泥土板模型进行网格划分
REAL,
ESYS,0
SECNUM,
VSEL,S,,,2
ASLV,S
LSLA,S
EPLOT      ！显示模型材料的有限元网格
FLST,5,12,4,ORDE,10
FITEM,5,2
FITEM,5,4
FITEM,5,-7
FITEM,5,15
FITEM,5,18
FITEM,5,23
FITEM,5,26
FITEM,5,34
FITEM,5,70
FITEM,5,-71
CM,_Y,LINE
LSEL,,,,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1,2,,,,,,,1    ！设置单元长度为2
CM,_Y,VOLU
VSEL,,,,2
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
EDPART,CREATE    ！创建PART
ALLSEL,ALL     ！选择所有实体
FLST,2,18,5,ORDE,14
FITEM,2,2
FITEM,2,-5
FITEM,2,28
FITEM,2,-29
FITEM,2,32
FITEM,2,-33
FITEM,2,37
FITEM,2,-38
FITEM,2,40
FITEM,2,42
FITEM,2,45
FITEM,2,-48
FITEM,2,50
FITEM,2,-51
DA,P51X,UZ,     ！对模型施加Z方向的约束
FLST,2,5,5,ORDE,5
FITEM,2,9
FITEM,2,11
FITEM,2,15
FITEM,2,19
FITEM,2,25
DA,P51X,UX,     ！对靠近炸药模型一侧施加X方向的约束
EDENERGY,1,0,1,1    ！设定分析选项
EDCPU,0,
EDBVIS,1.5,0.06,
TIME,100,     ！设置求解时间
EDCTS,0,0.6,     ！设置时间步
EDOPT,ADD,blank,LSDYNA   ！设置输出类型
EDRST,0,
EDHTIME,0,     ！设置时间间隔
EDDUMP,0,
ALLSEL,ALL     ！选择所有元素
EDWRITE,LSDYNA,explosion_undergroud,k ！生成K文件
/EXIT      ！退出ANSYS

6．K文件

*KEYWORD
*TITLE
EXPLOSION_UNDERGROUND.K
$                                                                   
*DATABASE_FORMAT
         0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                               NODE DEFINITIONS                          $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*NODE
       1 6.000000000E+00 6.000000000E+00 0.000000000E+00       0       0
       2 6.000000000E+00 6.000000000E+00 1.500000000E+00       0       0
       3 6.000000000E+00-6.000000000E+00 0.000000000E+00       0       0
       .................................................................
   91366 2.480000000E+02 8.400000000E+01 1.500000000E+00       0       0
   91367 2.480000000E+02 8.600000000E+01 1.500000000E+00       0       0
   91368 2.480000000E+02 8.800000000E+01 1.500000000E+00       0       0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             SECTION DEFINITIONS                        $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*SECTION_SOLID_ALE
1,11
*SECTION_SOLID_ALE
2,11
*SECTION_SOLID_ALE
3,11
*SECTION_SOLID
4,0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             ALE DEFINITIONS                             $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*ALE_MULTI-MATERIAL_GROUP
1,1
2,1
3,1
$
*CONSTRAINED_LAGRANGE_IN_SOLID
1,2,0,0,0,5,3,0
0,0,0.15
*SET_PART_LIST
         1
         4
*SET_PART_LIST
         2
         2         3
$
*CONTROL_ALE
2,1,2,-1.0000000,0.0000000,0.0000000,0.0000000
0.0000000,0.0000000,0.0000000
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             MATERIAL DEFINITIONS                        $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*MAT_HIGH_EXPLOSIVE_BURN
1,1.93,0.993,0.370,0.000E+00       
*EOS_JWL
1,3.71,7.430E-02,4.15,0.950,0.300,7.000E-02,1.00            
$
*INITIAL_DETONATION
1,0.00,0.00,0.00,0.00
*INITIAL_DETONATION
1,0.00,0.00,2.0,0.00
$
*MAT_SOIL_AND_FOAM
2,1.80,1.601E-04,1.328E+02,3.300E-03,1.310E-07,0.1232,0.0   
0.0,0.0
0.0,0.050,0.090,0.11,0.15,0.19,0.21,0.22
0.25,0.30
0.0,3.420E-02,4.530E-02,6.760E-02,1.270E-01,2.080E-01,2.710E-01,3.920E-01
5.660E-01,1.230E+00 
$  
*MAT_NULL
3,0.125E-02,0.00,0.00,0.00,0.00,0.00,0.00   
*EOS_GRUNEISEN
3,0.344,0.00,0.00,0.00,1.40,0.00,0.00   
0.00   
$
*MAT_PLASTIC_KINEMATIC
4,2.6500000,4.000E-01,0.3000000,1.000E-03,4.000E-02,0.5000000   
0.0000000,0.0000000,0.8000000   
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              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         0         0         0         0
$
*PART
Part          3 for Mat         3 and Elem Type         3
         3         3         3         3         0         0         0
$
*PART
Part          4 for Mat         4 and Elem Type         4
         4         4         4         0         0         0         0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             ELEMENT DEFINITIONS                         $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*ELEMENT_SOLID
 1      1       2       1      10      18      22      69      76      34
 2      1      18      10       9      17      34      76      75      35
 3      1      17       9       8      16      35      75      74      36
       .........................................................................
   45118       4   90484   90485   89614   89615   91358   91363   90621   90622
   45119       4   90485   90486   89613   89614   91363   91368   90620   90621
   45120       4   90486   89742   89612   89613   91368   90618   90613   90620
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              COORDINATE SYSTEMS                        $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                               LOAD DEFINITIONS                          $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                                RIGID BOUNDRIES                          $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                             BOUNDARY DEFINITIONS                        $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*SET_NODE_LIST
        1     0.000     0.000     0.000     0.000
        1       2         3         4         5         6         7         8
        9       10        11        12        13        14        15        16
       17       18        20        21        22        31        32        33
        ........................................................................
     91347     91348     91349     91350     91351     91352     91353     91354
     91355     91356     91357     91358     91359     91360     91361     91362
     91363     91364     91365     91366     91367     91368
*BOUNDARY_SPC_SET
        1       0         0         0         1         0         0         0
*SET_NODE_LIST
        2     0.000     0.000     0.000     0.000
       19       23        24        25        26        27        28        29
       30       55        59        60        61        62        63        64
        ........................................................................
     89478     89479     89480     89481     89605     89606     89607     89608
     89609     89610     89611     90487     90619     90744     90745     90746
     90747     90748
*BOUNDARY_SPC_SET
     2         0         1         0         1         0         0         0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              CONTACT DEFINITIONS                        $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              CONTROL OPTIONS                            $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*CONTROL_CPU
  0.00   
*CONTROL_ENERGY
         2         1         2         2
*CONTROL_SHELL
  20.0             1        -1         1         2         2         1
*CONTROL_TIMESTEP
    0.0000    0.6000         0  0.00      0.00   
*CONTROL_TERMINATION
 0.300E+04         0   0.00000   0.00000   0.00000
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                               TIME HISTORY                              $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_BINARY_D3PLOT
 20.00  
*DATABASE_BINARY_D3THDT
 4.000   
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              DATABASE OPTIONS                           $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_EXTENT_BINARY
         0       0         3         1         0         0         0         0
         0       0         4         0         0         0
*END