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

1、简支梁实体与预应力钢筋分析实例/COM,Structural/PREP7egjx=2e5!Eyagjx=140!单根钢绞线面积ehnt=4e4!Ehxzxs=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

2、,0.3!定义solid95单元的泊松系数blc4,100,200,3000!定义梁体/view,1,1,1,1!定义ISO查看/ang,1vplot!绘制梁体kwpave,6!工作平面移动到关键点6wpoff,-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,

3、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!将体划分单元的形状定位HEXmshkey,1!采用MAPPE捌分器vmesh,all!对体进行划分单元finish/soludl,3,all!对线line7施加名束（UX,UY,UZ）dl,16,all!对线line31施加名束（UX,UY,UZ）dl,23,all!

4、对线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/post1plnsol,s,z,0,1!绘制Z方向的应力etable,sigi,ls,1!定义钢筋单元数据表plls,sigi,sigi,1!绘制上述应力!f

5、inish!/exit,nosav模拟预应力损失的两种方法方法1:通过体分割的办法建立预应力钢束，然后通过定义降温荷载来模拟预应力效应。这种方法在张立明编的ALGORANSYS&桥梁工程中的应用方法与实例中有一个算例。命令流如下：/prep7/title,theanalysisofthesolidbeamandprestressedtendon/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,li

6、nk8!定义预应力钢筋单元类型为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,al

7、l,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.8solve

8、finish/post1pldispplnsol,s,1方法二：体线独立耦合法建立预应力钢束命令流如下，并附有命令说明：/prep7*SET,eg,2e5*SET,ag,140*SET,eh,4e4*SET,r0,9345*SET,yyl,200000et,1,link8et,2,solid95!钢束弹模!钢束面积!混凝土的弹模!曲线钢束的半径!钢束的预应力大小!定于预应力钢束单元类型为link8!定义混凝土单元类型为solid95r,1,ag,yyl/eg/ag*1.036258!定义单元实常数，对于link8单元为钢束面积、以及钢束的初始应变（除模拟预应力张拉还会用于应力刚度矩阵的计算）r

9、,2!solid95的单元实常数仅为x轴的调整，并且该实常数仅仅在单元坐标系统选择为1的时候填写。一般不填。mp,ex,1,egmp,prxy,1,0.3mp,ex,2,ehmp,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

10、中lsel,none!不选择所有的线llistk,kp0+1,50,160!定义关键点，kp0是已有的最大关键点编号k,kp0+2,50,160,3000!将定义的以上三点连成一条圆弧，半径为r0。k,kp0+3,50,800,1500Iarc,kp0+1,kp0+2,kp0+3,r0kp0+3用于定义圆弧所在面及曲线边中点kdele,kp0+3klistllist*get,line1,line,0,num,min先定义的预应力钢束lsel,s,loc,z,0加约束（该线为右侧截面的梁底线）lsel,r,loc,y,0dl,all,uylsel,s,loc,z,3000y=0的线）lsel,r

11、,loc,y,0dl,all,allUZallsel,alllsel,s,line1latt,1,1,1格化的线联系起来lesize,all,50lmesh,allvsel,allvlistvatt,2,2,2土）lsel,s,loc,z,0围内的所有的线lsel,r,loc,y,10,140lesize,all,8lsel,s,loc,z,0到140范围内的所有线lsel,u,loc,y,10,140lesize,all,4lsel,s,loc,y,0根线被选中）注意这根线是一根分网的扫掠线。!删除kp0+3这个关键点!将最小的线编号存于line1中，因为line1是最!选才Zz、y坐标均为

12、0的线，用于施!对线施加y方向的约束!选择左侧截面梁底线（即z=3000,而!对选择的线施加约束（UXUY!选择所有对象!选择线，最小值line1!将单元性质与选择的、尚未网!线单元划分的数目!对line1进行分网!选择所有的体!给体赋予相应的属性（混凝!选才?z=0,而y坐标从10到140范!分网数目8!再次选择z=0,而不选y坐标从10!分网数目为4!选才?x=0、y=0坐标处的线（仅有一lsel,r,loc,x,0lesize,all,50vsweep,all填充体!分网数目为50!通过扫掠网格的方法来用单元allsel,all!选择所有的对象lsel,s,line1!选才Iline1最

13、小值（预应力钢绞线）nsll,s,1!选择与所选择的线对应的所有的节点，即钢绞线上面的节点，因为网格数为50,所以节点总数为51cm,cmljnod,node!将钢绞线节点编组，其组名为cmljnod*get,max1,node,0,count!选择钢绞线中的节点单元数目，将其赋给变量maxi,此变量可以作为*do循环的上界来确定循环的次数*dim,ojd,max1为行的范围*dim,jd,max1为行的范围*get,nod1,node,0,num,min*do循环的上界*SET,ojd(1),nod1*do,i,2,max1*SET,ojd(i),ndnext(ojd(i-1)*enddoa

14、llsel,allnsel,allcmsel,u,cmljnod点*do,i,1,max1*SET,nod1,ojd(i)*SET,j,nnear(nod1)!定义数组序列。数组名为ojd,max1!定义数组序列。数组名为jd,max1!选择集中最小的节点号，赋给变量nod1,作为nod1的值赋给ojd(1)!do循环，从2到max1!将ndnext(ojd(i-1)赋给ojd(i)!结束循环!选择所有!选择所有节点!不选择cmljnod,即钢绞线对应的节!i从1循环开始到max1!将ojd(i)赋给nod1!将nnear(nod1)赋名nj,是指ojdi邻近的节点吗？nnear()是一个函数

15、，表示选择最近的节点。*SET,jd(i),j!将j赋给jd(i)*enddo!结束循环nsel,all*SET,ji,1*do,i,1,max1cp,ji,ux,ojd(i),jd(i)编号。ux为耦合节点的x方向的位移。cp,ji+1,uy,ojd(i),jd(i)cp,ji+2,uz,ojd(i),jd(i)*SET,ji,ji+3*enddoallsel,all*SET,ji,value为空格时标识删除)*SET,i,*SET,max1,*SET,nod1,*SET,ojd,*SET,jd,*SET,j,!选择所有节点!将1赋给ji!i从1到max1孑由度，ji只是随机设置的耦合!ji

16、+3赋给ji!结束循环!选择所有单元!将所有的参数删除(第三项!建模结束，进入求解!考虑自重*SET,ag,*SET,eg,*SET,eh,*SET,kp0,*SET,r0,*SET,yyl,*SET,line1,finish/soluallsel,allacel,9.8solveFINISH/POST1pldisp,1etable,sigi,ls,1plls,sigi,sigi,1plnsol,s,1连续梁桥ansys命令流建模求解!连续梁桥/prep7et,1,4!定义梁单元et,2,21!定义mass21单元!定义粱材料！!泊松比！!密度mp,ex,2,3.45e10!直线段梁材料和1M

17、段梁材料mp,nuxy,2,0.2mp,dens,2,3302.153125mp,ex,3,3.45e10mp,nuxy,3,0.2mp,dens,3,3301.658695mp,ex,4,3.45e10mp,nuxy,4,0.2mp,dens,4,3299.906778mp,ex,5,3.45e10mp,nuxy,5,0.2mp,dens,5,3298.327219mp,ex,6,3.45e10mp,nuxy,6,0.2mp,dens,6,3292.351605mp,ex,7,3.45e10mp,nuxy,7,0.2mp,dens,7,3284.137255mp,ex,8,3.45e10mp

18、,nuxy,8,0.2mp,dens,8,3271.802136mp,ex,9,3.45e10mp,nuxy,9,0.2mp,dens,9,3260.41903mp,ex,10,3.45e10mp,nuxy,10,0.2mp,dens,10,3248.193657mp,ex,11,3.45e10mp,nuxy,11,0.2mp,dens,11,3235.117644mp,ex,12,3.45e10mp,nuxy,12,0.2mp,dens,12,3221.585664mp,ex,13,3.45e10mp,nuxy,13,0.2mp,dens,13,3208.826871mp,ex,14,3.4

19、5e10mp,nuxy,14,0.2mp,dens,14,3194.279207mp,ex,15,3.45e10mp,nuxy,15,0.2mp,dens,15,3179.924673mp,ex,16,3.45e10mp,nuxy,16,0.2mp,dens,16,3166.445716mp,ex,17,3.45e10mp,nuxy,17,0.2mp,dens,17,3152.555731mp,ex,18,3.45e10mp,nuxy,18,0.2mp,dens,18,3138.312105mp,ex,19,3.45e10mp,nuxy,19,0.2mp,dens,19,3124.795334

20、mp,ex,20,3.45e10mp,nuxy,20,0.2mp,dens,20,3110.7135mp,ex,21,3.45e10mp,nuxy,21,0.2mp,dens,21,3097.080875mp,ex,22,3.45e10mp,nuxy,22,0.2mp,dens,22,3083.186268mp,ex,23,3.45e10mp,nuxy,23,0.2mp,dens,23,3068.968824mp,ex,24,3.45e10mp,nuxy,24,0.2mp,dens,24,3055.612436mp,ex,25,3.45e10mp,nuxy,25,0.2mp,dens,25,3

21、045.857147mp,ex,26,3.45e10mp,nuxy,26,0.2mp,dens,26,3035.174287mp,ex,27,3.45e10mp,nuxy,27,0.2mp,dens,27,3026.696551mp,ex,28,3.45e10mp,nuxy,28,0.2mp,dens,28,3015.795365mp,ex,29,3.45e10mp,nuxy,29,0.2mp,dens,29,3007.710181mp,ex,30,3.45e10mp,nuxy,30,0.2mp,dens,30,3000.513837mp,ex,31,3.45e10mp,nuxy,31,0.2

22、mp,dens,31,2978.611375mp,ex,32,3.45e10mp,nuxy,32,0.2mp,dens,32,2958.618861mp,ex,33,3.45e10mp,nuxy,33,0.2mp,dens,33,2937.888072mp,ex,34,3.45e10mp,nuxy,34,0.2mp,dens,34,2919.475751mp,ex,35,3.45e10mp,nuxy,35,0.2mp,dens,35,2903.359983!6700处mp,ex,36,3.45e10mp,nuxy,36,0.2mp,dens,36,3302.153125!可以不用mp,ex,3

23、7,3.45e10mp,nuxy,37,0.2mp,dens,37,3302.153125mp,ex,38,3.45e10!-700处梁mp,nuxy,38,0.2mp,dens,38,3180.578901mp,ex,39,3.45e10!边支点横隔板mp,nuxy,39,0.2mp,dens,39,2868.674818mp,ex,40,3.45e10!合拢段横隔板mp,nuxy,40,0.2mp,dens,40,2868.674818mp,ex,41,3.45e10!中支点横隔板mp,nuxy,41,0.2mp,dens,41,2757.470588mp,ex,1,3.25e10!定义墩

24、材料属性mp,nuxy,1,0.2mp,dens,1,2650!定义实常数编号，面积,IYY,IZZ,宽,高,RMORE,抗扭惯距!主梁截面r,1,11.851,28.52,215.151,16.95,4,!直线段rmore,56.8972r,2,11.862,28.57,215.458,16.95,4.001,!截面100处rmore,56.905r,3,11.871,28.91,215.622,16.95,4.009,rmore,57.345r,4,11.903,29.03,215.78,16.95,4.025,rmore,57.938r,5,11.932,29.626,215.88,16

25、.95,4.049,rmore,58.502r,6,12.043,30.402,217.093,16.95,4.081,rmore,59.768r,7,12.199,31.407,218.899,16.95,4.121,rmore,61.445r,8,12.441,32.716,221.806,16.95,4.169,rmore,63.652r,9,12.673,34.21,224.418,16.95,4.226,rmore,66.058r,10,12.932,35.94,227.279,16.95,4.290,rmore,68.801r,11,13,221,37,941,230,512,16

26、.95,4.362,rmore,71.948r,12,13.534,40.233,233.872,16.95,4.442,rmore,75.432r,13,13.843,42.747,237.133,16.95,4.530,rmore,79.173r,14,14.213,45.688,241.054,16.95,4.627,rmore,83.484r,15,14.598,48.964,245.004,16.95,4.731,rmore,88.155r,16,14.979,52.547,248.817,16.95,4.843,rmore,93.111r,17,15.393,56.562,252.

27、927,16.95,4.963,rmore,98.528r,18,15.842,61.069,257.347,16.95,5.092,rmore,104.47r,19,16.293,65.992,261.673,16.95,5.228,rmore,110.742r,20,16.791,71.484,266.426,16.95,5.372,rmore,117.648r,21,17.303,77.568,271.225,16.95,5.525,rmore,124.861r,22,17.858,84.282,276.363,16.95,5.685,rmore,132.631r,23,18.464,9

28、1.768,281.923,16.95,5.854,rmore,141.068r,24,19.072,99.936,287.364,16.95,6.03,rmore,149.875r,25,19.524,108.482,291.159,16.95,6.214,rmore,158.46r,26,20.084,118.046,295.597,16.95,6.407,rmore,167.857r,27,20.536,128.202,298.9,16.95,6.607,rmore,177.137r,28,21.148,139.669,303.86,16.95,6.816,rmore,187.495r,

29、29,21.626,151.737,307.186,16.95,7.032,rmore,197.507r,30,22.07,164.734,310.057,16.95,7.257,rmore,207.681r,31,23.541,182.225,323.433,16.95,7.490,rmore,223.527r,32,25.066,201.932,336.467,16.95,7.730,rmore,241.046r,33,26,871,223,700,351,448,16.95,7.979,rmore,257.765r,34,28.707,247.928,365.663,16.95,8.23

30、5,rmore,275.922r,35,30.533,274.36,378.934,16.95,8.500,rmore,294.461r,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

31、,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

32、,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,2

33、3,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-

34、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.63

35、1,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空心）5

36、M长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-2

37、4.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,0ns

38、ym,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*e

39、nddo*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

40、,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,306

41、mat,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,3

42、55mat,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,1