ANSYS摩擦接触理论知识详解！！

在有限元仿真分析中，我们常常会遇到计算结果与理论预期严重不符的情况，其中剪切自锁（Shear Locking）和体积自锁（Volumetric Locking）是两类典型的数值问题。这两种在仿真分析中的“锁死”现象会导致模型表现出异常的刚硬特性，使变形计算结果远小于实际值，常常造成非线性分析不收敛。本文将解析这两种锁死现象的产生机理、典型特征以及工程实践中常用的解决方案。—01—剪切自锁：弯曲变形“冻结”之手剪切自锁是有限元分析中最常见的数值问题之一，我们将理论上或实际上不应该产生剪切变形的场景中，由于单元原因在仿真中产生了剪切变形，或寄生了剪切能量，最终导致位移结果偏小的这种现象称为剪切自锁。这种现象的本质是单元位移场插值函数无法准确描述弯曲变形，导致本应不存在的“伪剪切应变”出现，从而吸收了大部分外力功，使结构表现出异常的过高刚度。从力学原理看，当材料受纯弯曲时，根据材料力学中的平截面假定，变形前后截面应保持平面且始终垂直于中性轴，理论上剪应力应为零。但在有限元离散模型中，完全积分的线性单元每条边仅由两个节点定义，其边的位移模式只支持线性变化，不能呈现单元边界弯曲形式，如图1所示。 (a) 单元正常弯曲变形 (b) 单元剪切自锁变形 图1 单元正常弯曲与剪切自锁当应用完全积分线性单元来模拟受纯弯曲的小块材料时，单元上表面拉伸，下表面压缩，这些都与纯受弯的小块材料一致，但是在每个积分点处，中间直线被迫“倾斜”，表现为图1所示的竖直虚线和水平虚线之间的夹角不再为90°，这说明了这点的剪应力不为零，即在积分点处产生了伴生剪切，这与纯受弯的情况不符。这种伪剪切变形消耗了变形能，而不产生所希望的弯曲变形，使得实际弯曲变形减小，即表现为单元“过硬”，计算得到的挠度可能比理论值低30%以上。 理论和实践表明，剪切自锁的产生通常与三个因素密切相关：① 单元阶次：线性单元比二次单元更容易锁死；② 受力状态：纯弯或弯主导工况；③ 网格形态：单元长宽比越大，锁死越明显。从以上剪切自锁的三个密切相关因素可以看出，剪力自锁主要影响受弯曲载荷的完全积分线性单元的行为。而二次单元(如下图)的边界可以弯曲，故它没有严重的剪力自锁问题。 图2 完全积分的线性单元与二次单元但是如果二次单元发生扭曲或它的弯曲应力有梯度(复杂受力状态)，那么也将有可能出现某种程度的剪切自锁，这两种情况在实际问题中是可能发生的。因此，对于完全积分的线性单元，我们只有当确信载荷只会在模型中产生很小的弯曲时，才可以采用完全积分的线性单元。如果对载荷产生的变形类型有所怀疑，可应采用不同类型的单元分析对比。对于完全积分的二次单元，在复杂应力下也有可能发生自锁，因此如果在模型中应用这类单元，也应仔细检查计算结果。那么在工程实践中，我们如何缓解剪切自锁呢？通常可通过以下四种方法缓解剪切自锁：① 采用缩减积分单元缩减积分的概念如图3所示，上图为完全积分，下图为缩减积分（注意在ANSYS软件中 ，Full Integration 实际是选择缩减积分，即将求解应变的矩阵分为体积项和偏差项，对体积项采用缩减积分，对偏差项采用完全积分）。此处所讲的缩减积分（Reduced Integration）是对体积项和偏差项均采用缩减积分。但该方法 会出现沙漏(变形的零能模式)。图3 完全积分单元与缩减积分单元值得注意的是，剪切自锁现象还受到单元长宽比的较大影响。如下表所示，有分析表明：当简支梁网格长宽比达到10:1（200mm×20mm单元）时，即使用缩减积分，跨中位移误差仍高达30%；而当采用近正方形单元（20mm×33.33mm）时，结果与理论解（1.042mm）基本吻合。② 采用增强应变单元（非协调单元）该方式适用于四边形或六面体低阶单元，该方法是在单元模型中增加附加自由度，以克服剪切自锁。单元质量越好，计算精度越高。计算原理如图4所示，在单元中附加了自由度相当于提高了形函数的阶数；同时允许网格中产生缝隙和重叠，故也称为非协调模式。图4 增强应变单元③ 采用高阶单元如前所述，剪切自锁主要是因为线性单元的边不能弯曲，因此采用二次或更高阶单元（如PLANE183），可大为缓解剪切自锁现象，但这会造成计算成本增加。④ 网格优化控制长宽比（建议<3）或在弯曲方向布置足够的单元层数（一般不小于3层）。—02—体积自锁：不可压缩材料的“刚性枷锁”与剪切自锁不同，体积自锁主要发生在处理近似不可压缩材料（如橡胶、生物软组织）的仿真中，当泊松比υ接近0.5时，完全积分单元会出现体积应变被过度约束的现象。从本构关系看，体积模量K的计算公式为：当υ→0.5时，分母趋近零，K→∞，材料表现出无限大的体积刚度，导致单元无法发生合理的体积变形。在超弹性材料的大变形分析中，这种锁死会使静水压应力呈现不真实的“棋盘式”振荡分布，甚至造成求解中止。体积自锁的产生需要两个条件：① 使用完全积分单元；② 材料泊松比接近0.5。值得注意的是，平面应力问题不会出现体积自锁，因为平面外应变可补偿体积变化。对于非平面应力问题，如何缓解体积自锁呢？通常可采取如下措施:①采用混合U-P单元(杂交单元)混合U-P（位移-压力）公式是解决体积自锁最有效的方法。其核心思想是将静水压应力(图5)作为独立自由度引入单元，与位移场分开求解。静水压力是一个点处应力张量的均匀部分，也可以理解为所有主应力的平均值，其数学表达为 ρ=−31(σ1+σ2+σ3)。其中 σ1、 σ2、 σ3是主应力。它是一种等向性应力状态，作用于材料的各个方向，主要引起体积变化。图5 静水压力下的单元在ANSYS中，采用U-P杂交单元可通过插入命令流实现。设置KEYOPT(6)=1启用该技术，同时建议配合体积相容公差（VTOL，默认1e-5）调整，如下: 例如对SOLID185单元，可插入command添加如下命令流：① 在橡胶材料下启用U-P杂交单元:KEYOPT,MATID,6,1② 在求解栏中插入一条 Command：/SOLUSOLC,,,,1E-3② 此外缓解体积自锁现象的其他措施还包括：B-Bar方法：对体积应变项单独降阶积分；增强应变单元技术：适合六面体/四边形单元；材料调整：υ设为0.495～0.499平衡精度与收敛；缩减积分：但需警惕沙漏效应。—03—结语:精准仿真的“解锁”之道在实际工程中，应对剪切自锁和体积自锁的“抗锁死”技术需权衡计算精度、稳定性和效率。需要特别注意的是，应对剪切自锁或体积自锁的这些技术往往会相互影响甚至冲突。例如同时采用缩减积分和增强应变可能导致过度软化。在应对这些问题时建议：① 先尝试网格细化，厚度方向应至少3层单元；② 对于橡胶密封件，首选混合U-P单元技术；③ 关注沙漏能比率，尽量控制在<5%，不应超过10%。剪切自锁和体积自锁作为有限元分析中的典型数值问题，其本质都是单元形函数与实际变形模式不匹配导致的过度约束现象。理解其背后的力学原理（如平截面假定、体积模量发散）和数值机理（积分方案、形函数阶次），才能合理选择单元技术。现代商业软件（如ANSYS、ABAQUS）虽然提供了多种高级单元选项，但作为工程师仍需掌握“问题诊断—方法匹配—结果验证”的完整闭环，特别是在处理非线性材料或复杂接触时，建议通过小规模试算验证单元技术的适用性。只有正确“解锁”这些数值陷阱，才能使仿真真正成为设计优化的可靠指南。来源：薛定谔的Cube