剪切锁定和体积锁定又叫剪切锁死、体积锁死，是有限元分析中常见的因单元选择造成的模型缺陷。

1、剪切锁定

关于剪切锁定的介绍可参考文章《有限元分析中的剪切锁定(shear lock)》，关于有限元分析中消除剪切锁定常用的方法可参考文章《消除有限元分析中剪切锁定(shear lock)的方法》，本文不再赘述，后面将详细讲解在ANSYS中处理剪切锁定的方法。

2、体积锁定

同剪切锁定相似，在模型出现体积锁定时，分析结果显示体积几乎不可压缩，体积应变表现为无穷小，体现为结构过硬，甚至导致非线性分析的不收敛。

在材料力学的诸多模量中，有一个与体积有关的模量，即“体积模量”，表达式如下：

体积模量可描述均质各向同性固体的弹性，可表示为单位面积的力，表示不可压缩性。当泊松比接近0.5时，上式中分母趋近于零，导致体积模量无穷大、体积应变无穷小。材料表现为不可压缩，在超弹性材料、塑性流动时出现这种不可压缩性的时候，会导致计算困难，产生单元伪应力。

3、ANSYS中剪切锁定和体积锁定的处理方法

对于单元过硬问题，一般采用不同的单元计算公式来消除。ANSYS 18X系列单元提供了B-Bar、URI、增强应变和混合U-P四种计算方法。下图所示为单元PLANE182所提供的四种消除剪切锁定和体积锁定的算法。

（1）选择缩减积分（B-Bar）

在PLANE182单元中选择缩减积分是在单元选项中设置关键字K1为Full Integration。此算法的思路是将容易发生体积锁定的单元积分项分解为体积项和偏差项两个部分，其中对体积项进行缩减积分，从而使单元软化。这种方法的优点是保持原模型分析方案不变，缺点是仍有发生体积锁定的可能。

（2）一致缩减积分（URI）

在PLANE182单元中将单元K1关键字设置为 reduced integration，用于处理体积锁定和剪切锁定问题。思路是对单元进行整体的缩减积分，使得单元软化。优点是计算性能好，速度快。缺点是可能会引起虚假的零能模式。

零能模式也叫沙漏模式，他在理论上是存在的一种变形模式，但是在实际模型中是不可能存在的。零能模式就是指有变形，但是不消耗能量。显然是一种伪变形模式，若不加以控制，计算模型会变得不稳定，并且计算出来的结果也是没有多大意义的。要抵制这种变形模式就得相应的消耗一定的能量，也就是沙漏能，如果这个比值太多，就说明模型和实际的变形有很大的差别，当然不正确。这也是缩减积分所付出的代价。

关于零能模式，张洪伟等编著《ANSYS非线性有限元分析方法及范例应用》一书37页的插图较好地做出了解释：

如上图所示，在单个积分点的低阶单元中，单元能得到的应变能信息只有横纵两个方向的长度，而对图中的变形“视而不见”。意味着单个积分点上所有应力分量都为零。

（3）增强应变积分

在PLANE182单元中设置K1关键字为 Enhanced strain，采用增强假定应变技术克服锁定问题。因为建立锁定是由于单元位移场函数不能模拟弯曲造成的，那么就给单元一个附加自由度——这个自由度专门用来解决剪切和体积锁定问题。这种方法在计算代价方面是一阶和二阶缩减积分的折衷效果。当其用于不可压缩材料时，需要与Mixed U-P方法一起应用，且只适用于四边形、六面体单元。

其简单形式，称为Simplified Enhanced Strain，主要用于解决剪切锁定问题，不适合体积锁定。

关于增强应变积分的更多原理，可参考这篇文献：ENHANCED STRAIN METHODS FOR ELASTICITY PROBLEMS。

（4）混合U-P

在ANSYS中用来处理体积锁定问题，适用于不可压缩材料(泊松比为0.5)和近似不可压缩材料(泊松比＞0.495)。将位移求解，从位移的函数导出改成将压力作为独立自由度求解，就可以处理大体积模量的问题。

4、注意事项

为了尽量减小剪切锁定和体积锁定对分析结果造成的影像，在ANSYS中尽量采用18X系列单元。其中plane182和solid185带有附加位移形状，plane183和solid186可使用一致缩减积分方案。模拟复杂结构、曲面结构，尽量采用高阶单元。对于实体结构和壳体结构，可采用低阶六面体和四节点四边形单元。

参考文献

【1】微信公众号：ANSYS学习与应用，关于ANSYS中的剪切锁定与体积锁定，2016-05-16