tixudaodao 发表于 2006-12-4 14:07

请教:能否用ansys ls-dyna模拟岩石(落石)冲击土体过程?

各位老大,请教:能否用ansys ls-dyna模拟岩石(落石)冲击土体过程?
能否得到冲击力、土体应力、变形等参数?模拟过程中需要注意哪些问题?
谢谢~~

wanglinfeng 发表于 2010-10-5 09:49

我也很想知道啊{:{13}:}

mbcy090101 发表于 2010-11-21 21:39

应该可以,这个问题和动力学侵彻问题一样,给石头块施加一个初速度,定义好石头和土坯的材料参数最重要。LS-DYNA的后处理器中可以直观的看出应力,压力云图以及土坯的破坏情况,如果你要看的比较细,建议你使用配置较好的电脑,选择小一点的动画步长。

Chelsea 发表于 2010-11-22 09:49

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

huizi42771 发表于 2011-9-20 14:54

学习了!!!!!!!
页: [1]
查看完整版本: 请教:能否用ansys ls-dyna模拟岩石(落石)冲击土体过程?