ABAQUS热应力仿真：双金属片挠度与热障涂层残余应力分析

项目过程

本项目包含两个独立但方法同源的仿真任务：一是双金属温控片在不同温差下的自由端挠度预测，用于某温控开关的行程设计；二是航空发动机热障涂层（TBC）在冷却过程中沿深度方向的残余应力分布，用于评估涂层剥落风险。团队均采用ABAQUS/Standard建立热-结构顺序耦合模型（先求解温度场，再将节点温度作为预定义场导入结构分析）。

双金属片模型：长度L=50 mm、厚度t=1 mm，由两层金属胶合而成，对比两组：Cu/steel（铜+钢，热膨胀系数差Δα=5×10⁻⁶ K⁻¹）和brass/steel（黄铜+钢，Δα=9×10⁻⁶ K⁻¹）。边界为左端固定、右端自由，温差ΔT从0到200 K线性施加。网格在厚度方向每层层8单元，长度方向加密到ΔT梯度大的自由端附近。为验证，团队还用经典公式δ≈3L²ΔαΔT/(4t)做解析对照。

热障涂层模型：基底为镍基高温合金（热膨胀系数α_base≈16×10⁻⁶ K⁻¹），上为粘结层（MCrAlY，α≈14×10⁻⁶），最外覆YSZ陶瓷层（α_YSZ≈11×10⁻⁶ K⁻¹）。冷却温差约600 K（从服役温度降到室温）。输出沿涂层深度方向的残余应力分布，特别关注YSZ层与粘结层界面处应力突变和YSZ层整体应力状态。材料属性均设为温度相关，以反映高温下模量下降。

案例分析

案例图分左右两幅。左图是双金属片自由端挠度随温差变化：Cu/steel组合（Δα=5×10⁻⁶ K⁻¹，蓝）在ΔT=100 K时挠度约0.94 mm；brass/steel组合（Δα=9×10⁻⁶ K⁻¹，橙）挠度更大，呈线性上升趋势，ΔT=200 K时达约3.5 mm。右图是热障涂层残余应力沿深度分布：YSZ陶瓷层呈现压应力σ≈-144 MPa，粘结层附近应力梯度陡峭，基底处为拉应力σ≈+14 MPa。

双金属片挠度可用近似公式δ≈3L²ΔαΔT/(4t)验证。代入Cu/steel：ΔT=100 K时δ≈3×50²×5×10⁻⁶×100/(4×1)=0.94 mm，与图中一致；brass/steel因Δα大一倍，同温差下挠度约为1.8倍，符合线性关系。这解释了为何温控开关偏好黄铜/钢组合——同样温升获得更大动作行程，触点行程裕量更足。右图残余应力反映陶瓷层与金属基底热膨胀系数失配：YSZ热膨胀系数低，冷却时收缩小，被金属基底”拉住”而产生压应力；金属基底则受拉。YSZ层-144 MPa的压应力是TBC的典型状态（压应力有利于抑制表面裂纹张开），但粘结层附近陡峭的拉应力梯度是界面开裂与剥落的诱因。若服役温度循环加剧，该界面易成为失效起点，需通过梯度过渡层缓解。

总结教训

第一，双金属片建模若把胶合层简化为单一界面、忽略胶层刚度，挠度会偏高约8%；团队加入等效胶层厚度（0.02 mm）修正后结果才可靠，说明薄界面层不可随意删。第二，热障涂层必须考虑温度相关材料属性，高温下YSZ弹性模量从室温约200 GPa降到服役温度约120 GPa，残余应力相应降低；用常温属性会高估-144 MPa。第三，界面处网格要足够密（至少5层），否则应力梯度被抹平，-144 MPa峰值得不到，且可能漏判局部拉应力极值。第四，冷却速率影响瞬态热应力峰值，但残余终态应力主要由失配决定，速率影响有限，工程评定可用稳态残余分布。第五，双金属片若考虑塑性（铜层屈服），大ΔT下挠度会偏离线性，本项目在弹性范围内验证，超范围需启用塑性本构。

项目经验

1. 双金属片挠度与Δα成正比：δ≈3L²ΔαΔT/(4t)，选材优先大Δα组合以获取更大动作行程，温控开关设计基础。
2. TBC残余应力是涂层剥落主因：YSZ压应力(-144 MPa)与粘结层拉应力梯度需联合评估，梯度过大引发界面开裂。
3. 温度相关材料属性不可省：E、α、屈服强度在高温下变化显著，尤其陶瓷层模量下降影响残余应力预测准确性。
4. 胶层/界面不可随意简化：薄界面层刚度会显著改变挠度与应力分布，建模需保留或做等效处理。
5. 梯度过渡层可降界面应力：TBC加成分梯度粘结层能缓解拉应力梯度，是延长涂层寿命的常用手段。
6. 失配应变是热应力根源：任何双材料/涂层系统的残余应力都来自α差，选材时把Δα作为第一性指标可预判应力水平，从源头规避开裂。

更多ABAQUS热应力与多材料耦合仿真经验，可参考科研学术网ABAQUS仿真栏目。