转角汇流条电阻仿真计算.docxVIP

转角汇流条电阻的计算【题目：转角汇流条电阻的计算】转角汇流条几何参数如下：a =1 m；d =1m；c= 0.5m；b= 0.5m；e = 10 mm按二维电流场考虑。设注入电流的两个端面的电压为1 V；采用有限差分方法，编写程序，计算汇流条电流分布，绘制电流线和等电位线；说明如何验证程序的正确性和提高计算准确度的方案；如不采用数值计算方法，请估算该汇流条电阻，比较估算方法结果和有限差分方法的结果，说明两种结果差异的原因。估算方法可参考附件中文章，或者自行提出估算方法。【解答】1、边值问题分析：1）、三维问题转换为二维问题：可建立如下坐标系。由于题中所给数据为：a =1 m；d =1m；c= 0.5m；b= 0.5m；e = 10 mm。a、b、c、d均至少为e的50倍，因此可忽略电势在导体厚度方面的变化，即可将此题放在图中黄色区域（x-y平面）的二维平面上求解。2）边界条件的分析：本题求解的是恒定电流场，考虑到导体本身的特性以及题中所给的电流流向，易得边界条件如下：（S2、S3、S5、S6上）；；；3）场域内部方程：求解区域是导体内部的区域，无自由电荷分布，故易得求解方程为： （导体内部区域）2、有限差分法求解：在用matlab求解该题中的电阻时，可以采用拉普拉斯方程的5点差分法格式进行有限差分求解。将上述的L型区域划分成网格，列写满足上述边值条件的差分方程即可求解。各种情况下的差分方程如下：1）在导体内部的情况：如上所示为导体内部的情况，此时，采用差分法进行计算，有：2）在边界上的情况：如上所示为在导体的边界上的情况，因为这样的点都是满足第二类边界条件的点，即位于平面S2、S3、S5、S6上的点，此时，采用有限差分法进行计算，有：3）在拐角处的情况：如上图所示即为在拐角处的情况，即图中L形导体的“凸点”，“凸点”的两条边都是边界，由于边界处的法向导数都为0，可知边界两侧的电势分布相同，则其5点差分方程为：3、利用matlab进行有限差分法求解：（1）直接使用矩阵求解：使用matlab来进行以上描述的有限差分法求解，我们可以使用如下的矩阵进行求解：其中，满足第一类边界条件时，设置A+I矩阵中相应的行向量为0，B矩阵为该对应点的电位值，易知这样的φ满足第一类条件；满足第二类边界条件时，（A+I）矩阵中对应的行向量应按照以上描述对应的设置为1/4和1/2，B矩阵对应的向量均为0；在导体内部的点，只需将（A+I）矩阵中将各行向量按以上描述设置为1/4，将B矩阵的行向量全设为0即可。这样，设置好A、B矩阵后，即可按照以上公式直接求出各点的φ。这样求解并做出电流的分布如下图：查看上端面的电场强度Ey，基本位于0.7948和0.8124之间，取各点的算术平均值，得Ey=0.7819V/m。降本题中所给的参数，代入，可得：（流入）（2）使用迭代的方法求解：迭代的方程正是在第2点“有限差分法求解”中所阐述的各种情况下的差分方程。根据电位点的位置使用相应的差分方程，进行迭代计算。计算时，取上边界点的电位初值为1V，其它各点电位初值为0。每迭代一次，求出各点的两次计算结果之差，取其中最大的值作为误差值，当误差值小于所设误差值上限时，认为精度已达要求，迭代结束。电阻的求法与矩阵求解中的相同，不再详述。下面为不同网格大小情况下的电流线和等电位线图，误差上限取为，如图所示：网格大小为，即划分为3030网格：求得电阻值为32.985网格大小为，即划分为5050网格：求得电阻值为32.521网格大小为，即划分为7070网格：求得电阻值为32.267对比以上各图可以看出，随着网格越来越小，划分的格点越来越多，等电位线越来越光滑，而且边界处等电位线垂直于边界的情况越来越好，实际上计算结果越来越接近于实际情况。而电阻值在逐渐减小，为了判断其趋势及最终的收敛情况，也如矩阵方法中进一步减小网格面积，增大网格数进行求解，结果如下：网格数303050507070100100200200300300电阻值32.98532.52132.26731.96930.876结果发现，在一定的范围内，电阻的计算值不收敛。分析有两个可能的原因，首先是网格的划分仍然不够细，但是迭代法计算较慢，在网格数为300300时计算了半个小时才得到结果，因此没有继续增大网格数，因此未能找到收敛值；另一个原因可能是迭代法本身的算法存在一定的误差，虽然可以控制误差的上限，但并不能完全消除误差，特别是当划分网格数较多时，计算结果较精确，此时误差的影响反而越来越大，因此在计算电阻值时造成了最终的结果不能收敛于某个固定值。为此取网格数为中间值的四种情况：5050，7070，100100，200200，计算其电阻的均值作为迭代法的电阻值，为31.91。4、结果分析：（1）首先进行等电位线与电流线的分析，根据所