2. 热循环载荷谱与工况解析
热循环载荷谱是把真实服役环境中的温度—时间历程,抽象为可用于试验、仿真与寿命评估的定量描述。它回答四个基本问题:温度在什么范围内变化(温度幅值)、变化有多快(升降温速率)、在高温或低温端停留多久(保温时间)、这样的过程重复多少次(循环次数)。在此基础上,还要考虑构件内部因温度分布不均产生的热梯度,以及多轴约束下形成的多轴热应力。只有把这几项参数同时刻画清楚,载荷谱才具备工程可用性。
一、载荷谱的五个核心参数
温度幅值通常用最高温度 T_max、最低温度 T_min 和温差 ΔT = T_max − T_min 表示。工程上更关注构件关键部位的实际金属温度,而非环境温度。例如涡轮盘轮缘与轮心在同一时刻可能相差 200 ℃ 以上,若只用燃气温度作为输入,会严重低估热应力。
升降温速率用 ℃/min 或 ℃/s 表示,直接决定瞬态热梯度的大小。速率越高,构件表面与内部温差越大,热应力越显著。航空发动机起飞—巡航—着陆过程中,涡轮叶片前缘的升温速率可达数百 ℃/s,而地面燃气轮机启停相对平缓,通常在几十 ℃/min 量级。
保温时间指在最高温度或最低温度附近的持续时长。保温阶段温度梯度趋于平缓,但高温端会叠加蠕变、氧化和扩散损伤,因此保温时间对总损伤的贡献不可忽略。
循环次数决定损伤累积的总量。工程上常把载荷谱分解为若干典型循环块,如“起动—运行—停机”为一个大循环,内部再叠加若干小循环,用雨流计数法统计各幅值对应的循环数。
热梯度与多轴热应力是载荷谱从“温度”走向“应力”的关键环节。温度梯度通过热膨胀系数 α、弹性模量 E 和约束条件转化为热应变与热应力。当构件同时承受离心力、气动载荷和热载荷时,主应力方向随时间和位置变化,形成多轴应力状态,此时需用等效应力或临界平面法进行损伤评估。
二、从温度谱到应力谱的转换
温度谱本身不是损伤的直接输入,必须经过热—结构耦合分析转换为应力—应变谱。以涡轮盘为例,典型做法是:先由 CFD 或经验公式给出轮缘、轮心等关键节点的温度—时间曲线;再建立有限元模型,赋予材料随温度变化的弹性模量、热膨胀系数、屈服强度和蠕变参数;然后施加离心力、气动载荷和位移约束,进行瞬态热—结构耦合求解,得到关键点的应力—时间历程。
在工程简化中,常用热应力估算式 σ_th ≈ E·α·ΔT/(1−ν) 判断热应力数量级,其中 ν 为泊松比。该式适用于完全约束、温度均匀变化的理想情况,实际构件需乘以约束系数。例如某高温合金 E = 180 GPa、α = 14×10⁻⁶ /℃、ν = 0.3,当 ΔT = 500 ℃ 时,σ_th ≈ 180×10³×14×10⁻⁶×500/0.7 ≈ 1800 MPa,已远超材料屈服强度,说明必须考虑塑性松弛和蠕变。这一估算结果属于教学示例,实际设计需以有限元分析为准。
三、真实工况谱的编制要点
编制载荷谱时,首先要明确任务剖面。以航空发动机为例,一个典型飞行循环包括地面起动、起飞、爬升、巡航、下降、着陆和地面停车。每个阶段对应不同的转速、燃气温度和持续时间。工程上通常将整个任务剖面离散为若干载荷步,每步给定温度、转速和持续时间,形成“载荷步—时间”序列表。
其次要处理多轴效应。涡轮盘轮缘部位同时承受径向离心应力、切向热应力和轴向气动弯矩,主应力方向随飞行阶段变化。此时不能简单用单轴应力幅值评估寿命,而应提取危险点的应力张量历程,计算等效应力幅值或采用临界平面法寻找最大损伤平面。
最后要进行循环计数与损伤累积。雨流计数法能把复杂的应力—时间历程分解为一系列完整循环,每个循环对应一个应力幅值和平均应力。再结合材料的高温疲劳数据、蠕变—疲劳交互作用模型,用线性累积损伤法则或非线性损伤模型计算总损伤。当总损伤 D ≥ 1 时,认为构件达到寿命极限。
四、四类典型工况的谱特征
航空发动机涡轮叶片的热循环以“快速升温—长时高温—快速降温”为特征,单次飞行循环的 ΔT 可达 600~1100 ℃,升温速率数百 ℃/s,主导损伤为低周热疲劳与蠕变交互。燃气轮机透平转子的启停相对平缓,ΔT 约 300~900 ℃,升降温速率几十 ℃/min,但启停次数多、保温时间长,氧化和热疲劳并重。热交换器管板承受频繁启停和介质温度波动,ΔT 约 100~400 ℃,热疲劳与腐蚀共同作用。电子封装焊点在环境温度循环和功率循环下工作,ΔT 约 −55~150 ℃,温变速率快、保温短,焊点因热膨胀失配产生剪切塑性应变,主导损伤为热疲劳与蠕变。
理解这些差异,才能在试验谱编制和寿命评估中选择正确的加速模型和损伤准则。例如电子封装常用 Coffin–Manson 模型关联塑性应变幅与寿命,而涡轮盘则需采用应变范围划分法或蠕变—疲劳交互损伤模型。载荷谱的每一个参数,最终都会影响寿命预测的保守程度和工程决策的可靠性。