功率谱密度响应分析
功率谱密度响应的理论基础
PSD响应分析是什么
老师，「PSD响应分析」是什么？
随机振动没有确定的波形吗？
是的。地震、道路面凹凸、喷气发动机音响、火箭发射振动…这些不是确定的时间历程，而是用统计特性（PSD）来描述。
输入PSD和响应PSD
随机振动的基本公式：
- $S_{in}(f)$ — 输入PSD（加速度 $g^2$/Hz、力 $N^2$/Hz等）
- $H(f)$ — 频率响应函数（FRF）
- $S_{out}(f)$ — 响应PSD（位移 $mm^2$/Hz、应力 $MPa^2$/Hz等）
FRF的平方×输入PSD = 输出PSD。很简单。
PSD响应分析使用频率响应分析的结果。首先求FRF $H(f)$ （模态法 或 直接法），然后乘以输入PSD $S_{in}(f)$。额外计算成本几乎为零。
RMS值
响应的RMS（Root Mean Square）值是响应PSD的积分：
PSD的面积的平方根是RMS。RMS是「振动大小」的指标。
在随机振动设计中，3σ（3×RMS）是最大响应的目安。假设正态分布，有99.7%的概率振幅在3σ以内。在设计中经常使用3σ作为安全系数。
Nastran
```
SOL 111 $ 模态法频率响应
CEND
METHOD = 10
RANDOM = 20
BEGIN BULK
RANDOM, 20, , FREQ, TABRND1, 100
TABRND1, 100, , ,
, 20., 0.01, 200., 0.01, 500., 0.04, 2000., 0.04, ENDT
```
Abaqus
```
*STEP
*RANDOM RESPONSE
20., 2000., 100, 1.
*PSD DEFINITION, NAME=base_psd, TYPE=BASE
20., 0.01
200., 0.01
500., 0.04
2000., 0.04
*END STEP
```
用表格（频率 vs. PSD值）定义输入PSD呢。
在航空航天中，MIL-STD-810和NASA-STD-7001的随机振动频谱用PSD规定。把这个PSD表输入到FEM。
总结
要点：
- $S_{out} = |H|^2 \cdot S_{in}$ — FRF的平方×输入PSD = 输出PSD
- RMS = $\sqrt{\int S_{out} df}$ — 响应的有效值
- 3σ是最大响应的目安 — 正态分布的99.7%
- 频率响应分析的PSD版本 — 额外成本几乎为零
- SOL 111 + RANDOM（Nastran）, *RANDOM RESPONSE（Abaqus）
PSD（功率谱密度）的定义
功率谱密度S(f)定义为单位频率带宽的平均平方值，单位为G²/Hz或Pa²/Hz。根据Wiener-Khinchin定理（1930年代Norbert Wiener和Alexander Khinchin独立证明），PSD等于自相关函数的傅立叶变换。从实测数据推估时，Welch重叠平均法（1967年提出）是标准方法，用窗函数宽度来控制频率分辨率和方差的权衡。
功率谱密度响应的数值计算方法
PSD分析的计算
PSD响应如何计算？
1. 固有振动数分析（SOL 103 / *FREQUENCY）
2. 频率响应分析（SOL 111 / *SSD）→ 获取FRF $H(f)$
3. PSD响应计算 — $S_{out}(f) = |H(f)|^2 \cdot S_{in}(f)$
4. RMS值的计算 — $x_{rms} = \sqrt{\int S_{out} df}$
步骤3和4是FRF的后处理，不需要额外的分析。
多自由度系的PSD响应
在多自由度系中也要考虑模态间的相关性（交叉模态项）：
相关项反映了「模态的重叠」。在密集模态中相关项很重要。
输入PSD的典型例
| 规格 | 频率范围 | PSD值 |
|---|---|---|
| MIL-STD-810 Method 514 | 20〜2000 Hz | 0.04 g²/Hz（平坦部分） |
| NASA-STD-7001 | 20〜2000 Hz | 等级相关 |
| IEC 60068-2-64 | 10〜500 Hz | 产品类别相关 |
总结
Welch法的PSD推估手步骤
Welch法PSD推估为①将时间序列数据50%重叠分割（典型512～4096点），②对各段施加Hann窗进行FFT，③将FFT的绝对值平方用采样间隔和段数标准化并平均化，的步骤实施。在Python中可用scipy.signal.welch(x, fs, nperseg=1024)一行计算。在振动试验规格ASTM E1049中推荐统计自由度至少120以上（约60个或更多段），根据此设定记录时间。
功率谱密度响应的实务应用
随机振动试验的仿真
用FEM事前仿真MIL-STD-810和IEC 60068的随机振动试验。输入PSD→响应PSD→RMS（3σ）应力来评估疲劳寿命。
宇宙器材的振动评估
火箭发射时的随机振动环境。各器材安装点的PSD输入→器材内部的响应PSD→电子部件的3σ加速度。
汽车的道路噪音
路面凹凸PSD→轮胎输入→车体响应PSD→车室内振动/噪音。
实务检查清单
3σ规则是设计的基本。用RMS×3推估最大响应。
对于非正态分布（有重裾部分的分布），3σ不足，有时使用4σ～6σ。要明确分布的假设。
人工卫星结构物的音响疲劳评估
欧洲航天局（ESA）的Arianne 6火箭发射时，整流罩内的声压级最高达144dBSPL（20～2000Hz）。卫星面板的PSD响应分析用Nassif & Kao的音响-结构耦合法进行，铝蜂窝面板的面内应力RMS用Nastran SOL 111计算。在2003年的火星快车探测器中，基于这个分析结果，面板厚从1.0mm改为1.2mm，降低了音响疲劳破坏风险的设计变更。
功率谱密度响应的软件比较
PSD响应的工具
选择指南
振动试验控制器的主要产品
随机振动试验控制器市场由Spectral Dynamics（现为Schenck Process旗下）、Crystal Instruments、Vibration Research、IMC Test & Measurement四强。Crystal Instruments的Spider-80X采样最高204.8kSPS、102.4kHz带宽，PSD控制精度±1.0dB。TEAM Corporation和MTS提供大型多轴振动台（MAST：Multi-Axis Simulation Table）系统，被Bosch、Toyota、NASA采用。
功率谱密度响应的先进研究
PSD响应的先进研究
非定常随机过程的Wigner-Ville分布
像火箭发射时的振动一样，统计特性随时间变化的非定常随机过程，用常规PSD不够，需要Wigner-Ville分布（时间-频率表现）。这是Eugene Wigner在1948年为量子力学引入的概念，后来应用于信号处理，可同时获得时间和频率的分辨率。缺点是产生交叉项（克罗斯项），所以实用上用Choi-Williams分布等平滑版本。Matlab Signal Processing Toolbox 2022a以后作为wvd()函数实现。
功率谱密度响应的故障处理
PSD响应的故障
频谱泄漏（Spectral Leakage）对策
PSD推估时不用FFT窗函数（矩形窗）的情况下，会出现频谱泄漏（Spectral Leakage），相邻频率的成分混入，峰值可能被伪装最多30dB以上。对策是应用Hann窗（Hanning窗）或Flat-Top窗。Flat-Top窗振幅精度达±0.1%，但频率分辨率降低，所以在振动试验的输入控制（振动台控制）中通常用Hann窗。B&K公司的LAN-XI系统默认窗函数设为Hann窗。
价值
更详细
错误