自由空间线圈产生的谐波磁场的分析模拟.docxVIP

自由空间线圈产生的谐波磁场的分析模拟郭宗林（哈工大物理系1411401班，1141140116）摘要：利用毕奥-萨伐尔定律给出圆电流所激发磁场磁感应强度的积分表达式，用ANSYS有限元分析对线圈产生的谐波磁场进行分析，模拟绘图。关键词：线圈 磁场 交变电压 ANSYS有限元分析 模拟绘图萨伐尔定律结合磁感应强度叠加原理，可以计算任意形状电流所产生的磁场。本文先利用毕奥-萨伐尔定律对圆电流激发磁场进行计算；再使用ANSYS有限元分析对线圈产生的谐波磁场进行分析模拟，并绘出三维图像。圆环电流的磁场分布对于P点，电流元所产生的磁感应强度为： d=用R表示电流半径：d=Rd=Rd=-=(x,y,z)-(,,0)=(x-+(y-)+zd==式中，r=d==d==d=d=上面的积分表达式被积函数都不是初等函数，无法得到其精确的解析解，可以利用MATLAB等计算软件近似计算器数值解。但考虑绘图、导出数据的便捷，下面采用ANSYS软件进行分析和模拟绘图。2、用ANSYS有限元分析对线圈谐波磁场进行模拟条件描述：一个载压线圈，电压为余弦交变电压，各材料参数如下表材料特性几何特性载荷线圈：相对磁导率=1.0线圈匝数n=500空气：相对磁导率=1.0S=0.02（线圈宽度）=12v电阻率（线圈）：=R=1.5*Sf =60Hz=2f该线圈为圆形对称，产生的电磁场在线圈的任一竖直截面（如下面草图）上是相同的，而截面对于电磁场是对称的，因此计算截面的1/4区域即可。假设大圆之外几乎没有磁场，而小圆和大圆之间的区域看成是远场区域，即磁场较小。于是得到如下模型，在r=6s到r=12s区域为远场区，r=12外的区域几乎无磁场，忽略不计。采用ANSYS进行电磁场分析，可分为三个阶段：前处理、划分网格和施加求解条件、求解和后处理。2.1前处理这一阶段包括创建物理环境、定义单元类型和定义各部分材料性质。具体步骤包括过滤图像界面，设置单位制以创建物理环境。然后定义单元类型，这里用到了3种单元类型：PLANE53模拟空气，带有CURR和AZ自由度的PLANE53模拟载压线圈，INFIN10模拟远场区域；再定义材料特性，输入材料的相对磁导率和电阻率。再绘制几何模型，利用GUI操作或命令流分别创建模拟线圈、空气、远场区域的几何模型，再叠加可以得到如左图所示结果。最后定义分析参数，在ScalarParamenters中输入一下参数：n=500s=0.02r=0.03 至此前处理阶段基本完成。2.2划分网格和施加求解条件为了便于划分网格，要把坐标系由全局笛卡尔坐标系改为全局柱坐标系。然后设置网格密度并划分网格，其中远场区域线上设定划分单元数为1，模型上所有区域网格数是8，线圈区域网格数为“s/4”。最终得到的网格划分结果如左图所示，加载求解条件主要包括加载边界条件和加载载荷，其中对Y轴上的节点加载磁力线平行边界条件，选择适当的分析类型，最后给线圈施加电压载荷并设置激励电压频率。2.3求解和后处理利用GUI操作或命令流执行求解命令，待求解成功后就可以通过后处理观察求解结果。首先查看磁力线分布，利用PLF命令可以快速绘出磁力线分布图（如下面左图所示）磁力线分布图扩展后的三维图像 这里绘制的磁力线（通量线）是一步载荷所产生的，谐性分析中，线性分析求解按一步进行，交变电流产生的是交变磁场，所以这里默认模式下绘制的是实部（叫变量用复数表示），可通过GUI操作查看虚部磁力线分布，结果与此图非常相似。利用扩展命令将上图所示的磁力线等值图绕Y轴旋转90度成一个三维实体（如上右图），再改变视角方向得到扩展后的节点磁流密度等值云图。除此之外，还可以用PRESOL命令快速打印节点磁通密度和磁场强度、磁能等量。（说明：节点磁通密度和磁场强度对于谐波分析是瞬时值：在WT=0和WT=-90时的实部/虚部，对于磁能除是以上瞬时值外还是单元总和），虽然调用的结果有诸多的限制，但对于磁场的研究仍有很大的价值。3.结果分析3.1图像分析从输出的二维磁力线分布可以看出，磁力线从线圈穿出，绕过其外部空间构成闭合曲线（图像只显示处四分之一的图像）。其次，就磁力线疏密而言，越靠近线圈的地方磁力线越密集，距线圈距离越大，磁力线越稀疏，说明磁感应强度在不断减小。3.2特殊位置节点处数值分析利用PRESOL操作打印的各节点数值是瞬时值，但也很好地刻画了周围空间各物理量的变化规律；节点也限于xoy平面上，但因为其高度对称性，结果仍不失一般性。下表（表3.2.1）是沿X轴各节点处的磁感应强度大小（B），磁感应强度在X,Y方向的分量大小（Bx,By）,以及节点所在网格单元的磁能总和。（因为