电磁场第三次数值仿真报告.docxVIP

电磁场第三次数值仿真题目叙述为了获得一定区域上的匀强磁场，可采用多组Helmholtz线圈结构。一种两对线圈的结构如图1所示。线圈半径a1，a2，线圈间距h1，h2，以及线圈中通过电流i1，i2可变化量，如图1（a）所示。为了定量衡量关注区域的磁场均压程度，过轴线做截面ABo1o2，取CD=0.8×AB和Eo1 = 0.8×o1o2，在CD和Eo1线段上每边均匀取20采样点，从而形成如图1（b）所示的采样节点，定义z方向B的不均压系数为：其中，为所有采样点的z方向磁感应强度平均值；为第n个采样点的z方向磁感应强度值。N为采样点总数。如给定a1 = 1 m， i1=1 A，且h1 0.2 m， 情况下，如果要使得?达到最小，线圈半径a2，线圈间距h1，h2，以及线圈中通过电流i2应如何取值？(a)线圈结构示意图（b）磁场采样节点示意图初步工作推导和计算直线电流元在任意观察点的磁感应强度电磁场课上曾讲过相关推导过程与运算结果。柱坐标系下，对于下图：在观察点A出的磁感应强度为：在直角坐标系下，相应的公式是：对于以上两个公式，在直角坐标系下，有：相应的MATLAB代码如下：（记得说只保留了Z方向分量）积分法计算线圈组在其内任一点的磁感应强度要算线圈组在其内任意一点的磁感应强度B，必须先知道一个载流圆线圈在空间任一点产生的磁场情况。建立空间直角坐标系，作图如下：对于空间任意观察点P，有以下推导：进而得出：相应matlab代码如下：*********************************************************************请简述你编的几个函数的功能每个函数都要说，每个函数用“******”隔开数值法计算线圈组在其内任一点的磁感应强度一方面，考虑到用积分法计算磁感应强度，运算步数多，运算时间长；另一方面，我们已经得出了直流电流元在任意观察点磁感应强度的表达式，基于相关表达式与函数，我们能够较为容易的得出正多边形线圈组在其内任一点的磁感应强度，进而用正多边形代替圆，在误差允许的范围内，进一步减少运算步数，缩短运算时间。在直角坐标系下，对于下图的正多边形（以正六边形代替）：它第i条边的起点坐标是（），终点坐标是（）。利用2.1中已经得出的calculate_b_line函数，我们循环算出每条边在观察点的B，并将其叠加，就得到边数为M的正多边形线圈组，在它们内部任意一点的磁感应强度。相应MATLAB代码如下：*********************************************************************function B=multi_solve_point(a,M,destinationpoint,h1,C) b=0;for i=1:M sita=2*pi/M; x1=a*cos((i-1)*sita); y1=a*sin((i-1)*sita); x2=a*cos((i)*sita); y2=a*sin((i)*sita); bn1=calculate_b_line([x1 y1 h1],[x2 y2 h1],destinationpoint,C); bn2=calculate_b_line([x1 y1 -h1],[x2 y2 -h1],destinationpoint,C); b=b+bn1+bn2;endB=b;*********************************************************************这种方法的优点是在精度要求不太高的情况下，能用M较小的正多边形来代替圆形，计算时间较短。缺点是本方法精度较小，如果要提高精度，就必须要增加边数M，进而计算时间会提高。现用单个线圈轴线上的B来讨论M取不同的值时，运算精度和运算时间时间的变化。理论计算公式是，来源于课件有关题目。运算结果如下：边数M10501001000数值结果B1/T5.6610e-0075.5268e-0075.5227e-0075.5213e-007理论结果B2/T5.5213e-0075.5213e-0075.5213e-0075.5213e-007相对误差2.5314%0.0991%0.0248%0%运算时间 t/s0.0032790.0038570.0052060.031530上表的结果验证了数值法的优缺点。由表中数据可得，在边数M=50时，运算精度和运算时间都较好。因此，在之后的运算中，均用M=50的正多边形来代替圆。不均压系数δ的计算。我们现在已经得到两个线圈中任一点磁感应强度B的计算公式，下面只需按照题目要求进行取点，再代入公式计算δ即可。利用对称法进行优化宏观的看，求