内容发布更新时间 : 2024/5/20 20:46:43星期一 下面是文章的全部内容请认真阅读。
梁壳体单元都包含了这 3 个平动自由度，并且具有相同的物理意义，因此杆单元与梁壳体单元的连接采用公共节点即可，无需建立约束方程。
实际工程结构如比较复杂的杆系结构，为简化计算用杆元模拟长细比很大的杆件，或者两端构造非常明确的铰销连接杆件，其余采用梁单元模拟。如桁架结构中的腹杆可用杆元模拟，而弦杆等采用梁单元模拟，从而形成“杆梁”单元连接。同样地，当杆件支承在刚度较大的实体上时，可用杆元模拟杆件，而用实体单元模拟几何实体时就形成“杆体”单元连接。
! EX6.25 桁架、梁和实体组合结构计算 finish $ /clear $ /prep7
l1=0.8 $ l2=2 $ l3=1 $ h1=3 $ h2=1.2 $ b=1.5 $ err=0.1 ! 定义参数 et,1,solid95 $ et,2,link8 $ et,3,beam4 ! 定义 3 种单元类型
mp,ex,1,3.3e10 $ mp,prxy,1,0.2 ! 定义材料 1 性质
mp,ex,2,2.1e11 $ mp,prxy,2,0.3 ! 定义材料 2 性质
r,1,0.25*acos(-1)*(0.16*0.16-0.144*0.144) ! 定义实常数 1
r,2,0.00908,2.14e-3,1.45e-2,0.3,0.2 $ rmore,,6.52e-5 ! 定义实常数 2 r,3,0.00908,1.45e-2,2.14e-3,0.2,0.3 $ rmore,,6.52e-5 ! 定义实常数 2 /view,1,1,1,1 $ /ang,1,-120,zs,1 ! 设置视图方向
blc4,,,l1,b,h1 $ wpoff,,,h1-h2 $ vsbw,all ! 创建长方体并切分为 2 个体
wpcsys,-1 ! 工作平面恢复，并创建关键点和线
k,20,l1+l2,,h1 $ k,21,l1+l2+l3,,h1 $ k,22,l1+l2+2*l3,,h1 $ k,23,l1+l2+3*l3,,h1 k,24,l1+l2,b,h1 $ k,27,l1+l2+3*l3,b,h1
k,28,l1+l2+3*l3,b/2,h1 $ k,29,l1+l2+3*l3,b/2
l,6,20 $ l,12,20 $ l,20,21 $ l,21,22 $ l,22,23 $ l,7,24 $ l,11,24
$
k,25,l1+l2+l3,b,h1
$
k,26,l1+l2+2*l3,b,h1
l,24,25 $ l,25,26 $ l,26,27 $ l,6,24 $ l,20,24 $ l,21,25 $ l,22,26
l,23,28 $ l,28,27 $ l,28,29 $ dk,29,all $ a,1,all ! 在关键点和面施加约束
ksel,s,,,20,23 $ fk,all,fz,-3e4 ! 对 4 个关键点各施加 30kN 载荷
ksel,s,,,24,27 $ fk,all,fz,-2e4 $ ksel,all ! 对 4 个关键点各施加 20kN 载荷
fk,23,fx,1e4 $ fk,27,my,1000 ! 对另外两个关键点施加载荷和扭矩
esize,0.2 $ mshape,0 $ mshkey,1 ! 定义单元网格尺寸、类型等
vatt,1,,1 $ vmesh,all ! 赋予几何体材料和单元属性，划分网格
lsel,s,loc,x,l1+err,l1+l2 ! 选择拟划分为桁元的线，并定义组件
lsel,a,loc,x,l1+l2+l3 $ lsel,a,loc,x,l1+l2+2*l3 $ cm,linkline,line
latt,2,1,2 ! 赋予线材料 2、实常数 1 和单元 2 属性
lesize,all,,,1 $ lmesh,all ! 定义每线划分一个单元，并划分之
lsel,s,loc,x,l1+l2+err,l1+l2+3*l3 ! 选择除几何体外的线 cmsel,u,linkline $ lsel,u,loc,z,0,h1-err ! 从中去掉桁元线和柱线 latt,2,2,3 $ lesize,all,,,4 $ lmesh,all ! 赋予线材料 2、实常数 2、单元 3 属性
lsel,s,loc,x,l1+l2+3*l3 ! 选择柱竖线，赋予属性，划分单元
lsel,r,loc,z,0,h1-err $ latt,2,3,3 $ lesize,all,,,4 $ lmesh,all $ allsel,all finish $ /solu $ solve $ /post1 $ pldisp,1 ! 求解并进入后处理