ansys算变截面连续梁的例子
ansys算变截面连续梁的例子/prep7
!上顶板关键点位置0m
k,1,
k,2,-2.1
k,3,-2.6,-0.125
k,4,-2.8,-0.125
k,5,-3,-0.125
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.4,0.025
!下底板关键点位置
k,13,-2.8,-1.85
k,14,,-1.85
!上顶板关键点位置49m
kgen,9,1,12,1,0,0,49/8,100
!上顶板关键点位置50m
kgen,2,1,12,1,,,50,900
!上顶板关键点位置50+34.5m
kgen,9,901,912,1,,,34.5/8,100
c1=0.000843399 !边跨二次抛物线系数
c2=0.001701323 !中跨二次抛物线系数
!yb＝1.85＋c1*x1^2 !边跨梁高变化曲线
!ym＝1.85＋c2*x2^2 !中跨梁高变化曲线
*dim,x1,array,8
x1(1)=49/8
x1(2)=2*49/8
x1(3)=3*49/8
x1(4)=4*49/8
x1(5)=5*49/8
x1(6)=6*49/8
x1(7)=7*49/8
x1(8)=8*49/8
*dim,x2,array,8
x2(1)=34.5/8
x2(2)=2*34.5/8
x2(3)=3*34.5/8
x2(4)=4*34.5/8
x2(5)=5*34.5/8
x2(6)=6*34.5/8
x2(7)=7*34.5/8
x2(8)=8*34.5/8
*dim,yb,array,8
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)
*dim,ym,array,8
ym(1)=-c2*(x2(1)-34.5)**2+2.2025 ym(2)=-c2*(x2(2)-34.5)**2+2.2025 ym(3)=-c2*(x2(3)-34.5)**2+2.2025 ym(4)=-c2*(x2(4)-34.5)**2+2.2025 ym(5)=-c2*(x2(5)-34.5)**2+2.2025 ym(6)=-c2*(x2(6)-34.5)**2+2.2025 ym(7)=-c2*(x2(7)-34.5)**2+2.2025 ym(8)=-c2*(x2(8)-34.5)**2+2.2025 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,813,814,1,0,0,1,100 kgen,2,913,914,1,0,ym(1),x2(1),100 kgen,2,913,914,1,0,ym(2),x2(2),200 kgen,2,913,914,1,0,ym(3),x2(3),300 kgen,2,913,914,1,0,ym(4),x2(4),400 kgen,2,913,914,1,0,ym(5),x2(5),500 kgen,2,913,914,1,0,ym(6),x2(6),600 kgen,2,913,914,1,0,ym(7),x2(7),700 kgen,2,913,914,1,0,ym(8),x2(8),800 kgen,2,1,14,1,,,-0.5,50
et,1,63
mp,dens,1,2500
mp,ex,1,3.5e10
mp,prxy,0.1667
r,1,0.25
*do,i,0,16,1
a,1+i*100,1+(i+1)*100,2+(i+1)*100,2+i*100
*enddo
a,1,51,52,2
aatt,1,1,1
asel,none
r,2,0.375
*do,i,0,16,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
r,3,0.5
*do,i,0,16,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,16,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
r,4,0.5
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
c3=6.2474e-5
c4=0.000126024
*dim,h1,array,8
h1(1)=0.25+c3*x1(1)*x1(1)
h1(2)=0.25+c3*x1(2)*x1(2)
h1(3)=0.25+c3*x1(3)*x1(3)
h1(4)=0.25+c3*x1(4)*x1(4)
h1(5)=0.25+c3*x1(5)*x1(5)
h1(6)=0.25+c3*x1(6)*x1(6)
h1(7)=0.25+c3*x1(7)*x1(7)
h1(8)=0.25+c3*x1(8)*x1(8)
*dim,h2,array,9
h2(1)=0.25+c4*x1(1)*x1(1)
h2(2)=0.25+c4*x1(2)*x1(2)
h2(3)=0.25+c4*x1(3)*x1(3)
h2(4)=0.25+c4*x1(4)*x1(4)
h2(5)=0.25+c4*x1(5)*x1(5)
h2(6)=0.25+c4*x1(6)*x1(6)
h2(7)=0.25+c3*x1(7)*x1(7)
h2(8)=0.25+c4*x1(8)*x1(8)
h2(9)=0.4
*do,i,1,8,1
r,10+i,h1(i)
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,9,1
r,20+i,h2(i)
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
a,1,51,64,14
a,801,901,914,814
aatt,1,50,1
alls
esize,0.4
mshape,0
mshkey,1
amesh,all
nsel,all
nsym,x,10000,1,80000,1 esym,,10000,1,10000,1 cskp,12,0,1714,1713,1701 nsym,z,100000,1,100000,1 esym,,100000,1,1000000,1 alls
nummrg,all
csys,0
nsel,s,loc,y,-1.85
nsel,r,loc,x,-3,-2.6
nsel,r,loc,z,-0.5,0
d,all,ux,,,,,uy,roty,rotz
alls
nsel,s,loc,y,-1.85
nsel,r,loc,x,2.6,3
nsel,r,loc,z,-0.5,0
d,all,ux,,,,,uy,roty,rotz
alls
nsel,s,loc,y,-1.85
nsel,r,loc,x,-3,-2.6
nsel,r,loc,z,169,169.5
d,all,ux,,,,,uy,roty,rotz
alls
nsel,s,loc,y,-1.85
nsel,r,loc,x,2.6,3
nsel,r,loc,z,169,169.5
d,all,ux,,,,,uy,roty,rotz
alls
nsel,s,loc,y,-3.875
nsel,r,loc,x,2.6,3
nsel,r,loc,z,49,50
d,all,ux,,,,,uy,roty,rotz
alls
nsel,s,loc,y,-3.875
nsel,r,loc,x,-3,-2.6
nsel,r,loc,z,49,50
d,all,ux,,,,,uy,roty,rotz
alls
nsel,s,loc,y,-3.875 nsel,r,loc,x,-3,-2.6 nsel,r,loc,z,119,120 d,all,ux,,,,,uy,roty,rotz alls
nsel,s,loc,y,-3.875 nsel,r,loc,x,2.6,3 nsel,r,loc,z,119,120 d,all,ux,,,,,uy,roty,rotz alls
!进行车道布载
!一恒
/solu
alls
acel,,9.8
solve
/post1
plnsol,s,1
nsel,s,loc,z,10 esln,s,,all
plnsol,s,1
finish