【原创】three span continus grider

cxy1225

three span continus grider

cxy1225

是一个变截面三跨连续箱梁问题。

yiby

帮忙贴出来，方便大家看！
  
/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

yiby

很好的一个东东，热烈支持一下。
  
顶！！！！！！！！！！！！！

scleilei

支持!!!

sunhang2004

支持你！

jacobi

看看先，谢谢。

wangji790512

真的好！