[转帖]斜拉桥的算例

AaronSpark

<P>/prep7</P>
<P>!*******************定义节点**********************</P>
<P>!*********************塔**************************<BR>!*******************塔上部************************<BR>*do,i,54.5,33.5,-0.5<BR>  k,110-i*2,-(82.1495+i)/7.403322,i,0<BR>*enddo</P>
<P>!*******************塔下部************************<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</P>
<P><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</P>
<P>!*********桥面(配跨为曲线)******</P>
<P>k,192,-2,0,-6<BR>k,193,0,0,-6<BR>k,194,2,0,-6<BR>k,195,9,0,-6</P>
<P>*do,i,196,205<BR>  k,i,(i-195)*9+9,0,-6<BR>*enddo<BR>k,206,105,0,-6</P>
<P><BR>k,207,-2,0,6<BR>k,208,0,0,6<BR>k,209,2,0,6<BR>k,210,9,0,6</P>
<P>*do,i,211,220<BR>  k,i,(i-210)*9+9,0,6<BR>*enddo<BR>k,221,105,0,6</P>
<P>!****************背索锚碇*******************<BR>k,222,-31,0,14<BR>k,223,-31,0,-14</P>
<P><BR>!***********************定义单元**********************<BR>!上部塔的截面几何特性<BR>*dim,BBS,,43    <BR>*dim,HHS,,43<BR>*dim,SSS,,43<BR>*dim,IIYYS,,43<BR>*dim,IIZZS,,43</P>
<P>*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</P>
<P>!**********平均**********<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</P>
<P>!下部塔的截面几何特性<BR>*dim,BBX,,80<BR>*dim,HHX,,80<BR>*dim,SSX,,80<BR>*dim,IIYYX,,80<BR>*dim,IIZZX,,80</P>
<P>*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</P>
<P>!**********平均**********<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</P>
<P>!******定义单元类型-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</P>
<P>!******定义单元类型-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</P>
<P><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,,,,,</P>
<P>!******定义单元类型-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,,,,,</P>
<P>!******定义单元类型-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,,,,,</P>
<P>!******定义单元类型-118-119---背索******<BR>*dim,cable_strain_bs,,2</P>
<P>cable_strain_bs(1)=3.978873577e-3    !2500<BR>cable_strain_bs(2)=3.978873577e-3    !2500</P>
<P>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)</P>
<P>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)</P>
<P><BR>!******定义单元类型-120-129---主索******<BR>*dim,cable_strain_zs,,11</P>
<P>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</P>
<P>*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</P>
<P><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</P>
<P><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</P>
<P><BR>!*************************塔下部-beam4-2***************</P>
<P>!*************************近朝*************************<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</P>
<P>!*************************远朝*************************<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</P>
<P>*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</P>
<P><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</P>
<P>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</P>
<P>!*******************主索*****************************<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</P>
<P>!****************定义桥面与塔的连接*************<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</P>
<P>!******************定义桥面-area-4********************</P>
<P>wproat,0,-90,0<BR>cyl4,-5.6695,-70,76,90,64,118.09989655,0<BR>wproat,0,90,0</P>
<P>a,188,192,207,191</P>
<P>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</P>
<P>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</P>
<P>*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</P>
<P><BR>!*****************活载****************<BR>!*****不计冲击力-验算最大竖向挠度*****<BR>!sfa,all,1,PRES,1834<BR>!f,196,fy,-114286<BR>!f,197,fy,-114286<BR>!***************************************</P>
<P><BR>f,196,fy,-160000<BR>f,197,fy,-160000</P>
<P>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</P>
<P>!***********重力加速度**********<BR>acel,,12<BR>sfa,all,1,PRES,2567<BR>!acel,,10<BR>solve<BR>finish</P>
<P><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></P>

cchenkaik

qiang!

tiger1980

厉害！！！

sysh320

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