平面应力问题
平面应力问题的理论基础
平面应力是什么
教授，"平面应力"是将3维问题降维为2维问题吗？
是的。3维弹性体求解的自由度非常多。但对于薄板结构，仅受面内荷载的情况，板厚方向的应力分量可视为零。这就是平面应力假设。
什么样的结构适用这个假设？
板厚与面内尺寸相比足够薄的结构。具体包括：
- 薄板受面内荷载（拉伸、压缩、剪切）
- 支架和冲压件的分析
- 受水压的薄壁
反过来，平面应力不适用的情况是？
板厚方向有约束的情况。例如厚实的大坝墙壁或长隧道的截面。这种情况应使用平面应变（板厚方向应变为零）假设。平面应力和平面应变看似相似但实际不同，混淆会导致较大误差。
控制方程
请说明平面应力的控制方程。
2维力的平衡方程：
应变适配条件（应变相容性）：
本构关系（Hooke定律）如何处理？
平面应力下的Hooke定律（矩阵形式）：
$E/(1-\nu^2)$ 比通常的 $E$ 大，对吗？泊松效应导致横向变形受约束，实际刚度增加。
好的观察。但在平面应力中，横向（$z$ 方向）变形是自由的。$E/(1-\nu^2)$ 出现是因为2个应力分量之间的耦合，与3维的 $E$ 含义略有不同。
重要的是：平面应力中 $\sigma_{zz} = 0$，但 $\varepsilon_{zz} \neq 0$：
板的厚度方向会产生应变。板在面内被拉伸时会变薄——这就是泊松效应。
平面应力 vs. 平面应变
请更详细地解释平面应变的区别。
平面应变的有效刚度比平面应力更大……即相同荷载下变形更小。
是的。$z$ 方向的变形完全受约束，导致3维约束效应使刚度增加。当 $\nu = 0.3$ 时，两者约相差10%。看似很小，但在应力评估中会产生不可忽视的差异。
Airy应力函数
听说平面问题可以用"应力函数"方法求解。
Airy应力函数 $\phi(x,y)$ 可以将平衡方程和相容条件合并为一个方程：
解出这个双调和方程后，应力分量为：
应力自动满足平衡方程！很方便。
经典弹性论的问题（有孔无限板、楔形、半平面接触问题）都可以用Airy函数求解。Kirsch问题（无限板圆孔周围应力集中）也是Airy函数解。这对FEM验证极其有用。
总结
整理一下平面应力的基础理论。
要点：
- $\sigma_{zz} = 0$ 假设 — 应用于薄板的面内问题
- 与平面应变的区别很关键 — 假设错误会导致10%以上误差
- 本构关系含 $E/(1-\nu^2)$ — 2维耦合效应
- $\varepsilon_{zz} \neq 0$ — 板厚方向应变仍然存在
- Airy函数可解古典问题 — 对FEM验证有用
要正确地将3维问题"降维"到2维，物理理解是不可或缺的。
完全同意。2维单元可以大幅降低计算成本，但前提是仅在假设成立的场景下使用。有疑虑时，可以用3维求解再与2维结果对比。
平面应力假设的理论依据
平面应力假设（σz=τyz=τxz=0）适用于"板厚相对于面内尺寸足够小"的薄板。1909年Kirchhoff整理了板弯曲与延伸问题的分离框架，成为后来薄板膜问题分析的基础。在飞机机翼桁腹板、汽车车门面板等板厚1～3mm的零件上，平面应力2D模型与3D壳体分析的偏差通常在1%以内。
平面应力问题的数值计算方法
FEM平面应力分析
用FEM求解平面应力问题时，使用什么样的单元？
2维平面应力单元。每个节点有2个自由度（$u_x, u_y$），约为3维固体单元自由度的1/3。
典型单元类型
| 单元 | 节点数 | 形状函数 | 精度 | 用途 |
|---|---|---|---|---|
| 3节点三角形（CST） | 3 | 线性 | 低（应变常数） | 自动网格填充 |
| 6节点三角形（LST） | 6 | 二次 | 高 | 复杂形状自动网格 |
| 4节点四边形（Q4） | 4 | 双线性 | 中等 | 规则网格 |
| 8节点四边形（Q8） | 8 | 二次 | 高 | 精密分析标准 |
CST（恒定应变三角形）如名所示，单元内应变是常数……这样就无法表示应力梯度。
完全正确。CST在应力集中部位（如孔周围）精度很低。教学演示或粗略筛选才用CST。实际工作中二次单元（6节点三角形或8节点四边形）是标准。
积分方案的影响
完全积分与低减积分的区别。
Q4（4节点四边形）的例子：
什么是剪切锁定？
弯曲变形时Q4的完全积分会产生有害的剪切应变，导致单元变得"刚硬"。实际上应该可以弯曲，但被锁定了。这就是剪切锁定。
对策：
1. 用低减积分 — 最简单，但需要沙漏控制
2. 非相容模式单元（incompatible modes） — Q4内加内部自由度。Abaqus的 CPS4I
3. 用二次单元 — Q8或T6不会出现剪切锁定
沙漏模式是什么样的变形？
低减积分的积分点（1个）应力为零，但单元以沙漏形状变形——这是应该不存在的零能量变形。数值上会自动出现但没有物理意义。沙漏控制（人工刚度）会抑制它。
求解器单元命名
各求解器的平面应力单元名称。
| 单元类型 | Nastran | Abaqus | Ansys |
|---|---|---|---|
| 4节点四边形 | CQUAD4（PSHELL） | CPS4, CPS4R | PLANE182 |
| 8节点四边形 | CQUAD8 | CPS8, CPS8R | PLANE183 |
| 3节点三角形 | CTRIA3 | CPS3 | PLANE182（退化） |
| 6节点三角形 | CTRIA6 | CPS6 | PLANE183（退化） |
Abaqus的CPS4R中"R"是低减积分的意思吗？
是的。"R = Reduced integration"。Abaqus中CPS4是完全积分，CPS4R是低减积分。实务中推荐用CPS4R（低减积分+沙漏控制）或CPS4I（非相容模式）。
Nastran的CQUAD4内置非相容模式，对剪切锁定有很强的鲁棒性。不需要特殊调整就能直接用。
网格生成要点
平面应力分析的网格需要特别注意什么？
板厚设置遗漏！这是盲点。3维固体单元的厚度是几何的一部分，但2维单元要当属性输入。
对。Nastran用PSHELL的板厚；Abaqus用 *SOLID SECTION 的厚度；Ansys用REAL常数或截面定义。板厚为1.0的分析结果力的量纲会不对，看结果时会显而易见。
总结
平面应力的数值方法，整理一下。
要点：
- 二次单元（Q8, T6）是工程标准 — CST精度不够
- 注意剪切锁定 — Q4的完全积分会出现。低减积分、非相容模式或二次单元可回避
- 沙漏模式须控制 — 低减积分的副作用
- 板厚设置必须做 — 平面应力单元的必要输入
- 应力集中部加密网格 — 精度关键部位要充分细化
CST与LST单元的精度对比
最简单的平面应力单元是3节点恒定应变三角形（CST），由Turner等人在1956年提出。CST内应力为常数，应力集中部需要大量单元才能达到精度。6节点二次三角形（LST）可以表示单元内线性应力变化，同样自由度下的精度是CST的4～8倍。1960年代Bogner的对比试验证实了这一点。
平面应力问题的实际应用
平面应力分析的实际应用
平面应力分析在实际中如何使用？
当今3维分析已普遍，平面应力仍有用处。
二维分析有效的场景
3维就能做，为什么还要用2维？
计算速度。2维成本约为3维的1/100以下。设计参数100个方案，3维逐个跑不现实，2维分钟级就完成。
应力集中系数评估
应力集中系数（$K_t$）评估可以用平面应力分析吗？
最典型的应用。无限板圆孔单轴拉伸的Kirsch问题，理论解 $K_t = 3.0$。用FEM再现这个就是平面应力的基础练习。
实务注意事项：
- 孔周网格足够密 — 二次单元至少16个围绕孔周
- 有限板的板宽效应 — 有限板 $K_t$ 与无限板不同。$d/W > 0.1$ 时需要修正
- 多孔干扰 — 孔间距小时应力集中会增加。FEM直接评估
Peterson's Stress Concentration Factors（应力集中系数手册）有大量表格。与FEM的关系？
Peterson手册是整理各种形状的 $K_t$ 的经典著作。FEM可以评估Peterson未收录的非标准形状（奇异切口、复合形状等）的 $K_t$。Peterson值与FEM结果的比对是基本的验证工作。
边界条件设置
2维平面应力问题的边界条件设置有难度吗？
3个要点：
1. 刚体运动约束 — 2维需要最少3个自由度的约束（平移2+旋转1）
2. 对称条件利用 — 形状和荷载都对称就用1/2或1/4模型。计算成本减少，网格细化更充分
3. Saint-Venant原理 — 荷载作用面对荷载细节敏感。关注部位距荷载面要足够远
Saint-Venant原理具体怎么用？
例如拉伸试样中央部分应力，端部约束条件精度差一些也没关系，中央应力仍正确。但"充分远离"的距离约为板幅的1～2倍是一个经验指南。
结果验证方法
平面应力分析的结果如何验证？
1. 反力检查 — 反力总和是否与荷载平衡
2. 应力连续性 — 单元间应力光滑过渡否。大的不连续表示精度不足
3. 网格收敛性 — 2～3个网格密度，应力5%以内收敛
4. 理论解比对 — 单轴拉伸、Kirsch问题等已知解比对
5. 位移量级确认 — 手算（如 $\delta = PL/AE$）验证位移数量级
实务检查清单
平面应力分析的检查清单。
"平面应力与平面应变混淆"是关键吧。混淆会改变基本物理。
空客A380机身外板的应力分析
空客A380（2005年首飞）机身外板采用厚1.5～4mm的铝锂合金，初期设计用平面应力有限元分析。单个外板面用约8000个单元建模，计算客舱增压荷载（0.6气压差）下的环向应力。精细的3D壳体分析与2D平面应力的主应力差为1.2%，设计文献中有记录。
平面应力问题的软件比较
平面应力分析工具
平面应力分析各求解器都能做吧。有差别吗？
平面应力是全部结构求解器的基本功能。差别在单元品质和易用性。
Nastran
Nastran的CQUAD4/CQUAD8单元被誉为高品质。特别CQUAD4虽然只有4节点，但内部含非相容模式，弯曲精度很高。用PSHELL属性定义板厚和MID3（面内刚性）。
Abaqus
Abaqus的CPS（Continuum Plane Stress）族丰富：
- CPS3 / CPS6 — 三角形（线性/二次）
- CPS4 / CPS4R / CPS4I — 四边形（完全积分/低减/非相容）
- CPS8 / CPS8R — 二次四边形
实务推荐为CPS8R（二次四边形低减积分）或CPS4I（非相容模式）。
Ansys
Ansys的PLANE182（线性）/PLANE183（二次）。用KEYOPT切换平面应力/平面应变/轴对称。KEYOPT设置错了会改变物理，非常关键。
开源软件
免费的FEM软件也能做平面应力分析吗？
有几个：
- FreeCAD + CalculiX — 开源。单元与Abaqus的CPS兼容
- Code_Aster + Salome-Meca — EDF开发。2维单元配置完备
- LISA — Windows用免费FEM。教学用最适合
- PrePoMax — CalculiX的GUI前端。易操作
教学的话LISA或PrePoMax很方便。
FEM基础学习，2维平面应力是最好的题材。用免费软件孔洞板的应力集中重现，对FEM理解效果非常好。
选型指南
总结？
平面应力是FEM的基础中的基础，工具差别比不上工程师理解程度对结果的影响大。
主要求解器平面应力单元比较
Abaqus的CPS4/CPS4R/CPS4I、ANSYS的PLANE182（KEYOPT(3)=0）、LS-DYNA的ELFORM=-1、NX Nastran的CPLSTS（平面应力）各使用不同的积分方式和稳定化技术。曲折基准试验（MacNeal-Harder 1985年标准问题集）中，Abaqus CPS4I与ANSYS PLANE182 Incompatible Modes分数最高，两者标准解误差在0.5%以内。
平面应力问题的前沿研究
平面应力前沿话题
平面应力是古典问题，但最前沿的研究有吗？
2维平面应力本身是成熟理论，但其拓展和应用在前沿研究中。
广义平面应力
"广义"平面应力是什么？
传统平面应力假设板厚方向应力均匀为零。广义平面应力（generalized plane stress）考虑板厚方向应力变化。厚板或板边缘附近 $\sigma_{zz} \neq 0$，常规平面应力会不准确。
用Asymptotic展开，板厚方向应力分布展开后用2维问题求解。板论（Kirchhoff、Mindlin）的3维补正也用到这种方法。
XFEM（扩展有限元法）
XFEM在平面应力问题上经常用吗？
2维裂纹问题是XFEM最基本的应用。常规FEM需要在裂纹尖端细化网格，XFEM的裂纹贯穿网格也可以。
Abaqus有XFEM功能吧。
Abaqus的 *ENRICHMENT + *CONTOUR INTEGRAL 可以计算2维裂纹的J积分、应力强度因子。网格生成工作大幅减少。
周期动力学
完全不同的方法也存在吗？
周期动力学（Peridynamics）用积分方程而非微分方程描述连续体。用2维平面应力问题时，裂纹产生与演进可用网格自由的形式自然表现。
常规FEM在裂纹尖端有奇点，周期动力学用"键"（bond）传递力，键断裂表示裂纹。2维脆性破坏模拟中研究快速进展。
相场法（Phase-Field）破坏
还有相场法，2维破坏问题中引人注目。损伤用连续场（0=健全、1=完全破坏）表示，裂纹分支、合并能自然追踪。在FEM框架可实现，Abaqus用户子程序（UEL/UMAT）可用。
2维平面应力是新破坏分析方法的"试验台"
正是。新方法总先在2维平面应力验证，再推广到3维。2维计算轻快，方法合理性可速查。平面应力是"已成熟理论"同时也是"最新研究支撑"。
总结
平面应力前沿话题，总结。
2维平面应力是"古"的"基本"最前沿研究支。
平面应力与混合单元定式
橡胶或金属大塑性变形时体积压缩性成问题，常规位移基平面应力单元会锁定。1983年Simo和Taylor发表增强假设应变（EAS）法，Abaqus CPS4I实装。对比CPS4R，曲折精度35%提升。薄板冲孔成形模拟（毛坯厚1mm）中，与实验的板厚分布误差从10%改善到3%。
平面应力问题的故障处理
平面应力分析故障
平面应力分析常见的故障。
平面应力是基础分析，初学者容易栽跟头。
平面应力与平面应变混淆
这是最多的错误吗？
间违最多。结果影响：
确认方法：
- Abaqus: 单元名首字 CPS = Plane Stress, CPE = Plane Strain
- Ansys: KEYOPT(3) = 0: Plane Stress, 1: Axisymmetric, 2: Plane Strain
- Nastran: PSHELL的MID参考隐含决定平面应力/应变
橡胶这样的非压缩材料用平面应变单元会发散？
$\nu = 0.5$ 时体积无变化（非压缩）的约束出现。平面应变的完全积分单元约束过剩，体积锁定发生。低减积分或混合单元（Abaqus的 CPE4H）需用。
应力集中值与理论不符
Kirsch问题（有孔板单轴拉伸）的 $K_t = 3.0$ 出不来。
原因排查：
1. 网格太粗 — 孔周围最少16个二次单元。一次单元需32个以上
2. 板为有限尺寸 — Kirsch理论解是无限板。$d/W < 0.1$ 以下才一致
3. 应力读取位置 — 最大应力在孔赤道（$\theta = 90°$）。单元积分点外推值或节点平均值读法不同结果不同
4. 单元类型 — CST无法表示应力梯度
应力读取方法会影响结果。
影响很大。FEM应力在积分点最正确。节点平均值（各单元值平均）光滑化但应力集中峰值低估。非平均化（unaveraged）应力云图看，单元间跳跃否检查必要。
位移与理论不符
均匀拉伸板的变位 $\delta = PL/AE$ 合不上。
检查项：
- 板厚设置 — 默认1.0还是？断面积 $A = t \times W$ 如果 $A$ 错位移会错
- 单位系统 — $E$ 单位与荷载、尺寸单位统一否
- 泊松效应 — 单轴拉伸也横收缩（$\varepsilon_y = -\nu \varepsilon_x$），$y$ 方向约束条件影响结果
单元扭曲精度下降
自动网格生成了扭曲单元。容许到什么程度？
2维单元品质指标：
| 指标 | 理想值 | 容许范围 | 影响 |
|---|---|---|---|
| 宽长比 | 1.0 | < 5.0 | 精度下降 |
| 歪斜度 | 0° | < 45° | 精度下降 |
| 雅各比行列式比 | 1.0 | > 0.3 | 负时单元退化 |
| 翘曲 | 0 | < 0.1 | 精度下降（3D问题） |
雅各比行列式为负什么意思？
单元节点顺序反了（单元翻过来了）。计算可能跑通但结果完全错。自动网格后必须品质查。
总结
平面应力故障处理，总结。
平面应力簡単見基本的結果直結分野。
FEM全分析基本的上成立。2次元平面応力基本徹底、複雑3次元分析誤防最良的。
平面应力中的厚縮小的見落
平面応力効果εz=−ν(εx+εy)/(1−ν)生、板面外変形。FEM求解器内部自動计算、後処理厚変化忘接触評価過小量見込誤起。成形仿真2010年代社的分析部門板厚変化的後処理整備、设计部門的情報連携標準化。