请教：能否用ansys ls-dyna模拟岩石（落石）冲击土体过程？

tixudaodao

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

wanglinfeng

我也很想知道啊

mbcy090101

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

Chelsea

ANSYS/LS-DYNA模拟落石冲击地面的案例

1. /title,Test No.A028
   
2. /PREP7
   
3. et,1,solid164                     !ANSYS/LS-DYNA采用显式求解，ANSYS采用的是隐式求解
   
4. et,2,solid164
   
5. et,3,shell163
   
6. 
   
7. H=2.95                                     !落石高度
   
8. Lside=0.30                              !落石边长
   
9. LX=5.0*Lside                           !模型X轴边界范围
   
10. LY=5.0*Lside                           !模型Y轴边界范围
    
11. LZ=5.0*Lside                           !模型Z轴边界范围
    
12. EL=20
    
13. vector=sqrt(2*9.8*H)                    !m/s
    
14. TT=0.006
    
15. DTT=1000                                 !SQRT(2*h/9.8)+3.0
    
16. DENS2=2390                              !落石的密度*9.8
    
17. FORC1=DENS2*Lside**3*9.8                !落石的自重
    
18. 
    
19. !地面岩土体--碎石土
    
20. MP,EX,1,1.0e9
    
21. MP,NUXY,1,0.24
    
22. MP,DENS,1,2100
    
23. TB,BISO
    
24. TBDATA,1,1.2e6,4.0e5              !屈服应力,越大回弹越高;切线模量,越大回弹越高
    
25. 
    
26. !落石
    
27. MP,EX,2,5.0e10
    
28. MP,NUXY,2,0.24
    
29. MP,DENS,2,DENS2
    
30. 
    
31. k,1,0.0,LY,-LZ
    
32. k,2,0.0,0,-LZ
    
33. k,3,2*LX,LY,-LZ
    
34. k,4,2*LX,0,-LZ
    
35. K,5,0.0,0.0,LZ
    
36. A,1,2,4,3
    
37. 
    
38. LESIZE,all,Lside/2,,,, , , ,1
    
39. amesh,all
    
40. l,2,5
    
41. TYPE,1   
    
42. EXTOPT,ESIZE,EL,1,  
    
43. EXTOPT,ACLEAR,1
    
44. EXTOPT,ATTR,0,0,0   
    
45. MAT,1   
    
46. REAL,1  
    
47. ESYS,0  
    
48. VDRAG,1,,,,,,5
    
49. cm,nsoil,node
    
50. cm,esoil,elem
    
51. ldele,5,,,1
    
52. allsel,all
    
53. !建立落石几何模型
    
54. block,LX-Lside/2,LX+Lside/2,LY,LY+Lside,-Lside/2,Lside/2
    
55. vsel,,,,2
    
56. vatt,2,2,2
    
57. aslv,s
    
58. lsla,s
    
59. LESIZE,ALL,,,2,1
    
60. vmesh,2
    
61. eslv,s
    
62. nsle,s
    
63. cm,nrock,node
    
64. cm,erock,elem
    
65. EDMP,RIGID,2,6,7      !用 EDMP命令定义刚性体,6指约束住x和z方向，7指约束住3个旋转自由度
    
66. aLLsel,aLL
    
67. *DIM,TIME,ARRAY,2,1,1,,,
    
68. *DIM,FORC,ARRAY,2,1,1,,,
    
69. *SET,FORC(1,1,1),-FORC1
    
70. *SET,FORC(2,1,1),-FORC1
    
71. *SET,TIME(2,1,1),TT
    
72. 
    
73. asel,,,,3
    
74. nsla,s,1
    
75. D,aLL,,,,,,ALL,,,,,
    
76. aLLsel,aLL
    
77. asel,,,,1,6,5
    
78. nsla,s,1
    
79. D,aLL,UZ,,,,,,,,,,
    
80. aLLsel,aLL
    
81. asel,,,,2,4,2
    
82. nsla,s,1
    
83. D,aLL,UX,,,,,,,,,,
    
84. aLLsel,aLL
    
85. EPLOT
    
86. EDPART,CREATE,1,esoil
    
87. EDPART,CREATE,2,erock
    
88. EDCGEN,STS,2,1,0.00,0.00,0.0,0,0.00,,,,,0,,0,0   !定义接触参数,粘滞阻尼系数VC=sy/sqrt(3),VDC  
    
89. !EDCONTACT,0.1,0,2,2,1,1,1,4,0                   !接触面控制
    
90. EDHGLS,0.10
    
91. EDCSC,ON
    
92. EDLOAD,ADD,RBFY,,2,TIME,FORC,0,,,,,
    
93. EDPVEL,VELO,2,vector,-vector,,0,0,0    !定义初始速度可以采用edvel或edpvel,X方向的速度没加载成功，不知道原因？？？？？
    
94. FINISH
    
95. 
    
96. /SOLU
    
97. antype,trans
    
98. TIMINT,ON
    
99. TRNOPT,FULL
    
100. TIME,TT,
     
101. EDRST,DTT,
     
102. EDHTIME,DTT,
     
103. EDDUMP,1,
     
104. !/STATUS,SOLU
     
105. SOLVE
     
106. FINISH
     
107. /DSCALE,1,OFF   
     
108. /REPLOT
     
109. 
     
110. vsel,,,,2
     
111. eslv,s
     
112. nsle,s
     
113. *get,Nrock,node,0,num,min
     
114. allsel,all
     
115. !落石运动特性与时间曲线
     
116. /POST26
     
117. NSOL,2,Nrock,U,Y,UY                           !块石Y方向位移的时程曲线
     
118. NSOL,3,Nrock,V,Y,VY                           !块石Y方向速度的时程曲线
     
119. NSOL,4,Nrock,A,Y,AY                           !块石Y方向加速度的时程曲线
     
120. NSOL,5,Nrock,U,X,UX                           !块石X方向位移的时程曲线
     
121. NSOL,6,Nrock,V,X,VX                           !块石X方向速度的时程曲线
     
122. NSOL,7,Nrock,A,X,AX                           !块石X方向加速度的时程曲线
     
123. STORE,MERGE  
     
124. XVAR,1                                        !定义X轴变量，1代表时间，n代表所定义的变量
     
125. PLVAR,2,5
     
126. PLVAR,3,6
     
127. PLVAR,4,7
     
128. PRVAR,2,3,4

复制代码

huizi42771

学习了！！！！！！！