2. 结构动力学建模基础
结构动力学分析的第一步不是打开软件,而是决定“用什么模型代表真实结构”。高端装备(如精密转台、卫星太阳翼、航空发动机机匣)的物理结构本质上是连续体:质量、刚度、阻尼在空间上连续分布,理论上具有无限多个自由度。直接求解连续体的偏微分方程在工程上几乎不可行,因此必须通过离散化把无限自由度问题转化为有限自由度问题。这一节的核心,就是理解连续体与离散模型之间的映射关系,并掌握自由度选择与边界条件处理这两项决定分析成败的关键决策。
连续体描述与离散化思想
连续体动力学的基本方程可写为 ρ(x)∂²u/∂t² + L[u] = f(x,t),其中 L 为微分算子,u 为位移场。它的解是空间连续函数,包含无穷多阶模态。工程上常用的离散化路线有三条:
- 集中质量法:把连续分布质量凝聚到若干节点上,用弹簧连接。适合梁、轴、桁架等细长结构,物理直观,但高频精度差。
- 有限元法(FEM):用形函数在单元内插值位移场,通过弱形式建立 K、M、C 矩阵。通用性强,是高端装备建模的主流手段。
- 模态综合法:先对子结构分别求解模态,再通过界面协调组装。适合大型复杂装备的分级建模。
离散化的本质是“用有限个广义坐标近似无限维位移场”。节点越多、形函数阶次越高,模型越接近真实连续体,但计算规模也越大。因此建模不是越细越好,而是在关心的频段内保证精度。
自由度选择:决定模型能否反映真实动力学
自由度(DOF)是描述结构运动所需的最少独立坐标数。对三维实体单元,每个节点有 3 个平动自由度;对壳、梁单元还包含转动自由度。自由度选择要回答三个问题:
- 保留哪些自由度:与激励方向、关注模态振型相关的自由度必须保留。例如转台的回转精度主要受绕轴扭转和径向弯曲影响,轴向平动可适当约束。
- 如何缩减:对大规模模型可用 Guyan 静力缩聚或 Craig-Bampton 模态缩聚,把不关心的内部自由度消去,保留界面自由度与低阶模态自由度。
- 自由度与频率上限的关系:经验上,要准确捕捉到频率 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。这些数值属于教学示例,实际应通过模态试验辨识。
面向高端装备的建模假设与简化原则
高端装备建模的难点在于:结构复杂、连接多样、工况严苛,但计算资源与周期有限。合理的简化原则可以归纳为四条:
- 按关注频段取舍细节:只保留对目标频段模态有贡献的特征。小孔、倒角、螺纹在低频段可忽略,但在高频局部应力分析中必须保留。
- 连接刚度显式化:螺栓、法兰、轴承不能简单用“绑定”代替。应通过试验或经验公式给出等效刚度,并做灵敏度分析。
- 对称性与周期性利用:对轴对称或周期结构,可用单扇区模型加循环对称边界,大幅降低自由度。
- 子结构分级建模:把装备拆成若干子结构,分别验证后组装,便于定位误差来源。
建模检查清单(可直接使用)
| 检查项 | 判断依据 | 示例阈值 |
|---|---|---|
| 网格收敛性 | 加密后前 6 阶频率变化 | < 3% |
| 自由度规模 | 与计算资源匹配 | 10⁵~10⁷ |
| 边界条件来源 | 手册 / 实测 / 经验 | 需标注 |
| 阻尼参数 | 由模态试验反算 | ζ 0.5%~5% |
| 连接刚度灵敏度 | ±30% 刚度变化对频率影响 | < 5% |