From: 啸百川的博客
http://sunjd2001.blog.163.com/blog/static/11129445201291695818657/
命令流:
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!梁与壳的偏置,使梁壳共节点!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
FINISH
/CLEAR
/PREP7
!板的尺寸
L=2000
B=600
T=20
!H型梁的尺寸
H=200
W1=200
W2=200
T0=8
T1=10
T2=10
!板面上的均布荷载
Q1=10
!单元及材料定义
ET,1,SHELL181
ET,2,BEAM188
MP,EX,1,2E5
MP,PRXY,1,0.3
!板的截面类型、截面数据及截面偏置
!注意:对板而言,并不是所有板都能设置截面数据,可以使用SECDATA(即可以进行截面偏置设置的板)的SHELL单元详SECDATA帮助文档
SECTYPE,1,SHELL
SECDATA,T,1,0.0,3
SECOFFSET,BOT !偏置位置在板顶
!梁的截面类型H型、相应的截面数据及截面偏置
SECTYPE,2,BEAM,I, H_Beam, 2
SECDATA,W1,W2,H,T1,T2,T0,0,0,0,0
SECOFFSET, USER, 0,H !偏置位置在梁顶(对于H型梁,截面原点在梁底,详SECDATA帮助文档)
!建立板
WPROTA,,90
RECTING,-B/2,B/2,0,L
areverse,all
!对板划分网格
WPROTA,,,90
ASBW,ALL
WPCSYS,-1
LSEL,S,TAN1,X
LESIZE,ALL,,,10
LSEL,S,TAN1,Z
LESIZE,ALL,,,5
AATT,1,,1,,1 !第一个1表示1#材料,第三个1表示定义的单元类型1,第四个1表示定义的截面类型1
AMESH,ALL
AREVERSE,ALL !令所选择的面法线反置
!对梁划分网格
K,100,,500,500 !设置梁的K点,用于定位H型梁腹板的方向
LSEL,S,LOC,X,0
LATT,1,,2,,100,,2 !第一个2表示定义的单元类型2,100表示梁的定位K点,第二个2表示定义的截面类型2
LMESH,ALL
!在板面上施加面荷载
LSEL,S,LOC,Z,0
DL,ALL,,ALL
LSEL,S,LOC,X,B1/2
SFL,ALL,PRES,Q2
SFA,ALL,1,PRES,Q1
ALLSEL,ALL
!求解
/SOLU
SOLVE
!后处理
/POST1
PLDISP,1
!显示实体(可以清楚的显示截面偏置)
/ESHAPE,1
GPLOT
结果如图所示: