2. 屈曲分析力学基础

屈曲不是一种“突然出现的失效”，而是结构平衡状态在载荷增加过程中发生性质改变的结果。要理解 ANSYS 中特征值屈曲、非线性屈曲与缺陷敏感性分析为什么给出不同答案，必须先把几个基本力学概念串起来：平衡路径描述结构随载荷变化的位移轨迹；分岔点与极限点标记路径上的临界事件；特征值问题给出理想结构的理论临界载荷；初始缺陷则把理想分岔问题转化为带扰动的非线性响应问题。

一、平衡路径：结构响应的“轨迹”

对一个受轴向压力 P 的细长杆，取横向挠度 w 为响应量，可以画出 P–w 平面上的平衡路径。理想直杆在 P < P_cr 时只有 w=0 的直线平衡路径，称为基本路径或主路径；当 P 达到 P_cr 时，路径出现分支，结构可以转向弯曲平衡状态，这就是分岔。平衡路径的斜率代表结构刚度：斜率为正表示稳定，斜率为零对应临界状态，斜率为负表示不稳定。

在有限元中，平衡路径由一系列平衡点组成。每个平衡点满足 K_T(u)·Δu = ΔF，其中 K_T 是切线刚度矩阵。当 K_T 奇异时，结构失去唯一平衡解的能力，对应临界点。因此，屈曲分析的核心就是追踪 K_T 随载荷和位移的变化。

二、分岔点与极限点：两类不同的临界事件

分岔点（bifurcation point）出现在理想对称结构中。以理想轴压杆为例，当 P=P_cr 时，基本路径与分支路径相交，结构存在两个可能的平衡方向。分岔点处的载荷就是经典欧拉临界载荷。分岔点又分为稳定分岔和不稳定分岔：前者在分支后仍能继续承载，后者一旦偏离基本路径就迅速失稳。

极限点（limit point）则出现在没有明显分岔的结构中，例如浅拱、球壳或受横向压力的曲板。载荷沿平衡路径上升至最大值 P_max 后开始下降，即使继续增加位移，载荷也无法再提高。极限点对应 dP/dw=0，此时结构发生跳跃失稳（snap-through）。与分岔点不同，极限点不需要理想对称条件，缺陷会直接降低极限载荷。

三、特征值问题与临界载荷

理想结构的线性屈曲分析归结为广义特征值问题：

(K₀ + λ·K_G) · φ = 0
即  K₀ · φ = -λ · K_G · φ

其中 K₀ 为小位移弹性刚度矩阵，K_G 为参考载荷下的几何刚度矩阵（初应力刚度矩阵），λ 为特征值，φ 为屈曲模态。最小特征值 λ₁ 乘以参考载荷 P_ref，即得到线性临界载荷 P_cr=λ₁·P_ref。在 ANSYS 中，这对应 ANTYPE,BUCKLE 与 BUCOPT,LANB 求解序列，输出的是特征值列表和对应的模态形状。

需要强调：特征值屈曲假设结构在屈曲前保持完美直线、材料线弹性、无初始缺陷。因此它给出的是理论上限，通常高于真实承载能力。对轴压杆，欧拉公式 P_cr=π²EI/(μL)² 与特征值解一致；对壳结构，特征值解可能比实测值高数倍。

四、初始缺陷与 imperfection 敏感性

真实结构不可避免地存在初始几何偏差、残余应力、载荷偏心与边界不理想。这些统称为初始缺陷。缺陷的作用相当于在基本路径上施加了一个横向扰动，使原本尖锐的分岔点被“抹平”，结构从加载开始就产生横向位移，平衡路径变为连续上升后趋于极限点的曲线。

缺陷敏感性可以用一个简单关系描述：设缺陷幅值为 δ，则极限载荷近似为

P_max(δ) ≈ P_cr · (1 - δ/δ₀)^(1/2)   （典型分岔结构）
或  P_max(δ) ≈ P_cr · (1 - a·δ)        （缺陷敏感结构）

其中 a 为敏感性系数。对缺陷敏感结构（如轴压圆柱壳、球壳），a 很大，微小缺陷即可导致承载能力大幅下降，工程中常用折减系数（knockdown factor）处理。对缺陷不敏感结构（如受横向压力的平板），缺陷影响较小，特征值解与非线性解差距有限。

五、从理论到 ANSYS 实现的关键对应

在 ANSYS 中，上述概念对应不同的分析类型与操作。特征值屈曲（ANTYPE,BUCKLE）求解 K₀φ = -λK_Gφ，得到 λ₁ 和模态；非线性屈曲（ANTYPE,STATIC + NLGEOM,ON + ARCLEN,ON）则追踪完整平衡路径，可识别极限点并越过下降段。初始缺陷的引入通常有两种方式：一是用特征值屈曲模态乘以一个小的比例因子作为几何扰动（UPGEOM），二是直接建立带有初始偏心的几何模型。

缺陷幅值的选取需要结合制造公差与结构敏感性。对缺陷敏感结构，建议取壁厚或跨度的千分之一到百分之一作为缺陷幅值，并做多组对比；对缺陷不敏感结构，可取更小值或直接采用特征值结果乘以折减系数。判断结果是否合理，应观察载荷–位移曲线是否出现极限点、是否收敛、以及缺陷幅值增大时极限载荷下降的速率是否符合预期。