3. ANSYS Workbench热分析模块入门

在ANSYS Workbench中,热分析主要通过两个独立系统完成:Steady-State Thermal(稳态热)与Transient Thermal(瞬态热)。二者共享同一套几何、材料与网格体系,但在物理假设、边界条件类型、求解器设置和结果时间属性上存在本质差异。理解这些差异,是正确选择分析类型并合理解读结果的前提。

一、两个系统的界面结构与数据连接

在Workbench项目原理图中,从Toolbox拖入Steady-State Thermal或Transient Thermal后,系统会生成标准的A/B/C…单元格序列:Engineering Data(材料)、Geometry(几何)、Model(模型)、Setup(求解设置)、Solution(结果)、Results(后处理)。其中Model单元格是共享入口,双击后进入Mechanical界面,热分析的全部边界条件、求解器设置和结果提取都在此完成。

两个系统最直观的区别在于时间属性。稳态热求解的是不随时间变化的温度场,其控制方程为∇·(k∇T)+q=0,结果中不存在“时间”维度;瞬态热求解的是温度随时间的演化过程,控制方程包含密度与比热项ρc∂T/∂t,因此必须定义初始温度场和总求解时间。若将稳态热结果作为瞬态分析的初始条件,可在项目原理图中用鼠标拖拽建立数据连接,这是工程中常见的“先稳态后瞬态”做法。

Workbench 热分析系统结构与数据流 Steady-State Thermal(稳态热) Engineering Data Geometry Model Setup Solution 控制方程:∇·(k∇T) + q = 0 结果:稳态温度场,无时间维度 典型边界:温度、对流、热流、辐射 Transient Thermal(瞬态热) Engineering Data Geometry Model Setup Solution 控制方程:ρc∂T/∂t = ∇·(k∇T) + q 结果:温度随时间变化历程 必须定义:初始温度、总时间、时间步 稳态结果 作为初始条件 热分析边界条件类型与适用性 边界条件类型 物理含义 稳态适用 瞬态适用 Temperature 指定温度值 ✓ 常用 ✓ 可随时间变化 Convection 对流换热系数h与流体温度 ✓ 常用 ✓ h可随时间变化 Heat Flow 热流率(W) ✓ 常用 ✓ 可随时间变化 Radiation 辐射换热,需定义发射率 ✓ 高温工况 ✓ 高温瞬态

二、边界条件类型与施加方式

在Mechanical界面中,右键单击Steady-State Thermal或Transient Thermal分支,可在Insert菜单中看到热边界条件列表。常用的四类边界条件如下:

  • Temperature(温度):在面、边或点上指定已知温度值。稳态分析中为定值;瞬态分析中可定义为随时间变化的表格数据或函数表达式。例如,将某表面温度设为25℃,模拟恒温环境。
  • Convection(对流):需要输入对流换热系数h(W/(m²·K))和流体参考温度T∞。h可设为常数,也可随温度或时间变化。例如,空气自然对流h≈5~25 W/(m²·K),强制风冷h≈25~250 W/(m²·K)。
  • Heat Flow(热流率):在面或体上施加总热流率(单位W),正号表示流入,负号表示流出。适用于已知发热功率的电子元件或加热器。
  • Radiation(辐射):需定义发射率ε和环境温度。当表面温度超过200℃时,辐射换热占比显著,不可忽略。

此外,还有Heat Flux(热流密度,W/m²)、Internal Heat Generation(内部生热,W/m³)等类型。施加边界条件时,选择对象可以是几何面、边或顶点,选择后会在模型树中生成对应分支,并可在Details面板中修改数值和作用范围。

三、求解设置的关键差异

稳态热分析的求解设置相对简单,主要关注非线性控制(如材料属性随温度变化时)和求解器类型。在Analysis Settings中,可设置是否开启大变形(热分析通常关闭)、是否使用自动时间步、以及非线性收敛准则。

瞬态热分析的求解设置则复杂得多,核心参数包括:

  • Number of Steps(载荷步数):定义载荷随时间变化的阶段数。例如,0~10s加热,10~60s冷却,可设为2个载荷步。
  • Step End Time(步结束时间):每个载荷步的终止时间,单位秒。
  • Time Integration(时间积分):推荐开启,使用隐式积分格式,保证数值稳定性。
  • Auto Time Stepping(自动时间步):设为On后,求解器根据温度变化率自动调整时间步长。初始时间步可设为总时间的1/100~1/50,最小时间步建议不小于1e-6s。
  • Initial Temperature(初始温度):在Initial Conditions中定义,可设为均匀温度或从稳态结果导入。
对于大多数工程瞬态热分析,建议先运行稳态热分析获得初始温度场,再将其作为瞬态分析的初始条件。这样可避免因初始温度假设不合理导致的计算振荡或物理失真。

四、结果查看与判读

求解完成后,右键单击Solution插入结果对象。稳态热分析常用结果包括:Temperature(温度场云图)、Total Heat Flux(总热流密度)、Directional Heat Flux(方向热流密度)。瞬态热分析除上述结果外,还可插入Temperature Probe(温度探针)和Chart(图表),用于提取特定节点温度随时间的变化曲线。

结果查看时需注意以下几点:

  • 温度云图的色标范围应手动设置,避免自动缩放掩盖局部高温区。例如,若关注点温度在80~120℃,可将色标范围设为20~150℃。
  • 瞬态分析中,可通过时间滑块或动画播放查看温度场演化过程。重点关注温度梯度最大区域和达到热平衡的时间点。
  • 热流密度结果可帮助识别散热路径和热阻集中位置。若某处热流密度异常高,可能意味着该区域需要加强散热或存在几何尖角。
  • 对于瞬态分析,建议提取关键点的温度-时间曲线,判断是否达到稳态。若在总时间结束时温度变化率小于1℃/s,可认为接近稳态。

五、一个完整的教学示例

以一块100mm×100mm×10mm的铝板为例,材料参数:密度2700kg/m³,比热900J/(kg·K),导热系数237W/(m·K)。初始温度25℃,底面施加100℃恒温,顶面与25℃空气对流,h=15W/(m²·K),四周绝热。

稳态分析:底面100℃,顶面对流散热,求解后顶面温度约为78℃(示例估算值,实际结果受网格和边界影响)。瞬态分析:总时间600s,初始时间步5s,自动时间步开启。求解后提取顶面中心点温度曲线,可见温度从25℃逐渐上升,约在400s后趋于稳定,最终接近稳态值。

通过对比稳态与瞬态结果,可直观理解热惯性对温度响应的影响:瞬态分析中,铝板需要一定时间才能达到热平衡,而稳态分析直接给出最终平衡状态。这一差异在电子散热、模具加热等工程场景中至关重要。

瞬态热分析中,若时间步长过大,可能遗漏温度快速变化阶段;若时间步长过小,计算时间显著增加。建议先以较粗时间步试算,观察温度曲线拐点,再在拐点附近细化时间步。