电缆散热设计这事看起来简单，实操起来全是坑。上个月给某变电站做电缆沟热仿真，甲方拿着计算器咔咔按公式说肯定没问题，结果实测温度超了十几度。后来用COMSOL重新建模才发现，土壤热阻和邻近电缆的热耦合完全没考虑进去。今天就带大家实操一把电热耦合仿真，手把手教你怎么避坑。
先看模型搭建要点。COMSOL里直接选"焦耳热"+"固体传热"多物理场耦合模块，重点是要把电缆的三层结构（导体、绝缘层、铠装）拆分开建模。这里有个偷懒技巧：用参数化扫描批量生成不同直径的同轴结构（如图1）。别傻乎乎手动画圆环，用下面这段代码自动生成：
for i = 1:3
model.geom('geom1').feature().create(...
strcat('circ',num2str(i)), 'Circle');
model.geom('geom1').feature(...
strcat('circ',num2str(i))).set('r', 0.02*i);
model.geom('geom1').feature(...
strcat('circ',num2str(i))).set('pos', {'0' '0'});
end
这段脚本循环生成三个同心圆，半径按0.02m递增。注意绝缘材料的热导率各向异性特性得单独设置，比如XLPE材料在径向和轴向的导热系数能差3倍。设置材料属性时别直接输数值，建议用变量表达式，方便后期调参：
model.param.set('k_radial', '0.35 [W/(m*K)]'); //径向导热
model.param.set('k_axial', '0.12 [W/(m*K)]'); //轴向导热
网格划分最容易踩雷。见过有人用自由四面体网格，结果导体边缘温度场全是锯齿。建议在导体表面做边界层网格（Boundary Layer），特别是电流趋肤效应明显的区域。分享个网格控制参数设置模板：
model.mesh('mesh1').feature.create('bndLayer1', 'BoundaryLayer');
model.mesh('mesh1').feature('bndLayer1').set('thickness', '0.002');
model.mesh('mesh1').feature('bndLayer1').set('numsublayers', '3');
跑完仿真别急着看温度云图，先检查焦耳热源分布是否合理。遇到过导体横截面上热源分布不均匀的情况，最后发现是材料电导率单位设成了S/m²（正确应该是S/m）。后处理阶段推荐用截面线图对比不同工况，比如这段代码生成沿电缆径向的温度分布曲线：
model.result().export().create('plot1', 'Line');
model.result().export('plot1').set('data', 'dset1');
model.result().export('plot1').set('expr', 'T');
model.result().export('plot1').set('line', {'0' '0' '0' '0.05'});
实际项目中最头疼的是土壤参数设置。某次仿真结果比实测低20℃，排查发现土壤含水率参数用的是夏季数据，实际冬季土壤干燥导致热阻增大。建议做参数化扫描，把土壤导热系数设为0.8~1.5 W/(m·K)范围动态分析。
最后给个实战案例：某地下管廊6回电缆并列敷设，仿真发现中间电缆温升比设计值高8℃。通过调整电缆间距从150mm增至200mm，同时将平面排列改为三角形排列，最终各回路温差控制在3℃以内。关键看温度场分布是否出现"热岛效应"，必要时得加装通风管道或采用导热回填材料。
记住，好的热仿真工程师都是被异常数据喂出来的。多对比实测数据，建立自己的材料参数库，比死磕仿真软件自带的材料库靠谱得多。下次碰到电缆载流量争议，直接甩出带环境耦合的仿真报告，比跟甲方扯皮热力学公式管用十倍。
