ESTS.rar_ansys_ansys lsdyna分析_ansys碰撞_lsdyna_全自动脚本程序

/FILNAM,ESTSanalysis /TITLE,Analysis of of ESTSanalysis /PREP7 *afun,deg !定义成角度制 !输入常数参数 Z1=0.4 !Z方向移动距离 !弹体模型选择参数 CHO1=1 !靶板模型选择参数 CHO2=15 !材料选择号 CHOMAT1=1 CHOMAT2=2 CHOMAT3=3 CHOMAT4=4 CHOMAT5=1 CHOMAT6=2 CHOMAT7=3 CHOMAT8=4 CHOMAT9=1 CHOMAT10=2 CHOMAT11=3 CHOMAT12=4 CHOMAT13=1 CHOMAT14=2 CHOMAT15=1 !模型凸台参数 D=1.3 D1=0.7 theat=45 HH1=2 !模型方块参数 D2=1.3 HH2=2 !模型圆柱参数 D3=1.3 HH3=2 !靶板尺寸 L1=8 H1=0.3 H2=0.3 H3=0.3 H4=0.3 H5=0.3 H6=0.3 H7=0.3 H8=0.3 H9=0.3 H10=0.3 H11=0.3 H12=0.3 H13=0.3 H14=0.3 H15=0.3 !网格尺寸,分别为弹头和靶板尺寸 size1=0.15 size2=0.1 !子弹速度 VEL1=0.13 !求解时间 TIME1=100 !求解步数 SUB1=50 !CPU数目 numcpus=4 !*钢板的材料参数 !密度 DENS1=7.83 !弹性模量 EX1=2 !泊松比 NUXY1=0.32 !屈服应力常数 TBDAT11=0.00792 !应变硬化常数 TBDAT12=0.0051 !应变硬化指数 TBDAT13=0.26 !应变率相关系数 TBDAT14=0.014 !温度相关系数 TBDAT15=1.03 !熔化温度 TBDAT16=1793 !室温 TBDAT17=294 !Ef P1 SR TBDAT18=1E-006 !Specific Heat(无单位) TBDAT19=4.77E-006 !失效应力 TBDAT110=-9 !第一断裂参数 TBDAT111=0.8 !第二断裂参数 TBDAT112=0 !第三断裂参数 TBDAT113=0 !第四断裂参数 TBDAT114=0 !第五断裂参数 TBDAT115=0 !波速 TBDAT116=0.4569 !状态方程参数S1 TBDAT117=1.49 !状态方程参数S2 TBDAT118=0 !状态方程参数S3 TBDAT119=0 !Gruneisen系数 TBDAT120=2.17 !一阶体积修正系数 TBDAT121=0.46 !E(0)(无单位) TBDAT122=0 !V(0)(无单位) TBDAT123=1 !*陶瓷的材料参数 !密度 DENS2=7.84 !弹性模量 EX2=2 !泊松比 NUXY2=0.32 !屈服应力常数 TBDAT21=0.00792 !应变硬化常数 TBDAT22=0.0051 !应变硬化指数 TBDAT23=0.26 !应变率相关系数 TBDAT24=0.014 !温度相关系数 TBDAT25=1.03 !熔化温度 TBDAT26=1793 !室温 TBDAT27=294 !Ef P1 SR TBDAT28=1E-006 !Specific Heat(无单位) TBDAT29=4.77E-006 !失效应力 TBDAT210=-9 !第一断裂参数 TBDAT211=0.8 !第二断裂参数 TBDAT212=0 !第三断裂参数 TBDAT213=0 !第四断裂参数 TBDAT214=0 !第五断裂参数 TBDAT215=0 !波速 TBDAT216=0.4569 !状态方程参数S1 TBDAT217=1.49 !状态方程参数S2 TBDAT218=0 !状态方程参数S3 TBDAT219=0 !Gruneisen系数 TBDAT220=2.17 !一阶体积修正系数 TBDAT221=0.46 !E(0)(无单位) TBDAT222=0 !V(0)(无单位) TBDAT123=1 !*带损伤复合材料参数 DENS3=1.85 EX3=0.238 EY3=0.09 EZ3=0.09 GXY3=0.0165 GYZ3=0.0302 GXZ3=0.0165 NUXY3=0.31 NUYZ3=0.49 NUXZ3=0.31 !压缩失效体积模量EBULK TBDAT31=0.204 !剪切强度S12 TBDAT32=0.0039 !纵向拉伸强度S1 TBDAT33=0.013 !横向拉伸强度S2 TBDAT34=0.0013 !横向压缩强度C2 TBDAT35=0.0065 !非线性剪应力参数A TBDAT36=0 !*另一种可描述金属损伤的模型参数 !密度 DENS4=8.33 !弹性模量 EX4=2 !泊松比 NUXY4=0.32 !屈服应力 TBDAT41=0.0025 !剪切模量 TBDAT42=0.015 !硬化参数 TBDAT43=0 !应变率(c) TBDAT44=40 !应变率(p) TBDAT45=5 !失效应变 TBDAT46=0.75 !*定义弹丸的材料模型参数 DENS5=8.93 EX5=1.17 NUXY5=0.35 M1=H1 M2=M1+H2 M3=M2+H3 M4=M3+H4 M5=M4+H5 M6=M5+H6 M7=M6+H7 M8=M7+H8 M9=M8+H9 M10=M9+H10 M11=M10+H11 M12=M11+H12 M13=M12+H13 M14=M13+H14 M15=M14+H15 !*定义单元类型 /PREP7 ET,1,SOLID164 !*定义钢板的材料模型 MP,DENS,1,DENS1 MP,EX,1,EX1 MP,NUXY,1,NUXY1 TB,EOS,1,,,1,2 TBDAT,1,TBDAT11 TBDAT,2,TBDAT12 TBDAT,3,TBDAT13 TBDAT,4,TBDAT14 TBDAT,5,TBDAT15 TBDAT,6,TBDAT16 TBDAT,7,TBDAT17 TBDAT,8,TBDAT18 TBDAT,9,TBDAT19 TBDAT,10,TBDAT110 TBDAT,11,TBDAT111 TBDAT,12,TBDAT112 TBDAT,13,TBDAT113 TBDAT,14,TBDAT114 TBDAT,15,TBDAT115 TBDAT,16,TBDAT116 TBDAT,17,TBDAT117 TBDAT,18,TBDAT118 TBDAT,19,TBDAT119 TBDAT,20,TBDAT120 TBDAT,21,TBDAT121 TBDAT,22,TBDAT122 TBDAT,23,TBDAT123 !*定义陶瓷的材料模型 MP,DENS,2,DENS2 MP,EX,2,EX2 MP,NUXY,2,NUXY2 TB,EOS,2,,,1,2 TBDAT,1,TBDAT21 TBDAT,2,TBDAT22 TBDAT,3,TBDAT23 TBDAT,4,TBDAT24 TBDAT,5,TBDAT25 TBDAT,6,TBDAT26 TBDAT,7,TBDAT27 TBDAT,8,TBDAT28 TBDAT,9,TBDAT29 TBDAT,10,TBDAT210 TBDAT,11,TBDAT211 TBDAT,12,TBDAT212 TBDAT,13,TBDAT213 TBDAT,14,TBDAT214 TBDAT,15,TBDAT215 TBDAT,16,TBDAT216 TBDAT,17,TBDAT217 TBDAT,18,TBDAT218 TBDAT,19,TBDAT219 TBDAT,20,TBDAT220 TBDAT,21,TBDAT221 TBDAT,22,TBDAT222 TBDAT,23,TBDAT123 !*定义带损伤复合材料的材料模型 MP,DENS,3,DENS3 MP,EX,3,EX3 MP,EY,3,EY3 MP,EZ,3,EZ3 MP,GXY,3,GXY3 MP,GYZ,3,GYZ3 MP,GXZ,3,GXZ3 MP,NUXY,3,NUXY3 MP,NUYZ,3,NUYZ3 MP,NUXZ,3,NUXZ3 TB,COMP,3,,,, TBDAT,1,TBDAT31 TBDAT,2,TBDAT32 TBDAT,3,TBDAT33 TBDAT,4,TBDAT34 TBDAT,5,TBDAT35 TBDAT,6,TBDAT36 !*另一种可描述金属损伤的模型材料模型 MP,DENS,4,DENS4 MP,EX,4,EX4 MP,NUXY,4,NUXY4 TB,PLAW,4,,,1 TBDATA,1,TBDAT41 TBDATA,2,TBDAT42 TBDATA,3,TBDAT43 TBDATA,4,TBDAT44 TBDATA,5,TBDAT45 TBDATA,6,TBDAT46 !*定义刚体弹丸的材料模型 EDMP,RIGI,5,4,7 MP,DENS,5,DENS5 MP,EX,5,EX5 MP,NUXY,5,NUXY5 !开始建模 !*建立弹丸实体模型 *IF,CHO1,EQ,1,THEN !凸台模型 !分别饶z,x,y旋转工作平面角度，正向为XY,YZ,ZX !LOCAL,11,0,0,0,0, ,90, ,1,1, CYL4,0,0,D/2, , , ,HH1 WPOFFS,D1/2,0,0 WPROTA,,,-theat VSBW,1 WPAVE,0,0,0 CSYS,0 WPCSYS,-1,0 WPOFFS,-D1/2,0,0 WPROTA,,,theat VSBW,3 VDELE,1,,,1 VDELE,2,,,1 WPAVE,0,0,0 CSYS,0 WPCSYS,-1,0 VGEN, ,ALL, , , , ,Z1, , ,1 NUMCMP,ALL NUMMRG,ALL *ELSEIF,CHO1,EQ,2,THEN !方块模型 BLOCK,-D2/2,D2/2,-D2/2,D2/2,0,HH2, VGEN, ,ALL, , , , ,Z1, , ,1 NUMCMP,ALL NUMMRG,ALL *ELSE,CHO1,EQ,3,THEN !模型圆柱 CYL4,0,0,D3/2 , , ,HH3 VGEN, ,ALL, , , , ,Z1, , ,1 NUMCMP,ALL NUMMRG,ALL *ENDIF !赋予材料，划分网格 !弹头网格划分 MSHKEY,0 !1采用映射划分，0自由，2先映射在自由 MSHAPE,1 !0生成六面体，1四面体 VSEL,S,,,1 !VATT,MAT,REAL,TYPE,ESYS VATT,5,,1 ESIZE,SIZE1 VMESH,ALL ALLSEL,ALL NUMCMP,ALL NUMMRG,ALL *IF,CHO2,EQ,1,THEN !*建立靶板实体模型1 BLOCK,-L1/2,L1/2,-L1/2,L1/2,0,-H1, NUMCMP,ALL NUMMRG,ALL MSHKEY,1 MSHAPE,0 VSEL,S,LOC,Z,-M1/2.0 !VATT,MAT,REAL,TYPE,ESYS VATT,CHOMAT1,,1 ESIZE,SIZE2 VMESH,ALL ALLSEL,ALL !设置接触信息 EDPART,CREATE EDCGEN,ESTS,2,1,0,0,0,0,0,0,1,0, ,0,10000000,0,0 EDCONTACT,10,0,2,0,1,1,1,4,0 *ELSEIF,CHO2,EQ,2,THEN !*建立靶板实体模型2 BLOCK,-L1/2,L1/2,-L1/2,L1/2,0,-M1, BLOCK,-L1/2,L1/2,-L1/2,L1/2,-M1,-M2, VGLUE,ALL NUMCMP,ALL NUMMRG,ALL MSHKEY,1 MSHAPE,0 VSEL,S,LOC,Z,-M1/2.0 !VATT,MAT,REAL,TYPE,ESYS VATT,CHOMAT1,,1 ESIZE,SIZE2 VMESH,ALL ALLSEL,ALL MSHKEY,1 MSHAPE,0 VSEL,S,LOC,Z,-(M1+H2/2.0) !VATT,MAT,REAL,TYPE,ESYS VATT,CHOMAT2,,1 ESIZE,SIZE2 VMESH,ALL ALLSEL,ALL !设置接触信息 EDPART,CREATE EDCGEN,ESTS,2,1,0,0,0,0,0,0,1,0, ,0,10000000,0,0 EDCONTACT,10,0,2,0,1,1,1,4,0 EDCGEN,ESTS,3,1,0,0,0,0,0,0,1,0, ,0,10000000,0,0 EDCONTACT,10,0,2,0,1,1,1,4,0 !设置层合板间接触 *ELSEIF,CHO2,EQ,3,THEN !*建立靶板实体模型3 BLOCK,-L1/2,L1/2,-L1/2,L1/2,0,-M1, BLOCK,-L1/2,L1/2,-L1/2,L1/2,-M1,-M2, BLOCK,-L1/2,L1/2,-L1/2,L1/2,-M2,-M3, VGLUE,ALL NUMCMP,ALL NUMMRG,ALL MSHKEY,1 MSHAPE,0 VSEL,S,LOC,Z,-M1/2.0 !VATT,MAT,REAL,TYPE,ESYS VATT,CHOMAT1,,1 ESIZE,SIZE2 VMESH,ALL ALLSEL,ALL MSHKEY,1 MSHAPE,0 VSEL,S,LOC,Z,-(M1+H2/2.0) !VATT,MAT,REAL,TYPE,ESYS VATT,CHOMAT2,,1 ESIZE,SIZE2 VMESH,ALL ALLSEL,ALL MSHKEY,1 MSHAPE,0 VSEL,S,LOC,Z,-(M2+H3/2.0) !VATT,MAT,REAL,TYPE,ESYS VATT,CHOMAT3,,1 ESIZE,SIZE2 VMESH,ALL ALLSEL,ALL !设置接触信息 EDPART,CREATE EDCGEN,ESTS,2,1,0,0,0,0,0,0,1,0, ,0,10000000,0,0 EDCONTACT,10,0,2,0,1,1,1,4,0 EDCGEN,ESTS,3,1,0,0,0,0,0,0,1,0, ,0,10000000,0,0 EDCONTACT,10,0,2,0,1,1,1,4,0 EDCGEN,ESTS,4,1,0,0,0,0,0,0,1,0, ,0,10000000,0,0 EDCONTACT,10,0,2,0,1,1,1,4,0 !设置层合板间接触 *ELSEIF,CHO2,EQ,4,THEN !*建立靶板实体模型4 BLOCK,-L1/2,L1/2,-L1/2,L1/2,0,-M1, BLOCK,-L1/2,L1/2,-L1/2,L1/2,-M1,-M2, BLOCK,-L1/2,L1/2,-L1/2,L1/2,-M2,-M3, BLOCK,-L1/2,L1/2,-L1/2,L1/2,-M3,-M4, VGLUE,ALL NUMCMP,ALL NUMMRG,ALL MSHKEY,1 MSHAPE,0 VSEL,S,LOC,Z,-M1/2.0 !VATT,MAT,REAL,TYPE,ESYS VATT,CHOMAT1,,1 ESIZE,SIZE2 VMESH,ALL ALLSEL,ALL MSHKEY,1 MSHAPE,0 VSEL,S,LOC,