如果你手头有个需要评估温度场的零件，或者刚打开ANSYS Workbench不知道从哪开始，E3-1这个案例就是个很典型的起步练习。它名叫“稳态热力学分析”，实际做起来并不复杂：建好模型、赋予材料、加发热和散热条件，让软件算出“热平衡后的温度分布”。它解决的问题很具体，比如一个铝合金壳体在固定热源下，最终哪个位置最烫、热量主要通过哪条路径散掉、散热片面积是否需要加大。刚入门热仿真的工程师，或者做结构设计顺带想验证散热方案的人，都能从中获得一套完整操作思路。
机械工程里常说的热分析，严格讲应该叫“传热分析”而不是“热力学分析”，但教程里沿用这个叫法太普遍了。E3-1的价值在于用一个极简模型，帮你把稳态传热的边界条件识别、材料参数设置、求解和后处理完整过一遍。我在带新人时经常先用这类金属块案例把流程跑通，再上完整装配体。大多数仿真结果误差不是软件不行，而是边界条件给错了，这个案例恰好能帮你建立正确的操作习惯。
1. 先吃透稳态热分析的基本逻辑
1.1 “稳态”到底在求什么
热分析和结构分析最大的差异在于：你关心的不再是“哪里应力最大”，而是“热量怎么流、温度怎么分布”。稳态和瞬态的区别，可以用开暖气来类比。刚开暖气的半小时，房间温度一直在上升，这是瞬态；等温度不再变化，室内与外界换热达到平衡，就进入了稳态。稳态分析不计算“需要多久达到平衡”，只计算“达到平衡后温度长什么样”。
数学上，稳态传热控制方程中时间项被置零，系统能量累积为零，进入系统的热量等于离开系统的热量。固体内部的导热方程退化成稳态拉普拉斯型方程，边界条件给的是温度、热流量、对流换热或辐射。对于电子设备外壳这类场景，芯片恒定功率发热，外壳表面通过对流把热量交给空气，一段时间后温度不再变化，这时工程上关心的最高温度、平均温度、热流路径，都是稳态分析能直接给出的。
一个容易出现的小误区是：有些人觉得“稳态”就是“均匀温度”。实际上稳态只表示温度场不随时间变化，不同位置照样可能存在几十度的温差。比如一根导热金属棒，一端接高温热源，另一端被风扇冷却，稳定后热端和冷端之间依然有明显温度梯度，热量持续从热端流向冷端。
1.2 E3-1案例想让你练什么
E3-1在各类教程里的模型并不花哨：一个带凸台或通孔的矩形金属块，某一表面或孔内壁施加固定温度，另外几个表面施加对流换热，也可能在某小块面积上施加热流密度。难度刻意控制在简单几何范围，为的是让你专注三件事：正确识别热边界、正确设置材料导热参数、正确解读温度云图和热通量图。
很多初学者一上来就研究网格要多细、求解精度要多高，反而忽略了边界条件是否反映真实物理。对同一个模型，把换热面积搞错，或者把对流系数单位弄混，结果差个几倍很常见。我在实际跑E3-1时，习惯先做手算估计：总功率多少瓦，自然对流散热能力大概多少瓦每开，预计温升多少开。心里有了数量级再上软件，一眼就能看出结果有没有离谱。
这个先手算后仿真、先粗估后细算的习惯，是E3-1最值得带走的东西。后面做复杂装配体，边界条件多到几十个，如果没有数量级概念兜底，算出一个温度云图根本不知道是对是错。
2. 从零搭建Workbench稳态热分析流程
2.1 分析系统选择与几何准备
打开ANSYS Workbench后，左侧工具箱里能看到一个叫“Steady-State Thermal”的系统。鼠标拖到项目概览表上，它会自动生成一个分析流程：工程数据、几何模型、模型、设置、求解、结果，每个环节都能展开。也有人习惯先建“Geometry”再手动关联热分析模块，但拖一个完整系统更方便，后续要接着做热应力分析时，可以直接右键复制结果链接到静力学模块，省去重新设置边界条件的时间。
几何来源一般有两种。简单教学模型可以在DesignModeler或SpaceClaim里直接建，导入后建议统一成米或毫米，并提前设置好单位。我踩过的坑是：默认显示单位毫米，但画草图时心里想着寸，一小块尺寸差很多，后面网格和换热面积全错。建议进入Workbench项目页面后，先确认单位制，再开始几何操作。
热分析对几何模型有两条额外要求：一是要保留真正的接触面，不要随手合并掉不同材料的交界面；二是尽量简化不必要的圆角和小孔。这些细节对刚度分析可能有影响，对温度场一般影响不大，但会拖慢网格生成和求解速度。对E3-1这种小型模型当然不用担心性能，但从一开始养成“按分析目的做几何清理”的习惯很重要。
2.2 材料参数设置的关键细节
进入Engineering Data后，默认材料一般有Structure Steel等。稳态热分析中真正参与计算的材料参数只有一个：导热系数，单位W/(m·K)。密度和比热容在纯稳态导热里不出现，因为温度不随时间变化，不需要存储热量。但如果后面要接着分析热应力，弹性模量、线膨胀系数和泊松比也要补齐，所以很多项目中干脆一开始就把材料参数配全，避免后续补设。
导热系数从哪来？最可靠是材料商手册或实测数据，其次是常用工程手册。以铝合金为例，6061-T6的导热系数大约为167 W/(m·K)，纯铜约385 W/(m·K)，结构钢约45 W/(m·K)。需要注意的是，导热系数随温度会变化，高温工况下使用常温数据会带来明显误差。Workbench支持定义随温度变化的分段线性表或多项式，E3-1入门阶段用常数没问题，但真实项目如果温度跨度超过两三百度，强烈建议配温度相关导热系数。
还有一个细节：软件可能会自动补全密度比热等参数，但导热系数缺失时，材料会显示黄色警告，计算要么失败要么用默认值。给所有实体赋予材料这一步，在Mechanical界面里通过选择实体再在材料面板里指定。我见过不少“求解成功但温度全是-273℃”的案例，最后发现就是材料没赋上或导热系数被误填成0。
2.3 网格划分与接触设置
热分析的网格要求和结构分析类似，但也有自己的特点。传热问题是扩散型方程，相比应力集中问题对网格不是极度苛刻，但热通量梯度和温度梯度大的区域需要加密，比如高功率热源附近、截面突变的颈部、散热齿顶端。E3-1模型通常用四面体网格就能拿到不错结果，五六个毫米的全局尺寸加局部细化，总计几万个单元，完全够用。
关键提醒是：网格不是越细越好。网格加密到一定程度，温度结果变化很小，而无脑加密只会增加求解时间。正确的做法是跑一次后，把网格尺寸减半再跑一次，对比关键点温度，差值在1%以内就说明网格基本收敛。这就是网格无关性检查，E3-1阶段应该主动做一遍，比背十遍教程都有用。
接触设置更要注意。如果模型是多零件装配体，Mechanical默认的接触类型是Bonded，它在热分析中表示理想无热阻接触，热量能连续穿过接触面。实际上两个固体表面贴合时存在界面热阻，温度在分界面会有一个跳变。如果模型里有导热垫、钎焊层或仅仅靠螺栓压紧的搭接面，就需要在接触属性里设置Thermal Contact Resistance，单位是m²·K/W。E3-1如果是一个单零件模型，这部分可以暂时跳过，但要记得有这么一回事。
3. 边界条件施加：决定结果正确性的第一道关
3.1 温度、对流、热流三类边界怎么选
稳态热分析里最常见的边界条件有四类：固定温度、对流换热、热流密度、内部热生成。固定温度对应“该表面的温度已知且保持不变”，比如设备外壳接触恒温夹具；对流对应“表面与流体之间存在换热”，需要给定换热系数和流体参考温度；热流密度对应“单位面积上流入固定热量”，比如贴在表面的加热片；内部热生成对应“体积内均匀发热”，比如电阻发热体或芯片的等效发热体。
选择原则其实很简单：你能确定的是温度就用Temperature，能确定的是热量流通量就用Heat Flux，能确定的是流体冷却条件就用Convection。千万不要在一个表面上又加固定温度、又加对流，这样会过约束，软件会报错或某一项被忽略。E3-1练习的一个目的就是让你体会这些边界在数学上都是“边界条件”，但物理意义完全不同。
3.2 对流换热系数的取值经验和单位换算
对流换热系数是E3-1里最容易出错的地方。常见量级是：自然对流空气约5到25 W/(m²·K)，强制风冷约25到250 W/(m²·K)，水冷约500到10000 W/(m²·K)。如果你把自然对流的系数填成100，结果会明显偏乐观，表面温度会比实际低一大截。这个系数和流体速度、表面粗糙度、形状方向都有关系，手册值只能用来做初步估算。
单位换算也得留意：几何用毫米时面积是平方毫米，热流密度如果按焦每平方毫米每秒输入，数字会非常小，容易看错。更安全的方法是把模型单位统一成米，或者先把热流密度换算成W/m²。固定温度边界用摄氏度还是开尔文，在Workbench里都可以，因为导热和对流只关心温差，但涉及到辐射计算时必须用绝对温度，否则结果会错得离谱。
3.3 辐射边界什么时候不要忽略
E3-1案例通常只涉及传导和对流，辐射默认忽略，因为计算量小、概念简单。但真实世界中，一个高温外壳在空气中自然冷却，辐射散热有时占总散热量的三成以上。判断要不要加辐射，有个经验：表面温度和环境温度差超过50℃时，辐射基本不能忽略；温差小或表面是低发射率抛光金属时，可以先不带。
Workbench里辐射边界有Surface-to-Surface和Surface-to-Ambient两种设置方式，需要给出发射率，常见值为0.9左右（氧化铝表面、黑色阳极氧化）或0.1左右（抛光铝）。注意辐射里所有温度都必须用开尔文，因为辐射定律是绝对零度以上四次方关系，用摄氏温度会产生错误。E3-1如果想做一个加辐射的进阶练习，模型不变、只增加一个辐射边界，正好能直观看到温度分布的差别。
4. 求解与后处理：别只看云图好看
4.1 求解设置与收敛判断
稳态传热如果线性材料参数恒定、边界条件不随温度变化，其实一次迭代就出结果，不存在“迭代不收敛”的问题。当一个模型出现“不收敛”时，多半是引入了非线性因素：温度相关导热系数、辐射边界、接触热阻随温度或压力变化。E3-1入门如果只算导热和对流，通常Very Stable。但为了以后能处理更复杂的项目，还是要知道非线性要打开Newton-Raphson残差控制，让求解器迭代到热通量残差小于设定容差。
在求解前可以插入一个温度探针，指定某个关键点，然后在求解过程中监视它的收敛曲线。如果该点温度在迭代中反复震荡，说明模型里可能存在过约束或辐射项过强，需要回查边界。稳态求解结束后，Mechanical会在结果树里给出能量平衡相关信息，把所有边界上的总热流量加起来，如果净热量远不为零，说明边界设置不一致，结果需要重新检查。
4.2 温度场和后处理量如何解读
结果后处理常用的云图有Temperature、Total Heat Flux、Directional Heat Flux。Temperature云图看最高温度分布、热点位置；Total Heat Flux云图显示热流密度大小，箭头能指出热量从哪个区域流到哪个区域；Directional Heat Flux能单独查看沿某个方向的热流分量，用来判断翅片垂直方向导热是否充分。
温度云图最重要的不是颜色好不好看，而是数值范围是否符合物理常识。比如一个自然对流条件下10W热源的铝块，最高温升只有零点几度，基本不可能明显异常。还可以插入探针，查看关键节点的具体温度值。导出温度结果列表也很简单：选中结果，右键Export，生成CSV文件，方便做后续数据处理或写报告。
需要注意：热通量云图在零件分界面处出现明显跳变，不一定是错误，可能是两种材料导热系数不同，热量密度在界面处连续但梯度不同。但如果同一材料内部热通量出现锯齿状分布，就是网格太粗或接触设置有问题，这时候先加密网格再看。
4.3 用理论解做结果验证
仿真结果必须经过验证才能用于工程判断，最简单的验证是用一维平壁导热理论解。比如取一块截面积为A、长度为L的金属杆，两端温差为ΔT，导热系数为k，那么热流量q=kAΔT/L。在Workbench里建同样模型，一端设高温，另一端设低温，让侧面绝热，算完得到的热通量应该和理论值几乎一致。
我实际用6061铝做过一个小练习：尺寸100mm长、20mm见方，两端温差50K，理论热流量约为167×0.0004×50/0.1=33.4W。仿真结果通常在33W左右，误差在1%以内。如果误差明显偏大，要查材料导热系数是否设错、接触界面是否有隐性热阻、单位是否一致。这个练习花不了十分钟，却能帮你确认整个仿真链路是可信的。
5. 我在实际项目中踩过的坑与排查清单
5.1 单位不统一：最常见也最坑
我做过的模拟项目X里，有一次几何模型标注为毫米，热源功率数据来自另一边的电子表格，用的是瓦特，但热流密度如果直接按W/mm²输入，数字会变成百万级别，温度结果直接爆表。解决方案是先统一到国际单位制，或者明确知道自己在用什么单位。Workbench的项目单位可以统一设，但几何内部的参数、材料数据和控制仿真变量往往各有各的来处，项目越复杂越容易混。
温度单位也是一个高频出错点。导热系数、对流换热系数跟K还是℃无关，但固定温度如果填成绝对温度而环境温度用摄氏温度，温差自然不对。辐射中的温度必须用K。建议在模型树里给每个边界命名时直接带上温度单位和换算关系，比如“入口温度60C”，减少后期混淆。
5.2 接触热阻、网格假热与结果跳变
接触热阻是最减弱“仿真很准”信心的地方。两个零件用螺栓连接，热量通过接触面时存在压紧应力、粗糙度、表面氧化膜的影响。理论上的界面热阻很难精确给定，实测数据也不易查到。E3-1不涉及多体接触，但如果你按现实装配体做，默认Bonded往往会让温度偏低，因为忽略了装配缝隙的实际热阻。
网格对热分析的影响虽然不像应力分析那样敏感，但碰到的怪异温度条纹多半还是网格问题。比如热流路径上出现长条状冷区，可能是有单元被拉得很长，几何质量下降。检查时看Element Quality和Skewness，优先处理小于0.1的单元。不要为了美观强行把全局网格设成0.1mm，那样求解时间拉满，结果却不一定有本质变化。
5.3 常见报错信息与对策清单
| 报错或异常现象 | 常见原因 | 对策 |
|---|---|---|
| 求解器提示没有求解位置 | 没设置任何类别，如没有温度边界或热生成 | 检查边界和内部热生成是否已施加 |
| 温度出现负值或绝对零度附近 | 材料导热系数为0，或单位制错误 | 核对材料参数和单位换算 |
| 非线性不收敛 | 辐射项过强或接触热阻反复振荡 | 降低辐射松弛因子，检查接触设置 |
| 结果云图全是同一种颜色 | 载荷太小导致温差在显示精度内 | 调整图例范围或检查边界载荷量级 |
| 热通量在体内出现锯齿状 | 网格过粗 | 在热源附近局部加密 |
| 导入模型后少了一个体 | 几何修复失败或零件被抑制 | 在DM或SpaceClaim检查体的显示和抑制状态 |
这张表我整理自多次实操记录，每次碰到类似问题都可以按这个顺序排查。尤其是“材料没赋值”这种低级错误，常常被云图上一大片异常颜色掩盖，先查单位再查材料，基本能解决一半问题。
6. 一些提升效率的建模习惯和扩展思路
6.1 命名选择和参数化：改尺寸后自动更新
在E3-1这个级别，命名选择可能看不出来有多大优势。但项目一旦有几十个表面和热源，每次选面都会变成体力活。我习惯在几何创建设阶段就给关键面命名，比如HotEnd、CoolSide、HeatSourceArea，这样到了Mechanical里选边界时直接下拉选择，不需要在三维视图里满屏点选。后续如果改几何或改边界，命名选择自动继承，也减少漏选错选。
参数化则是Workbench的另一大杀招。在DesignModeler或SpaceClaim里把关键尺寸设为参数，比如散热齿厚度、基板长度，回Workbench创建一个参数组，可以批量扫描不同尺寸下的最高温度。对于散热片选型这种问题，相当于不用反复手动改模型和重新设边界。建议在E3-1阶段就试着把一个尺寸参数化，跑完一次后修改参数更新几何，体会一下参数联动带来的便利。
6.2 仿真结果与热阻网络估算互验
传热仿真不是“算完就完事”，工程上常用热阻网络对结果做快速复核。把一个散热路径拆成串联热阻：芯片到外壳的扩散热阻、外壳到环境的对流热阻，再加上接触热阻。用总温差除以总热阻，就能估算出稳态温升。这个估算结果和Workbench结果差距在20%以内，说明模型合理；差距很大就要检查边界条件或者热阻拆分是否遗漏了并联路径。
我在做某跨平台系统的散热验证时，先用热阻公式算出预计最高温度75℃左右，Workbench算完是73.6℃，两者对得上，信心就足很多。如果仿真结果和手算差距大，通常不是手算错就是仿真边界漏了什么东西，值得多花时间查，而不是直接改报告数字。
6.3 后续还能怎么扩展
E3-1做完，最简单的扩展方向是把稳态温度场当作载荷，继续做热应力分析。在Workbench里复制一个Static Structural系统，把Thermal Analysis的温度结果作为体载荷或边界条件拖过去，就能看到热膨胀产生的应力，这是电子设备可靠性分析最常见的一条路。另一个方向是换成瞬态热分析，加一个时变功率曲线，观察温度过冲和热响应时间，适合评估芯片短时高负载场景。
无论怎么扩展，基础逻辑都是一样的：先把稳态传热看明白，知道温度梯度和热流路径，再谈结构强度和疲劳。E3-1虽然叫“热力学分析”，但它真正的教学目的是让你建立起对热边界条件、材料参数和求解验证的整体感觉。这个感觉一旦建立起来，后续遇到再复杂的热仿真模型，心里都不会发慌。
