3. ANSYS Workbench 屈曲分析环境
在 Workbench 中做屈曲分析，第一步不是急着点“Solve”，而是把分析系统搭对。屈曲问题在 Workbench 里通常由两个系统协同完成：Static Structural（静力结构）负责施加预载荷并求解结构的应力刚度状态，Eigenvalue Buckling（特征值屈曲）则读取前者的解，构造几何刚度矩阵并求解特征值问题。两者通过 Project Schematic 中的“Solution → Setup”连线传递数据，这条连线决定了屈曲分析是否真正建立在预应力基础上。
一、两个分析系统的分工与数据传递
Static Structural 系统承担的是“预加载”角色。它需要完成材料、几何、网格、边界条件、载荷的完整定义，并求解出一个包含应力状态的解。Eigenvalue Buckling 系统本身不重复定义这些内容，它从上游系统继承网格、材料与边界条件，只额外定义“屈曲载荷倍数”的求解设置。因此，如果上游静力分析没有收敛，或者载荷施加方式不合理，屈曲系统读到的几何刚度矩阵就是错的，特征值结果自然不可信。
在 Project Schematic 中，把 Eigenvalue Buckling 拖到 Static Structural 的 Solution 单元格上，会形成一条带箭头的连线。这条连线意味着：屈曲分析使用静力分析结束时的应力状态作为基准。若把屈曲系统单独放置、不连接上游，则它只能做无预应力的特征值分析，对大多数受压构件没有实际意义。
二、材料模型：弹性模量与泊松比是核心输入
特征值屈曲求解的是线性特征值问题 ([K] + λ[Kσ]){φ} = 0，其中 [K] 是弹性刚度矩阵，[Kσ] 是由预载荷产生的几何刚度矩阵，λ 是特征值（载荷倍数），φ 是屈曲模态。这个方程里只出现弹性刚度与几何刚度，因此特征值屈曲只依赖线弹性材料参数：弹性模量 E 和泊松比 ν。屈服强度、硬化曲线、塑性参数在特征值屈曲中不参与计算，它们只在后续的非线性屈曲（如弧长法、riks 法）中才起作用。
在 Engineering Data 中，对结构钢可给出示例参数：E = 200 GPa、ν = 0.3、密度 ρ = 7850 kg/m³。若使用铝合金，可取 E = 70 GPa、ν = 0.33。需要注意，Workbench 默认单位制下输入数值时必须与所选单位一致，否则 E 的量级错误会直接导致 λ 错误。例如在 mm-t-N-MPa 单位制下，E 应填 200000 MPa，而不是 200。
三、几何准备：细长受压构件要特别处理
屈曲分析对几何的敏感度远高于普通静力分析。对于梁、柱、桁架这类细长构件，直接用实体单元会带来两个问题：一是长细比大时单元数量激增，二是实体单元难以准确捕捉低阶弯曲屈曲模态。更合理的做法是使用梁单元（Beam）或壳单元（Shell）进行等效建模。在 Workbench 中，可通过“概念建模”或 SpaceClaim 抽取梁线、中面，再赋予截面与厚度。
几何准备还需注意：避免尖角、微小面、重叠面等导致网格畸变的特征；对薄壁结构，壁厚方向至少布置 2～3 层单元；对开口截面梁，要关注扭转屈曲模态是否被正确激发。若几何中存在对称性，可利用对称边界缩小模型规模，但必须确认最低阶屈曲模态是否破坏该对称性——若破坏，则不能使用对称简化。
四、单位制与求解器设置
Workbench 的单位制在 Project 层面统一设置，常用组合有两套：mm-t-N-MPa（毫米-吨-牛顿-兆帕）和 m-kg-N-Pa（米-千克-牛顿-帕）。前者适合中小型机械结构，后者适合大型土木或航空结构。无论选哪套，材料参数、几何尺寸、载荷数值必须自洽。一个常见错误是几何用 mm 建模，材料 E 却按 Pa 输入，导致 λ 偏差 10⁶ 量级。
Eigenvalue Buckling 的求解器设置中，关键项包括：
- Number of Modes：提取的屈曲模态阶数，一般取 4～10 阶，便于观察低阶模态是否被约束条件干扰。
- Load Multiplier：施加在预载荷上的倍数，用于构造几何刚度矩阵。通常保持默认 1.0，表示“以静力分析中的载荷为基准”。
- Solver Type：默认 Block Lanczos，适合大多数中小规模问题；模型很大时可考虑 PCG Lanczos。
- Search Range：可限定特征值搜索区间，避免求解不关心的模态。
求解完成后，在 Solution 中查看 Total Deformation 可显示各阶屈曲模态形状，查看 Eigenvalue 可得到 λ 值。临界载荷 = λ × 静力分析中施加的预载荷。例如静力分析施加 1000 N 压力，一阶特征值 λ₁ = 3.2，则临界屈曲载荷约为 3200 N。若 λ 小于 1，说明所施加的预载荷已经超过临界值，结构在当前载荷下就会失稳。
五、一个可复现的配置示例
以一根两端铰支的矩形截面钢柱为例，截面 20 mm × 20 mm，长度 500 mm，E = 200 GPa，ν = 0.3。在 mm-t-N-MPa 单位制下：几何长度 500 mm，截面 20 mm × 20 mm，材料 E = 200000 MPa，ν = 0.3。静力分析中在柱顶施加 1000 N 轴向压力，底端固定、顶端约束横向位移。连接 Eigenvalue Buckling 后，设置 Number of Modes = 6，求解得到一阶特征值。按欧拉公式估算，两端铰支柱的临界载荷约为 π²EI/L²，其中 I = 20⁴/12 ≈ 13333 mm⁴，代入可得理论值约 105 kN 量级。若 Workbench 一阶 λ 与理论值在同一量级且模态形状为半波弯曲，则说明系统搭建正确。
整个流程的关键在于：静力系统提供正确的预应力状态，屈曲系统只负责特征值求解，材料只取线弹性参数，几何要适合屈曲模态的捕捉，单位制必须自洽。把这五点确认清楚，后续的屈曲结果才有讨论价值。