二次变截面梁的分析

shajiabaiyang

finish
/clear
/prep7
*dim,sect,char,16
*set,sect(1),'1','2','3','4','5','6','7','8','9','10','11','12','13','14','15','16'
*dim,y,array,16
*do,i,1,16,1
y(i)=8+0.0011556*((i-1)*5)**2-0.173333*((i-1)*5)
*enddo
et,1,188
et,2,82
mp,dens,1,2600
mp,ex,1,3.45e10
mp,prxy,0.1667
*do,i,1,16,1
!定义15个截面的形状
k,1,0,0,0,  
k,2,0,-0.3,0,
k,3,12,0,0,
k,4,12,-0.3,0,  
k,5,2.5,-0.3,0,
k,6,3,-0.3,0,   
k,7,9,-0.3,0,   
k,8,9.5,-0.3,0,
k,9,9.5,-y(i),0,
k,10,9,-y(i)+0.5,0,  
k,11,3,-y(i)+0.5,0,  
k,12,2.5,-y(i),0,
k,13,2.5,-y(i)+0.5,0
k,14,9.5,-y(i)+0.5,0
!连接成线
lstr,1,2  
lstr,1,3  
lstr,3,4  
lstr,4,8  
lstr,8,7  
lstr,7,6  
lstr,6,5  
lstr,5,2  
lstr,5,13  
lstr,12,9  
lstr,14,8  
lstr,7,10  
lstr,10,11  
lstr,11,6  
lstr,10,14
lstr,11,13
lstr,9,14
lstr,12,13
!给各边划分网格
lsel,s,line,,1,8,1
al,all
allsel

lsel,s,line,,7,9,2
lsel,a,line,,14,16,2
al,all
allsel

lsel,s,line,,11,12,1
lsel,a,line,,5,15,10
al,all
allsel

lsel,s,line,,10,13,3
lsel,a,line,,15,16,1
lsel,a,line,,17,18,1
al,all
allsel

asel,s,,,1,4,1
aglue,all
allsel

lsel,s,line,,2,,
lesize,all,,,48,
allsel

lsel,s,line,,1,3,2
lesize,all,,,1,
allsel

lsel,s,length,,0.5
lesize,all,,,2,
allsel

lsel,s,line,,6,13,7
lesize,all,,,24,
allsel

lsel,s,length,,2.5
lesize,all,,,10,
allsel

lsel,s,line,,9,12,3
lsel,a,line,,11,14,3
lesize,all,,,24
allsel

lsel,s,line,,10,
lesize,all,,,28
allsel

lsel,s,line,,4,8,4
lesize,all,,,10
allsel

asel,all
aatt,,,2,,
allsel
!给截面画镜像网格
amap,1,1,2,3,4
amap,2,5,6,13,11
amap,3,7,8,10,14
amap,4,13,14,9,12
secwrite,sect(i),SECT,,2
sectype,sect(i), BEAM, mesh, , 0
ACLEAR,ALL
ADELE,ALL,
LDELE,ALL,
kdele,all,
*enddo
!定义变截面梁
*do,i,1,16
SECTYPE,i,BEAM,MESH,
SECOFFSET,ORIGIN,,,   
SECREAD,sect(i),SECT,,MESH
*enddo
*do,j,1,15
SECTYPE,j+20,TAPER, ,
SECDATA,j,(j-1)*5,0,0,
SECDATA,j+1,j*5,0,0,
*enddo

*do,i,1,16
n,i,(i-1)*5,0,0
*enddo

*do,i,1,15
secnum,20+i
e,i,i+1
*enddo
/eshape,1

/solu
d,1,all
f,16,fz,-2e6
ACEL,0,0,9.8,
solve