线性屈曲（特征值屈曲）分析
线性屈曲（特征值屈曲）的理论基础
概要
老师，"线性屈曲分析"和"特征值屈曲分析"是同一样东西吗？
是同一样东西。在有限元法的背景下，"线性屈曲分析"（linear buckling analysis）、"特征值屈曲分析"（eigenvalue buckling analysis）、"线性化前屈曲分析"（linearized pre-buckling analysis）都指同一种方法。以参考状态的应力构造几何刚性矩阵，求解特征值问题得到屈曲荷载系数。
这是我们在欧拉屈曲中学过的 $([K] + \lambda [K_\sigma])\{\phi\} = \{0\}$ 吗？
完全正确。欧拉屈曲是针对"柱"这种特定结构的讨论，而特征值屈曲分析是一种通用方法，可以应用于板、壳、框架，乃至这些结构的组合体。
特征值屈曲分析的数学结构
能更详细地解释一般公式吗？
分两个步骤。
步骤 1：以参考荷载 $\{F_{ref}\}$ 进行静力学分析，求得应力分布 $\{\sigma_{ref}\}$ ：
步骤 2：根据所得应力构造几何刚性矩阵 $[K_\sigma]$ ，求解特征值问题：
$\lambda_i$ 是屈曲荷载系数，$\{\phi_i\}$ 是模式形状。如果 $\lambda_1 = 3.5$ ，那就是在参考荷载的3.5倍时发生屈曲。
这就是为什么叫"线性"屈曲分析呢。
正是。这个"线性性假设"既是线性屈曲分析的本质局限，也是其计算速度快的源泉。
为什么是特征值问题
为什么屈曲可以定义为特征值问题呢？
这是个深刻的问题。屈曲是"荷载增加到某点后，变形状态突然从一种模式过渡到另一种模式"的现象。从数学上讲，这称为平衡路径的分岔（bifurcation）。
如果将荷载记为 $\lambda \{F_{ref}\}$ 的形式，那整体刚性矩阵就变成 $[K_0] + \lambda [K_\sigma]$ 。这里 $[K_\sigma]$ 是参考应力对应的几何刚性，其值与 $\lambda$ 成正比。当这个全局刚性矩阵变为奇异（行列式为零）时，就是分岔点，此时的 $\lambda$ 就是屈曲荷载系数。
$[K_0] + \lambda [K_\sigma]$ 的行列式为零的 $\lambda$ 就是要找的特征值。这就变成了特征值问题的形式！
完全正确。这与振动分析的 $([K] - \omega^2 [M])\{\phi\} = \{0\}$ 数学结构完全相同。振动中 $[M]$ 的位置换成了屈曲的几何刚性 $[K_\sigma]$ 。
几何刚性矩阵的物理意义
$[K_\sigma]$ 的物理意义是什么呢？
$[K_\sigma]$ 表示"当前应力状态对微小位移扰动做功有多大"。在压缩应力下，如果横向稍微弯曲一下，压缩力会朝着增大弯曲的方向做功。这表现为负刚性的形式。
具体来说，对于受轴力 $N$ 的梁单元，如果发生微小横向弯曲 $\delta v$ ，轴力做功约为 $N \cdot (\delta v')^2 / 2$ 。这种"应力引起的附加做功"被矩阵化后就成了 $[K_\sigma]$ 。
压缩（$N < 0$）时刚性减小，拉伸（$N > 0$）时刚性增大。这就是为什么拉伸构件不会屈曲。
直觉是对的。不过有一点注意：像预应力缆索那样的拉伸构件，在某些情况下（本质上是拉伸构件的横向屈曲，接近于横向振动现象）也可能出现屈曲问题。基本上要记住"压缩应力是屈曲的驱动力"就可以了。
线性屈曲分析的前提条件和限界
线性屈曲分析什么时候能给出"正确答案"呢？
满足以下前提条件时信度最高：
1. 屈曲前的变形很小 — 结构形状实际上没有改变
2. 材料处于弹性范围 — 屈曲点没有屈服
3. 荷载呈比例关系 — 所有荷载以相同比率增减
4. 初始不整很小 — 与理想形状的偏差可以忽略
不满足这些条件的情况会怎样？
线性屈曲分析的结果会是非保守估计（高估实际承载力）。具体高估多少取决于结构类型：
| 结构类型 | 线性屈曲的信度 | 实际失效荷载/线性屈曲荷载 |
|---|---|---|
| 柱（整体屈曲） | 高 | 0.85 ~ 1.0 |
| 平板（面内压缩） | 较高 | 0.7 ~ 0.95 |
| 圆筒壳（轴向压缩） | 极低 | 0.2 ~ 0.5 |
| 加肋板 | 中等 | 0.5 ~ 0.9 |
圆筒壳竟然会下降到五分之一？
这就是为什么对圆筒壳的屈曲设计，仅用线性屈曲分析是不够的。必须做包含初始不整的非线性分析。而柱的整体屈曲，线性分析能给出相当准确的估计。所以，要根据结构类型来理解"线性分析的可信度"。这对实务来说极其重要。
总结
总结一下线性屈曲分析的理论…
线性屈曲分析是个"初步筛选工具"而非"最终答案"，这个认识很重要呢。
正是这样。在设计初期快速找到危险部位时最好用。但最终设计是否合格，还要看结构的不整敏感性。
线性屈曲分析的起源与Hertz接触
线性屈曲分析（特征值屈曲分析）在FEM时代之前的1950～60年代就已作为矩阵特征值问题被理论化了。Turner和Clough（1960年）将FEM应用于结构分析之后，Martin（1965年）引入了几何刚性矩阵，线性屈曲FEM由此诞生。Nastran（1968年）首次将其商用化。
线性屈曲（特征值屈曲）的数值计算方法
特征值求解器的内部动作
具体怎么求解 $([K_0] + \lambda [K_\sigma])\{\phi\} = \{0\}$ 这个方程呢？
这是广义特征值问题。屈曲的情况，标准形变换为：
$[K_0]$ 是正定对称矩阵（在适当约束下），$[K_\sigma]$ 是对称但不定值矩阵。针对这种结构有几种高效算法。
Lanczos法
听说实务中Lanczos法是标准。为什么？
Lanczos法是针对大规模稀疏矩阵特征值问题的专门算法，不需要求所有特征值，只需要高效地提取指定范围内的特征值。屈曲分析只关心较小的特征值（下位数个模式），这正是Lanczos法的优势。
算法骨架是这样的：
1. 选定初始向量 $\{q_1\}$
2. 重复计算矩阵向量乘积 $[K_0]^{-1}[K_\sigma]\{q_j\}$
3. 用Gram-Schmidt正交化构造三对角矩阵 $[T]$
4. $[T]$ 的特征值是原问题的特征值近似
5. 迭代到收敛为止
直接计算 $[K_0]^{-1}$ 的话，不会变成巨大的稠密矩阵吗？
好问题。$[K_0]^{-1}$ 不是显式计算的。而是求解线性方程 $[K_0]\{x\} = [K_\sigma]\{q_j\}$ 。这样就能保持 $[K_0]$ 的稀疏性。事先对 $[K_0]$ 做LU分解（或LDLT分解），每次迭代就可以用前向后向代入来快速求解。
那LU分解是最耗时的步骤，之后的迭代相对较快，对吧。
完全正确。所以特征值问题的计算成本主要集中在矩阵分解。这与静力学分析的成本大致相同。也就是说线性屈曲分析通常只需静力学分析的1.5～2倍计算时间。
移位反演Lanczos法
如果特征值集中在某个特定值附近呢？
那就用移位反演Lanczos法来重点搜索：
如果把移位量 $\sigma$ 设为屈曲荷载的预估值，那么那附近的特征值就能高精度求出来。
实务中常见的做法是：先用无移位求出粗特征值，然后针对特定模式设定移位再精细求解。Nastran的EIGRL卡的V1、V2参数可以指定特征值搜索范围，Abaqus也有类似功能。
子空间迭代法
除了Lanczos法还有其他方法吗？
还有子空间迭代法（Bathe法）也在使用。特别是特征值密集时，有时比Lanczos法更稳健。
两种方法的特点对比：
| 特性 | Lanczos法 | 子空间迭代法 |
|---|---|---|
| 大规模问题效率 | 很高 | 较低 |
| 内存使用 | 少 | 多（需保持子空间） |
| 处理密集特征值 | 较差 | 较好 |
| 实现复杂度 | 高 | 中等 |
| 并行化容易度 | 中等 | 容易 |
Abaqus手册里有"subspace"和"lanczos"两个选项，通常lanczos就够了吧。
99%的情况下Lanczos是对的。改用子空间法通常是Lanczos收敛困难（少见）时的备选方案。
几何刚性矩阵的计算细节
$[K_\sigma]$ 的具体组建方式能细说吗？
以壳单元为例说。一个处于膜应力状态 $\{N_{xx}, N_{yy}, N_{xy}\}$ 的壳单元，其几何刚性是：
其中 $[G]$ 是将节点位移微分到空间来获得旋转增量的矩阵，$[S]$ 是应力分量矩阵。
$[S]$ 中是膜应力 $N_{xx}, N_{yy}, N_{xy}$，对吧。拉伸应力占主导的话 $[K_\sigma]$ 是正的，整体刚性增大…屈曲难…；反之压缩占主导就是负的，刚性减小…容易屈曲。
很好的理解。有趣的是，单纯的面内剪切应力 $N_{xy}$ 也能引发屈曲。这叫剪切屈曲，在飞机机翼腹板、I梁腹板等地常见。
计算效率优化技术
大规模模型的屈曲分析要快速完成，有什么技巧吗？
几个技巧：
- 超单元（部分结构法） — 对重复结构（如飞机肋骨）进行缩约减少自由度。Nastran SOL 105中常用
- 分量模式合成（CMS） — 用模态坐标表示部分结构。著名的有Hurty-Craig-Bampton法
- 对称条件活用 — 利用几何和荷载对称性只建1/2或1/4模型。但要注意反对称屈曲模式被遗漏的风险
反对称屈曲模式被遗漏？
比如对称面设置z方向位移为零，就只能出对称模式。反对称模式要在单独的对称面边界条件下（z自由、x旋转约束等）才能出现。用对称条件时必须对称、反对称都要解，这是原则。
总结
屈曲求解器的内部，有点清楚了。
要点总结：
- Lanczos法是实务标准 — 大规模稀疏矩阵最优，成本约为静力学分析的1.5～2倍
- 核心是矩阵分解 — LU/LDLT分解是最大瓶颈
- 移位反演来精细探测 — 处理特征值密集的情况
- $[K_\sigma]$ 依赖应力状态 — 膜力、弯矩、剪切都对屈曲有贡献
- 对称条件用时注意反对称模 — 必须两个都求
屈曲特征值的计算步骤和模式确认
线性屈曲分析分为三步：①进行线性静力分析取得应力②根据应力计算几何刚性矩阵Kg③求解(K+λKg)φ=0的特征值问题。特征值λ是安全系数。λ=1.0是临界点（设计荷载下屈曲），很多设计规范要求λ≥2.0。关键是不仅要看最小特征值，还要看第2、3模式，从模式形状判断实际的屈曲形态，这样才能做出设计判断。
线性屈曲（特征值屈曲）的实务应用
特征值屈曲分析的实务工作流程
在实际设计中怎样用特征值屈曲分析呢？
要根据设计阶段来用。差异很大。
| 设计阶段 | 屈曲分析的角色 | 精度要求 |
|---|---|---|
| 概念设计 | 候选断面的筛选 | 低（数量级确认） |
| 基本设计 | 屈曲荷载系数确认、危险部位特定 | 中（10%精度） |
| 详细设计 | 为非线性分析提供模式形状 | 高（网格收敛确认） |
| 设计验证 | 规范准合的屈曲检查 | 高（第三方可审核） |
概念设计要快速，详细设计要接近现实，对吧。
是的。线性屈曲分析有两种用法：直接用于设计决策，或者作为非线性分析的前阶段。后者典型做法是取屈曲模式形状作为初始不整来进行非线性分析。
参考荷载的设置
参考荷载 $\{F_{ref}\}$ 怎样确定？
初学者往往这里困惑。屈曲荷载系数 $\lambda$ 是相对于参考荷载的倍数，所以参考荷载换了 $\lambda$ 也就换了。但临界荷载 $P_{cr} = \lambda \times P_{ref}$ 始终是同一值。
比如参考荷载100 kN、$\lambda = 3.0$ 就是临界荷载300 kN。参考荷载改成300 kN的话 $\lambda = 1.0$，同样是300 kN。
完全正确。实务上最便利的做法是设计荷载（使用荷载）作为参考荷载。这样 $\lambda > 1.0$ 说明设计荷载不会屈曲， $\lambda < 1.0$ 说明会座屈。直接跟设计规范的安全系数对比。
要求安全系数3.0时，$\lambda > 3.0$ 才合格，对吧。
原则上是这样。但要记住线性屈曲是上界值的前提。不整敏感性高的结构还要乘以更多折减系数。
复合荷载的处理
实际结构是自重、外压、风荷等多种荷载同时作用。这怎样处理？
这是线性屈曲分析的重要限制。默认所有荷载等比例变化，不适用于独立变化的多个荷载。
有三种对策：
1. 多个荷载合并为一个参考荷载 — 最简单，但荷载比例固定
2. 按荷载分别进行屈曲分析 — 各荷载独立时
3. 改变荷载比例参数化求解 — 荷载比可变时最准确
"自重一定，只有风荷载变化"的情况呢？
那就用预加载+屈曲的方法。先做自重的静力分析（预加载），再基于那个应力状态分析风荷的屈曲。Nastran的SOL 105里STATSUB + BUCKLE子情形组合就是这个。Abaqus的话是Static步后加Buckle步。
预加载的应力反映在 $[K_0]$ 里，然后对新增荷载的屈曲进行评价…
模式形状的读取
屈曲荷载系数的数字只是一部分，模式形状也很重要吧。读什么？
从模式形状能看出屈曲的物理机制。检查点：
1. 屈曲类型的辨别
- 整体屈曲 — 部材整体横向倾倒。柱的整体屈曲、梁的横向屈曲等
- 局部屈曲 — 翼缘或腹板一部分波浪。板屈曲的一种
- 失稳屈曲（distortional） — 截面形状改变。薄壁冷弯型钢常见
2. 变形集中的部位
- 屈曲变形最大的地方 = 需要补强的地方
- 边界条件附近有屈曲的话，检查边界的固定度是否合理
3. 模式间的差异
- 1次和2次模式形状相似，则有模式相互作用可能
- 1次是局部屈曲、2次是整体屈曲，则局部屈曲对策优先
只看数字不看形状太危险了。屈曲荷载系数的数值、的物理的意味大事的。
太对了。$\lambda = 5.0$ 看起来很安全，但模式形状是局部板波，网格粗会遗漏真实的局部屈曲。网格细化后 $\lambda$ 急剧下降到1.2，这样的例子很常见。
设计规范的对应
建筑和土木设计规范怎样用屈曲分析结果？
规范处理异：
| 设计规范 | 屈曲的处理 | 固有值分析的位置 |
|---|---|---|
| 欧州基准3 (EN 1993) | 柱曲线(α系数)耐力低减 | 弹性屈曲荷重 $N_{cr}$ 的算定使用 |
| AISC 360 | 第E章(压缩材)屈曲应力评价 | 有效屈曲长的确认补助的使用 |
| AIJ钢结构设计规准 | 屈曲耐力式许容应力度算定 | FEM的固有值分析验证用 |
| ECSS-E-HB-32-24A (宇宙) | 适用 | 必须。理想屈曲荷重的出发点 |
欧州基准3固有值 $N_{cr}$ 出、柱曲线低减方法明示。
。欧州基准3的6.3.4「一般的方法」、FEM的弹性屈曲分析直接使屈曲耐力评价。最。
实务检查清单
屈曲分析的结果提出的检查清单。
「形状的可视化」最重要感。数字屈曲的本质。
的通。屈曲分析的价值$\lambda$ 值、形状读取物理。单「$\lambda$ 何倍」报告、分析意味半分。
集装箱船舱盖的屈曲设计
集装箱船舱盖（钢制，25×13m）受波浪、集装箱荷载、积雪的作用，线性屈曲分析是不可欠的。IACS（国际船级协会联盟）规定屈曲安全系数λ≥1.5，特局部屈曲（腹板屈曲）λ≥1.0。三菱重工业1990年代后半FEM线性屈曲舱盖全部材适用、设计工程30%短缩。
线性屈曲（特征值屈曲）的软件比较
主要求解器的特征值屈曲功能
各求解器屈曲分析的设定方法随分违？
基本的考方（静分析→固有值分析的2段阶）同、设定的书方结构违。求解器见。
Nastran（SOL 105）详解
Nastran的 SOL 105 少详教。
Nastran的线形屈曲SOL 105、Executive Control Case Control 的组合设定。
$ Executive Control
SOL 105
$ Case Control
SUBCASE 1
SUBTITLE = Static Preload (Dead Load)
SPC = 100
LOAD = 200
STRESS(PLOT) = ALL
SUBCASE 2
SUBTITLE = Buckling under Live Load
SPC = 100
LOAD = 300
METHOD = 10
STATSUB(PRELOAD) = 1STATSUB(PRELOAD) = 1 、案例1的应力引入。
付屈曲分析的。自重恒荷重案例1与、变动荷重（风荷重）的屈曲案例2评价。METHOD = 10 EIGRL 参照、固有值求解器的设定（数、方法）指定。
Nastran的屈曲分析注意：
- PARAM,BUCKLE,2 — 新有效化。古（BUCKLE,1）高精度
- EIGRL 的 ND 参数 — 求数。最低10推奖
- PARAM,POST,-1 — 结果（op2）形状书出
- 差分刚性法 vs. 荷重法 — Nastran的内部实装结果微妙变
Abaqus（*BUCKLE）详解
Abaqus的设定？
Abaqus2的顺设定：
** 1: （必要場合）
*
*STATIC
*CLOAD
node_set, 2, -1000.0
*END
**
** 2: 座屈分析
*STEP
*BUCKLE
15, , , ,
*CLOAD
node_set, 1, -500.0
*END STEP*BUCKLE 的最初的数字求数、的后的何？
eigensolver, number_of_vectors, max_iterations, blocksize 的顺。的 Lanczos 通常问题。法使明示的指定。
Abaqus的：
- 、Buckle的直前Static、的应力自动的使
- 出力 — 形状
.odb书出。$\lambda$.dat记载 *NODE FILE, GLOBAL=YES— 形状初期不整使场合.fil的出力必要- 负的固有值 — 正负两方的固有值出。正的最小值读
Ansys（Eigenvalue Buckling）详解
Ansys的Workbench中的设定手顺教。
Workbench的手顺直感的：
1. Static Structural 配置、荷重境界条件设定
2. Eigenvalue Buckling 追加、Static Structural
3. Eigenvalue Buckling 的设定 Max Modes to Find 10～20设定
4. 分析实行
5. Total Deformation 形状确认、Load Multiplier $\lambda$
WorkbenchGUI操作。APDL？
APDL的场合：
! 1: 静分析
/SOLU
ANTYPE, STATIC
PSTRES, ON ! 応力剛性有効化（重要！）
SOLVE
FINISH
! 2: 座屈分析
/SOLU
ANTYPE, BUCKLE
BUCOPT, LANB, 15 ! Lanczos法15
MXPAND, 15 ! 15結果展開
SOLVE
FINISH注意事项: APDL PSTRES, ON 忘应力刚性计算、屈曲分析无意味结果出。Workbench自动的ON、APDL手动指定必要。非常。
三求解器结果对比
同模型3的求解器解、同结果出？
理想的同、实际微妙异。要因：
| 差异的要因 | 影响度 | 说明 |
|---|---|---|
| 要素定式化的违 | 小～中 | 同「4节点」内部定式异 |
| 几何刚性的构成法 | 小 | 应力的内插方法点的违 |
| 固有值求解器的 | 极小 | Lanczos系差 |
| 网格的解释 | 小～中 | 同网格数据时差出 |
数%程度的差「正常」。
。基准测试问题（NAFEMS的标准问题）5%以内的一致得十分。10%以上、模型的解释差可能性疑。的际、形状同视觉的确认、的上 $\lambda$ 的值比较。
OptiStruct的屈曲优化
Altair OptiStruct屈曲何特别机能？
OptiStruct的最大的特征屈曲制约付优化对应。通常的优化刚性最大化目的、OptiStruct「屈曲荷重系数 $\lambda$ 指定值以上」制约追加。
屈曲考虑轻量化设计。
的通。薄肉结构的轻量化、材料削屈曲。屈曲制约付优化的自动解。航空宇宙自动车的设计特有用。、屈曲固有值的感度计算数值的不安定的、计算参数的调整必要。
选型指南
屈曲分析中心考、的求解器选？
判断的：
- 大规模航空宇宙 → Nastran SOL 105 的信赖性实绩压倒的
- Workbench完结设计业务 → Ansys Eigenvalue Buckling 最手轻
- 屈曲后的非线形分析连携 → Abaqus 的 BUCKLE → IMPERFECTION → *STATIC, RIKS
- 屈曲制约付优化 → OptiStruct
屈曲分析的求解器大差、的后的使方（非线形的连携、优化的连携）差。
的通。屈曲分析自体的求解器成熟技术。差别化要因流程全体的效率。
Simcenter Nastran的高速屈曲求解器
Siemens Simcenter Nastran大规模线形屈曲（SOL 105）100万DOF100固有值使用数时间以内计算。自动车车体的轻量薄板的屈曲分析500万DOF模型对固有值50个1夜出力2020年代确认。Lanczos法的并列化效率他的商用代码高设计者评价。
线性屈曲（特征值屈曲）的先进研究
特征值屈曲的先进应用方法
屈曲分析基本的方法思、发展的余地？
方法自体1960年代确立、使方解释的方法大研究的余地。特以下的3的方向活发。
连续体基屈曲（Continuum-Based Buckling）
从来的屈曲分析何违？
从来的FEM的屈曲分析、离散化行列的固有值问题解。一方、连续体基屈曲微分方程式屈曲条件记述、连续的屈曲场求。
具体的、的屈曲 Donnell 方程式 Sanders-Koiter 方程式数值的解。FEM的离散化误差「连续体的屈曲荷重」得。网格依存性点大利点。
、复杂形状适用？
。现状回转（圆筒、圆锥、球）的单纯形状对象。、不整敏感性高结构、FEM的网格依存性问题结构。补完的使效果的。
灵敏度分析和屈曲优化
屈曲荷重对设计参数的感度、计算？
屈曲的固有值 $\lambda$ 的设计变数 $s$ 对感度、随伴法（adjoint method）效率的计算：
形状 $\{\phi\}$ 、追加的固有值分析感度计算。
的通。屈曲优化的基盘。板厚0.1 mm变 $\lambda$ 变…的感度情报使勾配法优化回。
落穴。固有值近接（）、感度不连续。的入替（mode switching）起、优化振动的。
对处？
对策：
- MAC（Modal Assurance Criterion）追迹 — 优化的各反复的对应追迹
- 屈曲优化 — 最恶案例的固有值目的关数
- 集约关数（KS关数） — 复数的固有值滑1的关数总结：
KS关数「最小固有值微分可能」关数。巧妙。
GBT（Generalized Beam Theory）的联动
薄板冷间成形钢特殊屈曲闻？
歪屈曲（distortional buckling）。薄板沟形钢Z形钢、断面形状自体变屈曲。
一般化梁理论（GBT）、薄肉断面的屈曲「变形」分解分析方法：
- 1 — 轴变形
- 2～4 — 曲（古典的梁理论的）
- 5以上 — 歪（断面的变形）
FEM全混出的对、GBT分离理解。
。GBT的大利点屈曲的参加度（modal participation）。例「的屈曲60%歪屈曲、30%全体曲屈曲、10%屈曲」定量化。设计的判断非常有用。
的GBTUL（大学开发）CUFSM（大学开发、有限法）手轻试。设计代码（代码3、AISI S100）的DSM（Direct Strength Method）直结。
概率论屈曲评估
初期不整的考虑屈曲评价？
从来的决定论的屈曲评价（）代方法、确率论的屈曲评价注目。
手顺：
1. 实结构的初期不整测定（等）、统计的特性抽出
2. 仿真多数的不整生成
3. 各对非线形屈曲分析实施
4. 屈曲荷重的分布（平均标准偏差信赖区间）评价
何百回非线形分析回…计算。
屈曲分析使？
模型的入力变数设计段阶使。屈曲分析的形状初期不整的基底利用、各的振幅变数。屈曲分析「不整的」的道具活跃。
未来展望
屈曲分析的今后变？
3的方向见：
- 屈曲评价 — ROM（缩约模型）屈曲分析瞬时计算。结构健全性监视
- AI支援的屈曲判定 — FEM的屈曲结果不整敏感性的「补正」AI自动付与
- 等几何分析（IGA）的屈曲 — CAD完全一致形状中的屈曲分析。网格依存性的本质的排除
1960年代确立方法、AIIGA新命注入。
屈曲分析「古方法」侮。屈曲问题的、分类、优化的勾配情报、确率论的分析的基底…的前沿方法屈曲分析嵌入。基础、应用广。
确率论的线形屈曲感度分析
结构的不整（制造误差材料）屈曲荷重与影响确率论的评价「确率论的屈曲分析」2010年代确立。法1000案例的材料形状与FEM屈曲分析反复、5屈曲荷重（统计的下限值）设计荷重使。宇宙的轻量化设计重量15%削减达成事例报告。
线性屈曲（特征值屈曲）的故障对应
特征值屈曲分析的典型故障
屈曲分析回、结果的对处法教。
屈曲特有的故障整理。欧拉屈曲的故障违、屈曲分析的问题。
特征值为零或大量出现接近零的值
特征值出现了大量像 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）进行修正。仅凭计算值进行设计是危险的，原则上应进行考虑缺陷影响的非线性分析。