声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2318|回复: 5

[结构分析] [转帖]钢管混凝土拱桥建模实例

[复制链接]
发表于 2006-5-22 05:11 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
  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. n,1,0,2.025
  221. n,2,6.25,2.025
  222. n,3,6.25+5.1,2.025
  223. n,4,6.25+5.1+0.2+1.9/2,2.025
  224. n,5,6.25+5.1+2.3,2.025
  225. n,6,6.25+5.1+2.3+2.23,2.025
  226. nsym,x,10,2,6,1
  227. type,4
  228. r,1
  229. real,1
  230. mat,4
  231. secnum,7
  232. e,1,2
  233. secnum,6
  234. e,2,3
  235. secnum,5
  236. e,3,4
  237. e,4,5
  238. secnum,4
  239. e,5,6
  240. secnum,7
  241. e,1,12
  242. secnum,6
  243. e,12,13
  244. secnum,5
  245. e,13,14
  246. e,14,15
  247. secnum,4
  248. e,15,16
  249. egen,2,100,1,10,1,,,,,,0,4.7+1.25/2+1.5/2
  250. egen,17,100,11,20,1,,,,,,0,5.3+1.5
  251. egen,2,100,171,180,1,,,,,,0,4.7+1.5/2+1.25/2
  252. cm,hl,elem
  253. esel,none
  254. type,3
  255. r,2
  256. real,2
  257. mat,3
  258. sectype,10,beam,rect,zl
  259. secdata,0.4,1.2
  260. secnum,10
  261. e,1,101
  262. e,2,102
  263. e,3,103
  264. e,5,105
  265. e,12,112
  266. e,13,113
  267. e,15,115
  268. *do,i,0,15,1
  269. e,101+i*100,201+i*100
  270. e,102+i*100,202+i*100
  271. e,104+i*100,204+i*100
  272. e,112+i*100,212+i*100
  273. e,114+i*100,214+i*100
  274. *enddo
  275. e,1701,1801
  276. e,1702,1802
  277. e,1703,1803
  278. e,1705,1805
  279. e,1712,1812
  280. e,1713,1813
  281. e,1715,1815
  282. cm,zl,elem
  283. esel,none

  284. *dim,y,array,16
  285. *dim,cs,array,16
  286. *dim,z,array,16
  287. y(1)=8.1
  288. y(2)=14.9
  289. y(3)=21.7
  290. y(4)=28.5
  291. y(5)=35.3
  292. y(6)=42.1
  293. y(7)=48.9
  294. y(8)=55.7
  295. y(9)=62.5
  296. y(11)=21.7
  297. y(12)=28.5
  298. y(13)=42.1
  299. y(14)=48.9
  300. y(15)=59.1  
  301. y(16)=2.025
  302. f1=25
  303. m1=1.1
  304. k=log(m1+sqrt(m1*m1-1))
  305. *do,i,1,16,1
  306. cs(i)=1-y(i)*2/125
  307. z(i)=25-f1/(m1-1)*(cosh(k*cs(i))-1)
  308. *enddo
  309. *dim,ygd,array,23
  310. ygd(1) =0.0040004 +0.0129  
  311. ygd(2) =0.010685 +0.0222                        
  312. ygd(3) =0.019483 +0.0303  
  313. ygd(4) =0.029219 +0.0369  
  314. ygd(5) =0.039022 +0.0423  
  315. ygd(6) =0.048243 +0.0465  
  316. ygd(7) =0.05639 +0.0495         
  317. ygd(8) =0.06307 +0.0515         
  318. ygd(9) =0.068014 +0.0521  
  319. ygd(10) =0.071003 +0.0515  
  320. ygd(11) =0.071739 +0.0495  
  321. ygd(12) =0.069962 +0.0465  
  322. ygd(13) =0.065514 +0.423   
  323. ygd(14) =0.058289 +0.0369  
  324. ygd(15) =0.048248 +0.0304  
  325. ygd(16) =0.035507 +0.0224  
  326. ygd(17) =0.20415 +0.0123         
  327. ygd(18) =0.06577 +0.040234        
  328. ygd(19) =0.069769 +0.041215
  329. ngen,2,10000,104 ,114 ,10,0,0,z(1)+ygd(1 )
  330. ngen,2,10000,204 ,214 ,10,0,0,z(2)+ygd(2 )
  331. ngen,2,10000,304 ,314 ,10,0,0,z(3)+ygd(3 )
  332. ngen,2,10000,404 ,414 ,10,0,0,z(4)+ygd(4 )
  333. ngen,2,10000,504 ,514 ,10,0,0,z(5)+ygd(5 )
  334. ngen,2,10000,604 ,614 ,10,0,0,z(6)+ygd(6 )
  335. ngen,2,10000,704 ,714 ,10,0,0,z(7)+ygd(7 )
  336. ngen,2,10000,804 ,814 ,10,0,0,z(8)+ygd(8 )
  337. ngen,2,10000,904 ,914 ,10,0,0,z(9)+ygd(9 )
  338. ngen,2,10000,1004 ,1014 ,10,0,0,z(8)+ygd(10)
  339. ngen,2,10000,1104 ,1114 ,10,0,0,z(7)+ygd(11)
  340. ngen,2,10000,1204 ,1214 ,10,0,0,z(6)+ygd(12)
  341. ngen,2,10000,1304 ,1314 ,10,0,0,z(5)+ygd(13)
  342. ngen,2,10000,1404 ,1414 ,10,0,0,z(4)+ygd(14)
  343. ngen,2,10000,1504 ,1514 ,10,0,0,z(3)+ygd(15)
  344. ngen,2,10000,1604 ,1614 ,10,0,0,z(2)+ygd(16)
  345. ngen,2,10000,1704 ,1714 ,10,0,0,z(1)+ygd(17)
  346. ngen,2,3000,104 ,114 ,10,0,0,0
  347. ngen,2,3000,204 ,214 ,10,0,0,0
  348. ngen,2,3000,304 ,314 ,10,0,0,0
  349. ngen,2,3000,404 ,414 ,10,0,0,0
  350. ngen,2,3000,504 ,514 ,10,0,0,0
  351. ngen,2,3000,604 ,614 ,10,0,0,0
  352. ngen,2,3000,704 ,714 ,10,0,0,0
  353. ngen,2,3000,804 ,814 ,10,0,0,0
  354. ngen,2,3000,904 ,914 ,10,0,0,0
  355. ngen,2,3000,1004,1014 ,10,0,0,0
  356. ngen,2,3000,1104,1114 ,10,0,0,0
  357. ngen,2,3000,1204,1214 ,10,0,0,0
  358. ngen,2,3000,1304,1314 ,10,0,0,0
  359. ngen,2,3000,1404,1414 ,10,0,0,0
  360. ngen,2,3000,1504,1514 ,10,0,0,0
  361. ngen,2,3000,1604,1614 ,10,0,0,0
  362. ngen,2,3000,1704,1714 ,10,0,0,0
  363. ngen,2,6000,104 ,114 ,10,0,0,z(1)+ygd(1 )
  364. ngen,2,6000,204 ,214 ,10,0,0,z(2)+ygd(2 )
  365. ngen,2,6000,304 ,314 ,10,0,0,z(3)+ygd(3 )
  366. ngen,2,6000,404 ,414 ,10,0,0,z(4)+ygd(4 )
  367. ngen,2,6000,504 ,514 ,10,0,0,z(5)+ygd(5 )
  368. ngen,2,6000,604 ,614 ,10,0,0,z(6)+ygd(6 )
  369. ngen,2,6000,704 ,714 ,10,0,0,z(7)+ygd(7 )
  370. ngen,2,6000,804 ,814 ,10,0,0,z(8)+ygd(8 )
  371. ngen,2,6000,904 ,914 ,10,0,0,z(9)+ygd(9 )
  372. ngen,2,6000,1004,1014 ,10,0,0,z(8)+ygd(10)
  373. ngen,2,6000,1104,1114 ,10,0,0,z(7)+ygd(11)
  374. ngen,2,6000,1204,1214 ,10,0,0,z(6)+ygd(12)
  375. ngen,2,6000,1304,1314 ,10,0,0,z(5)+ygd(13)
  376. ngen,2,6000,1404,1414 ,10,0,0,z(4)+ygd(14)
  377. ngen,2,6000,1504,1514 ,10,0,0,z(3)+ygd(15)
  378. ngen,2,6000,1604,1614 ,10,0,0,z(2)+ygd(16)
  379. ngen,2,6000,1704,1714 ,10,0,0,z(1)+ygd(17)                                                            
  380. dgarea=109*3.1415926*0.25*0.007*0.007
  381. dgyl=0.3*1670*1e6
  382. dgl=dgyl*dgarea
  383. edg=1.95e11
  384. *dim,L0,array,17
  385. *dim,L,array,17
  386. *dim,deltL,array,17
  387. *dim,istrn,array,17
  388. *dim,fdg,array,17
  389. *do,i,1,9,1
  390. L0(i)=z(i)
  391. L0(18-i)=z(i)
  392. *enddo
  393. fdg(1 )=dgl
  394. fdg(2 )=dgl
  395. fdg(3 )=dgl
  396. fdg(4 )=dgl
  397. fdg(5 )=dgl
  398. fdg(6 )=dgl
  399. fdg(7 )=dgl
  400. fdg(8 )=dgl
  401. fdg(9 )=dgl
  402. fdg(10)=dgl
  403. fdg(11)=dgl
  404. fdg(12)=dgl
  405. fdg(13)=dgl
  406. fdg(14)=dgl
  407. fdg(15)=dgl
  408. fdg(16)=dgl
  409. fdg(17)=dgl
  410. type,6
  411. mat,6
  412. *do,i,1,17,1
  413. deltL(i)=fdg(i)*L0(i)/(edg*dgarea)
  414. L(i)=L0(i)+deltL(i)
  415. istrn(i)=deltL(i)/L(i)
  416. r,i+1000,dgarea,istrn(i)
  417. *enddo
  418. *do,i,1,17,1
  419. real,1+1000
  420. e,i*100+4+3000,i*100+4+6000
  421. e,i*100+10+4+3000,i*100+10+4+6000
  422. *enddo
  423. cm,dg,elem
  424. esel,none
  425. allsel
  426. real,1018
  427. n,20001, 0 ,21.7 ,z(11)        
  428. n,20002, 0 ,28.5 ,z(12)        
  429. n,20003, 0 ,42.1 ,z(13)        
  430. n,20004, 0 ,48.9 ,z(14)        
  431. n,20005, 0 ,59.1 ,z(15)        
  432. n,20006, 0 ,65.9 ,z(15)        
  433. n,20007, 0 ,76.1 ,z(14)        
  434. n,20008, 0 ,82.9 ,z(13)        
  435. n,20009, 0 ,96.5 ,z(12)        
  436. n,20010, 0 ,103.3 ,z(11)        
  437. n,20011, 12.5 ,21.7 ,z(11)        
  438. n,20012, 12.5 ,28.5 ,z(12)        
  439. n,20013, 12.5 ,42.1 ,z(13)        
  440. n,20014, 12.5 ,48.9 ,z(14)        
  441. n,20015, 12.5 ,59.1 ,z(15)+ygd(18)
  442. n,20016, 12.5 ,65.9 ,z(15)+ygd(19)
  443. n,20017, 12.5 ,76.1 ,z(14)        
  444. n,20018, 12.5 ,82.9 ,z(13)        
  445. n,20019, 12.5 ,96.5 ,z(12)        
  446. n,20020, 12.5 ,103.3 ,z(11)        
  447. n,20021, -12.5 ,21.7 ,z(11)        
  448. n,20022, -12.5 ,28.5 ,z(12)        
  449. n,20023, -12.5 ,42.1 ,z(13)        
  450. n,20024, -12.5 ,48.9 ,z(14)        
  451. n,20025, -12.5 ,59.1 ,z(15)+ygd(18)
  452. n,20026, -12.5 ,65.9 ,z(15)+ygd(19)
  453. n,20027, -12.5 ,76.1 ,z(14)        
  454. n,20028, -12.5 ,82.9 ,z(13)        
  455. n,20029, -12.5 ,96.5 ,z(12)        
  456. n,20030, -12.5 ,103.3 ,z(11)        
  457. n,20000, 0 ,62.5 ,z(15)
  458. type,5
  459. mat,5
  460. sectype,500,beam,ctube,windc
  461. secdata,0.486,0.50
  462. secnum,500
  463. e,20002,10404,90000
  464. e,20002,10414,90000
  465. e,20002,10304,90000
  466. e,20002,10314,90000
  467. e,20004,10704,90000
  468. e,20004,10714,90000
  469. e,20004,10604,90000
  470. e,20004,10614,90000
  471. e,20000,20015,90000
  472. e,20000,20016,90000
  473. e,20000,20025,90000
  474. e,20000,20026,90000
  475. e,20007,11104,90000
  476. e,20007,11114,90000
  477. e,20007,11204,90000
  478. e,20007,11214,90000
  479. e,20009,11404,90000
  480. e,20009,11414,90000
  481. e,20009,11504,90000
  482. e,20002,11514,90000
  483. cm,fc,elem
  484. esel,none
  485. n,30000,12.5,0,0
  486. n,30001,12.5,125,0
  487. n,30002,12.5,2.025,z(16)
  488. n,30003,12.5,125-2.025,z(16)
  489. n,40000,12.5,0,0
  490. n,40001,12.5,125,0
  491. n,40002,12.5,2.025,z(16)
  492. n,40003,12.5,125-2.025,z(16)
  493. *dim,zgzx,array,23
  494. *dim,ygzx,array,23
  495. zgzx(1 )=30000
  496. zgzx(2 )=30002
  497. zgzx(3 )=10104
  498. zgzx(4 )=10204
  499. zgzx(5 )=10304
  500. zgzx(6 )=10404
  501. zgzx(7 )=10504
  502. zgzx(8 )=10604
  503. zgzx(9 )=10704
  504. zgzx(10 )=10804
  505. zgzx(11 )=20015
  506. zgzx(12 )=10904
  507. zgzx(13 )=20016
  508. zgzx(14 )=11004
  509. zgzx(15 )=11104
  510. zgzx(16 )=11204
  511. zgzx(17 )=11304
  512. zgzx(18 )=11404
  513. zgzx(19 )=11504
  514. zgzx(20 )=11604
  515. zgzx(21 )=11704
  516. zgzx(22 )=30003
  517. zgzx(23 )=30001
  518. zgzx(1 )=40000
  519. zgzx(2 )=40002
  520. zgzx(3 )=10114
  521. zgzx(4 )=10214
  522. zgzx(5 )=10314
  523. zgzx(6 )=10414
  524. zgzx(7 )=10514
  525. zgzx(8 )=10614
  526. zgzx(9 )=10714
  527. zgzx(10 )=10814
  528. zgzx(11 )=20025
  529. zgzx(12 )=10914
  530. zgzx(13 )=20026
  531. zgzx(14 )=11014
  532. zgzx(15 )=11114
  533. zgzx(16 )=11214
  534. zgzx(17 )=11314
  535. zgzx(18 )=11414
  536. zgzx(19 )=11514
  537. zgzx(20 )=11614
  538. zgzx(21 )=11714
  539. zgzx(22 )=40003
  540. zgzx(23 )=40001
  541. n,90001,12.5,0,30
  542. n,90002,-12.5,0,30
  543. type,1
  544. mat,1
  545. r,2000
  546. secnum,1
  547. *do,i,1,22
  548. e,zgzx(i),zgzx(i+1),90001
  549. e,ygzx(i),ygzx(i+1),90002
  550. *enddo
  551. cm,gzgg,elem
  552. esel,none
  553. type,2
  554. mat,2
  555. r,2001
  556. secnum,2
  557. *do,i,1,22
  558. e,zgzx(i),zgzx(i+1),90001
  559. e,ygzx(i),ygzx(i+1),90002
  560. *enddo
  561. cm,gzhnt,elem
  562. esel,none
  563. type,9
  564. r,6
  565. real,6
  566. mat,9
  567. sectype,70,beam,rect,dizuo
  568. secdata,1.4,4
  569. secnum,70
  570. e,30001,30003,90001
  571. e,40001,40003,90002
  572. e,30000,30002,90001
  573. e,30000,40002,90002
  574. cm,gzdz,elem
  575. esel,none
  576. ngen,2,2000,1,16,1,0,-2.025,0
  577. ngen,2,100,1801,1816,1,0,2.025,0
  578. type,7
  579. mat,7
  580. r,7
  581. real,7
  582. sectype,60,beam,rect,gail
  583. secdata,4,1.6

  584. !!!!!!!!!!*****************************以上已经校核过!
复制代码
回复
分享到:

使用道具 举报

 楼主| 发表于 2006-5-22 05:11 | 显示全部楼层

回复:(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
发表于 2006-6-18 16:45 | 显示全部楼层
谢谢楼主,请问楼主有没有钢管相贯节点的命令流实例,谢谢
 楼主| 发表于 2006-6-19 03:42 | 显示全部楼层
命令流没有,gui得到是有一个,看帖子
http://forum.vibunion.com/forum/viewthread.php?tid=17479
发表于 2006-6-19 19:47 | 显示全部楼层
谢谢楼主
发表于 2008-2-25 20:57 | 显示全部楼层
谢谢!!非常
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-14 03:16 , Processed in 0.068640 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表