.
斜拉桥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 *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
*enddo
*do,i,KP_numone,kp_numone+diff1_elem-1 et,i,beam44 keyopt,i,9,0 mp,ex,i,3.5e10 mp,prxy,i,0.167 mp,dens,i,2600
mp,alpx,i,1e-05 !定义混凝土的线膨胀系数
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,0,0,0,0
rmore,0,0,BBX(i+1)/2,HHX(i+1)/2,BBX(i+2)/2,HHX(i+2)/2 *enddo
*do,i,kp_numtwo,kp_numthree-1 lsel,u,real,,1,150 lstr,i,i+1
latt,i-1,i-1,i-1 lesize,all,,,1,,1
精品
.
lmesh,all allsel,all
精品
斜拉桥ansys建模
