平面应变问题
平面应变问题的理论基础
平面应变是什么
在平面应力讲解中出现了"平面应变"，请详细说明。
平面应变（plane strain）是构造奥行方向（$z$ 方向）应变为零的假设：
不在奥行方向变形...这是什么情况？
奥行宽度相对于断面尺寸足够长的结构。在结构中央附近，端部影响消失，奥行方向变形受到约束。
典型例：
- 大坝的断面 — 奥行（河流方向）非常长
- 隧道断面 — 轴向截面一致
- 长堤坝 — 堤体纵向一致
- 压轧辊 — 宽度方向变形一致
- 地基滑面 — 假设奥行一致
平面应力是"薄板"，平面应变是"长柱状体的断面"对吧。
平面应变的本构关系
平面应变的胡克定律与平面应力有何不同？
平面应变的本构关系（矩阵形式）：
分母中出现了 $(1-2\nu)$！当 $\nu \to 0.5$ 时刚性变为无穷大。
这是平面应变最重要的特点。对于不可压缩材料（$\nu = 0.5$），体积变化为零，但同时还要求 $\varepsilon_{zz} = 0$，这严重限制了面内变形的自由度。这就是体积锁定的原因。
对于橡胶或近似不可压缩的材料，平面应变分析会很困难吧。
完全正确。当 $\nu > 0.49$ 时，普通单元失去作用。需要混合单元（将压力作为独立变量）或低减积分单元。
平面应变的应力
$\varepsilon_{zz} = 0$ 但 $\sigma_{zz} \neq 0$ 对吧。
是的。尽管 $z$ 方向应变为零，但由于泊松效应，$z$ 方向会产生应力：
如果 $\sigma_x + \sigma_y$ 是拉伸，那么 $\sigma_{zz}$ 也是拉伸...奥行方向会产生拉伸应力。
这是约束应力，在平面应变假设成立的条件下会自动产生。用3D分析求解同一问题时，在结构中央会看到这个 $\sigma_{zz}$，但在端部，$\sigma_{zz} \to 0$（趋向平面应力）。
地基力学中的平面应变
听说在地基工程中平面应变是标准的。
在挖掘、堤坝、挡土墙、隧道等纵向截面一致的地基问题中，平面应变事实上是标准做法。
但有一些要点要注意：
- 土的本构模型 — Mohr-Coulomb、Cam-Clay等在3D应力状态下定义，但在平面应变中，中间主应力 $\sigma_2 = \sigma_{zz} = \nu(\sigma_1 + \sigma_3)$ 会自动确定。这个 $\sigma_2$ 影响破坏判定
- 异向性 — 堆积土的水平和竖直刚性不同（横向各向同性）。即使在平面应变中也应该考虑这种异向性
Mohr-Coulomb破坏准则忽视中间主应力，所以在平面应变下是安全侧吗？
Mohr-Coulomb准则忽视中间主应力的影响，因此在平面应变下通常产生保守的（安全侧的）预测。使用考虑中间主应力的Drucker-Prager或Lade准则会得到更加现实的强度评估。
总结
整理平面应变的理论。
要点：
- $\varepsilon_{zz} = 0$ 的假设 — 适用于长柱状结构的断面分析
- $\sigma_{zz} = \nu(\sigma_x + \sigma_y)$ — 奥行方向产生约束应力
- $(1-2\nu)$ 在分母中 — $\nu \to 0.5$ 时刚性发散（体积锁定）
- 地基力学的标准假设 — 挖掘、隧道、堤坝
- 与平面应力混淆是禁忌 — 假设完全不同
平面应力和平面应变看起来相似，但物理意义根本不同。
是的。虽然都归结为2次元，但"什么是零"不同。应力为零（平面应力）还是应变为零（平面应变）。这个出发点的差异导致了一切的区别。
平面应变理论的历史背景
平面应变假设（εz=γyz=γxz=0）由Barré de Saint-Venant在1856年的剪应力分布理论中奠定了基础。它适用于"奥行相对于断面尺寸足够大"的结构，如隧道和长大坝，能将3D问题简化为2D。在地基工程中至今仍作为挡土墙、堤坝分析的设计标准。
平面应变问题的数值计算方法
FEM平面应变分析
平面应变的FEM单元与平面应力有什么区别？
单元形状和网格相同。区别在于本构矩阵 $[D]$。只需将平面应力的 $[D]$ 替换为平面应变的 $[D]$ 即可。
单元形状相同的话，就容易出现设置错误（选择平面应力/平面应变）而不被察觉。
这正是最大的陷阱。网格、荷载、边界条件完全相同，但单元类型的一个设置就会改变结果。
求解器的要素名称
| 单元 | Nastran | Abaqus | Ansys |
|---|---|---|---|
| 4节点四边形（平面应变） | CQUAD4 + PLPLANE | CPE4, CPE4R, CPE4H | PLANE182 (KEYOPT3=2) |
| 8节点四边形（平面应变） | CQUAD8 + PLPLANE | CPE8, CPE8R, CPE8RH | PLANE183 (KEYOPT3=2) |
Abaqus的"H"后缀是什么？
混合单元。将压力（静水压应力）作为独立变量处理，从而在 $\nu \to 0.5$ 时也能避免体积锁定。对于近不可压缩材料（橡胶、饱和土）的平面应变分析，CPE4H或CPE8RH是必须的。
体积锁定的问题
请详细说明体积锁定。
$\nu \to 0.5$ 的不可压缩材料中，体积变化 $\varepsilon_v = \varepsilon_x + \varepsilon_y + \varepsilon_z = 0$ 是约束条件。在平面应变下，$\varepsilon_z = 0$ 因此：
如果在每个积分点都施加这个约束，自由度会被过度约束。对于完全积分的4节点四边形，有4个积分点，每点1个约束，而只有8个自由度（4个节点×2DOF），自由度不足。
单元会"固化"...这就是锁定。
对策：
| 方法 | 原理 | 单元例（Abaqus） |
|---|---|---|
| 低减积分 | 减少积分点数以放松约束 | CPE4R |
| 混合单元 | 将压力作为独立变量分开处理体积约束 | CPE4H, CPE8RH |
| B-bar法 | 对体积应变进行平均化 | LS-DYNA的部分单元 |
| 二次单元 | 自由度较多，约束有余度 | CPE8R（完全积分CPE8仍有锁定） |
二次单元也会有锁定吗？
完全积分的二次单元（CPE8）容易发生锁定。应该使用低减积分（CPE8R）或混合单元（CPE8RH）。实务的铁则是"在平面应变中 $\nu > 0.45$ 就必须使用混合单元或低减积分"。
地基分析的设置
地基的平面应变分析有特殊设置吗？
有几个：
- 初始地压设置 — 地基在自重下具有初始应力。用 $K_0$ 法（$\sigma_h = K_0 \sigma_v$）给定初始应力
- 地下水位 — 孔隙水压影响有效应力。用耦合分析（Biot理论）还是非耦合处理
- 挖掘建模 — 用单元删除（element removal / kill）来表示挖掘
- 支护工程建模 — 用梁单元或壳单元表示喷射混凝土或岩石锚杆
挖掘时删除单元，应力平衡怎么处理？
删除的单元所承担的应力作为"释放力"加在剩余地基上。Abaqus中用 *MODEL CHANGE 删除单元，应力释放自动进行。Plaxis等地基专用软件能通过GUI直观设置。
总结
平面应变的数值方法整理一下。
要点：
- 单元类型的设置最关键 — 平面应力和平面应变只差一个设置值
- $\nu > 0.45$ 时混合单元必须 — 体积锁定对策
- 低减积分是通用对策 — 但要注意沙漏模式
- 地基分析需要初始地压·水压·挖掘步骤 — 专用软件更高效
- 验证结果中的 $\sigma_{zz}$ — 应与 $\nu(\sigma_x + \sigma_y)$ 一致
平面应变本构和厚方向应力
平面应变单元中厚方向应力σz不为零，而是 σz=ν(σx+σy) 计算。忽视这一点，仅用2D分量计算主应力会导致最大主应力被低估。1990年代一个配管法兰疲劳评估中，现场团队犯过这个错误，导致预测裂纹初生寿命约长30%。此事被记录在ASME PTC中。
平面应变问题的实务应用
平面应变的实务应用
请详细说明平面应变在实务中的应用。
有三个大的应用领域。
地基工程
平面应变分析的最大用户是地基工程。
| 问题 | 平面应变应用 | 注意事项 |
|---|---|---|
| 开挖 | 挖掘截面的变形·稳定性 | 挖掘端部效应需要3D |
| 盾构隧道 | 断面的环变形 | 掘进面是3D |
| 斜坡稳定性 | 滑面的安全系数 | 斜坡端部是3D |
| 挡土墙 | 土压和变形 | 标准应用 |
| 大坝 | 堤体断面应力 | 河谷部的3D效应另行评估 |
"端部效应需3D"的备注很多。
是的。平面应变是"无限长"的假设，所以端部（挖掘的起止点、斜坡的端部）会有3D效应。实务中通常先用平面应变做基本设计，再在重要位置用3D验证。这是两阶段方法。
金属加工（塑性加工）
金属加工模拟也用平面应变吗？
要点：
- 更新拉格朗日法或ALE法 — 应对大变形的定式化
- 摩擦模型 — 辊与材料间的库仑摩擦
- 温度耦合 — 热加工需要温度相关的材料特性
断裂力学
断裂力学中也用平面应变吗？
可以。裂纹前沿应力场取决于试件厚度：
- 薄试件 — 裂纹前沿接近平面应力。塑性区大
- 厚试件 — 裂纹前沿接近平面应变。塑性区小，更脆性
平面应变破壊靭性 $K_{IC}$ 是材料固有值，也是最保守的值。厚试件（ASTM E399）测定。用FEM模拟 $K_{IC}$ 试验时常用平面应变单元。
$K_{IC}$ 是平面应变值的原因是平面应变最苛刻（脆性）吗？
完全正确。平面应变有高三轴约束，塑性变形被抑制，裂纹容易扩展。所以 $K_{IC}$ 比平面应力破壞靭性低，用作设计的安全值。
实务检查清单
请给平面应变分析的检查清单。
"单位奥行"的结果是什么？
平面应变分析的是"奥行1单位"的切片。结果的力和力矩是"奥行1 m的"。要换算到实际结构全宽（如大坝宽100 m），需乘以100。这个换算常被初学者遗漏。
隧道挖掘分析的实践
东京外环道（外环道）大深度隧道（直径16.1m）设计中，SFRC与地山用平面应变单元建模。2014年检讨时用Abaqus分析挖掘各阶段应力重新分布，分段最大弯矩比3D分析误差在5%以内。相比3D，2D单个截面分析时间减少到1/80。
平面应变问题的软件比较
平面应变分析工具
有专门针对平面应变的工具吗？
通用FEM都支持，但地基工程有专用软件方便得多。
地基专用软件
| 软件 | 开发商 | 特点 |
|---|---|---|
| Plaxis 2D | Bentley Systems | 地基分析世界标准。GUI操作直观设置挖掘步骤 |
| GeoStudio (SIGMA/W) | Seequent | 浸透-变形耦合。与斜坡稳定性统合 |
| FLAC 2D | Itasca | 有限差分法。大变形·动力学强。岩石力学标准 |
| RS2 (Phase2) | Rocscience | 专门针对隧道·岩石挖掘 |
Plaxis 2D是世界标准吧。
Plaxis内置丰富的地基本构模型（Hardening Soil、Soft Soil Creep等），地基工程师无需编程就能处理复杂问题。挖掘步骤、水位变化、支护施工顺序都可通过GUI直观操作。
不过Plaxis仅限地基，结构部材（钢框架等）的耦合分析功能有限。需要详细分析结构和地基的联合效应时，得用Abaqus或Ansys这样的通用FEM。
通用FEM
用通用FEM做地基平面应变分析有什么好处？
Abaqus的 *SOILS 步（耦合浸透-变形分析）在地基工程中也被使用。
塑性加工软件
金属加工平面应变分析用什么？
选型指南
总结一下？
地基有地基、加工有加工，各有专用软件。平面应变不是"非主流"，而是各领域标准应用，这就清楚了。
平面应变绝不是"3D太重"时代的遗留，而是特定物理问题的最恰当建模。对奥行方向均匀的问题用3D没有必要。
求解器别平面应变单元对比
Abaqus的CPE4/CPE4R、ANSYS的PLANE182（KEYOPT(3)=2）、Nastran的CPLSTS（平面应变模式）、LS-DYNA的ELFORM=13都是4节点平面应变，但积分方式不同。Abaqus CPE4R的缩减积分+沙漏控制在地基弹塑性分析中收敛性高，2021年欧洲地基软件调查中被选为使用最多的求解器。
平面应变问题的前沿研究
平面应变的前沿课题
平面应变的前沿研究有哪些？
地基力学和断裂力学中有活跃展开。
MPM（物质点法）
MPM是什么？
物质点法（Material Point Method）是拉格朗日描述（追踪物质）和欧拉描述（固定网格）的混合方法。在大变形地基问题（滑坡、泥流、地基液化流动）中迅速普及。
2D平面应变的MPM：
- 后台网格（欧拉网格）求解方程
- 物质点（拉格朗日点）保持变形和历史
- 网格不变形 → 大变形也很稳定
解决FEM中网格变形的问题啊。
是的。土体滑动数百米的滑坡问题FEM不可能做到，但MPM可以追踪。Anura3D（代尔夫特工大）或CB-Geo MPM（剑桥大）可开源使用。
从平面应变到广义平面应变
"广义平面应变"是什么？
通常平面应变是 $\varepsilon_{zz} = 0$，而广义平面应变允许 $\varepsilon_{zz} = \text{const} \neq 0$。即奥行方向一致的应变（膨胀或收缩）可以有。
应用：
- 预应力混凝土的截面 — PC钢筋拉伸导致整个截面 $z$ 方向压缩
- 埋设管温度变化 — 管 $z$ 方向均匀膨胀
- 纤维复合材料微观结构 — RVE（代表体积单元）分析
Abaqus用 CPEG 单元（广义平面应变）直接建模。
孔隙水的连成分析
地基水耦合如何发展？
基于Biot压密理论的耦合浸透-变形分析正成为标准。
最新发展：
- 非饱和土耦合 — van Genuchten模型考虑含水量曲线
- 热-水-力学（THM）三耦合 — 放射性废弃物地层处置必须
- 化学-力学耦合 — 地基改良（药液注入等）建模
THM耦合是三种物理的耦合？
温度（Thermal）+水理（Hydraulic）+力学（Mechanical）三耦合。放射性废弃物产生的热推动地下水流动，地基应力状态改变，又影响透水性...。通常用2D平面应变截面分析。
总结
平面应变前沿研究总结。
平面应变是地基力学的基础，其上建立起最前沿的多物理耦合研究。
不完全平面应变与修正定式
长但有限的物体（如压轧辊）不能完全满足平面应变，端部有效应。1986年Seguin和Valliappan定式化了广义平面应变（Generalized Plane Strain），将端部挠度作为附加自由度纳入。Abaqus CGPE4单元实现此功能，用于压轧分析时端部应力预测精度比通常CPE4提高20～30%。
平面应变问题的故障排查
平面应变分析的故障
平面应变分析常见故障是什么？
除平面应力的故障外，平面应变有特有问题。
体积锁定
$\nu$ 较大的材料中位移异常偏小。
这是体积锁定，平面应变最代表性故障。
确认方法：
- 应力云图中有棋盘纹样（国际象棋棋盘纹）
- 改变 $\nu$，结果大幅变化（$\nu = 0.3$ 和 $\nu = 0.499$ 对比）
对策：
1. 换成低减积分单元 — CPE4R（Abaqus）
2. 换成混合单元 — CPE4H / CPE8RH
3. 二次元的低减积分 — CPE8R（最稳定）
4. $\nu = 0.5$ 改为 $\nu = 0.499$ — 避免完全不可压缩（使体积模量有限）
$\nu = 0.499$ 和 $\nu = 0.5$ 有差别吗？
数值上差别很大。$\nu = 0.5$ 时 $1/(1-2\nu)$ 无穷大，求解器无法正常工作。$\nu = 0.4999$ 时实质上不可压，但求解器稳定。不过混合单元即使 $\nu = 0.5$ 也没问题。
初始地压设置错误
地基分析中挖掘前的初始状态异常。
初始地压设置错误是地基分析最多见的故障。
检查项：
- $K_0$ 值正确吗 — 正常压密粘土 $K_0 = 1-\sin\phi'$（Jaky式）
- 自重与初始应力协调吗 — 初始步大位移出现就不协调
- 水压设置 — 全应力还是有效应力，$K_0$ 用法不同
初始步位移不为零是问题吗？
是问题。初始地压设置正确的话，初始步位移应为零（已经平衡）。出现位移说明设定的初始应力与自重产生的应力不一致。
Abaqus可用 *GEOSTATIC 步自动验证初始平衡。Plaxis用 $K_0$ procedure 自动生成初始状态。用通用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时）。
关联话题
了
详细
错误