前段时间利用ANSYS/LS-DYNA模拟落石冲击地面过程，能够通过加载初始竖向速度实现自由落体运动的模拟。在分析中，存在几点疑惑。希望能够得到大家的帮助。
       （1）无法模拟落石以任意角度入射地面的碰撞过程，原因在于只能加载初始竖向速度，不能加载初始水平速度，不知道为何？
       （2）模拟中假定落石为刚体，地面为弹塑性体，塑性部分采用的是BISO准则，其中相关塑性参数采用屈服应力和切线模量，不是我们常规使用的MC或者DP，这两个参数如何确定也是麻烦，或者与DP或MC准则参数是否存在转换关系？一头雾水！
       （3）ANSYS/LS-DYNA模拟落石冲击地面模拟中，如果地面是斜坡，程序不能模拟落石碰撞过程中的翻转，落石在整个过程中是平移，与实际存在差别。
     以下是我的程序，现附在后面，希望能够得到高手指点！
/title,Test No.A028
/PREP7
et,1,solid164                     !ANSYS/LS-DYNA采用显式求解，ANSYS采用的是隐式求解
et,2,solid164
et,3,shell163

H=2.95                                     !落石高度
Lside=0.30                              !落石边长
LX=5.0*Lside                           !模型X轴边界范围
LY=5.0*Lside                           !模型Y轴边界范围
LZ=5.0*Lside                           !模型Z轴边界范围
EL=20
vector=sqrt(2*9.8*H)                    !m/s
TT=0.006
DTT=1000                                 !SQRT(2*h/9.8)+3.0
DENS2=2390                              !落石的密度*9.8
FORC1=DENS2*Lside**3*9.8                !落石的自重

!地面岩土体--碎石土
MP,EX,1,1.0e9
MP,NUXY,1,0.24
MP,DENS,1,2100
TB,BISO
TBDATA,1,1.2e6,4.0e5              !屈服应力,越大回弹越高;切线模量,越大回弹越高

!落石
MP,EX,2,5.0e10
MP,NUXY,2,0.24
MP,DENS,2,DENS2

k,1,0.0,LY,-LZ
k,2,0.0,0,-LZ
k,3,2*LX,LY,-LZ
k,4,2*LX,0,-LZ
K,5,0.0,0.0,LZ
A,1,2,4,3

LESIZE,all,Lside/2,,,, , , ,1
amesh,all
l,2,5
TYPE,1   
EXTOPT,ESIZE,EL,1,  
EXTOPT,ACLEAR,1
EXTOPT,ATTR,0,0,0   
MAT,1   
REAL,1  
ESYS,0  
VDRAG,1,,,,,,5
cm,nsoil,node
cm,esoil,elem
ldele,5,,,1
allsel,all
!建立落石几何模型
block,LX-Lside/2,LX+Lside/2,LY,LY+Lside,-Lside/2,Lside/2
vsel,,,,2
vatt,2,2,2
aslv,s
lsla,s
LESIZE,ALL,,,2,1
vmesh,2
eslv,s
nsle,s
cm,nrock,node
cm,erock,elem
EDMP,RIGID,2,6,7      !用 EDMP命令定义刚性体,6指约束住x和z方向，7指约束住3个旋转自由度
aLLsel,aLL
*DIM,TIME,ARRAY,2,1,1,,,
*DIM,FORC,ARRAY,2,1,1,,,
*SET,FORC(1,1,1),-FORC1
*SET,FORC(2,1,1),-FORC1
*SET,TIME(2,1,1),TT

asel,,,,3
nsla,s,1
D,aLL,,,,,,ALL,,,,,
aLLsel,aLL
asel,,,,1,6,5
nsla,s,1
D,aLL,UZ,,,,,,,,,,
aLLsel,aLL
asel,,,,2,4,2
nsla,s,1
D,aLL,UX,,,,,,,,,,
aLLsel,aLL
EPLOT
EDPART,CREATE,1,esoil
EDPART,CREATE,2,erock
EDCGEN,STS,2,1,0.00,0.00,0.0,0,0.00,,,,,0,,0,0   !定义接触参数,粘滞阻尼系数VC=sy/sqrt(3),VDC  
!EDCONTACT,0.1,0,2,2,1,1,1,4,0                   !接触面控制
EDHGLS,0.10
EDCSC,ON
EDLOAD,ADD,RBFY,,2,TIME,FORC,0,,,,,
EDPVEL,VELO,2,vector,-vector,,0,0,0    !定义初始速度可以采用edvel或edpvel,X方向的速度没加载成功，不知道原因？？？？？
FINISH

/SOLU
antype,trans
TIMINT,ON
TRNOPT,FULL
TIME,TT,
EDRST,DTT,
EDHTIME,DTT,
EDDUMP,1,
!/STATUS,SOLU
SOLVE
FINISH
/DSCALE,1,OFF   
/REPLOT

vsel,,,,2
eslv,s
nsle,s
*get,Nrock,node,0,num,min
allsel,all
!落石运动特性与时间曲线
/POST26
NSOL,2,Nrock,U,Y,UY                           !块石Y方向位移的时程曲线
NSOL,3,Nrock,V,Y,VY                           !块石Y方向速度的时程曲线
NSOL,4,Nrock,A,Y,AY                           !块石Y方向加速度的时程曲线
NSOL,5,Nrock,U,X,UX                           !块石X方向位移的时程曲线
NSOL,6,Nrock,V,X,VX                           !块石X方向速度的时程曲线
NSOL,7,Nrock,A,X,AX                           !块石X方向加速度的时程曲线
STORE,MERGE  
XVAR,1                                        !定义X轴变量，1代表时间，n代表所定义的变量
PLVAR,2,5
PLVAR,3,6
PLVAR,4,7
PRVAR,2,3,4

转载本文请联系原作者获取授权，同时请注明本文来自章广成科学网博客。
链接地址：https://wap.sciencenet.cn/blog-304483-333136.html