可惜结果不好,希望哪位高手能指正,现在的小目标就是让钢弹不穿透铝板。急!
/com,弹为弹性。1/4
/com,**子弹侵彻靶板分析( 三维)
/filenam,3DERODING
/title,Projectile Penetrating Plate
/prep7
H=25.4*0.001
ZD=12.7*0.001 !子弹直径
ZL=47*0.001 !子弹长度
spac=8
high=0.0008 !单元厚度
Length=110*0.001 !板长度
/com,----单元及实常数
ET,1,solid164
ET,2,solid164
/com,子弹模型
MP,ex,1,200e9 ! Pa
MP,nuxy,1,.3 ! No units
MP,dens,1,7800 ! kg/m3
/com,Define materials 材料(纯铝)
MP,ex,2,72e9 ! Pa
MP,nuxy,2,.3 ! No units
MP,dens,2,2800 ! kg/m3
TB,PLAW,2,,,1
TBDATA,1,35e6 ! Yield stress (Pa)
TBDATA,2,140e6 ! Tangent modulus (Pa)
TBDATA,4,1.47e5 ! C (s-1)
TBDATA,5,2.8 ! P
TBDATA,6,0.8 ! Failure strain
/com,定义子弹外形/网格/初速度
K,1000,-2.5*0.001,0,0
K,1001,-2.5*0.001-ZD/2.0,ZD/2.0,0
K,1002,-2.5*0.001-ZD/2.0-ZL,ZD/2.0,0
K,1005,-2.5*0.001-ZD/2.0,0,0
K,1006,-2.5*0.001-ZD/2.0-ZL,0,0
KWPAVE,1000,1001,1002
LARC,1000,1001,1005,ZD/2.0
L,1001,1002
L,1002,1006
L,1006,1000
AL,1,2,3,4
esize,0.001
LESIZE,1,,,12
LESIZE,2,,,30
LESIZE,3,,,12
LESIZE,4,,,30
VROTAT,1,,,,,,1005,1006,90
type,2
mat,1
vmesh,all
ESEL,S,ENAME,,164 !选取单元类型164为Active,
NSLE,S !选取——选取单元上的——全部节点
CM,projectile,NODE !节点成组——名称projectile
EDIVELO,projectile,193,,,0,0,0,0,0,0,0, !指定初始速度,x方向193m/s
alls
/com,板模型
K,101,0 ,Length/2.0,0
K,102,H ,Length/2.0,0
K,103,0 ,0,0
K,104,H ,0,0
A,103,104,102,101
lesize,13,,,divi,1/spac, , , ,0
lesize,11,,,divi,spac, , , ,0
LESIZE,10,high, , , , , , ,1
LESIZE,12,high, , , , , , ,1
VROTAT,6,,,,,,1005,1006,90
mat,2
vmesh,2
/com,约束条件
local,11,1,0,0,0,0,0,-90
csys,11
nsel,s,loc,x,Length/2.0
d,all,all,0.0
csys,0
allsel,all
nsel,s,loc,y,0
DSYM,SYMM,Y
nsel,s,loc,z,0
dsym,symm,z
alls
/com,----时间控制参数
time,0.0004
edrst,100 !.RST文件的文件输出时间间隔
edhtime,100 !.HIS文件的文件输出时间间隔
edenergy,1,1,1,1 !指定能量耗散控制
save
/com,----创建PART
nsel,s,loc,x,0,H+0.1*0.001
CM,target,NODE
EDCGEN,ASTS,PROJECTILE,TARGET,0,0,0,0,0, , , , ,0,10000000,0,0
ALLSEL,ALL
FINISH
/SOLU
/STAT,SOLU
/REPLOT
EDOPT,ADD,blank,LSDYNA
/STATUS,SOLU
SOLVE
FINISH |