根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法
【专利摘要】根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法，属于航天航空领域，本发明为解决由静态磁场来确定霍尔推力器的线圈安匝变化百分率范围是不准确的问题。本发明方法：步骤一、将霍尔推力器的二维对称模型导入FEMM中，建立磁路模型，步骤二、三个线圈模拟通入初始电流产生静态磁场，获取静态磁场零坐标位置；步骤三、以20％为步长，逐渐改变其中一个线圈通入的电流值，获取线圈安匝变化百分率与零磁场位置变化百分率对应曲线图；步骤四、根据动态磁场零坐标位置偏离静态磁场零坐标位置的距离应小于通道特征尺寸的2％～2.5％的规定，及步骤三实现确定外线圈安匝变化百分率的范围。
【专利说明】根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法
【技术领域】
[0001]本发明涉及确定霍尔推力器运行参数的方法，属于航天航空领域。
【背景技术】
[0002]电推进装置以其高效率高比冲的优点取代化学推力器已经成为航天推进领域发展的一种趋势。霍尔推力器是电推进装置的一种，以其效率高、工作寿命长、功率密度高、t匕冲适中等优点成为卫星、探测器等航天飞行器的重要动力装置。霍尔推力器通过霍尔效应产生推力，磁场是产生霍尔效应的关键。因此，霍尔推力器必须设计出适当的磁场来提高电离率，增加等离子体密度，并有效地约束等离子体的行为。霍尔推力器的磁场通常是在静态磁路情况下设计的，利用适当的磁路结构和内线圈，外线圈以及附加线圈产生适当的出口区磁场强度、正梯度分布、凸向阳极的弯曲磁力线等优化的磁场位形。
[0003]由霍尔推力器工作原理可知，霍尔漂移电流作为霍尔推力器中一个重要的基础物理过程是必然存在的，同时闭环霍尔漂移电流存在宽谱的振荡现象，根据法拉第电磁感应定律可知，周向的霍尔漂移电流振荡感应出沿轴向随时间变化磁场，时变的磁场引起沿周向绕制励磁线圈中感生出波动的感应电动势，从而引起励磁线圈电流的波动，励磁电流波动耦合到放电回路中，引起放电电流在励磁线圈固有频率处的振荡现象，励磁电流振荡和放电电流振荡相互耦合，达到动态平衡，这就是电磁耦合振荡。电磁耦合振荡使推力器磁场在放电过程中处于波动状态，导致推力器实际放电过程中，通道内的磁场形貌和强度并不是静态原则设计出的优化磁场位形。因此，由静态磁场确定的线圈安匝变化百分率范围是不准确的。
【发明内容】
[0004]本发明目的是为了解决由静态磁场来确定霍尔推力器的线圈安匝变化百分率范围是不准确的问题，提供了一种根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法。
[0005]本发明包括三个方案。
[0006]第一个方案:根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法，所述霍尔推力器的二维对称模型包括外线圈、内线圈、附加线圈、内磁极、阳极、底板，其具有对称轴线，
[0007]线圈安匝变化百分率范围为外线圈的安匝变化百分率范围、内线圈的安匝变化百分率范围或附加线圈的安匝变化百分率范围；
[0008]确定外线圈的安匝变化百分率范围方法包括以下步骤:
[0009]步骤一、将霍尔推力器的二维对称模型导入电磁场有限元分析软件FEMM中，建立霍尔推力器的磁路模型，建立rz坐标系:以内磁极和底板的交点为原点坐标；以对称轴线为径向坐标z轴，以底板所在直线为轴向坐标r轴；[0010]步骤二、在电磁场有限元分析软件FEMM中，外线圈、内线圈和附加线圈模拟通入初始电流产生静态磁场，获取静态磁场零坐标位置OtlCrtl, z0)；
[0011]步骤三、以20%为步长，逐渐改变外线圈通入的电流值，改变外线圈通入的电流值时会产生一个新的磁场，该磁场为线圈稱合振荡引起的动态磁场与静态磁场的合成磁场，记录每次外线圈通入的电流值与对应的合成磁场的零磁场坐标O (r, Z)；
[0012]进而获取外线圈安匝变化百分率与零磁场位置变化百分率对应曲线图；
[0013]其中:外线圈安匝为外线圈通入的电流乘以外线圈的匝数；外线圈安匝变化百分率指每次外线圈通入电流后的外线圈安匝相对于外线圈通入初始电流的变化百分率；零磁场位置变化百分率指每次外线圈通入电流后的磁场零磁场坐标0(r，z)相对于静态磁场零坐标位置OtlO^ Z0)的变化百分率，包括轴向位置变化百分率和径向位置变化百分率；
[0014]步骤四、根据动态磁场零坐标位置0(r, z)偏尚静态磁场零坐标位置OtlCrtl, Ztl)的距离应小于通道特征尺寸的2%~2.5%的规定，及步骤三获取的外线圈安匝变化百分率与零磁场位置变化百分率对应曲线图，实现确定外线圈安匝变化百分率的范围。
[0015]第二个方案:根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法，所述霍尔推力器的二维对称模型包括外线圈、内线圈、附加线圈、内磁极、阳极、底板，其具有对称轴线，
[0016]线圈安匝变化百分率范围为外线圈的安匝变化百分率范围、内线圈的安匝变化百分率范围或附加线圈的安 匝变化百分率范围；
[0017]确定内线圈的安匝变化百分率范围方法包括以下步骤:
[0018]步骤一、将霍尔推力器的二维对称模型导入电磁场有限元分析软件FEMM中，建立霍尔推力器的磁路模型，建立rz坐标系:以内磁极和底板的交点为原点坐标；以对称轴线为径向坐标z轴，以底板所在直线为轴向坐标r轴；
[0019]步骤二、在电磁场有限元分析软件FEMM中，外线圈、内线圈和附加线圈模拟通入初始电流产生静态磁场，获取静态磁场零坐标位置OtlCrtl, z0)；
[0020]步骤三、以20%为步长，逐渐改变内线圈通入的电流值，改变内线圈通入的电流值时会产生一个新的磁场，该磁场为线圈稱合振荡引起的动态磁场与静态磁场的合成磁场，记录每次内线圈通入的电流值与对应的合成磁场的零磁场坐标0(r, z)；
[0021]进而获取内线圈安匝变化百分率与零磁场位置变化百分率对应曲线图；
[0022]其中:内线圈安匝为内线圈通入的电流乘以内线圈的匝数；内线圈安匝变化百分率指每次内线圈通入电流后的内线圈安匝相对于内线圈通入初始电流的变化百分率；零磁场位置变化百分率指每次内线圈通入电流后的磁场零磁场坐标0(r，z)相对于静态磁场零坐标位置OtlO^ Z0)的变化百分率，包括轴向位置变化百分率和径向位置变化百分率；
[0023]步骤四、根据动态磁场零坐标位置0(r, z)偏尚静态磁场零坐标位置OtlCrtl, Ztl)的距离应小于通道特征尺寸的2%~2.5%的规定，及步骤三获取的内线圈安匝变化百分率与零磁场位置变化百分率对应曲线图，实现确定内线圈安匝变化百分率的范围。
[0024]第三个方案:根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法，所述霍尔推力器的二维对称模型包括外线圈、内线圈、附加线圈、内磁极、阳极、底板，其具有对称轴线，
[0025]线圈安匝变化百分率范围为外线圈的安匝变化百分率范围、内线圈的安匝变化百分率范围或附加线圈的安匝变化百分率范围；
[0026]确定附加线圈的安匝变化百分率范围方法包括以下步骤:
[0027]步骤一、将霍尔推力器的二维对称模型导入电磁场有限元分析软件FEMM中，建立霍尔推力器的磁路模型，建立rz坐标系:以内磁极和底板的交点为原点坐标；以对称轴线为径向坐标z轴，以底板所在直线为轴向坐标r轴；
[0028]步骤二、在电磁场有限元分析软件FEMM中，外线圈、内线圈和附加线圈模拟通入初始电流产生静态磁场，获取静态磁场零坐标位置OtlCrtl, z0)；
[0029]步骤三、以20%为步长，逐渐改变附加线圈通入的电流值，改变附加线圈通入的电流值时会产生一个新的磁场，该磁场为线圈稱合振荡引起的动态磁场与静态磁场的合成磁场，记录每次附加线圈通入的电流值与对应的合成磁场的零磁场坐标0(r, z)；
[0030]进而获取附加线圈安匝变化百分率与零磁场位置变化百分率对应曲线图；
[0031]其中:附加线圈安匝为附加线圈通入的电流乘以附加线圈的匝数；附加线圈安匝变化百分率指每次附加线圈通入电流后的外线圈安匝相对于附加线圈通入初始电流的变化百分率；零磁场位置变化百分率指每次附加线圈通入电流后的磁场零磁场坐标0(r，z)相对于静态磁场零坐标位置Otl0-(ι，ζ0)的变化百分率，包括轴向位置变化百分率和径向位置变化百分率；
[0032]步骤四、根据动态磁场零坐标位置0(r, ζ)偏尚静态磁场零坐标位置OtlCrtl, Ztl)的距离应小于通道特征尺寸的2%~2.5%的规定，及步骤三获取的附加线圈安匝变化百分率与零磁场位置变化百分率对应曲线图，实现确定附加线圈安匝变化百分率的范围。
[0033]本发明的优点:本发明方法有效地分析振荡的励磁电流对通道内磁场设计的影响，克服静态磁场设计应用的局限性，合理的给出了考虑励磁/放电耦合振荡引起的动态磁场应对静态磁场的附加设计原则，给出准确的霍尔推力器的线圈安匝变化百分率范围。
【专利附图】
【附图说明】
[0034]图1是本发明方法涉及霍尔推力器的结构示意图；
[0035]图2是实施方式一所述根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法的流程图，此图为确定外线圈安匝变化百分率范围的方法；
[0036]图3是实施方式二所述根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法的流程图，此图为确定内线圈安匝变化百分率范围的方法；；
[0037]图4是实施方式三所述根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法的流程图，此图为确定附加线圈安匝变化百分率范围的方法；；
[0038]图5是霍尔推力器线圈通入初始电流，由FEMM仿真得到的静态磁场位形图；
[0039]图6是外线圈安阻变化百分率与零磁场位置变化百分率对应曲线图；
[0040]图7是内线圈安阻变化百分率与零磁场位置变化百分率对应曲线图；
[0041]图8是附加线圈安匝变化百分率与零磁场位置变化百分率对应曲线图。
【具体实施方式】
[0042]【具体实施方式】一:下面结合图1、图2、图5和图6说明本实施方式，本实施方式所述根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法，所述霍尔推力器的二维对称模型包括外线圈1、内线圈2、附加线圈3、内磁极4、阳极5、底板7，其具有对称轴线6，
[0043]线圈安匝变化百分率范围为外线圈I的安匝变化百分率范围、内线圈2的安匝变化百分率范围或附加线圈3的安匝变化百分率范围；
[0044]确定外线圈I的安匝变化百分率范围方法包括以下步骤:
[0045]步骤一、将霍尔推力器的二维对称模型导入电磁场有限元分析软件FEMM中，建立霍尔推力器的磁路模型，建立rz坐标系:以内磁极4和底板7的交点为原点坐标；以对称轴线6为径向坐标ζ轴，以底板7所在直线为轴向坐标r轴；
[0046]步骤二、在电磁场有限元分析软件FEMM中，外线圈1、内线圈2和附加线圈3模拟通入初始电流产生静态磁场，获取静态磁场零坐标位置OtlO^ Z0)；
[0047]步骤三、以20%为步长，逐渐改变外线圈I通入的电流值，改变外线圈I通入的电流值时会产生一个新的磁场，该磁场为线圈稱合振荡引起的动态磁场与静态磁场的合成磁场，记录每次外线圈I通入的电流值与对应的合成磁场的零磁场坐标O(r, Z)；
[0048]进而获取外线圈安匝变化百分率与零磁场位置变化百分率对应曲线图；
[0049]其中:外线圈安匝为外线圈I通入的电流乘以外线圈I的匝数；外线圈安匝变化百分率指每次外线圈I通入电流后的外线圈安匝相对于外线圈I通入初始电流的变化百分率；零磁场位置变化百分率指每次外线圈I通入电流后的磁场零磁场坐标0(r，ζ)相对于静态磁场零坐标位置Otl O^ ζ0)的变化百分率，包括轴向位置变化百分率和径向位置变化百分率；
[0050]步骤四、根据动态磁场零坐标位置0(r, ζ)偏尚静态磁场零坐标位置OtlCrtl, Ztl)的距离应小于通道特征尺寸的2%~2.5%的规定，及步骤三获取的外线圈安匝变化百分率与零磁场位置变化百分率对应曲线图，实现确定外线圈安匝变化百分率的范围。
[0051]改变外线圈I通入电流的范围为外线圈I通入初始电流的-60%~80%。
[0052]用安匝等效原则，在线圈匝数的前提下，改变线圈励磁电流。
[0053]外线圈1、内线圈2和附加线圈3模拟通入初始电流的值可以选择在霍尔推力器最优工作状态下的三个线圈通入的电流值。
[0054]步骤四中的通道特征尺寸是指霍尔推力器的固有参数，表示霍尔推力器外径，t匕如型号为HEP70的霍尔推力器，它的通道特征尺寸即为70mm。
[0055]在最优放电状态下的静态磁场如图5所示，其静态磁场零坐标位置OtlO^ Z0)，当改变外线圈I通入电流后，磁场位形会发生变化，磁场零坐标的位置也会随之改变，但它变化的范围是有规定的:动态磁场零坐标位置O(r, Z)偏离静态磁场零坐标位置OtlCrtl, Zci)的距离应小于通道特征尺寸的2%~2.5%。记录下每次外线圈I通入的电流值与对应的动态磁场的零磁场坐标0(r，ζ);进而获取外线圈安匝变化百分率与零磁场位置变化百分率对应曲线图；如图6所示，则根据这个对应关系，就能确定外线圈安匝变化百分率的范围。
[0056]【具体实施方式】二:下面结合图1、图3、图5和图7说明本实施方式，本实施方式所述霍尔推力器的二维对称模型包括外线圈1、内线圈2、附加线圈3、内磁极4、阳极5、底板7，其具有对称轴线6，
[0057]线圈安匝变化百分率范围为外线圈I的安匝变化百分率范围、内线圈2的安匝变化百分率范围或附加线圈3的安匝变化百分率范围；[0058]确定内线圈2的安匝变化百分率范围方法包括以下步骤:
[0059]步骤一、将霍尔推力器的二维对称模型导入电磁场有限元分析软件FEMM中，建立霍尔推力器的磁路模型，建立rz坐标系:以内磁极4和底板7的交点为原点坐标；以对称轴线6为径向坐标ζ轴，以底板7所在直线为轴向坐标r轴；
[0060]步骤二、在电磁场有限元分析软件FEMM中，外线圈1、内线圈2和附加线圈3模拟通入初始电流产生静态磁场，获取静态磁场零坐标位置OtlO^ Z0)；
[0061]步骤三、以20%为步长，逐渐改变内线圈2通入的电流值，改变内线圈2通入的电流值时会产生一个新的磁场，该磁场为线圈稱合振荡引起的动态磁场与静态磁场的合成磁场，记录每次内线圈2通入的电流值与对应的合成磁场的零磁场坐标0(r, ζ)；
[0062]进而获取内线圈安匝变化百分率与零磁场位置变化百分率对应曲线图；
[0063]其中:内线圈安匝为内线圈2通入的电流乘以内线圈2的匝数；内线圈安匝变化百分率指每次内线圈2通入电流后的内线圈安匝相对于内线圈2通入初始电流的变化百分率；零磁场位置变化百分率指每次内线圈2通入电流后的磁场零磁场坐标0(r，ζ)相对于静态磁场零坐标位置OtlO^ ζ0)的变化百分率，包括轴向位置变化百分率和径向位置变化百分率；
[0064]步骤四、根据动态磁场零坐标位置0(r, ζ)偏尚静态磁场零坐标位置OtlCrtl, Ztl)的距离应小于通道特征尺寸的2%~2.5%的规定，及步骤三获取的内线圈安匝变化百分率与零磁场位置变化百分 率对应曲线图，实现确定内线圈安匝变化百分率的范围。
[0065]改变内线圈2通入电流的范围为内线圈2通入初始电流的-60%~80%。
[0066]工作原理与实施方式一相同，不再赘述。
[0067]【具体实施方式】三:下面结合图1、图3、图5和图8说明本实施方式，本实施方式所述根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法，所述霍尔推力器的二维对称模型包括外线圈1、内线圈2、附加线圈3、内磁极4、阳极5、底板7，其具有对称轴线6，
[0068]线圈安匝变化百分率范围为外线圈I的安匝变化百分率范围、内线圈2的安匝变化百分率范围或附加线圈3的安匝变化百分率范围；
[0069]确定附加线圈3的安匝变化百分率范围方法包括以下步骤:
[0070]步骤一、将霍尔推力器的二维对称模型导入电磁场有限元分析软件FEMM中，建立霍尔推力器的磁路模型，建立rz坐标系:以内磁极4和底板7的交点为原点坐标；以对称轴线6为径向坐标ζ轴，以底板7所在直线为轴向坐标r轴；
[0071]步骤二、在电磁场有限元分析软件FEMM中，外线圈1、内线圈2和附加线圈3模拟通入初始电流产生静态磁场，获取静态磁场零坐标位置OtlO^ Z0)；
[0072]步骤三、以20%为步长，逐渐改变附加线圈3通入的电流值，改变附加线圈3通入的电流值时会产生一个新的磁场，该磁场为线圈稱合振荡引起的动态磁场与静态磁场的合成磁场，记录每次附加线圈3通入的电流值与对应的合成磁场的零磁场坐标0(r, ζ)；
[0073]进而获取附加线圈安匝变化百分率与零磁场位置变化百分率对应曲线图；
[0074]其中:附加线圈安匝为附加线圈3通入的电流乘以附加线圈3的匝数；附加线圈安匝变化百分率指每次附加线圈3通入电流后的外线圈安匝相对于附加线圈3通入初始电流的变化百分率；零磁场位置变化百分率指每次附加线圈3通入电流后的磁场零磁场坐标0(r, z)相对于静态磁场零坐标位置OtlO^ ζ0)的变化百分率，包括轴向位置变化百分率和径向位置变化百分率；
[0075] 步骤四、根据动态磁场零坐标位置0(r，ζ)偏离静态磁场零坐标位置OtlO^ Z0)的距离应小于通道特征尺寸的2%~2.5%的规定，及步骤三获取的附加线圈安匝变化百分率与零磁场位置变化百分率对应曲线图，实现确定附加线圈安匝变化百分率的范围。
[0076]改变附加线圈3通入电流的范围为附加线圈3初始电流的-60%~80%。
[0077]工作原理与实施方式一相同，不再赘述。
【权利要求】
1.根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法，所述霍尔推力器的二维对称模型包括外线圈(I)、内线圈⑵、附加线圈(3)、内磁极(4)、阳极(5)、底板(7)，其具有对称轴线(6)， 其特征在于，线圈安匝变化百分率范围为外线圈(I)的安匝变化百分率范围、内线圈(2)的安匝变化百分率范围或附加线圈(3)的安匝变化百分率范围； 确定外线圈(I)的安匝变化百分率范围方法包括以下步骤: 步骤一、将霍尔推力器的二维对称模型导入电磁场有限元分析软件FEMM中，建立霍尔推力器的磁路模型，建立rz坐标系:以内磁极(4)和底板(7)的交点为原点坐标；以对称轴线(6)为径向坐标z轴，以底板(7)所在直线为轴向坐标r轴； 步骤二、在电磁场有限元分析软件FEMM中，外线圈(I)、内线圈(2)和附加线圈(3)模拟通入初始电流产生静态磁场，获取静态磁场零坐标位置OtlO^ z0)； 步骤三、以20%为步长，逐渐改变外线圈(I)通入的电流值，改变外线圈(I)通入的电流值时会产生一个新的磁场，该磁场为线圈稱合振荡引起的动态磁场与静态磁场的合成磁场，记录每次外线圈(I)通入的电流值与对应的合成磁场的零磁场坐标0(r，z)； 进而获取外线圈安匝变化百分率与零磁场位置变化百分率对应曲线图； 其中:外线圈安匝为外线圈(I)通入的电流乘以外线圈(I)的匝数；外线圈安匝变化百分率指每次外线圈(I)通入电流后的外线圈安匝相对于外线圈(I)通入初始电流的变化百分率；零磁场位置变化百分率指每次外线圈(I)通入电流后的磁场零磁场坐标0(r，z)相对于静态磁场零坐标位置OtlO^ z0)的变化百分率，包括轴向位置变化百分率和径向位置变化百分率； 步骤四、根据动态磁场零坐标位置0(r，z)偏离静态磁场零坐标位置OtlO^ Ztl)的距离应小于通道特征尺寸的2%~2.5%的规定，及步骤三获取的外线圈安匝变化百分率与零磁场位置变化百分率对应曲线图，实现确定外线圈安匝变化百分率的范围。
2.根据权利要求1所述根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法，其特征在于，改变外线圈(I)通入电流的范围为外线圈(I)通入初始电流的-60%~80%。
3.根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法，所述霍尔推力器的二维对称模型包括外线圈(I)、内线圈⑵、附加线圈(3)、内磁极(4)、阳极(5)、底板(7)，其具有对称轴线(6)， 其特征在于，线圈安匝变化百分率范围为外线圈(I)的安匝变化百分率范围、内线圈(2)的安匝变化百分率范围或附加线圈(3)的安匝变化百分率范围； 确定内线圈(2)的安匝变化百分率范围方法包括以下步骤: 步骤一、将霍尔推力器的二维对称模型导入电磁场有限元分析软件FEMM中，建立霍尔推力器的磁路模 型，建立rz坐标系:以内磁极(4)和底板(7)的交点为原点坐标；以对称轴线(6)为径向坐标z轴，以底板(7)所在直线为轴向坐标r轴； 步骤二、在电磁场有限元分析软件FEMM中，外线圈(I)、内线圈(2)和附加线圈(3)模拟通入初始电流产生静态磁场，获取静态磁场零坐标位置OtlO^ z0)； 步骤三、以20%为步长，逐渐改变内线圈⑵通入的电流值，改变内线圈(2)通入的电流值时会产生一个新的磁场，该磁场为线圈稱合振荡引起的动态磁场与静态磁场的合成磁场，记录每次内线圈(2)通入的电流值与对应的合成磁场的零磁场坐标O(r, z)； 进而获取内线圈安匝变化百分率与零磁场位置变化百分率对应曲线图； 其中:内线圈安匝为内线圈(2)通入的电流乘以内线圈(2)的匝数；内线圈安匝变化百分率指每次内线圈(2)通入电流后的内线圈安匝相对于内线圈(2)通入初始电流的变化百分率；零磁场位置变化百分率指每次内线圈(2)通入电流后的磁场零磁场坐标0(r，z)相对于静态磁场零坐标位置OtlO^ z0)的变化百分率，包括轴向位置变化百分率和径向位置变化百分率； 步骤四、根据动态磁场零坐标位置0(r，z)偏离静态磁场零坐标位置OtlO^ Ztl)的距离应小于通道特征尺寸的2%~2.5%的规定，及步骤三获取的内线圈安匝变化百分率与零磁场位置变化百分率对应曲线图，实现确定内线圈安匝变化百分率的范围。
4.根据权利要求3所述根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法，其特征在于，改变内线圈(2)通入电流的范围为内线圈(2)通入初始电流的-60%~80%。
5.根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法，所述霍尔推力器的二维对称模型包括外线圈(I)、内线圈⑵、附加线圈(3)、内磁极(4)、阳极(5)、底板(7)，其具有对称轴线(6)， 其特征在于，线圈安匝变化百分率范围为外线圈(I)的安匝变化百分率范围、内线圈(2)的安匝变化百分率范围或附加线圈(3)的安匝变化百分率范围； 确定附加线圈(3)的安匝变化百分率范围方法包括以下步骤: 步骤一、将霍尔推力器的二维对称模型导入电磁场有限元分析软件FEMM中，建立霍尔推力器的磁路模型，建立rz坐标系:以内磁极(4)和底板(7)的交点为原点坐标；以对称轴线(6)为径向坐标z轴，以底板(7)所在直线为轴向坐标r轴； 步骤二、在电磁场有限元分析软件FEMM中，外线圈(I)、内线圈(2)和附加线圈(3)模拟通入初始电流产生静态磁场，获取静态磁场零坐标位置OtlO^ z0)； 步骤三、以20%为步长，逐渐改变附加线圈(3)通入的电流值，改变附加线圈(3)通入的电流值时会产生一个新的磁场，该磁场为线圈稱合振荡引起的动态磁场与静态磁场的合成磁场，记录每次附加线圈(3)通入的电流值与对应的合成磁场的零磁场坐标O(r, Z)； 进而获取附加线圈安匝变化百分率与零磁场位置变化百分率对应曲线图； 其中:附加线圈安匝为附加线圈(3)通入的电流乘以附加线圈(3)的匝数；附加线圈安匝变化百分率指每次附加线圈(3)通入电流后的外线圈安匝相对于附加线圈(3)通入初始电流的变化百分率；零磁场位置变化百分率指每次附加线圈(3)通入电流后的磁场零磁场坐标O (r, z)相对于静态磁场零坐标位置Otl (r0, z0)的变化百分率，包括轴向位置变化百分率和径向位置变化百分率； 步骤四、根据 动态磁场零坐标位置0(r，z)偏离静态磁场零坐标位置OtlO^ Ztl)的距离应小于通道特征尺寸的2%~2.5%的规定，及步骤三获取的附加线圈安匝变化百分率与零磁场位置变化百分率对应曲线图，实现确定附加线圈安匝变化百分率的范围。
6.根据权利要求5所述根据霍尔推力器中耦合振荡伴生的动态磁场确定线圈安匝变化百分率范围的方法，其特征在于，改变附加线圈(3)通入电流的范围为附加线圈(3)初始电流的-60%~80%。
【文档编号】G01R33/12GK103983927SQ201410257769
【公开日】2014年8月13日 申请日期:2014年6月11日 优先权日:2014年6月11日
【发明者】杨子怡, 魏立秋, 韩亮, 于达仁 申请人:哈尔滨工业大学