三维弹性体分析
三维弹性体理论基础
三维弹性论的基础
老师，有限元法的结构分析说到底就是在求解「三维弹性论」吧？
控制方程
请讲解三维弹性论的控制方程。
有三个基本方程。
1. 平衡方程（力的平衡）
2. 应变-位移关系（协调条件）
3. 本构关系（胡克定律）
6个应力分量（$\sigma_x, \sigma_y, \sigma_z, \tau_{xy}, \tau_{yz}, \tau_{xz}$）、6个应变分量、3个位移分量。共15个未知数对应15个方程，是吗？
理解完美。将这15个方程以位移为唯一未知数整合起来就是Navier方程（Lamé-Navier方程）：
其中 $\lambda, \mu$ 是Lamé常数。$\mu = G$（剪切弹性模量），$\lambda = \frac{E\nu}{(1+\nu)(1-2\nu)}$。
各向同性弹性体的D矩阵
有限元中使用的 $[D]$ 矩阵是什么样的？
各向同性弹性体的三维本构关系（Voigt表示法）：
这个6×6的矩阵就是一切的基础啊。平面应力和平面应变的 $[D]$ 是这个的约化形式吧。
是的。三维的 $[D]$ 通过平面应力的条件 $\sigma_z = 0$ 或平面应变的条件 $\varepsilon_z = 0$ 进行约化。壳体和梁也是根据各自的假设从三维 $[D]$ 导出的。
异向性材料
非各向同性材料的情况呢？
一般的异向性弹性体 $[D]$ 有21个独立常数（6×6对称矩阵）。但实用上以下特殊情况居多：
| 材料对称性 | 独立常数 | 例子 |
|---|---|---|
| 各向同性 | 2 ($E, \nu$) | 钢、铝 |
| 横向各向同性 | 5 | 单向强化CFRP层、沉积土 |
| 直交异向性 | 9 | 木材、织物CFRP |
| 一般异向性 | 21 | 晶体 |
CFRP（碳纤维复合材）是横向各向同性吗？
单向材（UD材）的一层是横向各向同性。纤维方向与其垂直的平面对称。多层积层后整体变成直交异向性或更复杂的对称性。
三维分析的必要场景
二维近似不适用、必须三维分析的情况有哪些？
「不用三维就解不了」的问题比想象中多啊。
计算机性能提高使三维分析成为常态。但三维分析不等于正确分析。网格质量、边界条件、材料模型如果有误，三维也会出垃圾结果。
总结
让我总结一下三维弹性体分析的理论。
要点：
- 三个基本方程 — 平衡、协调、本构关系。共15未知数15方程
- $[D]$ 矩阵是一切基础 — 二维单元是这个的约化版
- 各向同性时2个常数，异向性最多21个常数
- 三维分析的必要场景 — 3D应力集中、板厚应力、接触、复杂形状
- 三维≠精确 — 输入质量决定输出结果
有限元的所有单元理论都归结为三维弹性论。理解这一点，自然就知道该选哪种单元了。
正是如此。三维弹性论是「树干」，各单元类型是「枝叶」。理解树干后，选择枝叶就顺理成章了。
3D弹性体控制方程的历史
三维弹性体平衡方程在1820年代由Navier和Cauchy独立推导。九分量应力张量σij与六分量应变张量εij由通用Hooke定律（σ=Cε，C为四阶弹性张量）相连。C11〜C66最多有21个独立分量。对于各向同性材料，可以归结为两个Lamé参数λ和μ，通过Young模量E=μ(3λ+2μ)/(λ+μ)变换。
三维弹性体数值计算方法
三维实体单元
三维实体单元有哪些种类？
基本形状有三种：六面体（hex）、四面体（tet）、五面体（wedge/prism）。
| 单元 | 一阶 | 二阶 | 精度 | 网格生成 |
|---|---|---|---|---|
| 四面体 | TET4（4节点） | TET10（10节点） | TET4：低 / TET10：高 | 自动网格易 |
| 六面体 | HEX8（8节点） | HEX20（20节点） | HEX8：中 / HEX20：非常高 | 自动网格困难 |
| 五面体 | WEDGE6 | WEDGE15 | 中～高 | 用于hex和tet的过渡 |
TET4精度低吗？
TET4是常应变单元（CST的三维版）。单元内应变恒定，无法表示应变梯度。实务上不应使用TET4。用TET10（二阶四面体）精度就足够。
但自动网格生成的话TET4更容易啊。
确实如此。但用TET4达到精度需要HEX20的5～10倍单元数，计算成本反而更高。TET10自动网格已是实务标准。
HEX vs. TET 的讨论
六面体（HEX）和四面体（TET），该用哪个？
这是有限元中最热议的话题之一。
HEX的优点：
- 用更少单元数达到相同精度
- 完全积分或减次积分都稳健
- 接触问题的稳定性高
HEX的缺点：
- 复杂形状的自动网格困难（或不可能）
- 网格生成费时
- 坚持六面体易导致单元质量恶化
TET10的优点：
- 几乎任意形状可自动网格
- 从CAD直接网格生成容易
- 网格自适应（adaptive refinement）容易
TET10的缺点：
- 相同精度需HEX的2～5倍自由度
- 不可压材料易产生体积锁定
- 接触面稳定性略低
到底该怎么办？
实务的方法：
- 设计阶段的初步筛选 → TET10自动网格快速迭代
- 需要详细评估的部位 → HEX20精密网格
- 复杂形状＋高精度 → TET10高密度网格（网格生成时间计入总成本时，通常低于HEX）
求解器别的单元名称
| 单元 | Nastran | Abaqus | Ansys |
|---|---|---|---|
| TET4 | CTETRA（4节点） | C3D4 | SOLID185 |
| TET10 | CTETRA（10节点） | C3D10, C3D10M | SOLID187 |
| HEX8 | CHEXA（8节点） | C3D8, C3D8R, C3D8I | SOLID185 |
| HEX20 | CHEXA（20节点） | C3D20, C3D20R | SOLID186 |
Abaqus的C3D10M中的「M」是什么？
Modified（改进版）
网格质量管理
三维网格质量怎么管理？
自动网格全部满足这些指标很难吧。
所有单元达到理想品质基本不可能。关键是应力集中部或着眼部位的单元质量要保证。远处的质量稍差影响不大（Saint-Venant原理）。
总结
三维实体单元数值方法，总结一下。
要点：
- TET10是实务标准 — TET4禁用。自动网格优势大
- HEX20精度高但网格生成费时 — 用于详细评估部位
- 不可压材料用混合单元 — C3D8H, C3D10MH
- 网格质量优先保证着眼部位 — 整体完美不如局部精确
- 单元数目标 — 相同精度TET10需HEX的2～5倍自由度
三维实体单元的Gauss积分阶数
20节点六面体单元（HEX20）的完全积分采用3×3×3（27点）Gauss积分。2×2×2（8点）减次积分计算量快3.4倍，但可能产生最多12个沙漏模式。Barlow（1976）证明了Gauss积分点是应力的超收敛点，在此采样应力可达比单元内部更高阶精度。
三维弹性体的实务应用
三维分析的实务流程
三维实体分析的典型工作流程是什么？
1. CAD数据导入及形状简化 — 删除圆角、省略孔、利用对称
2. 网格生成 — TET10自动网格。着眼部位尺寸控制
3. 材料、边界条件、荷载设定
4. 分析执行
5. 结果评估 — 应力云图、位移、反力确认
6. 网格收敛性验证
7. 与理论解、简易计算比对
CAD形状的简化
CAD数据可以直接用于有限元吗？
直接用CAD形状会：
- 微小形状特征（倒角、刻印、小孔）破坏网格
- 不必要的精度导致自由度爆炸
- 网格生成失败
简化指南：
- 着眼部位外的小形状 → 删除
- 着眼部位的圆角 → 保留（影响应力集中）
- 螺栓孔 — 仅在着眼部位保留孔。远处省略
- 焊缝 — 疲劳评估保留；静强度可省
简化程度判断很难啊。
关键是「会不会影响着眼部位的应力」。形状特征离着眼部位板厚的2～3倍以上时，通常无影响（Saint-Venant原理）。
对称条件的活用
对称条件怎么用？
形状和荷载对称时，可用1/2、1/4、1/8模型。
| 对称类型 | 模型 | 自由度削减 |
|---|---|---|
| 单面对称 | 1/2模型 | 1/2 |
| 两面对称 | 1/4模型 | 1/4 |
| 三面对称 | 1/8模型 | 1/8 |
| 周期对称 | 1个扇区 | 1/N（N = 扇区数） |
对称面的边界条件：对称面垂直方向位移为零。对称面上节点法向位移拘束。
用对称条件时，反对称模式（座屈等）会漏掉吗？
静力分析没问题。座屈或固有振动分析时反对称模式会漏，需分别求对称、反对称或用全体模型。
结果评估方法
三维分析结果怎么评估？
应力平滑化
有限元的应力在积分点最精确。节点应力由多个单元的外推值平均化。注意：
- 平均化应力 — 云图光滑，但可能低估峰值
- 未平均应力 — 单元间不连续露出。不连续大=网格太粗
未平均应力的不连续越小，网格越充分，是吗？
对的。未平均应力不连续<5%的von Mises应力时，网格充分。这个应力误差指标也用于网格自适应的判据。
子模型
用全体模型的结果作为局部部位的边界条件再次分析的子模型是三维分析的必备技术。全体细致化代替局部精密化。
实务检查清单
三维分析的检查清单。
「TET4禁止」要加粗啊。
现在还有初学者用TET4「结果对不上」的烦恼。改成TET10就解决，实例太多了。
实体单元的发动机缸体分析
汽车发动机缸体的CAE分析用铸铁（弹性模量150GPa）和铝合金（E=70GPa）分区，采用10节点四面体（TET10）或20节点六面体（HEX20）实体单元。日产VQ型发动机热应力分析（2003技术报告）用HEX20约120万节点模型验证温度循环后塑性应变分布，提前评估衬套开裂风险。
三维弹性体软件比较
三维分析工具
三维实体分析哪个求解器强？
线性静力分析所有求解器精度都足够。差异出现在网格生成和大规模问题处理能力。
网格生成工具
三维网格生成很困难啊。
网格生成工具（前处理软件）选择左右整个工作效率。
| 工具 | 特征 |
|---|---|
| HyperMesh（Altair） | HEX网格生成强。手工控制灵活性最高 |
| ANSA（BETA CAE） | 汽车业标准。批处理强 |
| Abaqus/CAE | 与Abaqus一体。TET自动网格易用 |
| Ansys Meshing | 与Workbench一体。膨胀层生成优秀 |
| Gmsh | 开源。TET自动网格。Python API |
| Salome-Meca | 开源。Code_Aster前处理 |
HEX网格自动生成的工具有吗？
完全自动HEX网格生成仍是未解决的难题。研究中有自动HEX算法（paving、plastering等），但复杂形状仍需手工。HyperMesh用手工＋半自动组合成果最好。
求解器大规模性能
百万自由度以上的模型哪个求解器快？
Ansys的GPU加速有效吗？
线性静力分析的直接法求解器GPU加速报告2～5倍高速化。特别是稀疏直接求解器的LU分解很适合GPU。大规模模型效果更大。
选择指南
总结？
CAD融合看Ansys、非线性强Abaqus、大规模看Nastran，是这样的格局啊。
对的。不过三家年年改进彼此弱点，差异在缩小。关键仍是网格质量和边界条件的正确性。
三维实体单元主流实现比较
NastranCHEXA/CPENTA/CTETRA、AbaqusC3D8/C3D10/C3D20、AnsysSOLID185/186/187分别是代表性三维实体单元。LS-DYNAEL FORM=1（单点积分）在碰撞分析中高速但需要沙漏控制。2024年现在，Ansys Mechanical APDL求解器的HEX20单元内存效率比前代提升2倍以上，1亿自由度超大规模分析可在标准计算机集群上执行。
三维弹性体的前沿研究
三维分析的前沿课题
三维弹性体分析最前沿往什么方向发展？
三个大方向。
无网格法·粒子法
不用网格的方法有吗？
但现在的问题是：
- 精度有时不如FEM
- 边界条件处理复杂
- 计算成本高
所以实务普及有限。特殊问题（爆炸、冲击、大变形破坏）用。
等几何分析（IGA）
IGA经常听到啊。
等几何分析（Isogeometric Analysis）用CAD的形状表示（NURBS、T-样条）直接作分析基函数。网格生成这一步省了。革新的方法。
好处：
- CAD形状精确表示（形状近似误差零）
- 高阶连续性（$C^2$ 以上）让应力更光滑
- 网格生成不用（细化通过knot插入）
实用化进展如何？
LS-DYNA有部分IGA单元实现。Cimne的Kratos框架、Abaqus用户子程序也可用。但成为通用有限元主流还需时间。
AI与FEM的融合
AI用于有限元的研究有吗？
快速发展的领域：
- 代理模型 — 有限元结果用神经网络学习，实时预测应力
- Physics-Informed Neural Networks (PINN) — 支配方程作损失函数的神经网络。网格不用
- 自动网格优化 — AI学网格密度最优分配
- 异常检知 — AI检出有限元结果异常值
PINN取代有限元吗？
总结
三维分析的前沿课题，整理一下。
三维弹性体分析是「完成的技术」而非「新兴计算方法」与AI融合而不断进化。但基础（弹性论、有限元原理）不变。
实体单元的EAS法体积锁定对策
体积锁定（volumetric locking）是不可压材料（Poisson比ν→0.5）用完全积分实体求解时位移被低估的现象。1990年Simo和Rifai在Stanford提出增强假定应变（Enhanced Assumed Strain，EAS）法，追加内部应变模式解除锁定。AbaqusC3D8I实装此方法，橡胶弹性体（ν≈0.499）的大变形分析威力无穷。
三维弹性体故障排除
三维分析的故障
三维实体分析常见故障教一下。
三维自由度多，问题也多。
应力异常高（奇异点）
角部应力网格细化越细越高，一直上升。
应力奇异点。几何的角（90°边、凹口前沿）弹性论上应力无穷大。网格细化后「无穷大接近」而已。
对策：
- 応力奇異点実際存在 — 実结构圆角、材料降伏
- 奇異点的応力使 — 奇異点離位置評価
- 圆角模型化 — 圆角入応力有限
- 応力集中係数 $K_t$ 評価 — 理論的 $K_t$ 比较
「网格细化应力収束」部位奇異点的可能性、。
的通。奇異点网格収束。有限元的限界弾性論的特徴。実结构的圆角半径模型化奇異性消。
内存不足
三維分析足。
三维實体分析二维的100倍以上消費。対策：
| 对策 | 效果 | 缺点 |
|---|---|---|
| 对称条件活用 | 自由度1/2～1/8 | 对称问题限定 |
| 子模型 | 局所精密化 | 整体精度粗 |
| TET10→TET4多数要素 | 自由度削减 | 精度低（非推荐） |
| 迭代法求解器 | 内存大幅削减 | 不收敛的可能 |
| Out-of-core | 磁盘活用 | 计算时间增加 |
迭代法求解器直接法少済。
直接法（LU分解） $O(n^{1.5})$ 必要、迭代法（PCG等） $O(n)$ 足。100万自由度以上問題差桁違。迭代法収束、前処理選択重要。
位移·反力不符
反力的合計荷重一致。
確認项目：
1. 拘束不足 — 剛体移動残存。三維6自由度拘束必要
2. 荷重方向誤 — 座標座標混同
3. 単位系不整合 — mm/N/MPa m/N/Pa 混在
4. 対称面反力 — 対称面的反力含合計？
「单位系的問題」三维頻繁。
CADmm出力、有限元的材料MPa（= N/mm²）入力整合。CADm出力、有限元MPa10^6。荷重入力即座反力確認鉄則。
TET4精度不足
TET4要素应力全然合。
TET4使时点解決策明白：TET10变。TET4曲变形全表現。曲支配的问题（全的実构造）桁違的誤差。
TET4要素数増精度上。
h-refinement（网格細分化）精度上、TET10同等精度TET4的1/5要素必要。自由度 $(1/5)^3 \times 4 / 10 \approx 10$ 倍必要。TET10变方圧倒的効率的。
总结
三維分析的故障対処、整理。
「TET4TET10変」三維分析最費用対効果高改善。
間違。的1的变更精度劇的改善。三維分析始、TET10使徹底。
实体单元的应力不连续诊断
实体单元间的应力不连续（stress jump）是网格密度不足的指标。相邻单元von Mises应力差超过平均值的10%时认为精度不足。Nastran的「STRAIN ENERGY DENSITY」输出用能量梯度大的区域特定重新网格化很高效。HyperMesh的Error Estimation功能用Zienkiewicz-Zhu（ZZ）推定量自动可视化局所误差。
更详细
错误