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

第1页下一页