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

简支梁实体与预应力钢筋分析实例 /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,200000et,1,link8!定义预应力钢筋单元类型为Link8单元et,2,solid95 !r,1,agjx !定义实常数r,2mp,ex,1,egjxmp,prxy,1,0.3mp,dens,1,7800mp,dens,2,2500mp,alpx,1,xzxsmp,ex,2,ehntmp,prxy,2,0.3blc4,1,2,30/VIEW,1,1,1,1/ANG,1/AUTO,1/REP,FASTklist,allkwpave,6wpoff,-0.3wprot,0,0,90vsbw,1wpoff,0,0,-0.4vsbw,2wpoff,0,0.4wprot,0,90vsbw,allwpstylnummrg,all,lownumcmp,allesize,0.3klist,allllist,alllsel,s,28,38,10latt,1,1,1lmesh,allallsel,allgplotvsel,s,allvatt,2,2,2mshape,0,3dmshkey,1vmesh,allfinish/soludl,3,alldl,16,alldl,23,alldl,2,uydl,15,uydl,22,uy/AUTO,1/REP,FASTbfl,28,temp,-yjl/(xzxs*egjx*agjx)bfl,38,temp,-yjl/(xzxs*egjx*agjx)allselacel,9.8solvefinish/post1pldispplnsol,s,1方法二：体线独立耦合法建立预应力钢束命令流如下，并附有命令说明：/prep7*SET,eg,2e5!钢束弹模*SET,ag,140!钢束面积*SET,eh,4e4!混凝土的弹模*SET,r0,9345!曲线钢束的半径*SET,yyl,200000!钢束的预应力大小et,1,link8!定于预应力钢束单元类型为link8et,2,solid95!定义混凝土单元类型为solid95r,1,ag,yyl/eg/ag*1.036258!定义单元实常数，对于link8单元为钢束面积、以及钢束的初始应变（除模拟预应力张拉还会用于应力刚度矩阵的计算）r,2!solid95的单元实常数仅为x轴的调整，并且该实常数仅仅在单元坐标系统选择为1的时候填写。一般不填。mp,ex,1,eg!预应力钢束单元性质，包括之前定义的弹模，泊松比mp,prxy,1,0.3mp,ex,2,eh!混凝土单元性质，包括弹模、泊松比mp,prxy,2,0.2blc4,100,200,3000!定义长方体的长宽高（混凝土块体）/view,1,1,1,1!调整观察的角度/ang,1vplot!显示选择的体ksel,all!选择所有的关键点klist*get,kp0,kp,0,num,max !将最大的kp编号存于kp0中，因为事先不知道最大的关键点号是多少，所以先要存在kp0中lsel,none!不选择所有的线llistk,kp0+1,50,160!定义关键点，kp0是已有的最大关键点编号。k,kp0+2,50,160,3000k,kp0+3,50,800,1500larc,kp0+1,kp0+2,kp0+3,r0!将定义的以上三点连成一条圆弧，半径为r0。kp0+3用于定义圆弧所在面及曲线边中点kdele,kp0+3!删除kp0+3这个关键点klistllist*get,line1,line,0,num,min!将最小的线编号存于line1中，因为line1是最先定义的预应力钢束lsel,s,loc,z,0!选择z、y坐标均为0的线，用于施加约束（该线为右侧截面的梁底线）lsel,r,loc,y,0dl,all,uy!对线施加y方向的约束lsel,s,loc,z,3000!选择左侧截面梁底线（即z=3000，而y=0的线）lsel,r,loc,y,0dl,all,all!对选择的线施加约束（UX、UY、UZ）allsel,all!选择所有对象lsel,s,line1!选择线，最小值line1latt,1,1,1!将单元性质与选择的、尚未网格化的线联系起来lesize,all,50!线单元划分的数目lmesh,all!对line1进行分网vsel,all!选择所有的体vlistvatt,2,2,2!给体赋予相应的属性(混凝土)lsel,s,loc,z,0!选择z=0，而y坐标从10到140范围内的所有的线lsel,r,loc,y,10,140lesize,all,8!分网数目8lsel,s,loc,z,0!再次选择z=0，而不选y坐标从10到140范围内的所有线lsel,u,loc,y,10,140lesize,all,4!分网数目为4lsel,s,loc,y,0!选择x=0、y=0坐标处的线（仅有一根线被选中）注意这根线是一根分网的扫掠线。lsel,r,loc,x,0lesize,all,50!分网数目为50vsweep,all!通过扫掠网格的方法来用单元填充体allsel,all!选择所有的对象lsel,s,line1!选择line1最小值（预应力钢绞线）nsll,s,1!选择与所选择的线对应的所有的节点，即钢绞线上面的节点，因为网格数为50，所以节点总数为51cm,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到max1cp,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!建模结束，进入求解/soluallsel,allacel,9.8!考虑自重solveFINISH/POST1!进入后处理pldisp,1etable,sigi,ls,1plls,sigi,sigi,1plnsol,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.906778mp,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.117644 mp,ex,12,3.45e10 mp,nuxy,12,0.2 mp,dens,12,3221.585664 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.445716mp,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.45e10mp,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.578901 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.8972 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.899,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.964,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.648 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.464,91.768,281.923,16.95,5.854, rmore,141.068 r,24,19.072,99.936,287.364,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,164.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)空心2r,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*xln,3,6,6*xln,4,9,9*xln,5,11,11*xln,6,13,13*xl-0.0045n,7,15,15*xl-0.0135n,8,17,17*xl-0.023n,9,19,19*xl-0.038n,10,21,21*xl-0.062n,11,23,23*xl-0.092n,12,25,25*xl-0.1275n,13,27,27*xl-0.167n,14,29,29*xl-0.2115n,15,31,31*xl-0.261n,16,33,33*xl-0.314n,17,35,35*xl-0.3715n,18,37,37*xl-0.435n,19,39,39*xl-0.502n,20,41,41*xl-0.5725n,21,43,43*xl-0.648n,22,45,45*xl-0.728n,23,47,47*xl-0.8115n,24,49,49*xl-0.8995n,25,51,51*xl-0.9925n,26,53,53*xl-1.0895n,27,55,55*xl-1.1905n,28,57,57*xl-1.325n,29,59,59*xl-1.433n,30,61,61*xl-1.5155n,31,63,63*xl-1.6325n,32,65,65*xl-1.7535n,33,67,67*xl-1.879n,34,69,69*xl-2.007n,35,71,71*xl-2.1325n,36,73,73*xl-2.2635n,37,75,75*xl-2.399n,38,77,77*xl-2.533n,39,79,79*xl-2.631xl2=atan(0.7/140)local,11,0,80,80*xl-2.631,0,xl2nsym,x,35,5,39,1!复制粱结点从5到39结点编号增加35csys,0n,75,149.6,0.15852local,12,0,150,0.15852,0,xl2nsym,x,36,40,75,1local,13,0,220,-1.6027,0nsym,x,111,1,111,1csys,0n,500,80,80*xl-2.631n,501,220,-1.6027n,502,360,-1.6148csys,0!返回普通坐标!建立墩结点!PM112边墩截面n,300,0,-2.165 !实心)2米长n,301,0,-4.365n,302,0,-9.365!(1/2空心)5M长n,303,0,-13.365n,304,0,-17.365n,305,0,-21.365n,306,0,-24.93n,307,0,-28.93n,308,0,-30.63!PM113主敦截面n,309,80,0.5-6.665-0.8n,310,80,0.5-3.5-6.665!实心3.5mn,311,80,0.5-8.5-6.665!1/2(空心)5mn,312,80,0.5-12.5-6.665!(空心)4Mn,313,80,0.5-16.5-6.665!n,314,80,0.5-18.594-6.665n,315,80,0.5-22.594-6.665n,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空心5mn,352,220,0.5-8.5-6.665!(空心)5mn,353,220,0.5-12.5-6.665!空心)6.064Mn,354,220,0.5-16.5-6.665!1/2(空心)4Mn,355,220,0.5-19.294-6.665!(实心)2Mn,356,220,0.5-23.294-6.665n,357,220,0.5-25.294-6.665+0.4local,13,0,220,-1.6027,0nsym,x,17,300,316,1csys,0!生成边跨1mat,39real,36e,1,2mat,38real,46e,2,3mat,2real,1e,3,4e,4,5mat,2real,2e,5,6*do,i,3,35,1mat,ireal,ie,3+i,4+i*enddo*do,i,3,35,1mat,ireal,ie,38+i,39+i*enddo*do,i,3,35,1mat,ireal,ie,74+i,75+i*enddo*do,i,3,35,1mat,ireal,ie,114+i,115+i*enddo*do,i,3,35,1mat,ireal,ie,149+i,150+i*enddo*do,i,3,35,1mat,ireal,ie,185+i,186+i*enddomat,41real,38e,39,500mat,41real,38e,500,74mat,2real,2e,41,40mat,2real,1e,40,75mat,40real,37e,75,111mat,2real,1e,111,76mat,2real,2e,76,77mat,41real,38e,110,501mat,41real,38e,501,221mat,2real,2e,188,187mat,2real,1e,187,222mat,41real,37e,222,186mat,2real,1e,186,151mat,2real,2e,151,152mat,41real,38e,185,502mat,41real,38e,502,150mat,39real,36e,112,113mat,38real,46e,113,114mat,2real,1e,114,115e,115,116mat,2real,2e,116,117!建立桥墩!边墩1mat,1real,42e,300,301mat,1real,43e,301,302mat,1real,44e,302,303mat,1real,44e,303,304mat,1real,44e,304,305mat,1real,44e,305,306mat,1real,49e,306,307mat,1real,45e,307,308!主墩1mat,1real,39e,309,310mat,1real,40e,310,311mat,1real,41e,311,312mat,1real,41e,312,313mat,1real,41e,313,314mat,1real,40e,314,315mat,1real,39e,315,316!主墩2mat,1real,39e,350,351mat,1real,40e,351,352mat,1real,41e,352,353mat,1real,41e,353,354mat,1real,41e,354,355mat,1real,40e,355,356mat,1real,39e,356,357!主墩3mat,1real,39e,326,327mat,1real,40e,327,328mat,1real,41e,328,329mat,1real,41e,329,330mat,1real,41e,330,331mat,1real,40e,331,332mat,1real,39e,332,333!边墩2mat,1real,42e,317,318mat,1real,43e,318,319mat,1real,44e,319,320mat,1real,44e,320,321mat,1real,44e,321,322mat,1real,44e,322,323mat,1real,49e,323,324mat,1real,45e,324,325!添加质量惯距!两米段质量惯距r,100,1293348.4!100米处和直线段的质量惯距r,101,1296019.6r,102,1297493r,103,1301181.8r,104,1311723r,105,1326622r,106,1348967r,107,1370728r,108,1395061r,109,1422801r,110,1452757r,111,1483258r,112,1519733r,113,1558030r,114,1597229r,115,1640292r,116,1687605r,117,1736625r,118,1790923r,119,1848602r,120,1911418r,121,1980562r,122,2052690r,123,2118097r,124,2192307r,125,2263641r,126,2350704r,127,2432292r,128,2516392r,129,2679987r,130,2853514r,131,3048284r,132,3252032r,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,1064567!添加主墩质量惯性距从上到下!实心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,2real,135e,2type,2real,134e,3type,2real,148e,4type,2real,100e,5!右边边跨type,2real,135e,113type,2real,134e,114type,2real,148e,115type,2real,100e,116!对第一个主跨的循环*do,i,6,39,1type,2real,94+ie,i*enddo!对第二个主跨的循环*do,i,41,74,1type,2real,59+ie,i*enddo!对第3个循环*do,i,77,110,1type,2real,33+ie,i*enddo!对第4个循环*do,i,188,221,1type,2real,-88+ie,i*enddo!对第5个循环*d