ABAQUS中梁板壳单元的弯曲问题beam_plate_shell.pdfVIP

ABAQUS 中梁板壳单元的弯曲问题 曲哲 2007-4-3 一、Euler-Bernoulli 梁与Timoshenko 梁 在 ABAQUS 的单元库中，所有三次插值的梁单元（如 B23 ，B33 等），均为 Euler-Bernoulli 梁，而所 有线性和二次插值的梁单元（如 B21 ，B22 ，B31 ，B32 等），均为 Timoshenko 梁。 （1）细长梁与深梁 B23 为 2 结点三次插值的 Euler-Bernoulli 梁。由于在形成单元刚度矩阵时等效载荷项的被积函数至少 是 3 次的，所以至少需要 2 个积分点才能达到完全的高斯积分。而在 ABAQUS 中，B23 有 3 个积分点， 这意味着被积函数可以达到 5 次。总之 B23 是完全积分的单元。而 B21 和 B22 分别为 2 结点线性插值和 3 结点二次插值的 Timoshenko 梁，并且默认的采用减缩积分来避免剪切锁死。B22 只有 2 个积分点，B21 只有 1 个积分点，它们都只能达到 1 次的插值精度。 表 1：集中力作用下悬臂梁的自由端挠度（mm ） 细长梁（l/h=10 ） 深梁（l/h=3 ） 单元个数 1 2 4 材力解 1 2 4 材力解 B23 （E-B 梁） 4.000 4.000 4.000 4.000 0.1080 0.1080 0.1080 0.1080 B21 （Timoshenko 梁） 3.734 3.955 4.010 4.000 0.1086 0.1145 0.1160 0.1080 B22 （Timoshenko 梁） 4.028 4.028 4.028 4.000 0.1165 0.1165 0.1165 0.1080 表 1 比较了上述三种梁单元在应用于细长梁和深梁受弯时的表现。问题描述如图 1 所示，为端部受集 中载荷的悬臂梁。E-B 梁 B23 完全忠实于材料力解的解答，不考虑剪切应变的影响，并且只用 1 个 B23 单 元就可以得到与材力解一致的结果。B21 和 B22 考虑了梁的剪切变形，其分析得到的挠度略大于材力解。 同时可以看出，B21 和 B22 用于细长梁时并没有发生剪切自锁。 图 1：悬臂梁的构型图与截面图 图2 ：网格划分（2 个单元） （a ）B23 单元分析位移与应力结果 （b ）B21 单元分析位移与应力结果 （c ）B22 单元分析位移与应力结果 图3 ：悬臂梁在集中荷载作用下的位移与应力结果 图3 显示了用上述三种单元分析得到的梁的变形结果和应力结果，模型均划分为 2 个单元（图 2 ），且 均为图 1 中的细长梁。需要注意的是，ABAQUS 在后处理中并不显示梁内部的插值得到的位移，而只是以 直线连接梁的各结点，所以即使 B23 是三次插值，其变形图看上去仍然像是线性插值。同时从图 3 (b)中可 以看到，由于B21 单元只有一个积分点，所以单元内应力为常数。 （2 ）混合梁单元——轴力-弯矩耦合 ABAQUS 中还有一种用于模拟轴向刚度远大于弯曲刚度的梁的混合单元（Hybrid element ）。表2 比较 了混合单元 B21H 和一般的 Timoshenko 梁单元 B21 在压弯组合作用下的屈服荷载。构件与上节中一样， 也是端部受集中载荷的悬臂梁，只是端部添加了恒定的轴力（图 5 ）。这里采用理想弹塑性材料。网格划分 如图 4 所示。由表 2 可见，两种单元均能反映中轴力对于梁的屈服弯矩的影响。对于各向同性的对称截面， 拉、压对于屈服弯矩的影响也是相同的，表 2 中的结果也反映了这一点。由此可以看出 ABAQUS 中的 Timoshe