4. 结构力学基础与弹性波
在COMSOL中做声学—结构耦合，结构侧不是“配角”，它决定了声波如何被固体承载、反射、透射和辐射。本节把结构力学最核心的四件事串起来：应力应变关系、运动方程、模态分析，以及弹性波在固体中的传播；最后落到结构阻尼类型及其在COMSOL中的输入方法。理解这条链路，才能在“声学接口 + 固体力学接口”的耦合模型里判断：哪些参数必须给、哪些边界会改变波型、哪些阻尼设置会显著影响共振峰。
一、应力应变与材料本构
连续介质中，位移场 u = (u, v, w) 描述质点运动。小变形假设下，应变张量由位移梯度对称化得到：
ε_ij = 1/2 (∂u_i/∂x_j + ∂u_j/∂x_i)其中正应变 ε_xx、ε_yy、ε_zz 表示拉伸/压缩，剪应变 ε_xy 等表示角度畸变。应力张量 σ_ij 描述内力面密度，正应力引起体积变化，剪应力引起形状变化。对线弹性各向同性材料，本构关系为：
σ_ij = λ δ_ij ε_kk + 2μ ε_ij
λ = Eν / [(1+ν)(1-2ν)]
μ = G = E / [2(1+ν)]这里 E 是杨氏模量，ν 是泊松比，λ、μ 是拉梅常数，μ 即剪切模量 G。COMSOL固体力学接口中通常直接输入 E、ν 和密度 ρ，软件内部换算 λ、μ。需要注意：若材料接近不可压缩（ν→0.5），λ 急剧增大，体积模量主导，数值上容易出现体积锁定，此时应改用几乎不可压缩材料或混合格式。
二、运动方程与弹性动力学
把本构代入动量守恒，得到位移形式的弹性动力学方程（忽略体力）：
ρ ∂²u/∂t² = (λ + μ) ∇(∇·u) + μ ∇²u该式说明：固体中的扰动由两部分组成——体积膨胀/压缩部分以纵波（P波）传播，剪切畸变部分以横波（S波）传播。两种波速为：
c_P = sqrt((λ + 2μ)/ρ)
c_S = sqrt(μ/ρ)对典型铝材（E≈70 GPa，ν≈0.33，ρ≈2700 kg/m³），可估算 c_P≈6300 m/s、c_S≈3100 m/s；对环氧树脂（E≈3.5 GPa，ν≈0.35，ρ≈1180 kg/m³），c_P≈2200 m/s、c_S≈1050 m/s。这些数值直接决定网格尺寸：要解析波长，单元尺寸一般需小于最高关注频率对应波长的 1/6～1/10。若模型同时含空气和固体，声速差约 4～5 倍，网格与时间步必须按最快波速和最小波长共同约束。
三、模态分析与共振判据
忽略阻尼和外部载荷，令 u = U e^{jωt}，运动方程化为特征值问题：
(K - ω² M) U = 0K 为刚度矩阵，M 为质量矩阵，特征值 ω² 对应固有角频率，特征向量 U 为振型。COMSOL中“研究 → 特征频率”即求解该问题。模态分析的价值在于：它给出结构自身的共振频率和振型，帮助判断声学激励频率是否落在共振区、耦合界面是否会引起局部放大。实操中要注意：约束条件直接改变 K 的奇异性——自由—自由结构存在 6 个刚体模态（频率接近 0），必须识别并排除；固定约束会抬高频率；预应力和几何非线性会改变刚度，需在“研究”中启用相应设置。
四、阻尼类型与COMSOL输入方法
结构阻尼不是“可选项”，它决定共振峰的尖锐程度和声辐射的强弱。常见类型及在COMSOL中的对应输入如下：
- 粘性阻尼：阻尼力与速度成正比，时域运动方程为
Mü + C u̇ + K u = F。在“固体力学 → 阻尼”中可设为“粘性阻尼”，输入阻尼系数 c（单位 N·s/m）或阻尼比 ζ。适合描述流体膜、减振器。 - 滞后/结构阻尼：用复刚度
K(1 + jη)表示，η 为损耗因子。频域研究中最常用，在材料节点下设置“损耗因子 η”即可，适合金属、玻璃等小阻尼材料。 - 瑞利阻尼：
C = αM + βK，α 控制低频、β 控制高频。在“阻尼”中选择“瑞利阻尼”，输入 α、β。可由两个目标频率处的阻尼比反算：ζ_i = α/(2ω_i) + βω_i/2。 - 材料阻尼/各向异性阻尼：在材料属性中直接给各向异性损耗因子或阻尼矩阵，适合复合材料、层合板。
五、实操示例：铝梁—空气耦合的模态与阻尼设置
取一根 200 mm × 20 mm × 5 mm 的铝梁，一端固定，另一端自由，周围为空气。材料参数：E=70 GPa，ν=0.33，ρ=2700 kg/m³，损耗因子 η=0.002。操作路径：
- 建立“固体力学”接口，几何选梁，材料选铝，固定约束加在 x=0 端面。
- 添加“特征频率”研究，搜索 10 阶模态。预期前几阶为弯曲模态，一阶频率约在 1.3～1.6 kHz 量级（示例估算，实际以求解为准）。
- 在材料节点下设置损耗因子 η=0.002，或改用瑞利阻尼 α=0、β=2η/ω₁。
- 添加“压力声学，频域”接口，空气域包围梁，耦合面设为“声—结构边界”。
- 做频域扫频 500～5000 Hz，观察耦合面位移幅值和声压级。若某阶共振峰异常尖锐，说明阻尼偏小；若峰被抹平，检查是否重复设置了阻尼。
把上述关系记成一条判断链：材料参数（E、ν、ρ）→ 波速与波长 → 网格与时间步 → 模态频率与振型 → 阻尼类型与数值 → 频响峰形与声辐射。任何一步取值不合理，都会在耦合结果中以“峰位偏移、峰值失真或能量不守恒”的形式暴露出来。