基于磁场积分方程法的磁通门瞬态分析方法
基于磁场积分方程法的磁通门瞬态分析方法
【技术领域】
[0001 ]本发明设及磁通口输出信号的瞬态分析方法，具体设及利用磁场积分方程法对磁 通口输出信号进行瞬态分析的方法。
【背景技术】
[0002] 磁通口主要由高导磁率铁忍及两种线圈构成:在高导磁率铁忍外面总共有两组线 圈，一组沿两根铁忍(或环状铁忍两半)对称顺绕，称为激励线圈，供交变电流，另一组称测 量线圈或接收线圈，它们绕制在激励线圈外面;当存在外磁场时线圈两端出现感应电动势。 磁通口是磁通口传感器的重要组成部分，磁通口传感器是利用被测磁场中高导磁铁忍在交 变磁场的饱和激励下，其磁感应强度与磁场强度的非线性关系来测量弱磁场的一种传感 器。为了制作出性能优良的磁通口传感器，常常需要借助于各种数值分析方法来对磁通口 进行分析和优化。
[0003] 采用有限元法对平面磁通口进行分析，需要十分复杂的步骤设置W及大量的剖分 等处理环节，文南犬"An integrated micro-fluxgate magnetic sensor with front-end circuitry, IEEE Transactions on Instrumentation and Measurement, Vol. 58, 2009, P3269-3275."公开了一种采用有限元对平面磁通口输出信号进行分析的方法，进行 一次磁通口瞬态分析就需要约为120个小时，如果想获得磁通口线性范围、最佳激励条件等 性能指标，还要进行多次瞬态分析，则需要更长的时间。文献"Demagnetizing field in ferroma即etic sheet, Physica B, Vol. 306, 2001, pl72-177."公开了一种采用退磁 系数对磁通口进行分析的方法，虽然解决了有限元法计算时间长的问题，但是精度很差。文 南犬"About the use of numerical integral methods to simulate a fluxgate magnetometer: The ring-core example, Sensors & Actuators A Physical, Vol. 163, 2010, P48-53."公开了一种采用静态积分方程来分析磁通口的方法，虽然能够使得 计算时间大大减少，而且在精度方面也可W达到磁通口的要求，但是它只对磁通口进行了 准静态分析，并没有解决磁通口的瞬态分析问题。因此有必要设计一种磁通口瞬态分析方 法，满足耗时少、精度高的要求。
【发明内容】
[0004] 本发明所要解决的技术问题是克服现有技术的上述不足而提供一种基于磁场积 分方程法的磁通口瞬态分析方法，其能够对磁通口进行快速精确分析和优化。
[000引本发明解决上述问题所采用的技术方案为，基于磁场积分方程法的磁通口瞬态分 析方法，磁通口包括铁忍、激励线圈、接收线圈，铁忍为两条长条状薄片，磁通口铁忍数学模 型简化为平面二维模型，其特征在于，磁通口瞬态分析步骤包括： 首先计算磁通口铁忍的形状参数并测试铁忍材料的磁化曲线，建立准静态模型； 然后在准静态模型中引入激励电流频率的影响：在低频情况下可W忽略位移电流，电 流密度可分解为满流I:和给定的激励电流密度％，即!-?净I:，根据磁场叠加原理，考虑满 流影响后，磁通口铁忍内磁场强度1?可w分解成：
其中E是源磁场，包括激励电流产生的磁场强度和被测外磁场，11?是满流产生的 寄生磁场，?是铁忍材料产生的退磁场； 将满流产生的磁场强度移项到左端，并记堯红議哮;龜;，可得秦::二義束錢。裝自，铁忍内 磁场方程组为：
对于薄片状、圆柱状铁忍，可W用铁忍形状参数和磁感应强度时变率表示满流产生的 寄生磁场，其中铁忍材料数学模型??为：
最后建立磁通口瞬态分析模型， 在考虑频率的外磁场8中，铁忍材料内的磁场親与磁化强度夠的关系由W下方程组描 述：
其中3为形状参数，1?为磁滞损失参数，:S劝平均场参数，C为磁畴壁弯曲常数，为满 流产生的寄生磁场是与铁忍厚度和导电率有关的系数，为铁忍材料的饱和磁化 强度，鴻为不可逆磁化强度，該sb为非磁滞磁化强度，鋪为铁忍材料的磁化强度，B为磁感 应强度，Η为磁通口铁忍的总磁场，参数反。、：^、^1、％和^描述非磁滞磁化曲线的形状，线 性部分由Ra描述，而饱和部分由和、描述，S为方向参数，当地/dt〉0时裹取+1，当 地/化<0时瑟取-1; 联合铁忍材料的磁化曲线及W上公式、方程组，对磁通口进行瞬态分析。
[0006 ]建立准静态模型的步骤包括： 第一步，计算磁通口铁忍的形状参数并测试铁忍材料的磁化曲线； 第二步，建立磁场强度模型， 根据磁场叠加原理，磁通口铁忍任意一点的总磁场可W分为Ξ部分： 藉=:.裹;手變办難議..............................（1) 其中薇和顯为源场，囊为激励电流产生的激励磁场，蠢X为被测外磁场，龜为铁忍材 料产生的退磁场，揉为总磁场强度，旋度为零； 第Ξ步，对磁通口铁忍进行剖分，计算铁忍中每个剖分单元的磁化强度， 将磁通口铁忍剖分为η个单元，假设每个单元内的磁化强度是均匀的W及磁荷只在单 元之间的交界面和铁忍与空气的边界出现，可得：
其中Κ为第i个单元的总磁场强度，为激励电流在第i个单元产生的磁场，》棘为被 测外磁场在第i个单元产生的磁场，Sij为铁忍材料在第i个单元内产生的退磁场，苗为第i 个单元的磁化强度，下标X和y分别表示在铁忍平面X轴和y轴方向的磁场分量； 第四步，分析激励磁场强度， 根据毕奥-萨伐定律，激励电流在场点r处产生的激励磁场强度可W表述为：
其中为激励电流方向，i为激励电流密度，?为磁场的坐标，梦为电流的坐标，dV为激 励电流体积微元； 第五步，计算磁通口铁忍材料产生的退磁场， 根据：
退磁场旋度为零，存在一个标量满足:一聚藥，那么标量所满足的方程为：
铁忍材料所有剖分单元在第i个单元产生的磁场为：
其中霧为第j个单元对应的磁标势，r为源点和场点之间的距离，Sj为第j个单元的磁 化强度，苟为第j个单元的外表面，其中京为第i单元沿外表面法线单位矢量； 第六步，建立准静态模型， 将式(3)和式(6)代入式(2)可得：
其中雜?为退磁矩阵，将式(7)与磁通口铁忍材料的磁化曲线联合，建立磁通口的准静 态分析模型，求解在受到被测磁场和激励电流产生的激励磁场作用时铁忍内部的磁场和磁 化强度分布。
[0007]本发明的优点在于:采用本发明建立的磁通口瞬态分析模型不仅准确地综合考虑 了磁通口的结构、铁忍材料、激励电压、激励频率等因素的影响，而且准确地考虑了随着频 率增大而增强的磁感应强度滞后现象，从而能够准确的对磁通口进行瞬态分析，W便优化 磁通口的设计。
【附图说明】
[000引图1是磁通口结构示意图。
[0009] 图2是正弦电压激励的长条型磁通口电路。
[0010] 图3是被测磁场为30μΤ、激励信号为幅值为1.2V、激励频率为1曲Ζ的正弦电压时， 线圈电流仿真结果和实验结果比较图。
[0011] 图4是被测磁场为30μτ、激励信号为幅值为1.2V、激励频率为1曲Ζ的正弦电压时， 接收线圈输出电压的仿真结果和实验结果比较图。
[0012] 图5是激励频率为化化时的仿真磁滞回线和实验磁滞回线比较图。
[0013] 图6是被测磁场为30μΤ时和激励频率为化化时，随着激励电压增大，输出二次谐波 电压幅值与激励线圈电流幅值的关系。
【具体实施方式】
[0014] 下面结合附图、实施例对本发明进一步说明。
[0015] 实施例，如图1所示长条形磁通口，包括铁忍1、激励线圈2、接收线圈3,激励线圈2 和接收线圈3均匀绕制在铁忍1外，Ξ者组成磁通口。图2为正弦电压激励的长条型磁通口电 路，R为激励线圈的电阻。铁忍1为两条长条状薄片，厚度远小于其它尺寸，可W认为薄片内 部的磁场主要平行于其表面，垂直方向的磁化强度可近似忽略不计。磁通口接收线圈3中感 应到的是沿铁忍平面的磁场分量产生的磁通，而与之垂直的铁忍横截面内的磁场分量不会 对输出线圈电压产生影响，因此磁通口铁忍数学模型可W简化为平面二维模型。
[0016] 第一步，计算磁通口铁忍的形状参数并测试铁忍材料的磁化曲线 铁忍的形状参数、铁忍材料磁化曲线的计算测量方法已有现有技术公开，不再寶述。
[0017] 第二步，建立磁场强度模型 根据磁场叠加原理，对于平面二维模型，磁通口铁忍任意一点的总磁场可W分为Ξ部 分：
(1) 其中嚷和藏为源场，1|为激励电流产生的激励磁场，顯为被测外磁场，麵为铁忍材 料产生的退磁场，再为总磁场强度，旋度为零。
[0018] 第Ξ步，对磁通口铁忍进行剖分，计算铁忍中每个剖分单元的磁化强度 将磁通口铁忍剖分为η个单元，假设每个单元内的磁化强度是均匀的W及磁荷只在单 元之间的交界面和铁忍与空气的边界出现，可得：
其中K为第i个单元的总磁场强度，it为激励电流在第i个单元产生的磁场,??为被 测外磁场在第i个单元产生的磁场，为铁忍材料在第i个单元内产生的退磁场，瑪为第i 个单元的磁化强度，下标X和y分别表示在铁忍平面X轴和y轴方向的磁场分量。
[0019 ]第四步，分析激励磁场强度： 根据毕奥-萨伐定律，激励电流在场点r处产生的激励磁场强度可W表述为：
其中^为激励电流方向，i为激励电流密度J 为磁场的坐标，？f为电流的坐标，dV为激 励电流体积微元。
[0020] 第五步，计算磁通口铁忍材料产生的退磁场： 根据：
退磁场旋度为零，存在一个标量满足:1|:=心专缘，那么标量所满足的方程为：
铁忍材料所有剖分单元在第i个单元产生的磁场为：
其中嚷为第j个单元对应的磁标势，r为源点和场点之间的距离，为第j个单元的磁 化强度，覆为第j个单元的外表面，其中％.=敬i，京为第i单元沿外表面法线单位矢量。
[0021] 第六步，建立准静态模型 将式(3)和式(6)代入式(2)可得：
其中犧?为退磁矩阵。将式(7)与磁通口铁忍材料的磁化曲线联合，建立磁通口的准静 态分析模型，即可求解在受到被测磁场和激励电流产生的激励磁场作用时铁忍内部的磁场 和磁化强度分布，从而能够对磁通口进行准静态分析。
[0022] 第屯步，在准静态分析模型中引入激励电流频率的影响： 由于准静态分析模型是针对低频情况建立的模型，在低频情况下可W忽略位移电流， 电流密度可分解为满流拓和给定的激励电流密度％，即耗％料0。根据磁场叠加原理，考虑 满流影响后，磁通口铁忍内磁场强度寵（注：与前面用義表示铁忍内磁场不同，为方便起 见，用1表示铁忍内磁场)可W分解成：
其中ii。是源磁场，包括激励电流产生的磁场强度和被测外磁场。?暴S是满流产生的 寄生磁场，1%是铁忍材料产生的退磁场。
[0023] 将满流产生的磁场强度移项到左端，并记堯二％去赛a，可得藻二铁忍内 磁场方程组为：
对于薄片状、圆柱状铁忍，可W用铁忍形状参数和磁感应强度时变率表示满流产生的 寄生磁场，其中铁忍材料数学模型if!)为：
第八步，建立磁通口瞬态分析模型： 建立考虑频率的铁忍材料的数学模型，在考虑频率的外磁场採中，铁忍材料内的磁场 瑪与磁化强度級觀关系由W下方程组描述：
其中B为磁感应强度，3为形状参数，&为磁滞损失参数，α为平均场参数，;C为磁畴壁弯 曲常数。琴1为满流产生的寄生磁场晦iSyi，方程组的第一式可转换为i二錢S -鬚#歡，犹是 与铁忍厚度和导电率有关的系数。M。为铁忍材料的饱和磁化强度，絳《为不可逆磁化强度， 擁《^为非磁滞磁化强度，錄为铁忍材料的磁化强度。参数&。、而、麵、游和呜描述非磁滞磁 化曲线的形状。线性部分由觀描述，而饱和部分由觀巧鸣描述。S为方向参数，当 dH/dt〉0时瑟取+1，当dH/dt<0时送取-1，运样保证了不会出现不符合实际的总能量增加的现 象。
[0024] 联合铁忍材料的磁化曲线、式（10) W及式（11)所示的考虑频率的铁忍材料数学模 型，即可建立磁通口瞬态分析模型，能够准确考虑满流对铁忍材料的影响，从而能够对磁通 口输出信号进行瞬态分析。
[0025] W图1所示的磁通口采用瞬态分析模型进行验证。激励线圈2的应数为330,接收线 圈3的应数为1100。铁忍1为长度30mm，宽度1mm，厚度0.1mm的坡莫合金。铁忍材料由式(9)描 述，其模型参数为=1.0694 X 1〇-6，c=0.3，k=l. 7143，Ms=5.98 X 10日，Ru=0.3997，册=0.0001， Hl=1.1263Xl〇-2〇，aO=0.0827，al=0.0201，N=0.0014。图 2为正弦电压激励的长条型磁通口 电路。R为激励线圈的电阻，Κ=13.3Ω。
[0026] 采用正弦电压信号对磁通口进行激励，比较考虑频率的基于磁场积分方程法的磁 通口模型的仿真结果和实验结果，其中仿真结果是结合分析模型并通过仿真计算软件得 出，实验结果是通过建立的磁通口测量系统得到的，其中信号发生器和功率放大器产生正 弦电压激励信号，直流电源产生被测磁场，示波器用来检测磁通口输出信号。图3和图4分别 表示被测磁场为30μΤ、激励信号为幅值为1.2V和激励频率为1曲Ζ的正弦电压激励线圈电流 W及接收线圈输出电压的仿真结果和实验结果比较图。图5表示激励频率为Ik化时的仿真 磁滞回线和实验磁滞回线比较图。图6表示被测磁场为30μΤ时和激励频率为1曲Z时，随着激 励电压增大，输出二次谐波电压幅值与激励线圈电流幅值的关系。
[0027] 从图3和图4可W看出，因为实际中的激励电流需要一定的响应时间，在开始的一 段时间内，激励电流和输出信号的实验值滞后于仿真模型中的激励电流和输出信号。在后 续的时间条件下，仿真中的激励电流和输出信号和与实验中的激励电流和输出信号吻合良 好。从图5可W看出，仿真磁滞回线与实验测得的磁滞回线几乎一致。从图6可W看出，在小 于95mA的激励电流条件下，由仿真模型得到磁通口二次谐波幅值略低于实验测得的磁通口 二次谐波幅值;在大于95mA的激励电流条件下，由仿真模型得到磁通口二次谐波幅值略高 于实验测得的磁通口二次谐波幅值。从仿真结果与实验结果的对比可W看出，随着激励电 压，激励频率的改变，输出结果准确地反映了激励电流，输出电压等的变化，表明所建立的 磁通口瞬态分析模型可W准确的对磁通口输出信号进行瞬态分析。
【主权项】
1.基于磁场积分方程法的磁通门瞬态分析方法，磁通门包括铁芯、激励线圈、接收线 圈，铁芯为两条长条状薄片，磁通门铁芯数学模型简化为平面二维模型，其特征在于，磁通 门瞬态分析包括以下步骤： 首先计算磁通门铁芯的形状参数并测试铁芯材料的磁化曲线，建立准静态模型； 然后在准静态模型中引入激励电流频率的影响：在低频情况下可以忽略位移电流，电 流密度可分解为涡流U和给定的激励电流密度1￡5，即【产％1知，根据磁场叠加原理，考虑涡 流影响后，磁通门铁芯内磁场强度||可以分解成：其中!^是源磁场，包括激励电流产生的磁场强度和被测外磁场，.是涡流产生的寄 生磁场，是铁芯材料产生的退磁场； 将涡流产生的磁场强度移项到左端，并记铁芯内磁 场方程组为：对于薄片状、圆柱状铁芯，可以用铁芯形状参数和磁感应强度时变率表示涡流产生的 寄生磁场，其中铁芯材料数学模型義I?为：最后建立磁通门瞬态分析模型， 在考虑频率的外磁场Il中，铁芯材料内的磁场与磁化强度麵的关系由以下方程组描 述：其中獅状参数，咖滞损失参数，痛均场参数，畴壁弯曲常数， 流产生的寄生磁场1??，，辨是与铁芯厚度和导电率有关的系数，_为铁芯材料的饱和磁化 强度，Mirr为不可逆磁化强度,?涵为非磁滞磁化强度，議为铁芯材料的磁化强度，B为磁感 应强度，H为磁通门铁芯的总磁场，参数、:11、_、％和％描述非磁滞磁化曲线的形状，线 性部分由装。描述，而饱和部分由_、:_、和描述，§为方向参数，当dH/dt>0时S取+1，当 dH/dt〈0 时S 取-1; 联合铁芯材料的磁化曲线及以上公式、方程组，对磁通门进行瞬态分析。2.根据权利要求1所述的基于磁场积分方程法的磁通门瞬态分析方法，其特征在于，建 立准静态模型的步骤包括： 第一步，计算磁通门铁芯的形状参数并测试铁芯材料的磁化曲线； 第二步，建立磁场强度模型， 根据磁场叠加原理，磁通门铁芯任意一点的总磁场可以分为三部分：其中和1|^为源场，为激励电流产生的激励磁场，为被测外磁场，_为铁芯材 料产生的退磁场，M为总磁场强度，旋度为零； 第三步，对磁通门铁芯进行剖分，计算铁芯中每个剖分单元的磁化强度， 将磁通门铁芯剖分为η个单元，假设每个单元内的磁化强度是均匀的以及磁荷只在单 元之间的交界面和铁芯与空气的边界出现，可得：其中|丨为第i个单元的总磁场强度,?为激励电流在第i个单元产生的磁场，为被测 外磁场在第i个单元产生的磁场，1?...为铁芯材料在第i个单元内产生的退磁场，:轉为第i个 ·'；.￡ ... 单元的磁化强度，下标X和y分别表示在铁芯平面X轴和y轴方向的磁场分量； 第四步，分析激励磁场强度， 根据毕奥-萨伐定律，激励电流在场点r处产生的激励磁场强度可以表述为：其中_为激励电流方向I为激励电流密度J为磁场的坐标Ji为电流的坐标，dV为激励 电流体积微元； 第五步，计算磁通门铁芯材料产生的退磁场， 根据：I 退磁场旋度为零，存在一个标量满足：那么标量所满足的方程为：铁芯材料所有剖分单元在第i个单元产生的磁场为：其中％为第j个单元对应的磁标势，r为源点和场点之间的距离，为第j个单元的磁化 强度，sI为第j个单元的外表面，其中为第i单元沿外表面法线单位矢量； 第六步，建立准静态模型， 将式(3)和式(6)代入式(2)可得：其中辨k为退磁矩阵，将式(7)与磁通门铁芯材料的磁化曲线联合，建立磁通门的准静 态分析模型，求解在受到被测磁场和激励电流产生的激励磁场作用时铁芯内部的磁场和磁 化强度分布。
【专利摘要】本发明公开了一种基于磁场积分方程法的磁通门瞬态分析方法，根据磁通门铁芯特点，建立基于磁场积分方程法的磁通门准静态分析模型，并引入激励电流频率的影响，建立可以对磁通门进行瞬态分析的数学模型。采用本发明建立的磁通门瞬态分析模型不仅准确地综合考虑了磁通门的结构、铁芯材料、激励电压、激励频率等因素的影响，而且准确地考虑了随着频率增大而增强的磁感应强度滞后现象，从而能够准确的对磁通门进行瞬态分析。
【IPC分类】G01R33/00, G01R33/04
【公开号】CN105487025
【申请号】CN201510842978
【发明人】侯晓伟, 李菊萍, 郭俊杰, 倪大成, 王飞, 郑华雄, 郑良广
【申请人】宁波南车时代传感技术有限公司
【公开日】2016年4月13日
【申请日】2015年11月27日
【技术领域】
[0001 ]本发明设及磁通口输出信号的瞬态分析方法，具体设及利用磁场积分方程法对磁 通口输出信号进行瞬态分析的方法。
【背景技术】
[0002] 磁通口主要由高导磁率铁忍及两种线圈构成:在高导磁率铁忍外面总共有两组线 圈，一组沿两根铁忍(或环状铁忍两半)对称顺绕，称为激励线圈，供交变电流，另一组称测 量线圈或接收线圈，它们绕制在激励线圈外面;当存在外磁场时线圈两端出现感应电动势。 磁通口是磁通口传感器的重要组成部分，磁通口传感器是利用被测磁场中高导磁铁忍在交 变磁场的饱和激励下，其磁感应强度与磁场强度的非线性关系来测量弱磁场的一种传感 器。为了制作出性能优良的磁通口传感器，常常需要借助于各种数值分析方法来对磁通口 进行分析和优化。
[0003] 采用有限元法对平面磁通口进行分析，需要十分复杂的步骤设置W及大量的剖分 等处理环节，文南犬"An integrated micro-fluxgate magnetic sensor with front-end circuitry, IEEE Transactions on Instrumentation and Measurement, Vol. 58, 2009, P3269-3275."公开了一种采用有限元对平面磁通口输出信号进行分析的方法，进行 一次磁通口瞬态分析就需要约为120个小时，如果想获得磁通口线性范围、最佳激励条件等 性能指标，还要进行多次瞬态分析，则需要更长的时间。文献"Demagnetizing field in ferroma即etic sheet, Physica B, Vol. 306, 2001, pl72-177."公开了一种采用退磁 系数对磁通口进行分析的方法，虽然解决了有限元法计算时间长的问题，但是精度很差。文 南犬"About the use of numerical integral methods to simulate a fluxgate magnetometer: The ring-core example, Sensors & Actuators A Physical, Vol. 163, 2010, P48-53."公开了一种采用静态积分方程来分析磁通口的方法，虽然能够使得 计算时间大大减少，而且在精度方面也可W达到磁通口的要求，但是它只对磁通口进行了 准静态分析，并没有解决磁通口的瞬态分析问题。因此有必要设计一种磁通口瞬态分析方 法，满足耗时少、精度高的要求。
【发明内容】
[0004] 本发明所要解决的技术问题是克服现有技术的上述不足而提供一种基于磁场积 分方程法的磁通口瞬态分析方法，其能够对磁通口进行快速精确分析和优化。
[000引本发明解决上述问题所采用的技术方案为，基于磁场积分方程法的磁通口瞬态分 析方法，磁通口包括铁忍、激励线圈、接收线圈，铁忍为两条长条状薄片，磁通口铁忍数学模 型简化为平面二维模型，其特征在于，磁通口瞬态分析步骤包括： 首先计算磁通口铁忍的形状参数并测试铁忍材料的磁化曲线，建立准静态模型； 然后在准静态模型中引入激励电流频率的影响：在低频情况下可W忽略位移电流，电 流密度可分解为满流I:和给定的激励电流密度％，即!-?净I:，根据磁场叠加原理，考虑满 流影响后，磁通口铁忍内磁场强度1?可w分解成：
其中E是源磁场，包括激励电流产生的磁场强度和被测外磁场，11?是满流产生的 寄生磁场，?是铁忍材料产生的退磁场； 将满流产生的磁场强度移项到左端，并记堯红議哮;龜;，可得秦::二義束錢。裝自，铁忍内 磁场方程组为：
对于薄片状、圆柱状铁忍，可W用铁忍形状参数和磁感应强度时变率表示满流产生的 寄生磁场，其中铁忍材料数学模型??为：
最后建立磁通口瞬态分析模型， 在考虑频率的外磁场8中，铁忍材料内的磁场親与磁化强度夠的关系由W下方程组描 述：
其中3为形状参数，1?为磁滞损失参数，:S劝平均场参数，C为磁畴壁弯曲常数，为满 流产生的寄生磁场是与铁忍厚度和导电率有关的系数，为铁忍材料的饱和磁化 强度，鴻为不可逆磁化强度，該sb为非磁滞磁化强度，鋪为铁忍材料的磁化强度，B为磁感 应强度，Η为磁通口铁忍的总磁场，参数反。、：^、^1、％和^描述非磁滞磁化曲线的形状，线 性部分由Ra描述，而饱和部分由和、描述，S为方向参数，当地/dt〉0时裹取+1，当 地/化<0时瑟取-1; 联合铁忍材料的磁化曲线及W上公式、方程组，对磁通口进行瞬态分析。
[0006 ]建立准静态模型的步骤包括： 第一步，计算磁通口铁忍的形状参数并测试铁忍材料的磁化曲线； 第二步，建立磁场强度模型， 根据磁场叠加原理，磁通口铁忍任意一点的总磁场可W分为Ξ部分： 藉=:.裹;手變办難議..............................（1) 其中薇和顯为源场，囊为激励电流产生的激励磁场，蠢X为被测外磁场，龜为铁忍材 料产生的退磁场，揉为总磁场强度，旋度为零； 第Ξ步，对磁通口铁忍进行剖分，计算铁忍中每个剖分单元的磁化强度， 将磁通口铁忍剖分为η个单元，假设每个单元内的磁化强度是均匀的W及磁荷只在单 元之间的交界面和铁忍与空气的边界出现，可得：
其中Κ为第i个单元的总磁场强度，为激励电流在第i个单元产生的磁场，》棘为被 测外磁场在第i个单元产生的磁场，Sij为铁忍材料在第i个单元内产生的退磁场，苗为第i 个单元的磁化强度，下标X和y分别表示在铁忍平面X轴和y轴方向的磁场分量； 第四步，分析激励磁场强度， 根据毕奥-萨伐定律，激励电流在场点r处产生的激励磁场强度可W表述为：
其中为激励电流方向，i为激励电流密度，?为磁场的坐标，梦为电流的坐标，dV为激 励电流体积微元； 第五步，计算磁通口铁忍材料产生的退磁场， 根据：
退磁场旋度为零，存在一个标量满足:一聚藥，那么标量所满足的方程为：
铁忍材料所有剖分单元在第i个单元产生的磁场为：
其中霧为第j个单元对应的磁标势，r为源点和场点之间的距离，Sj为第j个单元的磁 化强度，苟为第j个单元的外表面，其中京为第i单元沿外表面法线单位矢量； 第六步，建立准静态模型， 将式(3)和式(6)代入式(2)可得：
其中雜?为退磁矩阵，将式(7)与磁通口铁忍材料的磁化曲线联合，建立磁通口的准静 态分析模型，求解在受到被测磁场和激励电流产生的激励磁场作用时铁忍内部的磁场和磁 化强度分布。
[0007]本发明的优点在于:采用本发明建立的磁通口瞬态分析模型不仅准确地综合考虑 了磁通口的结构、铁忍材料、激励电压、激励频率等因素的影响，而且准确地考虑了随着频 率增大而增强的磁感应强度滞后现象，从而能够准确的对磁通口进行瞬态分析，W便优化 磁通口的设计。
【附图说明】
[000引图1是磁通口结构示意图。
[0009] 图2是正弦电压激励的长条型磁通口电路。
[0010] 图3是被测磁场为30μΤ、激励信号为幅值为1.2V、激励频率为1曲Ζ的正弦电压时， 线圈电流仿真结果和实验结果比较图。
[0011] 图4是被测磁场为30μτ、激励信号为幅值为1.2V、激励频率为1曲Ζ的正弦电压时， 接收线圈输出电压的仿真结果和实验结果比较图。
[0012] 图5是激励频率为化化时的仿真磁滞回线和实验磁滞回线比较图。
[0013] 图6是被测磁场为30μΤ时和激励频率为化化时，随着激励电压增大，输出二次谐波 电压幅值与激励线圈电流幅值的关系。
【具体实施方式】
[0014] 下面结合附图、实施例对本发明进一步说明。
[0015] 实施例，如图1所示长条形磁通口，包括铁忍1、激励线圈2、接收线圈3,激励线圈2 和接收线圈3均匀绕制在铁忍1外，Ξ者组成磁通口。图2为正弦电压激励的长条型磁通口电 路，R为激励线圈的电阻。铁忍1为两条长条状薄片，厚度远小于其它尺寸，可W认为薄片内 部的磁场主要平行于其表面，垂直方向的磁化强度可近似忽略不计。磁通口接收线圈3中感 应到的是沿铁忍平面的磁场分量产生的磁通，而与之垂直的铁忍横截面内的磁场分量不会 对输出线圈电压产生影响，因此磁通口铁忍数学模型可W简化为平面二维模型。
[0016] 第一步，计算磁通口铁忍的形状参数并测试铁忍材料的磁化曲线 铁忍的形状参数、铁忍材料磁化曲线的计算测量方法已有现有技术公开，不再寶述。
[0017] 第二步，建立磁场强度模型 根据磁场叠加原理，对于平面二维模型，磁通口铁忍任意一点的总磁场可W分为Ξ部 分：
(1) 其中嚷和藏为源场，1|为激励电流产生的激励磁场，顯为被测外磁场，麵为铁忍材 料产生的退磁场，再为总磁场强度，旋度为零。
[0018] 第Ξ步，对磁通口铁忍进行剖分，计算铁忍中每个剖分单元的磁化强度 将磁通口铁忍剖分为η个单元，假设每个单元内的磁化强度是均匀的W及磁荷只在单 元之间的交界面和铁忍与空气的边界出现，可得：
其中K为第i个单元的总磁场强度，it为激励电流在第i个单元产生的磁场,??为被 测外磁场在第i个单元产生的磁场，为铁忍材料在第i个单元内产生的退磁场，瑪为第i 个单元的磁化强度，下标X和y分别表示在铁忍平面X轴和y轴方向的磁场分量。
[0019 ]第四步，分析激励磁场强度： 根据毕奥-萨伐定律，激励电流在场点r处产生的激励磁场强度可W表述为：
其中^为激励电流方向，i为激励电流密度J 为磁场的坐标，？f为电流的坐标，dV为激 励电流体积微元。
[0020] 第五步，计算磁通口铁忍材料产生的退磁场： 根据：
退磁场旋度为零，存在一个标量满足:1|:=心专缘，那么标量所满足的方程为：
铁忍材料所有剖分单元在第i个单元产生的磁场为：
其中嚷为第j个单元对应的磁标势，r为源点和场点之间的距离，为第j个单元的磁 化强度，覆为第j个单元的外表面，其中％.=敬i，京为第i单元沿外表面法线单位矢量。
[0021] 第六步，建立准静态模型 将式(3)和式(6)代入式(2)可得：
其中犧?为退磁矩阵。将式(7)与磁通口铁忍材料的磁化曲线联合，建立磁通口的准静 态分析模型，即可求解在受到被测磁场和激励电流产生的激励磁场作用时铁忍内部的磁场 和磁化强度分布，从而能够对磁通口进行准静态分析。
[0022] 第屯步，在准静态分析模型中引入激励电流频率的影响： 由于准静态分析模型是针对低频情况建立的模型，在低频情况下可W忽略位移电流， 电流密度可分解为满流拓和给定的激励电流密度％，即耗％料0。根据磁场叠加原理，考虑 满流影响后，磁通口铁忍内磁场强度寵（注：与前面用義表示铁忍内磁场不同，为方便起 见，用1表示铁忍内磁场)可W分解成：
其中ii。是源磁场，包括激励电流产生的磁场强度和被测外磁场。?暴S是满流产生的 寄生磁场，1%是铁忍材料产生的退磁场。
[0023] 将满流产生的磁场强度移项到左端，并记堯二％去赛a，可得藻二铁忍内 磁场方程组为：
对于薄片状、圆柱状铁忍，可W用铁忍形状参数和磁感应强度时变率表示满流产生的 寄生磁场，其中铁忍材料数学模型if!)为：
第八步，建立磁通口瞬态分析模型： 建立考虑频率的铁忍材料的数学模型，在考虑频率的外磁场採中，铁忍材料内的磁场 瑪与磁化强度級觀关系由W下方程组描述：
其中B为磁感应强度，3为形状参数，&为磁滞损失参数，α为平均场参数，;C为磁畴壁弯 曲常数。琴1为满流产生的寄生磁场晦iSyi，方程组的第一式可转换为i二錢S -鬚#歡，犹是 与铁忍厚度和导电率有关的系数。M。为铁忍材料的饱和磁化强度，絳《为不可逆磁化强度， 擁《^为非磁滞磁化强度，錄为铁忍材料的磁化强度。参数&。、而、麵、游和呜描述非磁滞磁 化曲线的形状。线性部分由觀描述，而饱和部分由觀巧鸣描述。S为方向参数，当 dH/dt〉0时瑟取+1，当dH/dt<0时送取-1，运样保证了不会出现不符合实际的总能量增加的现 象。
[0024] 联合铁忍材料的磁化曲线、式（10) W及式（11)所示的考虑频率的铁忍材料数学模 型，即可建立磁通口瞬态分析模型，能够准确考虑满流对铁忍材料的影响，从而能够对磁通 口输出信号进行瞬态分析。
[0025] W图1所示的磁通口采用瞬态分析模型进行验证。激励线圈2的应数为330,接收线 圈3的应数为1100。铁忍1为长度30mm，宽度1mm，厚度0.1mm的坡莫合金。铁忍材料由式(9)描 述，其模型参数为=1.0694 X 1〇-6，c=0.3，k=l. 7143，Ms=5.98 X 10日，Ru=0.3997，册=0.0001， Hl=1.1263Xl〇-2〇，aO=0.0827，al=0.0201，N=0.0014。图 2为正弦电压激励的长条型磁通口 电路。R为激励线圈的电阻，Κ=13.3Ω。
[0026] 采用正弦电压信号对磁通口进行激励，比较考虑频率的基于磁场积分方程法的磁 通口模型的仿真结果和实验结果，其中仿真结果是结合分析模型并通过仿真计算软件得 出，实验结果是通过建立的磁通口测量系统得到的，其中信号发生器和功率放大器产生正 弦电压激励信号，直流电源产生被测磁场，示波器用来检测磁通口输出信号。图3和图4分别 表示被测磁场为30μΤ、激励信号为幅值为1.2V和激励频率为1曲Ζ的正弦电压激励线圈电流 W及接收线圈输出电压的仿真结果和实验结果比较图。图5表示激励频率为Ik化时的仿真 磁滞回线和实验磁滞回线比较图。图6表示被测磁场为30μΤ时和激励频率为1曲Z时，随着激 励电压增大，输出二次谐波电压幅值与激励线圈电流幅值的关系。
[0027] 从图3和图4可W看出，因为实际中的激励电流需要一定的响应时间，在开始的一 段时间内，激励电流和输出信号的实验值滞后于仿真模型中的激励电流和输出信号。在后 续的时间条件下，仿真中的激励电流和输出信号和与实验中的激励电流和输出信号吻合良 好。从图5可W看出，仿真磁滞回线与实验测得的磁滞回线几乎一致。从图6可W看出，在小 于95mA的激励电流条件下，由仿真模型得到磁通口二次谐波幅值略低于实验测得的磁通口 二次谐波幅值;在大于95mA的激励电流条件下，由仿真模型得到磁通口二次谐波幅值略高 于实验测得的磁通口二次谐波幅值。从仿真结果与实验结果的对比可W看出，随着激励电 压，激励频率的改变，输出结果准确地反映了激励电流，输出电压等的变化，表明所建立的 磁通口瞬态分析模型可W准确的对磁通口输出信号进行瞬态分析。
【主权项】
1.基于磁场积分方程法的磁通门瞬态分析方法，磁通门包括铁芯、激励线圈、接收线 圈，铁芯为两条长条状薄片，磁通门铁芯数学模型简化为平面二维模型，其特征在于，磁通 门瞬态分析包括以下步骤： 首先计算磁通门铁芯的形状参数并测试铁芯材料的磁化曲线，建立准静态模型； 然后在准静态模型中引入激励电流频率的影响：在低频情况下可以忽略位移电流，电 流密度可分解为涡流U和给定的激励电流密度1￡5，即【产％1知，根据磁场叠加原理，考虑涡 流影响后，磁通门铁芯内磁场强度||可以分解成：其中!^是源磁场，包括激励电流产生的磁场强度和被测外磁场，.是涡流产生的寄 生磁场，是铁芯材料产生的退磁场； 将涡流产生的磁场强度移项到左端，并记铁芯内磁 场方程组为：对于薄片状、圆柱状铁芯，可以用铁芯形状参数和磁感应强度时变率表示涡流产生的 寄生磁场，其中铁芯材料数学模型義I?为：最后建立磁通门瞬态分析模型， 在考虑频率的外磁场Il中，铁芯材料内的磁场与磁化强度麵的关系由以下方程组描 述：其中獅状参数，咖滞损失参数，痛均场参数，畴壁弯曲常数， 流产生的寄生磁场1??，，辨是与铁芯厚度和导电率有关的系数，_为铁芯材料的饱和磁化 强度，Mirr为不可逆磁化强度,?涵为非磁滞磁化强度，議为铁芯材料的磁化强度，B为磁感 应强度，H为磁通门铁芯的总磁场，参数、:11、_、％和％描述非磁滞磁化曲线的形状，线 性部分由装。描述，而饱和部分由_、:_、和描述，§为方向参数，当dH/dt>0时S取+1，当 dH/dt〈0 时S 取-1; 联合铁芯材料的磁化曲线及以上公式、方程组，对磁通门进行瞬态分析。2.根据权利要求1所述的基于磁场积分方程法的磁通门瞬态分析方法，其特征在于，建 立准静态模型的步骤包括： 第一步，计算磁通门铁芯的形状参数并测试铁芯材料的磁化曲线； 第二步，建立磁场强度模型， 根据磁场叠加原理，磁通门铁芯任意一点的总磁场可以分为三部分：其中和1|^为源场，为激励电流产生的激励磁场，为被测外磁场，_为铁芯材 料产生的退磁场，M为总磁场强度，旋度为零； 第三步，对磁通门铁芯进行剖分，计算铁芯中每个剖分单元的磁化强度， 将磁通门铁芯剖分为η个单元，假设每个单元内的磁化强度是均匀的以及磁荷只在单 元之间的交界面和铁芯与空气的边界出现，可得：其中|丨为第i个单元的总磁场强度,?为激励电流在第i个单元产生的磁场，为被测 外磁场在第i个单元产生的磁场，1?...为铁芯材料在第i个单元内产生的退磁场，:轉为第i个 ·'；.￡ ... 单元的磁化强度，下标X和y分别表示在铁芯平面X轴和y轴方向的磁场分量； 第四步，分析激励磁场强度， 根据毕奥-萨伐定律，激励电流在场点r处产生的激励磁场强度可以表述为：其中_为激励电流方向I为激励电流密度J为磁场的坐标Ji为电流的坐标，dV为激励 电流体积微元； 第五步，计算磁通门铁芯材料产生的退磁场， 根据：I 退磁场旋度为零，存在一个标量满足：那么标量所满足的方程为：铁芯材料所有剖分单元在第i个单元产生的磁场为：其中％为第j个单元对应的磁标势，r为源点和场点之间的距离，为第j个单元的磁化 强度，sI为第j个单元的外表面，其中为第i单元沿外表面法线单位矢量； 第六步，建立准静态模型， 将式(3)和式(6)代入式(2)可得：其中辨k为退磁矩阵，将式(7)与磁通门铁芯材料的磁化曲线联合，建立磁通门的准静 态分析模型，求解在受到被测磁场和激励电流产生的激励磁场作用时铁芯内部的磁场和磁 化强度分布。
【专利摘要】本发明公开了一种基于磁场积分方程法的磁通门瞬态分析方法，根据磁通门铁芯特点，建立基于磁场积分方程法的磁通门准静态分析模型，并引入激励电流频率的影响，建立可以对磁通门进行瞬态分析的数学模型。采用本发明建立的磁通门瞬态分析模型不仅准确地综合考虑了磁通门的结构、铁芯材料、激励电压、激励频率等因素的影响，而且准确地考虑了随着频率增大而增强的磁感应强度滞后现象，从而能够准确的对磁通门进行瞬态分析。
【IPC分类】G01R33/00, G01R33/04
【公开号】CN105487025
【申请号】CN201510842978
【发明人】侯晓伟, 李菊萍, 郭俊杰, 倪大成, 王飞, 郑华雄, 郑良广
【申请人】宁波南车时代传感技术有限公司
【公开日】2016年4月13日
【申请日】2015年11月27日
最新回复(0)