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

AaronSpark

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. 
    
15. *set,kp_yy1,0 !定义桥塔上部的索锚固点竖向距离（从塔顶算起）
    
16. *set,kp_yy2,2.5185
    
17. *set,kp_yy3,3.5788
    
18. 
    
19. *set,kp_yy4,4.6469
    
20. *set,kp_yy5,5.7248
    
21. *set,kp_yy6,6.8151
    
22. *set,kp_yy7,7.9211
    
23. *set,kp_yy8,9.0479
    
24. *set,kp_yy9,10.2027
    
25. 
    
26. *set,kp_yy10,11.3965
    
27. *set,kp_yy11,12.6470
    
28. *set,kp_yy12,13.9848
    
29. *set,kp_yy13,15.7143
    
30. *set,kp_yy14,17.7041
    
31. *set,kp_yy15,22.0000
    
32. 
    
33. k,1,-x1,y1,
    
34. k,6,-x1+kp_yy2*tan(alfa1*pi/180),y1-kp_yy2
    
35. k,8,-x1+kp_yy3*tan(alfa1*pi/180),y1-kp_yy3
    
36. k,10,-x1+kp_yy4*tan(alfa1*pi/180),y1-kp_yy4
    
37. k,12,-x1+kp_yy5*tan(alfa1*pi/180),y1-kp_yy5
    
38. k,14,-x1+kp_yy6*tan(alfa1*pi/180),y1-kp_yy6
    
39. k,16,-x1+kp_yy7*tan(alfa1*pi/180),y1-kp_yy7
    
40. k,18,-x1+kp_yy8*tan(alfa1*pi/180),y1-kp_yy8
    
41. k,20,-x1+kp_yy9*tan(alfa1*pi/180),y1-kp_yy9
    
42. k,22,-x1+kp_yy10*tan(alfa1*pi/180),y1-kp_yy10
    
43. k,24,-x1+kp_yy11*tan(alfa1*pi/180),y1-kp_yy11
    
44. k,26,-x1+kp_yy12*tan(alfa1*pi/180),y1-kp_yy12
    
45. k,28,-x1+kp_yy13*tan(alfa1*pi/180),y1-kp_yy13
    
46. k,30,-x1+kp_yy14*tan(alfa1*pi/180),y1-kp_yy14
    
47. k,38,-x1+kp_yy15*tan(alfa1*pi/180),y1-kp_yy15
    
48. 
    
49. kfill,1,6 !在已建关键点内内插关键点
    
50. kfill,6,8
    
51. kfill,8,10
    
52. kfill,10,12
    
53. kfill,12,14
    
54. kfill,14,16
    
55. kfill,16,18
    
56. kfill,18,20
    
57. kfill,20,22
    
58. kfill,22,24
    
59. kfill,24,26
    
60. kfill,26,28
    
61. kfill,28,30
    
62. kfill,30,38
    
63. *set,kp_numone,38!定义桥塔上部的最后一个关键点号
    
64. 
    
65. /com,define and mesh the above part of tower
    
66. *dim,BBS,,40
    
67. *dim,HHS,,40
    
68. *dim,SSS,,40
    
69. *dim,IIYYS,,40
    
70. *dim,IIZZS,,40
    
71. *set,length1,2.9546108
    
72. *set,length2,4.9251168
    
73. *set,width1,3
    
74. *set,width2,5
    
75. *set,diff1_length,length2-length1 !桥塔上部两个截面的长度之差（纵桥向）
    
76. *set,diff1_width,width2-width1 !桥塔上部两个截面的宽度之差（横桥向）
    
77. 
    
78. 
    
79. *dim,yy,,40 !定义桥塔上部的关键点竖向距离数组（从塔顶计算）
    
80. *do,i,1,38
    
81. *set,yy(i),y1-ky(i)
    
82. *set,hhs(i),diff1_length*yy(i)/(y1-y2)+length1
    
83. *set,bbs(i),diff1_width*yy(i)/(y1-y2)+width1
    
84. SSS(i)=BBS(i)*HHS(i) !按照实心截面
    
85. IIZZS(i)=BBS(i)*HHS(i)**3/12
    
86. IIyyS(i)=HHS(i)*BBS(i)**3/12
    
87. *enddo
    
88. 
    
89. 
    
90. *do,i,1,kp_numone-1
    
91. et,i,beam44
    
92. keyopt,i,9,0
    
93. mp,ex,i,3.5e10
    
94. mp,prxy,i,0.167
    
95. mp,dens,i,2.6e3
    
96. mp,alpx,i,1e-05 !定义混凝土的线膨胀系数
    
97. r,i,SSs(i),IIZZS(i),IIyyS(i),BBS(i)/2,HHS(i)/2,IIZZS(i)+IIyyS(i)
    
98. rmore,SSs(i+1),IIZZS(i+1),IIyyS(i+1),BBS(i+1)/2,HHS(i+1)/2,IIZZS(i+1)+IIyyS(i+1)
    
99. rmore,0,0,0,0,0,0
    
100. rmore,0,0,BBS(i)/2,HHS(i)/2,BBS(i+1)/2,HHS(i+1)/2
     
101. lsel,u,real,,1,kp_numone
     
102. lstr,i,i+1
     
103. latt,i,i,i
     
104. lesize,all,,,1,,1
     
105. lmesh,all,all
     
106. allsel,all
     
107. *enddo
     
108. 
     
109. cm,uptower,elem
     
110. 
     
111. *get,emin_ts,elem,0,num,min
     
112. *get,emax_ts,elem,0,num,max
     
113. allsel
     
114. kp_ts=emax_ts-emin_ts+1 !提取出来桥塔上部的关键点数目
     
115. 
     
116. /com,createt the kps of tower down
     
117. *set,zfirst,1.25 !三个控制点的横桥向坐标
     
118. *set,zsecond,7.4839
     
119. *set,kp_numtwo,kp_ts+2 !定义桥塔下部的第一个关键点
     
120. *set,diff1_elem,67 !桥塔下部的单元数目
     
121. *set,kp_numthree,kp_numtwo+diff1_elem !桥塔下部的最后一个关键点
     
122. k,kp_numtwo,-y2/tan(alfa2*pi/180),y2 ,-zfirst
     
123. k,kp_numthree,,,-zsecond
     
124. kfill,kp_numtwo,kp_numthree
     
125. 
     
126. *dim,BBX,,200
     
127. *dim,HHX,,200
     
128. *dim,SSX,,200
     
129. *dim,IIYYX,,200
     
130. *dim,IIZZX,,200
     
131. *set,HHXQ,0.8*sin(alfa2*pi/180) !定义砍掉部分的长度
     
132. *set,BBXQ,0.4*sin(alfa3*pi/180)
     
133. *set,SSXQ,HHXQ*BBXQ
     
134. *dim,IIYYXQ,,200
     
135. *dim,IIZZXQ,,200
     
136. !*set,w4,2.6875*2 !在桥塔折角处单元的宽度
     
137. !*set,w5,4.3682*2 !在坐标原点处单元的宽度
     
138. *set,w4,4.5361558
     
139. *set,w5,8.1609208
     
140. 
     
141. *do,i,1,diff1_elem+1
     
142. BBX(i+kp_numone)=2.5*sin(alfa3*3.1415926/180)
     
143. HHX(i+kp_numone)=w4+(i-1)*(w5-w4)/67
     
144. SSX(i+kp_numone)=BBX(i+kp_numone)*HHX(i+kp_numone)-SSXQ*2
     
145. IIYYXQ(i+kp_numone)=HHXQ*BBXQ**3/12+SSXQ*BBX(i+kp_numone)**2/4
     
146. IIZZXQ(i+kp_numone)=BBXQ*HHXQ**3/12+SSXQ*HHX(i+kp_numone)**2/4
     
147. IIYYX(i+kp_numone)=HHX(i+kp_numone)*BBX(i+kp_numone)**3/12-IIYYXQ(i+kp_numone)*2
     
148. IIZZX(i+kp_numone)=BBX(i+kp_numone)*HHX(i+kp_numone)**3/12-IIZZXQ(i+kp_numone)*2
     
149. *enddo
     
150. 
     
151. *do,i,KP_numone,kp_numone+diff1_elem-1
     
152. et,i,beam44
     
153. keyopt,i,9,0
     
154. mp,ex,i,3.5e10
     
155. mp,prxy,i,0.167
     
156. mp,dens,i,2600
     
157. mp,alpx,i,1e-05 !定义混凝土的线膨胀系数
     
158. 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)
     
159. rmore,SSX(i+2),IIZZX(i+2),IIYYX(i+2),BBX(i+2)/2,HHX(i+2)/2,IIZZX(i+2)+IIYYX(i+2)
     
160. rmore,0,0,0,0,0,0
     
161. rmore,0,0,BBX(i+1)/2,HHX(i+1)/2,BBX(i+2)/2,HHX(i+2)/2
     
162. *enddo
     
163. *do,i,kp_numtwo,kp_numthree-1
     
164. lsel,u,real,,1,150
     
165. lstr,i,i+1
     
166. latt,i-1,i-1,i-1
     
167. lesize,all,,,1,,1
     
168. lmesh,all
     
169. allsel,all
     
170. *enddo
     
171. 
     
172. esel,u,real,,1,kp_numtwo-2
     
173. cm,downtower,elem
     
174. allsel,all
     
175. cmsel,s,downtower
     
176. *get,emin_tx,elem,0,num,min
     
177. *get,emax_tx,elem,0,num,max
     
178. allsel
     
179. kp_tx=emax_tx-emin_tx+1
     
180. 
     
181. esel,all,all
     
182. cm,tower,elem
     
183. 
     
184. lsel,s,real,,kp_numone,kp_numthree-2
     
185. lsymm,z,all,,,100
     
186. allsel,all
     
187. *set,kp_numfour,kp_numone+2*kp_tx+100+1 !主梁的第一个关键点号
     
188. !esel,s,ename,,beam44
     
189. !tunif,0
     
190. !tref,-30
     
191. !allsel,all
     
192. 
     
193. 
     
194. /com,couple the tower upside and tower down
     
195. cerig,node(kx(38),ky(38),kz(38)),node(kx(39),ky(39),kz(39)),all,
     
196. cerig,node(kx(38),ky(38),kz(38)),node(kx(139),ky(139),kz(139)),all,
     
197. 
     
198. *set,beam_height,1.2725-0.3　 !主梁节点即锚固点到原点的距离
     
199. !*set,kp_numfour,kp_numone+2*(kp_tx+1)+100+1 !主梁的第一个关键点号
     
200. *set,kp_numfour_inc,210 !主梁的关键点数目-1
     
201. *set,kp_numfive,kp_numfour+kp_numfour_inc
     
202. k,kp_numfour,,beam_height
     
203. k,kp_numfive,105,beam_height
     
204. kfill,kp_numfour,kp_numfive
     
205. *do,i,1,kp_numfour_inc
     
206. lstr,kp_numfour+i-1,kp_numfour+i
     
207. *enddo
     
208. *set,enum_beam,emax_tx+1
     
209. et,enum_beam,beam188
     
210. mp,ex,enum_beam,3.5e10
     
211. mp,prxy,enum_beam,0.0.167
     
212. mp,dens,enum_beam,3038.8 !考虑到二期恒载后的换算密度
     
213. mp,alpx,enum_beam,1e-05 !定义混凝土的线膨胀系数
     
214. KEYOPT,enum_beam,7,2
     
215. keyopt,enum_beam,8,3
     
216. keyopt,enum_beam,9,3
     
217. SECTYPE,1,BEAM,MESH,sec1
     
218. SECOFFSET,user,,-1.40 !截面读入时主梁的平移
     
219. SECREAD,'main_beam','SECT',' ',MESH
     
220. k,5000,,1000000
     
221. lsel,u,real,,1,200
     
222. latt,enum_beam,enum_beam,enum_beam,,5000,,1
     
223. lesize,all,,,1,,1
     
224. lmesh,all
     
225. allsel,all
     
226. esel,s,ename,,beam188
     
227. cm,main_beam,elem
     
228. allsel,all
     
229. 
     
230. /com,create the cable element
     
231. *set,enum_link,enum_beam+1 !拉索的开始单元号
     
232. *dim,cable_area,,13 !定义拉索单元的面积数组
     
233. *dim,cable_istrain,,13 !定义拉索单元的初始应变数组
     
234. *dim,cable_dens,,13
     
235. 
     
236. /com,define the angle of all cable
     
237. *set,cable_area1,1.668E-03
     
238. *set,cable_area2,1.668E-03
     
239. *set,cable_area3,2.6410E-03
     
240. *set,cable_area4,2.6410E-03
     
241. *set,cable_area5,2.6410E-03
     
242. *set,cable_area6,2.6410E-03
     
243. *set,cable_area7,3.0580E-03
     
244. *set,cable_area8,3.0580E-03
     
245. *set,cable_area9,3.0580E-03
     
246. *set,cable_area10,3.7530E-03
     
247. *set,cable_area11,3.7530E-03
     
248. *set,cable_area12,3.7530E-03
     
249. *set,cable_area13,3.7530E-03
     
250. *set,cable_area_back,2.0155E-02
     
251. 
     
252. 
     
253. *set,cable_area(1),cable_area1
     
254. *set,cable_area(2),cable_area2
     
255. *set,cable_area(3),cable_area3
     
256. *set,cable_area(4),cable_area4
     
257. *set,cable_area(5),cable_area5
     
258. *set,cable_area(6),cable_area6
     
259. *set,cable_area(7),cable_area7
     
260. *set,cable_area(8),cable_area8
     
261. *set,cable_area(9),cable_area9
     
262. *set,cable_area(10),cable_area10
     
263. *set,cable_area(11),cable_area11
     
264. *set,cable_area(12),cable_area12
     
265. *set,cable_area(13),cable_area13
     
266. 
     
267. *set,cable_dens(1),13.2/cable_area1
     
268. *set,cable_dens(2),13.2/cable_area2
     
269. *set,cable_dens(3),20.9/cable_area3
     
270. *set,cable_dens(4),20.9/cable_area4
     
271. *set,cable_dens(5),20.9/cable_area5
     
272. *set,cable_dens(6),20.9/cable_area6
     
273. *set,cable_dens(7),24.2/cable_area7
     
274. *set,cable_dens(8),24.2/cable_area8
     
275. *set,cable_dens(9),24.2/cable_area9
     
276. *set,cable_dens(10),29.7/cable_area10
     
277. *set,cable_dens(11),29.7/cable_area11
     
278. *set,cable_dens(12),29.7/cable_area12
     
279. *set,cable_dens(13),29.7/cable_area13
     
280. *set,cable_dens_back,159.5/cable_area_back
     
281. 
     
282. *set,cable_istrain1,0.26032E-02
     
283. *set,cable_istrain2,0.25568E-02
     
284. 
     
285. *set,cable_istrain3,0.23210E-02
     
286. 
     
287. *set,cable_istrain4,0.23456E-02
     
288. *set,cable_istrain5,0.23892E-02
     
289. *set,cable_istrain6,0.24412E-02
     
290. 
     
291. *set,cable_istrain7,0.28199E-02
     
292. *set,cable_istrain8,0.28719E-02
     
293. *set,cable_istrain9,0.29143E-02
     
294. 
     
295. *set,cable_istrain10,0.28321E-02
     
296. *set,cable_istrain11,0.28559E-02
     
297. *set,cable_istrain12,0.28743E-02
     
298. *set,cable_istrain13,0.28926E-02
     
299. 
     
300. cable_back_istrain1=0.32891E-02
     
301. cable_back_istrain2=0.33661E-02
     
302. 
     
303. *set,cable_istrain(1),cable_istrain1
     
304. *set,cable_istrain(2),cable_istrain2
     
305. *set,cable_istrain(3),cable_istrain3
     
306. *set,cable_istrain(4),cable_istrain4
     
307. *set,cable_istrain(5),cable_istrain5
     
308. *set,cable_istrain(6),cable_istrain6
     
309. *set,cable_istrain(7),cable_istrain7
     
310. *set,cable_istrain(8),cable_istrain8
     
311. *set,cable_istrain(9),cable_istrain9
     
312. *set,cable_istrain(10),cable_istrain10
     
313. *set,cable_istrain(11),cable_istrain11
     
314. *set,cable_istrain(12),cable_istrain12
     
315. *set,cable_istrain(13),cable_istrain13

复制代码

AaronSpark

<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>

haobaobaos

:lol :lol :lol :lol 有没有模型啊