ANSYS实例命令流程序_接触分析例题_岩土方面

Chelsea

例子1是为了获得初始应力，计算可以完成。命令流为：

1. !以便获得初始应力，收敛！！考虑土与岩石之间的摩擦，点面接触
   
2. fini !考虑土体的相互作用。自由划分网格，局部细化,不收敛！！！！！！！！！！！
   
3. /cle
   
4. /filname,diji11
   
5. /title,地基应力、应变的分布规律
   
6. /units,si
   
7. /PREP7
   
8. !ET,1,PLANE2
   
9. ET,1,PLANE42
   
10. KEYOPT,1,3,2
    
11. KEYOPT,1,5,0
    
12. KEYOPT,1,6,0
    
13. MP,EX,1,8e9 !岩石性质
    
14. MP,PRXY,1,0.29
    
15. MP,DENS,1,2600
    
16. 
    
17. MP,EX,2,5.0e7 !风化槽土体性质?
    
18. MP,PRXY,2,0.25
    
19. TB,dp,2
    
20. TBDATA,1,15e3,20,0
    
21. MP,DENS,2,2000
    
22. 
    
23. et,2,contac48,,,1
    
24. mp,mu,3,2.5
    
25. r,1,5.0e7,5.0e6,0.01,1.0,0.001,0.001
    
26. r,2,5.0e7,5.0e6,0.01,1.0,0.001,0.001
    
27. K,1,0,0,0,
    
28. K,2,2.5,-12,0,
    
29. K,3,6.5,-12,0,
    
30. K,4,13.5,-26.5,0,
    
31. K,5,20,-30,0,
    
32. K,6,24.5,-27,0,
    
33. K,7,31.5,-12,0
    
34. K,8,35.5,-12,0,
    
35. K,9,38,0,0,
    
36. k,10,-40,0,0,
    
37. k,11,-40,-60,0,
    
38. k,12,80,-60,0,
    
39. k,13,80,0,0,
    
40. K,14,6.5,-12,0,
    
41. K,15,13.5,-26.5,0,
    
42. K,16,20,-30,0,
    
43. K,17,24.5,-27,0,
    
44. K,18,31.5,-12,0
    
45. 
    
46. *do,i,1,8
    
47. lstr,i,i+1
    
48. *enddo
    
49. LSTR,1,10
    
50. 
    
51. *do,i,10,12
    
52. lstr,i,i+1
    
53. *enddo
    
54. LSTR,13,9
    
55. 
    
56. *do,i,14,17
    
57. lstr,i,i+1
    
58. *enddo
    
59. LARC,14,18,16,14.57
    
60. 
    
61. FLST,2,13,4
    
62. *do,i,1,13
    
63. FITEM,2,i
    
64. *enddo
    
65. AL,P51X
    
66. 
    
67. FLST,2,5,4
    
68. *do,i,14,18
    
69. FITEM,2,i
    
70. *enddo
    
71. AL,P51X
    
72. 
    
73. asel,s,area,,1
    
74. aatt,1,,1
    
75. asel,s,area,,2
    
76. aatt,2,,1
    
77. esize,4
    
78. asel,all
    
79. amesh,all
    
80. 
    
81. lsel,s,line,,3,4,1
    
82. nsll,s,1
    
83. CM,_TARGET1,node
    
84. lsel,s,line,,14,15,1
    
85. nsll,s,1
    
86. CM,_CONTACT1,node
    
87. lsel,s,line,,5,6,1
    
88. nsll,s,1
    
89. CM,_TARGET2,node
    
90. lsel,s,line,,16,17,1
    
91. nsll,s,1
    
92. CM,_CONTACT2,node
    
93. TYPE,2
    
94. mat,3
    
95. real,1
    
96. gcgen,_CONTACT1,_TARGET1
    
97. real,2
    
98. gcgen,_CONTACT2,_TARGET2
    
99. fini
    
100. FINISH
     
101. 
     
102. /SOLU
     
103. antype,static
     
104. nlgeom,on
     
105. nropt,full,,off
     
106. neqit,40
     
107. pred,on
     
108. autots,on
     
109. time,1
     
110. nsubst,100,1000,10
     
111. alls
     
112. LSEL,s,line,,10,12,2
     
113. CM,_Y,LINE
     
114. CMSEL,,_Y
     
115. DL,_Y, ,ux,0,0
     
116. LSEL,s,line,,11
     
117. CM,_Y,LINE
     
118. CMSEL,,_Y
     
119. DL,_Y, ,uy,0,0
     
120. alls
     
121. acel,,9.8
     
122. iswrite,on
     
123. lnsrch,on
     
124. outres,all,all
     
125. solve

复制代码

Chelsea

例2：

1. !读入初始应力，收敛！！考虑土和混凝土与岩石之间的摩擦，点面接触
   
2. fini !考虑土体和混凝土之间的相互作用。不收敛！！！！！！！！！！！
   
3. /cle
   
4. /filname,diji
   
5. /title,地基应力、应变的分布规律
   
6. /units,si
   
7. /PREP7
   
8. ET,1,PLANE42
   
9. KEYOPT,1,3,2
   
10. KEYOPT,1,5,0
    
11. KEYOPT,1,6,0
    
12. MP,EX,1,8e9 !岩石性质
    
13. MP,PRXY,1,0.29
    
14. MP,DENS,1,2600
    
15. 
    
16. MP,EX,2,5.0e7 !风化槽土体性质?
    
17. MP,PRXY,2,0.25
    
18. TB,dp,2
    
19. TBDATA,1,15e3,20,0
    
20. MP,DENS,2,2000
    
21. 
    
22. MP,EX,3,2.55e10 !混凝土性质
    
23. MP,PRXY,3,0.167
    
24. MP,DENS,3,2400
    
25. 
    
26. et,2,contac48,,,1
    
27. mp,mu,4,2.5
    
28. mp,mu,5,4.5
    
29. mp,mu,6,3.0
    
30. r,1,5.0e7,5.0e6,0.01,1.0,0.001,0.001
    
31. r,2,5.0e7,5.0e6,0.01,1.0,0.001,0.001
    
32. r,3,8.0e9,8.0e8,0.01,1.0,0.001,0.001
    
33. r,4,8.0e9,8.0e8,0.01,1.0,0.001,0.001
    
34. r,5,5.0e7,5.0e6,0.01,1.0,0.001,0.001
    
35. 
    
36. K,1,0,0,0,
    
37. K,2,2.5,-12,0,
    
38. K,3,6.5,-12,0,
    
39. K,4,13.5,-26.5,0,
    
40. K,5,20,-30,0,
    
41. K,6,24.5,-27,0,
    
42. K,7,31.5,-12,0
    
43. K,8,35.5,-12,0,
    
44. K,9,38,0,0,
    
45. k,10,-40,0,0,
    
46. k,11,-40,-60,0,
    
47. k,12,80,-60,0,
    
48. k,13,80,0,0,
    
49. K,14,6.5,-12,0,
    
50. K,15,13.5,-26.5,0,
    
51. K,16,20,-30,0,
    
52. K,17,24.5,-27,0,
    
53. K,18,31.5,-12,0
    
54. K,19,0,0,0,
    
55. K,20,2.5,-12,0,
    
56. K,21,6.5,-12,0,
    
57. K,22,31.5,-12,0
    
58. K,23,35.5,-12,0,
    
59. K,24,38,0,0,
    
60. K,25,20,0,0,
    
61. *do,i,1,8
    
62. lstr,i,i+1
    
63. *enddo
    
64. LSTR,1,10
    
65. 
    
66. *do,i,10,12
    
67. lstr,i,i+1
    
68. *enddo
    
69. LSTR,13,9
    
70. 
    
71. *do,i,14,17
    
72. lstr,i,i+1
    
73. *enddo
    
74. LARC,14,18,16,14.57
    
75. 
    
76. *do,i,19,20
    
77. lstr,i,i+1
    
78. *enddo
    
79. LARC,21,22,16,14.57
    
80. 
    
81. *do,i,22,24
    
82. lstr,i,i+1
    
83. *enddo
    
84. LSTR,25,19
    
85. 
    
86. FLST,2,13,4
    
87. *do,i,1,13
    
88. FITEM,2,i
    
89. *enddo
    
90. AL,P51X
    
91. 
    
92. FLST,2,5,4
    
93. *do,i,14,18
    
94. FITEM,2,i
    
95. *enddo
    
96. AL,P51X
    
97. 
    
98. FLST,2,7,4
    
99. *do,i,19,25
    
100. FITEM,2,i
     
101. *enddo
     
102. AL,P51X
     
103. 
     
104. asel,s,area,,1
     
105. aatt,1,,1
     
106. asel,s,area,,2
     
107. aatt,2,,1
     
108. esize,4
     
109. asel,s,area,,1,2
     
110. amesh,all
     
111. 
     
112. 
     
113. lsel,s,line,,3,4,1
     
114. nsll,s,1
     
115. CM,_TARGET1,node
     
116. lsel,s,line,,14,15,1
     
117. nsll,s,1
     
118. CM,_CONTACT1,node
     
119. lsel,s,line,,5,6,1
     
120. nsll,s,1
     
121. CM,_TARGET2,node
     
122. lsel,s,line,,16,17,1
     
123. nsll,s,1
     
124. CM,_CONTACT2,node
     
125. TYPE,2
     
126. mat,4
     
127. real,1
     
128. gcgen,_CONTACT1,_TARGET1
     
129. real,2
     
130. gcgen,_CONTACT2,_TARGET2
     
131. 
     
132. asel,s,area,,3
     
133. aatt,3,,1
     
134. esize,4
     
135. amesh,all
     
136. 
     
137. lsel,s,line,,1,8,7
     
138. nsll,s,1
     
139. CM,_TARGET3,node
     
140. lsel,s,line,,19,23,4
     
141. nsll,s,1
     
142. CM,_CONTACT3,node
     
143. lsel,s,line,,2,7,5
     
144. nsll,s,1
     
145. CM,_TARGET4,node
     
146. lsel,s,line,,20,22,2
     
147. nsll,s,1
     
148. CM,_CONTACT4,node
     
149. TYPE,2
     
150. mat,5
     
151. real,3
     
152. gcgen,_CONTACT3,_TARGET3
     
153. real,4
     
154. gcgen,_CONTACT4,_TARGET4
     
155. 
     
156. 
     
157. lsel,s,line,,18
     
158. nsll,s,1
     
159. CM,_TARGET5,node
     
160. lsel,s,line,,21
     
161. nsll,s,1
     
162. CM,_CONTACT5,node
     
163. TYPE,2
     
164. mat,6
     
165. real,5
     
166. gcgen,_CONTACT5,_TARGET5
     
167. fini
     
168. 
     
169. /SOLU
     
170. antype,static
     
171. nlgeom,on
     
172. nropt,full,,off
     
173. neqit,40
     
174. pred,on
     
175. autots,on
     
176. time,1
     
177. esel,s,,,1,451
     
178. nsle,s
     
179. isfile,read,diji11,ist,,1
     
180. nsubst,100,1000,10
     
181. LSEL,s,line,,10,12,2
     
182. CM,_Y,LINE
     
183. CMSEL,,_Y
     
184. DL,_Y, ,ux,0,0
     
185. LSEL,s,line,,11
     
186. CM,_Y,LINE
     
187. CMSEL,,_Y
     
188. DL,_Y, ,uy,0,0
     
189. acel,,9.8
     
190. alls
     
191. LSEL,s,line,,25
     
192. CM,_Y1,LINE
     
193. CMSEL,,_Y1
     
194. SFL,_Y1,PRES,4.5e5,6e5,
     
195. LSEL,s,line,,24
     
196. CM,_Y1,LINE
     
197. CMSEL,,_Y1
     
198. SFL,_Y1,PRES,4.5e5,4.5e5
     
199. alls
     
200. iswrite,on
     
201. lnsrch,on
     
202. outres,all,all
     
203. solve

复制代码

李子

支持一下，看了你好多帖子，希望以后可以看到你更多的帖子

笑话大哥

太感谢了，呵呵

sjwluobo

回复 2 # Chelsea 的帖子

楼主你好，我在做一个柔柔接触的时候，分析结果又很大的刚体位移，不知如何解决，能请您帮助看看吧
   
我固定筒中的内核的一端，然后在内核的另一端施加位移荷载。几个荷载步运行后发现DMX很大，可是我是假的位移荷载只有10个单位，查阅资料说是刚体位移，请问这个问题如何解决啊。
不知道是不是我的命令流有问题啊？