轴对称解析
轴对称的理论基础
轴对称问题是什么
老师，什么是「轴对称」解析？能用2次元解决3次元问题吗？
是的。当形状、材料、荷载、边界条件都相对回转轴对称时，可将3次元问题归结为2次元截面解析。在圆柱坐标 $(r, \theta, z)$ 中，$\theta$ 方向没有变化。
什么样的结构是轴对称的？
实际上非常普遍：
- 压力容器 — 圆筒体、镜板、喷嘴（周向均匀情况下）
- 法兰·螺栓 — 紧固体的轴向力解析
- 轴承 — 内圈、外圈、滚动体的接触
- 活塞·气缸 — 内燃机压力解析
- 管道内压 — 厚壁圆筒的Lamé问题
- 旋转体 — 圆盘、飞轮、涡轮转子的离心力
压力容器设计中必须的吧。
压力容器在ASME规范中也认可轴对称解析。与其用3次元建模整个压力容器，轴对称的截面解析用百分之一以下的自由度就能获得同等或更高的精度。
支配方程
轴对称问题的应力成分有哪些？
圆柱坐标系中存在4个应力成分：
平面应力和平面应变是3个成分，但轴对称是4个成分吗？
是的。轴对称的特点是存在环向应力 $\sigma_\theta$。$r$ 方向的变位 $u_r$ 会产生周向应变：
$u_r / r$！只需将变位除以半径就得到应变。这是几何效应吗？
完全正确的理解。半径 $r$ 的圆周长是 $2\pi r$，当半径增加 $u_r$ 时，圆周长变为 $2\pi(r + u_r)$。应变为 $\Delta L / L = u_r / r$。这个$1/r$ 的项是轴对称解析的核心，是平面应力/应变所没有的特性。
本构关系
轴对称的Hooke定律：
与3次元的Hooke定律形式相同呢。虽然4个成分，但本质上是3次元。
是的。轴对称是「外观是2次元，内部是3次元」。应力状态完全是3次元的，$\sigma_\theta$ 不会变为零（除了 $r = 0$ 处）。
厚壁圆筒（Lamé问题）
受内压的厚壁圆筒有理论解吗？
Lamé问题是轴对称最基本的理论解。当内半径 $a$、外半径 $b$、内压 $p$ 时：
内表面（$r = a$）处环向应力最大，向外半径处逐渐减小。内表面是拉伸，向外变小。
内表面的最大环向应力：
薄壁情况（$b \approx a$）下接近 $\sigma_\theta \approx pD/(2t)$（薄壁圆筒公式）。这是验证FEM轴对称解析的最佳基准问题。
小结
我来总结轴对称解析的理论。
要点：
- 形状、荷载、材料相对回转轴对称就能用2次元求解 — 自由度降低至百分之一以下
- 4个应力成分 — $\sigma_r, \sigma_\theta, \sigma_z, \tau_{rz}$。环向应力 $\sigma_\theta$ 是独有的
- $\varepsilon_\theta = u_r / r$ — 几何学周向应变
- 本质是3次元应力状态 — 外观是2次元但内部是3次元
- 用Lamé问题验证 — 厚壁圆筒的内压是最佳基准
压力容器设计中不用轴对称就没有理由了。成本只有3次元的百分之一而精度相当。
完全同意。但当有非轴对称荷载（喷嘴荷载、地震荷载）或非轴对称形状（孔、不连续部分）时就需要3次元解析。先用轴对称进行筛选，然后只对必要部分进行3次元详细解析是高效的方法。
轴对称的历史起源
轴对称解析的理论基础可以追溯到1850年代G.B.Airy引入应力函数。1960年代Clough & Rashid定式化了环形单元，将浮点运算降至3次元解析的约1/10。在压力容器和原子反应堆安全壳的设计中至今仍是现役手法。
轴对称的数值计算方法
FEM的轴对称解析
轴对称的FEM单元有什么特点？
轴对称单元是 $(r, z)$ 平面的2次元单元，但内部定式化中有$2\pi r$ 的权重。这对应于单元绕回转轴旋转一周的体积。
单元的定式化
轴对称单元的刚度矩阵积分：
平面问题中的 $t$（板厚）被替换为 $2\pi r$（周长）。
$r = 0$（回转轴）处 $2\pi r = 0$，这没问题吗？
非常好的问题。$r = 0$ 处的积分具有奇异性。实现中采用以下对策：
- Gauss积分点不设置在 $r = 0$ 上（积分点在单元内部）
- $r = 0$ 上的节点处 $\varepsilon_\theta = u_r / r$ 变为 $0/0$，但用洛必达法则得 $\varepsilon_\theta = \varepsilon_r$
- 实际上没有问题，但 $r = 0$ 附近用极细网格时会数值不稳定
求解器单元名称对照
| 单元 | Nastran | Abaqus | Ansys |
|---|---|---|---|
| 4节点四边形 | CQUADX4 / CTRIAX | CAX4, CAX4R, CAX4H | PLANE182 (KEYOPT3=1) |
| 8节点四边形 | CQUADX8 | CAX8, CAX8R | PLANE183 (KEYOPT3=1) |
| 3节点三角形 | CTRIAX6 | CAX3 | PLANE182（退化） |
Abaqus的CAX是「Continuum AXisymmetric」的缩写吧。CPS是平面应力，CPE是平面应变，CAX是轴对称。
完全正确。Ansys用同一PLANE182/183单元通过KEYOPT来切换。Nastran是分别的CQUADX系单元。
边界条件的设置
轴对称解析需要注意的边界条件是什么？
1. 回转轴（$r = 0$）的约束
- $r = 0$ 上的节点必须设 $u_r = 0$（径向变位为零）
- $u_z$ 自由（轴向可以移动）
- 遗漏此条件回转轴会横向移动，结果非物理
2. 对称面的处理
- 轴向对称面（如圆筒中心截面）处 $u_z = 0$
- 这样可以用半模型
3. 荷载输入
- 集中力是「断面上的线荷载」而非「$2\pi r$ 除以的值」
- 内压作为面压直接输入（$2\pi r$ 自动考虑）
集中力输入最容易出错吧？
是的。轴对称解析中在 $r = 50$ mm 的节点输入「100 kN的集中荷载」时，实际上相当于周向分布的 $100 / (2\pi \times 0.05) = 318$ kN/m 的线荷载。输入值变成「环形荷载」，不理解这一点荷载量级会对不上。
非轴对称荷载的处理（傅里叶展开）
荷载不是轴对称的情况怎么办？
用傅里叶展开可以分离非轴对称成分。将荷载按周向高调波（$\cos n\theta, \sin n\theta$）分解，对每个高调波进行2次元解析，最后叠加。
Nastran可用CQUADX + SOL 101 进行傅里叶展开解析（谐波单元）。Abaqus用 CAXA 系单元（非轴对称-轴对称）处理非轴对称荷载。
喷嘴荷载这样的局部非轴对称荷载也能用傅里叶展开？
原理上可以，但局部荷载用傅里叶级数表示需要大量高调波，计算成本接近3次元。喷嘴荷载这样的局部问题用3次元子建模更高效。
小结
轴对称的数值方法，总结一下。
要点：
- $2\pi r$ 权重进入刚度矩阵 — 与平面单元的根本不同
- $r = 0$ 上的节点 $u_r = 0$ — 别忘回转轴约束
- 荷载作为环形荷载 — 集中荷载的意义与平面问题不同
- 用傅里叶展开处理非轴对称荷载 — 但局部荷载用3D更高效
- 混合单元 — 非压缩材料（橡胶）的轴对称解析也需要
轴对称单元的数值积分
轴对称单元的刚度矩阵含有r方向的半径权重，需要与平面单元不同的Gauss积分法则。Nastran的CTRAX单元采用2×2 Gauss点，从1972年版本起就实装了原点（r=0）附近的数值不稳定回避特殊处理。
轴对称的实务应用
轴对称解析的实务应用
轴对称解析最主要的应用领域是什么？
压力容器的设计。ASME BPVC Section VIII Div. 2中，Design by Analysis方法中广泛使用FEM轴对称解析。
压力容器的轴对称解析
压力容器的哪些部分用轴对称解析？
| 部位 | 可轴对称解析否 | 注意事项 |
|---|---|---|
| 圆筒体 | ○ | 基本情况 |
| 半球镜板 | ○ | 与圆筒接续部的不连续应力 |
| 半椭圆镜板 | ○ | 滑口部的高应力 |
| 法兰 | ○ | 螺栓荷载用等价环形荷载表示 |
| 喷嘴（轴向） | △ | 轴对称则○，偏心喷嘴则× |
| 支撑脚·鞍座 | × | 非轴对称。需要3D |
镜板与圆筒的接续部很重要呢。
是的。半椭圆或碟形镜板的滑口部（曲率变化的部分）产生不连续应力。这种应力用膜理论算不出来，只有用FEM轴对称解析才能准确评估。
实务中的要点：
- 滑口部网格应小于板厚的1/4
- 板厚方向最少4个单元（二阶单元则2个）
- 焊接部的错边（目差）作为初始不整引入
O形圈沟的接触解析
轴对称也能解接触问题吗？
O形圈的压缩变形是典型的轴对称接触问题。橡胶O形圈被压入沟内产生变形。
要点：
- 橡胶的超弹性本构 — Mooney-Rivlin、Ogden等。线性弹性不准
- 混合单元（CAX4H） — 非压缩材的锁定对策
- 接触面定义 — O形圈与沟面接触。含摩擦
- 大变形解析 — O形圈压缩率10～30%较大
2次元能解接触就比3次元快太多了呢。
O形圈设计时往往要评估100多个参数组合（沟深、沟宽、O形圈径、材料硬度）。3次元无法做到，但轴对称每个情况数秒就算完。
旋转体的离心力解析
涡轮盘的离心力解析也是轴对称吗？
盘本身是轴对称，但叶片是非轴对称。实务中：
1. 盘用轴对称解析 — 将离心力 $\omega^2 r \rho$ 作为体积力输入
2. 叶片离心力作为等价环形荷载输入 — 叶片数×离心力 / $2\pi r$
3. 关键部位（螺栓孔等）用3次元子模型详细解析
叶片作等价环形荷载是近似，盘的全体应力有足够精度吗？
盘的轴孔部（中心孔周围）和腹板部的应力用轴对称就足够准确。叶片安装部的榫槽是局部非轴对称，需要3次元子模型。
实务检查表
请给出轴对称解析的检查表。
「结果是全周还是截面」是什么意思？
轴对称解析的反力是截面值，要换成全周时乘以 $2\pi r$。比如 $r = 100$ mm 处节点反力 $F_r = 50$ N，全周力是 $50 \times 2\pi \times 0.1 = 31.4$ N…不对，轴对称单元通常反力已含全周值。求解器不同处理也不同，要查手册确认。
涡轮泵解析的实务
H-IIA火箭LE-7A发动机涡轮泵解析中，轴对称模型在初期设计阶段大量使用。在3次元模型构筑前，用轴对称对板厚、压力、旋转荷载进行试算，据JAXA内部报告可削减设计工时约40%。
轴对称的软件比较
轴对称解析的工具
轴对称解析特别强的工具有吗？
通用FEM都支持轴对称，但压力容器设计用专用工具更方便。
压力容器专用工具
| 软件 | 特点 |
|---|---|
| PVElite | 压力容器规范计算（ASME、EN 13445等）。不是FEM而是规则库 |
| Compress | 符合ASME Section VIII的设计计算 |
| Nozzle Pro | 喷嘴接续部应力解析。内置轴对称FEM |
| FEPipe | 配管柔性解析和局部应力评估 |
Nozzle Pro内置轴对称FEM呢。
Nozzle Pro除了WRC 297/537计算，还内置轴对称FEM自动计算不连续应力。还自动进行ASME Div. 2的应力分类（膜应力、弯曲应力）。对压力容器工程师比通用FEM更高效。
通用FEM
| 观点 | Nastran | Abaqus | Ansys |
|---|---|---|---|
| 轴对称单元 | CQUADX4/8 | CAX4/8, CAXA | PLANE182/183 |
| 傅里叶展开 | 谐波单元 | CAXA（有限制） | PLANE25 |
| 非压缩材轴对称 | 有限 | CAX4H, CAX8RH | 混合对应 |
| 接触（轴对称） | 对应 | 对应（最灵活） | 对应 |
Abaqus的傅里叶展开「有限制」吗？
Abaqus的CAXA单元傅里叶展开的高调波数有限制。Nastran和Ansys的谐波单元傅里叶展开更灵活。不过多数情况非轴对称荷载用3次元子建模更实际。
选型指南
总结一下？
压力容器用专用工具比通用FEM与ASME规范的整合更容易呢。
ASME Div. 2的应力分类（一般膜、局部膜、一次弯曲、二次应力）需从通用FEM结果手动分类。专用工具自动化了这个，大大减轻设计者负担。
主要求解器的轴对称对应
Abaqus从Rev 3.0（1984年）开始同时提供轴对称单元CAX4R的全积分和缩减积分。MSC Nastran的CTRIAX6和Abaqus CAX6M采用中点节点的6边形，在应力集中部精度不同，已成为现场的标准配置选择。
轴对称的先端研究
轴对称解析的先端课题
轴对称解析的先端研究有什么？
轴对称本身是成熟技术，但应用面有发展。
2.5次元解析（傅里叶展开的高度化）
「2.5次元」解析是什么？
例如轴承接触问题。内外圈是轴对称，但荷载（径向荷载）非轴对称。用傅里叶展开分离非轴对称接触力，各高调波分别进行轴对称接触解析再叠加。
如果能用2.5次元代替3次元接触，计算成本大幅下降呢。
能期待1～2个数量级的加速。旋转体接触如轴承和齿轮，这种方法收益大。
从轴对称到3次元的子建模
最近实务中轴对称全局模型 → 3次元局部模型的子建模成为标准。
步骤：
1. 用轴对称解析压力容器全体
2. 抽取喷嘴接续部等非轴对称部位的变位
3. 作为3次元子模型的边界条件使用
4. 用子模型详细评估局部应力
Abaqus和Ansys的子建模功能能自动化吗？
Abaqus的 *SUBMODEL 关键字支持轴对称→3次元子建模。自动将轴对称模型的 $(r, z)$ 坐标映射到3次元的 $(x, y, z)$。Ansys的Workbench子建模功能也能做类似的事。
氢脆化·扩散耦合
压力容器相关的先端课题？
高压储氢容器的设计是注目领域。氢在材料中扩散引起脆化（氢脆化）这一现象与应力解析耦合。
应力场→氢的应力驱动扩散→脆化导致材料劣化→应力重分配…这样的连锁反应用FEM追踪。压力容器大部分是轴对称，所以这个耦合解析能在2次元高效实施。
是与氢社会实现直接相关的研究呢。
燃料电池汽车（FCV）的氢罐（700气压）和管道的氢输送中，轴对称FEM耦合解析正作为设计工具使用。
小结
轴对称的先端课题总结一下。
轴对称不是「古老的手法」而是最先端的耦合解析的基础。
傅里叶级数轴对称扩展
对有非轴对称荷载的结构，将荷载傅里叶级数展开，各调和成分分别作为独立轴对称问题求解的手法由1970年代Wilson等人确立。ADINA-A求解器以此手法进行圆筒形储罐地震应答解析，计算时间相比3次元解析缩短至1/20。
轴对称的故障排除
轴对称解析的故障
轴对称解析的常见故障请教一下。
轴对称特有的盲点有几个。
回转轴上的约束遗漏
$r = 0$ 上的节点忘记设 $u_r = 0$ 会怎样？
回转轴会横向移动，应力分布变成非物理的。有时还会出现奇异刚度矩阵错误（Nastran的FATAL 2012等）。
很多求解器对轴对称单元的 $r = 0$ 上的节点自动设 $u_r = 0$，但有的求解器需要手动设置。一定要确认。
环向应力的遗漏
结果只看 $\sigma_r$ 和 $\sigma_z$ 忘了看 $\sigma_\theta$ 的情况很多。
轴对称中 $\sigma_\theta$（环向应力）往往是最大应力。内压受压的圆筒中 $\sigma_\theta > \sigma_z > \sigma_r$。遗漏 $\sigma_\theta$ 是严重的，会导致最大应力的重大低估。
检查von Mises应力时 $\sigma_\theta$ 自动含入，但应力分类（ASME Div. 2的膜/弯曲分离）需分别确认。
荷载输入单位错误
轴对称荷载输入要注意什么？
求解器对荷载的解释不同：
- Abaqus — 集中力是全周值（$2\pi r$ 分）。输入 $F = 100$ N 则全周100 N
- Nastran — FORCE卡是全周值。$F = 100$ N = 全周100 N
- Ansys — 节点集中力是全周值
所有求解器都是全周值呢。
「环形荷载」的意思。内压100 MPa作为面压给时所有求解器都处理正确，但给集中力时以全周的合计力输入要记住。
$r = 0$ 附近的应力集中
回转轴附近应力异常大。
$r \to 0$ 时 $\varepsilon_\theta = u_r / r$ 分母接近零，数值精度下降。特别是网格粗糙时 $r = 0$ 附近的单元应力易波动。
对策：
- $r = 0$ 附近网格要足够细
- $r = 0$ 上的应力作为「参考值」，不用于设计判断
- 需要时用3次元模型验证 $r = 0$ 附近
与Lamé理论解的对比
FEM结果与Lamé理论解不符时？
检查项目：
1. 开放端还是密闭端 — Lamé问题假设平面应变（无限长圆筒）。密闭端时 $\sigma_z$ 会变
2. 端部条件 — FEM模型端部（$z$ 方向边界条件）是否与理论一致
3. 网格密度 — 肉厚方向最少4个单元（二阶单元则2个）
4. 薄壁近似的混淆 — 薄壁公式 $\sigma_\theta = pD/2t$ 是Lamé的近似。$D/t < 20$ 时差异大
薄壁公式与Lamé式的差在 $D/t$ 多少时出现？
$D/t < 20$ 时差异超过5%。$D/t < 10$（厚壁）时必须用Lamé式。FEM当然都能准确解，但验证时别搞错用哪个理论解。
小结
轴对称的故障处理总结一下。
环向应力的遗漏最可怕呢。$\sigma_\theta$ 是最大应力还没发现的话，设计就危险。
压力容器事故多是环向应力引起的破裂。$\sigma_\theta$ 是轴对称解析中最要关注的成分，要时刻铭记。
轴线奇异点的故障
轴对称解析中在对称轴r=0上配置节点时，环向应力σθ有时不与σr一致而产生「伪应力」。Ansys PLANE25单元中，不对对称轴全节点强制Ux=0的边界条件时，1980年代初期版本报告过最大15%误差。
相关课题
价值
更详细
错误