铜芯电缆温度分布MATLAB计算.docVIP

PAGE \* MERGEFORMAT1 题目 如图1所示铜芯电缆，电流为5000A，内径为10mm，外包材料聚氯乙烯的厚度为2mm，导热系数为0.15+0.00013{t}。电缆左半边为绝热边界条件，右半边为第三类边界条件，空气温度为20℃，绝缘层表面与环境间的复合表面传热系数为10。铜的电阻率为，，，t的单位为摄氏度。试通过数值方法求解温度分布。 图 SEQ 图 \* ARABIC 1 编程计算 控制方程 根据题意，本题为二维稳态导热问题，其控制方程为： 边界条件： ： ： 其中：。 方程离散 为建立通用方程，考虑非稳态项的控制方程为： 采用全隐格式，在时间内，对控制容积积分，整理后可得： 其中： ，， ，， ，， 采用通用表达式，各表达式如下表： 表1 坐标及系数表达式 坐标系 极坐标 通用表达式 东西坐标 南北坐标 半径 东西尺度系数 东西节点间距 南北节点间距 东西导热面积 南北导热面积 控制体体积 边界条件处理 对于北边界，采用附加源项法处理。由于北边界()为第三类边界条件，则最靠近边界的控制容积加入以下附加源项： 其中： 将附加源项加到相应控制容积后，再令相应的。 对于南边界，可认为定温边界条件，由于其导热面积为零，。 对于东西边界，计算时取计算区域，故东西边界重合，可认为为定温边界条件，温度为上一层相邻控制容积的温度。 导热系数与计算 取铜导热系数为常数，。 每个控制容积各界面对应导热系数分别为、、、。对于铜芯或保温层内部控制容积，各导热系数均为常数。两者交接界面的导热系数用调和平均法计算。 方程求解 方程采用ADI-TDMA方法求解，首先在Y方向进行隐式计算，X方向采用显式计算。各方向对应方程为三对角矩阵，使用TDMA法求解。然后再在X方向进行隐式计算，Y方向采用显式计算。  结果输出与分析 计算结果 程序中温度T为二维数组，采用坐标变换方法，将温度表示在极坐标系中。 设定温度初场为23℃，循环结束判定条件为，网格数为条件下，输出结果如图2： 图 SEQ 图 \* ARABIC 2 网格独立性考察 保持迭代精度不变： 网格数为时，计算结果为： 图 SEQ 图 \* ARABIC 3 网格数为时，计算结果为： 图 SEQ 图 \* ARABIC 4 网格数为时，计算结果为： 图 SEQ 图 \* ARABIC 5 结论：从以上各图可以看出，程序运行结果与网格划分无关，程序具有较好的网格独立性。 收获与体会 通过这次matlab编程作业，我对二维扩散问题有了更加深刻的理解，对网格划分、通用离散形式、边界条件处理等有了进一步的认识。在编写Matlab程序过程中，我为了直接求解三对角矩阵还曾编写一个Solution.m文件，经过对比后发现此文件相比于TDMA方法在速度上稍微快一点，结果基本相同。 通过编程，我更加深刻的认识到只有亲自动手才能加深对问题理解，才能真正获得属于自己的知识。  程序语句 程序采用Matlab编写，主要分为4部分，分别是主程序，用于给定题目条件，调用其他函数，循环求解等；网格划分函数Grid.m，用于划分网格；SolutionTDMA.m，用于执行交替隐式计算；TDMA.m，用于求解三对角矩阵。 Main.m clear all clear global format long global X global Y global dX global dY global DX global DXn global DXs global DY global Cv global CV global T global T0 global Tf global nodX global nodY nodX=200; nodY=350; %给出题目参数 X=2*pi; Y=7E-3; Grid; %划分网格 Tf=20; h=10; %计算导热系数 for i=1:5/7*nodY Le(i)=400;%假设铜的导热系数为400W/(m.K) Lw(i)=400; end for i=5/7*nodY+1:nodY Le(i)=0.15; Lw(i)=0.15; end for i=1:5/7*nodY-1 Ln(i)=400; end for i=5/7*nodY+1:nodY Ln(i)=0.15; end Ln(5/7*nodY)=2/(1/Ln(5/7*nodY-1)+1/Ln(5/7*nodY+1)); for i=1:5/7*nodY Ls(i)=400; end for i=5/7*nodY+2:nodY Ls(i)=0.15; end Ls(5/7*nodY+1)=2/(1/Ls(5/7