1. 热应力与热结构耦合分析概述
在工程结构中,温度变化几乎无处不在:发动机缸体在冷启动与全负荷之间经历上百摄氏度的温差,电子封装在回流焊后从两百多摄氏度冷却到室温,管道在输送高温介质时被支架约束而无法自由伸长。这些场景的共同点是——温度场并不"安分"地只停留在热学范畴,它会通过材料的热膨胀转化为变形,而变形一旦被约束或分布不均,就会在结构内部生成应力。这就是热应力,也是热结构耦合分析要解决的核心问题。
热应力产生的机理
热应力的根源可以归结为一条基本关系:材料受热要膨胀,遇冷要收缩,而膨胀或收缩受到阻碍时才会产生应力。用一维情形说明最清楚。一根长度为 L、热膨胀系数为 α 的杆,温度均匀升高 ΔT,若两端完全自由,其自由伸长量为 ΔL = α·L·ΔT,此时杆内应力为零。若两端被刚性固定,伸长被完全抑制,则杆内产生压应力 σ = −E·α·ΔT,其中 E 为弹性模量。若约束是部分弹性约束,实际应力介于零与全约束值之间,取决于约束刚度与杆件轴向刚度的比值。
由此可以提炼出热应力产生的三个必要条件:存在温度变化、材料具有非零热膨胀系数、变形受到约束。约束的来源有三类,理解它们的区别对后续建模至关重要。第一类是外部约束,如固定支座、螺栓连接、与刚性件的接触;第二类是内部约束,即结构内部温度分布不均匀,高温区想多膨胀、低温区不让,彼此牵制,即使整体没有任何外部支撑也会产生应力;第三类是材料约束,多材料结构中各组分热膨胀系数不同,同一温升下变形量不一致,界面处必然产生应力,典型如双金属片、陶瓷封装与金属基板、复合材料层合板。工程中这三类约束往往同时存在,因此热应力问题通常不能靠单一公式估算,需要数值方法求解。
还有一点容易被忽略:热应力与温度场的时间历程密切相关。急剧升温时,表面先热、内部仍冷,表面膨胀受内部牵制而产生压应力,内部则受拉;缓慢升温时温度趋于均匀,应力水平显著降低。这说明热应力分析往往不能只取一个稳态温度,还要考虑瞬态过程。
热结构耦合的两条技术路线
按温度场与位移场是否在同一方程组中求解,热结构耦合分为顺序耦合与直接耦合。
顺序耦合是先做热分析得到节点温度,再把温度作为"体载荷"施加到结构分析中计算应力与变形。两个物理场各自使用独立的单元类型、材料属性和求解器,通过温度结果的传递建立联系。它的优点是计算规模小、每步都可独立检查、便于排查问题,是绝大多数工程热应力问题的首选。顺序耦合又分单向与双向:若变形对温度场的影响可以忽略(大多数情况如此),单向传递即可;若变形显著改变传热路径,例如接触间隙随热变形张开或闭合、流体通道截面随热膨胀变化,则需要双向迭代,把结构变形反馈回热模型重新计算温度。
直接耦合则使用同时具有温度自由度和位移自由度的耦合单元,在一个矩阵中同时求解温度与位移,天然包含双向效应。它适合强耦合场景,如压电热耦合、热-机-电多场问题、接触热阻随压力显著变化的情况。代价是矩阵规模大、非线性强、收敛困难、调试成本高。
选择原则可以概括为:先判断变形是否反过来显著影响温度场。若影响小,用顺序耦合单向;若影响明显但可以分步迭代,用顺序耦合双向;若两场强耦合且难以分离,才用直接耦合。对课程所覆盖的绝大多数工程问题,顺序耦合已经足够。
Ansys Mechanical 在热分析中的定位
Ansys Mechanical 是 Workbench 环境下的结构与热分析模块,它本身并不替代 Fluent、CFX 这类 CFD 软件去求解复杂流场,而是承担"固体域传热 + 结构响应"这一段的求解。它的热分析能力包括稳态热、瞬态热、热-结构顺序耦合,以及通过耦合场单元或系统耦合实现的直接耦合。温度载荷可以来自自身热求解,也可以从外部读入(如从 CFD 映射过来的对流换热系数与流体温度、从电路仿真导入的功耗分布)。
在 Workbench 项目流程中,典型做法是建立"Steady-State Thermal 或 Transient Thermal → Static Structural 或 Transient Structural"的链接,把热求解的 Temperature 结果直接拖入结构分析作为载荷。材料数据方面,热分析需要导热系数 k、密度 ρ、比热容 c;结构分析需要弹性模量 E、泊松比 ν、热膨胀系数 α,若考虑高温蠕变或塑性还需补充相应本构。需要特别提醒的是,α 必须与参考温度(Reference Temperature)配套使用,参考温度通常取装配或应力自由状态的温度,取错会直接导致整体刚体膨胀被误算成应力。
典型工程应用场景
热应力分析在工程中覆盖面很广,常见的有:电子封装与 PCB 在回流焊及功率循环下的焊点应力与疲劳寿命评估;涡轮叶片、燃烧室在高温燃气与冷却气流作用下的热梯度应力;压力容器与管道在启停、暖管过程中的热冲击与热棘轮;制动盘在反复制动中的热裂纹倾向;模具在注塑或压铸循环中的热疲劳;光学镜筒与精密结构在温度波动下的热变形导致的光轴偏移。这些场景的共同特征是:温度场往往不均匀或随时间剧烈变化,而结构又受到装配、接触或材料差异的约束。
课程学习路线
本课程按"机理—工具—实操—进阶"的顺序展开。先建立热应力与耦合方式的概念框架(本章),再依次掌握稳态热分析、瞬态热分析、热载荷与边界条件的施加、材料热物性设置、顺序耦合工作流的搭建与结果传递、热应力结果的判读与常见错误排查,随后进入接触热阻、热疲劳、多载荷步与参数化等进阶主题,最后通过若干完整工程案例把流程串起来。学习过程中建议始终围绕一条主线:温度场从哪来、约束在哪里、应力由什么主导、结果是否合理。