图 A 为一方板,在方板中心位置焊接一圆柱杆件,圆柱杆件顶端固结,方板四角作用有
集中力。图 B 为梁与壳简化计算的几何模型,图 C 为有限元模型,图 D 为梁与壳连接点
局部的单元的节点编号,按图中坐标系,节点 2 自由度 ROTZ 与其余节点自由度之间的关
系为:
写成标准方程形式如下:
! EX6.31 3D 梁单元与板壳单元刚接
finish $ /clear $ /prep7
l1=1.4 $ t=0.02 $ l2=1 $ r=0.1 $ p=20000 $ q=300 !板宽、板厚、柱高、柱半径、荷载
et,1,shell63,,,2 $ et,2,beam189 ! 定义两类单元:壳和梁单元
mp,ex,1,2.1e11 $ mp,prxy,1,0.3 $ r,1,t ! 定义材料性质及实常数(壳厚度)
sectype,1,beam,csolid $ secdata,r ! 定义梁截面及数据:实心圆柱
blc5,,,l1,l1 $ wprota,,90 $ asbw,all ! 创建方板,并切分为 4 部分
wprota,,,90 $ asbw,all $ wpcsys,-1
k,50,,,l2+t/2 $ kp0=kp(0,0,0) $ l,50,kp0 ! 创建点 50、获取 0,0,0 处点号、连线
lsel,s,loc,z,0.1,l2 ! 选择(柱)线
latt,1,,2,,,,1 $ lesize,all,,,4 $ lmesh,all ! 赋予属性、定义划分个数、划分网格
aatt,1,1,1 $ esize,0.1 ! 赋予面属性、定义单元尺寸
mshape,0,2d $ mshkey,1 $ amesh,all ! 定义网格形状、划分类型、划分网格
dk,50,all $ fk,1,fy,p $ fk,4,fx,p ! 在几何模型上施加约束和荷载
fk,3,fy,-p $ fk,2,fx,-p $ fk,4,fz,-q $ allsel,all
ce,1,0,143,ux,1,23,ux,-1,2,rotz,ny(143)-ny(23) ! 建立约束方程 1
ce,2,0,92,uy,1,30,uy,-1,2,rotz,-(nx(92)-nx(30)) ! 建立约束方程 2
/solu $ solve $ /post1 $ pldisp,1 ! 求解并进入后处理等
上述示例中,梁单元与壳单元共用节点,如不共用节点或各自节点独立,也可采用约束
方程。如梁单元节点位于某个壳单元的某位置,需要编写除 ROTZ 外的其余自由度的约束
方程;如果梁单元节点与壳单元节点位置重合,但各自独立,则需要将 ROTZ 外的自由度
耦合。因此建议采用共用节点,只需编写关于 ROTZ 的约束方程,以减少工作量。
如果采用刚性区法,可在节点 2 附近创建一刚性区(自动生成约束方程),即将该小区
域的 SHELL 视为刚性,这样势必就增加了结构的刚性。当采用 MPC184 建立几个刚性梁
单元时,情况与之类似,也增加了结构的刚性。因此,这两种方法的结果不如编写约束方
程合理。