变截面梁实例二

变截面梁实例二

finish
/CLEAR,NOSTART
/UIS,MSGPOP,3
/dev,font,1,Courier*New,400,0,-14,0,0,,,
/dev,font,2,Courier*New,400,0,-14,0,0,,,
!!This Macro was made for VIP entrance canopy of DongChang Fortune Plaza Project located in shanghai china
!!!!!!!!!!Parameter to define the model
L_x=6.4
W_y=3.3
L_div=5
W_div=4
U0_x=1
U0_Y=1
U0_Z=0
Dis=0.05
H_z=1.5
ediv_same=3
ediv_taper=4
Thi=20!!!mm
Gcoe=1.2
Wk=2.154!!!KN

!!!!!!!!!!!!!Customrized element type and section and materials
/prep7
et,1,189
et,2,44
keyopt,2,7,11
keyopt,2,8,11
mp,ex,1,2.06e11
mp,prxy,1,0.3
mp,dens,1,7850

sectype,1,beam,Hrec,,0
secoffset,cent
secdata,0.06,0.12,0.006,0.006,0.006,0.006

sectype,2,beam,ctube,,0
secoffset,cent
secdata,0.05,0.06,

sectype,3,beam,ctube,Strut,0
secoffset,cent
secdata,0.025,0.03,

sectype,4,beam,Hrec,,0
secoffset,user,0.06,0.06
secdata,0.12,0.2,0.006,0.006,0.006,0.006

sectype,11,beam,Hrec,,0
secoffset,cent
secdata,0.06,0.3,0.006,0.006,0.006,0.006

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Automatic Process below!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
k,1,u0_x,U0_y,U0_z
kgen,2,1,,,L_x
kgen,2,1,,,,W_y
kgen,2,1, , , , ,H_z
kgen,2,kp(u0_x,U0_y,U0_z+H_z), , ,L_x
l,1,kp(U0_x+L_x,U0_y,U0_Z)
l,1,kp(U0_x,U0_y+W_y,U0_Z)
lgen,W_div+1,1,,,,W_y/W_div
lgen,L_div+1,2,,,L_x/L_div
lsel,s,loc,x,U0_x
lgen,2,all,,,dis
lsel,all
lsel,s,loc,x,U0_x+L_x
lgen,2,all,,,-dis
lsel,all
/UIS,MSGPOP,3
lovlap,all
/UIS,MSGPOP,2
lsel,s,loc,x,U0_x-0.001,U0_x+0.001
ldele,all
lsel,all
lsel,s,loc,x,U0_x+L_x-0.001,U0_x+L_x+0.001
ldele,all
lsel,s,loc,y,u0_y
lsel,a,loc,x,U0_x+0.5*L_x
ksll,s
kgen,2,all,,,,,-1
lsel,all
ksel,all
gplot
l,kp(U0_x,U0_y,U0_Z+H_z),kp(U0_x,U0_y+2/3*W_y,U0_Z)
l,kp(U0_x+L_x,U0_y,U0_Z+H_z),kp(U0_x+L_x,U0_y+2/3*W_y,U0_Z)
!!!!!!!!!!!!!!!!!!!!!!!!
*do,i,1,2,1
lsel,s,loc,y,U0_y+(i-1)*W_y
ksel,s,loc,z,U0_z-1
latt,1, ,1, ,kp(u0_x+0.5*L_x,U0_y+(i-1)*W_y,U0_z-1) , ,2
lesize,all,,,ediv_same
lmesh,all
*enddo
*do,i,1,W_div-1,1
lsel,s,loc,y,U0_y+i/W_div*W_y
ksel,s,loc,z,U0_z-1
latt,1, ,1, ,kp(u0_x+0.5*L_x,U0_y+i/W_div*W_y,U0_z-1)

第1页下一页