[转帖]一个变截面三跨连续箱梁问题

MVH

　　/prep7

　　/title,three span continus grider

　　!上顶板关键点位置0m位置

　　k,1,0.0000,0.0000

　　k,2,-2.1000,0.0000 !定义关键点

　　k,3,-2.6000,-0.1250

　　k,4,-2.8000,-0.1250

　　k,5,-3.0000,-0.1250

　　k,6,-3.4857,-0.1036

　　k,7,-3.9714,-0.0821

　　k,8,-4.4571,-0.0607

　　k,9,-4.9429,-0.0393

　　k,10,-5.4286,-0.0179

　　k,11,-5.9143,0.0036

　　k,12,-6.4000,0.0250

　　!下底板关键点位置

　　k,13,-2.8000,-1.775

　　k,14,0.0000,-1.775

　　!上顶板关键点位置54m位置

　　kgen,11,1,12,1,0,0,54/10,100

　　!上顶板关键点位置55m位置

　　kgen,2,1,12,1,0,0,55,1100

　　!上顶板关键点位置55+39.5m位置

　　kgen,11,1101,1112,1,0,0,39.5/10,100

　　!上顶板关键点位置55+40m位置

　　!kgen,2,1101,1112,1,0,0,40,1100

　　!板的高度由1.775变化到3.875

　　!55m

　　c1=0.000720165 !54m的边跨二次抛物线系数

　　!39.5m

　　c2=0.0013459 !39.5m的中跨二次抛物线系数

　　!yb=1.775+c1*X1*X1 !边跨梁高变化曲线

　　!ym=1.775+c2*X2*X2 !中跨梁高变化曲线

　　*dim,x1,array,10 !定义边跨的10分点位置

　　x1(1)=54/10

　　x1(2)=2*54/10

　　x1(3)=3*54/10

　　x1(4)=4*54/10

　　x1(5)=5*54/10

　　x1(6)=6*54/10

　　x1(7)=7*54/10

　　x1(8)=8*54/10

　　x1(9)=9*54/10

　　x1(10)=10*54/10

　　*dim,x2,array,10 !定义中跨的10分点位置

　　x2(1)=39.5/10

　　x2(2)=2*39.5/10

　　x2(3)=3*39.5/10

　　x2(4)=4*39.5/10

　　x2(5)=5*39.5/10

　　x2(6)=6*39.5/10

　　x2(7)=7*39.5/10

　　x2(8)=8*39.5/10

　　x2(9)=9*39.5/10

　　x2(10)=10*39.5/10

　　!下底板边跨(54m)10分点关键点的y坐标位置

　　*dim,yb,array,10

　　yb(1)=c1*x1(1)*x1(1)

　　yb(2)=c1*x1(2)*x1(2)

　　yb(3)=c1*x1(3)*x1(3)

　　yb(4)=c1*x1(4)*x1(4)

　　yb(5)=c1*x1(5)*x1(5)

　　yb(6)=c1*x1(6)*x1(6)

　　yb(7)=c1*x1(7)*x1(7)

　　yb(8)=c1*x1(8)*x1(8)

　　yb(9)=c1*x1(9)*x1(9)

　　yb(10)=c1*x1(10)*x1(10)

　　!下底板中跨(39.5m)10分点关键点的y坐标位置

　　*dim,ym,array,10

　　ym(1)=-c2*(x2(1)-39.5)**2+2.1

　　ym(2)=-c2*(x2(2)-39.5)**2+2.1

　　ym(3)=-c2*(x2(3)-39.5)**2+2.1

　　ym(4)=-c2*(x2(4)-39.5)**2+2.1

　　ym(5)=-c2*(x2(5)-39.5)**2+2.1

　　ym(6)=-c2*(x2(6)-39.5)**2+2.1

　　ym(7)=-c2*(x2(7)-39.5)**2+2.1

　　ym(8)=-c2*(x2(8)-39.5)**2+2.1

　　ym(9)=-c2*(x2(9)-39.5)**2+2.1

　　ym(10)=-c2*(x2(10)-39.5)**2+2.1

　　!下底板边跨(55m)10分点关键点的坐标

　　kgen,2,13,14,1,0,-yb(1),x1(1),100

　　kgen,2,13,14,1,0,-yb(2),x1(2),200

　　kgen,2,13,14,1,0,-yb(3),x1(3),300

　　kgen,2,13,14,1,0,-yb(4),x1(4),400

　　kgen,2,13,14,1,0,-yb(5),x1(5),500

　　kgen,2,13,14,1,0,-yb(6),x1(6),600

　　kgen,2,13,14,1,0,-yb(7),x1(7),700

　　kgen,2,13,14,1,0,-yb(8),x1(8),800

　　kgen,2,13,14,1,0,-yb(9),x1(9),900

　　kgen,2,13,14,1,0,-yb(10),x1(10),1000

　　!下底板关键点位置55m位置

　　kgen,2,1013,1014,1,0,0,1,100

　　!下底板关键点位置55m+39.5m位置的关键点坐标

　　kgen,2,1113,1114,1,0,ym(1),x2(1),100

　　kgen,2,1113,1114,1,0,ym(2),x2(2),200

　　kgen,2,1113,1114,1,0,ym(3),x2(3),300

　　kgen,2,1113,1114,1,0,ym(4),x2(4),400

　　kgen,2,1113,1114,1,0,ym(5),x2(5),500

　　kgen,2,1113,1114,1,0,ym(6),x2(6),600

　　kgen,2,1113,1114,1,0,ym(7),x2(7),700

　　kgen,2,1113,1114,1,0,ym(8),x2(8),800

　　kgen,2,1113,1114,1,0,ym(9),x2(9),900

　　kgen,2,1113,1114,1,0,ym(10),x2(10),1000

　　!边墩处的横隔板位置

　　!1，2，3，……14为0m处横隔板位置的所有关键点

　　!51，52，，，，，64为-0.5m处横隔板位置的所有关键点

　　kgen,2,1,14,1,0,0,-0.5,50 !为-0.5m处横隔板位置

　　!中墩处的横隔板位置

　　!1001，1002，1003，……1014为54m处横隔板位置的所有关键点

　　!顶板关键点位置50m位置

　　!1101，1102，1103，……1114为55m处横隔板位置的所有关键点

　　!连成板面

　　!上顶板板厚0.25m

　　et,1,shell63 !定义箱梁的单元类型

　　mp,dens,1,2500 !定义质量密度

　　mp,ex,1,3.5e10 !定义弹性模量

　　mp,prxy,0.1667 !定义泊松比

　　r,1,0.25 !定义上顶板的厚度

　　*do,i,0,20,1 !通过循环生成顶板板面

　　a,1+i*100,1+(i+1)*100,2+(i+1)*100,2+i*100

　　*enddo !定义加速度，为y方向，大小9.8

　　a,1,51,52,2 !生成端横隔板的顶板板面

　　aatt,1,1,1 !赋予顶板以特性

　　asel,none !使面的选择集为空

　　!上顶板的加掖，板厚0.375m

　　r,2,0.375 !定义实常数(上顶板的加掖，板厚0.375m)

　　*do,i,0,20,1

　　a,2+i*100,2+(i+1)*100,3+(i+1)*100,3+i*100 !通过循环生成顶板加掖的板面

　　*enddo

　　a,2,52,53,3 !生成端横隔板顶板加掖的板面

　　aatt,1,2,1 !赋予顶板加掖以特性

　　asel,none !使面的选择集为空

　　!肋板与翼缘交接处，板厚0.5m

　　r,3,0.375 !定义实常数(肋板与翼缘交接处，板厚0.5m)

　　*do,i,0,20,1 !通过循环生成肋板与翼缘交接处的板面

　　a,3+i*100,3+(i+1)*100,4+(i+1)*100,4+i*100

　　a,4+i*100,4+(i+1)*100,5+(i+1)*100,5+i*100

　　*enddo

　　a,3,53,54,4 !生成端横隔板位置肋板与翼缘交接处的板面

　　a,4,54,55,5

　　aatt,1,3,1 !赋予顶板加掖以特性

　　asel,none !使面的选择集为空

　　!翼缘板，厚度渐变

　　!定义渐变的翼缘厚度

　　*dim,hd,array,9

　　hd(9)=0.2

　　hd(8)=(3.4/7)*0.3/3.4+0.2

　　hd(7)=2*(3.4/7)*0.3/3.4+0.2

　　hd(6)=3*(3.4/7)*0.3/3.4+0.2

　　hd(5)=4*(3.4/7)*0.3/3.4+0.2

　　hd(4)=5*(3.4/7)*0.3/3.4+0.2

　　hd(3)=6*(3.4/7)*0.3/3.4+0.2

　　hd(2)=7*(3.4/7)*0.3/3.4+0.2

　　hd(1)=0.5

　　*do,i,1,8,1

　　r,30+i,hd(i+1),hd(i+1),hd(i),hd(i) !定义实常数(渐变的翼缘厚度)

　　*enddo

　　*do,i,0,20,1 !通过循环生成翼缘板面

　　*do,k,5,11,1

　　a,k+i*100 ,k+(i+1)*100,k+1+(i+1)*100,k+1+i*100

　　a,k,k+50,k+1+50,k+1

　　aatt,1,k+26,1 !赋予翼缘以特性

　　asel,none !使面的选择集为空

　　*enddo

　　*enddo

　　asel,none !使面的选择集为空

　　!肋板厚度0.5

　　r,4,0.5 !定义实常数(肋板厚度)

　　*do,i,0,20,1 !通过循环生成肋板板面

　　a,4+i*100,4+(i+1)*100,13+(i+1)*100,13+i*100

　　*enddo

　　a,4,54,63,13 !生成端横隔板位置肋板的板面

　　aatt,1,4,1 !赋予肋板以特性

　　asel,none !使面的选择集为空

　　!下底板厚度0.4

　　*do,i,1,10,1 !通过循环生成边跨的底板板面

　　r,10+i,0.4 !定义底板的实常数

　　a,13+(i-1)*100,13+i*100,14+i*100,14+(i-1)*100

　　aatt,1,10+i,1 !赋予边跨底板以特性

　　asel,none !使面的选择集为空

　　*enddo

　　*do,i,1,13,1 !通过循环生成中跨的底板板面 ?13

　　r,20+i,0.4 !定义底板的实常数

　　a,13+(i+7)*100,13+(i+8)*100,14+(i+8)*100,14+(i+7)*100

　　a,13,63,64,14

　　aatt,1,20+i,1 !赋予中跨底板以特性

　　asel,none !使面的选择集为空

　　*enddo

　　r,50,2 !横隔板厚度2m

　　a,1,51,64,14 !生成端横隔板

　　a,1001,1101,1114,1014

　　aatt,1,50,1 !赋予横隔板以特性

　　asel,none !使面的选择集为空

　　allsel !全选

　　esize,0.4 !定义单元尺寸

　　mshape,0 !以四边形单元形状划分

　　mshkey,1 !以映射网格划分

　　amesh,all !划分所有面

　　nsel,all !全选所有节点

　　nsym,x,100000,1,200000,1 !镜像复制所有节点

　　esym,,100000,1,200000,1 !镜像复制所有单元

　　cskp,12,0,2114,2113,2101 !定义坐标平面

　　nsym,z,500000,1,1000000,1 !沿z方向镜像复制所有的节点

　　esym,,500000,1,1000000,1 !镜像复制所有的单元

　　allsel !全选

　　nummrg,all !进行编号合并

　　csys,0 !恢复整体坐标系

　　!施加支座约束

　　!支座位置选择在肋板处(端部支座)

　　nsel,s,loc,y,-1.775

　　nsel,r,loc,x,-3,-2.6

　　nsel,r,loc,z,-0.5,0

　　d,all,ux,,,,,uy,roty,rotz

　　allsel

　　nsel,s,loc,y,-1.775

　　nsel,r,loc,x,2.6,3

　　nsel,r,loc,z,-0.5,0

　　d,all,ux,,,,,uy,roty,rotz !(端部支座-0.5到0m位置)

　　allsel

　　nsel,s,loc,y,-1.775

　　nsel,r,loc,x,-3,-2.6

　　nsel,r,loc,z,189,189.5

　　d,all,ux,,,,,uy,roty,rotz

　　allsel

　　nsel,s,loc,y,-1.775

　　nsel,r,loc,x,2.6,3

　　nsel,r,loc,z,189,189.5

　　d,all,ux,,,,,uy,roty,rotz !(端部支座189到189.5m位置)

　　allsel

　　!中部支座

　　nsel,s,loc,y,-3.875

　　nsel,r,loc,x,2.6,3

　　nsel,r,loc,z,54,55

　　d,all,ux,,,,,uy,roty,rotz

　　allsel

　　nsel,s,loc,y,-3.875

　　nsel,r,loc,x,-3,-2.6

　　nsel,r,loc,z,54,55

　　d,all,ux,,,,,uy,roty,rotz !(中部支座54到55m位置)

　　allsel

　　nsel,s,loc,y,-3.875

　　nsel,r,loc,x,-3,-2.6

　　nsel,r,loc,z,134,135

　　d,all,ux,,,,,uy,roty,rotz

　　allsel

　　nsel,s,loc,y,-3.875

　　nsel,r,loc,x,2.6,3

　　nsel,r,loc,z,134,135

　　d,all,ux,,,,,uy,roty,rotz !(中部支座134到135m位置)

　　allsel !全选

　　!进行车道布载

　　!工况1一期横载

　　/solu !进入求解模块

　　allsel !全选

　　acel,,9.8 !定义重力加速度y方向

　　solve !求解

　　/post1 !进入后处理

　　plnsol,s,1 !显示主应力变化图

　　Nsel,s,loc,z,27 !选择10m处的所有节点

　　Esln,s,,all !选择包含这些节点的所有单元

　　Plnsol,s,1 !绘制主应力云图

　　Finish

　　Save

　　!挂-80荷载布置在边跨最不利位置

　　allsel

　　fdele,all,all

　　acel,,,

　　nsel,s,loc,x,-2.60000000000

　　nsel,r,loc,y,-0.125000000000

　　nsel,r,loc,z,21.9857142857

　　f,all,fy,-200000 !挂-80车轮位置1

　　allsel

　　nsel,s,loc,x,-3.48570000000

　　nsel,r,loc,y,-0.103600000000

　　nsel,r,loc,z,21.9857142857

　　f,all,fy,-200000 !挂-80车轮位置2

　　allsel

　　nsel,s,loc,x,-4.45710000000

　　nsel,r,loc,y,-0.607000000000E-01

　　nsel,r,loc,z,21.9857142857

　　f,all,fy,-200000 !挂-80车轮位置3

　　allsel

　　nsel,s,loc,x,-5.18575000000

　　nsel,r,loc,y, -0.286000000000E-01

　　nsel,r,loc,z,21.9857142857

　　f,all,fy,-200000 !挂-80车轮位置4

　　allsel

　　nsel,s,loc,x,-2.60000000000

　　nsel,r,loc,y,-0.125000000000

　　nsel,r,loc,z,23.1428571429

　　f,all,fy,-200000 !挂-80车轮位置5

　　allsel

　　nsel,s,loc,x,-3.48570000000

　　nsel,r,loc,y,-0.103600000000

　　nsel,r,loc,z,23.1428571429

　　f,all,fy,-200000 !挂-80车轮位置6

　　allsel

　　nsel,s,loc,x,-4.45710000000

　　nsel,r,loc,y,-0.607000000000E-01

　　nsel,r,loc,z,23.1428571429

　　f,all,fy,-200000 !挂-80车轮位置7

　　allsel

　　nsel,s,loc,x,-5.18575000000

　　nsel,r,loc,y,-0.286000000000E-01

　　nsel,r,loc,z,23.1428571429

　　f,all,fy,-200000 !挂-80车轮位置8

　　allsel

　　nsel,s,loc,x,-2.60000000000

　　nsel,r,loc,y,-0.125000000000

　　nsel,r,loc,z,27.0000000000

　　f,all,fy,-200000 !挂-80车轮位置9

　　allsel

　　nsel,s,loc,x,-3.48570000000

　　nsel,r,loc,y,-0.103600000000

　　nsel,r,loc,z,27.0000000000

　　f,all,fy,-200000 !挂-80车轮位置10

　　allsel

　　nsel,s,loc,x,-4.45710000000

　　nsel,r,loc,y,-0.607000000000E-01

　　nsel,r,loc,z,27.0000000000

　　f,all,fy,-200000 !挂-80车轮位置11

　　allsel

　　nsel,s,loc,x,-5.18575000000

　　nsel,r,loc,y,-0.286000000000E-01

　　nsel,r,loc,z,27.0000000000 !挂-80车轮位置12

　　f,all,fy,-200000

　　allsel

　　nsel,s,loc,x,-2.60000000000

　　nsel,r,loc,y,-0.125000000000

　　nsel,r,loc,z,28.5428571429

　　f,all,fy,-200000 !挂-80车轮位置13

　　allsel

　　nsel,s,loc,x,-3.48570000000

　　nsel,r,loc,y,-0.103600000000

　　nsel,r,loc,z,28.5428571429

　　f,all,fy,-200000 !挂-80车轮位置14

　　allsel

　　nsel,s,loc,x,-4.45710000000

　　nsel,r,loc,y,-0.607000000000E-01

　　nsel,r,loc,z,28.5428571429

　　f,all,fy,-200000 !挂-80车轮位置15

　　allsel

　　nsel,s,loc,x,-5.18575000000

　　nsel,r,loc,y,-0.286000000000E-01

　　nsel,r,loc,z,28.5428571429

　　f,all,fy,-200000 !挂-80车轮位置16

　　allsel

　　/solu

　　solve

　　/post1

　　plnsol,u,y

　　!挂-80布置在中跨最不利位置

　　/prep7

　　allsel

　　acel,,,

　　fdele,all,all

　　nsel,s,loc,x,-2.60000000000

　　nsel,r,loc,y,-0.125000000000

　　nsel,r,loc,z,89.3650000000

　　f,all,fy,-200000

　　allsel

　　nsel,s,loc,x,-3.48570000000

　　nsel,r,loc,y,-0.103600000000

　　nsel,r,loc,z,89.3650000000

　　f,all,fy,-200000

　　allsel

　　nsel,s,loc,x,-4.45710000000

　　nsel,r,loc,y,-0.607000000000E-01

　　nsel,r,loc,z,89.3650000000

　　f,all,fy,-200000

　　allsel

　　nsel,s,loc,x,-5.18575000000

　　nsel,r,loc,y,-0.286000000000E-01

　　nsel,r,loc,z,89.3650000000

　　f,all,fy,-200000

　　allsel

　　nsel,s,loc,x,-2.60000000000

　　nsel,r,loc,y,-0.125000000000

　　nsel,r,loc,z,90.5500000000

　　f,all,fy,-200000

　　allsel

　　nsel,s,loc,x,-3.48570000000

　　nsel,r,loc,y,-0.103600000000

　　nsel,r,loc,z,90.5500000000

　　f,all,fy,-200000

　　allsel

　　nsel,s,loc,x,-4.45710000000

　　nsel,r,loc,y,-0.607000000000E-01

　　nsel,r,loc,z,90.5500000000

　　f,all,fy,-200000

　　allsel

　　nsel,s,loc,x,-5.18575000000

　　nsel,r,loc,y,-0.286000000000E-01

　　nsel,r,loc,z,90.5500000000

　　f,all,fy,-200000

　　allsel

　　nsel,s,loc,x,-2.60000000000

　　nsel,r,loc,y,-0.125000000000

　　nsel,r,loc,z,94.5000000000

　　f,all,fy,-200000

　　allsel

　　nsel,s,loc,x,-3.48570000000

　　nsel,r,loc,y,-0.103600000000

　　nsel,r,loc,z,94.5000000000

　　f,all,fy,-200000

　　allsel

　　nsel,s,loc,x,-4.45710000000

　　nsel,r,loc,y,-0.607000000000E-01

　　nsel,r,loc,z,94.5000000000

　　f,all,fy,-200000

　　allsel

　　nsel,s,loc,x,-5.18575000000

　　nsel,r,loc,y,-0.286000000000E-01

　　nsel,r,loc,z,94.5000000000

　　f,all,fy,-200000

　　allsel

　　nsel,s,loc,x,-2.60000000000

　　nsel,r,loc,y,-0.125000000000

　　nsel,r,loc,z,96.0800000000

　　f,all,fy,-200000

　　allsel

　　nsel,s,loc,x,-3.48570000000

　　nsel,r,loc,y, -0.103600000000

　　nsel,r,loc,z,96.0800000000

　　f,all,fy,-200000

　　allsel

　　nsel,s,loc,x,-4.45710000000

　　nsel,r,loc,y,-0.607000000000E-01

　　nsel,r,loc,z,96.0800000000

　　f,all,fy,-200000

　　allsel

　　nsel,s,loc,x,-5.18575000000

　　nsel,r,loc,y,-0.286000000000E-01

　　nsel,r,loc,z,96.0800000000

　　f,all,fy,-200000

　　allsel

　　/solu

　　solve

　　/post1

　　phnsol,u,y

tonan_888@163.c