線形座屈（固有値座屈）解析
理论与物理
概述
老师，“线性屈曲分析”和“特征值屈曲分析”是同一个东西吗？
是的。在FEM的语境中，“线性屈曲分析”（linear buckling analysis）、“特征值屈曲分析”（eigenvalue buckling analysis）、“线性化前屈曲分析”（linearized pre-buckling analysis）都指的是同一种方法。这是一种基于基准状态的应力构建几何刚度矩阵，然后求解特征值问题以获得屈曲载荷系数的方法。
就是欧拉屈曲那部分学过的 $([K] + \lambda [K_\sigma])\{\phi\} = \{0\}$ 对吧？是把那个应用到一般结构上吗？
没错。欧拉屈曲是针对“柱”这种特定结构的讨论，而特征值屈曲分析是一种可以应用于板、壳、框架、以及它们的组合结构的通用方法。
特征值屈曲分析的数学结构
能再详细地讲解一下一般公式吗？
分为两个步骤。
Step 1: 执行参考载荷 $\{F_{ref}\}$ 下的静力分析，求得应力分布 $\{\sigma_{ref}\}$：
Step 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$ 的比例缩放。当这个整体刚度变为奇异（det = 0）时，就是分岔点，此时的 $\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. 初始缺陷小 — 与理想形状的偏差可忽略不计
如果不满足这些条件呢？
特征值屈曲分析的结果将是上界值（unconservative estimate）。也就是说，实际的破坏载荷低于特征值屈曲载荷。具体低多少取决于结构类型：
| 结构类型 | 特征值屈曲的可靠度 | 实际破坏载荷/特征值屈曲载荷 |
|---|---|---|
| 柱（整体屈曲） | 高 | 0.85 〜 1.0 |
| 平板（面内压缩） | 较高 | 0.7 〜 0.95 |
| 圆柱壳（轴压） | 非常低 | 0.2 〜 0.5 |
| 带加筋板 | 中等 | 0.5 〜 0.9 |
圆柱壳有时会降到五分之一啊…。
所以对于圆柱壳的屈曲设计，仅靠特征值分析是不够的，必须结合引入初始缺陷的非线性分析。另一方面，柱的整体屈曲可以通过特征值分析得到相当好的估计。在实际工作中，根据结构类型把握“特征值分析的可靠度”是极其重要的。
总结
总结一下特征值屈曲分析的理论…
认识到特征值屈曲分析是“筛选工具”而非“最终答案”，这一点很重要呢。
没错。在设计的初期阶段，要快速找出“哪里危险”，这是最合适的方法。但是，请不要忘记，仅凭这个能否给出设计OK的判断，取决于结构对缺陷的敏感性。
特征值屈曲分析的起源与Hertz接触
特征值屈曲分析（线性屈曲分析）在FEM时代之前的1950〜60年代就已经被公式化为矩阵特征值问题。Turner和Clough（1960年）将FEM应用于结构分析后不久，Martin（1965年）引入了几何刚度矩阵，线性屈曲FEM由此诞生。NASTRAN（1968年）首次将其商业化实现。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。您有过急刹车时身体被向前甩出去的经历吗？那种“被带走的感觉”正是惯性力。物体越重越难启动，一旦启动也越难停止。地震时建筑物摇晃，也是因为地面突然移动，而建筑物的质量“被落下”。静力分析中此项设为零，这是基于“缓慢施加载荷，加速度可忽略”的假设。对于冲击载荷或振动问题，此项绝对不能省略。
- 刚度项（弹性恢复力）：$Ku$ 或 $\nabla \cdot \sigma$。拉弹簧时会感觉到“想要恢复原状的力”对吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——铁棒和橡皮筋，用相同的力拉，哪个伸得更长？当然是橡皮筋。这种“不易伸长”的性质就是杨氏模量 $E$，它决定了刚度。常见的误解：“刚度高=强度高”是不对的。刚度是“不易变形”，强度是“不易破坏”，是不同的概念。
- 外力项（载荷项）：体积力 $f_b$（重力等）和表面力 $f_s$（压力、接触力等）。可以这样想——桥上卡车的重量是“作用在整个内部上的力”（体积力），轮胎压路面的力是“只作用在表面上的力”（表面力）。风压、水压、螺栓紧固力…这些都是外力。这里容易犯的错误：弄错载荷方向。本想“拉伸”却成了“压缩”——听起来像笑话，但在3D空间中坐标系发生旋转时，确实会发生。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。弹一下吉他的弦试试。声音会一直响吗？不，会逐渐变小对吧。这是因为振动能量通过空气阻力和弦的内部摩擦转化成了热。汽车的减震器也是同样原理——特意吸收振动能量来改善乘坐舒适性。如果阻尼为零会怎样？建筑物会在地震后一直摇晃不停。实际上不会这样，所以设定适当的阻尼非常重要。
假设条件与适用局限
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 以mm输入时，载荷·弹性模量也需统一为MPa/N系 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意工程应变与对数应变的区别（大变形时） |
| 弹性模量 $E$ | Pa | 钢: 约210 GPa，铝: 约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm系中为tonne/mm³（钢约为 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm系和m系都统一用N |
数值解法与实现
特征值求解器的内部运作
求解 $([K_0] + \lambda [K_\sigma])\{\phi\} = \{0\}$ 的算法，具体是怎么运作的呢？
结构分析中的特征值问题通常是广义特征值问题。对于屈曲，转换为标准形式为：
$[K_0]$ 是正定对称矩阵（如果约束正确的话），$[K_\sigma]$ 是对称但不定矩阵。利用这种结构，存在几种高效的算法。
Lanczos法
听说在实际工作中Lanczos法是标准，为什么呢？
なった
詳しく
報告