（2）输出土壤鼓包运动的速度及位移时程曲线 1）运行LS-PREPOST，读入d3plot结果数据文件。

图9-13爆炸后土壤的鼓包及混凝土板运动过程

2）选择主菜单SelPar，在Part Selection中选择H4（只显示Part4）。

3）选择主单History→Nodal→Y-velocity命令，在图形对话框中选择H4对称面上最下方的一个节点（节点号为89606），如图9-14所示。单击Plot按钮，将弹出Plot Window-1对话框，输出混凝土板顶点的Y方向速度时程曲线，如图9-15所示。

图9-14节点89606

4）选择菜单History→Nodal→Y-displacement命令，在图形对话框中仍选择节点89606，单击Plot按钮，将弹出Plot Window-1对话框，输出混凝土板顶点的Y方向位移时程曲线，如图9-16所示。

图9-15鼓包运动的速度时程曲线

图9-16鼓包运动的位移时程曲线

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