二次变截面梁的分析
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
页:
[1]