关于斜拉桥ansys建模

斜拉桥ansys建模

/com,new model of linjiang cable_stayed bridge,2004.2.09

/prep7

/title, cable_stayed bridge,author is Sunhang

/com,define the keypoints

*set,alfa1,10 !angle of tower upside

*set,alfa2,65 !angle of tower downside

*set,alfa3,79.04594 !angle of tower with bridge surface

*set,y1,55.5 !桥塔顶面到原点的距离

*set,y2,33.5 !桥塔中部的Y轴向长度

*set,pi,3.1415926

*set,x3,y2/tan(alfa2*pi/180) !桥塔中部的X轴向长度

*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,define and mesh the above part of tower

*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

*enddo

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,createt the kps of tower down

*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