ls-dyna的一个问题
最近做一个练习，模拟一个物体从一米高的地方掉下来，但是我的结果总是坠落体穿透了固定板，帮忙看看。finish
/clear
/filname,droping_test
/PREP7
*SET,m_length,99e-3 !坠落体形状
*SET,m_width,45e-3
*SET,m_thick,23e-3
*SET,high,1 !坠落高度
*SET,g_length,0.2 !固定板尺寸
*SET,g_width,0.1
*SET,g_thick,0.05
ET,1,SOLID164 !选择单元
MP,DENS,1,7800 !坠落体材料
MP,EX,1,2e11
MP,NUXY,1,0.3
EDMP,RIGI,2,7,7 !固定板材料
MP,DENS,2,7800
MP,EX,2,2e11
MP,NUXY,2,0.3
block,-m_length/2,m_length/2,-m_width/2,m_width/2,high,high+m_thick !坠落体建模
block,-g_length/2,g_length/2,-g_width/2,g_width/2,-g_thick,0 !固定板
/vup,1,z
/view,1,1,1,1 !单元大小控制
lsel,s,loc,x,0
lsel,r,loc,z,high,high+m_thick
lesize,all,,,10
lsel,s,loc,y,0
lsel,r,loc,z,high,high+m_thick
lesize,all,,,5
lsel,s,loc,z,m_thick/2+high
lesize,all,,,3
allsel
lsel,s,loc,z,-g_thick,0
lesize,all,9e-3
allsel
gplot
type,1 !划分单元
mat,1
real,1
vmesh,1
type,1
mat,2
real,1
vmesh,2
*dim,time,array,2 !构建加载数组
*dim,accelerat_g,array,2
time(1)=0,1
accelerat_g(1)=-9.8,-9.8
nsel,s,loc,z,high,high+m_thick
/MREP,EPLOT
CM,droping,NODE
nsel,s,loc,z,-g_thick,0 !创建坠落体节点组
cm,groud,node
allsel
EDPART,CREATE !creat parts
EDBX,Add,1,-m_length/2,m_length/2,-m_width/2,m_width/2,high,high+m_thick,
EDBX,Add,2,-g_length/2,g_length/2,-g_width/2,g_width/2,-g_thick,0,
ALLSEL
EDCGEN,ASTS, 2, 1,0.35,0.29,0,0,0, , , , ,0,10000000,2,1 !定义接触？？
EDCONTACT,0.9,0,2,0,1,1,1,0.4,1
/SOL
EDLOAD,ADD,ACLZ,0,DROPING,TIME,ACCELERAT_G, 0, , , , ,
TIME,0.6,
EDCTS,0,0.9,
EDRST,100,
EDHTIME,1000,
EDDUMP,1,
EDENERGY,1,1,1,1
EDBVIS,1.5,0.06,
SOLVE 我也初学，感觉是不是定义接触有问题啊 我也觉得是接触的问题，如果用自DTM做的话就没什么问题，但是自己定义接触就出问题了。
页:
[1]