马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
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
|