斜拉桥ansys建模Word下载.doc

斜拉桥ansys建模Word下载.doc

《斜拉桥ansys建模Word下载.doc》由会员分享，可在线阅读，更多相关《斜拉桥ansys建模Word下载.doc（10页珍藏版）》请在冰豆网上搜索。

*set,x2,（y1-y2）*tan（alfa1*pi/180）!

桥塔上部的X轴向长度

*set,x1,x2+x3!

桥塔的X轴向长度

*set,kp_yy1,0!

定义桥塔上部的索锚固点竖向距离（从塔顶算起）

*set,kp_yy2,2.5185

*set,kp_yy3,3.5788

*set,kp_yy4,4.6469

*set,kp_yy5,5.7248

*set,kp_yy6,6.8151

*set,kp_yy7,7.9211

*set,kp_yy8,9.0479

*set,kp_yy9,10.2027

*set,kp_yy10,11.3965

*set,kp_yy11,12.6470

*set,kp_yy12,13.9848

*set,kp_yy13,15.7143

*set,kp_yy14,17.7041

*set,kp_yy15,22.0000

k,1,-x1,y1,

k,6,-x1+kp_yy2*tan（alfa1*pi/180）,y1-kp_yy2

k,8,-x1+kp_yy3*tan（alfa1*pi/180）,y1-kp_yy3

k,10,-x1+kp_yy4*tan（alfa1*pi/180）,y1-kp_yy4

k,12,-x1+kp_yy5*tan（alfa1*pi/180）,y1-kp_yy5

k,14,-x1+kp_yy6*tan（alfa1*pi/180）,y1-kp_yy6

k,16,-x1+kp_yy7*tan（alfa1*pi/180）,y1-kp_yy7

k,18,-x1+kp_yy8*tan（alfa1*pi/180）,y1-kp_yy8

k,20,-x1+kp_yy9*tan（alfa1*pi/180）,y1-kp_yy9

k,22,-x1+kp_yy10*tan（alfa1*pi/180）,y1-kp_yy10

k,24,-x1+kp_yy11*tan（alfa1*pi/180）,y1-kp_yy11

k,26,-x1+kp_yy12*tan（alfa1*pi/180）,y1-kp_yy12

k,28,-x1+kp_yy13*tan（alfa1*pi/180）,y1-kp_yy13

k,30,-x1+kp_yy14*tan（alfa1*pi/180）,y1-kp_yy14

k,38,-x1+kp_yy15*tan（alfa1*pi/180）,y1-kp_yy15

kfill,1,6!

在已建关键点内内插关键点

kfill,6,8

kfill,8,10

kfill,10,12

kfill,12,14

kfill,14,16

kfill,16,18

kfill,18,20

kfill,20,22

kfill,22,24

kfill,24,26

kfill,26,28

kfill,28,30

kfill,30,38

*set,kp_numone,38!

定义桥塔上部的最后一个关键点号

/com,defineandmeshtheabovepartoftower

*dim,BBS,,40

*dim,HHS,,40

*dim,SSS,,40

*dim,IIYYS,,40

*dim,IIZZS,,40

*set,length1,2.9546108

*set,length2,4.9251168

*set,width1,3

*set,width2,5

*set,diff1_length,length2-length1!

桥塔上部两个截面的长度之差（纵桥向）

*set,diff1_width,width2-width1!

桥塔上部两个截面的宽度之差（横桥向）

*dim,yy,,40!

定义桥塔上部的关键点竖向距离数组（从塔顶计算）

*do,i,1,38

*set,yy（i）,y1-ky（i）

*set,hhs（i）,diff1_length*yy（i）/（y1-y2）+length1

*set,bbs（i）,diff1_width*yy（i）/（y1-y2）+width1

SSS（i）=BBS（i）*HHS（i）!

按照实心截面

IIZZS（i）=BBS（i）*HHS（i）**3/12

IIyyS（i）=HHS（i）*BBS（i）**3/12

*enddo

*do,i,1,kp_numone-1

et,i,beam44

keyopt,i,9,0

mp,ex,i,3.5e10

mp,prxy,i,0.167

mp,dens,i,2.6e3

mp,alpx,i,1e-05!

定义混凝土的线膨胀系数

r,i,SSs（i）,IIZZS（i）,IIyyS（i）,BBS（i）/2,HHS（i）/2,IIZZS（i）+IIyyS（i）

rmore,SSs（i+1）,IIZZS（i+1）,IIyyS（i+1）,BBS（i+1）/2,HHS（i+1）/2,IIZZS（i+1）+IIyyS（i+1）

rmore,0,0,0,0,0,0

rmore,0,0,BBS（i）/2,HHS（i）/2,BBS（i+1）/2,HHS（i+1）/2

lsel,u,real,,1,kp_numone

lstr,i,i+1

latt,i,i,i

lesize,all,,,1,,1

lmesh,all,all

allsel,all

cm,uptower,elem

*get,emin_ts,elem,0,num,min

*get,emax_ts,elem,0,num,max

allsel

kp_ts=emax_ts-emin_ts+1!

提取出来桥塔上部的关键点数目

/com,createtthekpsoftowerdown

*set,zfirst,1.25!

三个控制点的横桥向坐标

*set,zsecond,7.4839

*set,kp_numtwo,kp_ts+2!

定义桥塔下部的第一个关键点

*set,diff1_elem,67!

桥塔下部的单元数目

*set,kp_numthree,kp_numtwo+diff1_elem!

桥塔下部的最后一个关键点

k,kp_numtwo,-y2/tan（alfa2*pi/180）,y2,-zfirst

k,kp_numthree,,,-zsecond

kfill,kp_numtwo,kp_numthree

*dim,BBX,,200

*dim,HHX,,200

*dim,SSX,,200

*dim,IIYYX,,200

*dim,IIZZX,,200

*set,HHXQ,0.8*sin（alfa2*pi/180）!

定义砍掉部分的长度

*set,BBXQ,0.4*sin（alfa3*pi/180）

*set,SSXQ,HHXQ*BBXQ

*dim,IIYYXQ,,200

*dim,IIZZXQ,,200

!

*set,w4,2.6875*2!

在桥塔折角处单元的宽度

*set,w5,4.3682*2!

在坐标原点处单元的宽度

*set,w4,4.5361558

*set,w5,8.1609208

*do,i,1,diff1_elem+1

BBX（i+kp_numone）=2.5*sin（alfa3*3.1415926/180）

HHX（i+kp_numone）=w4+（i-1）*（w5-w4）/67

SSX（i+kp_numone）=BBX（i+kp_numone）*HHX（i+kp_numone）-SSXQ*2

IIYYXQ（i+kp_numone）=HHXQ*BBXQ**3/12+SSXQ*BBX（i+kp_numone）**2/4

IIZZXQ（i+kp_numone）=BBXQ*HHXQ**3/12+SSXQ*HHX（i+kp_numone）**2/4

IIYYX（i+kp_numone）=HHX（i+kp_numone）*BBX（i+kp_numone）**3/12-IIYYXQ（i+kp_numone）*2

IIZZX（i+kp_numone）=BBX（i+kp_numone）*HHX（i+kp_numone）**3/12-IIZZXQ（i+kp_numone）*2

*do,i,KP_numone,kp_numone+diff1_elem-1

mp,dens,i,2600

r,i,SSX（i+1）,IIZZX（i+1）,IIYYX（i+1）,BBX（i+1）/2,HHX（i+1）/2,IIZZX（i+1）+IIYYX（i+1）

rmore,SSX（i+2）,IIZZX（i+2）,IIYYX（i+2）,BBX（i+2）/2,HHX（i+2）/2,IIZZX（i+2）+IIYYX（i+2）

rmore,0,0,BBX（i+1）/2,HHX（i+1）/2,BBX（i+2）/2,HHX（i+2）/2

*do,i,kp_numtwo,kp_numthree-1

lsel,u,real,,1,150

latt,i-1,i-1,i-1

lmesh,all

esel,u,real,,1,kp_numtwo-2

cm,downtower,elem

cmsel,s,downtower

*get,emin_tx,elem,0,num,min

*get,emax_tx,elem,0,num,max

kp_tx=emax_tx-emin_tx+1

esel,all,all

cm,tower,elem

lsel,s,real,,kp_numone,kp_numthree-2

lsymm,z,all,,,100

*set,kp_numfour,kp_numone+2*kp_tx+100+1!

主梁的第一个关键点号

esel,s,ename,