|
/PREP7
ET,1,BEAM4
et,2,link10 !缆索
et,3,mass21,,,0 !附加质量惯距
!et,4,shell63 !板壳单元
!et,5,mass21,,,0 !附加质量惯距
KEYOPT,1,6,1
keyopt,1,2,1
keyopt,2,3,0 !只受拉
keyopt,2,2,2 !很小横向刚度
et,6,link10!缆索
!* 材料
MP,EX,1,3.5e10
MP,PRXY,1,0.167 ! 塔/柱/加劲梁/横隔板
mp,dens,1,2600
mp,ALPX,1,0.00001
MP,EX,2,2.0e11 ! 主缆、吊索、斜索
MP,PRXY,2,0.3
mp,dens,2,7850
mp,ALPX,2,0.000012
MP,EX,3,2.0e14 ! 钢臂
MP,PRXY,3,0.000002
MP,EX,4,3.5e10
MP,PRXY,4,0.167 ! 主梁
mp,dens,4,2600
mp,ALPX,4,0.00001
MP,EX,5,3.5e10
MP,PRXY,5,0.167 ! 主梁
mp,dens,5,2600
mp,ALPX,5,0.00001
!主梁
!r,1,0.39,0.0137,0.0117,0.6,0.65,,,0.0216
!r,1,5.1788,3.837,93.1022,14.5,2.23
r,1,7.903,3.837,93.1022,14.5,2.23
rmore,,7.415
!横隔板刚臂 +加劲肋
!r,2,0,1.0e3,1.0e3,0.9,0.9,0,
!rmore,0,1.0e3,,,,
r,2,0.769,0.225,0.00468,13.3,1.98
rmore,0,0,0,0,0,!
!柱
r,3,4*2,2*4*4*4/12,4*2*2*2/12,2,4,,
rmore,0,0,0,0,0,!tower_up
!塔
r,4,5*3,3*5*5*5/12,5*3*3*3/12,3,5,,
rmore,0,,0,0,0,!tower_up
!柱横梁
r,5,3*1.5,3*1.5*1.5*1.5/12,1.5*3*3*3/12,3,1.5,,
rmore,0,0,0,0,0,!tower_up
!塔横梁
r,6,3*2,3*2*2*2/12,2*3*3*3/12,3,2,,
rmore,0,0,0,0,0,!tower_up
!索鞍&索梁刚臂
r,7,1000,1.0e3,1.0e3,0.9,0.9,0,
rmore,0,1.0e3,,,,!刚愎
r,8,19*2489.2/1000000,0.0046! 0.00323! !主缆1# 跨中c=40600300/(19*137.2)/2e5
!!!!!斜拉索(end)
r,51,19*137.2/1000000,0.0025!1156000/(19*137.2)/2e5 !1# 中跨 623
r,52,19*137.2/1000000,0.0025!934000/(19*137.2)/2e5 !2# 627
r,53,19*137.2/1000000,0.0025!1124000/(19*137.2)/2e5 !3# 631
r,54,19*137.2/1000000,0.0025!1052000/(19*137.2)/2e5 !4# 635
r,55,19*137.2/1000000,0.0025!1105000/(19*137.2)/2e5 !5# 639
r,56,19*137.2/1000000,0.0025!1185000/(19*137.2)/2e5 !6# 643
r,57,19*137.2/1000000,0.0025!1235000/(19*137.2)/2e5 !7# 647
r,58,19*137.2/1000000,0.0025!1251000/(19*137.2)/2e5 !8# 651
r,59,19*137.2/1000000,0.0025!1330000/(19*137.2)/2e5 !9# 655
r,60,19*137.2/1000000,0.0025!1403000/(19*137.2)/2e5 !10# 659
r,71,19*137.2/1000000,0.0025!1382000/(19*137.2)/2e5 !1# 边跨 663
r,72,19*137.2/1000000,0.0025!1347000/(19*137.2)/2e5 !2# 667
r,73,19*137.2/1000000,0.0025!1300000/(19*137.2)/2e5 !3# 671
r,74,19*137.2/1000000,0.0025!1223000/(19*137.2)/2e5 !4# 675
r,75,19*137.2/1000000,0.0025!1196000/(19*137.2)/2e5 !5# 679
r,76,19*137.2/1000000,0.0025!1120000/(19*137.2)/2e5 !6# 683
r,77,19*137.2/1000000,0.0025!1040000/(19*137.2)/2e5 !7# 687
r,78,19*137.2/1000000,0.0025!983000/(19*137.2)/2e5 !8# 691
r,79,19*137.2/1000000,0.0025!875000/(19*137.2)/2e5 !9# 695
r,80,19*137.2/1000000,0.0025!1176000/(19*137.2)/2e5 !10# 699
!!!!!板单元
!r,61,0.14,0.14,0.14,0.14,
!r,62,0.10,0.10,0.10,0.10,
!r,63,0.12,0.12,0.12,0.12,
!!!!!吊杆(1#~15#)
!r,101,0.000001,0.001646! 1#
r,102,0.00167,0.0! 2#
r,103,0.00167,0.00! 3#
r,104,0.00167,0.00! 4#
r,105,0.00167,0.00! 5#
r,106,0.00167,0.0! 6#
r,107,0.00167,0.0! 7#
r,108,0.00167,0.0! 8#
r,109,0.00167,0.0! 9#
r,110,0.00167,0.0! 10#
r,111,0.00167,0.0! 11#
r,112,0.00167,0.0! 12#
r,113,0.00167,0.0! 13#
r,114,0.00167,0.0! 14#
r,115,0.00167,0.0! 15#
!中跨主缆
q1=371.4
H=49241.7
b=H/q1*2.75
a=288/(2*b)
k=(exp(a)+exp(-a))/2
*do,i,1,145
xi=2*(i-1)
yi=-b*(k-(exp(1/b*xi-a)+exp(-1/b*xi+a))/2)+30.308
n,i,xi,yi,-7
n,i+1000,xi,yi,7
*enddo
!加劲梁节点
n,2001,-219.1+144,0,0
*do,i,1,71
n,i+2001,-210+6*(i-1)+144,0,0
*enddo
n,2073,219.1+144,0,0
!左横梁节点
n,3001,-219.1+144,0,-7
*do,i,1,71
n,i+3001,-210+6*(i-1)+144,0,-7
*enddo
n,3073,219.1+144,0,-7
!右横梁节点
n,4001,-219.1+144,0,7
*do,i,1,71
n,i+4001,-210+6*(i-1)+144,0,7
*enddo
n,4073,219.1+144,0,7
!塔节点
n,8001,0,-26,-7
n,9001,0,-26,7
n,10001,288,-26,-7
n,11001,288,-26,7
*do,i,1,13
n,i+8001,0,-26+2*i,-7
n,i+9001,0,-26+2*i,7
n,i+10001,288,-26+2*i,-7
n,i+11001,288,-26+2*i,7
*enddo
!柱节点
*do,i,1,14
n,i+8014,0,2*i,-7
n,i+9014,0,2*i,7
n,i+10014,288,2*i,-7
n,i+11014,288,2*i,7
*enddo
n,8029,0,30.308,-7
n,9029,0,30.308,7
n,10029,288,30.308,-7
n,11029,288,30.308,7
!生成单元
!!!!主梁单元
type,1
mat,5
real,1
*do,i,1,22
e,2001+i,i+2002
*enddo
type,1
mat,5
real,1
*do,i,23,48
e,2001+i,i+2002
*enddo
type,1
mat,5
real,1
*do,i,49,70
e,2001+i,i+2002
*enddo
!!!!横隔板单元
type,1
mat,3
real,2
*do,i,1,71
e,2001+i,3001+i
e,2001+i,4001+i
*enddo
!!!!塔单元
type,1
mat,1
real,4
*do,i,1,13
e,8000+i,8001+i
e,9000+i,9001+i
e,10000+i,10001+i
e,11000+i,11001+i
*enddo
!!!!柱单元
type,1
mat,1
real,3
*do,i,1,15
e,8013+i,8014+i
e,9013+i,9014+i
e,10013+i,10014+i
e,11013+i,11014+i
*enddo
!!!!塔横梁单元
type,1
mat,1
real,6
n,8050,0,-4,0
n,10050,288,-4,0
e,8012,8050
e,8050,9012
e,10012,10050
e,10050,11012
!!!!柱横梁单元
type,1
mat,1
real,5
e,8026,9026
e,10026,11026
!主缆单元
n,8030,0,30.308,-7
n,9030,0,30.308,7
n,10030,288,30.308,-7
n,11030,288,30.308,7
type,2
mat,2
real,8
e,3001,8030
e,4001,9030
e,3073,10030
e,4073,11030
*do,i,1,144
e,i,i+1
e,1000+i,1001+i
*enddo
!斜拉索单元（中跨1#~10#）
*do,i,1,10
type,2
mat,2
real,50+i
e,3024-i,8029-i
e,4024-i,9029-i
e,3050+i,10029-i
e,4050+i,11029-i
*enddo
!斜拉索单元（边跨1#~10#）
*do,i,1,10
type,2
mat,2
real,70+i
e,3002+i,8029-i
e,4002+i,9029-i
e,3072-i,10029-i
e,4072-i,11029-i
*enddo
!吊杆单元（1#~10#）
*do,i,2,14
type,2
mat,2
real,100+i
e,3022+i,28+3*i
e,3052-i,118-3*i
e,4022+i,1028+3*i
e,4052-i,1118-3*i
*enddo
e,73,3037
e,1073,4037
!!!!塔顶索鞍与塔耦合
cp,next,ux,1,8030
!cp,next,ux,8029,1,8030
cp,next,uy,8029,1,8030
cp,next,uz,8029,1,8030
cp,next,ux,1001,9030
!cp,next,ux,9029,1001,9030
cp,next,uy,9029,1001,9030
cp,next,uz,9029,1001,9030
cp,next,ux,145,10030
!cp,next,ux,10029,145,10030
cp,next,uy,10029,145,10030
cp,next,uz,10029,145,10030
cp,next,ux,1145,11030
!cp,next,ux,11029,1145,11030
cp,next,uy,11029,1145,11030
cp,next,uz,11029,1145,11030
allsel
/solu
D,3001,all !塔，墩底固结
D,3073,all
D,4001,all
D,4073,all
D,8001,all
D,9001,all
D,10001,all
D,11001,all
D,3002,ux
D,3002,uy
D,3002,uz
D,4002,ux
D,4002,uy
D,4002,uz
D,3072,ux
D,3072,uy
D,3072,uz
D,4072,ux
D,4072,uy
D,4072,uz
allsel
f,34,fy,-2700000
f,1034,fy,-2700000
f,112,fy,-2700000
f,1112,fy,-2700000
f,35,fy,-300000
f,1035,fy,-300000
f,111,fy,-300000
f,1111,fy,-300000
*do,i,3,25
f,34+3*(i-1),fy,110000
f,1034+3*(i-1),fy,110000
f,34+3*(i-1),fy,110000
f,1034+3*(i-1),fy,110000
*enddo
/solu
acel,,9.8 !施加重力加速度
allsel
esel,s,,,1,70
!sfbeam,all,2,pres,26670,26670
!sfbeam,all,2,pres,26670+500+10530*3,26670+500+10530*3
sfbeam,all,2,pres,26670+500,26670+500
!sfbeam,all,2,pres,10530*3,10530*3
allsel
*do,i,28,46,1
f,2000+i,fy,-63180*3
*enddo
f,2037,fy,-(360000)*3
allsel
eqslv,-1
nsubst,2
nlgeom,on
nropt,full
sstif,on
allsel
allsel
solve |
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
显身卡
