很经典的钢管混凝土拱桥ANSYS命令流

Chelsea

1. fini
   
2. /clear
   
3. /title,concrete-filled tube arch bridge
   
4. /prep7
   
5. 
   
6. !!!截面1
   
7. et,1,82
   
8. cyl4,0.6,0.6,0.6,,0.586
   
9. cyl4,0.6,2.4,0.6,,0.586
   
10. rectng,0.35,0.364,0.8,2.2
    
11. rectng,0.85,0.836,0.8,2.2
    
12. allsel
    
13. aptn,all
    
14. adele,5,8,1
    
15. allsel
    
16. 
    
17. aadd,all
    
18. smrtsize,5
    
19. amesh,all
    
20. 
    
21. secwrite,gg,sect,,1
    
22. sectype,1,beam,mesh
    
23. secoffset,cent,,,
    
24. secread,'gg','sect',,mesh
    
25. asel,all
    
26. aclear,all
    
27. adele,all,,,1
    
28. 
    
29. !!!截面2
    
30. cyl4,0.6,0.6,0.586
    
31. cyl4,0.6,2.4,0.586
    
32. rectng,0.364,0.836,0.8,2.2
    
33. allsel
    
34. aadd,all
    
35. 
    
36. smrtsize,5
    
37. amesh,all
    
38. 
    
39. secwrite,hnt,sect,,1
    
40. sectype,2,beam,mesh
    
41. secoffset,cent,,,
    
42. secread,'hnt','sect',,mesh
    
43. asel,all
    
44. aclear,all
    
45. adele,all,,,1
    
46. 
    
47. !!!!!!建立横梁标准截面a-a
    
48. k,1,-0.5
    
49. k,2,-0.75,0.25
    
50. k,3,-0.75,0.5
    
51. k,4,-0.55,0.5
    
52. k,5,-0.55,0.85
    
53. k,6,0.55,0.85
    
54. k,7,0.55,0.5
    
55. k,8,0.75,0.5
    
56. k,9,0.75,0.25
    
57. k,10,0.5,0
    
58. a,1,2,3,4,5,6,7,8,9,10
    
59. 
    
60. smrtsize,5
    
61. amesh,all
    
62. 
    
63. secwrite,hl1,sect,,1
    
64. sectype,4,beam,mesh
    
65. secoffset,cent,,,
    
66. secread,'hl1','sect',,mesh
    
67. asel,all
    
68. aclear,all
    
69. adele,all,,,1
    
70. 
    
71. !!!!!!建立横梁标准截面b-b
    
72. 
    
73. k,1,-0.5
    
74. k,2,-0.5,1
    
75. k,3,-0.75,1.25
    
76. k,4,-0.75,1.5
    
77. k,5,0.75,1.5
    
78. k,6,0.75,1.25
    
79. k,7,0.5,1
    
80. k,8,0.5
    
81. a,1,2,3,4,5,6,7,8
    
82. 
    
83. smrtsize,5
    
84. amesh,all
    
85. 
    
86. secwrite,hl2,sect,,1
    
87. sectype,5,beam,mesh
    
88. secoffset,cent,,,
    
89. secread,'hl2','sect',,mesh
    
90. asel,all
    
91. aclear,all
    
92. adele,all,,,1
    
93. 
    
94. 
    
95. !!!!!!建立横梁标准截面c-c
    
96. 
    
97. k,1,-0.5
    
98. k,2,-0.5,1.04
    
99. k,3,-0.75,1.04+0.25
    
100. k,4,-0.75,1.04+0.5
     
101. k,5,-0.55,1.04+0.5
     
102. k,6,-0.55,1.04+0.5+0.35
     
103. k,7,0.55,1.04+0.5+0.35
     
104. k,8,0.55,1.04+0.5
     
105. k,9,0.75,1.04+0.5
     
106. k,10,0.75,1.04+0.25
     
107. k,11,0.5,1.04
     
108. k,12,0.5
     
109. k,13,-0.2,0.58
     
110. k,14,-0.2,1.04+0.5-0.3
     
111. k,15,0.2,1.04+0.5-0.3
     
112. k,16,0.2,0.58
     
113. a,1,2,3,4,5,6,7,8,9,10,11,12
     
114. a,13,14,15,16
     
115. asba,1,2
     
116. 
     
117. smrtsize,5
     
118. amesh,all
     
119. 
     
120. secwrite,hl3,sect,,1
     
121. sectype,6,beam,mesh
     
122. secoffset,cent,,,
     
123. secread,'hl3','sect',,mesh
     
124. asel,all
     
125. aclear,all
     
126. adele,all,,,1
     
127. 
     
128. !!!!!!建立横梁标准截面d-d
     
129. ldele,all
     
130. kdele,all
     
131. 
     
132. k,1,-0.5
     
133. k,2,-0.5,1.23
     
134. k,3,-0.75,1.23+0.25
     
135. k,4,-0.75,1.23+0.5
     
136. k,5,-0.55,1.23+0.5
     
137. k,6,-0.55,1.23+0.5+0.35
     
138. k,7,0.55,1.23+0.5+0.35
     
139. k,8,0.55,1.23+0.5
     
140. k,9,0.75,1.23+0.5
     
141. k,10,0.75,1.23+0.25
     
142. k,11,0.5,1.23
     
143. k,12,0.5
     
144. k,13,-0.2,0.38
     
145. k,14,-0.3,0.38+0.2
     
146. k,15,-0.3,1.23+0.5-0.1-0.2
     
147. k,16,-0.2,1.23+0.5-0.1
     
148. k,17,0.2,1.23+0.5-0.1
     
149. k,18,0.3,1.23+0.5-0.1-0.2
     
150. k,19,0.3,0.38+0.2
     
151. k,20,0.2,0.38
     
152. 
     
153. a,1,2,3,4,5,6,7,8,9,10,11,12
     
154. a,13,14,15,16,17,18,19,20
     
155. asba,1,2
     
156. 
     
157. smrtsize,5
     
158. amesh,all
     
159. 
     
160. secwrite,hl4,sect,,1
     
161. sectype,7,beam,mesh
     
162. secoffset,cent,,,
     
163. secread,'hl4','sect',,mesh
     
164. asel,all
     
165. aclear,all
     
166. adele,all,,,1
     
167. 
     
168. 
     
169. !!!!!!建立横梁标准截面e-e
     
170. 
     
171. k,1,-0.5
     
172. k,2,-0.5,1.23
     
173. k,3,-0.75,1.23+0.25
     
174. k,4,-0.75,1.23+0.5
     
175. k,5,-0.55,1.23+0.5
     
176. k,6,-0.55,1.23+0.5+0.35
     
177. k,7,0.55,1.23+0.5+0.35
     
178. k,8,0.55,1.23+0.5
     
179. k,9,0.75,1.23+0.5
     
180. k,10,0.75,1.23+0.25
     
181. k,11,0.5,1.23
     
182. k,12,0.5
     
183. 
     
184. a,1,2,3,4,5,6,7,8,9,10,11,12
     
185. 
     
186. smrtsize,5
     
187. amesh,all
     
188. 
     
189. secwrite,hl5,sect,,1
     
190. sectype,8,beam,mesh
     
191. secoffset,cent,,,
     
192. secread,'hl5','sect',,mesh
     
193. asel,all
     
194. aclear,all
     
195. adele,all,,,1
     
196. 
     
197. etdele,1         !!将辅助单元删除
     
198. 
     
199. !!!***********************
     
200. !定义所有材料特性
     
201. et,1,beam44          !!钢管特性
     
202. mp,ex,1,2.1e11
     
203. mp,dens,1,7800
     
204. mp,prxy,1,0.3
     
205. 
     
206. n,90000,0,0,30       !!参考点
     
207. 
     
208. et,2,beam44          !!钢管内50＃混凝土特性
     
209. mp,ex,2,3.5e10
     
210. mp,dens,2,2600
     
211. mp,prxy,2,0.1667
     
212. 
     
213. et,3,beam44          !!纵梁30＃混凝土钢管特性
     
214. mp,ex,3,3.0e10
     
215. mp,dens,3,2600
     
216. mp,prxy,3,0.1667
     
217. 
     
218. et,4,beam44          !!横梁30＃混凝土钢管特性
     
219. mp,ex,4,3.0e10
     
220. mp,dens,4,2600
     
221. mp,prxy,4,0.1667
     
222. 
     
223. et,5,beam44          !!风撑特性
     
224. mp,ex,5,2.1e11
     
225. mp,dens,5,7800
     
226. mp,prxy,5,0.3
     
227. 
     
228. et,6,link10          !!吊杆特性（钢绞线）
     
229. mp,ex,6,1.9e11
     
230. mp,dens,6,7800
     
231. mp,prxy,6,0.3
     
232. keyopt,6,3,0             !只拉吊杆
     
233. 
     
234. et,7,beam44          !!盖梁30＃混凝土特性
     
235. mp,ex,7,3.0e10
     
236. mp,dens,7,2600
     
237. mp,prxy,7,0.1667
     
238. 
     
239. et,8,beam44          !!墩柱40＃混凝土特性
     
240. mp,ex,8,3.3e10
     
241. mp,dens,8,2600
     
242. mp,prxy,8,0.1667
     
243. 
     
244. et,9,beam44          !!横梁30＃混凝土特性
     
245. mp,ex,9,3.0e10
     
246. mp,dens,9,2600
     
247. mp,prxy,9,0.1667
     
248. 
     
249. et,10,link10          !!系杆特性（钢绞线）
     
250. mp,ex,10,2.0e11
     
251. mp,dens,10,7800
     
252. mp,prxy,10,0.3
     
253. keyopt,10,3,0             !只拉系杆
     
254. 
     
255. et,11,beam44          !!承台桩基基础25＃混凝土特性
     
256. mp,ex,11,2.85e10
     
257. mp,dens,11,2500
     
258. mp,prxy,11,0.1667
     
259. 
     
260. !!!!!!!!##############################
     
261. !开始建立模型
     
262. !0号横梁（端横梁）
     
263. n,1,,2.025
     
264. n,2,6.25,2.025
     
265. n,3,6.25+5.1,2.025
     
266. n,4,6.25+5.1+0.2+1.9/2,2.025   !12.5m位置
     
267. n,5,6.25+5.1+2.3,2.025
     
268. n,6,6.25+5.1+2.3+2.23,2.025
     
269. nsym,x,10,2,6,1
     
270. type,4
     
271. r,1
     
272. real,1
     
273. mat,4
     
274. secnum,7
     
275. !单元1到10号为0号横梁
     
276. e,1,2
     
277. secnum,6
     
278. e,2,3
     
279. secnum,5
     
280. e,3,4
     
281. e,4,5
     
282. secnum,4
     
283. e,5,6
     
284. secnum,7
     
285. e,1,12
     
286. secnum,6
     
287. e,12,13
     
288. secnum,5
     
289. e,13,14
     
290. e,14,15
     
291. secnum,4
     
292. e,15,16
     
293. !!1号端横梁
     
294. egen,2,100,1,10,1,,,,,,0,4.7+1.25/2+1.5/2
     
295. 
     
296. !!其余端横梁2～17号
     
297. egen,17,100,11,20,1,,,,,,0,5.3+1.5
     
298. 
     
299. !!18号端横梁
     
300. egen,2,100,171,180,1,,,,,,0,4.7+1.25/2+1.5/2
     
301. 
     
302. cm,hl,elem
     
303. esel,none
     
304. 
     
305. !!********************************
     
306. !纵梁单元
     
307. type,3
     
308. r,2
     
309. real,2
     
310. mat,3
     
311. sectype,10,beam,rect,zl
     
312. secdata,0.4,1.2
     
313. secnum,10
     
314. 
     
315. !0号纵梁
     
316. e,1,101
     
317. e,2,102
     
318. e,3,103
     
319. e,5,105
     
320. e,12,112
     
321. e,13,113
     
322. e,15,115
     
323. 
     
324. !中间纵梁1～16号
     
325. *do,i,0,15,1
     
326. e,101+i*100,201+i*100
     
327. e,102+i*100,202+i*100
     
328. e,104+i*100,204+i*100
     
329. e,112+i*100,212+i*100
     
330. e,114+i*100,214+i*100
     
331. *enddo
     
332. 
     
333. !17号纵梁（端横梁）
     
334. e,1701,1801
     
335. e,1702,1802
     
336. e,1703,1803
     
337. e,1705,1805
     
338. e,1712,1812
     
339. e,1713,1813
     
340. e,1715,1815
     
341. cm,zl,elem
     
342. esel,none
     
343. 
     
344. !!!****************************
     
345. !定义拱轴线
     
346. *dim,y,array,16   !将纵桥向距离定义为数组
     
347. *dim,cs,array,16  !定义分值cs＝（1－2x/l）
     
348. *dim,z,array,16   !将拱高度定义为数组
     
349. 
     
350. !吊杆上的点
     
351. !首先定义吊杆的y向（顺桥向）坐标
     
352. y(1)=8.1
     
353. y(2)=14.9
     
354. y(3)=21.7
     
355. y(4)=28.5
     
356. y(5)=35.3
     
357. y(6)=42.1
     
358. y(7)=48.9
     
359. y(8)=55.7
     
360. y(9)=62.5
     
361. !风撑上的点
     
362. !首先定义风撑的y向（顺桥向）坐标
     
363. y(11)=21.7
     
364. y(12)=28.5
     
365. y(13)=42.1
     
366. y(14)=48.9
     
367. y(15)=59.1
     
368. 
     
369. !0号横梁的y向坐标
     
370. y(16)=2.025
     
371. 
     
372. fl=25              !矢高
     
373. ml=1.1
     
374. k=log(ml+sqrt(ml*ml-1))   !k=ch-1(m)
     
375. *do,i,1,16,1
     
376. cs(i)=1-y(i)*2/125
     
377. z(i)=25-fl/(ml-1)*(cosh(k*cs(i))-1)   !拱轴高度
     
378. *enddo
     
379. !上述定义的拱轴线可参照姚玲森《桥梁工程》p1299
     
380. 
     
381. *dim,ygd,array,23  !定义预工拱度（是在后续工况中反复调算的结果）
     
382.                    !初始运行中可将预拱度全部设为0
     
383. ygd(1)=0.4000e-02+1.29e-02
     
384. ygd(2)=0.10685e-01+2.22e-02
     
385. ygd(3)=0.19483e-01+3.03e-02
     
386. ygd(4)=0.29219e-01+3.69e-02
     
387. ygd(5)=0.39022e-01+4.23e-02
     
388. ygd(6)=0.48243e-01+4.65e-02
     
389. ygd(7)=0.56390e-01+4.95e-02
     
390. ygd(8)=0.63070e-01+5.15e-02
     
391. ygd(9)=0.68014e-01+5.21e-02
     
392. ygd(10)=0.71003e-01+5.15e-02
     
393. ygd(11)=0.71739e-01+4.95e-02
     
394. ygd(12)=0.69962e-01+4.65e-02
     
395. ygd(13)=0.65514e-01+4.23e-02
     
396. ygd(14)=0.58289e-01+3.69e-02
     
397. ygd(15)=0.48248e-01+3.04e-02
     
398. ygd(16)=0.35507e-01+2.24e-02
     
399. ygd(17)=0.20415e-01+1.23e-02
     
400. ygd(18)=0.65770e-01+0.40234e-01
     
401. ygd(19)=0.69769e-01+0.41215e-01
     
402. 
     
403. !生成各吊杆的竖向节点（z向）
     
404. ngen,2,10000,104,114,10,0,0,z(1)+ygd(1)
     
405. ngen,2,10000,204,214,10,0,0,z(2)+ygd(2)
     
406. ngen,2,10000,304,314,10,0,0,z(3)+ygd(3)
     
407. ngen,2,10000,404,414,10,0,0,z(4)+ygd(4)
     
408. ngen,2,10000,504,514,10,0,0,z(5)+ygd(5)
     
409. ngen,2,10000,604,614,10,0,0,z(6)+ygd(6)
     
410. ngen,2,10000,704,714,10,0,0,z(7)+ygd(7)
     
411. ngen,2,10000,804,814,10,0,0,z(8)+ygd(8)
     
412. ngen,2,10000,904,914,10,0,0,z(9)+ygd(9)
     
413. ngen,2,10000,1004,1014,10,0,0,z(8)+ygd(10)
     
414. ngen,2,10000,1104,1114,10,0,0,z(7)+ygd(11)
     
415. ngen,2,10000,1204,1214,10,0,0,z(6)+ygd(12)
     
416. ngen,2,10000,1304,1314,10,0,0,z(5)+ygd(13)
     
417. ngen,2,10000,1404,1414,10,0,0,z(4)+ygd(14)
     
418. ngen,2,10000,1504,1514,10,0,0,z(3)+ygd(15)
     
419. ngen,2,10000,1604,1614,10,0,0,z(2)+ygd(16)
     
420. ngen,2,10000,1704,1714,10,0,0,z(1)+ygd(17)
     
421. 
     
422. ngen,2,3000,104,114,10,0,0,0  
     
423. ngen,2,3000,204,214,10,0,0,0  
     
424. ngen,2,3000,304,314,10,0,0,0  
     
425. ngen,2,3000,404,414,10,0,0,0  
     
426. ngen,2,3000,504,514,10,0,0,0  
     
427. ngen,2,3000,604,614,10,0,0,0  
     
428. ngen,2,3000,704,714,10,0,0,0  
     
429. ngen,2,3000,804,814,10,0,0,0  
     
430. ngen,2,3000,904,914,10,0,0,0  
     
431. ngen,2,3000,1004,1014,10,0,0,0
     
432. ngen,2,3000,1104,1114,10,0,0,0
     
433. ngen,2,3000,1204,1214,10,0,0,0
     
434. ngen,2,3000,1304,1314,10,0,0,0
     
435. ngen,2,3000,1404,1414,10,0,0,0
     
436. ngen,2,3000,1504,1514,10,0,0,0
     
437. ngen,2,3000,1604,1614,10,0,0,0
     
438. ngen,2,3000,1704,1714,10,0,0,0
     
439. 
     
440. ngen,2,6000,104,114,10,0,0,z(1)+ygd(1)
     
441. ngen,2,6000,204,214,10,0,0,z(2)+ygd(2)
     
442. ngen,2,6000,304,314,10,0,0,z(3)+ygd(3)
     
443. ngen,2,6000,404,414,10,0,0,z(4)+ygd(4)
     
444. ngen,2,6000,504,514,10,0,0,z(5)+ygd(5)
     
445. ngen,2,6000,604,614,10,0,0,z(6)+ygd(6)
     
446. ngen,2,6000,704,714,10,0,0,z(7)+ygd(7)
     
447. ngen,2,6000,804,814,10,0,0,z(8)+ygd(8)
     
448. ngen,2,6000,904,914,10,0,0,z(9)+ygd(9)
     
449. ngen,2,6000,1004,1014,10,0,0,z(8)+ygd(10)
     
450. ngen,2,6000,1104,1114,10,0,0,z(7)+ygd(11)
     
451. ngen,2,6000,1204,1214,10,0,0,z(6)+ygd(12)
     
452. ngen,2,6000,1304,1314,10,0,0,z(5)+ygd(13)
     
453. ngen,2,6000,1404,1414,10,0,0,z(4)+ygd(14)
     
454. ngen,2,6000,1504,1514,10,0,0,z(3)+ygd(15)
     
455. ngen,2,6000,1604,1614,10,0,0,z(2)+ygd(16)
     
456. ngen,2,6000,1704,1714,10,0,0,z(1)+ygd(17)
     
457. 
     
458. !*do,i,1,9,1
     
459. !ngen,2,10000,i*100+4,i*100+14,10,0,0,z(i)+ygd(i)    !吊杆节点
     
460. !*enddo
     
461. !
     
462. !jj=1
     
463. !*do,i,10,17,1
     
464. !ngen,2,10000,i*100+4,i*100+14,10,0,0,z(i-2*jj)+ygd(i)    !吊杆节点
     
465. !jj=jj+1
     
466. !*enddo
     
467. 
     
468. !!为了对以后吊杆转动进行放松，重复吊杆节点
     
469. !*do,i,1,9,1
     
470. !ngen,2,3000,i*100+4,i*100+14,10,0,0,z(i)+ygd(i)   !吊杆下节点
     
471. !*enddo
     
472. !
     
473. !jj=1
     
474. !*do,i,10,17,1
     
475. !ngen,2,3000,i*100+4,i*100+14,10,0,0,z(i-2*jj)+ygd(i)    !吊杆节点
     
476. !jj=jj+1
     
477. !*enddo
     
478. !
     
479. !*do,i,1,10,1
     
480. !ngen,2,6000,i*100+4,i*100+14,10,0,0,z(i)+ygd(i)   !吊杆节点
     
481. !*enddo
     
482. !
     
483. !jj=1
     
484. !*do,i,10,17,1
     
485. !ngen,2,6000,i*100+4,i*100+14,10,0,0,z(i-2*jj)+ygd(i)    !吊杆节点
     
486. !jj=jj+1
     
487. !*enddo
     
488. 
     
489. !定义吊杆的参数（每根吊杆均采用109根直径7mm高强钢丝）
     
490. dgarea=109*3.141593*0.25*0.007*0.007
     
491. dgyl=0.3*1670*1e6      !吊杆除应力采用0.3倍抗拉强度（n/mm*mm）
     
492. 
     
493. dgl=dgyl*dgarea
     
494. edg=1.95e11            !吊杆弹模
     
495. 
     
496. 
     
497. *dim,L0,array,17       !定义吊杆的初始0应力长度数组
     
498. *dim,L,array,17       !定义吊杆的受力后长度数组
     
499. 
     
500. *dim,deltL,array,17       !定义吊杆的受力后长度增量数组
     
501. *dim,istrn,array,17       !定义吊杆的应力数组
     
502. *dim,fdg,array,17       !定义在吊杆上施加的力数组
     
503. 
     
504. *do,i,1,9,1
     
505. L0(i)=z(i)           !0应力时单元长度
     
506. L0(18-i)=z(i)        !吊杆沿中轴对称
     
507. *enddo
     
508. 
     
509. *do,i,1,17,1
     
510. fdg(i)=dgl        !各吊杆施加的预应力
     
511. *enddo
     
512. 
     
513. type,6
     
514. mat,6
     
515. *do,i,1,17,1
     
516. deltL(i)=fdg(i)*L0(i)/(edg*dgarea)  !吊杆受力后的增长量
     
517. L(i)=L0(i)+deltL(i)
     
518. istrn(i)=deltL(i)/L(i)              !吊杆受力后的应变
     
519. r,i+1000,dgarea,istrn(i)            !定义吊杆的实常数特性
     
520. *enddo
     
521. 
     
522. *do,i,1,17,1
     
523. real,i+1000                   !赋予吊杆的实常数特性
     
524. e,i*100+4+3000,i*100+4+6000    !生成吊杆单元
     
525. e,i*100+10+4+3000,i*100+10+4+6000
     
526. *enddo
     
527. cm,dg,elem
     
528. esel,none
     
529. allsel
     
530. 
     
531. !!*********************************
     
532. real,1018      !定义一个全为零的实常数便于利用截面特性时，不产生冲突
     
533. !!*********************************
     
534. !生成风撑上的节点
     
535. n,20001,0,21.7,z(11)
     
536. n,20002,0,28.5,z(12)
     
537. n,20003,0,42.1,z(13)
     
538. n,20004,0,48.9,z(14)
     
539. n,20005,0,59.1,z(15)
     
540. n,20006,0,65.9,z(15)
     
541. n,20007,0,76.1,z(14)
     
542. n,20008,0,82.9,z(13)
     
543. n,20009,0,96.5,z(12)
     
544. n,20010,0,103.3,z(11)
     
545. n,20011,12.5,21.7,z(11)
     
546. n,20012,12.5,28.5,z(12)
     
547. n,20013,12.5,42.1,z(13)
     
548. n,20014,12.5,48.9,z(14)
     
549. n,20015,12.5,59.1,z(15)+ygd(18)
     
550. n,20016,12.5,65.9,z(15)+ygd(19)
     
551. n,20017,12.5,76.1,z(14)
     
552. n,20018,12.5,82.9,z(13)
     
553. n,20019,12.5,96.5,z(12)
     
554. n,20020,12.5,103.3,z(11)
     
555. n,20021,-12.5,21.7,z(11)
     
556. n,20022,-12.5,28.5,z(12)
     
557. n,20023,-12.5,42.1,z(13)
     
558. n,20024,-12.5,48.9,z(14)
     
559. n,20025,-12.5,59.1,z(15)+ygd(18)
     
560. n,20026,-12.5,65.9,z(15)+ygd(19)
     
561. n,20027,-12.5,76.1,z(14)
     
562. n,20028,-12.5,82.9,z(13)
     
563. n,20029,-12.5,96.5,z(12)
     
564. n,20030,-12.5,103.3,z(11)
     
565. n,20000,0,62.5,z(15)       !中撑中点
     
566. 
     
567. !风撑1（倒k撑）
     
568. type,5
     
569. mat,5
     
570. sectype,500,beam,ctube,windc          !定义风撑形状为空心圆
     
571. secdata,0.486,0.50
     
572. secnum,500                            !调用截面
     
573. 
     
574. !生成风撑单元
     
575. !风撑1（倒k撑）
     
576. e,20002,10404,90000
     
577. e,20002,10414,90000
     
578. e,20002,10304,90000
     
579. e,20002,10314,90000
     
580. 
     
581. !风撑2（倒k撑）
     
582. e,20004,10704,90000
     
583. e,20004,10714,90000
     
584. e,20004,10604,90000
     
585. e,20004,10614,90000
     
586. 
     
587. !风撑3（中撑x撑）
     
588. e,20000,20015,90000
     
589. e,20000,20016,90000
     
590. e,20000,20025,90000
     
591. e,20000,20026,90000
     
592. 
     
593. !风撑4（倒k撑）
     
594. e,20007,11104,90000
     
595. e,20007,11114,90000
     
596. e,20007,11204,90000
     
597. e,20007,11214,90000
     
598. 
     
599. !风撑5（倒k撑）
     
600. e,20009,11404,90000
     
601. e,20009,11414,90000
     
602. e,20009,11504,90000
     
603. e,20009,11514,90000
     
604. 
     
605. esel,s,real,,1018
     
606. cm,fc,elem
     
607. esel,none
     
608. 
     
609. !!*******************************************
     
610. !连接拱轴线
     
611. !定义拱轴上的坐标点
     
612. n,30000,12.5                     !左拱轴起点
     
613. n,30001,12.5,125,0               !左拱轴端点
     
614. n,30002,12.5,2.025,z(16)         !左横梁起点
     
615. n,30003,12.5,125-2.025,z(16)     !左横梁端点
     
616. 
     
617. n,40000,-12.5                     !右拱轴起点
     
618. n,40001,-12.5,125,0               !右拱轴端点
     
619. n,40002,-12.5,2.025,z(16)         !右横梁起点
     
620. n,40003,-12.5,125-2.025,z(16)     !右横梁端点
     
621. 
     
622. *dim,zgzx,array,23                !定义拱轴线，左轴线
     
623. *dim,ygzx,array,23                !定义拱轴线，右轴线
     
624. zgzx(1)=30000                     !以下为拱轴线上的节点号
     
625. zgzx(2)=30002
     
626. *do,i,3,10,1
     
627. zgzx(i)=10004+(i-2)*100
     
628. *enddo
     
629. zgzx(11)=20015
     
630. zgzx(12)=10904
     
631. zgzx(13)=20016
     
632. *do,i,14,21,1
     
633. zgzx(i)=11004+(i-14)*100
     
634. *enddo
     
635. zgzx(22)=30003
     
636. zgzx(23)=30001
     
637. ygzx(1)=40000
     
638. ygzx(2)=40002
     
639. *do,i,3,10,1
     
640. ygzx(i)=10014+(i-2)*100
     
641. *enddo
     
642. ygzx(11)=20025
     
643. ygzx(12)=10914
     
644. ygzx(13)=20026
     
645. *do,i,14,21,1
     
646. ygzx(i)=11014+(i-14)*100
     
647. *enddo
     
648. ygzx(22)=40003
     
649. ygzx(23)=40001
     
650. n,90001,12.5,0,30         !钢管左轴线参考点
     
651. n,90002,-12.5,0,30        !钢管右轴线参考点
     
652. type,1
     
653. mat,1
     
654. r,2000
     
655. secnum,1
     
656. *do,i,1,22                 !生成拱轴单元（钢管）
     
657. e,zgzx(i),zgzx(i+1),90001
     
658. e,ygzx(i),ygzx(i+1),90002
     
659. *enddo
     
660. cm,gzgg,elem
     
661. esel,none
     
662. type,2
     
663. mat,2
     
664. secnum,2
     
665. *do,i,1,22
     
666. e,zgzx(i),zgzx(i+1),90001             !生成拱轴单元（混凝土）
     
667. e,ygzx(i),ygzx(i+1),90002
     
668. *enddo
     
669. cm,gzhnt,elem
     
670. esel,none  

复制代码

Chelsea

1. !!!*****************************************
   
2. !建立拱轴底座
   
3. type,9
   
4. r,6
   
5. real,6
   
6. mat,9
   
7. sectype,70,beam,rect,dizuo             !定义截面形状
   
8. secdata,1.4,4
   
9. secnum,70
   
10. e,30001,30003,90001
    
11. e,40001,40003,90002
    
12. e,30000,30002,90001
    
13. e,40000,40002,90002
    
14. cm,gzdz,elem
    
15. esel,none
    
16. !!!*****************************************
    
17. !建立盖梁
    
18. ngen,2,2000,1,16,1,0,-2.025,0       !复制节点到拱轴线起点位置
    
19. ngen,2,100,1801,1816,1,0,2.025,0    !复制节点到拱轴线终点位置
    
20. type,7
    
21. mat,7
    
22. r,7
    
23. real,7
    
24. sectype,60,beam,rect,gail
    
25. secdata,4,1.6
    
26. secnum,60
    
27. n,90005,0,0,20
    
28. n,90006,0,125,20
    
29. e,1901,1902,90006         !盖梁单元125米处
    
30. e,1902,1903,90006
    
31. e,1903,1904,90006
    
32. e,1904,1905,90006
    
33. e,1905,1906,90006
    
34. e,1901,1912,90006
    
35. e,1912,1913,90006
    
36. e,1913,1914,90006
    
37. e,1914,1915,90006
    
38. e,1915,1916,90006
    
39. e,2001,2002,90005
    
40. e,2002,2003,90005
    
41. e,2003,2004,90005
    
42. e,2004,2005,90005
    
43. e,2005,2006,90005
    
44. e,2001,2012,90005
    
45. e,2012,2013,90005
    
46. e,2013,2014,90005
    
47. e,2014,2015,90005
    
48. e,2015,2016,90005
    
49. cm,gailiang,elem
    
50. esel,none
    
51. !*****************************
    
52. !建立墩柱
    
53. n,90003,14.5              ! 参考点
    
54. n,90004,-14.5
    
55. 
    
56. type,8
    
57. mat,8
    
58. r,8
    
59. real,8
    
60. sectype,50,beam,hrec,dun              !定义截面形状
    
61. secdata,3,4,0.8,0.8,0.8,0.8
    
62. ngen,2,10000,2004,2014,10,0,0,-10
    
63. ngen,2,10000,1904,1914,10,0,0,-10      !墩柱节点，距离桥面-10m
    
64. 
    
65. secnum,50
    
66. e,2004,12004,90003
    
67. e,2014,12014,90004
    
68. e,1904,11904,90003
    
69. e,1914,11914,90004
    
70. cm,dunzhu,elem
    
71. esel,none
    
72. !!!*****************************************
    
73. !承台单元
    
74. type,11
    
75. mat,11
    
76. r,1
    
77. sectype,40,beam,rect,chengtai
    
78. secdata,12.7,9
    
79. r,11
    
80. real,11
    
81. secnum,40
    
82. secdata,3,4,0.8,0.8,0.8,0.8
    
83. ngen,2,20000,2004,2014,10,0,0,-11.5
    
84. ngen,2,20000,1904,1914,10,0,0,-11.5
    
85. ngen,2,30000,2004,2014,10,0,0,-13
    
86. ngen,2,30000,1904,1914,10,0,0,-13
    
87. 
    
88. e,12004,22004,90003         !生成承台单元
    
89. e,11904,21904,90003
    
90. e,12014,22014,90004
    
91. e,11914,21914,90004
    
92. 
    
93. e,22004,32004,90003        
    
94. e,21904,31904,90003
    
95. e,22014,32014,90004
    
96. e,21914,31914,90004
    
97. 
    
98. cm,chtai,elem
    
99. esel,none
    
100. 
     
101. !!**********************************
     
102. !横系梁
     
103. type,9
     
104. mat,9
     
105. r,9
     
106. real,9
     
107. n,90007,12.5,0,10                !参考点
     
108. n,90008,12.5,0,10
     
109. sectype,90,beam,rect,hxl
     
110. secdata,3,3
     
111. secnum,90
     
112. e,22004,22014,90007
     
113. e,21904,21914,90008
     
114. cm,hxl,elem
     
115. esel,none
     
116. 
     
117. !!**********************************
     
118. !桩
     
119. n,50000,6.15+1.85,2.7,-13          !定义桩节点
     
120. n,50001,6.15+1.85,-2.7,-13
     
121. n,50002,6.15+1.85+4.5,2.7,-13
     
122. n,50003,6.15+1.85+4.5,-2.7,-13
     
123. n,50004,6.15+12.7-1.85,2.7,-13
     
124. n,50005,6.15+12.7-1.85,-2.7,-13
     
125. nsel,s,,,50000,50005,1
     
126. nsym,x,10,50000,50005,1
     
127. allsel
     
128. ngen,2,100,50000,50015,1,0,125,0          !首先生成平面内的24个节点
     
129. 
     
130. !每3m生成一层节点（共8层）
     
131. ngen,10,1000,50000,51000,1,0,0,-2.5
     
132. !生成桩单元
     
133. esel,none
     
134. type,11
     
135. mat,11
     
136. r,11
     
137. real,11
     
138. sectype,20,beam,csolid,pier
     
139. secdata,0.9
     
140. secnum,20
     
141. *do,i,0,8,1
     
142. *do,k,0,5,1
     
143. e,50000+k+i*1000,50000+k+(i+1)*1000,90003
     
144. e,50000+k+10+i*1000,50000+k+10+(i+1)*1000,90003
     
145. e,50000+k+100+i*1000,50000+k+100+(i+1)*1000,90003
     
146. e,50000+k+110+i*1000,50000+k+110+(i+1)*1000,90003
     
147. *enddo
     
148. *enddo
     
149. cm,zhuang,elem
     
150. esel,none
     
151. 
     
152. !!!***********************************
     
153. !在桩上施加弹簧单元
     
154. !第一层距离承台-3m   317-340
     
155. !计算弹簧的弹性系数   
     
156. d=1.8
     
157. kf=0.9
     
158. k0=1+1/d
     
159. b1=0.6       !和每排桩有关的系数，当为2排时为0.6
     
160. 
     
161. L1=3.6       !桩间净距
     
162. h1=3*(d+1)   !桩在地面或冲刷线下的计算深度
     
163. 
     
164. k=b1+(1-b1)/0.6*L1/h1
     
165. B=kf*k0*k*d
     
166. !第1层土粉质量粘土m=8厚度5
     
167. !第2层土粉粘土（细砂）m=15厚度5
     
168. !第3层土中砂m=20厚度5
     
169. !第4层土粘土m=25厚度5
     
170. !第5层土基岩c0=300厚度10
     
171. *dim,m,array,5
     
172. m(1)=8                     
     
173. m(2)=15
     
174. m(3)=20
     
175. m(4)=25
     
176. m(5)=4000                 !基岩系数
     
177. *dim,h,array,5
     
178. h(1)=5                    !第一层土的厚度                  
     
179. h(2)=5
     
180. h(3)=5
     
181. h(4)=5
     
182. h(5)=10
     
183. *dim,c,array,4
     
184. c(1)=m(1)*h(1)
     
185. c(2)=m(2)*(h(1)+h(2))
     
186. c(3)=m(3)*(h(1)+h(2)+h(3))
     
187. c(4)=m(4)*(h(1)+h(2)+h(3)+h(4))
     
188. 
     
189. *dim,e,array,4            !弹簧弹性系数
     
190. e(1)=0.5*c(1)*h(1)*B*1E6
     
191. *do,i,2,4,1
     
192. e(i)=0.5*(c(i-1)+c(i))*h(i)*B*1E6
     
193. *enddo
     
194. et,20,combin14                  !弹簧单元
     
195. r,20,e(1)                       !第一层土的弹性系数
     
196. r,30,e(2)                       !第二层土的弹性系数
     
197. r,40,e(3)                       !第三层土的弹性系数
     
198. r,50,e(4)                       !第四层土的弹性系数
     
199. r,60,m(5)*3.1415926*B*B/4*1E6   !基岩弹性系数
     
200. !第1层底部距离承台－5m
     
201. ngen,2,10000,51000,51115,1,1,0,0  !x方向
     
202. ngen,2,20000,51000,51115,1,1,0,0  !y方向
     
203. type,20
     
204. real,20
     
205. *do,i,0,5,1
     
206. e,51000+i,51000+i+10000
     
207. e,51000+i+10,51000+i+10+10000
     
208. e,51100+i,51100+i+10000
     
209. e,51100+i+10,51100+i+10+10000
     
210. e,51000+i,51000+i+20000
     
211. e,51000+i+10,51000+i+10+20000
     
212. e,51100+i,51100+i+20000
     
213. e,51100+i+10,51100+i+10+20000
     
214. *enddo
     
215. cm,spring1,elem
     
216. esel,none
     
217. 
     
218. 
     
219. !第2层底部距离承台－10m
     
220. ngen,2,10000,53000,53115,1,1,0,0  !x方向
     
221. ngen,2,20000,53000,53115,1,1,0,0  !y方向
     
222. type,20
     
223. real,30
     
224. *do,i,0,5,1
     
225. e,53000+i,53000+i+10000
     
226. e,53000+i+10,53000+i+10+10000
     
227. e,53100+i,53100+i+10000
     
228. e,53100+i+10,53100+i+10+10000
     
229. e,53000+i,53000+i+20000
     
230. e,53000+i+10,53000+i+10+20000
     
231. e,53100+i,53100+i+20000
     
232. e,53100+i+10,53100+i+10+20000
     
233. *enddo
     
234. cm,spring2,elem
     
235. esel,none
     
236. 
     
237. !第3层底部距离承台－15m
     
238. ngen,2,10000,55000,55115,1,1,0,0  !x方向
     
239. ngen,2,20000,55000,55115,1,1,0,0  !y方向
     
240. type,20
     
241. real,40
     
242. *do,i,0,5,1
     
243. e,55000+i,55000+i+10000
     
244. e,55000+i+10,55000+i+10+10000
     
245. e,55100+i,55100+i+10000
     
246. e,55100+i+10,55100+i+10+10000
     
247. e,55000+i,55000+i+20000
     
248. e,55000+i+10,55000+i+10+20000
     
249. e,55100+i,55100+i+20000
     
250. e,55100+i+10,55100+i+10+20000
     
251. *enddo
     
252. cm,spring3,elem
     
253. esel,none
     
254. 
     
255. 
     
256. !第4层底部距离承台－20m
     
257. ngen,2,10000,57000,57115,1,1,0,0  !x方向
     
258. ngen,2,20000,57000,57115,1,1,0,0  !y方向
     
259. type,20
     
260. real,50
     
261. *do,i,0,5,1
     
262. e,57000+i,57000+i+10000
     
263. e,57000+i+10,57000+i+10+10000
     
264. e,57100+i,57100+i+10000
     
265. e,57100+i+10,57100+i+10+10000
     
266. e,57000+i,57000+i+20000
     
267. e,57000+i+10,57000+i+10+20000
     
268. e,57100+i,57100+i+20000
     
269. e,57100+i+10,57100+i+10+20000
     
270. *enddo
     
271. cm,spring4,elem
     
272. esel,none
     
273. 
     
274. !第5层底部距离承台－22.5m
     
275. ngen,2,10000,58000,58115,1,1,0,0  !x方向
     
276. ngen,2,20000,58000,58115,1,1,0,0  !y方向
     
277. type,20
     
278. real,60
     
279. *do,i,0,5,1
     
280. e,58000+i,58000+i+10000
     
281. e,58000+i+10,58000+i+10+10000
     
282. e,58100+i,58100+i+10000
     
283. e,58100+i+10,58100+i+10+10000
     
284. e,58000+i,58000+i+20000
     
285. e,58000+i+10,58000+i+10+20000
     
286. e,58100+i,58100+i+20000
     
287. e,58100+i+10,58100+i+10+20000
     
288. *enddo
     
289. cm,spring5,elem
     
290. esel,none
     
291. 
     
292. !第6层底部距离承台－25m
     
293. ngen,2,10000,59000,59115,1,1,0,0  !x方向
     
294. ngen,2,10000,59000,59115,1,1,0,0  
     
295. ngen,2,20000,59000,59115,1,0,1,0  !y方向
     
296. ngen,2,20000,59000,59115,1,0,1,0
     
297. ngen,2,30000,59000,59115,1,0,0,-1  !z方向
     
298. ngen,2,30000,59000,59115,1,0,0,-1
     
299. type,20
     
300. real,60
     
301. *do,i,0,5,1
     
302. e,59000+i,59000+i+10000
     
303. e,59000+i+10,59000+i+10+10000
     
304. e,59100+i,59100+i+10000
     
305. e,59100+i+10,59100+i+10+10000
     
306. e,59000+i,59000+i+20000
     
307. e,59000+i+10,59000+i+10+20000
     
308. e,59100+i,59100+i+20000
     
309. e,59100+i+10,59100+i+10+20000
     
310. e,59000+i,59000+i+30000      
     
311. e,59000+i+10,59000+i+10+30000
     
312. e,59100+i,59100+i+30000
     
313. e,59100+i+10,59100+i+10+30000
     
314. *enddo
     
315. cm,spring6,elem
     
316. allsel
     
317. 
     
318. !!!*******************************************************
     
319. !剔除不需要进行对称复制的单元
     
320. nsel,s,loc,z,-50,-0.01
     
321. nsel,r,loc,y,-10,10
     
322. cm,zdnode,node
     
323. esln,s,,all
     
324. cm,zdelem,elem
     
325. allsel
     
326. nsel,u,,,zdnode
     
327. nsym,y,100000,all                      !节点对称
     
328. allsel
     
329. esel,u,,,zdelem
     
330. esym,y,100000,all
     
331. allsel
     
332. 
     
333. 
     
334. !!!!****************************************************
     
335. !定义系杆单元（全桥通长）
     
336. !定义系杆的参数（每根系杆都采用8束27×7×d5高强钢丝）
     
337. xgarea=8*27*7*3.1415926*0.25*0.005*0.005     !系杆集束面积
     
338. exg=1.95e11
     
339. *dim,fxg,array,2
     
340. *dim,Llast,array,2
     
341. *dim,Listrn,array,2
     
342. *dim,detL,array,2
     
343. xgyl=0.3*1860*1e6
     
344. xgl=xgyl*xgarea
     
345. fxg(1)=xgl              !施加的预应力
     
346. fxg(2)=xgl              !施加的预应力
     
347. 
     
348. Lzero=250               !0应力时单元长度
     
349. *do,i,1,2,1
     
350. detL(i)=Lzero*fxg(i)/(exg*xgarea)
     
351. Llast(i)=Lzero+detL(i)              !单元最终长度
     
352. Listrn(i)=detl(i)/Llast(i)
     
353. r,i+100,xgarea,Listrn(i)
     
354. *enddo
     
355. 
     
356. type,10
     
357. mat,10
     
358. real,101
     
359. e,30001,130001                  !系杆1  250m
     
360. real,102
     
361. e,40001,140001                  !系杆2  250m
     
362. cm,xigan,elem
     
363. esel,none
     
364. allsel
     
365. 
     
366. !!*******************************************
     
367. !将盖梁和拱轴交界点进行节点耦合，放松uy，和rotx
     
368. cp,1,ux,2004,30000
     
369. cp,1,uy,2004,30000
     
370. cp,next,uz,2004,30000
     
371. cp,next,roty,2004,30000
     
372. cp,next,rotz,2004,30000
     
373. cp,next,ux,1904,30001
     
374. cp,next,uy,1904,30001
     
375. cp,next,uz,1904,30001
     
376. cp,next,roty,1904,30001
     
377. cp,next,rotz,1904,30001
     
378. cp,next,ux,2014,40000
     
379. cp,next,uy,2014,40000
     
380. cp,next,uz,2014,40000
     
381. cp,next,roty,2014,40000
     
382. cp,next,rotz,2014,40000
     
383. cp,next,ux,1914,40001
     
384. cp,next,uy,1914,40001
     
385. cp,next,uz,1914,40001
     
386. cp,next,roty,1914,40001
     
387. cp,next,rotz,1914,40001
     
388. cp,next,ux,101914,140001
     
389. cp,next,uy,101914,140001
     
390. cp,next,uz,101914,140001
     
391. cp,next,roty,101914,140001
     
392. cp,next,rotz,101914,140001
     
393. cp,next,ux,101904,130001
     
394. cp,next,uy,101904,130001
     
395. cp,next,uz,101904,130001
     
396. cp,next,roty,101904,130001
     
397. cp,next,rotz,101904,130001
     
398. allsel
     
399. nsel,s,,,30000,130000,100000
     
400. nsel,a,,,40000,140000,100000
     
401. nummrg,node
     
402. allsel
     
403. 
     
404. !!!**********************************************
     
405. *do,i,1,3,1
     
406. cp,next,all,1800+i,1900+i               !盖梁和0号横梁耦合
     
407. *enddo
     
408. allsel
     
409. *do,i,5,6,1
     
410. cp,next,all,1800+i,1900+i              
     
411. *enddo
     
412. allsel
     
413. *do,i,1,2,1
     
414. cp,next,all,1800+i+11,1900+i+11
     
415. cp,next,all,i+11,2000+i+11,100000+i+11              
     
416. *enddo
     
417. allsel
     
418. *do,i,4,5,1
     
419. cp,next,all,1800+i+11,1900+i+11
     
420. cp,next,all,i+11,2000+i+11,100000+i+11            
     
421. *enddo
     
422. allsel
     
423. *do,i,1,3,1
     
424. cp,next,all,i,2000+i,100000+i              
     
425. *enddo
     
426. allsel
     
427. *do,i,5,6,1
     
428. cp,next,all,i,2000+i,100000+i              
     
429. *enddo
     
430. allsel
     
431. *do,i,1,3,1
     
432. cp,next,all,101800+i,101900+i              
     
433. *enddo
     
434. allsel
     
435. *do,i,5,6,1
     
436. cp,next,all,101800+i,101900+i              
     
437. *enddo
     
438. allsel
     
439. *do,i,1,2,1
     
440. cp,next,all,101800+i+11,101900+i+11              
     
441. *enddo
     
442. allsel
     
443. *do,i,4,5,1
     
444. cp,next,all,101800+i+11,101900+i+11              
     
445. *enddo
     
446. allsel
     
447. 
     
448. !!!*******************************************
     
449. !!桩节点和承台节点耦合
     
450. nsel,s,loc,z,-13
     
451. nsel,r,loc,y,-3,3
     
452. nsel,r,loc,x,8,20
     
453. cp,next,all,all
     
454. allsel
     
455. nsel,s,loc,z,-13
     
456. nsel,r,loc,y,-3,3
     
457. nsel,r,loc,x,-8,-20
     
458. cp,next,all,all
     
459. allsel
     
460. nsel,s,loc,z,-13
     
461. nsel,r,loc,y,122,128
     
462. nsel,r,loc,x,-8,-20
     
463. cp,next,all,all
     
464. allsel
     
465. nsel,s,loc,z,-13
     
466. nsel,r,loc,y,122,128
     
467. nsel,r,loc,x,8,20
     
468. cp,next,all,all
     
469. allsel
     
470. nsel,s,loc,z,-13
     
471. nsel,r,loc,y,-122,-128
     
472. nsel,r,loc,x,8,20
     
473. cp,next,all,all
     
474. allsel
     
475. nsel,s,loc,z,-13
     
476. nsel,r,loc,y,-122,-128
     
477. nsel,r,loc,x,-8,-20
     
478. cp,next,all,all
     
479. allsel
     
480. 
     
481. !!!*****************************************
     
482. !将吊杆转动放松
     
483. *do,i,1,17,1
     
484. cp,next,ux,i*100+4,i*100+4+3000
     
485. cp,next,ux,i*100+4+6000,i*100+4+10000
     
486. cp,next,uy,i*100+4,i*100+4+3000
     
487. cp,next,uy,i*100+4+6000,i*100+4+10000
     
488. cp,next,uz,i*100+4,i*100+4+3000
     
489. cp,next,uz,i*100+4+6000,i*100+4+10000
     
490. cp,next,ux,i*100+4+10,i*100+4+3000+10
     
491. cp,next,ux,i*100+4+6000+10,i*100+4+10000+10
     
492. cp,next,uy,i*100+4+10,i*100+4+3000+10      
     
493. cp,next,uy,i*100+4+6000+10,i*100+4+10000+10
     
494. cp,next,uz,i*100+4+10,i*100+4+3000+10      
     
495. cp,next,uz,i*100+4+6000+10,i*100+4+10000+10
     
496. cp,next,ux,i*100+4+100000,i*100+4+3000+100000      
     
497. cp,next,ux,i*100+4+6000+100000,i*100+4+10000+100000
     
498. cp,next,uy,i*100+4+100000,i*100+4+3000+100000      
     
499. cp,next,uy,i*100+4+6000+100000,i*100+4+10000+100000
     
500. cp,next,uz,i*100+4+100000,i*100+4+3000+100000      
     
501. cp,next,uz,i*100+4+6000+100000,i*100+4+10000+100000
     
502. cp,next,ux,i*100+4+100000+10,i*100+4+3000+100000+10      
     
503. cp,next,ux,i*100+4+6000+100000+10,i*100+4+10000+100000+10
     
504. cp,next,uy,i*100+4+100000+10,i*100+4+3000+100000+10      
     
505. cp,next,uy,i*100+4+6000+100000+10,i*100+4+10000+100000+10
     
506. cp,next,uz,i*100+4+100000+10,i*100+4+3000+100000+10      
     
507. cp,next,uz,i*100+4+6000+100000+10,i*100+4+10000+100000+10
     
508. *enddo
     
509. allsel
     
510. !!*****************************************
     
511. !将土弹簧的一端约束
     
512. nsel,s,loc,x,13.5
     
513. nsel,a,loc,x,9
     
514. nsel,a,loc,x,18
     
515. nsel,a,loc,x,-11.5
     
516. nsel,a,loc,x,-7
     
517. nsel,a,loc,x,-16
     
518. nsel,a,loc,y,3.7
     
519. nsel,a,loc,y,-1.7
     
520. nsel,a,loc,y,128.7
     
521. nsel,a,loc,y,123.3
     
522. nsel,a,loc,y,-128.7
     
523. nsel,a,loc,y,-123.3
     
524. nsel,a,loc,y,-128.7
     
525. nsel,a,loc,y,-123.3
     
526. nsel,r,loc,z,-40,-13
     
527. d,all,all
     
528. allsel
     
529. nsel,s,loc,z,-36.5
     
530. d,all,all
     
531. /color,pbak,off
     
532. /rep
     
533. allsel
     
534. nsel,s,loc,y,0
     
535. nsel,r,loc,z,0
     
536. nummrg,node,all
     
537. esel,none
     
538. allsel
     
539. fini
     
540. 
     
541. !!!******************************
     
542. !工况1 不包括桥面板时，求横梁和纵梁的总重量
     
543. /solu
     
544. allsel
     
545. mpdele,dens,all
     
546. mp,dens,3,2600
     
547. mp,dens,4,2600
     
548. acel,,,9.8
     
549. pstres,on
     
550. solve
     
551. save
     
552. !总重力：0.65064e＋08N
     
553. !*********************************************
     
554. !工况2 恒载下调节预拱度
     
555. !首先计算只有拱肋自重作用下拱的变位
     
556. /solu
     
557. allsel
     
558. mpdele,dens,3,11,1
     
559. mp,dens,1,7800
     
560. mp,dens,2,2600
     
561. *do,i,1,17,1
     
562. r,1000+i,dgarea,0
     
563. *enddo
     
564. r,101,xgarea,0
     
565. r,102,xgarea,0
     
566. acel,,,9.8
     
567. solve
     
568. 
     
569. !!********************************************
     
570. !工况3 一期恒载（包括桥面板）作用下的无预应力状态
     
571. /solu
     
572. allsel
     
573. *do,i,1,17,1
     
574. r,1000+i,dgarea,0
     
575. *enddo
     
576. r,101,xgarea,0
     
577. r,102,xgarea,0
     
578. qmb=(250/1.5)*32/6.8*5*1000*9.8    !桥面板的总重量38431372。5N
     
579. acel,,,9.8
     
580. mpdele,dens,3,4,1
     
581. mp,dens,1,7800
     
582. mp,dens,2,2600
     
583. mp,dens,3,4167.7        !加上桥面板后的换算质量密度
     
584. mp,dens,4,4167.7
     
585. mp,dens,5,7800          !风撑特性
     
586. mp,dens,6,7800          !吊杆特性
     
587. mp,dens,7,2600
     
588. mp,dens,8,2600
     
589. mp,dens,9,2600
     
590. mp,dens,10,7800         !系杆质量密度
     
591. mp,dens,11,2600
     
592. solve
     
593. !拱脚位移0.05m
     
594. !*******************************************
     
595. !工况4
     
596. /solu
     
597. allsel
     
598. asel,,,9.8
     
599. r,101,xgarea,Listrn(1)*1.3
     
600. r,102,xgarea,Listrn(2)*1.3        !采用单控，即先调节系杆
     
601. pstres,on
     
602. solve
     
603. !0.3倍应力极限时的拉力为16565931.9N
     
604. !拱脚的位移为-0.001487（合乎要求）
     
605. !****************************************
     
606. !工况5 一期恒载下计算调平
     
607. /solu
     
608. allsel
     
609. r,1001,dgarea,1.6787E-03+7.64E-05
     
610. r,1002,dgarea,1.7538E-03+8.48E-07
     
611. r,1003,dgarea,1.7457E-03+1.03E-06
     
612. r,1004,dgarea,1.7400E-03+2.71E-06
     
613. r,1005,dgarea,1.7292E-03+6.17E-06
     
614. r,1006,dgarea,1.7303E-03+6.37E-06
     
615. r,1007,dgarea,1.7267E-03+6.47E-06
     
616. r,1008,dgarea,1.7279E-03+7.69E-06
     
617. r,1009,dgarea,1.7365E-03+7.91E-06
     
618. r,10010,dgarea,1.7277E-03+7.81E-06
     
619. r,10011,dgarea,1.7263E-03+6.56E-06
     
620. r,10012,dgarea,1.7299E-03+6.24E-06
     
621. r,10013,dgarea,1.7290E-03+5.53E-06
     
622. r,10014,dgarea,1.7412E-03+2.03E-06
     
623. r,10015,dgarea,1.7504E-03+3.03E-06
     
624. r,10016,dgarea,1.7587E-03+1.28E-05
     
625. r,10017,dgarea,1.6143E-03+4.90E-05
     
626. PSTRES,ON
     
627. solve
     
628. finish
     
629. 
     
630. !******************************************
     
631. /solu
     
632. antype,buckle
     
633. bucopt,subsp,4,
     
634. outpr,nsol,all
     
635. solve
     
636. save
     
637. fini
     
638. /solu
     
639. expass,on
     
640. mxpand,4,,,
     
641. outpr,all
     
642. outres
     
643. solve
     
644. save
     
645. finish

复制代码

kakakey

有没有关于桁架门式起重机的ansys建模及分析呀？

happyyikuahai

你好  我也研究钢管砼  想请教你在ansys中钢管砼怎么建模   望指导  麻烦加一下547767612    先谢谢了