3. 结构动力学与频响函数
随机振动疲劳寿命仿真的核心任务，是把输入端的加速度功率谱密度（PSD）经过结构动力学系统，映射为危险点的应力功率谱密度，再结合材料S-N曲线与疲劳损伤模型估算寿命。这条链路能否成立，取决于一个关键中间环节：结构动力学模型是否真实反映了从基础激励到局部应力的传递关系。模态分析、频响函数（FRF）与传递函数，正是描述这一传递关系的三件核心工具。
一、从物理空间到模态空间：模态分析在做什么
一个连续结构经有限元离散后，动力学方程为：
[M]{ẍ} + [C]{ẋ} + [K]{x} = {F(t)}其中 [M]、[C]、[K] 分别为质量、阻尼、刚度矩阵，{F(t)} 为外载荷向量。当结构自由度达到数十万甚至数百万时，直接求解该方程在频域内代价极高。模态分析通过求解无阻尼特征值问题：
([K] − ω²[M]){φ} = 0得到固有频率 ωr 与振型 {φ}r。振型具有关于质量矩阵和刚度矩阵的正交性，利用这一性质可将耦合的物理方程解耦为一组独立的单自由度方程。这一步的意义在于：结构在任意激励下的响应，可以表示为各阶模态响应的线性叠加。对于随机振动疲劳而言，这意味着我们只需关注那些被激励频带覆盖、且对危险点应力贡献显著的模态，而不必求解全部自由度。
建模要点上，模态提取阶数应覆盖激励PSD有效频带的上限，通常建议提取到分析频带上限的1.5～2倍。若只提取到200 Hz而激励在500 Hz仍有能量，高频模态的缺失会直接导致应力响应被低估。此外，模态分析阶段应检查是否有局部模态（如薄板、支架、连接件）落在激励频带内，这类模态往往对应力热点贡献极大，却容易被整体网格疏漏。
二、频响函数与传递函数：从激励到应力的映射算子
频响函数 H(ω) 定义为系统在单位简谐激励下的稳态响应，是结构动力学特性的频域表达。对于基础加速度激励 a(ω) 到某点应力 σ(ω) 的映射，可写为：
σ(ω) = H_σa(ω) · a(ω)这里 H_σa(ω) 就是加速度到应力的传递函数。它既可以通过模态叠加法由模态参数合成，也可以由有限元频响分析（如直接法或模态法）直接计算。在随机振动中，由于输入是功率谱而非确定信号，映射关系作用于PSD：
S_σσ(ω) = |H_σa(ω)|² · S_aa(ω)这一式子揭示了随机振动疲劳仿真的核心逻辑：应力PSD等于加速度PSD乘以传递函数模的平方。因此，传递函数的幅频特性——尤其是峰值频率、峰值高度和带宽——直接决定了应力响应的量级与分布。若传递函数在某个频率出现尖锐峰值，而激励PSD在该频率也有能量，应力响应会被显著放大，疲劳损伤随之集中。
三、传递函数的建模要点与常见陷阱
传递函数的准确性取决于三个层面：模态参数、模态截断和应力恢复。
模态参数层面，固有频率的误差会直接导致传递函数峰值偏移。若仿真频率比实际偏低，而激励谱在对应频段有能量，应力响应会被高估；反之则低估。阻尼比则控制峰值高度和带宽，随机振动中结构阻尼通常取1%～3%，但装配界面、螺栓连接处的摩擦阻尼可能显著更高，需要根据试验模态或经验数据修正。阻尼比从1%增大到3%，共振峰幅值约下降至原来的1/3，对应应力PSD峰值下降近一个数量级，对寿命影响极大。
模态截断层面，模态叠加法只保留有限阶模态。若截断频率以下模态不完整，或高频模态对危险点应力有不可忽略的贡献，传递函数会在截断频率附近出现虚假衰减。工程上常用模态有效质量参与系数来判断截断是否充分：在激励方向上，累计有效质量达到总质量的80%～90%以上，通常认为截断合理。对于局部应力热点，还应单独检查该区域附近模态是否被保留。
应力恢复层面，频响分析得到的是节点位移或模态坐标，需要经过应力恢复矩阵才能得到单元应力。对于壳单元和实体单元，表面应力与内部应力可能差异显著，疲劳危险点通常位于表面。若直接使用节点平均应力，可能平滑掉应力集中。建议在危险区域细化网格，并提取单元积分点或表面层的应力分量，再合成等效应力。
四、从传递函数到应力PSD的实操检查
完成频响分析后，不要急于进入疲劳计算，应先检查传递函数本身是否合理。具体可执行以下检查：
- 峰值频率核对：将传递函数峰值频率与模态分析固有频率对比，确认峰值对应的是哪一阶模态，是否落在激励PSD的有效频带内。
- 峰值幅值量级：在共振频率处，传递函数幅值近似为 1/(2ζ)，若阻尼比取2%，峰值约为25倍。若仿真得到数百倍的峰值，需检查阻尼设置或模态耦合是否异常。
- 低频段趋势：在远低于第一阶固有频率的频段，基础加速度到应力的传递函数应趋于一个常数（准静态刚度比）。若低频段出现异常上升或下降，可能是约束条件或单位制不一致。
- 应力PSD积分：对应力PSD在分析频带内积分得到应力均方值，结合峰值因子可估算应力幅值分布，与经验判断对比是否在合理范围。
综上，模态分析提供解耦基础，频响函数与传递函数建立从加速度到应力的频域映射，二者共同决定了应力PSD的形状与量级。建模时应围绕激励频带覆盖、阻尼合理性、模态截断充分性和应力恢复精度四个要点展开检查，确保后续疲劳损伤估算建立在可信的动力学传递关系之上。