4.1 结构力学上手：悬臂梁的静力与模态

4.1 结构力学上手：悬臂梁的静力与模态

本节摘要：固体力学接口的首秀选了最经典的对账题：100 mm 悬臂梁端部受 10 N 集中力，挠度与材料力学公式对账；再加一次特征频率研究，读出四百赫兹出头的一阶模态。约束是结构仿真的头号陷阱，本节给你一套防呆配置。

从挠度公式开始倒推模型

传热那套五步工作流原封不动搬过来，先在纸上过第一步。悬臂梁：长 L = 100 mm、宽 b = 10 mm、厚 h = 5 mm，材料结构钢（取 E = 200 GPa、泊松比 0.3、密度 7850），自由端向下集中力 F = 10 N。手算先走一步：截面惯性矩 I = b·h³/12 = 10 × 5³/12 ≈ 104 mm⁴，端部挠度 δ = F·L³/(3·E·I) = 10 × 0.1³/(3 × 200e9 × 1.042e-10)，分子 0.01、分母约 62.5，得 δ ≈ 0.16 mm。固定端最大弯曲应力 σ = M·c/I = F·L·(h/2)/I ≈ 24 MPa，远低于结构钢几百兆帕的屈服强度，线弹性假设成立。这两个数字待会儿就是仿真答案的裁判。

模型向导：三维、固体力学（solid）、稳态研究。几何一个长方体 100 × 10 × 5 mm。材料库里找 Steel AISI 4340（E 约 205 GPa，与手算的 200 略有出入，对账时心里有数即可，或者干脆在 3.1 学过的空材料里手填 200 GPa，对得更干净）。网格自由四面体常规档，注意一点：弯曲问题厚度方向至少要有三层单元，5 mm 的厚度配常规四面体勉强够，稳妥起见对整个域加"尺寸"子节点、最大单元 1.5 mm。

约束与载荷：结构仿真的主战场

右键"固体力学 (solid)"→ 固定约束，选 x=0 的端面——这是全模型唯一消除刚体位移的地方，漏了它，求解器会报"奇异矩阵"然后拒绝工作，这条报错在结构仿真里的含义几乎恒等于"约束不足"。再右键 → 边界载荷，选 x=100 mm 的端面，力栏填：Fx=0、Fy=0、Fz=-10[N]。注意边界载荷默认选项是"总计"，10 N 均摊在整个端面上——这与材料力学里"集中力"的经典解有微小系统差，精确对账时可改用"单位面积"施加于端面中线附近的窄条，工程上 0.16 mm 这个量级不必较真。

计算，后处理：结果→派生值→最大值，表达式选位移场的 z 分量（w 或 solid.w），结果约 -0.16 mm，与手算仅差尾数——梁的端部效应、泊松比取值、钢材实际弹性模量把这点差异解释得干干净净。再看 von Mises 应力面图：最大值出现在固定端上表面棱角处，数值却可能离谱地大，比如几百兆帕——这不是梁要断了，是应力奇点：直角尖角加理想点约束在数学上产生无穷应力。工程处理办法是取"离开尖角一个单元尺寸处"的应力做判据，或把固定端约束改为弹性的螺栓连接模型。认出奇点、不惊慌、换位置读数，这三步是结构新手的必修礼。

图：悬臂梁的静力结果与模态振型

模态研究：同一模型换一步研究

静力模型原封不动，右键"研究 1"复制出"研究 2"，把研究步骤从稳态换成"特征频率"，要找的模态数填 6。计算结果是一张固有频率表：一阶约 408 Hz，正是悬臂一阶弯曲；对照公式 f1 = (1.875²/2π)·√(EI/(ρA·L⁴))，代入算得约 400 Hz 上下，对得上。二阶频率约在一阶的六倍出头（悬臂梁的经典倍率），约 2.6 kHz，振型变成一个波峰。模态研究的妙处是"零载荷"：它不问结构受多大劲，只问结构天生喜欢怎么振动，激励源（风扇、电机、地震）只要凑近这些频率，共振就会把微小输入放大到灾难。设计循环因此清晰：静力保证不坏，模态保证不振，两步都过才算结构过关。

本节要点回顾

• 手算先行：惯性矩、挠度、名义应力三笔账在建模前完成，仿真结果落进预期区间才算通过。
• 奇异矩阵报错 ≈ 约束不足：固定约束要消除全部六个刚体自由度，这是结构仿真第一报错。
• 应力奇点认得出：尖角与理想约束处的应力无意义，读数离开奇点一个单元尺寸。
• 边界载荷分清"总计"与"单位面积"，载荷施加方式与经典解的微小系统差要心里有数。
• 静力与模态是同一模型的两步研究，复制研究、换步骤即可，一阶弯曲 408 Hz 与公式对账。

补一个结构力学仿真的网格与边界条件实例。以散热器的翅片在气流冲击下的结构变形为例：翅片根部固定（固定约束）、翅片表面承受风压（面载荷，比如 50 Pa）。COMSOL 的设置步骤：加固体力学接口、选铝的材料参数（杨氏模量 70 GPa、泊松比 0.33）、翅片根部加固定约束、翅片表面加面载荷、网格物理场控制、稳态求解。结果关注两个量：最大变形量（如果超过 1 mm 可能碰到相邻部件）和最大 von Mises 应力（如果超过铝的屈服强度 276 MPa 会有永久变形）。这个案例展示了结构力学仿真的最小完整流程——约束定义"哪里固定"、载荷定义"受了什么力"、材料定义"什么材质"、求解给出"变形多大会不会坏"。把这四步记熟，你就掌握了结构仿真的骨架——后续只是根据实际工况替换约束和载荷的类型与数值。

结构力学仿真中的网格要求与传热有一个关键差异：结构力学关心的是应力，而应力在几何突变处（孔边、根部、截面变化处）会集中，因此这些区域的网格必须足够细才能准确捕捉应力集中系数。经验法则：应力集中区域的最小单元尺寸应为圆角半径的 1/5 到 1/10。对于散热器翅片根部的圆角（通常 0.5 mm），最小单元应为 0.05-0.1 mm——这个尺寸比传热分析的网格细了一个数量级。这解释了为什么同一个模型在传热分析和应力分析中可能需要不同的网格——传热看的是整体温度分布（粗网格够用），应力看的是局部集中（细网格必须）。COMSOL 的解决方案是使用不同的物理场网格或在同一次求解中对不同物理场使用不同网格（非保形网格），但后者需要用一致对（identity pair）连接界面。网格策略的差异是"传热简单、结构复杂"的根本原因。

一个从零到出结果的结构力学实操案例

以散热器翅片的弯曲分析为例走一遍完整流程。模型准备：从第 2 章的散热器几何出发，只保留一片翅片做简化二维平面应力分析。步骤一，物理场选择：固体力学，平面应力假设（翅片厚度远小于其他尺寸）。步骤二，材料：铝，杨氏模量 70 GPa、泊松比 0.33。步骤三，边界条件：翅片底部加固定约束（模拟与基板的焊接），翅片一侧表面加 50 Pa 的均布压力（模拟风阻）。步骤四，网格：物理场控制常规，在根部圆角处局部加密。步骤五，求解后查看：von Mises 应力云图（看峰值在哪里）、y 方向变形量（翅片尖端偏了多少）。结果解读的参考值：铝的屈服强度 276 MPa，如果 von Mises 峰值只有几 MPa，安全系数上百，说明结构绰绰有余；如果变形超过 1 mm，可能碰到相邻部件。这个案例从建模到出结果不超过十五分钟——结构力学之所以是 COMSOL 入门的首选，就是因为它的方程最简、边界条件最直观、结果最容易用常识验证。