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