[转帖]斜拉桥的算例

AaronSpark · 发表于 2006-5-22 05:16

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

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

x

/prep7
!*******************定义节点**********************
!*********************塔************************** !*******************塔上部************************ *do,i,54.5,33.5,-0.5 k,110-i*2,-(82.1495+i)/7.403322,i,0 *enddo
!*******************塔下部************************ !*远朝 !*上 *do,i,33,0,-0.5 k,110-i*2,-15.6213*i/33.5,i,7.4839*i/33.5-7.4839 *enddo !*下 *do,i,-0.5,-2.5,-0.5 k,110-i*2,-2.7978*i/6,i,-i/2-7.4839 *enddo
 !*近朝 !*上 *do,i,33,0,-0.5 k,182-i*2,-15.6213*i/33.5,i,-7.4839*i/33.5+7.4839 *enddo !*下 *do,i,-0.5,-2.5,-0.5 k,182-i*2,-2.7978*i/6,i,i/2+7.4839 *enddo
!*********桥面(配跨为曲线)******
k,192,-2,0,-6 k,193,0,0,-6 k,194,2,0,-6 k,195,9,0,-6
*do,i,196,205 k,i,(i-195)*9+9,0,-6 *enddo k,206,105,0,-6
 k,207,-2,0,6 k,208,0,0,6 k,209,2,0,6 k,210,9,0,6
*do,i,211,220 k,i,(i-210)*9+9,0,6 *enddo k,221,105,0,6
!****************背索锚碇******************* k,222,-31,0,14 k,223,-31,0,-14
 !***********************定义单元********************** !上部塔的截面几何特性 *dim,BBS,,43 *dim,HHS,,43 *dim,SSS,,43 *dim,IIYYS,,43 *dim,IIZZS,,43
*do,i,1,43 BBS(i)=1.5*(i-1)/42+3 HHS(i)=(i-1)/21+3 SSS(i)=BBS(i)*HHS(i) IIzzS(i)=BBS(i)*HHS(i)**3/12 IIyyS(i)=HHS(i)*BBS(i)**3/12 *enddo
!**********平均********** *dim,BS,,42 *dim,HS,,42 *dim,SS,,42 *dim,IYYS,,42 *dim,IZZS,,42 *do,i,1,42 hS(i)=(hhS(i)+hhS(i+1))/2 bS(i)=(bbS(i)+bbS(i+1))/2 SS(i)=(SSS(i)+SSS(i+1))/2 IYYS(i)=(IIYYS(i)+IIYYS(i+1))/2 IZZS(i)=(IIZZS(i)+IIZZS(i+1))/2 *enddo
!下部塔的截面几何特性 *dim,BBX,,80 *dim,HHX,,80 *dim,SSX,,80 *dim,IIYYX,,80 *dim,IIZZX,,80
*do,i,1,80 BBX(i)=2.5 HHX(i)=(i-1)*4.7164/79+5 SSX(i)=BBX(i)*HHX(i) IIzzX(i)=BBX(i)*HHX(i)**3/12 IIyyX(i)=HHX(i)*BBX(i)**3/12 *enddo
!**********平均********** *dim,BX,,79 *dim,HX,,79 *dim,SX,,79 *dim,IYYX,,79 *dim,IZZX,,79 *do,i,1,79 BX(i)=(BBX(i)+BBX(i+1))/2 HX(i)=(HHX(i)+HHX(i+1))/2 SX(i)=(SSX(i)+SSX(i+1))/2 IYYX(i)=(IIYYX(i)+IIYYX(i+1))/2 IZZX(i)=(IIZZX(i)+IIZZX(i+1))/2 *enddo
!******定义单元类型-1-42---塔上部****** *do,i,1,42 et,i,beam4 mp,ex,i,3.4e10 MP,DENS,i,2.5e3 mp,prxy,i,0.3 r,i,SS(i),IZZS(i),IYYS(i),bS(i),hS(i) *enddo
!******定义单元类型-43-115---塔下部****** *do,i,43,114 et,i,beam4 mp,ex,i,3.4e10 MP,DENS,i,2.5e3 mp,prxy,i,0.3 r,i,SX(i-42),IZZX(i-42),IYYX(i-42),bX(i-42),hX(i-42) *enddo
 !******定义单元类型-115---桥面-1****** et,115,shell63 MP,EX,115,3.4e10 MP,DENS,115,2.5e3 mp,prxy,115,0.3 r,115,1,,,,,
!******定义单元类型-116---桥面-2****** et,116,shell63 MP,EX,116,3.4e10 MP,DENS,116,2.5e3 mp,prxy,116,0.3 r,116,4,,,,,
!******定义单元类型-117---桥面-3****** et,117,shell63 MP,EX,117,3.4e10 MP,DENS,117,2.5e3 mp,prxy,117,0.3 r,117,1,,,,,
!******定义单元类型-118-119---背索****** *dim,cable_strain_bs,,2
cable_strain_bs(1)=3.978873577e-3 !2500 cable_strain_bs(2)=3.978873577e-3 !2500
et,118,link10 MP,EX,118,2e11 MP,DENS,118,7.85e3 mp,prxy,118,0.3 R,118,0.031416,cable_strain_bs(1)
et,119,link10 MP,EX,119,2e11 MP,DENS,119,7.85e3 mp,prxy,119,0.3 R,119,0.031416,cable_strain_bs(2)
 !******定义单元类型-120-129---主索****** *dim,cable_strain_zs,,11
cable_strain_zs(1)=2.419155135e-3 !380 cable_strain_zs(2)=2.419155135e-3 !380 cable_strain_zs(3)=2.2918311805e-3 !360 cable_strain_zs(4)=2.2918311805e-3 !360 cable_strain_zs(5)=2.228169203e-3 !350 cable_strain_zs(6)=2.228169203e-3 !350 cable_strain_zs(7)=1.909859317e-3 !300 cable_strain_zs(8)=1.65521141e-3 !260 cable_strain_zs(9)=1.65521141e-3 !260 cable_strain_zs(10)=1.40056340e-3 !220 cable_strain_zs(11)=1.273239545e-3 !200
*do,i,120,130 et,i,link10 MP,EX,i,2e11 MP,DENS,i,7.85e3 mp,prxy,i,0.3 R,i,0.007853982,cable_strain_zs(i-119) !定义索面积直径=0.1米，初始应变 *enddo
 !********定义桥面与塔的连接*********** et,131,beam4 mp,ex,131,3.4e10 MP,DENS,131,2.5e3 mp,prxy,131,0.3 r,131,16,21.3333,21.3333,4,4
 !*************************定义塔上部-beam4-1************** *do,i,1,42 l,i,i+1 lsel,u,type,,1,131 latt,i,i,i,,,, lesize,all,,,1,,1 lmesh,all,all alls *enddo
 !*************************塔下部-beam4-2***************
!*************************近朝************************* *do,i,43,114 l,i,i+1 lsel,u,type,,1,131 latt,i,i,i,,,, lesize,all,,,1,,1 lmesh,all,all alls *enddo
!*************************远朝************************* l,43,116 lsel,u,type,,1,131 latt,43,43,43,,,, lesize,all,,,1,,1 lmesh,all,all alls
*do,i,116,186 l,i,i+1 lsel,u,type,,1,131 latt,i-72,i-72,i-72,,,, lesize,all,,,1,,1 lmesh,all,all alls *enddo
 !*******************定义索-link10******************** !*******************背索***************************** l,9,222 lsel,u,type,,1,131 latt,118,118,118,,,, lesize,all,,,1,,1 lmesh,all,all alls
l,9,223 lsel,u,type,,1,131 latt,119,119,119,,,, lesize,all,,,1,,1 lmesh,all,all alls
!*******************主索***************************** *do,i,1,11 l,i*4+1,206-i l,i*4+1,221-i lsel,u,type,,1,131 latt,i+119,i+119,i+119,,,, lesize,all,,,1,,1 lmesh,all,all alls *enddo
!****************定义桥面与塔的连接************* l,110,193 l,182,208 lsel,u,type,,1,131 latt,131,131,131,,,, lesize,all,,,1,,1 lmesh,all,all alls
!******************定义桥面-area-4********************
wproat,0,-90,0 cyl4,-5.6695,-70,76,90,64,118.09989655,0 wproat,0,90,0
a,188,192,207,191
lsel,u,type,,1,131 aatt,115,115,115,,,, !lesize,all,2,,,,2 lesize,213,2,,,,2 lesize,214,2,,,,2 lesize,215,2,,,,2 lesize,216,2,,,,2 lesize,217,,,,4,2 lesize,218,2,,,,2 lesize,219,,,,4,2 alls amesh,all
a,192,193,208,207 a,193,194,209,208 lsel,u,type,,1,131 aatt,116,116,116,,,, lesize,all,2,,,,1 alls amesh,all
*do,i,1,12 a,i+193,i+194,i+209,i+208 *enddo lsel,u,type,,1,131 aatt,117,117,117,,,, lesize,all,2,,,,1 alls amesh,all
 !*****************活载**************** !*****不计冲击力-验算最大竖向挠度***** !sfa,all,1,PRES,1834 !f,196,fy,-114286 !f,197,fy,-114286 !***************************************
 f,196,fy,-160000 f,197,fy,-160000
finish !***********************加约束************************* /SOLU d,115,all d,187,all d,188,all d,189,all d,235,uy d,239,uy d,280,uy d,292,uy d,735,uy d,739,uy
!***********重力加速度********** acel,,12 sfa,all,1,PRES,2567 !acel,,10 solve finish
 /POST1 !进入后处理模块POST1 esel,s,type,,118,130 etable,link-fx,smisc,1 esel,s,type,,1,114 etable,MXI,smisc,4 esel,s,type,,1,114 etable,MXJ,smisc,10 esel,s,type,,1,114 etable,MYI,smisc,5 esel,s,type,,1,114 etable,MYJ,smisc,11 esel,s,type,,1,114 etable,MZI,smisc,6 esel,s,type,,1,114 etable,MZJ,smisc,12 esel,s,type,,all

cchenkaik · 发表于 2006-5-22 22:12

qiang!

tiger1980 · 发表于 2006-8-24 16:35

厉害！！！

sysh320 · 发表于 2006-8-24 20:55

这个资料确实是好啊，对斜拉桥的分析真的是一个好的引导

账号		自动登录	找回密码
密码			我要加入

[转帖]斜拉桥的算例

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

评分