2. 稳定性理论基础

在结构稳定性分析中，我们关心的核心问题不是“结构在给定载荷下应力多大”，而是“当载荷继续增加时，结构是否还能维持原有的平衡形态”。要回答这个问题，必须把结构的平衡状态放到一个随载荷参数变化的路径上去考察。这条路径就是平衡路径，而路径上的转折、分岔与极值点，正是稳定性丧失的数学标志。

一、平衡路径：把“平衡”看成一条曲线

设结构受到一个按比例加载的载荷向量 P，引入载荷因子 λ，则总载荷为 λP。结构的位移向量记为 u。有限元离散后的平衡方程可写为：

R(u, λ) = F_int(u) − λP = 0

其中 F_int(u) 是内力向量，R 是残差向量。对于给定的 λ，满足该方程的 u 就是平衡解。把所有 (λ, u) 组合画在“载荷—位移”空间中，就得到一条或多条曲线，即平衡路径。路径上每一点的切线刚度矩阵为：

K_T(u, λ) = ∂F_int/∂u

当 K_T 正定时，结构处于稳定平衡；当 K_T 奇异时，平衡路径出现临界点；当 K_T 出现负特征值时，结构已进入不稳定分支。因此，稳定性分析在数学上等价于研究 K_T 随载荷变化的性质。

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

临界点分为两类。第一类是分岔点，典型代表是理想压杆的欧拉屈曲。在分岔点之前，结构只有一条基本平衡路径（例如压杆保持直线）；到达临界载荷后，路径分叉为两条：一条是原路径的延续，另一条是新的弯曲平衡路径。分岔点的数学条件是：

det K_T = 0，且存在非零的屈曲模态 φ 使 K_T φ = 0

第二类是极限点，典型代表是浅拱或球壳的失稳。此时结构没有明显的分叉，而是载荷—位移曲线先上升后下降，在最高点处切线水平，即 dλ/du = 0。极限点之后，结构必须依靠位移控制才能追踪到下降段。极限点的数学条件是：

det K_T = 0，且 dλ/du = 0

两类临界点的区别在于：分岔点处存在多个可能的平衡方向，而极限点处平衡路径自身发生折返。在有限元中，判断属于哪一类，需要观察特征值屈曲分析得到的模态是否与基本路径正交，以及非线性追踪时路径是否出现分叉。

三、特征值屈曲：线性化框架下的临界载荷估计

特征值屈曲分析是最常用的稳定性评估手段。它的基本假设是：在临界点附近，结构处于中性平衡，总刚度由弹性刚度 K_e 和几何刚度 K_g 组成。几何刚度与当前应力状态成正比，而应力又与外载荷成正比，因此可写成：

(K_e + λ K_g) φ = 0

这是一个广义特征值问题。求解后得到一系列特征值 λ₁ ≤ λ₂ ≤ … 和对应的屈曲模态 φ₁, φ₂, …。最小的特征值 λ₁ 乘以参考载荷 P，就是经典的特征值屈曲载荷。在 Ansys Mechanical 中，这对应 Eigenvalue Buckling 分析模块：先做一次静力预载分析得到应力状态，再求解上述特征值问题。

特征值屈曲的优点是计算快、能给出屈曲模态形状，适合初步筛选和理想化结构。但它有两个明显局限：一是假设材料线弹性、结构无初始缺陷，因此结果通常偏高；二是它只能给出临界载荷，不能描述后屈曲行为。对于实际工程结构，特征值屈曲载荷往往需要乘以一个折减系数（经验上常取 0.2～0.5，具体取决于结构类型和缺陷敏感性）才能作为设计参考。

四、非线性屈曲：追踪真实平衡路径

要获得真实的屈曲载荷和后屈曲响应，必须进行非线性屈曲分析。其数学基础是增量形式的平衡方程。将载荷分成若干增量步，在第 n 步到第 n+1 步之间，对平衡方程做泰勒展开并保留一阶项，得到：

K_T Δu = λ_{n+1} P − F_int(u_n)

其中 K_T 是切线刚度矩阵，Δu 是位移增量。求解这个方程得到位移修正量，再迭代直到残差收敛。当 K_T 接近奇异时，普通载荷控制会失败，此时需要改用弧长法（Arc-Length Method），把载荷因子 λ 也作为未知量，同时约束位移增量和载荷增量的组合：

ΔuᵀΔu + Δλ² = Δs²

弧长法能够自动越过极限点，追踪到下降段。在 Ansys Mechanical 中，这通过开启大变形（Large Deflection）、设置非线性自适应载荷步、并激活弧长法来实现。判断是否发生屈曲的关键观察量包括：载荷—位移曲线是否出现水平切线或负斜率、切线刚度矩阵是否出现零或负特征值、以及求解器是否报告“刚度矩阵奇异”或“弧长法收敛困难”。

五、有限元实现逻辑：从数学到求解器

把上述理论落到有限元实现上，核心逻辑可以概括为三条主线。第一，几何刚度矩阵的组装：在每次迭代中，根据当前应力状态计算 K_g，它反映了轴向力对横向刚度的影响——压力降低刚度，拉力提高刚度。第二，切线刚度的更新：非线性分析中 K_T 由材料切线模量、几何刚度和接触刚度共同组成，每次迭代都要重新组装和分解。第三，临界点的识别：通过监测 K_T 的特征值符号变化，或者监测载荷—位移曲线的切线斜率，来判断是否接近或越过临界点。

在实际操作中，一个典型的非线性屈曲分析流程是：先做特征值屈曲得到临界载荷和模态，再以该模态为基础施加初始几何缺陷（幅值通常取结构特征尺寸的 1/1000 到 1/100），然后进行非线性静力分析，开启大变形和弧长法，逐步加载直到追踪到后屈曲路径。这样得到的极限载荷才是更接近真实结构的承载能力。