简支梁+连续梁建模+预应力钢筋分析实例-APDL

简支梁实体与预应力钢筋分析实例

/COM, Structural /PREP7

egjx=2e5 !Ey

agjx=140 !单根钢绞线面积 ehnt=4e4 !Eh

xzxs=1.0e-5 !线胀系数 yjl=200000 !定义预加力

et,1,link8 !定义link8单元 et,2,solid95 !定义solid95单元 r,1,agjx !!定义link8单元的面积 r,2 !定义第2种实常数

mp,ex,1,egjx !定义link8单元的弹性模量 mp,prxy,1,0.3 !定义link8单元的泊松系数 mp,alpx,1,1.0e-5 !定义线膨胀系数

mp,ex,2,ehnt !定义solid95单元的弹性模量 mp,prxy,2,0.3 !定义solid95单元的泊松系数 blc4, , ,100,200,3000 !定义梁体 /view,1,1,1,1 !定义ISO查看 /ang,1

vplot !绘制梁体

kwpave,6 !工作平面移动到关键点6 wpoff,-30 !工作平面移动-30mm(X) wprot,0,0,90 !工作平面旋转 vsbw,1 !分割梁体

wpoff,0,0,-40 !工作平面移动-40mm(Z) vsbw,2 !分割梁体

wpoff,0,40 !工作平面移动40mm(Y) wprot,0,90 !工作平面旋转 vsbw,all !分割梁体

wpstyl !关闭工作平面显示 nummrg,all,,,,low !整理 numcmp,all !压缩编号

esize,30 !定义网分时边长控制

lsel,s,,,28,38,10 !定义line28和38为新的选择集 latt,1,1,1 !定义选择集的属性 lmesh,all !对线划分单元

allsel,all !新的选择集为所有的实体 gplot !绘制所有的实体

vsel,s,,,all !定义所有体为选择集

vatt,2,2,2 !定义选择集的属性

mshape,0,3d !将体划分单元的形状定位HEX mshkey,1 !采用MAPPED划分器 vmesh,all !对体进行划分单元 finish /solu

dl,3,,all !对线line7施加约束(UX,UY,UZ) dl,16,,all !对线line31施加约束(UX,UY,UZ) dl,23,,all !对线line23施加约束(UX,UY,UZ) dl,2,,uy !对线line4施加约束(UY) dl,15,,uy !对线line30施加约束(UY) dl,22,,uy !对线line23施加约束(UY) dk,2,,,,,ux,uy !对关键点2约束(UX,UY)

bfl,28,temp,-yjl/(xzxs*egjx*agjx) !对钢绞线施加温度 bfl,38,temp,-yjl/(xzxs*egjx*agjx) !对钢绞线施加温度 solve !求解 finish /post1

plnsol,s,z,0,1 !绘制Z方向的应力 etable,sigi,ls,1 !定义钢筋单元数据表 plls,sigi,sigi,1 !绘制上述应力 !finish

!/exit,nosav

模拟预应力损失的两种方法

方法1：通过体分割的办法建立预应力钢束，然后通过定义降温荷载来模拟预应力效应。这种方法在张立明编的《ALGOR、ANSYS在桥梁工程中的应用方法与实例》中有一个算例。 命令流如下： /prep7

/title,the analysis of the solid beam and prestressed tendon /DIST,1,1.08222638492,1 /REP,FAST

*SET,egjx,2e11 !定义弹性模量 *SET,agjx,60e-4 ! *SET,ehnt,3.5e10 *SET,xzxs,1.2e-5 *SET,yjl,200000

et,1,link8 !定义预应力钢筋单元类型为Link8单元 et,2,solid95 !

r,1,agjx !定义实常数 r,2

mp,ex,1,egjx

mp,prxy,1,0.3 mp,dens,1,7800 mp,dens,2,2500 mp,alpx,1,xzxs mp,ex,2,ehnt

mp,prxy,2,0.3 blc4,,,1,2,30 /VIEW,1,1,1,1 /ANG,1 /AUTO,1

/REP,FAST klist,all kwpave,6

wpoff,-0.3 wprot,0,0,90 vsbw,1

wpoff,0,0,-0.4 vsbw,2 wpoff,0,0.4 wprot,0,90 vsbw,all wpstyl

nummrg,all,,,,low numcmp,all esize,0.3 klist,all llist,all

lsel,s,,,28,38,10 latt,1,1,1 lmesh,all allsel,all gplot vsel,s,,,all vatt,2,2,2 mshape,0,3d mshkey,1

vmesh,all finish /solu dl,3,,all dl,16,,all dl,23,,all dl,2,,uy

dl,15,,uy

dl,22,,uy /AUTO,1

/REP,FAST

bfl,28,temp,-yjl/(xzxs*egjx*agjx) bfl,38,temp,-yjl/(xzxs*egjx*agjx) allsel

acel,,9.8 solve finish /post1 pldisp plnsol,s,1

方法二：体线耦合法建立预应力钢束 命令流如下，并附有命令说明： /prep7

*SET,eg,2e5 !钢束弹模 *SET,ag,140 !钢束面积 *SET,eh,4e4 !混凝土的弹模

*SET,r0,9345 !曲线钢束的半径 *SET,yyl,200000 !钢束的预应力大小

et,1,link8 !定于预应力钢束单元类型为link8 et,2,solid95 !定义混凝土单元类型为solid95 r,1,ag,yyl/eg/ag*1.036258 !定义单元实常数，对于link8单元为钢束面积、以及钢束的初始应变（除模拟预应力张拉还会用于应力刚度矩阵的计算） r,2 !solid95的单元实常数仅为x轴的调整，并且该实常数仅仅在单元坐标系统选择为1的时候填写。一般不填。

mp,ex,1,eg !预应力钢束单元性质，包括之前定义的弹模，泊松比 mp,prxy,1,0.3

mp,ex,2,eh !混凝土单元性质，包括弹模、泊松比 mp,prxy,2,0.2

blc4,,,100,200,3000 !定义长方体的长宽高（混凝土块体） /view,1,1,1,1 !调整观察的角度 /ang,1

vplot !显示选择的体 ksel,all !选择所有的关键点 klist

*get,kp0,kp,0,num,max !将最大的kp编号存于kp0中，因为事先不知道最大的关键点号是多少，所以先要存在kp0中

lsel,none !不选择所有的线 llist

k,kp0+1,50,160 !定义关键点，kp0是已有的最大关键点编号。 k,kp0+2,50,160,3000 k,kp0+3,50,800,1500

larc,kp0+1,kp0+2,kp0+3,r0 !将定义的以上三点连成一条圆弧，半径为r0。kp0+3用于定义圆弧所在面及曲线边中点

kdele,kp0+3 !删除kp0+3这个关键点 klist llist

*get,line1,line,0,num,min !将最小的线编号存于line1中，因为line1是最先定义的预应力钢束

lsel,s,loc,z,0 !选择z、y坐标均为0的线，用于施加约束（该线为右侧截面的梁底线） lsel,r,loc,y,0

dl,all,,uy lsel,s,loc,z,3000 y=0的线） lsel,r,loc,y,0

dl,all,,all UZ）

allsel,all lsel,s,,,line1 latt,1,1,1 格化的线联系起来

lesize,all,,,50 lmesh,all vsel,all vlist

vatt,2,2,2 土)

lsel,s,loc,z,0 围内的所有的线 lsel,r,loc,y,10,140

lesize,all,,,8 lsel,s,loc,z,0 到140范围内的所有线 lsel,u,loc,y,10,140

lesize,all,,,4 lsel,s,loc,y,0 根线被选中）注意这根线是一根分网的扫掠线。 lsel,r,loc,x,0

lesize,all,,,50 vsweep,all 填充体 allsel,all lsel,s,,,line1 线）

nsll,s,1 !对线施加y方向的约束 !选择左侧截面梁底线（即z=3000，而 !对选择的线施加约束（UX、UY、 !选择所有对象 !选择线，最小值line1

!将单元性质与选择的、尚未网 !线单元划分的数目

!对line1进行分网 !选择所有的体 !给体赋予相应的属性(混凝 !选择z=0，而y坐标从10到140范 !分网数目8

!再次选择z=0，而不选y坐标从10 !分网数目为4 !选择x=0、y=0坐标处的线（仅有一 !分网数目为50

!通过扫掠网格的方法来用单元

!选择所有的对象

!选择line1最小值（预应力钢绞 !选择与所选择的线对应的所

有的节点，即钢绞线上面的节点，因为网格数为50，所以节点总数为51

cm,cmljnod,node !将钢绞线节点编组，其组名为cmljnod

*get,max1,node,0,count !选择钢绞线中的节点单元数目，将其赋给变量max1，此变量可以作为*do循环的上界来确定循环的次数

*dim,ojd,,max1 !定义数组序列。数组名为ojd，max1为行的范围

*dim,jd,,max1 !定义数组序列。数组名为jd，max1为行的范围

*get,nod1,node,0,num,min !选择集中最小的节点号，赋给变量nod1，作为*do循环的上界

*SET,ojd(1),nod1 !将nod1的值赋给ojd(1) *do,i,2,max1 !do循环，从2到max1 *SET,ojd(i),ndnext(ojd(i-1)) !将ndnext(ojd(i-1))赋给ojd(i)

*enddo !结束循环 allsel,all !选择所有 nsel,all !选择所有节点

cmsel,u,cmljnod !不选择cmljnod，即钢绞线对应的节点

*do,i,1,max1 !i从1循环开始到max1 *SET,nod1,ojd(i) !将ojd(i)赋给nod1

*SET,j,nnear(nod1) !将nnear(nod1)赋给j，是指ojdi邻近的节点吗？nnear()是一个函数，表示选择最近的节点。

*SET,jd(i),j !将j赋给jd(i)

*enddo !结束循环 nsel,all !选择所有节点 *SET,ji,1 !将1赋给ji *do,i,1,max1 !i从1到max1

cp,ji,ux,ojd(i),jd(i) !定义耦合自由度，ji只是随机设置的耦合编号。ux为耦合节点的x方向的位移。 cp,ji+1,uy,ojd(i),jd(i) cp,ji+2,uz,ojd(i),jd(i)

*SET,ji,ji+3 !ji+3赋给ji *enddo !结束循环 allsel,all !选择所有单元

*SET,ji, !将所有的参数删除（第三项value为空格时标识删除） *SET,i,

*SET,max1, *SET,nod1, *SET,ojd, *SET,jd, *SET,j, *SET,ag,

*SET,eg, *SET,eh,

*SET,kp0, *SET,r0,

*SET,yyl, *SET,line1,

finish !建模结束，进入求解 /solu

allsel,all

acel,,9.8 solve FINISH

/POST1 pldisp,1

etable,sigi,ls,1 plls,sigi,sigi,1 plnsol,s,1

连续梁桥ansys 命令流建模求解

!!连续梁桥 /prep7

et,1,4定义梁单元

et,2,21 定义mass21 单元 !!定义粱材料 !!泊松比!!密度

mp,ex,2,3.45e10 !!直线段梁材料 和1M 段梁材料 mp,nuxy,2,0.2

mp,dens,2,3302.153125 mp,ex,3,3.45e10 mp,nuxy,3,0.2

mp,dens,3,3301.658695 mp,ex,4,3.45e10 mp,nuxy,4,0.2

mp,dens,4,3299.906778 mp,ex,5,3.45e10 mp,nuxy,5,0.2

mp,dens,5,3298.327219 mp,ex,6,3.45e10 mp,nuxy,6,0.2

考虑自重 !进入后处理 ! mp,dens,6,3292.351605 mp,ex,7,3.45e10 mp,nuxy,7,0.2

mp,dens,7,3284.137255 mp,ex,8,3.45e10 mp,nuxy,8,0.2

mp,dens,8,3271.802136 mp,ex,9,3.45e10 mp,nuxy,9,0.2

mp,dens,9,3260.41903 mp,ex,10,3.45e10 mp,nuxy,10,0.2

mp,dens,10,3248.193657 mp,ex,11,3.45e10 mp,nuxy,11,0.2

mp,dens,11,3235.1174 mp,ex,12,3.45e10 mp,nuxy,12,0.2

mp,dens,12,3221.5856 mp,ex,13,3.45e10 mp,nuxy,13,0.2

mp,dens,13,3208.826871 mp,ex,14,3.45e10 mp,nuxy,14,0.2

mp,dens,14,3194.279207 mp,ex,15,3.45e10 mp,nuxy,15,0.2

mp,dens,15,3179.924673 mp,ex,16,3.45e10 mp,nuxy,16,0.2

mp,dens,16,3166.445716 mp,ex,17,3.45e10 mp,nuxy,17,0.2

mp,dens,17,3152.555731 mp,ex,18,3.45e10 mp,nuxy,18,0.2

mp,dens,18,3138.312105 mp,ex,19,3.45e10 mp,nuxy,19,0.2

mp,dens,19,3124.795334 mp,ex,20,3.45e10 mp,nuxy,20,0.2 mp,dens,20,3110.7135 mp,ex,21,3.45e10

mp,nuxy,21,0.2

mp,dens,21,3097.080875 mp,ex,22,3.45e10 mp,nuxy,22,0.2

mp,dens,22,3083.186268 mp,ex,23,3.45e10 mp,nuxy,23,0.2

mp,dens,23,3068.968824 mp,ex,24,3.45e10 mp,nuxy,24,0.2

mp,dens,24,3055.612436 mp,ex,25,3.45e10 mp,nuxy,25,0.2

mp,dens,25,3045.857147 mp,ex,26,3.45e10 mp,nuxy,26,0.2

mp,dens,26,3035.174287 mp,ex,27,3.45e10 mp,nuxy,27,0.2

mp,dens,27,3026.696551 mp,ex,28,3.45e10 mp,nuxy,28,0.2

mp,dens,28,3015.795365 mp,ex,29,3.45e10 mp,nuxy,29,0.2

mp,dens,29,3007.710181 mp,ex,30,3.45e10 mp,nuxy,30,0.2

mp,dens,30,3000.513837 mp,ex,31,3.45e10 mp,nuxy,31,0.2

mp,dens,31,2978.611375 mp,ex,32,3.45e10 mp,nuxy,32,0.2

mp,dens,32,2958.618861 mp,ex,33,3.45e10 mp,nuxy,33,0.2

mp,dens,33,2937.888072 mp,ex,34,3.45e10 mp,nuxy,34,0.2

mp,dens,34,2919.475751 mp,ex,35,3.45e10 mp,nuxy,35,0.2

mp,dens,35,2903.359983!!6700 处

mp,ex,36,3.45e10 mp,nuxy,36,0.2

mp,dens,36,3302.153125!!可以不用 mp,ex,37,3.45e10 mp,nuxy,37,0.2

mp,dens,37,3302.153125 mp,ex,38,3.45e10!!-700 处梁 mp,nuxy,38,0.2

mp,dens,38,3180.5701 mp,ex,39,3.45e10!!边支点横隔板 mp,nuxy,39,0.2

mp,dens,39,2868.674818 mp,ex,40,3.45e10 !!合拢段横隔板 mp,nuxy,40,0.2

mp,dens,40,2868.674818 mp,ex,41,3.45e10 !!中支点横隔板 mp,nuxy,41,0.2

mp,dens,41,2757.470588 mp,ex,1,3.25e10 定义墩材料属性 mp,nuxy,1,0.2 mp,dens,1,2650

!!定义实常数 编号 ,面积,IYY,IZZ,宽,高,,RMORE,,抗扭惯距 !!主梁截面

r,1,11.851,28.52,215.151,16.95,4,,!!直线段 rmore,,56.72

r,2,11.862,28.57,215.458,16.95,4.001,, !!截面100 处 rmore,,56.905

r,3,11.871,28.91,215.622,16.95,4.009,, rmore,,57.345

r,4,11.903,29.03,215.78,16.95,4.025,, rmore,,57.938

r,5,11.932,29.626,215.88,16.95,4.049,, rmore,,58.502

r,6,12.043,30.402,217.093,16.95,4.081,, rmore,,59.768

r,7,12.199,31.407,218.9,16.95,4.121,, rmore,,61.445

r,8,12.441,32.716,221.806,16.95,4.169,, rmore,,63.652

r,9,12.673,34.21,224.418,16.95,4.226,, rmore,,66.058

r,10,12.932,35.94,227.279,16.95,4.290,, rmore,,68.801

r,11,13.221,37.941,230.512,16.95,4.362,, rmore,,71.948

r,12,13.534,40.233,233.872,16.95,4.442,, rmore,,75.432

r,13,13.843,42.747,237.133,16.95,4.530,, rmore,,79.173

r,14,14.213,45.688,241.054,16.95,4.627,, rmore,,83.484

r,15,14.598,48.9,245.004,16.95,4.731,, rmore,,88.155

r,16,14.979,52.547,248.817,16.95,4.843,, rmore,,93.111

r,17,15.393,56.562,252.927,16.95,4.963,, rmore,,98.528

r,18,15.842,61.069,257.347,16.95,5.092,, rmore,,104.47

r,19,16.293,65.992,261.673,16.95,5.228,, rmore,,110.742

r,20,16.791,71.484,266.426,16.95,5.372,, rmore,,117.8

r,21,17.303,77.568,271.225,16.95,5.525,, rmore,,124.861

r,22,17.858,84.282,276.363,16.95,5.685,, rmore,,132.631

r,23,18.4,91.768,281.923,16.95,5.854,, rmore,,141.068

r,24,19.072,99.936,287.3,16.95,6.03,, rmore,,149.875

r,25,19.524,108.482,291.159,16.95,6.214,, rmore,,158.46

r,26,20.084,118.046,295.597,16.95,6.407,, rmore,,167.857

r,27,20.536,128.202,298.9,16.95,6.607,, rmore,,177.137

r,28,21.148,139.669,303.86,16.95,6.816,, rmore,,187.495

r,29,21.626,151.737,307.186,16.95,7.032,, rmore,,197.507

r,30,22.07,1.734,310.057,16.95,7.257,, rmore,,207.681

r,31,23.541,182.225,323.433,16.95,7.490,, rmore,,223.527

r,32,25.066,201.932,336.467,16.95,7.730,, rmore,,241.046

r,33,26.871,223.700,351.448,16.95,7.979,, rmore,,257.765

r,34,28.707,247.928,365.663,16.95,8.235,, rmore,,275.922

r,35,30.533,274.36,378.934,16.95,8.500,, rmore,,294.461

r,46,14.58,35.599,238.702,16.95,4,, rmore,,238.702 !!横隔板截面 !!边支点

r,36,35.376,54.393,339.84,16.95,4.00,, rmore,,128.747 !!合拢段

r,37,35.376,53.405,339.84,16.95,4.00,, rmore,,127.015 !!中支点

r,38,71.981,454.871,543.871,16.95,8.5,, rmore,,684.565 !!主墩截面 !!截面1(实心)

r,39,39.932,82.876,212.138,8,5,, rmore,,203.709 !!截面(1/2 空心)

r,40,23.863,72.571,168.024,8,5,, rmore,,168.491 !!截面3 空心)

r,41,14.236,52.366,114.529,8,5,, rmore,,116.391 !!边墩截面 !!截面1(实心1)

r,42,36.113,66.217,177.703,7.7,4.7,, rmore,,165.119 !!截面2(空心1/2)

r,43,20.753,35.743,99.711,7,4,, rmore,,88.486 !!截面3(空心)

r,44,9.103,16.082,45.911,6.5,3.5,, rmore,,39.39 !!截面4(实心2)

r,45,22.673,23.005,79.361,6.5,3.5,, rmore,,61.527 !!截面5(1/2)空心2

r,49,16,21.958,68.708,6.5,3.5,, rmore,,55.682

!!建立第一跨梁节点 xl=atan(1.029/80) n,1,0,

n,2,2.25,2.25*xl n,3,6,6*xl n,4,9,9*xl n,5,11,11*xl n,6,13,13*xl-0.0045 n,7,15,15*xl-0.0135 n,8,17,17*xl-0.023 n,9,19,19*xl-0.038 n,10,21,21*xl-0.062 n,11,23,23*xl-0.092 n,12,25,25*xl-0.1275 n,13,27,27*xl-0.167 n,14,29,29*xl-0.2115 n,15,31,31*xl-0.261 n,16,33,33*xl-0.314 n,17,35,35*xl-0.3715 n,18,37,37*xl-0.435 n,19,39,39*xl-0.502 n,20,41,41*xl-0.5725 n,21,43,43*xl-0.8 n,22,45,45*xl-0.728 n,23,47,47*xl-0.8115 n,24,49,49*xl-0.95 n,25,51,51*xl-0.9925 n,26,53,53*xl-1.05 n,27,55,55*xl-1.1905 n,28,57,57*xl-1.325 n,29,59,59*xl-1.433 n,30,61,61*xl-1.5155 n,31,63,63*xl-1.6325 n,32,65,65*xl-1.7535 n,33,67,67*xl-1.879 n,34,69,69*xl-2.007 n,35,71,71*xl-2.1325 n,36,73,73*xl-2.2635 n,37,75,75*xl-2.399 n,38,77,77*xl-2.533 n,39,79,79*xl-2.631

xl2=atan(0.7/140)

local,11,0,80,80*xl-2.631,0,xl2

nsym,x,35,5,39,1!!复制粱结点从5到39结点编号增加35 csys,0

n,75,149.6,0.15852

local,12,0,150,0.15852,0,xl2 nsym,x,36,40,75,1 local,13,0,220,-1.6027,0 nsym,x,111,1,111,1 csys,0

n,500,80,80*xl-2.631 n,501,220,-1.6027 n,502,360,-1.6148

csys,0!!返回普通坐标 !!建立墩结点

!!PM112边墩截面

n,300,0,-2.165 !!实心)2米长 n,301,0,-4.365

n,302,0,-9.365!!(1/2空心)5M长 n,303,0,-13.365 n,304,0,-17.365 n,305,0,-21.365 n,306,0,-24.93 n,307,0,-28.93 n,308,0,-30.63

!!PM113主敦截面 n,309,80,0.5-6.665-0.8

n,310,80,0.5-3.5-6.665!!实心3.5m n,311,80,0.5-8.5-6.665!!1/2(空心)5m n,312,80,0.5-12.5-6.665!!(空心)4M n,313,80,0.5-16.5-6.665!! n,314,80,0.5-18.594-6.665 n,315,80,0.5-22.594-6.665

n,316,80,0.5-24.594-6.665+0.4!(实心)2m

!!pm114主墩截面 n,350,220,0.5-6.665-0.8!

n,351,220,0.5-3.5-6.665!!1/2空心5m n,352,220,0.5-8.5-6.665!!(空心)5m n,353,220,0.5-12.5-6.665!!空心)6.0M n,354,220,0.5-16.5-6.665!!1/2(空心)4M n,355,220,0.5-19.294-6.665!!(实心)2M n,356,220,0.5-23.294-6.665

n,357,220,0.5-25.294-6.665+0.4

local,13,0,220,-1.6027,0 nsym,x,17,300,316,1 csys,0 !!生成边跨1 mat,39 real,36 e,1,2 mat,38 real,46 e,2,3 mat,2 real,1 e,3,4 e,4,5 mat,2 real,2 e,5,6

*do,i,3,35,1 mat,i real,i e,3+i,4+i *enddo

*do,i,3,35,1 mat,i real,i e,38+i,39+i *enddo

*do,i,3,35,1 mat,i real,i e,74+i,75+i *enddo

*do,i,3,35,1 mat,i

real,i

e,114+i,115+i *enddo

*do,i,3,35,1 mat,i real,i

e,149+i,150+i *enddo

*do,i,3,35,1 mat,i real,i

e,185+i,186+i *enddo mat,41 real,38 e,39,500 mat,41 real,38 e,500,74 mat,2 real,2 e,41,40 mat,2 real,1 e,40,75 mat,40 real,37 e,75,111 mat,2 real,1 e,111,76 mat,2 real,2 e,76,77

mat,41 real,38 e,110,501 mat,41 real,38 e,501,221 mat,2 real,2 e,188,187 mat,2 real,1 e,187,222 mat,41 real,37 e,222,186 mat,2 real,1 e,186,151 mat,2 real,2 e,151,152 mat,41 real,38 e,185,502 mat,41 real,38 e,502,150 mat,39 real,36 e,112,113 mat,38 real,46 e,113,114

mat,2 real,1 e,114,115 e,115,116 mat,2 real,2 e,116,117

!!建立桥墩 !!边墩1 mat,1 real,42 e,300,301 mat,1 real,43 e,301,302 mat,1 real,44 e,302,303 mat,1 real,44 e,303,304 mat,1 real,44 e,304,305 mat,1 real,44 e,305,306 mat,1 real,49 e,306,307 mat,1 real,45 e,307,308 !!主墩1 mat,1 real,39 e,309,310 mat,1 real,40 e,310,311 mat,1

real,41 e,311,312 mat,1 real,41 e,312,313 mat,1 real,41 e,313,314 mat,1 real,40 e,314,315 mat,1 real,39 e,315,316 !!主墩2 mat,1 real,39 e,350,351 mat,1 real,40 e,351,352 mat,1 real,41 e,352,353 mat,1 real,41 e,353,354 mat,1 real,41 e,354,355 mat,1 real,40 e,355,356 mat,1 real,39 e,356,357 !!主墩3 mat,1 real,39 e,326,327 mat,1 real,40

e,327,328 mat,1 real,41 e,328,329 mat,1 real,41 e,329,330 mat,1 real,41 e,330,331 mat,1 real,40 e,331,332 mat,1 real,39 e,332,333 !!边墩2 mat,1 real,42 e,317,318 mat,1 real,43 e,318,319 mat,1 real,44 e,319,320 mat,1 real,44 e,320,321 mat,1 real,44 e,321,322 mat,1 real,44 e,322,323 mat,1 real,49 e,323,324 mat,1 real,45 e,324,325

!!添加质量惯距

!!两米段质量惯距

r,100,,,,1293348.4!!100米处和直线段的质量惯距 r,101,,,,1296019.6 r,102,,,,1297493 r,103,,,,1301181.8 r,104,,,,1311723 r,105,,,,1326622 r,106,,,,13467 r,107,,,,1370728 r,108,,,,1395061 r,109,,,,1422801 r,110,,,,1452757 r,111,,,,1483258 r,112,,,,1519733 r,113,,,,1558030 r,114,,,,1597229 r,115,,,,10292 r,116,,,,1687605 r,117,,,,1736625 r,118,,,,1790923 r,119,,,,1848602 r,120,,,,1911418 r,121,,,,1980562 r,122,,,,2052690 r,123,,,,2118097 r,124,,,,2192307 r,125,,,,22631 r,126,,,,2350704 r,127,,,,2432292 r,128,,,,2516392 r,129,,,,2679987 r,130,,,,2853514 r,131,,,,3048284 r,132,,,,3252032 r,133,,,,3462506

!!边跨直线3米段 r,148,,,,1293348.4*3/2

!!-700处3.75米段的质量惯距 r,134,,,,2725866.2

!!边支点横隔

r,135,,,,2246142.5

!!合龙段横隔 r,136,,,,833679/2

!!中支点横隔 r,137,,,,5292760

!!添加边墩质量惯距从上到下

!!变截面实心2米段 r,138,,,,1163498

!!1/2空心5米段 r,139,,,,1794766

!!空心段4米 r,140,,,,657211

!!1/2空心4米段 r,141,,,,1085079

!!实心2米段 r,142,,,,10567

!!添加主墩质量惯性距从上到下

!!实心3米段 r,143,,,,2345361

!!1/2空心5米段 r,144,,,,2869095

!!空心4米段 r,145,,,,1769087

!!1/2空心4米段 r,146,,,,2869095*4/5

!!实心2米段 r,147,,,,2345361*2/3

!!在结点处加入质量惯性距

!!边跨合龙段 !左边边跨 type,2 real,135 e,2 type,2 real,134 e,3 type,2 real,148 e,4 type,2 real,100 e,5

!!右边边跨 type,2 real,135 e,113 type,2 real,134 e,114 type,2 real,148 e,115 type,2 real,100 e,116

!!对第一个主跨的循环 *do,i,6,39,1 type,2 real,94+i e,i *enddo

!!对第二个主跨的循环 *do,i,41,74,1 type,2 real,59+i e,i *enddo

!!对第3个循环

*do,i,77,110,1 type,2 real,33+i e,i *enddo !!对第4个循环 *do,i,188,221,1 type,2 real,-88+i e,i *enddo

!!对第5个循环 *do,i,152,185,1 type,2 real,i-52 e,i *enddo

!!对第6个循环 *do,i,117,150,1 type,2 real,-17+i e,i *enddo

!!添加质量惯距到桥墩 !!左边墩 type,2 real,138 e,301 type,2 real,139 e,302 type,2 real,140 e,303 type,2 real,140 e,304 type,2 real,140 e,305

type,2 real,140 e,306 type,2 real,141 e,307 type,2 real,142 e,308 !!右边墩 type,2 real,138 e,318 type,2 real,139 e,319 type,2 real,140 e,320 type,2 real,140 e,321 type,2 real,140 e,322 type,2 real,140 e,323 type,2 real,141 e,324 type,2 real,142 e,325

!!左边主墩 type,2 real,143 e,310 type,2 real,144 e,311 type,2

real,145 e,312 type,2 real,145 e,313 type,2

r,400,,,,1769087*3/4 e,314 type,2 real,146 e,315 type,2 real,147 e,316

!!中间主墩 type,2 real,143 e,351 type,2 real,144 e,352 type,2 real,145 e,353 type,2 real,145 e,354 type,2

r,400,,,,1769087*3/4 e,355 type,2 real,146 e,356 type,2 real,147 e,357

!!右边主墩 type,2 real,143 e,327 type,2 real,144

e,328 type,2 real,145 e,329 type,2 real,145 e,330 type,2

r,400,,,,1769087*3/4 e,331 type,2 real,146 e,332 type,2 real,147 e,333

!!添加横隔板处质量惯性距 !!中支点 type,2 real,137 e,500 type,2 real,137 e,501 type,2 real,137 e,502

!!合龙段1 type,2 real,136 e,75 type,2 real,136 e,111 type,2

r,411,,,,1293348.4*3/10 e,40 type,2

r,412,,,,1293348.4*3/10 e,76

!!合龙段2 type,2 real,136 e,222 type,2 real,136 e,186 type,2

r,411,,,,1293348.4*3/10 e,187 type,2

r,412,,,,1293348.4*3/10 e,151

!!偶合结点 cp,1,uy,1,300 cp,2,uz,1,300 cp,3,rotx,1,300 cp,4,uy,309,500 cp,5,uz,309,500 cp,6,rotx,309,500 cp,7,ux,350,501 cp,8,uy,350,501 cp,9,uz,350,501 cp,10,rotx,350,501 cp,11,roty,350,501 cp,12,uy,326,502 cp,13,uz,326,502 cp,14,rotx,326,502 cp,15,uy,317,112 cp,16,uz,317,112 cp,17,rotx,317,112

加约束 d,308,all,all,0 d,316,all,all,0 d,357,all,all,0 d,333,all,all,0 d,325,all,all,0 allsel finish 求解 /solu

antype,2 !模态分析 acel,,9.8 !加载重力

modopt,subsp,20 !提取前30阶模态 mxpand,20,,,0 LUMPM,1 solve finish /post1 set,list