偶正在做一个东东,就是炮弹冲击装甲材料(本例为混凝土)
可是在ansys6.1上运行时总出错,我得命令流可能有错,麻烦各位大虾或管理员老哥帮忙看看,不胜感激
k文件见附件:
命令流:
!炮弹冲击混凝土材料
finish
/cle
/filename,bullet impacts concrete plate
/title,Impact of a Bullet into a Concrete Plate
/PREP7
/view,,1,2,3
!*
! --- Enter steel Material Properties For Bullet
et,1,mesh200
keyopt,1,1,6 !好像落了个1
mp,ex,1,210e9
mp,nuxy,1,.27
mp,dens,1,7800
!*
et,2,solid164 !定义单元类型(3D)
mp,ex,2,210e9
mp,nuxy,2,.27
mp,dens,2,7800
!*
! --- Create and Mesh Bullet
k,50000,0,0,9e-2
k,50001,0,10,9e-2
circle,50000,8e-2,50001,,90 ! 以第 1 个keypoin为圆心画弧线 ,半径为8e-2,画出1/4个圆
kdele,50000,50001
/replot
k,3,8e-2,0,81e-2
k,4,0,,81e-2
k,5,0,0,9e-2
l,1,3
l,3,4
l,4,5
l,5,2
l,1,5
al,1,5,6
al,2,3,4,6
esize,0.015
mshape,0,2d
mshkey,1
amesh,all
vrotat,1,2,,,,,2,4,360,4
/pnum,volu,1
/number,0
vplot,all
mat,2
type,2
vmesh,all
cm,nbullet,node
cm,ebullet,elem
!*
! --- 对炮弹定义初始速度
edvel,vgen,nbullet,,,-1000.0
!*
! --- 混凝土板的材料参数
et,3,solid164 !定义单元类型
mp,ex,3,20e9
mp,nuxy,3,.3
mp,dens,3,2800
! --- 混凝土板建立模型
block,2,-2,2,-2,-0.1,-1.1
esize,.1
mat,3
type,3
vmesh,9
esel,s,mat,,3
nsle
cm,nplate,node
cm,eplate,elem
nsel,all
esel,all
! --- 对混凝土板的底边施加约束
nsel,s,loc,z,-1.1
nsel,u,loc,x,-1.99,1.99
d,all,ux,0,,,,uy,uz
nsel,s,loc,z,-1.1
nsel,u,loc,y,-1.99,1.99
d,all,ux,0,,,,uy,uz
nsel,all
/replot
eplot
! --- 定义接触 Eroding Surface-To-Surface Contact Between Bullet and Plate
edcgen,ests,nbullet,nplate,,,,,,1,1,1 ! Set V1=V2=V3=1
edcontact,1.0,,,,,,2,,0 ! SFSI, ORIE, and PENCHK explicitly set
! --- 对后处理POST26定义集合
nsel,s,loc,x,-2,2
nsel,r,loc,z,-1.1,0.81
nsel,r,loc,y,-2,2
esln
nsle
cm,elems,elem
cm,nodes,node
nsel,all
esel,all
eplot
finish
/solu
time,.002 ! Endtime
edrst,10 ! Time Increment For Results Output
edopt,,,lsdyna
edhtime,100
edhist,nodes
edhist,elems
eplot
save
! tb材料定义
finish
/prep7
!*
!*
tb,plaw,3,,,1
tbdata,1,289e6,
tbdata,2,500e6,
tbdata,3,0,
tbdata,4,40,
tbdata,5,5,
tbdata,6,0.125,
!*
save
finish
!求解
/solu
/status,solu
solve |