【Abaqus有限元模拟】小片试验
1 问题描述
本算例采用矩形壳单元开展小片试验。模型由四个矩形单元组成，共包含九个节点。通过在外边界八个节点上施加由给定位移场计算得到的位移与转角，求解中心节点的响应，并将有限元结果与位移场直接计算结果进行对比，以检验所选壳单元对该位移模式的再现能力。
给定的横向位移场为：
w = 1 + 3x + 4y - 2x² + 4y² - 5xy
模型节点坐标及编号如图3-1所示。矩形区域宽度为4 m、高度为3 m，中心节点5位于坐标（2，1.5）处。试验的核心是利用边界节点的已知位移和转角反求中心节点响应，并据此判断单元是否能够通过小片试验。
图1-1 试验模型图
2 模型建立及边界条件的设置
在Part模块中选择3D、Deformable、Shell、Planar建立矩形壳模型。模型高度设置为3 m，宽度设置为4 m。建模时在x=2 m和y=1.5 m处对几何区域进行分割，使模型形成2×2的规则单元布局，并保证中心交点与节点5的位置一致。
材料采用各向同性线弹性模型，弹性模量E=3000 GPa，泊松比μ=0.3；壳截面厚度设置为0.1 m。上述参数均沿用原模型设置，本算例主要考察单元插值与位移场再现能力，材料参数不作为对比变量。
图1-2 所建模型
在Mesh模块中将模型划分为4个矩形单元。规则网格能够保证各单元几何形状一致，避免由于畸变或不规则网格对试验结果造成额外影响。网格划分后的模型如图3-3所示。
图1-3 划分网格后的模型
将各边界节点坐标代入给定的位移函数，求得四条边上8个节点沿z方向的位移。位移函数为：
w = 1 + 3x + 4y - 2x² + 4y² - 5xy
根据该位移场的一阶导数，可得到节点绕x轴和y轴的转角表达式：
UR1 = 4 + 8y - 5x
UR2 = -(3 - 4x - 5y)
将边界八个节点的坐标分别代入上述表达式，得到对应的U3、UR1和UR2，并将这些数值作为位移/转角边界条件施加于模型。随后创建Static, General通用静力分析步，输出中心节点5的位移与转角历史，提交Job进行求解。该设置能够避免额外外载荷干扰，使中心节点响应完全由边界位移场决定。
图1-4 施加边界条件
3 单元类型的选择
（a）单元类型为S4R
首先选用S4R四节点减缩积分壳单元进行计算。S4R单元计算效率较高，适用于一般薄壳问题。由于本模型仅包含四个规则单元，计算结果能够直接反映该单元对给定位移场的再现情况。结构变形结果如图3-5所示。
图1-5 单元类型为S4R时的结构变形图
小片试验共包含九个节点，其中外边界八个节点的位移与转角作为已知条件输入，中心节点5作为检验点。由有限元结果可得，中心节点沿z方向位移为-1 m，UR1=6 rad，UR2=12.5 rad。
将中心节点坐标（2，1.5）直接代入位移场及转角表达式，可得到完全相同的理论结果。有限元计算曲线与理论值吻合，说明S4R单元能够准确恢复中心节点的位移与转角。
图1-6 单元类型为S4R时的中间点位移和转角
（b）单元类型为S4
随后将单元类型更换为S4四节点完全积分壳单元，其余模型尺寸、材料参数、网格布置、分析步和边界条件保持不变，以保证两种单元结果具有可比性。S4单元对应的结构变形如图3-7所示。
图1-7 单元类型为S4时的结构变形图
提取中心节点5的结果可得：沿z方向位移为-1 m，UR1=6 rad，UR2=12.5 rad。该结果同样与位移函数直接计算值一致，表明S4单元也能准确再现给定的位移模式。
图1-8 单元类型为S4时的中间点位移和转角
4 结果分析和心得
通过对有限元计算结果与给定位移场理论结果进行对比可知，S4R和S4单元计算得到的中心节点位移及转角均为U3=-1 m、UR1=6 rad、UR2=12.5 rad，与理论值一致。因此，两种单元均能够通过本算例的小片试验，可认为其对该位移模式具有正确的插值和传递能力。
本算例也说明，小片试验的关键不在于模型规模，而在于边界位移和转角的准确施加。节点坐标、转角符号或自由度方向一旦设置错误，中心节点结果就会出现明显偏差。因此，建模时应逐点核对边界数据，并确保壳单元局部方向与转角定义一致。
对于非协调单元，能够通过小片试验是收敛的充分条件。需要说明的是，小片试验主要用于检验单元的基本完备性和协调性，并不能替代复杂工程问题中的网格收敛性分析。实际应用中仍应结合网格尺寸、单元畸变程度和具体受力形式对计算精度进行进一步验证。
