[转帖]钢管混凝土拱桥建模实例

AaronSpark

1. finish
   
2. /clear
   
3. /prep7
   
4. /title,the analyse of steel arch bridge
   
5. et,1,82
   
6. cyl4,0.6,0.6,0.6,,0.586
   
7. cyl4,0.6,2.4,0.6,,0.586
   
8. rectng,0.35,0.364,0.8,2.2
   
9. rectng,0.85,0.836,0.8,2.2
   
10. asel,all
    
11. aptn,all
    
12. adele,5,8,1
    
13. allsel
    
14. aadd,all
    
15. smrtsize,5
    
16. amesh,all
    
17. secwrite,gg,sect,,1
    
18. sectype,1,beam,mesh,
    
19. secoffset,cent,,,
    
20. secread,'gg','sect',,mesh
    
21. asel,all
    
22. aclear,all
    
23. adele,all,,,1
    
24. /replot
    
25. cyl4,0.6,0.6,0.586
    
26. cyl4,0.6,2.4,0.586
    
27. rectng,0.364,0.836,0.8,2.2
    
28. allsel
    
29. aadd,all
    
30. smrtsize,5
    
31. amesh,all
    
32. secwrite,hnt,sect,,1
    
33. sectype,2,beam,mesh,
    
34. secoffset,cent,,,
    
35. secread,'hnt','sect',,mesh
    
36. asel,all
    
37. aclear,all
    
38. adele,all,,,1
    
39. /replot
    
40. k,1,-0.5,0
    
41. k,2,-0.75,0.25
    
42. k,3,-0.75,0.5
    
43. k,4,-0.55,0.5
    
44. k,5,-0.55,0.85
    
45. k,6,0.55,0.85
    
46. k,7,0.55,0.5
    
47. k,8,0.75,0.5
    
48. k,9,0.75,0.25
    
49. k,10,0.5,0
    
50. a,1,2,3,4,5,6,7,8,9,10
    
51. smrtsize,5
    
52. amesh,all
    
53. secwrite,hl1,sect,,1
    
54. sectype,4,beam,mesh,
    
55. secoffset,cent,,,
    
56. secread,'hl1','sect',,mesh
    
57. asel,all
    
58. aclear,all
    
59. adele,all,,,1
    
60. /replot
    
61. k,1,-0.5,0
    
62. k,2,-0.5,1
    
63. k,3,-0.75,1.25
    
64. k,4,-0.75,1.5
    
65. k,5,0.75,1.5
    
66. k,6,0.75,1.25
    
67. k,7,0.5,1
    
68. k,8,0.5
    
69. a,1,2,3,4,5,6,7,8
    
70. smrtsize,5
    
71. amesh,all
    
72. secwrite,hl2,sect,,1
    
73. sectype,5,beam,mesh,
    
74. secoffset,cent,,,
    
75. secread,'hl2','sect',,mesh
    
76. asel,all
    
77. aclear,all
    
78. adele,all,,,1
    
79. /replot
    
80. k,1,-0.5,0
    
81. k,2,-0.5,1.04
    
82. k,3,-0.75,1.04+0.25
    
83. k,4,-0.75,1.04+0.5
    
84. k,5,-0.55,1.04+0.5
    
85. k,6,-0.55,1.04+0.5+0.35
    
86. k,7,0.55,1.04+0.5+0.35
    
87. k,8,0.55,1.04+0.5
    
88. k,9,0.75,1.04+0.5
    
89. k,10,0.75,1.04+0.25
    
90. k,11,0.5,1.04
    
91. k,12,0.5,0
    
92. k,13,-0.2,0.58
    
93. k,14,-0.2,1.04+0.5-0.3
    
94. k,15,0.2,1.04+0.5-0.3
    
95. k,16,0.2,0.58
    
96. a,1,2,3,4,5,6,7,8,9,10,11,12
    
97. a,13,14,15,16
    
98. asba,1,2
    
99. smrtsize,5
    
100. amesh,all
     
101. secwrite,hl3,sect,,1
     
102. sectype,6,beam,mesh,
     
103. secoffset,cent,,,
     
104. secread,'hl3','sect',,mesh
     
105. asel,all
     
106. aclear,all
     
107. adele,all,,,1
     
108. /replot  
     
109. ldele,all
     
110. kdele,all
     
111. k,1,-0.5,0
     
112. k,2,-0.5,1.23                     
     
113. k,3,-0.75,1.23+0.25               
     
114. k,4,-0.75,1.23+0.5               
     
115. k,5,-0.55,1.23+0.5               
     
116. k,6,-0.55,1.23+0.5+0.35           
     
117. k,7,0.55,1.23+0.5+0.35            
     
118. k,8,0.55,1.23+0.5                 
     
119. k,9,0.75,1.23+0.5                 
     
120. k,10,0.75,1.23+0.25               
     
121. k,11,0.5,1.23                     
     
122. k,12,0.5,0                        
     
123. k,13,-0.2,0.38                    
     
124. k,14,-0.3,0.38+0.2            
     
125. k,15,-0.3,1.23+0.5-0.3            
     
126. k,16,-0.2,1.23+0.5-0.1
     
127. k,17,0.2,1.23+0.5-0.1
     
128. k,18,0.3,1.23+0.5-0.1-0.2
     
129. k,19,0.3,0.38+0.2
     
130. k,20,0.2,0.38
     
131.                   
     
132. a,1,2,3,4,5,6,7,8,9,10,11,12      
     
133. a,13,14,15,16,17,18,19,20                     
     
134. asba,1,2                          
     
135. smrtsize,5                        
     
136. amesh,all                        
     
137. secwrite,hl4,sect,,1              
     
138. sectype,7,beam,mesh,              
     
139. secoffset,cent,,,                 
     
140. secread,'hl4','sect',,mesh        
     
141. asel,all                          
     
142. aclear,all                        
     
143. adele,all,,,1                     
     
144.                                  
     
145. /replot  
     
146. k,1,-0.5,0
     
147. k,2,-0.5,1.23                     
     
148. k,3,-0.75,1.23+0.25               
     
149. k,4,-0.75,1.23+0.5               
     
150. k,5,-0.55,1.23+0.5               
     
151. k,6,-0.55,1.23+0.5+0.35           
     
152. k,7,0.55,1.23+0.5+0.35            
     
153. k,8,0.55,1.23+0.5                 
     
154. k,9,0.75,1.23+0.5                 
     
155. k,10,0.75,1.23+0.25               
     
156. k,11,0.5,1.23                     
     
157. k,12,0.5,0                        
     
158. a,1,2,3,4,5,6,7,8,9,10,11,12
     
159. smrtsize,5                        
     
160. amesh,all                        
     
161. secwrite,hl5,sect,,1              
     
162. sectype,8,beam,mesh,              
     
163. secoffset,cent,,,                 
     
164. secread,'hl5','sect',,mesh        
     
165. asel,all                          
     
166. aclear,all                        
     
167. adele,all,,,1                     
     
168.                                  
     
169. /replot     
     
170. etdele,1
     
171. /replot
     
172. et,1,beam44
     
173. mp,ex,1,2.1e11
     
174. mp,dens,1,7800
     
175. mp,prxy,1,0.3
     
176. n,90000,0,0,30
     
177. et,2,beam44
     
178. mp,ex,2,3.5e10
     
179. mp,dens,2,2600
     
180. mp,prxy,2,0.1667
     
181. et,3,beam44
     
182. mp,ex,3,3.0e10
     
183. mp,dens,3,2600
     
184. mp,prxy,3,0.1667
     
185. et,4,beam44
     
186. mp,ex,4,3.0e10
     
187. mp,dens,4,2600
     
188. mp,prxy,4,0.1667
     
189. et,5,beam44
     
190. mp,ex,5,2.1e11
     
191. mp,dens,5,7800
     
192. mp,prxy,5,0.3
     
193. et,6,link10
     
194. mp,ex,6,1.9e11
     
195. mp,dens,6,7800
     
196. mp,prxy,6,0.3
     
197. keyopt,6,3,0
     
198. et,7,beam44
     
199. mp,ex,7,3.0e10
     
200. mp,dens,7,2600
     
201. mp,prxy,7,0.1667
     
202. et,8,beam44
     
203. mp,ex,8,3.3e10
     
204. mp,dens,8,2600
     
205. mp,prxy,8,0.1667
     
206. et,9,beam44
     
207. mp,ex,9,3.0e10
     
208. mp,dens,9,2600
     
209. mp,prxy,9,0.1667
     
210. et,10,link10
     
211. mp,ex,10,2.0e11
     
212. mp,dens,10,7800
     
213. mp,prxy,10,0.3
     
214. keyopt,10,3,0
     
215. et,11,beam44
     
216. mp,ex,11,2.85e10
     
217. mp,dens,11,2500
     
218. mp,prxy,11,0.1667
     
219. n,90000,0,0,30
     
220. 
     
221. n,1,0,2.025
     
222. n,2,6.25,2.025
     
223. n,3,6.25+5.1,2.025
     
224. n,4,6.25+5.1+0.2+1.9/2,2.025
     
225. n,5,6.25+5.1+2.3,2.025
     
226. n,6,6.25+5.1+2.3+2.23,2.025
     
227. nsym,x,10,2,6,1
     
228. type,4
     
229. r,1
     
230. real,1
     
231. mat,4
     
232. secnum,7
     
233. e,1,2
     
234. secnum,6
     
235. e,2,3
     
236. secnum,5
     
237. e,3,4
     
238. e,4,5
     
239. secnum,4
     
240. e,5,6
     
241. secnum,7
     
242. e,1,12
     
243. secnum,6
     
244. e,12,13
     
245. secnum,5
     
246. e,13,14
     
247. e,14,15
     
248. secnum,4
     
249. e,15,16
     
250. egen,2,100,1,10,1,,,,,,0,4.7+1.25/2+1.5/2
     
251. egen,17,100,11,20,1,,,,,,0,5.3+1.5
     
252. egen,2,100,171,180,1,,,,,,0,4.7+1.5/2+1.25/2
     
253. cm,hl,elem
     
254. esel,none
     
255. type,3
     
256. r,2
     
257. real,2
     
258. mat,3
     
259. sectype,10,beam,rect,zl
     
260. secdata,0.4,1.2
     
261. secnum,10
     
262. e,1,101
     
263. e,2,102
     
264. e,3,103
     
265. e,5,105
     
266. e,12,112
     
267. e,13,113
     
268. e,15,115
     
269. *do,i,0,15,1
     
270. e,101+i*100,201+i*100
     
271. e,102+i*100,202+i*100
     
272. e,104+i*100,204+i*100
     
273. e,112+i*100,212+i*100
     
274. e,114+i*100,214+i*100
     
275. *enddo
     
276. e,1701,1801
     
277. e,1702,1802
     
278. e,1703,1803
     
279. e,1705,1805
     
280. e,1712,1812
     
281. e,1713,1813
     
282. e,1715,1815
     
283. cm,zl,elem
     
284. esel,none
     
285. 
     
286. *dim,y,array,16
     
287. *dim,cs,array,16
     
288. *dim,z,array,16
     
289. y(1)=8.1
     
290. y(2)=14.9
     
291. y(3)=21.7
     
292. y(4)=28.5
     
293. y(5)=35.3
     
294. y(6)=42.1
     
295. y(7)=48.9
     
296. y(8)=55.7
     
297. y(9)=62.5
     
298. y(11)=21.7
     
299. y(12)=28.5
     
300. y(13)=42.1
     
301. y(14)=48.9
     
302. y(15)=59.1  
     
303. y(16)=2.025
     
304. f1=25
     
305. m1=1.1
     
306. k=log(m1+sqrt(m1*m1-1))
     
307. *do,i,1,16,1
     
308. cs(i)=1-y(i)*2/125
     
309. z(i)=25-f1/(m1-1)*(cosh(k*cs(i))-1)
     
310. *enddo
     
311. *dim,ygd,array,23
     
312. ygd(1) =0.0040004 +0.0129  
     
313. ygd(2) =0.010685 +0.0222                        
     
314. ygd(3) =0.019483 +0.0303  
     
315. ygd(4) =0.029219 +0.0369  
     
316. ygd(5) =0.039022 +0.0423  
     
317. ygd(6) =0.048243 +0.0465  
     
318. ygd(7) =0.05639 +0.0495         
     
319. ygd(8) =0.06307 +0.0515         
     
320. ygd(9) =0.068014 +0.0521  
     
321. ygd(10) =0.071003 +0.0515  
     
322. ygd(11) =0.071739 +0.0495  
     
323. ygd(12) =0.069962 +0.0465  
     
324. ygd(13) =0.065514 +0.423   
     
325. ygd(14) =0.058289 +0.0369  
     
326. ygd(15) =0.048248 +0.0304  
     
327. ygd(16) =0.035507 +0.0224  
     
328. ygd(17) =0.20415 +0.0123         
     
329. ygd(18) =0.06577 +0.040234        
     
330. ygd(19) =0.069769 +0.041215
     
331. ngen,2,10000,104 ,114 ,10,0,0,z(1)+ygd(1 )
     
332. ngen,2,10000,204 ,214 ,10,0,0,z(2)+ygd(2 )
     
333. ngen,2,10000,304 ,314 ,10,0,0,z(3)+ygd(3 )
     
334. ngen,2,10000,404 ,414 ,10,0,0,z(4)+ygd(4 )
     
335. ngen,2,10000,504 ,514 ,10,0,0,z(5)+ygd(5 )
     
336. ngen,2,10000,604 ,614 ,10,0,0,z(6)+ygd(6 )
     
337. ngen,2,10000,704 ,714 ,10,0,0,z(7)+ygd(7 )
     
338. ngen,2,10000,804 ,814 ,10,0,0,z(8)+ygd(8 )
     
339. ngen,2,10000,904 ,914 ,10,0,0,z(9)+ygd(9 )
     
340. ngen,2,10000,1004 ,1014 ,10,0,0,z(8)+ygd(10)
     
341. ngen,2,10000,1104 ,1114 ,10,0,0,z(7)+ygd(11)
     
342. ngen,2,10000,1204 ,1214 ,10,0,0,z(6)+ygd(12)
     
343. ngen,2,10000,1304 ,1314 ,10,0,0,z(5)+ygd(13)
     
344. ngen,2,10000,1404 ,1414 ,10,0,0,z(4)+ygd(14)
     
345. ngen,2,10000,1504 ,1514 ,10,0,0,z(3)+ygd(15)
     
346. ngen,2,10000,1604 ,1614 ,10,0,0,z(2)+ygd(16)
     
347. ngen,2,10000,1704 ,1714 ,10,0,0,z(1)+ygd(17)
     
348. ngen,2,3000,104 ,114 ,10,0,0,0
     
349. ngen,2,3000,204 ,214 ,10,0,0,0
     
350. ngen,2,3000,304 ,314 ,10,0,0,0
     
351. ngen,2,3000,404 ,414 ,10,0,0,0
     
352. ngen,2,3000,504 ,514 ,10,0,0,0
     
353. ngen,2,3000,604 ,614 ,10,0,0,0
     
354. ngen,2,3000,704 ,714 ,10,0,0,0
     
355. ngen,2,3000,804 ,814 ,10,0,0,0
     
356. ngen,2,3000,904 ,914 ,10,0,0,0
     
357. ngen,2,3000,1004,1014 ,10,0,0,0
     
358. ngen,2,3000,1104,1114 ,10,0,0,0
     
359. ngen,2,3000,1204,1214 ,10,0,0,0
     
360. ngen,2,3000,1304,1314 ,10,0,0,0
     
361. ngen,2,3000,1404,1414 ,10,0,0,0
     
362. ngen,2,3000,1504,1514 ,10,0,0,0
     
363. ngen,2,3000,1604,1614 ,10,0,0,0
     
364. ngen,2,3000,1704,1714 ,10,0,0,0
     
365. ngen,2,6000,104 ,114 ,10,0,0,z(1)+ygd(1 )
     
366. ngen,2,6000,204 ,214 ,10,0,0,z(2)+ygd(2 )
     
367. ngen,2,6000,304 ,314 ,10,0,0,z(3)+ygd(3 )
     
368. ngen,2,6000,404 ,414 ,10,0,0,z(4)+ygd(4 )
     
369. ngen,2,6000,504 ,514 ,10,0,0,z(5)+ygd(5 )
     
370. ngen,2,6000,604 ,614 ,10,0,0,z(6)+ygd(6 )
     
371. ngen,2,6000,704 ,714 ,10,0,0,z(7)+ygd(7 )
     
372. ngen,2,6000,804 ,814 ,10,0,0,z(8)+ygd(8 )
     
373. ngen,2,6000,904 ,914 ,10,0,0,z(9)+ygd(9 )
     
374. ngen,2,6000,1004,1014 ,10,0,0,z(8)+ygd(10)
     
375. ngen,2,6000,1104,1114 ,10,0,0,z(7)+ygd(11)
     
376. ngen,2,6000,1204,1214 ,10,0,0,z(6)+ygd(12)
     
377. ngen,2,6000,1304,1314 ,10,0,0,z(5)+ygd(13)
     
378. ngen,2,6000,1404,1414 ,10,0,0,z(4)+ygd(14)
     
379. ngen,2,6000,1504,1514 ,10,0,0,z(3)+ygd(15)
     
380. ngen,2,6000,1604,1614 ,10,0,0,z(2)+ygd(16)
     
381. ngen,2,6000,1704,1714 ,10,0,0,z(1)+ygd(17)                                                            
     
382. dgarea=109*3.1415926*0.25*0.007*0.007
     
383. dgyl=0.3*1670*1e6
     
384. dgl=dgyl*dgarea
     
385. edg=1.95e11
     
386. *dim,L0,array,17
     
387. *dim,L,array,17
     
388. *dim,deltL,array,17
     
389. *dim,istrn,array,17
     
390. *dim,fdg,array,17
     
391. *do,i,1,9,1
     
392. L0(i)=z(i)
     
393. L0(18-i)=z(i)
     
394. *enddo
     
395. fdg(1 )=dgl
     
396. fdg(2 )=dgl
     
397. fdg(3 )=dgl
     
398. fdg(4 )=dgl
     
399. fdg(5 )=dgl
     
400. fdg(6 )=dgl
     
401. fdg(7 )=dgl
     
402. fdg(8 )=dgl
     
403. fdg(9 )=dgl
     
404. fdg(10)=dgl
     
405. fdg(11)=dgl
     
406. fdg(12)=dgl
     
407. fdg(13)=dgl
     
408. fdg(14)=dgl
     
409. fdg(15)=dgl
     
410. fdg(16)=dgl
     
411. fdg(17)=dgl
     
412. type,6
     
413. mat,6
     
414. *do,i,1,17,1
     
415. deltL(i)=fdg(i)*L0(i)/(edg*dgarea)
     
416. L(i)=L0(i)+deltL(i)
     
417. istrn(i)=deltL(i)/L(i)
     
418. r,i+1000,dgarea,istrn(i)
     
419. *enddo
     
420. *do,i,1,17,1
     
421. real,1+1000
     
422. e,i*100+4+3000,i*100+4+6000
     
423. e,i*100+10+4+3000,i*100+10+4+6000
     
424. *enddo
     
425. cm,dg,elem
     
426. esel,none
     
427. allsel
     
428. real,1018
     
429. n,20001, 0 ,21.7 ,z(11)        
     
430. n,20002, 0 ,28.5 ,z(12)        
     
431. n,20003, 0 ,42.1 ,z(13)        
     
432. n,20004, 0 ,48.9 ,z(14)        
     
433. n,20005, 0 ,59.1 ,z(15)        
     
434. n,20006, 0 ,65.9 ,z(15)        
     
435. n,20007, 0 ,76.1 ,z(14)        
     
436. n,20008, 0 ,82.9 ,z(13)        
     
437. n,20009, 0 ,96.5 ,z(12)        
     
438. n,20010, 0 ,103.3 ,z(11)        
     
439. n,20011, 12.5 ,21.7 ,z(11)        
     
440. n,20012, 12.5 ,28.5 ,z(12)        
     
441. n,20013, 12.5 ,42.1 ,z(13)        
     
442. n,20014, 12.5 ,48.9 ,z(14)        
     
443. n,20015, 12.5 ,59.1 ,z(15)+ygd(18)
     
444. n,20016, 12.5 ,65.9 ,z(15)+ygd(19)
     
445. n,20017, 12.5 ,76.1 ,z(14)        
     
446. n,20018, 12.5 ,82.9 ,z(13)        
     
447. n,20019, 12.5 ,96.5 ,z(12)        
     
448. n,20020, 12.5 ,103.3 ,z(11)        
     
449. n,20021, -12.5 ,21.7 ,z(11)        
     
450. n,20022, -12.5 ,28.5 ,z(12)        
     
451. n,20023, -12.5 ,42.1 ,z(13)        
     
452. n,20024, -12.5 ,48.9 ,z(14)        
     
453. n,20025, -12.5 ,59.1 ,z(15)+ygd(18)
     
454. n,20026, -12.5 ,65.9 ,z(15)+ygd(19)
     
455. n,20027, -12.5 ,76.1 ,z(14)        
     
456. n,20028, -12.5 ,82.9 ,z(13)        
     
457. n,20029, -12.5 ,96.5 ,z(12)        
     
458. n,20030, -12.5 ,103.3 ,z(11)        
     
459. n,20000, 0 ,62.5 ,z(15)
     
460. type,5
     
461. mat,5
     
462. sectype,500,beam,ctube,windc
     
463. secdata,0.486,0.50
     
464. secnum,500
     
465. e,20002,10404,90000
     
466. e,20002,10414,90000
     
467. e,20002,10304,90000
     
468. e,20002,10314,90000
     
469. e,20004,10704,90000
     
470. e,20004,10714,90000
     
471. e,20004,10604,90000
     
472. e,20004,10614,90000
     
473. e,20000,20015,90000
     
474. e,20000,20016,90000
     
475. e,20000,20025,90000
     
476. e,20000,20026,90000
     
477. e,20007,11104,90000
     
478. e,20007,11114,90000
     
479. e,20007,11204,90000
     
480. e,20007,11214,90000
     
481. e,20009,11404,90000
     
482. e,20009,11414,90000
     
483. e,20009,11504,90000
     
484. e,20002,11514,90000
     
485. cm,fc,elem
     
486. esel,none
     
487. n,30000,12.5,0,0
     
488. n,30001,12.5,125,0
     
489. n,30002,12.5,2.025,z(16)
     
490. n,30003,12.5,125-2.025,z(16)
     
491. n,40000,12.5,0,0
     
492. n,40001,12.5,125,0
     
493. n,40002,12.5,2.025,z(16)
     
494. n,40003,12.5,125-2.025,z(16)
     
495. *dim,zgzx,array,23
     
496. *dim,ygzx,array,23
     
497. zgzx(1 )=30000
     
498. zgzx(2 )=30002
     
499. zgzx(3 )=10104
     
500. zgzx(4 )=10204
     
501. zgzx(5 )=10304
     
502. zgzx(6 )=10404
     
503. zgzx(7 )=10504
     
504. zgzx(8 )=10604
     
505. zgzx(9 )=10704
     
506. zgzx(10 )=10804
     
507. zgzx(11 )=20015
     
508. zgzx(12 )=10904
     
509. zgzx(13 )=20016
     
510. zgzx(14 )=11004
     
511. zgzx(15 )=11104
     
512. zgzx(16 )=11204
     
513. zgzx(17 )=11304
     
514. zgzx(18 )=11404
     
515. zgzx(19 )=11504
     
516. zgzx(20 )=11604
     
517. zgzx(21 )=11704
     
518. zgzx(22 )=30003
     
519. zgzx(23 )=30001
     
520. zgzx(1 )=40000
     
521. zgzx(2 )=40002
     
522. zgzx(3 )=10114
     
523. zgzx(4 )=10214
     
524. zgzx(5 )=10314
     
525. zgzx(6 )=10414
     
526. zgzx(7 )=10514
     
527. zgzx(8 )=10614
     
528. zgzx(9 )=10714
     
529. zgzx(10 )=10814
     
530. zgzx(11 )=20025
     
531. zgzx(12 )=10914
     
532. zgzx(13 )=20026
     
533. zgzx(14 )=11014
     
534. zgzx(15 )=11114
     
535. zgzx(16 )=11214
     
536. zgzx(17 )=11314
     
537. zgzx(18 )=11414
     
538. zgzx(19 )=11514
     
539. zgzx(20 )=11614
     
540. zgzx(21 )=11714
     
541. zgzx(22 )=40003
     
542. zgzx(23 )=40001
     
543. n,90001,12.5,0,30
     
544. n,90002,-12.5,0,30
     
545. type,1
     
546. mat,1
     
547. r,2000
     
548. secnum,1
     
549. *do,i,1,22
     
550. e,zgzx(i),zgzx(i+1),90001
     
551. e,ygzx(i),ygzx(i+1),90002
     
552. *enddo
     
553. cm,gzgg,elem
     
554. esel,none
     
555. type,2
     
556. mat,2
     
557. r,2001
     
558. secnum,2
     
559. *do,i,1,22
     
560. e,zgzx(i),zgzx(i+1),90001
     
561. e,ygzx(i),ygzx(i+1),90002
     
562. *enddo
     
563. cm,gzhnt,elem
     
564. esel,none
     
565. type,9
     
566. r,6
     
567. real,6
     
568. mat,9
     
569. sectype,70,beam,rect,dizuo
     
570. secdata,1.4,4
     
571. secnum,70
     
572. e,30001,30003,90001
     
573. e,40001,40003,90002
     
574. e,30000,30002,90001
     
575. e,30000,40002,90002
     
576. cm,gzdz,elem
     
577. esel,none
     
578. ngen,2,2000,1,16,1,0,-2.025,0
     
579. ngen,2,100,1801,1816,1,0,2.025,0
     
580. type,7
     
581. mat,7
     
582. r,7
     
583. real,7
     
584. sectype,60,beam,rect,gail
     
585. secdata,4,1.6
     
586. 
     
587. !!!!!!!!!!*****************************以上已经校核过!

复制代码

AaronSpark

secnum,60
n,90005,0,0,20
n,90006,0,125,20
e,1901, 1902, 90006
e,1902, 1903, 90006
e,1903, 1904, 90006
e,1904, 1905, 90006
e,1905, 1906, 90006
e,1906, 1907, 90006
e,1907, 1908, 90006
e,1908, 1909, 90006
e,1909, 1910, 90006
e,1910, 1911, 90006
e,1911, 1912, 90006
e,1912, 1913, 90006
e,1913, 1914, 90006
e,1914, 1915, 90006
e,1915, 1916, 90006
e,2001, 2002, 90005
e,2002, 2003, 90005
e,2003, 2004, 90005
e,2004, 2005, 90005
e,2005, 2006, 90005
e,2006, 2007, 90005
e,2007, 2008, 90005
e,2008, 2009, 90005
e,2009, 2010, 90005
e,2010, 2011, 90005
e,2011, 2012, 90005
e,2012, 2013, 90005
e,2013, 2014, 90005
e,2014, 2015, 90005
e,2015, 2016, 90005
cm,gailiang,elem
esel,none
n,90003,14.5,0,0
n,90004,-14,5,0,0
type,8
r,8
real,8
sectype,50,beam,hrec,dun
secdata,3,4,0.8,0.8,0.8,0.8
ngen,2,10000,2004,2014,10,0,0,-10
ngen,2,10000,1904,1914,10,0,0,-10
secnum,50
e,2004,12004,90003
e,2014,12014,90004
e,1904,11904,90003
e,1914,11914,90004
cm,dunzhu,elem
esel,none
type,11
mat,11
r,1
sectype,40,beam,rect,chengtai
secdata,12.7,9
r,11
real,11
secnum,40
ngen,2,20000,2004,2014,10,0,0,-11.5
ngen,2,20000,1904,1914,10,0,0,-11.5
ngen,2,30000,2004,2014,10,0,0,-13
ngen,2,30000,1904,1914,10,0,0,-13  
e,12004,22004,90003
e,11904,21904,90003
e,12014,22014,90004
e,11904,21904,90004
e,22004,32004,90003
e,21904,31904,90003
e,22014,32014,90004
e,21904,31904,90004
cm,chtai,elem
esel,none
type,9
mat,9
r,9
real,9
n,90007,12.5,0,10
n,90008,-12.5,0,10
cm,hxl,elem
esel,none

n,50000,6.15+1.85,2.7,-13
n,50001,6.15+1.85,-2.7,-13
n,50002,6.15+1.85+4.5,2.7,-13
n,50003,6.15+1.85+4.5,-2.7,-13
n,50004,6.15+1.27-1.85,2.7,-13
n,50005,6.15+1.27-1.85,-2.7,-13
nsel,s,,,50000,50005,1
nsym,x,10,50000,50005,1
allsel
ngen,2,100,50000,50015,1,0,125,-2.5
ngen,10,1000,50000,50015,1,0,0,-2.5
esel,none
type,11
mat,11
r,11
real,11
sectype,20,beam,csolid,pier
secdata,0.9
secnum,20
*do,i,0,8,1
*do,k,0,5,1
e,50000+k+i*1000,50000+k+(i+1)*1000,90003
e,50000+k+10+i*1000,50000+k+10+(i+1)*1000,90003
e,50000+k+100+i*1000,50000+k+100+(i+1)*1000,90003
e,50000+k+110+i*1000,50000+k+110+(i+1)*1000,90003
*enddo
*enddo
cm,zhuang,elem
esel,none
d=1.8
kf=0.9
k0=1+1/d
b1=0.6
L1=3.6
h1=3*(d+1)
k=b1+(1-b1)/0.6*L1/h1
B=Kf*k0*k*d

*dim,m,array,5
m(1)=8
m(2)=15
m(3)=20
m(4)=25
m(5)=4000
*dim,h,array,5
h(1)=5
h(2)=5
h(3)=5
h(4)=5
h(5)=10
*dim,c,array,4
c(1)=m(1)*h(1)
c(2)=m(2)*(h(1)+h(2))
c(3)=m(3)*(h(1)+h(2)+h(3))
c(4)=m(4)*(h(1)+h(2)+h(3)+h(4))
*dim,e,array,4
e(1)=0.5*c(1)*h(1)*B*1e6
*do,i,2,4,1
e(i)=0.5*(c(i-1)+c(i))*h(1)*B*1e6
*enddo
et,20,combin14
r,20,e(1)
r,30,e(2)
r,40,e(3)
r,50,e(4)
r,60,m(5)*3.1415926*B*B/4*1e6
ngen,2,10000,51000,51115,1,1,0,0
ngen,2,20000,51000,51115,1,0,1,0
ngen,2,30000,51000,51115,1,0,0,1
type,20
real,20
*do,i,0,5,1
e,51000+i,51000+i+10000
e,51000+i+10,51000+i+10+10000
e,51100+i,51100+i+10000
e,51100+i+10,51100+i+10+10000
e,51000+i,51000+i+20000
e,51000+i+10,51000+i+10+20000
e,51100+i,51100+i+20000
e,51100+i+10,51100+i+10+20000
*enddo
cm,spring1,elem
esel,none
ngen,2,10000,53000,53115,1,1,0,0
ngen,2,20000,53000,53115,1,0,1,0
ngen,2,30000,53000,53115,1,0,0,1
type,20
real,30
*do,i,0,5,1
e,53000+i,53000+i+10000
e,53000+i+10,53000+i+10+10000
e,53100+i,53100+i+10000
e,53100+i+10,53100+i+10+10000
e,53000+i,53000+i+20000
e,53000+i+10,53000+i+10+20000
e,53100+i,53100+i+20000
e,53100+i+10,53100+i+10+20000
*enddo
cm,spring2,elem
esel,none

ngen,2,10000,55000,55115,1,1,0,0  
ngen,2,20000,55000,55115,1,0,1,0  
                                 
type,20                           
real,40                           
*do,i,0,5,1                       
e,55000+i,55000+i+10000           
e,55000+i+10,55000+i+10+10000     
e,55100+i,55100+i+10000           
e,55100+i+10,55100+i+10+10000     
e,55000+i,55000+i+20000           
e,55000+i+10,55000+i+10+20000     
e,55100+i,55100+i+20000           
e,55100+i+10,55100+i+10+20000     
*enddo                           
cm,spring3,elem                  
esel,none                        
ngen,2,10000,57000,57115,1,1,0,0     
ngen,2,20000,57000,57115,1,0,1,0     
                                    
type,20                              
real,50                              
*do,i,0,5,1                          
e,57000+i,57000+i+10000              
e,57000+i+10,57000+i+10+10000        
e,57100+i,57100+i+10000              
e,57100+i+10,57100+i+10+10000        
e,57000+i,57000+i+20000              
e,57000+i+10,57000+i+10+20000        
e,57100+i,57100+i+20000              
e,57100+i+10,57100+i+10+20000        
*enddo                              
cm,spring4,elem                     
esel,none
ngen,2,10000,58000,58115,1,1,0,0                                
ngen,2,20000,58000,58115,1,0,1,0     
                                    
type,20                              
real,60                              
*do,i,0,5,1                          
e,58000+i,58000+i+10000              
e,58000+i+10,58000+i+10+10000        
e,58100+i,58100+i+10000              
e,58100+i+10,58100+i+10+10000        
e,58000+i,58000+i+20000              
e,58000+i+10,58000+i+10+20000        
e,58100+i,58100+i+20000              
e,58100+i+10,58100+i+10+20000        
*enddo                              
cm,spring5,elem                     
esel,none
ngen,2,10000,59000,59115,1,1,0,0                                 
ngen,2,10000,59000,59115,1,1,0,0
ngen,2,20000,59000,59115,1,0,1,0
ngen,2,20000,59000,59115,1,0,1,0
ngen,2,30000,59000,59115,1,0,0,-1
ngen,2,30000,59000,59115,1,0,0,-1
   
                                    
type,20                              
real,60                              
*do,i,0,5,1                          
e,59000+i,59000+i+10000              
e,59000+i+10,59000+i+10+10000        
e,59100+i,59100+i+10000              
e,59100+i+10,59100+i+10+10000        
e,59000+i,59000+i+20000              
e,59000+i+10,59000+i+10+20000        
e,59100+i,59100+i+20000              
e,59100+i+10,59100+i+10+20000
e,59000+i,59000+i+30000      
e,59000+i+10,59000+i+10+30000
e,59100+i,59100+i+30000     
e,59100+i+10,59100+i+10+30000
      
*enddo                              
cm,spring6,elem                     
allsel

allsel
nsel,s,loc,z,-50,-0.01
nsel,r,loc,y,-10,10
cm,zdnode,node
esln,s,,all
cm,zdelem,elem
allsel
nsel,u,,,zdnode
nsym,y,100000,all
allsel
esel,u,,,zdelem
esym,,100000,all
allsel

xgarea=8*27*7*3.1415926*0.25*0.005*0.005
exg=1.95e11
*dim,fxg,array,2
*dim,Llast,array,2
*dim,Listrn,array,2
*dim,detL,array,2
xgyl=0.3*1860*1e6
xgl=xgyl*xgarea
fxg(1)=xgl
fxg(2)=xgl
Lzero=250
*do,i,1,2,1
detL(i)=Lzero*Fxg(i)/(exg*xgarea)
Llast(i)=Lzero+detL(i)
Listrn(i)=detl(i)/Llast(i)
r,i+100,xgarea,Listrn(i)
*enddo
type,10
mat,10
real,101
e,30001,130001
real,102
e,40001,140001
cm,xigan,elem
esel,none
allesl
cp, 1, ux ,2004 ,30000
cp, 1, uy ,2004 ,30000                     
cp, next, uz ,2004 ,30000
cp, next, roty ,2004 ,30000
cp, next, rotz ,2004 ,30000
cp, next, ux ,1904 ,30001
cp, next, uy ,1904 ,30001
cp, next, uz ,1904 ,30001
cp, next, roty ,1904 ,30001
cp, next, rotz ,1904 ,30001
cp, next, ux ,2014 ,40000
cp, next, uy ,2014 ,40000
cp, next, uz ,2014 ,40000
cp, next, roty ,2014 ,40000
cp, next, rotz ,2014 ,40000
cp, next, ux ,1914 ,40001
cp, next, uy ,1914 ,40001
cp, next, uz ,1914 ,40001
cp, next, roty ,1914 ,40001
cp, next, rotz ,1914 ,40001
cp, next, ux ,101914 ,140001
cp, next, uy ,101914 ,140001
cp, next, uz ,101914 ,140001
cp, next, roty ,101914 ,140001
cp, next, rotz ,101914 ,140001
cp, next, ux ,101904 ,130001
cp, next, uy ,101904 ,130001
cp, next, uz ,101904 ,130001
cp, next, roty ,101904 ,130001
cp, next, rotz ,101904 ,130001
allsel
nsel,s,,,30000,130000,100000
nsel,a,,,40000,140000,100000
nummrg,node
allsel
*do,i,1,3,1
cp,next,all,1800+i,1900+i
*enddo
allsel
*do,i,5,6,1
cp,next,all,1800+i,1900+i
*enddo
allsel
*do,i,1,2,1              
cp,next,all,1800+i+11,1900+i+11
cp,next,all,i+11,2000+i+11,i+100000+11
*enddo                  
allsel
*do,i,4,5,1                                                
cp,next,all,1800+i+11,1900+i+11         
cp,next,all,i+11,2000+i+11,i+100000+11   
*enddo                                   
allsel   
*do,i,1,3,1                                                              
cp,next,all,i,2000+i,100000+i
*enddo                                   
allsel  
*do,i,5,6,1                                                              
cp,next,all,i,2000+i,100000+i
*enddo                                   
allsel
*do,i,1,3,1                                                              
cp,next,all,101800+i,101900+i
*enddo                                   
allsel   
*do,i,5,6,1                                                              
cp,next,all,101800+i,101900+i
*enddo                                   
allsel
*do,i,1,2,1                                                              
cp,next,all,101800+i+11,101900+i+11
*enddo                                   
allsel
*do,i,4,5,1                                                              
cp,next,all,101800+i+11,101900+i+11
*enddo                                   
allsel
nsel,s,loc,z,-13
nsel,r,loc,y,-3,3
nsel,r,loc,x,8,20
cp,next,all,all
allsel   
nsel,s,loc,z,-13
nsel,r,loc,y,-3,3
nsel,r,loc,x,-8,-20
cp,next,all,all
allsel
nsel,s,loc,z,-13
nsel,r,loc,y,122,128
nsel,r,loc,x,-8,-20
cp,next,all,all
allsel
nsel,s,loc,z,-13
nsel,r,loc,y,122,128
nsel,r,loc,x,8,20
cp,next,all,all
allsel
nsel,s,loc,z,-13
nsel,r,loc,y,-122,-128
nsel,r,loc,x,8,20
cp,next,all,all
allsel
nsel,s,loc,z,-13
nsel,r,loc,y,-122,-128
nsel,r,loc,x,-8,-20
cp,next,all,all
allsel   
*do,i,1,17,1
cp, next, ux, i*100+4                 ,i*100+4+3000
cp, next, ux, i*100+4+6000         ,i*100+4+10000
cp, next, uy, i*100+4                 ,i*100+4+3000
cp, next, uy, i*100+4+6000         ,i*100+4+10000
cp, next, uz, i*100+4                 , i*100+4+3000
cp, next, uz, i*100+4+6000         ,i*100+4+10000
cp, next, ux, i*100+4+10         ,i*100+4+3000+10
cp, next, ux, i*100+4+6000+10         ,i*100+4+10000+10
cp, next, uy, i*100+4+10         ,i*100+4+3000+10
cp, next, uy, i*100+4+6000+10         ,i*100+4+10000+10
cp, next, uz, i*100+4+10         ,i*100+4+3000+10
cp, next, uz, i*100+4+6000+10         ,i*100+4+10000+10
cp, next, ux, i*100+4+100000         ,i*100+4+3000+100000
cp, next, ux, i*100+4+6000+100000 ,i*100+4+10000+100000
cp, next, uy, i*100+4+100000         , i*100+4+3000+100000
cp, next, uy, i*100+4+6000+100000 ,i*100+4+10000+100000
cp, next, uz, i*100+4+100000         ,i*100+4+3000+100000
cp, next, uz, i*100+4+6000+100000 ,i*100+4+10000+100000
cp, next, ux, i*100+4+100000+10 ,i*100+4+3000+100000+10
cp, next, ux, i*100+4+6000+100000 ,i*100+4+10000+100000+10
cp, next, uy, i*100+4+100000+10 ,i*100+4+3000+100000+10
cp, next, uy, i*100+4+6000+100000 ,i*100+4+10000+100000+10
cp, next, uz, i*100+4+100000+10 ,i*100+4+3000+100000+10
cp, next, uz, i*100+4+6000+100000 ,i*100+4+10000+100000+10
*enddo
allsel
allsel
nsel,u,,,erase
nsel,s,loc,x,13.5
nsel,a,loc,x,9
nsel,a,loc,x,18
nsel,a,loc,x,-11.5
nsel,a,loc,x,-7
nsel,a,loc,x,-16
nsel,a,loc,y,3.7
nsel,a,loc,y,-1.7
nsel,a,loc,y,128.7
nsel,a,loc,y,123.3
nsel,a,loc,y,-128.7
nsel,a,loc,y,-123.3
nsel,a,loc,z,-40,-13
d,all,all
allsel
nsel,s,loc,z,-36,5
d,all,all
/color,pbak,off
/replot
allsel
nsel,s,loc,y,0
nsel,r,loc,z,0
nummrg,node,all
esel,none
allsel
/solu
allsel
mpdele,dens,all
mp,dens,3,2600
mp,dens,4,2600
acel,,,9.8
pstres,on
solve
save

/solu
allsel
mpdele,dens,3,11,1
mp,dens,1,7800
mp,dens,2,2600
*do,i,1,17,1
r,1000+i,dgarea,0
*enddo
r,101,xgarea,0
r,102,xgarea,0
acel,,,9.8
solve
/solu
allsel
*do,i,1,17,1
r,1000+1,dgarea,0
*enddo
r,101,xgarea,0
r,102,xgarea,0
r,101,xgarea,0
r,102,xgarea,0
qmb=(250/1.5)*32/6.8*5*1000*9.8
acel,,,9.8
mpdele,dens,3,4,1
mp,dens,1,7800
mp,dens,2,2600
mp,dens,3,4167.7
mp,dens,4,4167.7
mp,dens,5,7800
mp,dens,6,7800
mp,dens,7,2600
mp,dens,8,2600
mp,dens,9,2600
mp,dens,10,7800
mp,dens,11,2600
solve
/solu
allsel
acel,,,9.8
r,101,xgarea,Listrn(1)*1.3
r,102,xgarea,Listrn(2)*1.3
pstres,on
solve

/solu
allsel
r, 1001, dgarea,1.6787e-03+7.64e-05
r, 1002, dgarea,1.7538e-03+8.48e-07
r, 1003, dgarea,1.7547e-03+1.03e-06
r, 1004, dgarea,1.7400e-03+2.71e-06
r, 1005, dgarea,1.7292e-03+6.17e-06
r, 1006, dgarea,1.7303e-03+6.37e-06
r, 1007, dgarea,1.7267e-03+6.47e-06
r, 1008, dgarea,1.7279e-03+7.69e-06
r, 1009, dgarea,1.7365e-03+7.91e-06
r, 1010, dgarea,1.7277e-03+7.81e-06
r, 1011, dgarea,1.7263e-03+6.56e-06
r, 1012, dgarea,1.7299e-03+6.24e-06
r, 1013, dgarea,1.7290e-03+5.53e-06
r, 1014, dgarea,1.7412e-03+2.03e-06
r, 1015, dgarea,1.7504e-03+3.03e-06
r, 1016, dgarea,1.7587e-03+1.28e-05
r, 1017, dgarea,1.6143e-03+4.90e-05
pstres,on
solve
finish
/solu
antype,buckle
bucopt,subsp,4
outpr,nsol,all
save
solve
finish
/solu
expass,on
mxpand,4,,,
outpr,all
outres
solve
finish

Crystalfly

谢谢楼主，请问楼主有没有钢管相贯节点的命令流实例，谢谢

AaronSpark

命令流没有，gui得到是有一个，看帖子
http://forum.vibunion.com/forum/viewthread.php?tid=17479

yinyi

谢谢楼主

wuwei8696

谢谢！！非常