MVH 发表于 2005-7-5 10:39

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

本帖最后由 wdhd 于 2016-3-7 14:20 编辑

  /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 发表于 2023-11-7 16:15

{:3_53:}{:3_53:}{:3_53:}
页: [1]
查看完整版本: [转帖]一个变截面三跨连续箱梁问题