2. 结构动力学建模基础

结构动力学分析的第一步不是打开软件,而是决定“用什么模型代表真实结构”。高端装备(如精密转台、卫星太阳翼、航空发动机机匣)的物理结构本质上是连续体:质量、刚度、阻尼在空间上连续分布,理论上具有无限多个自由度。直接求解连续体的偏微分方程在工程上几乎不可行,因此必须通过离散化把无限自由度问题转化为有限自由度问题。这一节的核心,就是理解连续体与离散模型之间的映射关系,并掌握自由度选择与边界条件处理这两项决定分析成败的关键决策。

连续体描述与离散化思想

连续体动力学的基本方程可写为 ρ(x)∂²u/∂t² + L[u] = f(x,t),其中 L 为微分算子,u 为位移场。它的解是空间连续函数,包含无穷多阶模态。工程上常用的离散化路线有三条:

  • 集中质量法:把连续分布质量凝聚到若干节点上,用弹簧连接。适合梁、轴、桁架等细长结构,物理直观,但高频精度差。
  • 有限元法(FEM):用形函数在单元内插值位移场,通过弱形式建立 K、M、C 矩阵。通用性强,是高端装备建模的主流手段。
  • 模态综合法:先对子结构分别求解模态,再通过界面协调组装。适合大型复杂装备的分级建模。

离散化的本质是“用有限个广义坐标近似无限维位移场”。节点越多、形函数阶次越高,模型越接近真实连续体,但计算规模也越大。因此建模不是越细越好,而是在关心的频段内保证精度。

自由度选择:决定模型能否反映真实动力学

自由度(DOF)是描述结构运动所需的最少独立坐标数。对三维实体单元,每个节点有 3 个平动自由度;对壳、梁单元还包含转动自由度。自由度选择要回答三个问题:

  1. 保留哪些自由度:与激励方向、关注模态振型相关的自由度必须保留。例如转台的回转精度主要受绕轴扭转和径向弯曲影响,轴向平动可适当约束。
  2. 如何缩减:对大规模模型可用 Guyan 静力缩聚或 Craig-Bampton 模态缩聚,把不关心的内部自由度消去,保留界面自由度与低阶模态自由度。
  3. 自由度与频率上限的关系:经验上,要准确捕捉到频率 f_max,单元尺寸应满足每个波长至少 6~10 个节点。例如钢梁中弯曲波速约 5000 m/s,若关心 2000 Hz,波长约 2.5 m,单元长度应控制在 0.25~0.4 m 量级。

自由度选择不当的典型后果是:模态遗漏、局部振型被“锁死”、或出现虚假高频模态污染结果。工程上应通过网格收敛性检查(逐步加密网格,观察前若干阶频率变化小于 2%~5%)来确认自由度是否足够。

边界条件处理:约束、连接与阻尼

边界条件是动力学模型与真实装备之间最容易出偏差的环节。常见类型包括:

  • 固定约束:如地脚螺栓连接的基础,可近似为全约束。但要注意实际螺栓连接存在局部柔度,完全刚性约束会高估整体刚度。
  • 弹性约束:用接地弹簧模拟隔振器、轴承支承。弹簧刚度需由产品手册或实测得到,不能随意取值。
  • 自由边界:用于模态试验的悬挂状态,要求悬挂频率远低于结构一阶频率(通常低于 1/5~1/10)。
  • 连接界面:螺栓、焊缝、胶接等连接刚度对整体模态影响显著。高端装备中常用“刚性连接 + 局部弹簧”混合建模。

阻尼处理同样属于边界条件范畴。工程中常用 Rayleigh 阻尼 C = αM + βK,α 控制低频、β 控制高频。α、β 可由两阶参考模态的阻尼比反算:给定 ζ₁、ζ₂ 与 ω₁、ω₂,则 β = 2(ζ₂ω₂ − ζ₁ω₁)/(ω₂² − ω₁²),α = 2ζ₁ω₁ − βω₁²。示例:若 ζ₁=ζ₂=0.02,ω₁=2π×50,ω₂=2π×500,可算得 β≈1.1×10⁻⁵,α≈12.0。这些数值属于教学示例,实际应通过模态试验辨识。

连续体 → 离散模型 → 边界条件 → 求解 的建模链路 连续体 ρ(x), E(x), 无限自由度 离散化 FEM / 集中质量 / 模态综合 K, M, C 有限维矩阵 自由度选择 平动 ux, uy, uz 转动 θx, θy, θz 缩聚:Guyan / C-B 网格收敛检查 边界条件与阻尼处理 固定约束 地脚 / 焊接 弹性约束 隔振器 / 轴承 自由边界 悬挂模态试验 Rayleigh 阻尼 C = αM + βK 特征值求解 (K − ω²M)φ = 0 → 频率与振型 模型修正 / 参数辨识

面向高端装备的建模假设与简化原则

高端装备建模的难点在于:结构复杂、连接多样、工况严苛,但计算资源与周期有限。合理的简化原则可以归纳为四条:

  1. 按关注频段取舍细节:只保留对目标频段模态有贡献的特征。小孔、倒角、螺纹在低频段可忽略,但在高频局部应力分析中必须保留。
  2. 连接刚度显式化:螺栓、法兰、轴承不能简单用“绑定”代替。应通过试验或经验公式给出等效刚度,并做灵敏度分析。
  3. 对称性与周期性利用:对轴对称或周期结构,可用单扇区模型加循环对称边界,大幅降低自由度。
  4. 子结构分级建模:把装备拆成若干子结构,分别验证后组装,便于定位误差来源。
建模假设必须记录在案:哪些部位做了刚性化、哪些连接用了弹簧、弹簧刚度来源是什么。这些信息是后续模型修正与试验对比的基础。

建模检查清单(可直接使用)

检查项判断依据示例阈值
网格收敛性加密后前 6 阶频率变化< 3%
自由度规模与计算资源匹配10⁵~10⁷
边界条件来源手册 / 实测 / 经验需标注
阻尼参数由模态试验反算ζ 0.5%~5%
连接刚度灵敏度±30% 刚度变化对频率影响< 5%