多自由度系统的随机振动
多自由度系统的随机振动的理论基础
多自由度系统的随机振动
老师，多自由度系统的随机振动如何处理？
多自由度系统中模态间的相关性（交叉模态项）很重要。响应PSD为：
若模态充分分离，则SRSS合成已足够。对于密集模态，则需要CQC（完全二次组合）。
多点输入
多个支撑点具有相关性的随机输入：
$[S_{in}]$ 是交叉谱矩阵。对角线=自PSD，非对角线=交叉PSD。
总结
要点：
- 交叉模态项 — 密集模态时重要。使用CQC合成
- 交叉谱矩阵 — 描述多点输入的相关性
- SRSS vs. CQC — 模态分离时用SRSS，密集时用CQC
- FEM的PSD分析自动包含交叉模态项
多自由度随机振动的基础理论
多自由度（MDOF）随机振动，通过模态叠加法分解为各固有模态，求得各模态坐标下的随机响应后进行叠加。考虑模态间相关性（交叉PSD）的方法有SRSS（平方和的平方根）法和CQC（完全二次组合）法，当固有频率接近时（比率<10%）必须使用CQC。该方法由E.L. Wilson等人在1970年代于加州大学伯克利分校完成公式化。
数值解法与实现
FEM中的多自由度PSD
Nastran的SOL 111 + RANDOM自动计算所有模态的自项+交叉项。
多点输入：
```
RANDPS, 1, 1, 1, 100, 0.0 $ 输入1的自PSD
RANDPS, 1, 2, 1, 101, 0.5 $ 输入1-2的交叉PSD
RANDPS, 2, 2, 1, 102, 0.0 $ 输入2的自PSD
```
若相关性未知，则用完全相关+不相关两种工况评估范围。
总结
使用交叉PSD矩阵的解法
MDOF随机振动的精确解法是，从输入交叉功率谱密度矩阵S_FF(ω)和频率响应函数矩阵H(ω)计算输出交叉PSD矩阵S_XX(ω)＝H(ω)S_FF(ω)H*(ω)^T，并通过频率积分求得响应方差。当存在多个输入点时（例如：4点悬架），此公式化是必须的。可以使用Python的scipy库中的signal.coherence或结合numpy.fft进行数值实现。
多自由度系统的随机振动多自由度系统的随机振动实践指南
实务检查清单
在汽车底盘上的实际应用
汽车行业将路面输入作为多个轮胎（4点输入）的随机载荷处理，进行底盘的疲劳耐久评估。福特公司在1990年代基于比利时Pave测试路面的PSD数据，制定了道路模拟试验台（LST）试验标准。分析中使用MSC Nastran的随机振动SOL 111求解最大10,000自由度模型，并通过应力RMS值，用雨流计数法评估疲劳损伤量。
多自由度系统的随机振动软件与求解器比较
工具
选型指南
各求解器的MDOF支持比较
MDOF 随机振动主要求解器比较：Ansys Mechanical支持大规模模型的GPU并行化（使用RTX 4090可达传统速度8倍），MSC Nastran拥有50多年经验和丰富的NASA验证案例，Abaqus不支持随机振动（用瞬态隐式法+时程分析替代），SIMcenter Nastran可与DAKOTA联动进行概率设计优化。在HPC环境中，Nastran和AnSys Mechanical具有几乎同等的扩展性能。
尖端技术
前沿研究
非线性MDOF的等效线性化
对于包含橡胶衬套或摩擦阻尼器的非线性MDOF随机振动，等效线性化法（Equivalent Linearization, EL法）有效。1954年由Bogoliubov提出的随机EL法，用等效线性刚度和等效阻尼替换非线性恢复力，通过迭代计算收敛。Simulia研究团队在2018年发表的论文中报告，在汽车悬架非线性模型中，将响应RMS的误差控制在5%以内。
多自由度系统的随机振动常见问题与调试
故障
模态截断导致的低估
1997年发射的NASA卡西尼号探测器在振动分析中使用了200个以上模态，但在截断为150个模态的初期模型中，特定频带的响应被低估了18%。多自由度随机分析中，MIL-STD-810G也推荐包含目标频率2倍以内的固有值。