4. 屈曲特征值分析
问题定位：线性摄动屈曲能回答什么
薄壁结构在面内压力下，失稳往往发生在材料屈服之前。线性摄动屈曲分析（Linear Perturbation Buckling）把问题拆成两步：先做一次静力预载，得到结构的应力刚度；再求解广义特征值问题
(K₀ + λᵢ Kσ) φᵢ = 0
其中 K₀ 为弹性刚度矩阵，Kσ 为预载引起的几何（应力）刚度矩阵，λᵢ 为第 i 阶特征值（即载荷乘子），φᵢ 为对应屈曲模态。临界载荷 = 预载 × λᵢ。本例以四边简支矩形薄板受面内均布压力为对象，提取前 6 阶模态与临界载荷因子。
建模与求解设置
在 Abaqus/Standard 中完成该分析，关键设置如下：
- 单元类型：S4R（4 节点减缩积分壳单元），厚度 2 mm，沿厚度方向 5 个积分点。
- 材料：线弹性，E = 210000 MPa，ν = 0.3。屈曲特征值分析只用到弹性刚度，不涉及塑性。
- 分析步：先建立
Static, General预载步，施加面内压力 Nx = 1 N/mm；再建立Buckle步，请求前 6 阶特征值。 - 边界条件：四边简支，即板边 U3 = 0，面内位移不约束（允许泊松收缩）。
- 网格：全局种子 10 mm，得到 60×40 网格，共 2400 个 S4R 单元。
关键操作：在 Buckle 步中，
Number of eigenvalues requested 设为 6；Minimum eigenvalue 和 Maximum eigenvalue 留空，让求解器自动搜索。若只关心最低阶，可设 Minimum = 0，Maximum = 5。
结果判读与验证
提交计算后，在 Job Monitor 中确认 THE ANALYSIS HAS COMPLETED SUCCESSFULLY。进入 Visualization 模块，切换到 Buckle 步，逐阶查看模态形状与对应的 EigenValue。
本例中一阶载荷因子 λ₁ = 3.82，意味着当预载放大到 3.82 N/mm 时结构发生屈曲。由于实际施加的预载为 1 N/mm，远小于临界值，结构在当前载荷下稳定。若把预载改为 5 N/mm，则 λ₁ = 3.82/5 = 0.764 < 1，结构将失稳。
验证方法：用经典解析解校核。四边简支矩形板在单轴压力下的临界应力 σ_cr = kπ²E/[12(1-ν²)(b/t)²]，其中 k 由边长比 a/b 决定。本例 a/b = 1.5，查表 k ≈ 4.0。代入得 σ_cr ≈ 4.0 × π² × 210000 / [12 × (1-0.09) × (400/2)²] ≈ 3.85 N/mm²。与有限元 λ₁ = 3.82 相比，误差约 0.8%，在网格收敛范围内。
常见问题与排查
| 现象 | 可能原因 | 处理方式 |
|---|---|---|
| 特征值为负 | 预载方向与压应力符号相反 | 检查预载步中压力方向，确保为面内压缩 |
| 前几阶 λ 非常接近 | 结构存在重根模态或几何对称 | 属正常现象，查看模态形状确认是否为不同变形模式 |
| λ₁ 远大于解析解 | 网格过粗或边界条件过刚 | 细化网格至 5 mm，检查简支边是否约束了面内位移 |
| 求解不收敛 | 预载步存在刚体位移 | 检查边界条件是否约束了所有刚体模态 |
注意：线性摄动屈曲分析假设结构在屈曲前处于线性弹性状态，且忽略几何缺陷。若结构在屈曲前已进入塑性，或需要追踪后屈曲下降段，应改用 Riks 弧长法并引入初始缺陷。