3. 疲劳分析所需的光机载荷谱

光机结构在轨或运输过程中承受的载荷并不是单一形式,而是随机振动、正弦振动、冲击响应谱、温度循环和装配预紧等多种载荷的叠加。疲劳寿命分析的第一步,就是把这些分散在环境规范、试验条件和装配工艺中的载荷,统一转化为可用于寿命评估的载荷时间历程——即随时间变化的应力、应变或加速度序列。没有这一步,后续的雨流计数、S-N 曲线查取和 Miner 累积损伤都无从谈起。

本节的核心任务可以概括为:把频域或统计域的载荷描述,还原或等效为时域历程。不同载荷的转化路径不同,下面逐一说明。

随机振动:从功率谱密度到时域历程

随机振动通常以加速度功率谱密度(PSD)给出,单位 g²/Hz。它描述的是统计意义上的能量分布,本身不含相位信息,因此不能直接用于时域疲劳计算。工程上有两条常用路径:

  • 时域合成法:假设随机振动为平稳高斯过程,用 PSD 作为目标谱,通过随机相位叠加生成时域信号。典型做法是取 N 个频率点,幅值 A(f)=√(2·PSD(f)·Δf),相位在 [0,2π) 内随机取值,叠加后得到加速度时程 a(t)。为保证统计稳定性,样本时长一般不少于 60 s,并做多次样本平均。
  • 频域疲劳法:直接由 PSD 计算应力响应的谱矩,再用 Dirlik 或 Lalanne 等经验公式估算损伤。此法不需要显式时程,但要求结构响应近似线性、应力过程平稳。

实际工程中,若结构存在明显非线性(如间隙、摩擦、预紧松弛),优先采用时域合成法。合成时需注意:PSD 的斜率(如 +3 dB/oct、-3 dB/oct)决定了低频和高频能量分配,斜率越陡,时程中的低频大位移成分越突出,对镜面面形和支撑结构的影响也越大。

正弦振动:从扫频条件到等效时程

正弦振动一般给出频率范围、量级(加速度或位移)和扫频速率。它本质上是确定性载荷,转化相对直接:

a(t) = A·sin(2π·f(t)·t)
f(t) = f0 + (f1 - f0)·t / T   (线性扫频)

其中 A 为加速度幅值,f0、f1 为起止频率,T 为扫频总时长。对于定频振动,f(t) 为常数。需要注意的是,扫频过程中结构响应会随频率变化,若某阶模态被激励,局部应力可能远大于输入量级。因此正弦振动的时程必须与结构频响函数结合,得到关注点的应力时程,而不是直接用输入加速度做疲劳计算。

冲击响应谱:从 SRS 到瞬态时程

冲击响应谱(SRS)描述的是不同固有频率的单自由度系统对同一冲击的最大响应。它丢失了相位和持续时间信息,因此从 SRS 反推时程是一个欠定问题。常用做法是构造一个基波合成瞬态:

  • 选择若干频率点,使其覆盖 SRS 的主要峰值区间;
  • 为每个频率点分配一个衰减正弦子波,幅值由 SRS 反算,相位随机或按经验设定;
  • 叠加所有子波,得到候选时程,再计算其 SRS 并与目标谱对比,迭代修正幅值和相位,直到包络误差在 ±3 dB 以内。

冲击疲劳通常关注的是最大应力幅值和循环次数。对于光机结构,冲击往往发生在运输或火工品解锁瞬间,持续时间短(几毫秒到几十毫秒),但峰值高,容易在支撑点、螺钉连接处产生局部塑性变形或微动疲劳。

温度循环:从温度剖面到热应力历程

温度循环载荷以温度-时间曲线给出,例如 -20℃ → +60℃ → -20℃,升降温速率 3℃/min,高低温各保持 2 h。要转化为应力历程,需要经过热-结构耦合分析:

  1. 将温度剖面离散为若干时间步;
  2. 对每个时间步施加温度场,计算结构的热变形和热应力;
  3. 提取关注点的应力分量,形成应力时间历程。

对于光机结构,温度循环的疲劳风险主要来自材料热膨胀系数不匹配。例如铝合金镜框与玻璃镜体粘接,CTE 差异可达 20×10⁻⁶/K 以上,温度变化 80℃ 时界面剪应变可能超过胶层许用值。此时载荷谱中应包含界面剪应力和正应力的时程,而不仅仅是整体温度。

装配预紧:从拧紧力矩到初始应力

装配预紧不是交变载荷,但它决定了结构的平均应力水平,直接影响疲劳极限。螺纹连接中,预紧力 F0 与拧紧力矩 T 的关系为:

T = K · F0 · d
K ≈ 0.2(钢-钢,无润滑)

预紧力在螺杆中产生拉应力,在被连接件中产生压应力。疲劳分析时,预紧应力作为平均应力 σm 叠加到交变应力上,需用 Goodman 或 Gerber 修正。若预紧力过大,局部进入塑性;若过小,连接在振动下容易松动,导致载荷重新分配,反而加剧疲劳。

载荷谱的统一与合成

上述五类载荷最终要合成到同一时间轴上。工程上常用两种策略:

  • 分别计算损伤后叠加:对每类载荷单独生成时程、单独计数、单独查 S-N 曲线,最后用 Miner 准则累加损伤。适用于载荷之间相互独立、耦合较弱的情况。
  • 时域叠加后统一计数:将随机振动、正弦、冲击的加速度时程按时间顺序拼接,温度应力和预紧应力作为慢变分量叠加,形成一条完整的应力时程,再做雨流计数。适用于多载荷耦合明显、平均应力变化不可忽略的情况。

无论哪种策略,都必须保证应力方向一致、参考坐标系统一。例如随机振动给出的是 X、Y、Z 三向加速度,需分别投影到关注点的最大主应力方向;温度应力则需明确是拉应力还是压应力,正负号不能混用。

光机载荷谱转化路径 随机振动 PSD (g²/Hz) 正弦振动 f, A, 扫频速率 冲击响应谱 SRS (g) 温度循环 T(t), 速率, 保持 装配预紧 T, F0, K 随机相位叠加 / 频域谱矩 样本 ≥ 60 s,多次平均 a(t)=A·sin(2πf(t)t) 结合频响函数取应力 衰减正弦子波合成 迭代包络误差 ±3 dB 热-结构耦合分析 CTE 不匹配 → 界面应力 T=K·F0·d → σm Goodman/Gerber 修正 载荷时间历程 σ(t) 或 a(t) 统一坐标系 方向一致 平均应力叠加 雨流计数 → 损伤 D → Miner 累积 关键约定: ① 所有载荷统一到同一参考坐标系,应力正负号按“拉为正、压为负”定义; ② 随机振动与正弦振动取关注点最大主应力方向投影;温度应力需区分升降温阶段;预紧力作为平均应力叠加。

上图概括了五类载荷的转化路径与汇聚关系。需要强调的是,载荷谱的精度直接决定寿命评估的可信度。工程实践中,应优先采用实测数据或经过验证的仿真结果,而不是直接套用通用规范。对于光机结构,尤其要关注镜体与支撑结构连接界面、螺钉连接处和胶粘界面这三个疲劳薄弱位置,它们的载荷历程往往与整体结构响应存在明显差异,需要单独提取。

在载荷谱合成时,若随机振动和正弦振动同时存在,建议先分别计算损伤,再按 Miner 准则叠加。因为两者的循环计数分布差异很大,强行拼接时程可能导致雨流计数结果失真。
温度循环与振动载荷的耦合不可忽略。高温下材料弹性模量和屈服强度下降,相同振动量级产生的应力可能增大 10%~20%。若温度循环与振动试验同时进行,应在载荷谱中体现温度对材料参数的修正。