ls-dyna的一个问题

liyun99

最近做一个练习，模拟一个物体从一米高的地方掉下来，但是我的结果总是坠落体穿透了固定板，帮忙看看。
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

jpy2008

我也初学，感觉是不是定义接触有问题啊

liyun99

我也觉得是接触的问题，如果用自DTM做的话就没什么问题，但是自己定义接触就出问题了。