电缆中间接头的电 - 热 - 力多物理场耦合仿真之旅Comsol 6.3 实战干过电缆附件设计的朋友都有这种体会中间接头是整个输电回路里最容易“翻车”的环节。绝缘击穿、局部放电、过热老化十起电缆故障里中间接头能占一半以上。这不是制造端不认真而是中间接头天然就是个多物理场矛盾集中地——电压加在绝缘上电流在导体里发热发热又让橡胶和绝缘层膨胀膨胀又反过来改变界面压力界面压力直接决定电场分布是否均匀。这些因素纠缠在一块单看电场或者单看温度怎么看都是对的但放在一起就可能出问题。我这次用 COMSOL 6.3 对一个典型的中压电缆中间接头做了完整的电 - 热 - 力多物理场耦合仿真把电场、温度场、应力场放在同一个模型里同时解出来跑了稳态工况下的全耦合分析。整个过程下来最大的感受是模型本身不难搭真正的功夫全在“耦合关系怎么设”“参数怎么定”“结果怎么解读”这三件事上。这篇文章把我这次从建模到后处理的全过程、关键参数、踩过的坑都记录下来适合正在做电缆附件仿真、或者准备入手多物理场耦合分析的朋友参考。1. 为什么中间接头非要“电-热-力”一起算1.1 中间接头的结构特点和应力锥的作用中间接头的工作环境和本体电缆完全不同。本体的绝缘是均匀挤包成型的电场在绝缘层里基本是径向分布接头不一样——它要把两段电缆的导体接起来绝缘恢复要靠预制件套上去这就带来了两个先天问题一是导体连接处形状突变电场会集中二是预制绝缘和电缆本体绝缘的界面是贴合的贴合压力一旦松了层间就可能出现微观气隙气隙里的电场强度远高于固体绝缘局部放电就从这里开始。为了把电场理顺接头设计里必然有应力锥。应力锥本质上是一个半导电材料做成的锥形结构它把绝缘层表面原本很高的轴向电场逐步过渡到低场强区域让电场分布尽量平缓。从仿真角度看应力锥的尖端就是一个电场“防火墙”它的锥面角度、贴合长度、材料电导率直接决定了整个接头的绝缘裕度。1.2 单向物理场仿真的盲区如果只做静电场分析模型会画得很漂亮但工程上不够用。电场分布依赖的是界面几何而界面几何在真实运行中不是设计图纸上那个固定形状——电缆一载流导体发热铜的线性膨胀系数虽然不高但绝缘和半导电橡胶的热膨胀系数要大一两个数量级。橡胶受热往外顶应力锥和绝缘界面的压力就变了界面压力变小接触变差局放风险就上去了。这就是电 - 力耦合前提是先有温度场。反过来温度场也不是独立的一门课。绝缘层里的介质损耗、导体和压接区的焦耳热主要取决于电压和电流但导体电阻率随温度升高而增大温度越高发热越厉害这是一个正反馈。同时热传导路径上如果材料变形了、接触界面叠差变了热阻也会跟着变。所以温度场本身就和电场、力场绑定在一起。1.3 三条耦合路径是怎么串起来的实际工程里最关心的三条耦合路径一句话就能串起来电压和电流输入电能一部分变成热量热量让材料温度升高进而改变材料电导率、导热系数和热膨胀热膨胀产生机械应力机械应力又反过来改变接触状态和电场分布。用 COMSOL 的术语来说就是“电流接口 → 电磁热多物理场 → 固体传热接口”以及“固体传热 → 热膨胀多物理场 → 固体力学接口”两条链路。介质损耗放在绝缘域里作为热源应力场求解结果可以回填到接触界面条件里再参与电场求解。这三者循环迭代直到收敛。有些刚接触耦合仿真的朋友喜欢把三个物理场拆开单独算然后手动相互赋值迭代。这么做不是不行但接口之间的手动传递容易忘掉“温度依赖材料参数”这一步误差也会在迭代中积累。用 COMSOL 的多物理场节点把三个接口绑定在一个研究里求解器会自动处理交叉耦合项省心得多。2. 仿真前的思路梳理模型、参数、物理场接口怎么定2.1 几何简化的取舍仿真和论文不一样论文恨不得把所有细节都画出来仿真首先要考虑算得动。我这次建的是二维轴对称模型以接头轴向为对称轴取一半截面。10kV 等级的电气设备结构尺寸相对紧凑二维轴对称精度足够三维模型留给后续做弯头、水下敷设等特殊工况再考虑。几何上我保留了这些关键部件导体压接管、导体屏蔽层、XLPE 绝缘、绝缘屏蔽层、半导电应力锥、硅橡胶增强绝缘、外半导电层、铜丝屏蔽层、防水密封壳。有些细节故意省略了比如压接管的齿痕、密封胶填缝等因为它们的尺寸在毫米级以下对宏观电场和温度分布影响可以忽略却会在网格剖分时制造大量畸形单元。一个容易忽视的细节是几何中的尖角。应力锥末端和导体压接管端部都是天然的高场强集中区这些地方必须真实保留绝不能为了网格好画就倒个圆角。我看着网格剖分器在这些位置报错时也犹豫过要不要简化后来发现只要局部网格足够密尖角部分是可以算稳的。2.2 材料参数收集与单位陷阱材料参数是整个仿真里最脏最累的活也是最容易出错的地方。我这次用到的核心材料有四类铜导体、XLPE 绝缘、半导电硅橡胶、硅橡胶增强绝缘。参数我从厂家公开数据以及通用材料库中汇总整理成表供参考。材料相对介电常数介质损耗因数 tanδ导热系数 W/(m·K)热膨胀系数 1/K弹性模量 MPa备注铜导体--4001.7e-5110000电导率 5.998e7 S/m温度系数采用默认XLPE 绝缘2.30.0010.2852.0e-4250线性弹性近似实际可选粘弹性模型半导电硅橡胶视频率/配方0.005 左右0.352.5e-42力学参数按超弹性预处理硅橡胶增强绝缘3.20.0050.352.5e-42同上参数单位是一个经典坑。比如导热系数工程手册里经常给 W/(m·K)但有些供应商提供的是 W/(cm·K)一换算就差了 100 倍热膨胀系数有的给 10⁻⁶/K有的给 10⁻⁴/K差一个数量级直接导致应力结果完全失真。我这次把所有参数先在表格里统一换算好再填进“全局定义”中的材料参数表避免在物理场设置界面里零散修改带来的前后不一致。比热容和密度建议按体积单位填因为 COMSOL 内建的固体传热接口默认使用体积热容。质量热容和密度要同时正确否则瞬态温升曲线会明显偏慢或偏快。我第一次跑瞬态时温升稳定时间异常长查了半天发现是密度用了 g/cm³ 而比热容用了 J/(kg·K)两个单位体系混在一起体积热容算出来整整差了 1000 倍。2.3 物理场接口与耦合节点的选择COMSOL 6.3 里针对电场计算有两个常用接口静电接口和电流接口。工频 50Hz 的导体系统波长以千公里计电气尺寸远远小于波长静电场近似足够但中间接头有泄漏电流和介质损耗靠静电接口没法算发热。所以我选了“电流”接口它求解的是电流守恒方程既能得到电场分布又能算出焦耳热和介质损耗热。热场我用“固体传热”接口。力学场用“固体力学”接口材料模型先用各向同性线弹性初步评价应力分布足够如果后续要做蠕变、回弹等精细化分析再升级为超弹性模型。耦合节点在 6.3 的“多物理场”选项卡里直接选“电磁热”和“热膨胀”程序会自动把电流场的损耗密度作为传热源项把温度场作为力学载荷源项。我最初以为还需要手动添加一个“电流→热”的耦合实际上 6.3 这两个内置节点已经把主要耦合项覆盖了能少绕很多弯。3. 一步步把模型搭起来COMSOL 6.3 实操记录3.1 全局定义与几何构建打开新模型后第一步先设“选择空间维度”为二维轴对称并添加三个物理场接口电流、固体传热、固体力学。然后在“全局定义”里建好所有参数包括几何尺寸、载荷电流、环境温度、材料参数。把参数集中写在这里后面改方案时只需要动一处。几何绘制我用的是内建 CAD 建模功能把接头按轴向分段画矩形再合并成复合域。这里要注意二维轴对称模型里 X 轴代表径向、Y 轴代表轴向画图时别把方向看反了。我一开始画完检查几何怎么都对不上后来才发现把径向和轴向搞反了。这个错误隐蔽性强因为云图出来还是有形状的只是数值全乱。关键几何数值我按典型中等截面设置了以下参考值实际项目中请以具体产品图纸为准导体压接管外径40 mm导体内孔直径15 mm应力锥锥面长度约 60 mm锥面末端距屏蔽层开口约 10 mm硅橡胶增强绝缘层厚度约 8 mm整体模型轴向长度500 mm其中端部预留 50 mm 电缆本体几何建完后我用“合并”操作把所有相邻域组合成几个大域便于后面分配合适的材料和边界条件。管线类的模型我一般不追求把每个小垫圈、密封圈都画出来CAD 布尔运算生成的碎片域不仅难赋属性还会拖慢网格剖分速度。3.2 边界条件与载荷设置几何搞定后接下来把三个物理场的边界条件逐个设好顺序建议按“电→热→力”来每一步思路都清晰。电流接口里我在一端导体端面设置“电势”边界条件电压幅值取相电压峰值另一端端面设为“接地”。实际电缆运行时三相都带电但接头局部模型的电位参考点可以设在远端屏蔽层这样结果更接近工程实测。绝缘和半导电层的材料属性里要填相对介电常数和电导率绝缘电导率取极小值模拟纯容性电流介质损耗靠 tanδ 折算的等效电导率来表达。固体传热接口里导体端面设为温度边界 25℃ 近似连接电缆本体温度外护套表面设为“对流传热”边界换热系数给 8 W/(m²·K)环境温度 25℃。如果你仿真是直埋或管廊敷设这个边界就要换成土壤热阻模型数值上会有明显差异。固体力学接口里我在接头中截面施加“对称”约束两端面施加“固定约束”模拟接头被两端电缆锚固、无法自由伸缩的实际情况。热膨胀产生的应力主要就体现在轴向上如果两端不加约束模型会变成自由膨胀应力结果趋近于零等于白算。这个点我特别提醒一下很多人第一次做热力耦合漏了固定约束最后出来应力云图稀稀拉拉原因就在这。3.3 网格剖分策略网格是整个耦合仿真里决定成败的一步。我的经验是结构性网格比自由三角形网格好控制但复杂几何做起来太费时间所以折中方案是整体用自由三角形网格在关键边界上加边界层。具体的剖分策略分三层第一层是全局尺寸。整个模型的最大单元尺寸设为 5 mm最小单元尺寸 0.2 mm这样保证网格数量不至于太夸张。第二层是局部加密。应力锥端部、压接管端部、绝缘和半导电层界面处我手动添加了尺寸控制域最大单元尺寸设为 0.5 mm。电场和应力的最大值都出现在这些几何突变点网格加密之后数值会更可信。第三层是边界层。在导体外表面、绝缘内外表面这些有热流和电场梯度的边界上添加了 4 层边界层第一层厚度取 0.05 mm增长率 1.2。边界层能明显改善边界附近的温度梯度和电场梯度分辨率特别是在评估绝缘表面切向电场时不加边界层的解会偏保守加多了计算量陡增4 层是个性价比不错的选择。网格剖分完成后记得看单元质量和单元数量。我这次全模型大约 35 万单元在普通工作站上稳态全耦合求解一个工况需要十几分钟属于可以接受的节奏。单元质量最小不要低于 0.1低于这个值的地方大概率是几何里出现了异常小的尖角缝隙建议返回几何处理而不是硬着头皮求解。3.4 求解配置与收敛技巧求解设置里我先建了一个“稳态”研究把三个物理场全部加入。COMSOL 6.3 默认的全耦合求解器会自动处理多物理场耦合但直接硬解容易不收敛。我的经验是先给“固体力学”单独跑一个纯力学的预分析给后续耦合迭代提供一个合理的初始位移场然后切到全耦合设置阻尼因子为 0.6最大迭代次数 25 次。实际跑下来大概 8 次迭代就收敛了稳定得超出预期。如果你遇到全耦合不收敛还有个思路是把“电磁热”和“热膨胀”两个耦合节点暂时停用先算一个纯电场传热稳态把温度场结果作为初始值再开启全部耦合节点继续求解。这种方法本质上是给求解器一个更接近真解的起步点而不是让它在黑暗中摸索。我在调试初期反复发散用这个方法后基本一次通过。还有一个小细节研究配置里要把“电流接口”的求解变量改为“稳态初始值”。如果不改耦合求解时电场会在第一个迭代步突然加载传热方程来不及响应局部温度出现跳变求解器很容易直接报错。4. 后处理与结果解读从云图到工程结论4.1 温度场怎么看稳态求解完先看温度场。典型的中间接头温度分布是导体压接管温度最高向外逐渐降低绝缘层内部有很明显的径向温度梯度。因为介质损耗在绝缘里也存在其温度曲线不会是一条直线的传导衰减形态而是略带抛物状的总趋势。我这次稳态工况下导体压接管表面温度约 72℃外护套表面约 42℃温升 17℃在合理范围内。XLPE 绝缘最高温度只要不超过 90℃ 长期耐受线这个接头设计就有足够的散热裕度。在看温度云图时我最建议做两件事一是沿径向多取几条探针线看温度梯度是否光滑二是把导体温度单独剖出来核对它与输入的载流量之间是否符合预期的阻抗发热关系。如果导体温度明显偏低但绝缘温升很高大概率是介质损耗参数设置过大或者相反。4.2 电场与应力场对应关系电场结果要看应力锥端部区域。理论上应力锥端部的峰值场强是整个接头里最高的地方也是绝缘设计裕度的控制点。我这次算出来峰值电场出现在应力锥尖端和硅橡胶增强绝缘交界附近场强约 3.5 kV/mm对于硅橡胶材料来说留有充足裕度。应力场方面最大等效应力出现在压接管端部因为这里几何刚度突变热膨胀后产生楔入效应。应力量级在 1.5 MPa 左右而硅橡胶的拉伸强度至少十几个兆帕短期看没有问题。但要注意聚合物材料在长期热载荷下存在应力松弛稳态应力低不代表寿命内不出事我后面用派生值做了一条沿轴向的应力分布曲线去看变化趋势。最有意思的结果是电场和应力场的高度相关性应力集中区和最大场强区基本处于同一位置。这就验证了工程上“把应力锥做光滑、做长”的经验——几何突变越缓界面压力分布越均匀电场和应力就越不容易同时冲到高点。4.3 用探针和派生值评估接头寿命余量云图直观归直观工程上做结论还得靠数值。我在模型中添加了三个“全局探针”导体最高温度、绝缘域最大电场、应力锥表面最小接触压力。这三个值分别对应热老化、绝缘击穿、界面微放电三个失效模式。尤其是接触压力这个参数很多人会忽略。预制接头靠界面压力保证贴合紧密度如果运行中接触压力降到零界面就会分离产生气隙局放随之而来。我通过派生值计算了应力锥和 XLPE 绝缘之间的界面正应力分布结果显示在满载发热状态下界面压力相比初始装配值下降了约 18%还在可控范围但这个下降量在 1.5 倍过载工况下会更严重值得进一步优化。把三个探针的数值和材料许用值、寿命模型结合起来就能粗略估计接头的长期运行可靠性。比如阿伦尼乌斯寿命模型需要输入热点温度把探针测出的导体温度填入即可得到相对寿命损耗率最大场强和材料击穿场强一对比就得到绝缘裕度系数。这些指标都是设计阶段可以直接用来把关的。5. 常见问题与排错实录5.1 不收敛的常见诱因全耦合稳态求解最容易爆“不收敛”的地方通常是电磁热耦合和接触力学交叉的位置。如果是初次做先做一个“降温试验”把电流载荷降低到额定值的 10%如果这样都不收敛说明模型本身有毛病大概率是边界条件或材料属性填错如果降温后能收敛就逐步提高载荷用连续延拓法逼近目标工况。另一个常见诱因是力学材料参数数量级失真。I 单位制下弹性模量必须用 Pa而不是 MPa。把 2 MPa 填成 2 Pa 会导致材料极软位移结果会飞上天反过来把 250 MPa 填成 250E6 Pa 其实是对的但很多人会因为看到“Pa”这个单位就怀疑自己填错了。我的建议是每个材料参数填完后先在纯力学研究中单独跑一次自重/热膨胀测试看变形是否合理再继续。5.2 结果震荡与突变温度场或电场分布不光滑通常是网格和边界条件不匹配的信号。我遇到过一次绝缘层中场强出现了明显的齿状跳变复查发现是绝缘域被几何布尔运算切成了两个域两个域的材料属性虽一样但网格在边界处不连续导致场量在跨界面处出现数值振荡。解决办法是用“联结”特征把两个相邻域的边界设为一致或者重新建模避免多余的分域。周期性震荡还有一个来源是求解器阻尼过大。如果看到残差曲线像心跳一样上下跳动把阻尼因子调到 0.2或者改用“分离”式求解器让电场先单独收敛再热-力迭代能快速定位是哪一组变量在互相拖后腿。5.3 网格敏感性评估网格密度直接影响尖角处的峰值场强所以必须做网格敏感性分析。我这次把应力锥端部网格从 0.5 mm 逐步加密到 0.2 mm再看峰值电场的变化量。如果加密后峰值场强变化小于 5%说明解已经进入网格不敏感区如果变化很大就得继续加密或者接受保守值。网格加密时也要注意峰值场强只集中在很小一个区域全局网格都加密会导致计算时间陡增。正确做法是利用 COMSOL 的“自适应网格细化”功能或者手动在峰值区域画小尺寸控制域。工程上不需要追求绝对收敛判断标准是网格翻倍后结果变化能否接受。5.4 排错速查表现象可能原因处理办法全耦合不收敛初始值距离真解太远先单场预求解再加耦合温度异常高导体电阻率未设温度依赖打开材料温度依赖选项温度异常低历史遗留的单位制错误核对密度与比热容量纲应力结果接近零缺固定约束两端加固定约束中间加对称电场出现锯齿跳变域被分割、边界网格不一致使用联结特征或重新建模结果随网格剧烈变化网格未收敛加密局部关键区并做敏感性分析瞬态分析奇慢时间步设得过细初试每 60 秒一步稳定后加密6. 仿真成果落地与工程扩展6.1 设计验证与结构优化仿真跑完了不能停在“云图好看”这个层面。我通常会把仿真结果和设计规范里的参考数值对齐比如验证最大场强与材料短时击穿场强之间保留至少 2~3 倍安全系数验证长期运行温度不超过 90℃ 的使用限制。满足的话设计就有数据支撑不满足就要对着几何和材料做优化。优化方向一般有两个。一是结构参数优化比如加大应力锥长度、调整锥面角度、增加绝缘层厚度每个方案重跑一次多物理场计算用探针数值做对比。二是材质优化比如换用导热系数更高的硅橡胶填料能直接把热点温度压下来。几次迭代之后我最后确定的方案在额定载流量下热点温度比初始设计低 6℃同时应力锥端部场强下降了约 12%效果相当明显。有一点我必须提醒仿真优化出来的几何不能直接当作最终产品设计还要考虑制造工艺、公差、装配应力等因素。比如应力锥锥面角度优化后角度更小了但模压脱模可能成问题。所以仿真更多是提供“方向性结论”具体尺寸还要工艺工程师拍板。6.2 扩展到暂态、寿命与序列化分析稳态工况是所有耦合分析的地基但它只是第一层。工程里还有短路冲击、负荷突变、冷启动等暂态工况这时候需要把研究类型从“稳态”切换成“瞬态”在载荷里加上时间函数模拟电流的突变过程。瞬态模型可以直接复用稳态模型的所有设置只需要在电流接口的求解配置里为时域分析调整载荷圆滑因子避免阶跃载荷产生非物理振荡。寿命预测方面阿伦尼乌斯模型是工程上最常用的热老化模型其核心输入就是热点温度。把稳态和暂态仿真得到的温升数据代入模型就能估算接头的加速老化因子和剩余寿命。再加上电场强度对寿命的协同影响可以做电-热联合老化评估比单看温度或者单看场强要符合实际得多。最后一个建议是把自己的模型参数、网格设置、求解配置整理成一个标准模板。中间接头结构在项目间变化不大大部分几何尺寸和物理场设置可以直接套用。我这边已经把整个模型打包成库函数新项目只需要改几何参数和材料参数一个白天就能完成计算和分析。这个习惯帮我从重复劳动里解放出来后面再做其他电缆附件仿真都是一条成熟流水线。我个人看完这轮电 - 热 - 力耦合结果后最大的体会是多个物理场耦合不是三个仿真绑在一起而是要把自己切换到“系统视角”去看问题。单纯降低场强可能牺牲散热单纯降低温度又可能引入应力集中只有在一个模型里同时看它们彼此制约才能判断最终设计是不是真的更优。这也是为什么我建议不要跳过耦合仿真直接做样品实测——实测只能告诉你结果好坏耦合仿真能告诉你为什么好、为什么坏。最后分享一个小技巧每次跑新工况前我都在“结果”里保存一个基准快照包含温升、最大场强、最大应力三个探针值。方案改完就跑一遍对比快照三个数字的变化一目了然。这个习惯帮我追回了好几次因为参数误调整带来的方向性偏差也让我能快速在多种结构方案里选出更优的那一个。仿真做到一定程度比的不是软件操作谁更熟练而是对物理机理的理解深不深。电 - 热 - 力耦合仿真不过是个工具真正值钱的是你能不能用这套工具去回答“这个接头在真实运行条件下到底会发生什么”。如果你正准备做类似的仿真希望这份记录能让你少走几步弯路。