找回密码
 注册
Simdroid-非首页
查看: 256|回复: 0

【讨论】自已写的钢弹穿甲铝板命令流

[复制链接]
发表于 2003-12-5 23:33:42 | 显示全部楼层 |阅读模式 来自 陕西西安
可惜结果不好,希望哪位高手能指正,现在的小目标就是让钢弹不穿透铝板。急!
/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  

评分

1

查看全部评分

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-28 18:30 , Processed in 0.048971 second(s), 15 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表