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