平面应变问题 — 故障排除指南
更丰富的内容可在 plane-strain.html 中查看。
平面应变分析的故障
平面应变分析中常见的故障有哪些？
除了平面应力的问题外，还有平面应变特有的问题。
体积锁定
在泊松比$\nu$较大的材料中，位移异常小。
体积锁定是平面应变最具代表性的故障。
确认方法：
- 应力等值线中是否出现了棋盘模式（国际象棋棋盘图样）
- 改变$\nu$时结果是否大幅变化（比较$\nu = 0.3$和$\nu = 0.499$）
对策：
1. 改用低减积分要素 — CPE4R（Abaqus）
2. 改用混合要素 — CPE4H / CPE8RH
3. 二阶要素的低减积分 — CPE8R（最稳定）
4. 使用$\nu = 0.499$而非$\nu = 0.5$ — 避免完全非压缩（保持体积弹性模量有限）
ν = 0.499和ν = 0.5有差异吗？
数值上有很大的差异。当ν = 0.5时，$1/(1-2\nu)$变为无穷大，求解器无法正常工作。如果设置ν = 0.4999左右，实际上是非压缩的，但求解器保持稳定。但是如果使用混合要素，ν = 0.5也没有问题。
初始地应力设定错误
地基分析中开挖前的初始状态不正常。
初始地应力的设定错误是地基分析中最常见的故障。
确认项：
- $K_0$的值是否正确 — 正常固结粘土采用$K_0 = 1-\sin\phi'$（Jaky公式）
- 自重与初始应力是否一致 — 初始步骤出现大位移说明不一致
- 孔隙水压力的设置 — 总应力还是有效应力，会改变$K_0$的适用对象
初始步骤的位移不为零是问题吗？
这是一个问题。如果初始地应力设置正确，初始步骤的位移应该为零（已处于平衡状态）。位移出现说明自重产生的应力与设定的初始应力不一致。
在Abaqus中，可以通过*GEOSTATIC步骤自动检验初始平衡。在Plaxis中，可以通过$K_0$过程自动生成初始状态。对于手动设定的通用FEM，这个检查尤为重要。
开挖分析中应力跳跃
开挖步骤中移除要素后，应力异常增大。
开挖面的应力释放一次性进行时，特别是在弹塑性地基中会出现超调。
对策：
- 分阶段开挖 — 每步开挖1 m等
- 应力释放率分阶段 — 使用Plaxis中的$\beta$参数分阶段释放应力
- 增加增量数 — 每个开挖步骤的增量更细致
结果单位注意事项
平面应变结果的单位中常见的错误有什么？
最常见的是力的单位错误。平面应变是基于深度方向1单位长度的分析，因此：
- 反力为「kN/m」（深度方向每1 m）
- 集中荷载应作为「kN/m」输入（等同于线荷载）
- 转换到实际结构的全宽度需乘以深度长度
原意是施加「1 kN的集中荷载」，结果实际是「1 kN/m的线荷载」……
正是这样。平面应变中的「集中荷载」实际上是深度方向均匀分布的线荷载。面荷载也一样，如果输入值单位为kN/m²，结果则为kN/m。必须在报告中明确说明这个转换。
总结
整理一下平面应变的故障处理。
体积锁定和初始地应力是两大故障啊。
前者是FEM的数值问题，后者是地基力学的物理问题。两者都理解才能做好平面应变地基分析。
平面应变边界条件设置错误
平面应变分析中，对所有节点固定z方向位移的错误在初学者中很常见。正确做法是在构成律层面强制εz=0，不施加位移边界条件。在ANSYS PLANE182中，如果没有正确设置单元类型的KEYOPT(3)=2（平面应变），会被作为平面应力计算。这个设置错误可能导致应力差异达到ν/(1-ν)≈43%（ν=0.3时）。
相关主题
价值
详细
错误