/filname,cable-stayed bridge,1
keyw,pr_struc,1
/prep7
!定义单元类型
et,1,beam4
et,2,link10
!定义材料属性
mp,ex,1,3.5e10
mp,prxy,1,0.17
mp,dens,1,2500
mp,ex,2,10e15
mp,prxy,2,0
mp,dens,2,0
mp,ex,3,1.9e10
mp,prxy,3,0.25
mp,dens,3,1200
mp,damp,3,0.5
!定义实常数
!定义实常数
r,1,25.6,20,546.133,16,1.6
r,2,16,29.417,15.394,3.4,4.7
r,3,54,364.5,162,6,9
r,4,40,213.3,83.3,5,8
r,5,1,1/12,1/12,1,1
r,6,0.012,0.012 !索的
!创建节点和单元
!建立主梁节点
/view,1,1,1,1
/angle,1,270,xm,0
/replot
*do,i,1,59 !此循环用于建立主梁的半跨节点
x=-174*2+(i-1)*6 !最左端x=174*2,x=0左边的节点x坐标值,间距为6
y1=-14 !桥面宽28米,故左边节点为-14
y2=14 !桥面宽28米,故右边节点为-14
n,3*(i-1)+1,x !建立主梁节点 3*(i-1)+1为节点号
n,3*(i-1)+2,x,y1 !以下两行建立桥面两边节点
n,3*i,x,y2 !能想出这种建模命令的绝对是编程高手,哈哈
*enddo !完全可以先建立端部的三个节点,然后用这三个节点在x方向上复制59份,间距为6
!建立主梁单元
type,1
real,1
mat,1
*do,i,1,58,1 !以下循环建立建立桥面中线主梁单元
j=3*(i-1)+1
e,j,j+3
*enddo
!建立鱼刺刚横梁
type,1
real,5
mat,2
*do,i,1,59,1 !以下循环用于建立桥面鱼刺横梁的节点
j=3*(i-1)+1
j1=3*(i-1)+2
j2=3*i
e,j,j1
e,j,j2
*enddo
!建立半跨主塔
i=59*3 !变量用于记录桥面的节点数,即至此已经建立了59*3个节点了,用于指导以后设定节点的编号
n,i+1,-174,-10,-30 !以下两行记录塔脚节点
n,i+2,-174,10,-30
n,i+3,-174,-15 !以下两行用于建立与桥面齐高的主塔节点
n,i+4,-174,15
*do,j,1,5,1 !以下循环用于建立索塔在桥面以上的节点
k=i+4+j
n,k,-174,0,60+(j-1)*18
*enddo
!建立下索塔单元
type,1
real,4
mat,1
e,i+1,i+3 !以下用于建立主塔在桥面以下的两根塔柱单元
e,i+2,i+4
!建立中索塔单元
type,1
real,3
mat,1
e,i+3,i+5 !以下用于建立倒Y分叉点到桥面间的两根塔柱单元
e,i+4,i+5
!建立上索塔单元
type,1
real,2
mat,1
*do,j,1,4,1 !以下用于建立倒Y分叉点以上的塔柱单元
k=i+4+j
e,k,k+1
*enddo
!建立与塔的倒Y分叉点链接的索单元
type,2
real,6
mat,3
e,i+5,89
e,i+5,90
!建立主塔倒Y分叉点以上第一个张拉点连接的索单元
*do,j,1,8,1
!此循环用于建立主塔倒Y分叉点以上第一个张拉点连接的所有索单元,共32个
e,i+6,89+3*j
e,i+6,89-3*j
e,i+6,90+3*j
e,i+6,90-3*j
*enddo
!建立与主塔的其他三个张拉点连接的单元
*do,k,1,3,1
*do,j,1,7,1
e,i+6+k,113+(k-1)*21+3*j !一共有28个索单元连接在每个张拉点上
e,i+6+k,65-(k-1)*21-3*j
e,i+6+k,114+(k-1)*21+3*j
e,i+6+k,66-(k-1)*21-3*j
*enddo
*enddo
!生成全桥模型节点
i=i+9 !记录半跨的所有节点数
nsym,x,i,all !用映射法直接建立另半跨节点
esym,,i,all !用映射法直接建立另半跨单元
nummrg,all !合并所有节点和单元
!建立索塔连接横梁单元
type,1
real,5
mat,2
j=i
i=i-9
n,1000,-174
e,1000,i+3
e,1000,i+4
n,2000,174
e,2000,i+3+j
e,2000,i+4+j
!施加主塔的四个脚上的全约束
nsel,s,loc,z,-30
d,all,all
allsel
!在左桥端施加y,z约束
nsel,s,loc,x,-348 !仅给左端主梁施加约束
nsel,r,loc,y,0
d,all,uy
d,all,uz
allsel
!在右桥端施加y约束
nsel,s,loc,x,348 !仅给右端主梁施加约束
nsel,r,loc,y,0
d,all,uy
allsel
numcmp,all
!施加重力场
acel,,,9.8
!耦合节点,耦合跨中由于对称而重复的单元节点以及两主塔上塔横梁和主梁的重合节点,
cpintf,uy
cpintf,uz
cpintf,rotx
cpintf,rotz