设为首页收藏本站联系我们
切换到宽版

声振论坛

 找回密码
 我要加入

QQ登录

只需一步，快速开始

声振论坛»论坛 › 工程软件 › Ansys讨论区 › 如何模拟钢管混凝土
返回列表 发新帖
查看: 2940|回复: 6

[材料与接触] 如何模拟钢管混凝土

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

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

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

x
如何模拟钢管和混凝土之间的接触，用什么单元比较好呢！
网格划分的时候应该注意些什么？
回复
分享到：

使用道具 举报

AaronSpark
发表于 2006-11-11 07:34 | 显示全部楼层
这个问题比较复杂，很多参数选择比较困难，除非能够做实验来获取，比如摩擦系数等

下面转两个例子你可以参考一下

[ 本帖最后由 AaronSpark 于 2006-11-11 07:37 编辑 ]

评分

1

查看全部评分

回复 支持 反对

使用道具 举报

AaronSpark
发表于 2006-11-11 07:36 | 显示全部楼层
《Algor、Ansys在桥梁工程中的应用方法与实例》中最后一个例子

  1. finish
  2. /clear
  3. /prep7
  4. /title,the analyse of steel arch bridge

  6. et,1,82
  7. 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. asel,all
  13. aptn,all
  14. adele,5,8,1
  15. allsel
  16. aadd,all
  17. smrtsize,5
  18. amesh,all
  19. secwrite,gg,sect,,1
  20. sectype,1,beam,mesh,
  21. secoffset,cent,,,
  22. secread,'gg','sect',,mesh
  23. asel,all
  24. aclear,all
  25. adele,all,,,1

  27. /replot

  29. 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. smrtsize,5
  36. amesh,all
  37. secwrite,hnt,sect,,1
  38. sectype,2,beam,mesh,
  39. secoffset,cent,,,
  40. secread,'hnt','sect',,mesh
  41. asel,all
  42. aclear,all
  43. adele,all,,,1

  45. /replot

  47. k,1,-0.5,0
  48. k,2,-0.75,0.25
  49. k,3,-0.75,0.5
  50. k,4,-0.55,0.5
  51. k,5,-0.55,0.85
  52. k,6,0.55,0.85
  53. k,7,0.55,0.5
  54. k,8,0.75,0.5
  55. k,9,0.75,0.25
  56. k,10,0.5,0
  57. a,1,2,3,4,5,6,7,8,9,10
  58. smrtsize,5
  59. amesh,all
  60. secwrite,hl1,sect,,1
  61. sectype,4,beam,mesh,
  62. secoffset,cent,,,
  63. secread,'hl1','sect',,mesh
  64. asel,all
  65. aclear,all
  66. adele,all,,,1

  68. /replot

  70. k,1,-0.5,0
  71. k,2,-0.5,1
  72. k,3,-0.75,1.25
  73. k,4,-0.75,1.5
  74. k,5,0.75,1.5
  75. k,6,0.75,1.25
  76. k,7,0.5,1
  77. k,8,0.5
  78. a,1,2,3,4,5,6,7,8
  79. smrtsize,5
  80. amesh,all
  81. secwrite,hl2,sect,,1
  82. sectype,5,beam,mesh,
  83. secoffset,cent,,,
  84. secread,'hl2','sect',,mesh
  85. asel,all
  86. aclear,all
  87. adele,all,,,1

  89. /replot

  91. k,1,-0.5,0
  92. k,2,-0.5,1.04
  93. k,3,-0.75,1.04+0.25
  94. k,4,-0.75,1.04+0.5
  95. k,5,-0.55,1.04+0.5
  96. k,6,-0.55,1.04+0.5+0.35
  97. k,7,0.55,1.04+0.5+0.35
  98. k,8,0.55,1.04+0.5
  99. k,9,0.75,1.04+0.5
  100. k,10,0.75,1.04+0.25
  101. k,11,0.5,1.04
  102. k,12,0.5,0
  103. k,13,-0.2,0.58
  104. k,14,-0.2,1.04+0.5-0.3
  105. k,15,0.2,1.04+0.5-0.3
  106. k,16,0.2,0.58
  107. a,1,2,3,4,5,6,7,8,9,10,11,12
  108. a,13,14,15,16
  109. asba,1,2
  110. smrtsize,5
  111. amesh,all
  112. secwrite,hl3,sect,,1
  113. sectype,6,beam,mesh,
  114. secoffset,cent,,,
  115. secread,'hl3','sect',,mesh
  116. asel,all
  117. aclear,all
  118. adele,all,,,1

  120. /replot  
  121. ldele,all
  122. kdele,all

  124. k,1,-0.5,0
  125. k,2,-0.5,1.23                     
  126. k,3,-0.75,1.23+0.25               
  127. k,4,-0.75,1.23+0.5               
  128. k,5,-0.55,1.23+0.5               
  129. k,6,-0.55,1.23+0.5+0.35           
  130. k,7,0.55,1.23+0.5+0.35            
  131. k,8,0.55,1.23+0.5                 
  132. k,9,0.75,1.23+0.5                 
  133. k,10,0.75,1.23+0.25               
  134. k,11,0.5,1.23                     
  135. k,12,0.5,0                        
  136. k,13,-0.2,0.38                    
  137. k,14,-0.3,0.38+0.2            
  138. k,15,-0.3,1.23+0.5-0.3            
  139. k,16,-0.2,1.23+0.5-0.1
  140. k,17,0.2,1.23+0.5-0.1
  141. k,18,0.3,1.23+0.5-0.1-0.2
  142. k,19,0.3,0.38+0.2
  143. k,20,0.2,0.38
  144.                   
  145. a,1,2,3,4,5,6,7,8,9,10,11,12      
  146. a,13,14,15,16,17,18,19,20                     
  147. asba,1,2                          
  148. smrtsize,5                        
  149. amesh,all                        
  150. secwrite,hl4,sect,,1              
  151. sectype,7,beam,mesh,              
  152. secoffset,cent,,,                 
  153. secread,'hl4','sect',,mesh        
  154. asel,all                          
  155. aclear,all                        
  156. adele,all,,,1                     
  157.                                  
  158. /replot  

  160. k,1,-0.5,0
  161. k,2,-0.5,1.23                     
  162. k,3,-0.75,1.23+0.25               
  163. k,4,-0.75,1.23+0.5               
  164. k,5,-0.55,1.23+0.5               
  165. k,6,-0.55,1.23+0.5+0.35           
  166. k,7,0.55,1.23+0.5+0.35            
  167. k,8,0.55,1.23+0.5                 
  168. k,9,0.75,1.23+0.5                 
  169. k,10,0.75,1.23+0.25               
  170. k,11,0.5,1.23                     
  171. k,12,0.5,0                        
  172. a,1,2,3,4,5,6,7,8,9,10,11,12
  173. smrtsize,5                        
  174. amesh,all                        
  175. secwrite,hl5,sect,,1              
  176. sectype,8,beam,mesh,              
  177. secoffset,cent,,,                 
  178. secread,'hl5','sect',,mesh        
  179. asel,all                          
  180. aclear,all                        
  181. adele,all,,,1                     
  182.                                  
  183. /replot     

  185. etdele,1
  186. /replot

  188. et,1,beam44
  189. mp,ex,1,2.1e11
  190. mp,dens,1,7800
  191. mp,prxy,1,0.3
  192. n,90000,0,0,30

  194. et,2,beam44
  195. mp,ex,2,3.5e10
  196. mp,dens,2,2600
  197. mp,prxy,2,0.1667

  199. et,3,beam44
  200. mp,ex,3,3.0e10
  201. mp,dens,3,2600
  202. mp,prxy,3,0.1667

  204. et,4,beam44
  205. mp,ex,4,3.0e10
  206. mp,dens,4,2600
  207. mp,prxy,4,0.1667

  209. et,5,beam44
  210. mp,ex,5,2.1e11
  211. mp,dens,5,7800
  212. mp,prxy,5,0.3

  214. et,6,link10
  215. mp,ex,6,1.9e11
  216. mp,dens,6,7800
  217. mp,prxy,6,0.3
  218. keyopt,6,3,0

  220. et,7,beam44
  221. mp,ex,7,3.0e10
  222. mp,dens,7,2600
  223. mp,prxy,7,0.1667

  225. et,8,beam44
  226. mp,ex,8,3.3e10
  227. mp,dens,8,2600
  228. mp,prxy,8,0.1667

  230. et,9,beam44
  231. mp,ex,9,3.0e10
  232. mp,dens,9,2600
  233. mp,prxy,9,0.1667

  235. et,10,link10
  236. mp,ex,10,2.0e11
  237. mp,dens,10,7800
  238. mp,prxy,10,0.3
  239. keyopt,10,3,0

  241. et,11,beam44
  242. mp,ex,11,2.85e10
  243. mp,dens,11,2500
  244. mp,prxy,11,0.1667
  245. n,90000,0,0,30


  248. n,1,0,2.025
  249. n,2,6.25,2.025
  250. n,3,6.25+5.1,2.025
  251. n,4,6.25+5.1+0.2+1.9/2,2.025
  252. n,5,6.25+5.1+2.3,2.025
  253. n,6,6.25+5.1+2.3+2.23,2.025
  254. nsym,x,10,2,6,1
  255. type,4
  256. r,1
  257. real,1
  258. mat,4
  259. secnum,7
  260. e,1,2
  261. secnum,6
  262. e,2,3
  263. secnum,5
  264. e,3,4
  265. e,4,5
  266. secnum,4
  267. e,5,6
  268. secnum,7
  269. e,1,12
  270. secnum,6
  271. e,12,13
  272. secnum,5
  273. e,13,14
  274. e,14,15
  275. secnum,4
  276. e,15,16

  278. egen,2,100,1,10,1,,,,,,0,4.7+1.25/2+1.5/2
  279. egen,17,100,11,20,1,,,,,,0,5.3+1.5
  280. egen,2,100,171,180,1,,,,,,0,4.7+1.5/2+1.25/2
  281. cm,hl,elem
  282. esel,none

  284. type,3
  285. r,2
  286. real,2
  287. mat,3
  288. sectype,10,beam,rect,zl
  289. secdata,0.4,1.2
  290. secnum,10

  292. e,1,101
  293. e,2,102
  294. e,3,103
  295. e,5,105
  296. e,12,112
  297. e,13,113
  298. e,15,115

  300. *do,i,0,15,1
  301. e,101+i*100,201+i*100
  302. e,102+i*100,202+i*100
  303. e,104+i*100,204+i*100
  304. e,112+i*100,212+i*100
  305. e,114+i*100,214+i*100
  306. *enddo

  308. e,1701,1801
  309. e,1702,1802
  310. e,1703,1803
  311. e,1705,1805
  312. e,1712,1812
  313. e,1713,1813
  314. e,1715,1815
  315. cm,zl,elem
  316. esel,none


  319. *dim,y,array,16
  320. *dim,cs,array,16
  321. *dim,z,array,16

  323. y(1)=8.1
  324. y(2)=14.9
  325. y(3)=21.7
  326. y(4)=28.5
  327. y(5)=35.3
  328. y(6)=42.1
  329. y(7)=48.9
  330. y(8)=55.7
  331. y(9)=62.5

  333. y(11)=21.7
  334. y(12)=28.5
  335. y(13)=42.1
  336. y(14)=48.9
  337. y(15)=59.1  

  339. y(16)=2.025
  340. f1=25
  341. m1=1.1
  342. k=log(m1+sqrt(m1*m1-1))
  343. *do,i,1,16,1
  344. cs(i)=1-y(i)*2/125
  345. z(i)=25-f1/(m1-1)*(cosh(k*cs(i))-1)
  346. *enddo

  348. *dim,ygd,array,23

  350. ygd(1)        =0.0040004        +0.0129  
  351. ygd(2)        =0.010685        +0.0222                        
  352. ygd(3)        =0.019483        +0.0303  
  353. ygd(4)        =0.029219        +0.0369  
  354. ygd(5)        =0.039022        +0.0423  
  355. ygd(6)        =0.048243        +0.0465  
  356. ygd(7)        =0.05639        +0.0495         
  357. ygd(8)        =0.06307        +0.0515         
  358. ygd(9)        =0.068014        +0.0521  
  359. ygd(10)        =0.071003        +0.0515  
  360. ygd(11)        =0.071739        +0.0495  
  361. ygd(12)        =0.069962        +0.0465  
  362. ygd(13)        =0.065514        +0.423   
  363. ygd(14)        =0.058289        +0.0369  
  364. ygd(15)        =0.048248        +0.0304  
  365. ygd(16)        =0.035507        +0.0224  
  366. ygd(17)        =0.20415        +0.0123         
  367. ygd(18)        =0.06577        +0.040234        
  368. ygd(19)        =0.069769        +0.041215

  370. ngen,2,10000,104        ,114        ,10,0,0,z(1)+ygd(1 )
  371. ngen,2,10000,204        ,214        ,10,0,0,z(2)+ygd(2 )
  372. ngen,2,10000,304        ,314        ,10,0,0,z(3)+ygd(3 )
  373. ngen,2,10000,404        ,414        ,10,0,0,z(4)+ygd(4 )
  374. ngen,2,10000,504        ,514        ,10,0,0,z(5)+ygd(5 )
  375. ngen,2,10000,604        ,614        ,10,0,0,z(6)+ygd(6 )
  376. ngen,2,10000,704        ,714        ,10,0,0,z(7)+ygd(7 )
  377. ngen,2,10000,804        ,814        ,10,0,0,z(8)+ygd(8 )
  378. ngen,2,10000,904        ,914        ,10,0,0,z(9)+ygd(9 )
  379. ngen,2,10000,1004        ,1014        ,10,0,0,z(8)+ygd(10)
  380. ngen,2,10000,1104        ,1114        ,10,0,0,z(7)+ygd(11)
  381. ngen,2,10000,1204        ,1214        ,10,0,0,z(6)+ygd(12)
  382. ngen,2,10000,1304        ,1314        ,10,0,0,z(5)+ygd(13)
  383. ngen,2,10000,1404        ,1414        ,10,0,0,z(4)+ygd(14)
  384. ngen,2,10000,1504        ,1514        ,10,0,0,z(3)+ygd(15)
  385. ngen,2,10000,1604        ,1614        ,10,0,0,z(2)+ygd(16)
  386. ngen,2,10000,1704        ,1714        ,10,0,0,z(1)+ygd(17)

  388. ngen,2,3000,104        ,114        ,10,0,0,0
  389. ngen,2,3000,204        ,214        ,10,0,0,0
  390. ngen,2,3000,304        ,314        ,10,0,0,0
  391. ngen,2,3000,404        ,414        ,10,0,0,0
  392. ngen,2,3000,504        ,514        ,10,0,0,0
  393. ngen,2,3000,604        ,614        ,10,0,0,0
  394. ngen,2,3000,704        ,714        ,10,0,0,0
  395. ngen,2,3000,804        ,814        ,10,0,0,0
  396. ngen,2,3000,904        ,914        ,10,0,0,0
  397. ngen,2,3000,1004,1014        ,10,0,0,0
  398. ngen,2,3000,1104,1114        ,10,0,0,0
  399. ngen,2,3000,1204,1214        ,10,0,0,0
  400. ngen,2,3000,1304,1314        ,10,0,0,0
  401. ngen,2,3000,1404,1414        ,10,0,0,0
  402. ngen,2,3000,1504,1514        ,10,0,0,0
  403. ngen,2,3000,1604,1614        ,10,0,0,0
  404. ngen,2,3000,1704,1714        ,10,0,0,0

  406. ngen,2,6000,104        ,114        ,10,0,0,z(1)+ygd(1 )
  407. ngen,2,6000,204        ,214        ,10,0,0,z(2)+ygd(2 )
  408. ngen,2,6000,304        ,314        ,10,0,0,z(3)+ygd(3 )
  409. ngen,2,6000,404        ,414        ,10,0,0,z(4)+ygd(4 )
  410. ngen,2,6000,504        ,514        ,10,0,0,z(5)+ygd(5 )
  411. ngen,2,6000,604        ,614        ,10,0,0,z(6)+ygd(6 )
  412. ngen,2,6000,704        ,714        ,10,0,0,z(7)+ygd(7 )
  413. ngen,2,6000,804        ,814        ,10,0,0,z(8)+ygd(8 )
  414. ngen,2,6000,904        ,914        ,10,0,0,z(9)+ygd(9 )
  415. ngen,2,6000,1004,1014        ,10,0,0,z(8)+ygd(10)
  416. ngen,2,6000,1104,1114        ,10,0,0,z(7)+ygd(11)
  417. ngen,2,6000,1204,1214        ,10,0,0,z(6)+ygd(12)
  418. ngen,2,6000,1304,1314        ,10,0,0,z(5)+ygd(13)
  419. ngen,2,6000,1404,1414        ,10,0,0,z(4)+ygd(14)
  420. ngen,2,6000,1504,1514        ,10,0,0,z(3)+ygd(15)
  421. ngen,2,6000,1604,1614        ,10,0,0,z(2)+ygd(16)
  422. ngen,2,6000,1704,1714        ,10,0,0,z(1)+ygd(17)                                                            

  424. dgarea=109*3.1415926*0.25*0.007*0.007
  425. dgyl=0.3*1670*1e6

  427. dgl=dgyl*dgarea
  428. edg=1.95e11
  429. *dim,L0,array,17
  430. *dim,L,array,17
  431. *dim,deltL,array,17
  432. *dim,istrn,array,17
  433. *dim,fdg,array,17
  434. *do,i,1,9,1
  435. L0(i)=z(i)
  436. L0(18-i)=z(i)
  437. *enddo

  439. fdg(1 )=dgl
  440. fdg(2 )=dgl
  441. fdg(3 )=dgl
  442. fdg(4 )=dgl
  443. fdg(5 )=dgl
  444. fdg(6 )=dgl
  445. fdg(7 )=dgl
  446. fdg(8 )=dgl
  447. fdg(9 )=dgl
  448. fdg(10)=dgl
  449. fdg(11)=dgl
  450. fdg(12)=dgl
  451. fdg(13)=dgl
  452. fdg(14)=dgl
  453. fdg(15)=dgl
  454. fdg(16)=dgl
  455. fdg(17)=dgl

  457. type,6
  458. mat,6
  459. *do,i,1,17,1
  460. deltL(i)=fdg(i)*L0(i)/(edg*dgarea)
  461. L(i)=L0(i)+deltL(i)
  462. istrn(i)=deltL(i)/L(i)
  463. r,i+1000,dgarea,istrn(i)
  464. *enddo

  466. *do,i,1,17,1
  467. real,1+1000
  468. e,i*100+4+3000,i*100+4+6000
  469. e,i*100+10+4+3000,i*100+10+4+6000
  470. *enddo
  471. cm,dg,elem
  472. esel,none
  473. allsel

  475. real,1018

  477. n,20001,        0        ,21.7        ,z(11)               
  478. n,20002,        0        ,28.5        ,z(12)               
  479. n,20003,        0        ,42.1        ,z(13)               
  480. n,20004,        0        ,48.9        ,z(14)               
  481. n,20005,        0        ,59.1        ,z(15)               
  482. n,20006,        0        ,65.9        ,z(15)               
  483. n,20007,        0        ,76.1        ,z(14)               
  484. n,20008,        0        ,82.9        ,z(13)               
  485. n,20009,        0        ,96.5        ,z(12)               
  486. n,20010,        0        ,103.3        ,z(11)               
  487. n,20011,        12.5        ,21.7        ,z(11)               
  488. n,20012,        12.5        ,28.5        ,z(12)               
  489. n,20013,        12.5        ,42.1        ,z(13)               
  490. n,20014,        12.5        ,48.9        ,z(14)               
  491. n,20015,        12.5        ,59.1        ,z(15)+ygd(18)
  492. n,20016,        12.5        ,65.9        ,z(15)+ygd(19)
  493. n,20017,        12.5        ,76.1        ,z(14)               
  494. n,20018,        12.5        ,82.9        ,z(13)               
  495. n,20019,        12.5        ,96.5        ,z(12)               
  496. n,20020,        12.5        ,103.3        ,z(11)               
  497. n,20021,        -12.5        ,21.7        ,z(11)               
  498. n,20022,        -12.5        ,28.5        ,z(12)               
  499. n,20023,        -12.5        ,42.1        ,z(13)               
  500. n,20024,        -12.5        ,48.9        ,z(14)               
  501. n,20025,        -12.5        ,59.1        ,z(15)+ygd(18)
  502. n,20026,        -12.5        ,65.9        ,z(15)+ygd(19)
  503. n,20027,        -12.5        ,76.1        ,z(14)               
  504. n,20028,        -12.5        ,82.9        ,z(13)               
  505. n,20029,        -12.5        ,96.5        ,z(12)               
  506. n,20030,        -12.5        ,103.3        ,z(11)               
  507. n,20000,        0        ,62.5        ,z(15)

  509. type,5
  510. mat,5
  511. sectype,500,beam,ctube,windc
  512. secdata,0.486,0.50

  514. secnum,500

  516. e,20002,10404,90000
  517. e,20002,10414,90000
  518. e,20002,10304,90000
  519. e,20002,10314,90000

  521. e,20004,10704,90000
  522. e,20004,10714,90000
  523. e,20004,10604,90000
  524. e,20004,10614,90000

  526. e,20000,20015,90000
  527. e,20000,20016,90000
  528. e,20000,20025,90000
  529. e,20000,20026,90000

  531. e,20007,11104,90000
  532. e,20007,11114,90000
  533. e,20007,11204,90000
  534. e,20007,11214,90000

  536. e,20009,11404,90000
  537. e,20009,11414,90000
  538. e,20009,11504,90000
  539. e,20002,11514,90000
  540. cm,fc,elem
  541. esel,none

  543. n,30000,12.5,0,0
  544. n,30001,12.5,125,0
  545. n,30002,12.5,2.025,z(16)
  546. n,30003,12.5,125-2.025,z(16)

  548. n,40000,12.5,0,0
  549. n,40001,12.5,125,0
  550. n,40002,12.5,2.025,z(16)
  551. n,40003,12.5,125-2.025,z(16)

  553. *dim,zgzx,array,23
  554. *dim,ygzx,array,23
  555. zgzx(1        )=30000
  556. zgzx(2        )=30002
  557. zgzx(3        )=10104
  558. zgzx(4        )=10204
  559. zgzx(5        )=10304
  560. zgzx(6        )=10404
  561. zgzx(7        )=10504
  562. zgzx(8        )=10604
  563. zgzx(9        )=10704
  564. zgzx(10        )=10804
  565. zgzx(11        )=20015
  566. zgzx(12        )=10904
  567. zgzx(13        )=20016
  568. zgzx(14        )=11004
  569. zgzx(15        )=11104
  570. zgzx(16        )=11204
  571. zgzx(17        )=11304
  572. zgzx(18        )=11404
  573. zgzx(19        )=11504
  574. zgzx(20        )=11604
  575. zgzx(21        )=11704
  576. zgzx(22        )=30003
  577. zgzx(23        )=30001
  578. zgzx(1        )=40000
  579. zgzx(2        )=40002
  580. zgzx(3        )=10114
  581. zgzx(4        )=10214
  582. zgzx(5        )=10314
  583. zgzx(6        )=10414
  584. zgzx(7        )=10514
  585. zgzx(8        )=10614
  586. zgzx(9        )=10714
  587. zgzx(10        )=10814
  588. zgzx(11        )=20025
  589. zgzx(12        )=10914
  590. zgzx(13        )=20026
  591. zgzx(14        )=11014
  592. zgzx(15        )=11114
  593. zgzx(16        )=11214
  594. zgzx(17        )=11314
  595. zgzx(18        )=11414
  596. zgzx(19        )=11514
  597. zgzx(20        )=11614
  598. zgzx(21        )=11714
  599. zgzx(22        )=40003
  600. zgzx(23        )=40001
  601. n,90001,12.5,0,30
  602. n,90002,-12.5,0,30
  603. type,1
  604. mat,1
  605. r,2000
  606. secnum,1
  607. *do,i,1,22
  608. e,zgzx(i),zgzx(i+1),90001
  609. e,ygzx(i),ygzx(i+1),90002
  610. *enddo
  611. cm,gzgg,elem
  612. esel,none
  613. type,2
  614. mat,2
  615. r,2001
  616. secnum,2
  617. *do,i,1,22
  618. e,zgzx(i),zgzx(i+1),90001
  619. e,ygzx(i),ygzx(i+1),90002
  620. *enddo
  621. cm,gzhnt,elem
  622. esel,none

  624. type,9
  625. r,6
  626. real,6
  627. mat,9
  628. sectype,70,beam,rect,dizuo
  629. secdata,1.4,4
  630. secnum,70
  631. e,30001,30003,90001
  632. e,40001,40003,90002
  633. e,30000,30002,90001
  634. e,30000,40002,90002
  635. cm,gzdz,elem
  636. esel,none


  639. ngen,2,2000,1,16,1,0,-2.025,0
  640. ngen,2,100,1801,1816,1,0,2.025,0
  641. type,7
  642. mat,7
  643. r,7
  644. real,7
  645. sectype,60,beam,rect,gail
  646. secdata,4,1.6


  649. !!!!!!!!!!*****************************以上已经校核过!

  651. secnum,60
  652. n,90005,0,0,20
  653. n,90006,0,125,20
  654. e,1901,        1902, 90006
  655. e,1902,        1903, 90006
  656. e,1903,        1904, 90006
  657. e,1904,        1905, 90006
  658. e,1905,        1906, 90006
  659. e,1906,        1907, 90006
  660. e,1907,        1908, 90006
  661. e,1908,        1909, 90006
  662. e,1909,        1910, 90006
  663. e,1910,        1911, 90006
  664. e,1911,        1912, 90006
  665. e,1912,        1913, 90006
  666. e,1913,        1914, 90006
  667. e,1914,        1915, 90006
  668. e,1915,        1916, 90006
  669. e,2001,        2002, 90005
  670. e,2002,        2003, 90005
  671. e,2003,        2004, 90005
  672. e,2004,        2005, 90005
  673. e,2005,        2006, 90005
  674. e,2006,        2007, 90005
  675. e,2007,        2008, 90005
  676. e,2008,        2009, 90005
  677. e,2009,        2010, 90005
  678. e,2010,        2011, 90005
  679. e,2011,        2012, 90005
  680. e,2012,        2013, 90005
  681. e,2013,        2014, 90005
  682. e,2014,        2015, 90005
  683. e,2015,        2016, 90005
  684. cm,gailiang,elem
  685. esel,none

  687. n,90003,14.5,0,0
  688. n,90004,-14,5,0,0
  689. type,8
  690. r,8
  691. real,8
  692. sectype,50,beam,hrec,dun
  693. secdata,3,4,0.8,0.8,0.8,0.8
  694. ngen,2,10000,2004,2014,10,0,0,-10
  695. ngen,2,10000,1904,1914,10,0,0,-10

  697. secnum,50
  698. e,2004,12004,90003
  699. e,2014,12014,90004
  700. e,1904,11904,90003
  701. e,1914,11914,90004
  702. cm,dunzhu,elem
  703. esel,none

  705. type,11
  706. mat,11
  707. r,1
  708. sectype,40,beam,rect,chengtai
  709. secdata,12.7,9
  710. r,11
  711. real,11
  712. secnum,40
  713. ngen,2,20000,2004,2014,10,0,0,-11.5
  714. ngen,2,20000,1904,1914,10,0,0,-11.5
  715. ngen,2,30000,2004,2014,10,0,0,-13
  716. ngen,2,30000,1904,1914,10,0,0,-13  

  718. e,12004,22004,90003
  719. e,11904,21904,90003
  720. e,12014,22014,90004
  721. e,11904,21904,90004

  723. e,22004,32004,90003
  724. e,21904,31904,90003
  725. e,22014,32014,90004
  726. e,21904,31904,90004
  727. cm,chtai,elem
  728. esel,none

  730. type,9
  731. mat,9
  732. r,9
  733. real,9
  734. n,90007,12.5,0,10
  735. n,90008,-12.5,0,10
  736. cm,hxl,elem
  737. esel,none


  740. n,50000,6.15+1.85,2.7,-13
  741. n,50001,6.15+1.85,-2.7,-13
  742. n,50002,6.15+1.85+4.5,2.7,-13
  743. n,50003,6.15+1.85+4.5,-2.7,-13
  744. n,50004,6.15+1.27-1.85,2.7,-13
  745. n,50005,6.15+1.27-1.85,-2.7,-13
  746. nsel,s,,,50000,50005,1
  747. nsym,x,10,50000,50005,1
  748. allsel
  749. ngen,2,100,50000,50015,1,0,125,-2.5
  750. ngen,10,1000,50000,50015,1,0,0,-2.5
  751. esel,none
  752. type,11
  753. mat,11
  754. r,11
  755. real,11
  756. sectype,20,beam,csolid,pier
  757. secdata,0.9
  758. secnum,20
  759. *do,i,0,8,1
  760. *do,k,0,5,1
  761. e,50000+k+i*1000,50000+k+(i+1)*1000,90003
  762. e,50000+k+10+i*1000,50000+k+10+(i+1)*1000,90003
  763. e,50000+k+100+i*1000,50000+k+100+(i+1)*1000,90003
  764. e,50000+k+110+i*1000,50000+k+110+(i+1)*1000,90003
  765. *enddo
  766. *enddo
  767. cm,zhuang,elem
  768. esel,none

  770. d=1.8
  771. kf=0.9
  772. k0=1+1/d
  773. b1=0.6

  775. L1=3.6
  776. h1=3*(d+1)

  778. k=b1+(1-b1)/0.6*L1/h1
  779. B=Kf*k0*k*d


  782. *dim,m,array,5
  783. m(1)=8
  784. m(2)=15
  785. m(3)=20
  786. m(4)=25
  787. m(5)=4000
  788. *dim,h,array,5
  789. h(1)=5
  790. h(2)=5
  791. h(3)=5
  792. h(4)=5
  793. h(5)=10
  794. *dim,c,array,4
  795. c(1)=m(1)*h(1)
  796. c(2)=m(2)*(h(1)+h(2))
  797. c(3)=m(3)*(h(1)+h(2)+h(3))
  798. c(4)=m(4)*(h(1)+h(2)+h(3)+h(4))
  799. *dim,e,array,4
  800. e(1)=0.5*c(1)*h(1)*B*1e6
  801. *do,i,2,4,1
  802. e(i)=0.5*(c(i-1)+c(i))*h(1)*B*1e6
  803. *enddo
  804. et,20,combin14
  805. r,20,e(1)
  806. r,30,e(2)
  807. r,40,e(3)
  808. r,50,e(4)
  809. r,60,m(5)*3.1415926*B*B/4*1e6

  811. ngen,2,10000,51000,51115,1,1,0,0
  812. ngen,2,20000,51000,51115,1,0,1,0
  813. ngen,2,30000,51000,51115,1,0,0,1
  814. type,20
  815. real,20
  816. *do,i,0,5,1
  817. e,51000+i,51000+i+10000
  818. e,51000+i+10,51000+i+10+10000
  819. e,51100+i,51100+i+10000
  820. e,51100+i+10,51100+i+10+10000
  821. e,51000+i,51000+i+20000
  822. e,51000+i+10,51000+i+10+20000
  823. e,51100+i,51100+i+20000
  824. e,51100+i+10,51100+i+10+20000
  825. *enddo
  826. cm,spring1,elem
  827. esel,none

  829. ngen,2,10000,53000,53115,1,1,0,0
  830. ngen,2,20000,53000,53115,1,0,1,0
  831. ngen,2,30000,53000,53115,1,0,0,1

  833. type,20
  834. real,30
  835. *do,i,0,5,1
  836. e,53000+i,53000+i+10000
  837. e,53000+i+10,53000+i+10+10000
  838. e,53100+i,53100+i+10000
  839. e,53100+i+10,53100+i+10+10000
  840. e,53000+i,53000+i+20000
  841. e,53000+i+10,53000+i+10+20000
  842. e,53100+i,53100+i+20000
  843. e,53100+i+10,53100+i+10+20000
  844. *enddo
  845. cm,spring2,elem
  846. esel,none


  849. ngen,2,10000,55000,55115,1,1,0,0  
  850. ngen,2,20000,55000,55115,1,0,1,0  
  851.                                  
  852. type,20                           
  853. real,40                           
  854. *do,i,0,5,1                       
  855. e,55000+i,55000+i+10000           
  856. e,55000+i+10,55000+i+10+10000     
  857. e,55100+i,55100+i+10000           
  858. e,55100+i+10,55100+i+10+10000     
  859. e,55000+i,55000+i+20000           
  860. e,55000+i+10,55000+i+10+20000     
  861. e,55100+i,55100+i+20000           
  862. e,55100+i+10,55100+i+10+20000     
  863. *enddo                           
  864. cm,spring3,elem                  
  865. esel,none                        

  867. ngen,2,10000,57000,57115,1,1,0,0     
  868. ngen,2,20000,57000,57115,1,0,1,0     
  869.                                     
复制代码


未完......
回复 支持 反对

使用道具 举报

AaronSpark
发表于 2006-11-11 07:37 | 显示全部楼层
续上....

  1. type,20                              
  2. real,50                              
  3. *do,i,0,5,1                          
  4. e,57000+i,57000+i+10000              
  5. e,57000+i+10,57000+i+10+10000        
  6. e,57100+i,57100+i+10000              
  7. e,57100+i+10,57100+i+10+10000        
  8. e,57000+i,57000+i+20000              
  9. e,57000+i+10,57000+i+10+20000        
  10. e,57100+i,57100+i+20000              
  11. e,57100+i+10,57100+i+10+20000        
  12. *enddo                              
  13. cm,spring4,elem                     
  14. esel,none

  16. ngen,2,10000,58000,58115,1,1,0,0                                
  17. ngen,2,20000,58000,58115,1,0,1,0     
  18.                                     
  19. type,20                              
  20. real,60                              
  21. *do,i,0,5,1                          
  22. e,58000+i,58000+i+10000              
  23. e,58000+i+10,58000+i+10+10000        
  24. e,58100+i,58100+i+10000              
  25. e,58100+i+10,58100+i+10+10000        
  26. e,58000+i,58000+i+20000              
  27. e,58000+i+10,58000+i+10+20000        
  28. e,58100+i,58100+i+20000              
  29. e,58100+i+10,58100+i+10+20000        
  30. *enddo                              
  31. cm,spring5,elem                     
  32. esel,none

  34. ngen,2,10000,59000,59115,1,1,0,0                                 
  35. ngen,2,10000,59000,59115,1,1,0,0
  36. ngen,2,20000,59000,59115,1,0,1,0
  37. ngen,2,20000,59000,59115,1,0,1,0
  38. ngen,2,30000,59000,59115,1,0,0,-1
  39. ngen,2,30000,59000,59115,1,0,0,-1
  40.    
  41.                                     
  42. type,20                              
  43. real,60                              
  44. *do,i,0,5,1                          
  45. e,59000+i,59000+i+10000              
  46. e,59000+i+10,59000+i+10+10000        
  47. e,59100+i,59100+i+10000              
  48. e,59100+i+10,59100+i+10+10000        
  49. e,59000+i,59000+i+20000              
  50. e,59000+i+10,59000+i+10+20000        
  51. e,59100+i,59100+i+20000              
  52. e,59100+i+10,59100+i+10+20000
  53. e,59000+i,59000+i+30000      
  54. e,59000+i+10,59000+i+10+30000
  55. e,59100+i,59100+i+30000     
  56. e,59100+i+10,59100+i+10+30000
  57.       
  58. *enddo                              
  59. cm,spring6,elem                     
  60. allsel


  63. allsel
  64. nsel,s,loc,z,-50,-0.01
  65. nsel,r,loc,y,-10,10
  66. cm,zdnode,node
  67. esln,s,,all

  69. cm,zdelem,elem
  70. allsel
  71. nsel,u,,,zdnode
  72. nsym,y,100000,all
  73. allsel
  74. esel,u,,,zdelem
  75. esym,,100000,all
  76. allsel


  79. xgarea=8*27*7*3.1415926*0.25*0.005*0.005
  80. exg=1.95e11
  81. *dim,fxg,array,2
  82. *dim,Llast,array,2
  83. *dim,Listrn,array,2
  84. *dim,detL,array,2
  85. xgyl=0.3*1860*1e6
  86. xgl=xgyl*xgarea
  87. fxg(1)=xgl
  88. fxg(2)=xgl
  89. Lzero=250
  90. *do,i,1,2,1
  91. detL(i)=Lzero*Fxg(i)/(exg*xgarea)
  92. Llast(i)=Lzero+detL(i)
  93. Listrn(i)=detl(i)/Llast(i)
  94. r,i+100,xgarea,Listrn(i)
  95. *enddo
  96. type,10
  97. mat,10
  98. real,101
  99. e,30001,130001
  100. real,102
  101. e,40001,140001
  102. cm,xigan,elem
  103. esel,none
  104. allesl

  106. cp,        1,        ux        ,2004        ,30000
  107. cp,        1,        uy        ,2004        ,30000                     
  108. cp,        next,        uz        ,2004        ,30000
  109. cp,        next,        roty        ,2004        ,30000
  110. cp,        next,        rotz        ,2004        ,30000
  111. cp,        next,        ux        ,1904        ,30001
  112. cp,        next,        uy        ,1904        ,30001
  113. cp,        next,        uz        ,1904        ,30001
  114. cp,        next,        roty        ,1904        ,30001
  115. cp,        next,        rotz        ,1904        ,30001
  116. cp,        next,        ux        ,2014        ,40000
  117. cp,        next,        uy        ,2014        ,40000
  118. cp,        next,        uz        ,2014        ,40000
  119. cp,        next,        roty        ,2014        ,40000
  120. cp,        next,        rotz        ,2014        ,40000
  121. cp,        next,        ux        ,1914        ,40001
  122. cp,        next,        uy        ,1914        ,40001
  123. cp,        next,        uz        ,1914        ,40001
  124. cp,        next,        roty        ,1914        ,40001
  125. cp,        next,        rotz        ,1914        ,40001
  126. cp,        next,        ux        ,101914        ,140001
  127. cp,        next,        uy        ,101914        ,140001
  128. cp,        next,        uz        ,101914        ,140001
  129. cp,        next,        roty        ,101914        ,140001
  130. cp,        next,        rotz        ,101914        ,140001
  131. cp,        next,        ux        ,101904        ,130001
  132. cp,        next,        uy        ,101904        ,130001
  133. cp,        next,        uz        ,101904        ,130001
  134. cp,        next,        roty        ,101904        ,130001
  135. cp,        next,        rotz        ,101904        ,130001
  136. allsel
  137. nsel,s,,,30000,130000,100000
  138. nsel,a,,,40000,140000,100000
  139. nummrg,node
  140. allsel

  142. *do,i,1,3,1
  143. cp,next,all,1800+i,1900+i
  144. *enddo
  145. allsel
  146. *do,i,5,6,1
  147. cp,next,all,1800+i,1900+i
  148. *enddo
  149. allsel
  150. *do,i,1,2,1              
  151. cp,next,all,1800+i+11,1900+i+11
  152. cp,next,all,i+11,2000+i+11,i+100000+11
  153. *enddo                  
  154. allsel
  155. *do,i,4,5,1                                                
  156. cp,next,all,1800+i+11,1900+i+11         
  157. cp,next,all,i+11,2000+i+11,i+100000+11   
  158. *enddo                                   
  159. allsel   
  160. *do,i,1,3,1                                                              
  161. cp,next,all,i,2000+i,100000+i
  162. *enddo                                   
  163. allsel  
  164. *do,i,5,6,1                                                              
  165. cp,next,all,i,2000+i,100000+i
  166. *enddo                                   
  167. allsel
  168. *do,i,1,3,1                                                              
  169. cp,next,all,101800+i,101900+i
  170. *enddo                                   
  171. allsel   
  172. *do,i,5,6,1                                                              
  173. cp,next,all,101800+i,101900+i
  174. *enddo                                   
  175. allsel
  176. *do,i,1,2,1                                                              
  177. cp,next,all,101800+i+11,101900+i+11
  178. *enddo                                   
  179. allsel
  180. *do,i,4,5,1                                                              
  181. cp,next,all,101800+i+11,101900+i+11
  182. *enddo                                   
  183. allsel

  185. nsel,s,loc,z,-13
  186. nsel,r,loc,y,-3,3
  187. nsel,r,loc,x,8,20
  188. cp,next,all,all
  189. allsel   
  190. nsel,s,loc,z,-13
  191. nsel,r,loc,y,-3,3
  192. nsel,r,loc,x,-8,-20
  193. cp,next,all,all
  194. allsel
  195. nsel,s,loc,z,-13
  196. nsel,r,loc,y,122,128
  197. nsel,r,loc,x,-8,-20
  198. cp,next,all,all
  199. allsel
  200. nsel,s,loc,z,-13
  201. nsel,r,loc,y,122,128
  202. nsel,r,loc,x,8,20
  203. cp,next,all,all
  204. allsel
  205. nsel,s,loc,z,-13
  206. nsel,r,loc,y,-122,-128
  207. nsel,r,loc,x,8,20
  208. cp,next,all,all
  209. allsel
  210. nsel,s,loc,z,-13
  211. nsel,r,loc,y,-122,-128
  212. nsel,r,loc,x,-8,-20
  213. cp,next,all,all
  214. allsel   

  216. *do,i,1,17,1
  217. cp,        next,        ux,        i*100+4                        ,i*100+4+3000
  218. cp,        next,        ux,        i*100+4+6000                ,i*100+4+10000
  219. cp,        next,        uy,        i*100+4                        ,i*100+4+3000
  220. cp,        next,        uy,        i*100+4+6000                ,i*100+4+10000
  221. cp,        next,        uz,        i*100+4                        , i*100+4+3000
  222. cp,        next,        uz,        i*100+4+6000                ,i*100+4+10000
  223. cp,        next,        ux,        i*100+4+10                ,i*100+4+3000+10
  224. cp,        next,        ux,        i*100+4+6000+10                ,i*100+4+10000+10
  225. cp,        next,        uy,        i*100+4+10                ,i*100+4+3000+10
  226. cp,        next,        uy,        i*100+4+6000+10                ,i*100+4+10000+10
  227. cp,        next,        uz,        i*100+4+10                ,i*100+4+3000+10
  228. cp,        next,        uz,        i*100+4+6000+10                ,i*100+4+10000+10
  229. cp,        next,        ux,        i*100+4+100000                ,i*100+4+3000+100000
  230. cp,        next,        ux,        i*100+4+6000+100000        ,i*100+4+10000+100000
  231. cp,        next,        uy,        i*100+4+100000                , i*100+4+3000+100000
  232. cp,        next,        uy,        i*100+4+6000+100000        ,i*100+4+10000+100000
  233. cp,        next,        uz,        i*100+4+100000                ,i*100+4+3000+100000
  234. cp,        next,        uz,        i*100+4+6000+100000        ,i*100+4+10000+100000
  235. cp,        next,        ux,        i*100+4+100000+10        ,i*100+4+3000+100000+10
  236. cp,        next,        ux,        i*100+4+6000+100000        ,i*100+4+10000+100000+10
  237. cp,        next,        uy,        i*100+4+100000+10        ,i*100+4+3000+100000+10
  238. cp,        next,        uy,        i*100+4+6000+100000        ,i*100+4+10000+100000+10
  239. cp,        next,        uz,        i*100+4+100000+10        ,i*100+4+3000+100000+10
  240. cp,        next,        uz,        i*100+4+6000+100000        ,i*100+4+10000+100000+10
  241. *enddo
  242. allsel

  244. allsel
  245. nsel,u,,,erase
  246. nsel,s,loc,x,13.5
  247. nsel,a,loc,x,9
  248. nsel,a,loc,x,18
  249. nsel,a,loc,x,-11.5
  250. nsel,a,loc,x,-7
  251. nsel,a,loc,x,-16
  252. nsel,a,loc,y,3.7
  253. nsel,a,loc,y,-1.7
  254. nsel,a,loc,y,128.7
  255. nsel,a,loc,y,123.3
  256. nsel,a,loc,y,-128.7
  257. nsel,a,loc,y,-123.3
  258. nsel,a,loc,z,-40,-13
  259. d,all,all
  260. allsel
  261. nsel,s,loc,z,-36,5
  262. d,all,all
  263. /color,pbak,off
  264. /replot
  265. allsel
  266. nsel,s,loc,y,0
  267. nsel,r,loc,z,0
  268. nummrg,node,all
  269. esel,none
  270. allsel

  272. /solu
  273. allsel
  274. mpdele,dens,all
  275. mp,dens,3,2600
  276. mp,dens,4,2600
  277. acel,,,9.8
  278. pstres,on
  279. solve
  280. save


  283. /solu
  284. allsel
  285. mpdele,dens,3,11,1
  286. mp,dens,1,7800
  287. mp,dens,2,2600
  288. *do,i,1,17,1
  289. r,1000+i,dgarea,0
  290. *enddo
  291. r,101,xgarea,0
  292. r,102,xgarea,0
  293. acel,,,9.8
  294. solve

  296. /solu
  297. allsel
  298. *do,i,1,17,1
  299. r,1000+1,dgarea,0
  300. *enddo
  301. r,101,xgarea,0
  302. r,102,xgarea,0
  303. r,101,xgarea,0
  304. r,102,xgarea,0
  305. qmb=(250/1.5)*32/6.8*5*1000*9.8
  306. acel,,,9.8
  307. mpdele,dens,3,4,1
  308. mp,dens,1,7800
  309. mp,dens,2,2600
  310. mp,dens,3,4167.7
  311. mp,dens,4,4167.7

  313. mp,dens,5,7800
  314. mp,dens,6,7800
  315. mp,dens,7,2600
  316. mp,dens,8,2600
  317. mp,dens,9,2600
  318. mp,dens,10,7800
  319. mp,dens,11,2600
  320. solve

  322. /solu
  323. allsel
  324. acel,,,9.8
  325. r,101,xgarea,Listrn(1)*1.3
  326. r,102,xgarea,Listrn(2)*1.3
  327. pstres,on
  328. solve


  331. /solu
  332. allsel
  333. r,        1001,        dgarea,1.6787e-03+7.64e-05
  334. r,        1002,        dgarea,1.7538e-03+8.48e-07
  335. r,        1003,        dgarea,1.7547e-03+1.03e-06
  336. r,        1004,        dgarea,1.7400e-03+2.71e-06
  337. r,        1005,        dgarea,1.7292e-03+6.17e-06
  338. r,        1006,        dgarea,1.7303e-03+6.37e-06
  339. r,        1007,        dgarea,1.7267e-03+6.47e-06
  340. r,        1008,        dgarea,1.7279e-03+7.69e-06
  341. r,        1009,        dgarea,1.7365e-03+7.91e-06
  342. r,        1010,        dgarea,1.7277e-03+7.81e-06
  343. r,        1011,        dgarea,1.7263e-03+6.56e-06
  344. r,        1012,        dgarea,1.7299e-03+6.24e-06
  345. r,        1013,        dgarea,1.7290e-03+5.53e-06
  346. r,        1014,        dgarea,1.7412e-03+2.03e-06
  347. r,        1015,        dgarea,1.7504e-03+3.03e-06
  348. r,        1016,        dgarea,1.7587e-03+1.28e-05
  349. r,        1017,        dgarea,1.6143e-03+4.90e-05
  350. pstres,on
  351. solve
  352. finish

  354. /solu
  355. antype,buckle
  356. bucopt,subsp,4
  357. outpr,nsol,all
  358. save
  359. solve
  360. finish
  361. /solu
  362. expass,on
  363. mxpand,4,,,
  364. outpr,all
  365. outres
  366. solve
  367. finish

  369. save
  370. finish
复制代码


代码提供者：maple520
回复 支持 反对

使用道具 举报

AaronSpark
发表于 2006-11-11 07:39 | 显示全部楼层
某铁路钢管混凝土拱桥

  1. /prep7  
  2. /title,analysis of a cftb
  3.   
  4. !主梁节点
  5.   
  6. /prep7
  7. z=-3.0124   !(=-6.0248/2)
  8. y=-22.9207  ! (=-45.3653/2-0.2380)
  9. b=0
  10.   
  11. n,b+1,-84,y,z  $  n,b+2,-69,y,z
  12. n,b+3,-57,y,z   $  n,b+4,-45,y,z
  13.   
  14. *do,i,b+5,b+14,1
  15. n,i,-45+9*(i-4),y,z
  16. *enddo
  17.   
  18. n,b+15,56,y,z
  19.   
  20. ngen,2,15,all,,,,,5.72  
  21.   
  22. !主梁左侧节点由对称产生，这里不写
  23. !吊杆底部节点放在后面
  24.   
  25. ！拱肋关键点
  26.   
  27. c=30  !c的值可变
  28. n,c+1,0,1.3454,-0.2130   $  n,c+26,0,-1.3454,0.2130
  29. n,c+2,-3,1.2736,-0.2016   $  n,c+27,-3,-1.4182,0.2245
  30. n,c+3,-6,1.0591,-0.1677   $  n,c+28,-6,-1.6374,0.2592
  31. n,c+4,-9,0.7014,-0.1110   $  n,c+29,-9,-2.0027,0.3170
  32. n,c+5,-12,0.2006,-0.0316  $  n,c+30,-12,-2.5146,0.3981
  33. n,c+6,-15,-0.4438,0.0702 $  n,c+31,-15,-3.1728,0.5023
  34. n,c+7,-18,-1.2311,0.1949  $  n,c+32,-18,-3.9785,0.6298
  35. n,c+8,-21,-2.1621,0.3423  $  n,c+33,-21,-4.9315,0.7807
  36. n,c+9,-24,-3.2372,0.5125  $  n,c+34,-24,-6.0321,0.9549
  37. n,c+10,-27,-4.4564,0.7055  $  n,c+35,-27,-7.2811,1.1526
  38. n,c+11,-30,-5.8205,0.9214  $  n,c+36,-30,-8.6785,1.3738
  39. n,c+12,-33,-7.3287,1.1602  $  n,c+37,-33,-10.2262,1.6188
  40. n,c+13,-36,-8.9825,1.4219  $  n,c+38,-36,-11.9238,1.8876
  41. n,c+14,-39,-10.7820,1.7068  $  n,c+39,-39,-13.7722,2.1802
  42. n,c+15,-42,-12.7303,2.0152  $  n,c+40,-42,-15.7699,2.4964
  43. n,c+16,-45,-14.8188,2.3459  $  n,c+41,-45,-17.9269,2.8379
  44. n,c+17,-48,-17.0570,2.7002  $  n,c+42,-48,-20.2351,3.2033
  45. n,c+18,-51,-19.4426,3.0778  $  n,c+43,-51,-22.6984,3.5932
  46. n,c+19,-54,-21.9766,3.4790  $  n,c+44,-54,-25.3178,4.0079
  47. n,c+20,-57,-24.6590,3.9036  $  n,c+45,-57,-28.0949,4.4475
  48. n,c+21,-60,-27.4903,4.3518  $  n,c+46,-60,-31.0313,4.9123
  49. n,c+22,-63,-30.4717,4.8238  $  n,c+47,-63,-34.1279,5.4025
  50. n,c+23,-66,-33.6029,5.3194  $  n,c+48,-66,-37.3871,5.9185
  51. n,c+24,-69,-36.8825,5.8386  $  n,c+49,-69,-40.8127,6.4608
  52. n,c+25,-71,-39.1646,6.1998  
  53.   
  54. nsel,s,,,31,79,1
  55. ngen,2,49,all,,,,0.2345/2,1.4815/2   
  56. !（49＝79+1-31）
  57.   
  58. nsel,s,,,80,128,1
  59. ngen,2,49,all,,,,0.2345/2,1.4815/2
  60.   
  61. nsel,s,,,80,128,1
  62. ngen,2,49,all,,,,0.2345/2,1.4815/2
  63.   
  64. nsel,all
  65. nsel,s,,,32,49,1
  66. nsym,x,146 ,all,,,
  67. nsel,s,,,57,74,1
  68. nsym,x,139 ,all,,,
  69.   
  70. nsel,s,,,81,98,1
  71. nsym,x,133 ,all,,,
  72.   
  73. nsel,s,,,106,123,1
  74. nsym,x,126 ,all,,,
  75.   
  76. nsel,s,,,130,147,1
  77. nsym,x,120 ,all,,,
  78.   
  79. nsel,s,,,155,172,1
  80. nsym,x,113,all,,,
  81.   
  82. nsel,all
  83. nplot
  84. !补右拱脚节点
  85. n,286,55,-26.221,4.1509  
  86. n,287,55,-26.1038,4.8917
  87. n,288,55,-25.9865,5.6324
  88. n,289,57,-24.6590,3.9036   
  89. n,290,57,-24.4245,4.64435
  90. n,291,57,-24.4245,5.3851
  91.   
  92. ！横撑中心节点
  93. nsel,s,,,5
  94. ngen,2,287,all,,,0,10.9969,0
  95. nsel,s,,,5
  96. ngen,2,288,all,,,0,13.9382,0
  97.   
  98. nsel,s,,,7
  99. ngen,2,287,all,,,0,18.9422,0
  100. nsel,s,,,7
  101. ngen,2,288,all,,,0,21.6896,0
  102.   
  103. nsel,s,,,9
  104. ngen,2,287,all,,,0,21.5753,0
  105. nsel,s,,,9
  106. ngen,2,288,all,,,0,24.2661,0
  107.   
  108. nsel,s,,,11
  109. ngen,2,287,all,,,0,18.9422,0
  110. nsel,s,,,11
  111. ngen,2,288,all,,,0,21.6896,0
  112.   
  113. nsel,s,,,13
  114. ngen,2,287,all,,,0,10.9969,0
  115. nsel,s,,,13
  116. ngen,2,288,all,,,0,13.9382,0
  117.   
  118. nsel,all
  119. nplot
  120.   
  121. *do,i,1,11,1
  122. n,301+i,-45+9*(i-1),-22.6827,2.6876
  123. *enddo
  124.   
  125. *do,i,1,6,1
  126. nsel,s,,,105+3*(i-1)   
  127. ngen,2,208-2*(i-1),all,,,0,-0.4981,0.01358
  128. *enddo
  129.   
  130. *do,i,1,5,1
  131. nsel,s,,,234+3*(i-1)   
  132. ngen,2,85-2*(i-1),all,,,0,-0.4981,0.01358
  133. *enddo
  134.   
  135. nsel,all
  136. nplot
  137.   
  138. !立柱
  139. !324，2要耦合
  140. !325，3要耦合
  141. !326，17要耦合
  142. !327，18要耦合
  143. nsel,s,,,2,3,1
  144. ngen,2,322,all,,,0,-(1.6-0.238),0
  145. nsel,s,,,17,18,1
  146. ngen,2,309,all,,,0,-(1.6-0.238),0
  147.   
  148. nsel,s,,,326
  149. ngen,2,2,all,,,0,0,3.131
  150. nsel,s,,,327
  151. ngen,2,2,all,,,0,0,1.196
  152.   
  153. nsel,all
  154. nplot
  155.   
  156. !定义单元类型
  157. et,1,beam4
  158. et,2,link10
  159. keyopt,2,2,2   
  160. keyopt,2,3,0   
  161. et,3,shell63
  162.   
  163. !定义材料
  164. !拱肋弦杆
  165. MP,EX,1,0.527e11
  166. MP,nuxy,1,0.25
  167. MP,dens,1,3188
  168. !拱肋腹杆，横撑强、弱构件
  169. MP,EX,2,2.06e11
  170. MP,nuxy,2,0.3
  171. MP,dens,2,7850
  172. !拱肋上下平联
  173. MP,EX,3,0.355e11
  174. MP,nuxy,3,0.2
  175. MP,dens,3,2600
  176. !吊杆
  177. MP,EX,4,1.95e11
  178. MP,nuxy,4,0.3
  179. MP,dens,4,7850
  180. !立柱、立柱横梁
  181. MP,EX,5,0.3e11
  182. MP,nuxy,5,0.2
  183. MP,dens,5,2500
  184. !鱼骨梁主梁
  185. MP,EX,6,0.355e11
  186. MP,nuxy,6,0.2
  187. MP,dens,6,4119
  188. !鱼骨梁刺梁
  189. MP,EX,7,3.55e11     !弹性模量输入主梁的十倍
  190. MP,nuxy,7,0.2
  191. MP,dens,7,0.1
  192.   
  193. !定义实常数（面积除了拱肋弦杆外，其他输入实际面积以和密度对应）
  194. !拱肋弦杆
  195. r,1,0.37,0.03956,0.03956,0.7,0.7
  196. rmore,,0.02357
  197. !拱肋腹杆、横撑强构件
  198. r,2,0.01225,0.000233,0.000233,0.4,0.4
  199. !拱肋横撑弱构件
  200. r,3,0.004825,0.0000233,0.0000233,0.2,0.2
  201. !立柱
  202. r,4,1,0.083333,0.083333,1,1
  203. rmore,,0.1426
  204. !立柱横梁
  205. r,5,1.44,0.1728,0.1728,1.2,1.2
  206. rmore,,0.2957
  207. !鱼骨梁主梁
  208.   
  209. r,6,6.718,2.061,103.4,2.5,1.6  !注意输入
  210. rmore,0.02357
  211. !鱼骨梁刺梁
  212. r,7,0.01,1030.4,20.016,0.1,0.1 !抗弯惯性矩为主梁的十倍
  213.   
  214. !吊杆
  215. r,8,0.00535,1e-8
  216. !上下平联板
  217. r,9,0.5，0.5，0.5，0.5，
  218.   
  219. !主梁
  220.   
  221. ！刺梁
  222. type,1  
  223. mat,7
  224. real,7
  225.   
  226. *do,i,1,15,1
  227. e,i,i+15
  228. *enddo
  229.   
  230. ！主梁
  231. type,1  
  232. mat,6  
  233. real,6
  234.   
  235. *do,i,1,14,1
  236. e,i,i+1
  237. !e,e+15,i+16  !三主梁的边梁
  238. *enddo
  239.   
  240. !拱肋
  241.   
  242. ！弦杆
  243. type,1  
  244. mat,1  
  245. real,1
  246.   
  247. *do,i,1,24,1
  248. e,i+30,i+31
  249. e,i+128,i+129
  250. *enddo
  251.   
  252. *do,i,1,23,1
  253. e,i+55,i+56
  254. e,i+153,i+154
  255. *enddo
  256.   
  257. *do,i,1,17,1
  258. e,i+177,i+178
  259. e,i+195,i+196
  260. e,i+249,i+250
  261. e,i+267,i+268
  262. *enddo
  263.   
  264. e,31,178  $ e,129,250
  265. e,56,196  $ e,154,268
  266. e,213,286  $ e,285,288
  267. e,195,289  $ e,267,291
  268.   
  269. ！竖腹杆
  270. type,1  
  271. mat,2  
  272. real,2
  273.   
  274. *do,i,1,24,1
  275. e,i+30,i+55
  276. e,i+128,i+153
  277. *enddo
  278.   
  279. *do,i,1,18,1
  280. e,i+177,i+195
  281. e,i+249,i+267
  282. *enddo
  283.   
  284. ！斜腹杆
  285. type,1  
  286. mat,2  
  287. real,2
  288.   
  289. *do,i,1,23,1
  290. e,i+31,i+55
  291. e,i+129,i+153
  292. *enddo
  293.   
  294. *do,i,1,17,1
  295. e,i+178,i+195
  296. e,i+250,i+267
  297. *enddo
  298.   
  299. e,178,56  $  e,250,154
  300. e,213,289  $ e,285,291
  301.   
  302. !横撑
  303. ！横撑强构件
  304. type,1  
  305. mat,2  
  306. real,2
  307.   
  308. e,292,293  $   e,292,43   
  309. e,294,295  $    e,294,37
  310. e,296,297  $  e,296,31
  311. e,298,299  $  e,298,183
  312. e,300,301  $ e,300,189
  313.   
  314. ！横撑弱构件
  315. type,1  
  316. mat,2  
  317. real,3
  318.   
  319. e,293,43   $ e,292,68
  320. e,295,37   $  e,294,62
  321. e,297,31  $  e,296,56
  322. e,299,183  $ e,298,201
  323. e,301,189  $ e,300,207
  324.   
  325. !吊杆
  326. type,2
  327. mat,4  
  328. real,8
  329.   
  330. e,302,318  $ e,303,317  $ e,304,316
  331. e,305,315  $ e,306,314  $ e,307,313
  332. e,308,319  $ e,309,320  $ e,310,321
  333. e,311,322  $ e,312,323
  334.   
  335. e,95,120  $ e,92,117  $ e,89,114
  336. e,86,111  $ e,83,108  $ e,80,105
  337. e,216,234  $ e,219,237  $ e,222,240
  338. e,225,243  $ e,228,246
  339.   
  340. !立柱
  341. ！立柱横梁
  342. type,1  
  343. mat,5  
  344. real,5
  345. e,324,326  $ e,326,328   
  346. e,325,327  $ e,327,329   
  347.   
  348. ！竖立柱
  349. type,1  
  350. mat,5  
  351. real,4
  352.   
  353. e,328,54   $ e,329,50
  354.   
  355. !上下平联板
  356. type,3  
  357. mat,3  
  358. real,9
  359.   
  360. ！左平联
  361. *do,i,1,24,1
  362. e,i+30,i+31,i+80,i+79
  363. e,i+79,i+80,i+129,i+128
  364. *enddo
  365.   
  366. *do,i,1,23,1
  367. e,i+55,i+56,i+105,i+104
  368. e,i+104,i+105,i+154,i+153
  369. *enddo
  370.   
  371. e,31,178,214,80  $ e,80,214,250,129
  372. e,56,196,232,105  $ e,105,232,268,154
  373.   
  374. ！右平联
  375. *do,i,1,17,1
  376. e,i+177,i+178,i+214,i+213
  377. e,i+213,i+214,i+250,i+249
  378. e,i+195,i+196,i+232,i+231
  379. e,i+231,i+232,i+268,i+267
  380. *enddo
  381.   
  382. e,195,289,290,231  $ e,231,290,291,267
  383. e,213,286,287,249  $ e,249,287,288,285
  384.   
  385. /eshape,1
  386. eplot
  387.   
  388. cs,11,0,297,296,298,1,1
  389. csys,11
  390. dsys,11
  391.   
  392. nsym,z,329,all  !前面一共three two nine 个节点   
  393. ensym,585,,329,all !前面一共five eight five个单元   
  394. eplot
  395.   
  396. csys,0   
  397. dsys,0  
  398. /replo
  399.   
  400. nummrg,all, , , ,low
  401. numcmp,all
  402.   
  403. *do,i,1,6,1
  404. cp,next,ux,312+i,102+3*i  
  405. cp,next,uy,312+i,102+3*i  
  406. cp,next,uz,312+i,102+3*i  
  407. *enddo
  408.   
  409. *do,i,1,5,1
  410. cp,next,ux,318+i,231+3*i  
  411. cp,next,uy,318+i,231+3*i  
  412. cp,next,uz,318+i,231+3*i  
  413. *enddo
  414.   
  415. *do,i,1,11,1
  416. cp,next,ux,301+i,18+i  
  417. cp,next,uy,301+i,18+i  
  418. cp,next,uz,301+i,18+i  
  419. *enddo
  420.   
  421. *do,i,1,6,1
  422. cp,next,ux,616+i,416+3*i  
  423. cp,next,uy,616+i,416+3*i   
  424. cp,next,uz,616+i,416+3*i  
  425. *enddo
  426.   
  427. *do,i,1,5,1
  428. cp,next,ux,622+i,545+3*i  
  429. cp,next,uy,622+i,545+3*i   
  430. cp,next,uz,622+i,545+3*i  
  431. *enddo
  432.   
  433. *do,i,1,11,1
  434. cp,next,ux,605+i,332+i  
  435. cp,next,uy,605+i,332+i  
  436. cp,next,uz,605+i,332+i  
  437. *enddo
  438.   
  439. *do,i,1,2,1
  440. cp,next,ux,16+i,325+i  
  441. cp,next,uy,16+i,325+i  
  442. cp,next,uz,16+i,325+i
  443. cp,next,rotx,16+i,325+i  
  444. cp,next,roty,16+i,325+i  
  445. cp,next,rotz,16+i,325+i
  446. *enddo
  447.   
  448. *do,i,1,2,1
  449. cp,next,ux,1+i,323+i  
  450. cp,next,uy,1+i,323+i  
  451. cp,next,uz,1+i,323+i
  452. cp,next,rotx,1+i,323+i  
  453. cp,next,roty,1+i,323+i  
  454. cp,next,rotz,1+i,323+i
  455. *enddo
  456.   
  457. *do,i,1,2,1
  458. cp,next,ux,330+i,627+i  
  459. cp,next,uy,330+i,627+i  
  460. cp,next,uz,330+i,627+i
  461. cp,next,rotx,330+i,627+i  
  462. cp,next,roty,330+i,627+i  
  463. cp,next,rotz,330+i,627+i
  464. *enddo
  465.   
  466. /solu  
  467.    
  468. antype,2  
  469. modopt,lanb,10   
  470. mxpand,10, , ,0  
  471. lumpm,0  
  472. pstres,0
  473. modopt,lanb,10,0,100, ,off
  474.   
  475. d,286, , , , 288,2 ,all
  476. d,289, , , , 291,2 ,all
  477. d,600, , , , 602,2 ,all
  478. d,603, , , , 605,2 ,all
  479. d,467, , , , , ,all
  480. d,369, , , , , ,all
  481. d,491, , , , , ,all
  482. d,393, , , , , ,all
  483. d,55, , , , , ,all
  484. d,153, , , , , ,all
  485. d,79, , , , , ,all
  486. d,177, , , , , ,all
  487. d,344,,,,,,ux,uy,uz
  488. d,30,,,,,,ux,uy,uz
  489. d,15,,,,,,ux,uy,uz
  490. d,330,,,,,,uy,uz
  491. d,1,,,,,,uy,uz
  492. d,16,,,,,,uy,uz
  493.   
  494. solve
  495.   
  496. finish
复制代码


代码提供者：西王赏功
回复 支持 反对

使用道具 举报

valen_zt
发表于 2012-2-19 11:47 | 显示全部楼层
回复 5 # AaronSpark 的帖子

这两个都是钢管混凝土拱桥的例子吧。。如果想要分别用壳、实体单元去模拟钢管混凝土柱。。。应该怎么模拟?特别是他们之间的接触应该怎么去设置？望高手不吝解答。。。耐心等待中。。感谢！
回复 支持 反对

使用道具 举报

ymcheng123
发表于 2012-4-17 16:23 | 显示全部楼层
看不懂，好复杂啊
回复 支持 反对

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-12-30 01:28 , Processed in 0.057695 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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