连续钢构桥三维仿真分析【转】

Chelsea

1. 
   
2. 
   
3. /filname,rigid bridge,1
   
4. keyw,pr_struc,1
   
5. /prep7
   
6. !定义单元类型
   
7. et,1,shell63
   
8. et,2,link8
   
9. !定义实常数
   
10. r,1,0.25            !顶层与底层
    
11. r,2,0.6             !腹板
    
12. r,3,1               !v形支撑的厚度
    
13. r,4,0.02,-0.005
    
14. !定义材料属性
    
15. mp,ex,1,3.25e10   !混凝土
    
16. mp,prxy,1,0.2
    
17. mp,dens,1,2700
    
18. mp,ex,2,207e9     !钢筋
    
19. mp,prxy,2,0.3
    
20. mp,dens,2,7800
    
21. 
    
22. 
    
23. !建模
    
24. /view,1,1,1,1    !以正等侧方式观看
    
25. k,1              !建立主桥边界面上的关键点
    
26. k,2,4
    
27. k,3,10
    
28. k,4,16
    
29. k,5,20
    
30. k,6,24
    
31. k,7,30
    
32. k,8,36
    
33. k,9,40
    
34. k,10,4,-1.8
    
35. k,11,10,-1.8
    
36. k,12,16,-1.8
    
37. k,13,24,-1.8
    
38. k,14,30,-1.8
    
39. k,15,36,-1.8
    
40. kgen,2,1,15,1,,,-16      !建立梁高为1.8m处的关键点
    
41. kgen,2,1,15,1,,,-85
    
42. kgen,2,1,15,1,,,-154
    
43. kgen,2,1,15,1,,,-170
    
44. k,101,,,-36              !建立梁高为3m处的关键点
    
45. k,102,4,,-36
    
46. k,103,10,,-36
    
47. k,104,16,,-36
    
48. k,105,20,,-36
    
49. k,106,24,,-36
    
50. k,107,30,,-36
    
51. k,108,36,,-36
    
52. k,109,40,,-36
    
53. kgen,2,102,104,1,,-3,,8
    
54. kgen,2,106,108,1,,-3,,7
    
55. kgen,2,101,115,1,,,-12,15
    
56. kgen,2,116,130,1,,,-12,15
    
57. kgen,2,131,145,1,,,-50,15
    
58. kgen,3,146,160,1,,,-12,15
    
59. !建立箱形梁的顶层并划分网格
    
60. *do,i,1,8,1
    
61. a,i,i+15,i+16,i+1
    
62. *enddo
    
63. *do,i,16,23,1
    
64. a,i,i+85,i+86,i+1
    
65. *enddo
    
66. *do,i,101,108,1
    
67. a,i,i+15,i+16,i+1
    
68. *enddo
    
69. *do,i,116,123,1
    
70. a,i,i+15,i+16,i+1
    
71. *enddo
    
72. *do,i,31,38,1
    
73. a,i,i+100,i+101,i+1
    
74. *enddo
    
75. *do,i,31,38,1
    
76. a,i,i+115,i+116,i+1
    
77. *enddo
    
78. *do,i,146,153,1
    
79. a,i,i+15,i+16,i+1
    
80. *enddo
    
81. *do,i,161,168,1
    
82. a,i,i+15,i+16,i+1
    
83. *enddo
    
84. *do,i,46,53,1
    
85. a,i,i+130,i+131,i+1
    
86. *enddo
    
87. *do,i,46,53,1
    
88. a,i,i+15,i+16,i+1
    
89. *enddo
    
90. lesize,all,4            !设置单元长度为4
    
91. aatt,1,1,1
    
92. amesh,all
    
93. 
    
94. !建立箱形梁的底层并划分网格
    
95. asel,u,,,all
    
96. kplot
    
97. *do,i,10,11,1
    
98. a,i,i+15,i+16,i+1
    
99. *enddo
    
100. *do,i,13,14,1
     
101. a,i,i+15,i+16,i+1
     
102. *enddo
     
103. *do,i,25,26,1
     
104. a,i,i+85,i+86,i+1
     
105. *enddo
     
106. *do,i,28,29,1
     
107. a,i,i+85,i+86,i+1
     
108. *enddo
     
109. *do,i,110,111,1
     
110. a,i,i+15,i+16,i+1
     
111. *enddo
     
112. *do,i,113,114,1
     
113. a,i,i+15,i+16,i+1
     
114. *enddo
     
115. *do,i,125,126,1
     
116. a,i,i+15,i+16,i+1
     
117. *enddo
     
118. *do,i,128,129,1
     
119. a,i,i+15,i+16,i+1
     
120. *enddo
     
121. *do,i,40,41,1
     
122. a,i,i+100,i+101,i+1
     
123. *enddo
     
124. *do,i,43,44,1
     
125. a,i,i+100,i+101,i+1
     
126. *enddo
     
127. *do,i,40,41,1
     
128. a,i,i+115,i+116,i+1
     
129. *enddo
     
130. *do,i,43,44,1
     
131. a,i,i+115,i+116,i+1
     
132. *enddo
     
133. *do,i,155,156,1
     
134. a,i,i+15,i+16,i+1
     
135. *enddo
     
136. *do,i,158,159,1
     
137. a,i,i+15,i+16,i+1
     
138. *enddo
     
139. *do,i,170,171,1
     
140. a,i,i+15,i+16,i+1
     
141. *enddo
     
142. *do,i,173,174,1
     
143. a,i,i+15,i+16,i+1
     
144. *enddo
     
145. *do,i,55,56,1
     
146. a,i,i+130,i+131,i+1
     
147. *enddo
     
148. *do,i,58,59,1
     
149. a,i,i+130,i+131,i+1
     
150. *enddo
     
151. *do,i,55,56,1
     
152. a,i,i+15,i+16,i+1
     
153. *enddo
     
154. *do,i,58,59,1
     
155. a,i,i+15,i+16,i+1
     
156. *enddo
     
157. lesize,all,4
     
158. aatt,1,1,1
     
159. amesh,all
     
160. !建立腹板并划分网格
     
161. asel,u,,,all
     
162. *do,i,2,4,1
     
163. a,i,i+8,i+23,i+15
     
164. *enddo
     
165. 
     
166. *do,i,6,8,1
     
167. a,i,i+7,i+22,i+15
     
168. *enddo
     
169. *do,i,17,19,1
     
170. a,i,i+8,i+93,i+85
     
171. *enddo
     
172. *do,i,21,23,1
     
173. a,i,i+7,i+92,i+85
     
174. *enddo
     
175. *do,i,102,104,1
     
176. a,i,i+8,i+23,i+15
     
177. *enddo
     
178. *do,i,106,108,1
     
179. a,i,i+7,i+22,i+15
     
180. *enddo
     
181. *do,i,117,119,1
     
182. a,i,i+8,i+23,i+15
     
183. *enddo
     
184. *do,i,121,123,1
     
185. a,i,i+7,i+22,i+15
     
186. *enddo
     
187. *do,i,32,34,1
     
188. a,i,i+8,i+108,i+100
     
189. *enddo
     
190. *do,i,36,38,1
     
191. a,i,i+7,i+107,i+100
     
192. *enddo
     
193. *do,i,32,34,1
     
194. a,i,i+8,i+123,i+115
     
195. *enddo
     
196. *do,i,36,38,1
     
197. a,i,i+7,i+122,i+115
     
198. *enddo
     
199. *do,i,147,149,1
     
200. a,i,i+8,i+23,i+15
     
201. *enddo
     
202. *do,i,151,153,1
     
203. a,i,i+7,i+22,i+15
     
204. *enddo
     
205. *do,i,162,164,1
     
206. a,i,i+8,i+23,i+15
     
207. *enddo
     
208. *do,i,166,168,1
     
209. a,i,i+7,i+22,i+15
     
210. *enddo
     
211. *do,i,47,49,1
     
212. a,i,i+8,i+138,i+130
     
213. *enddo
     
214. *do,i,51,53,1
     
215. a,i,i+7,i+137,i+130
     
216. *enddo
     
217. *do,i,47,49,1
     
218. a,i,i+8,i+23,i+15
     
219. *enddo
     
220. *do,i,51,53,1
     
221. a,i,i+7,i+22,i+15
     
222. *enddo
     
223. lsel,s,loc,y,-0.9
     
224. lesize,all,,,1
     
225. lsel,s,loc,y,-1.5
     
226. lesize,all,,,2
     
227. esize,all,4
     
228. aatt,1,2,1
     
229. amesh,all
     
230. !预应力钢筋布置及网格划分
     
231. lsel,s,loc,x,3.9,4.1
     
232. lsel,a,loc,x,9,11
     
233. lsel,a,loc,x,15,17
     
234. lsel,a,loc,x,23,25
     
235. lsel,a,loc,x,29,31
     
236. lsel,a,loc,x,35,37
     
237. lsel,u,loc,y,-1.7,-0.1
     
238. lesize,all,4             !单元长度为4
     
239. latt,2,4,2               !单元为2，实常数为4，材料属性为2
     
240. lmesh,all
     
241. !建立v形支撑模型并划分网格
     
242. asel,u,loc,y,-3.1,0.1
     
243. k,1001,4,-15,-48
     
244. k,1002,16,-15,-48
     
245. k,1003,24,-15,-48
     
246. k,1004,36,-15,-48
     
247. k,2001,4,-15,-122
     
248. k,2002,16,-15,-122
     
249. k,2003,24,-15,-122
     
250. k,2004,36,-15,-122
     
251. a,110,112,1002,1001
     
252. a,112,113,1003,1002
     
253. a,113,115,1004,1003
     
254. a,140,142,1002,1001
     
255. a,142,143,1003,1002
     
256. a,143,145,1004,1003
     
257. a,155,157,2002,2001
     
258. a,157,158,2003,2002
     
259. a,158,160,2004,2003
     
260. a,185,187,2002,2001
     
261. a,187,188,2003,2002
     
262. a,188,190,2004,2003
     
263. lesize,all,4
     
264. aatt,1,3,1
     
265. amesh,all
     
266. nummrg,all
     
267. numcmp,all
     
268. allsel
     
269. !加载及求解
     
270. !施加位移约束
     
271. nsel,s,loc,z,-0.1,0.1
     
272. nsel,a,loc,z,-171,-169
     
273. nsel,r,loc,y,-1.9,-1.7
     
274. d,all,uy
     
275. allsel
     
276. nsel,s,loc,y,-16,-14
     
277. d,all,all
     
278. allsel
     
279. !预应力静力求解
     
280. /solu
     
281. acel,,9.8
     
282. pstres,on
     
283. allsel
     
284. solve
     
285. finish
     
286. !预应力模态求解
     
287. /solu
     
288. antype,2
     
289. upcoord,1,on
     
290. pstres,on
     
291. modopt,subsp,10
     
292. mxpand,10,,,1
     
293. allsel
     
294. solve
     
295. finish
     

复制代码

hancunliang

很给力  真的