线性座屈（特征值座屈）分析 — 故障排除指南
可在 linear-buckling.html 查看更充实的内容。
特征值座屈分析的典型故障
特征值座屈分析运行时结果异常时的处理方法？
让我们整理特征值座屈特有的故障。这些与欧拉座屈故障排除不同，是特征值分析特有的问题。
特征值为零或接近零出现大量值
特征值出现0.001、0.0002等接近零的大量值…
这是前阶段静解析中没有产生有意义的应力时发生的情况。$[K_\sigma]$ 几乎是零矩阵，所以特征值会呈现极端的大或小的值。
原因检查：
1. 参考荷载为零或极小 — 确认荷载的值和方向
2. 荷载全部作用在约束点上 — 仅反力平衡，结构内部没有应力产生
3. 刚体位移未被约束 — 静解析变成奇异的
4. Ansys APDL中忘记了 PSTRES, ON — 应力刚度没有被计算
第4点是盲点。这是Ansys特有的问题吗？
是的。Nastran的SOL 105或Abaqus的*BUCKLE会自动计算应力刚度，但Ansys APDL中，PSTRES, ON 必须明确指定，否则几何刚度矩阵不会被计算。Workbench是自动的，但APDL脚本中容易忘记。
所有特征值都为负
所有模式的特征值都是负数。这意味着什么？
参考荷载的方向与意图相反的可能性很高。例如，想要压缩但给了拉伸荷载的情况。
将荷载符号反转就能得到正特征值吗？
可以的。或者保持参考荷载不变，将 $\lambda$ 解释为负也可以。$\lambda = -3.0$ 表示"反向荷载的3倍时座屈"。但这样容易引起混淆，最好从一开始就确保荷载方向正确。
第一阶和第二阶模式的特征值非常接近
$\lambda_1 = 3.42$、$\lambda_2 = 3.44$，两个模式的值几乎相同。
这就是特征值聚类。物理上意味着"两个座屈模式几乎同时发生"，存在模式相互作用导致耐力下降的风险。
对策：
- 比较模式形状 — 形状相同则是网格对称性导致的重复（无问题）。形状不同则需要注意
- 对称模式和反对称模式 — 对称结构中对称/反对称模式容易成对出现
- 非线性分析验证 — 将两个模式作为初始不完整叠加，确认相互作用的影响
在非线性分析中，如何重叠两个模式？
设定 $\{u_{impf}\} = \alpha_1 \{\phi_1\} + \alpha_2 \{\phi_2\}$，改变 $\alpha_1, \alpha_2$ 的比率进行多个工况分析。例如 $\alpha_1 : \alpha_2 = 1:0, 1:1, 0:1$ 至少需要3个工况。将最低的崩溃荷载作为设计值。
局部模式大量出现找不到整体模式
按特征值从低到高排列时，前20个模式全是翼缘局部座屈，找不到整体座屈。
这是薄壁结构中常见的问题。翼缘或腹板的板单元局部座屈模式可能具有比整体座屈更低的特征值。
对策：
1. 大幅增加求解的模式数（50～100个模式） — 整体模式被埋没在上位
2. 用位移反演法探索特定特征值域 — 在整体座屈的概算值附近设置位移
3. 子建模 — 全局模型（梁单元）评估整体座屈，局部模型（壳单元）单独评估局部座屈
子建模看起来是最清晰的方法。
完全同意。对于大规模结构（整栋建筑、桥梁全体等），将整体用详细的壳模型建模在实际上不可行。梁模型确认整体座屈，仅在关键部材用壳模型确认局部座屈的多尺度方法是实务性的做法。
预载的有无导致结果大幅变化
加入预载（自重）后座屈荷载下降了两成。为什么？
预载产生的应力会反映在 $[K_\sigma]$ 中。如果自重在柱中产生压缩，对追加荷载的座屈余度当然会变小。
反之，预载为拉伸（例如：预应力），座屈荷载会上升。预应力混凝土柱耐座屈性强的原因之一就是这个。
忽视预载进行座屈分析就会导致危险方向的过大评估吗？
不一定都是过大评估。因为有些部材的预载方向是拉伸方向。但自重这类永久荷载没有理由忽视，所以应该总是包含预载在内。这样更安全。
基准检验
有没有基准问题来验证自己的设置是否正确？
NAFEMS公开发布的基准问题最具可信度：
| 基准 | 结构 | 理论解 | 出处 |
|---|---|---|---|
| NAFEMS BM-1 | 正方形板面内压缩座屈 | $\sigma_{cr} = 4.0 \times \frac{\pi^2 D}{b^2 t}$ | NAFEMS R0015 |
| NAFEMS BM-3 | 圆筒壳外压座屈 | $p_{cr}$ 的理论式 | NAFEMS R0015 |
| 欧拉柱 | 两端铰支柱 | $P_{cr} = \pi^2 EI / L^2$ | 古典理论 |
先用欧拉柱调整，再用NAFEMS的板座屈确认是稳健的做法吧。
是的。使用新的求解器或新的单元类型时，一定要先进行基准验证再开始实际问题。座屈分析是设置错误直接影响结果的领域，不能跳过验证步骤。
终于系统地理解了座屈分析的故障排除。关键是"物理性地解释特征值的含义"吧。
完全同意。不要只看数字判断，要养成常常思考"为什么会出现这个特征值"的习惯。这样的话，问题的原因自然就能看清。
线性座屈分析结果非保守的情况
薄壁壳体（圆筒·球壳）的线性座屈荷载实验值往往大2～8倍，非常非保守（锯齿因子问题）。初始不完美和残余应力的影响极其显著，需要乘以NASA SP-8007的经验式锯齿因子（圆筒为0.6～0.9）进行修正。仅用计算值进行设计很危险，以考虑不完美影响的非线性分析为原则。
相关主题
价值
更详细
错误