声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3693|回复: 2

[结构分析] [分享]斜拉桥模型命令流

[复制链接]
发表于 2005-8-27 16:49 | 显示全部楼层 |阅读模式

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

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

x
  1. /com,new model of linjiang cable_stayed bridge,2004.2.09
  2. /prep7
  3. /title, cable_stayed bridge,author is Sunhang
  4. /com,define the keypoints
  5. *set,alfa1,10 !angle of tower upside
  6. *set,alfa2,65 !angle of tower downside
  7. *set,alfa3,79.04594 !angle of tower with bridge surface
  8. *set,y1,55.5 !桥塔顶面到原点的距离
  9. *set,y2,33.5 !桥塔中部的Y轴向长度
  10. *set,pi,3.1415926
  11. *set,x3,y2/tan(alfa2*pi/180) !桥塔中部的X轴向长度
  12. *set,x2,(y1-y2)*tan(alfa1*pi/180) !桥塔上部的X轴向长度
  13. *set,x1,x2+x3 !桥塔的X轴向长度

  14. *set,kp_yy1,0 !定义桥塔上部的索锚固点竖向距离(从塔顶算起)
  15. *set,kp_yy2,2.5185
  16. *set,kp_yy3,3.5788

  17. *set,kp_yy4,4.6469
  18. *set,kp_yy5,5.7248
  19. *set,kp_yy6,6.8151
  20. *set,kp_yy7,7.9211
  21. *set,kp_yy8,9.0479
  22. *set,kp_yy9,10.2027

  23. *set,kp_yy10,11.3965
  24. *set,kp_yy11,12.6470
  25. *set,kp_yy12,13.9848
  26. *set,kp_yy13,15.7143
  27. *set,kp_yy14,17.7041
  28. *set,kp_yy15,22.0000

  29. k,1,-x1,y1,
  30. k,6,-x1+kp_yy2*tan(alfa1*pi/180),y1-kp_yy2
  31. k,8,-x1+kp_yy3*tan(alfa1*pi/180),y1-kp_yy3
  32. k,10,-x1+kp_yy4*tan(alfa1*pi/180),y1-kp_yy4
  33. k,12,-x1+kp_yy5*tan(alfa1*pi/180),y1-kp_yy5
  34. k,14,-x1+kp_yy6*tan(alfa1*pi/180),y1-kp_yy6
  35. k,16,-x1+kp_yy7*tan(alfa1*pi/180),y1-kp_yy7
  36. k,18,-x1+kp_yy8*tan(alfa1*pi/180),y1-kp_yy8
  37. k,20,-x1+kp_yy9*tan(alfa1*pi/180),y1-kp_yy9
  38. k,22,-x1+kp_yy10*tan(alfa1*pi/180),y1-kp_yy10
  39. k,24,-x1+kp_yy11*tan(alfa1*pi/180),y1-kp_yy11
  40. k,26,-x1+kp_yy12*tan(alfa1*pi/180),y1-kp_yy12
  41. k,28,-x1+kp_yy13*tan(alfa1*pi/180),y1-kp_yy13
  42. k,30,-x1+kp_yy14*tan(alfa1*pi/180),y1-kp_yy14
  43. k,38,-x1+kp_yy15*tan(alfa1*pi/180),y1-kp_yy15

  44. kfill,1,6 !在已建关键点内内插关键点
  45. kfill,6,8
  46. kfill,8,10
  47. kfill,10,12
  48. kfill,12,14
  49. kfill,14,16
  50. kfill,16,18
  51. kfill,18,20
  52. kfill,20,22
  53. kfill,22,24
  54. kfill,24,26
  55. kfill,26,28
  56. kfill,28,30
  57. kfill,30,38
  58. *set,kp_numone,38!定义桥塔上部的最后一个关键点号

  59. /com,define and mesh the above part of tower
  60. *dim,BBS,,40
  61. *dim,HHS,,40
  62. *dim,SSS,,40
  63. *dim,IIYYS,,40
  64. *dim,IIZZS,,40
  65. *set,length1,2.9546108
  66. *set,length2,4.9251168
  67. *set,width1,3
  68. *set,width2,5
  69. *set,diff1_length,length2-length1 !桥塔上部两个截面的长度之差(纵桥向)
  70. *set,diff1_width,width2-width1 !桥塔上部两个截面的宽度之差(横桥向)


  71. *dim,yy,,40 !定义桥塔上部的关键点竖向距离数组(从塔顶计算)
  72. *do,i,1,38
  73. *set,yy(i),y1-ky(i)
  74. *set,hhs(i),diff1_length*yy(i)/(y1-y2)+length1
  75. *set,bbs(i),diff1_width*yy(i)/(y1-y2)+width1
  76. SSS(i)=BBS(i)*HHS(i) !按照实心截面
  77. IIZZS(i)=BBS(i)*HHS(i)**3/12
  78. IIyyS(i)=HHS(i)*BBS(i)**3/12
  79. *enddo


  80. *do,i,1,kp_numone-1
  81. et,i,beam44
  82. keyopt,i,9,0
  83. mp,ex,i,3.5e10
  84. mp,prxy,i,0.167
  85. mp,dens,i,2.6e3
  86. mp,alpx,i,1e-05 !定义混凝土的线膨胀系数
  87. r,i,SSs(i),IIZZS(i),IIyyS(i),BBS(i)/2,HHS(i)/2,IIZZS(i)+IIyyS(i)
  88. rmore,SSs(i+1),IIZZS(i+1),IIyyS(i+1),BBS(i+1)/2,HHS(i+1)/2,IIZZS(i+1)+IIyyS(i+1)
  89. rmore,0,0,0,0,0,0
  90. rmore,0,0,BBS(i)/2,HHS(i)/2,BBS(i+1)/2,HHS(i+1)/2
  91. lsel,u,real,,1,kp_numone
  92. lstr,i,i+1
  93. latt,i,i,i
  94. lesize,all,,,1,,1
  95. lmesh,all,all
  96. allsel,all
  97. *enddo

  98. cm,uptower,elem

  99. *get,emin_ts,elem,0,num,min
  100. *get,emax_ts,elem,0,num,max
  101. allsel
  102. kp_ts=emax_ts-emin_ts+1 !提取出来桥塔上部的关键点数目

  103. /com,createt the kps of tower down
  104. *set,zfirst,1.25 !三个控制点的横桥向坐标
  105. *set,zsecond,7.4839
  106. *set,kp_numtwo,kp_ts+2 !定义桥塔下部的第一个关键点
  107. *set,diff1_elem,67 !桥塔下部的单元数目
  108. *set,kp_numthree,kp_numtwo+diff1_elem !桥塔下部的最后一个关键点
  109. k,kp_numtwo,-y2/tan(alfa2*pi/180),y2 ,-zfirst
  110. k,kp_numthree,,,-zsecond
  111. kfill,kp_numtwo,kp_numthree

  112. *dim,BBX,,200
  113. *dim,HHX,,200
  114. *dim,SSX,,200
  115. *dim,IIYYX,,200
  116. *dim,IIZZX,,200
  117. *set,HHXQ,0.8*sin(alfa2*pi/180) !定义砍掉部分的长度
  118. *set,BBXQ,0.4*sin(alfa3*pi/180)
  119. *set,SSXQ,HHXQ*BBXQ
  120. *dim,IIYYXQ,,200
  121. *dim,IIZZXQ,,200
  122. !*set,w4,2.6875*2 !在桥塔折角处单元的宽度
  123. !*set,w5,4.3682*2 !在坐标原点处单元的宽度
  124. *set,w4,4.5361558
  125. *set,w5,8.1609208

  126. *do,i,1,diff1_elem+1
  127. BBX(i+kp_numone)=2.5*sin(alfa3*3.1415926/180)
  128. HHX(i+kp_numone)=w4+(i-1)*(w5-w4)/67
  129. SSX(i+kp_numone)=BBX(i+kp_numone)*HHX(i+kp_numone)-SSXQ*2
  130. IIYYXQ(i+kp_numone)=HHXQ*BBXQ**3/12+SSXQ*BBX(i+kp_numone)**2/4
  131. IIZZXQ(i+kp_numone)=BBXQ*HHXQ**3/12+SSXQ*HHX(i+kp_numone)**2/4
  132. IIYYX(i+kp_numone)=HHX(i+kp_numone)*BBX(i+kp_numone)**3/12-IIYYXQ(i+kp_numone)*2
  133. IIZZX(i+kp_numone)=BBX(i+kp_numone)*HHX(i+kp_numone)**3/12-IIZZXQ(i+kp_numone)*2
  134. *enddo

  135. *do,i,KP_numone,kp_numone+diff1_elem-1
  136. et,i,beam44
  137. keyopt,i,9,0
  138. mp,ex,i,3.5e10
  139. mp,prxy,i,0.167
  140. mp,dens,i,2600
  141. mp,alpx,i,1e-05 !定义混凝土的线膨胀系数
  142. r,i,SSX(i+1),IIZZX(i+1),IIYYX(i+1),BBX(i+1)/2,HHX(i+1)/2,IIZZX(i+1)+IIYYX(i+1)
  143. rmore,SSX(i+2),IIZZX(i+2),IIYYX(i+2),BBX(i+2)/2,HHX(i+2)/2,IIZZX(i+2)+IIYYX(i+2)
  144. rmore,0,0,0,0,0,0
  145. rmore,0,0,BBX(i+1)/2,HHX(i+1)/2,BBX(i+2)/2,HHX(i+2)/2
  146. *enddo
  147. *do,i,kp_numtwo,kp_numthree-1
  148. lsel,u,real,,1,150
  149. lstr,i,i+1
  150. latt,i-1,i-1,i-1
  151. lesize,all,,,1,,1
  152. lmesh,all
  153. allsel,all
  154. *enddo

  155. esel,u,real,,1,kp_numtwo-2
  156. cm,downtower,elem
  157. allsel,all
  158. cmsel,s,downtower
  159. *get,emin_tx,elem,0,num,min
  160. *get,emax_tx,elem,0,num,max
  161. allsel
  162. kp_tx=emax_tx-emin_tx+1

  163. esel,all,all
  164. cm,tower,elem

  165. lsel,s,real,,kp_numone,kp_numthree-2
  166. lsymm,z,all,,,100
  167. allsel,all
  168. *set,kp_numfour,kp_numone+2*kp_tx+100+1 !主梁的第一个关键点号
  169. !esel,s,ename,,beam44
  170. !tunif,0
  171. !tref,-30
  172. !allsel,all


  173. /com,couple the tower upside and tower down
  174. cerig,node(kx(38),ky(38),kz(38)),node(kx(39),ky(39),kz(39)),all,
  175. cerig,node(kx(38),ky(38),kz(38)),node(kx(139),ky(139),kz(139)),all,

  176. *set,beam_height,1.2725-0.3  !主梁节点即锚固点到原点的距离
  177. !*set,kp_numfour,kp_numone+2*(kp_tx+1)+100+1 !主梁的第一个关键点号
  178. *set,kp_numfour_inc,210 !主梁的关键点数目-1
  179. *set,kp_numfive,kp_numfour+kp_numfour_inc
  180. k,kp_numfour,,beam_height
  181. k,kp_numfive,105,beam_height
  182. kfill,kp_numfour,kp_numfive
  183. *do,i,1,kp_numfour_inc
  184. lstr,kp_numfour+i-1,kp_numfour+i
  185. *enddo
  186. *set,enum_beam,emax_tx+1
  187. et,enum_beam,beam188
  188. mp,ex,enum_beam,3.5e10
  189. mp,prxy,enum_beam,0.0.167
  190. mp,dens,enum_beam,3038.8 !考虑到二期恒载后的换算密度
  191. mp,alpx,enum_beam,1e-05 !定义混凝土的线膨胀系数
  192. KEYOPT,enum_beam,7,2
  193. keyopt,enum_beam,8,3
  194. keyopt,enum_beam,9,3
  195. SECTYPE,1,BEAM,MESH,sec1
  196. SECOFFSET,user,,-1.40 !截面读入时主梁的平移
  197. SECREAD,'main_beam','SECT',' ',MESH
  198. k,5000,,1000000
  199. lsel,u,real,,1,200
  200. latt,enum_beam,enum_beam,enum_beam,,5000,,1
  201. lesize,all,,,1,,1
  202. lmesh,all
  203. allsel,all
  204. esel,s,ename,,beam188
  205. cm,main_beam,elem
  206. allsel,all

  207. /com,create the cable element
  208. *set,enum_link,enum_beam+1 !拉索的开始单元号
  209. *dim,cable_area,,13 !定义拉索单元的面积数组
  210. *dim,cable_istrain,,13 !定义拉索单元的初始应变数组
  211. *dim,cable_dens,,13

  212. /com,define the angle of all cable
  213. *set,cable_area1,1.668E-03
  214. *set,cable_area2,1.668E-03
  215. *set,cable_area3,2.6410E-03
  216. *set,cable_area4,2.6410E-03
  217. *set,cable_area5,2.6410E-03
  218. *set,cable_area6,2.6410E-03
  219. *set,cable_area7,3.0580E-03
  220. *set,cable_area8,3.0580E-03
  221. *set,cable_area9,3.0580E-03
  222. *set,cable_area10,3.7530E-03
  223. *set,cable_area11,3.7530E-03
  224. *set,cable_area12,3.7530E-03
  225. *set,cable_area13,3.7530E-03
  226. *set,cable_area_back,2.0155E-02


  227. *set,cable_area(1),cable_area1
  228. *set,cable_area(2),cable_area2
  229. *set,cable_area(3),cable_area3
  230. *set,cable_area(4),cable_area4
  231. *set,cable_area(5),cable_area5
  232. *set,cable_area(6),cable_area6
  233. *set,cable_area(7),cable_area7
  234. *set,cable_area(8),cable_area8
  235. *set,cable_area(9),cable_area9
  236. *set,cable_area(10),cable_area10
  237. *set,cable_area(11),cable_area11
  238. *set,cable_area(12),cable_area12
  239. *set,cable_area(13),cable_area13

  240. *set,cable_dens(1),13.2/cable_area1
  241. *set,cable_dens(2),13.2/cable_area2
  242. *set,cable_dens(3),20.9/cable_area3
  243. *set,cable_dens(4),20.9/cable_area4
  244. *set,cable_dens(5),20.9/cable_area5
  245. *set,cable_dens(6),20.9/cable_area6
  246. *set,cable_dens(7),24.2/cable_area7
  247. *set,cable_dens(8),24.2/cable_area8
  248. *set,cable_dens(9),24.2/cable_area9
  249. *set,cable_dens(10),29.7/cable_area10
  250. *set,cable_dens(11),29.7/cable_area11
  251. *set,cable_dens(12),29.7/cable_area12
  252. *set,cable_dens(13),29.7/cable_area13
  253. *set,cable_dens_back,159.5/cable_area_back

  254. *set,cable_istrain1,0.26032E-02
  255. *set,cable_istrain2,0.25568E-02

  256. *set,cable_istrain3,0.23210E-02

  257. *set,cable_istrain4,0.23456E-02
  258. *set,cable_istrain5,0.23892E-02
  259. *set,cable_istrain6,0.24412E-02

  260. *set,cable_istrain7,0.28199E-02
  261. *set,cable_istrain8,0.28719E-02
  262. *set,cable_istrain9,0.29143E-02

  263. *set,cable_istrain10,0.28321E-02
  264. *set,cable_istrain11,0.28559E-02
  265. *set,cable_istrain12,0.28743E-02
  266. *set,cable_istrain13,0.28926E-02

  267. cable_back_istrain1=0.32891E-02
  268. cable_back_istrain2=0.33661E-02

  269. *set,cable_istrain(1),cable_istrain1
  270. *set,cable_istrain(2),cable_istrain2
  271. *set,cable_istrain(3),cable_istrain3
  272. *set,cable_istrain(4),cable_istrain4
  273. *set,cable_istrain(5),cable_istrain5
  274. *set,cable_istrain(6),cable_istrain6
  275. *set,cable_istrain(7),cable_istrain7
  276. *set,cable_istrain(8),cable_istrain8
  277. *set,cable_istrain(9),cable_istrain9
  278. *set,cable_istrain(10),cable_istrain10
  279. *set,cable_istrain(11),cable_istrain11
  280. *set,cable_istrain(12),cable_istrain12
  281. *set,cable_istrain(13),cable_istrain13
复制代码
回复
分享到:

使用道具 举报

 楼主| 发表于 2005-8-27 16:50 | 显示全部楼层

回复:(gaby2033)[分享]斜拉桥模型命令流

<P>/com,create the kps of croSSbeams<BR>*set,w5,5.35  !横梁到主梁中心的距离<BR>*do,i,1,13<BR>k,kp_numfour+i*14+400+12,i*7+6,beam_height,-w5<BR>k,kp_numfour+i*14+800+12,i*7+6,beam_height,w5<BR>*enddo<BR>!定义拉索单元</P>
<P>*do,i,1,13<BR>et,enum_link+i-1,link10 <BR>mp,ex,enum_link+i-1,2.0e11<BR>mp,prxy,enum_link+i-1,0.3<BR>mp,dens,enum_link+i-1,cable_dens(i)<BR>!mp,alpx,enum_link+i-1,1.5e-05  !定义索(钢材)的线膨胀系数<BR>lsel,u,real,,1,1000<BR>r,i+enum_link-1,cable_area(i),cable_istrain(i)<BR>lstr,30-2*(i-1),kp_numfour+i*14+400+12<BR>lstr,30-2*(i-1),kp_numfour+i*14+800+12<BR>latt,enum_link+i-1,enum_link+i-1,enum_link+i-1<BR>lesize,all,,,1,,1<BR>lmesh,all<BR>lsel,all<BR>*enddo</P>
<P>/com,create the back cables<BR>*set,enum_back_cable,enum_link+13<BR>et,enum_back_cable,link10<BR>mp,ex,enum_back_cable,2.0e11<BR>mp,prxy,enum_back_cable,0.3<BR>mp,dens,enum_back_cable,cable_dens_back<BR>r,enum_back_cable,cable_area_back,cable_back_istrain1<BR>lsel,u,real,,1,enum_back_cable-1<BR>k,3001,-45,<BR>k,3002,-45+2.956<BR>lstr,3001,10<BR>latt,enum_back_cable,enum_back_cable,enum_back_cable<BR>lesize,all,,,1,,1<BR>lmesh,all<BR>allsel,all</P>
<P>et,enum_back_cable+1,link10<BR>mp,ex,enum_back_cable+1,2.0e11<BR>mp,prxy,enum_back_cable+1,0.3<BR>mp,dens,enum_back_cable+1,cable_dens_back<BR>r,enum_back_cable+1,cable_area_back,cable_back_istrain2<BR>lsel,u,real,,1,200<BR>lstr,3002,18<BR>latt,enum_back_cable+1,enum_back_cable+1,enum_back_cable+1<BR>lesize,all,,,1,,1<BR>lmesh,all<BR>allsel,all</P>
<P>esel,s,ename,,link10<BR>cm,cable,elem<BR>alls</P>
<P>/com,create the croSSbeams<BR>*set,enum_croSSbeam,enum_back_cable+2 !横梁的单元号<BR>et,enum_croSSbeam,beam4<BR>mp,ex,enum_croSSbeam,3.5e20<BR>mp,prxy,enum_croSSbeam,0.167<BR>mp,dens,enum_croSSbeam,2500<BR>r,enum_croSSbeam,0.5,0.5**3/12,0.5/12,0.5,1,<BR>lsel,u,real,,1,150<BR>*do,i,1,13<BR>lstr,kp_numfour+i*14+400+12,kp_numfour+i*14+12<BR>lstr,kp_numfour+i*14+12,kp_numfour+14*i+800+12<BR>lesize,all,,,3,,1<BR>latt,enum_croSSbeam,enum_croSSbeam,enum_croSSbeam,<BR>lmesh,all<BR>*enddo<BR>allsel,all</P>
<P>dk,kp_numthree,all,<BR>dk,kp_numthree+100,all,<BR>dk,kp_numfour,all,<BR>dk,3001,all,<BR>dk,3002,all,<BR>dk,kp_numfive,uy,</P>
<P>/solu<BR>antype,0<BR>acel,,10<BR>alls<BR>solve<BR>fini<BR>save<BR><BR>来自仿真在线</P>
发表于 2009-2-21 15:21 | 显示全部楼层
:lol :lol :lol :lol 有没有模型啊
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-30 15:19 , Processed in 0.056509 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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