OptFuture | 随机振动实战深度解析

在工程实践中，随机振动无处不在——车辆在路面飞驰、火箭发射升空、飞机穿越湍流、电子设备在发动机旁工作。这些真实工况下的振动并非规则的简谐运动，而是呈现出统计意义上的随机性。随机振动不同于我们熟悉的正弦运动，它的瞬时值无法用确定的数学函数描述，必须用统计方法处理。其激励通常以功率谱密度的形式定义——这并非通常意义上的“载荷”，而是一种“能量谱密度”，描述了振动能量在频域上的分布。

第一次看到PSD的单位，很多人都会疑惑：为什么是加速度的平方？为什么还要除以频率？

用“降雨量”来打个比方，RMS（均方根）就像总降雨量——它告诉你一共下了多少水。而PSD则像降雨强度（密度）——它告诉你“每平方米下了多少毫米的雨”。

分子g²：代表“能量”。在统计学里，信号的强弱用幅值的平方来衡量，对应到降雨就是“雨量”。

分母/Hz：代表“密度”。描述的是单位频率带宽，对应到降雨就是“单位面积/m²”

两者结合在一起，g²/Hz就是“单位频率范围内的能量密度”，好比“每平方米下了多少毫米的雨”，所以PSD的本质就是一张能量分布图：曲线越高，说明该频率附近的能量越“密集”。你不可能问“100Hz这个点有多少能量”，就像你不能问“第三粒盐有多咸”——你只能问“在100Hz附近这一小段频率带宽内，能量密度是多少”。

能量密度知道以后，总能量显而易见：把所有频率下的密度加起来，就是总能量——PSD曲线下的面积，就是总能量。

   

在实际工程实践中，结构不会分毫不差地传递振动——它会放大某些频率的能量（共振），也会隔离某些频率的能量（衰减），这种结构对振动信号的“过滤”特性一般用传递函数H(f)描述。在频域中，每个频率f对应一个复数H(f)，表示结构对该频率振动的放大或衰减倍数（幅值的放大倍数），相位部分则表示振动经过后的时间延迟。

由PSD单位含义可知，在随机振动中输入与输出存在平方关系，随机振动分析最核心的公式出现了：

   

   

这意味着：在随机振动中，共振的破坏力被平方级放大了。举个例子，假设输入PSD在某频率附近是0.01g²/Hz，结构在该频率的传递率是10倍，那么输出PSD在该频率附近就变成了

   

这个数值可能比输入大两个数量级，而输出响应的RMS值是对PSD曲线下的面积开方得到的，所以最终响应很可能远超过直观预期。这也解释了为什么在随机振动测试中，稍微改变一点阻尼或刚度，响应可能发生剧烈变化——因为你动的不是线性项，是平方项。

对于单自由度系统，在平直谱（白噪声）激励下，响应的均方根值可近似为：

   

其中：

   

该公式首次发表于 John W. Miles 于1954年的论文 “On Structural Fatigue Under Random Loading”（Journal of the Aeronautical Sciences），至今仍是工程手算的标准工具。

工程随机振动通常符合高斯分布：

1σ(RMS)：振动幅值有68.27%的时间在±1·RMS范围内

2σ(RMS)：振动幅值有95.45%的时间在±1·RMS范围内

3σ(RMS)：振动幅值有99.73%的时间在±1·RMS范围内

因此，在99.73%的时间里，振动的瞬间峰值都不会超过3倍的RMS，当我们对结构进行强度校核时，通常认为结构可能承受的最大瞬时加速度是3·Grms，3σ准则将动力学概率问题转化为了一个确定性的静力学设计指标。工程师可以直接用“最大可能峰值”进行强度校核，不必纠结于复杂的概率分布。

共振的“平方放大”，直接指导了工程中的两条设计原则：

控制共振是首要任务：共振点哪怕只放大了5倍，能量也会放大25倍。所以随机振动分析中，第一件事就是找出共振频率，并尽可能通过刚度调整、阻尼增加来压低传递率。

阻尼的价值被放大：Miles方程告诉我们，响应与阻尼的平方根成反比——阻尼提升带来的收益是指数级的。把阻尼提高一倍，响应仅降低约29%；但如果将阻尼提高至原来四倍，响应就能直接减半。换句话说，在随机振动中，阻尼的“投资回报率”远超直觉：花在增加阻尼上的努力，比同等程度的刚度调整更能有效控制共振。

   

图1几何模型

在消费电子领域，主板（PCB）作为核心部件，其可靠性直接影响整个系统的寿命，主板在服役过程中承受的振动通常是随机的，必须采用随机振动分析方法，该模型如图1所示，主板在三个全局方向（X、Y、Z）分别施加加速度激励，计算元器件的加速度响应及校核连接区域的结构完整性。以下为由CPU、GPU、散热器、内存芯片、圆柱形电容器、连接器组成的主板，材料参数详见表1。

表1 材料参数

   

OptFuture中执行随机振动分析的标准流程分为三个步骤：

随机振动分析建立在模态叠加法的基础上，首先要提取结构的固有频率和振型，作为后续分析的基础。在该案例中，工程师需要在模态分析中设置边界条件（如图2所示，通常是螺栓孔处的固定约束），提取足够的模态（需要覆盖频率激励范围）。

   

图2 边界条件示意图

在该案例中，根据带元器件PCB材料特性设置结构阻尼系数为0.16。

   

三个方向（X、Y、Z）分别施加加速度PSD谱，每个方向独立计算，输入加速度PSD谱参数详见图3。

   

图3 加速度PSD谱

RMS云图对比

   

   

由上述仿真结果对比可得，OptFuture与某国外商软各项仿真结果云图趋势一致，最大值误差不超过10%。基于OptFuture仿真结果对该案例进行分析（参考响应最大的方向），具体情况如下（单独元器件云图及最大值可在软件内查看）：

   

根据仿真可得，当前设计基本满足随机振动要求，PCB基板、散热器、圆柱电容器三项指标均合格，内存加速度RMS响应处于临界状态，连接器位移RMS响应超限，存在接触瞬断风险。后续优化可优先对连接器进行加固，将位移控制在0.01mm以内；将内存加速度降至15g以下，提高可靠性裕量。