地下电力电缆温度场及载流量的数值计算.pdfVIP

东 北 电 力 大 学 学 报 第28卷第4期 JournalOfNortheastDianliUniversity Vo1．28．N0．4 2008年 8月 NaturalScienceEdition Aug．，2008 文章编号：1005-2992(2008)04—0062—04 地下电力电缆温度场及载流量的数值计算 于建立，常树生，牛远方，葛丽婷 (东北电力大学 电气工程学院，吉林 吉林 132012) 摘 要：对地下电力电缆温度场和载流量进行了研究。基于有限差分法，采用直角坐标和极坐标 相结合的思想，对计算区域进行不均匀网格的划分。对电缆及其周围土壤的温度场进行迭代，并进行两 种坐标系之间的插值计算 ，得到一个收敛的温度场。多次循环后电缆缆芯达到上限温度值，可得电缆的 载流量。 关 键 词：电力电缆；温度场；载流量；有限差分；坐标组合 中图分类号：TM247 文献标识码：A 在地下电缆运行过程中，为避免缆芯通过的电流过大导致电缆绝缘因超温而损坏，通常电缆的运行 电流都小于额定载流量，这样就造成了缆芯材料的浪费。对电缆的温度场和载流量进行准确计算具有 重要意义。受电缆初始状态、环境温度、埋设区域的几何特点，及土壤复杂程度等因素的影响，而国际电 工委员会(IEC)的标准只能解决均匀土壤的计算问题，已经不能满足实际工程的需要。为了尽可能准 确的模拟电缆埋设区域的温度场，必须借助数值方法。一直以来，解决这一问题的数值方法有：有限元 法…、边界元法 J、有限差分法 3等。本文采取的数值计算方法是在有限差分法基础上进行的坐标组 合方法。这种方法的优点是可充分考虑电缆的内部结构对传热的影响。电缆区域采用极坐标，电缆外 的土壤区域采用直角坐标，这样能够将一些参数 (如媒质的导热系数等)考虑为温度的变量，不只是作 为简单的常数来处理。该模型适用于不同的铺设形式和土壤情况比较复杂的工况，且计算精度高。为 了减少计算时间，本文采用不均匀网格．将土壤区域离散的传热方程采用交替方向块方法计算，加快求 解速度。在求得一个达到收敛标准的温度场时，计算出此时缆芯的温度，转而进入载流量的求解。本方 法编程简单，计算速度快，均匀土壤情况下与IEC-60287 标准 比较，证 明本方法正确。 1 物理模型和数学模型 1．1 物理模型 根据电缆发热的特点，可以把此问题看做是具有内热源的闭域导热问题。有限差分法在计算导热 问题时边界条件很关键 ，若根据电缆的物理形状单纯采用极坐标计算，土壤区域的外边界温度无法确 定；若只采用直角坐标计算，电缆区域比较难处理。故本文采用直角坐标网格与极坐标网格组合的方法 对电缆埋设区域的温度场进行数值计算，称为坐标组合法。土壤区域采用直角坐标，电缆区域采用极 坐标。 收稿 日期：2008一o4—25 作者简介：于建立(1983一)，男，东北电力大学电气工程学院在读硕士研究生 第4期 于建立等：地下电力电缆温度场及载流量的数值计算 63 1．2 数学模型 根据物理模型，本文采用通用的导热微分方程 J： pc警=未(A)+(A)+亳(A)+s (1) 式中，JD，C和A分别为密度 (kg／m)，比热 (J／kg·cc)和导热系数 (W／m ·k)；T为温度(℃)；t为时问 (s)；S为内热源。电缆的轴向长度远远大于其径向长度，在长度方向温度基本没有变化，可将本文模型 由三维问题简化为二维问题。 1．2．1 直角坐标导热方程 稳态情况下pc0i1=0， 土壤区域无内热源S=O，直角坐标导热方程为： 击(A)+(A)：0 (2) 计算过程中，采用以控制容积中心点(以下称节点)温度代替 (6x) (8x) 该控制容积温度的方法。为了提高计算精