一个斜拉桥的计算空例

一个斜拉桥的计算空例[转帖]!******************2003.10.27-13:00*********************
!*************以shell单元模拟横梁-3-恒载****************
!******************配跨为曲线**************************
/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,,,,,

!******定义单元类型-1

17---桥面-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
志于学而独善其身，行其道而兼济天下