CN105844095B - 一种基于性能退化的离子推力器小子样可靠性评估方法 - Google Patents

一种基于性能退化的离子推力器小子样可靠性评估方法 Download PDF

Info

Publication number
CN105844095B
CN105844095B CN201610165279.8A CN201610165279A CN105844095B CN 105844095 B CN105844095 B CN 105844095B CN 201610165279 A CN201610165279 A CN 201610165279A CN 105844095 B CN105844095 B CN 105844095B
Authority
CN
China
Prior art keywords
ion thruster
performance degradation
performance
time
reliability
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610165279.8A
Other languages
English (en)
Other versions
CN105844095A (zh
Inventor
李军星
王治华
张勇波
傅惠民
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ningbo Zhizheng Weiying Information Technology Co ltd
Original Assignee
Beihang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN201610165279.8A priority Critical patent/CN105844095B/zh
Publication of CN105844095A publication Critical patent/CN105844095A/zh
Application granted granted Critical
Publication of CN105844095B publication Critical patent/CN105844095B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass

Landscapes

  • Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)

Abstract

一种基于性能退化的离子推力器小子样可靠性评估方法,包括以下步骤:一:对离子推力器的结构特点和关键失效模式进行分析,选取影响其寿命的关键性能参数;二:根据关键性能参数y随时间的退化规律,建立离子推力器小子样性能退化模型;三:确定离子推力器平均失效时间tMTTF和可靠度函数R(t);四:对第i台离子推力器,记录在时间tij时的测试性能退化值yij,i=1,2,...,m,j=1,2,...,ni,m为投入试验的离子推力器台数,ni为第i台离子推力器性能测试的总次数;五:估计性能退化模型中的未知参数;六:利用可靠度函数对离子推力器进行小子样可靠性评估。本发明开发利用不同时刻性能测试数据间的纵向信息,使信息量大幅增加。提高分析精度,同时在精度相同的情况下,则可节省大量试样。

Description

一种基于性能退化的离子推力器小子样可靠性评估方法
【技术领域】
本发明涉及一种基于性能退化的离子推力器极小子样可靠性分析方法,尤其涉及离子推力器单机寿命试验且试验数据为无失效的可靠性分析方法。属于航天可靠性分析技术领域。
【背景技术】
电推进是未来十大尖端技术之一,与传统化学推力器相比,具有高比冲,高效率,长寿命,推力小等特点。可以用来执行南北位保、轨道提升及深空探测等空间任务。离子推力器是电推进的一种,由于其推力比较小,这就要求必须运行较长的时间才能达到总冲量的要求。因此,离子推力器作为长寿命设备一般要求能够在轨可靠运行数千小时甚至上万小时。所以,对离子推力器服役期间的动态可靠性进行实时分析具有重要意义。
离子推力器子系统的基本组件包括空心阴极、放电室、栅极系统和中和器等,如图1所示。在离子源内由空心阴极发射的电子碰撞推进剂原子使之电离,进入放电室,工作介质氙气体在空心阴极发射的电子作用下在离子腔内离子化,被电离了的离子在含有电势差的离子光学系统(屏栅和加速栅)加速下以非常高的速度喷出发动机产生推力。离子被加速到所要求的排气速度后,由于物质是以离子形式喷出,携带净正电荷,中和器将发射等量的电子到离子束中确保电荷均衡。
为了确定离子推力器的关键失效模式及其失效机理,对其进行可靠性分析,国内外已经做了大量的地面寿命试验。美国NASA应用于深空一号的30cm NSTAR(SolarElectric Propulsion Technology Application Readiness)推力器地面寿命扩展试验累计时间达到30352h,而在NSTAR的基础上改进的应用于黎明号的离子推力器NEXT(NASA’sEvolutionary Xenon Thruster)创下了连续运行超50000小时的新世界纪录;德国的射频离子推力器RIT-10最终验证的寿命大于20000h;日本的10cm微波离子推力器地面试验寿命到2003年已经达18000h。而我国的LIPS-200离子推力器的寿命要求为10000~15000小时。离子推力器长寿命、高可靠性的特点给其可靠性分析带来一定的困难。
由于离子推力器的制造成本昂贵,同时受试验条件的制约,每次只能单机进行寿命试验,导致可靠性试验子样有限,属于极小子样范畴;同时受到研制周期紧迫的影响,其寿命试验往往无法长时间地进行,试验结果多为无失效数据。另外,由于对电推进技术的研究尚不够成熟,还不明确影响其寿命的加速因子,因此无法对其进行加速寿命试验。针对这种极小子样、长寿命和高可靠的特点,本发明提出一种基于性能退化的离子推力器小子样可靠性评估方法,来解决上述问题。
【发明内容】
本发明的目的就是针对上述离子推力器可靠性评估的特点和难点而提供一种方便有效的离子推力器小子样可靠性评估方法。
本发明一种基于性能退化的离子推力器小子样可靠性评估方法,可以通过以下步骤来实现:
步骤一:对离子推力器的结构特点和关键失效模式进行分析,选取影响其寿命的关键性能参数;
步骤二:根据关键性能参数y随时间的退化规律,建立离子推力器小子样性能退化模型;
步骤三:确定离子推力器平均失效时间tMTTF和可靠度函数R(t);
步骤四:对第i台离子推力器,记录在时间tij时的测试性能退化值yij,i=1,2,…,m,j=1,2,…,ni,m为投入试验的离子推力器台数,ni为第i台离子推力器性能测试的总次数;
步骤五:估计性能退化模型中的未知参数;
步骤六:利用可靠度函数对离子推力器进行小子样可靠性评估;
其中,在步骤一中所述的“选取影响其寿命的关键性能特征”,其选取的具体实现过程为:
加速栅结构失效是制约离子推力器寿命的最关键的失效模式之一,其工作原理如图2(a)和图2(b)所示。屏栅具有高电势,加速栅为低电势,两者结构相同并且产生电势差使放电室内的正离子高速喷出产生推力。加速栅的结构失效主要是由于加速栅的下游表面与中和区的电荷交换导致加速栅下游表面溅射腐蚀出凹槽和坑,如图3所示。交换电荷离子对加速栅相邻两孔之间的厚边溅射腐蚀出凹槽,相邻三孔之间腐蚀出坑,当凹槽和坑的腐蚀深度达到加速栅厚度时就会出现结构失效。
试验研究表明,中心区域凹槽最大腐蚀深度可以作为加速栅结构失效的判别依据,当最大腐蚀深度达到加速栅厚度时认为离子推力器寿命终止。因此,本发明选取加速栅中心区域凹槽最大腐蚀深度作为影响离子推力器寿命的关键性能参数y。
其中,在步骤二中所述的“建立离子推力器小子样性能退化模型”,其建立的具体实现过程为:
2.1.离子推力器性能退化路径的确定。试验表明,离子推力器中心区域凹槽最大腐蚀深度与时间的关系可以近似为线性关系,则对第i台离子推力器,有
yi(t)=ai+bit,i=1,2,…,m (1)
2.2.样本方差的确定。在t时刻,关键性能参数的方差的无偏估计为
其中,则,
因此,离子推力器性能退化过程的方差可用二次函数表示。
2.3.建立离子推力器小子样性能退化模型。考虑到同一台离子推力器不同测试时刻测得的性能参数之间的相关性,本发明提出一种独立增量线性过程性能退化模型,可以表示为
式中,a,b,d1,d2和d0均为未知参数,其中,d0≥0,d2≥0。根据独立增量过程定义可知,令yij=y(tij),i=1,2,…,m,j=1,2,...,ni,则对于测试时间yi1与性能退化增量之间相互独立,其中,
其中,在步骤三中所述的“确定离子推力器平均失效时间tMTTF和可靠度函数R(t)”,其确定的具体实现过程为:
3.1.假设离子推力器中心区域凹槽最大腐蚀深度达到Df时发生失效,则其平均失效时间公式为
3.2.可靠度函数R(t)可由以下公式给出
其中,在所述步骤五中,估计性能退化模型中的未知参数的具体实现过程如下:
5.1.重新参数化。令
式中,t1为各样本的第一个测量时间。则性能退化过程y(t)的方差项变为
5.2.定义似然函数。令Δyij=yij-yi(j-1),Δtij=tij-ti(j-1)则有
根据离子推力器性能退化数据,得到样本的对数似然函数为
5.3.确定边缘似然函数。固定参数θ,分别对未知参数a,b,求对数,同时考虑到参数估计的无偏性,可得
将式(17)—(20)带入式(16)得到关于θ的边缘似然函数为
5.4.未知参数估计。对于利用二分法对式(21)进行求解得到θ的估计值然后将带入式(17)~(19)可得到模型参数a,b,的估计值参数可由式(20)得到。
5.5.性能退化模型方差项未知参数估计。将参数估计值带入式(9)~(12)可得到模型方差项的未知参数估计值,进而进行可靠性评估。
本发明的优点与积极效果在于:
1.针对离子推力器寿命试验难以得到足够失效数据的情况建立了一种基于性能退化的小子样可靠性分析方法,解决了离子推力器小子样可靠性评估这一难题。
2.针对小样本退化试验的情况,本发明方法能够开发利用不同时刻性能测试数据间的纵向信息,使信息量大幅增加。所以,在试样数相同的情况下,本发明方法可以提高分析精度;而在精度相同的情况下,则可节省大量试样。
3.本发明给出了离子推力器给定退化失效阈值下的可靠度函数和平均失效时间公式,为以后离子推力器的可靠性评估提供了依据。
【附图说明】
图1为离子推力器结构示意图。
图2(a)为加速栅结构示意图。
图2(b)为加速栅工作原理示意图。
图3为加速栅凹槽和坑腐蚀示意图。
图4为本发明方法流程示意图。
【符号说明】
本发明中符号说明如下:
t为性能退化测试时刻;
yi(t)为第i台离子推力器在t时刻测得的性能退化量;
ai,bi为第i台离子推力器的退化模型的待求参数;
S2(t)为性能退化量在t时刻的方差无偏估计;
为性能退化量在t时刻的均值;
为退化模型参数的均值;
ε(t)为退化模型在t时刻的方差项;
d1,d2,d0为退化模型方差项中的待求参数;
N()为正态分布;
Cov()为协方差;
yij为第i台离子推力器在第j次测得的性能退化量;
tij为第i台离子推力器在第j次时的测试时刻;
μ1为离子推力器第1次测得的性能退化量的均值;
为离子推力器第1次测得的性能退化量的方差;
Df为离子推力器性能退化量的阈值;
Φ()为标准正态分布函数;
I()为退化模型重新参数化后的参数公式;
θ,为退化模型重新参数化后方差项中的待求参数;
Δyij为第i台离子推力器在第j次和第j-1次测得的性能退化量增量;
Δtij为第i台离子推力器在第j次和第j-1次测试时刻增量;
为第i台离子推力器在第j次和第j-1次测试时刻的均值;
【具体实施方式】
下面将结合附图对本发明作进一步的详细说明。
本发明提供一种基于性能退化的离子推力器极小子样可靠性分析方法,其流程图如图4所示,本发明方法的具体实现包括以下步骤:
步骤一:对离子推力器的结构特点及主要失效模式进行分析,选取影响其寿命的关键性能参数,具体实施方式如下:
1.离子推力器结构复杂,失效模式多,影响离子推力器运行寿命的关键部组件是栅极系统和阴极组件,其主要失效模式有电子返流、加速栅结构失效、发射体耗尽、加热丝熔断和触持极磨损等。研究显示,加速栅结构失效是制约离子推力器寿命的最关键的失效模式之一,其工作原理如图2所示。屏栅具有高电势,加速栅为低电势,两者结构相同并且产生电势差使放电室内的正离子高速喷出产生推力。加速栅的结构失效主要是由于加速栅的下游表面与中和区的电荷交换导致加速栅下游表面溅射腐蚀出凹槽和坑,如图3所示。交换电荷离子对加速栅相邻两孔之间的厚边溅射腐蚀出凹槽,相邻三孔之间腐蚀出坑,当凹槽和坑的腐蚀深度达到加速栅厚度时就会出现结构失效。
2.试验研究表明,中心区域凹槽最大腐蚀深度可以作为加速栅结构失效的判别依据,当最大腐蚀深度达到加速栅厚度时认为离子推力器寿命终止。因此,本发明选取加速栅中心区域凹槽最大腐蚀深度作为影响离子推力器寿命的关键性能参数y。
步骤二:采集离子推力器关键性能参数退化数据。假设有m台离子推力器进行地面寿命试验,对第i台离子推力器,i=1,2,…,m,分别在时刻处测量加速栅中心区域凹槽最大腐蚀深度并记录相应的性能退化值,ni为第i台离子推力器性能测试的时刻数。
步骤三:建立离子推力器小子样性能退化模型。本发明采用一种独立增量线性过程模型进行离子推力器小子样可靠性评估,该模型可以表示为
步骤四:模型参数估计。具体实施过程如下:
1.重新参数化。首先对模型参数做如下变换:
式中,t1为各样本的第一个测量时间。则性能退化过程y(t)的方差项变为
2.定义似然函数。令Δyij=yij-yi(j-1),△tij=tij-ti(j-1)根据离子推力器性能退化数据,得到样本的对数似然函数为
3.确定边缘似然函数。固定参数θ,分别对未知参数a,b,求对数,同时考虑到参数估计的无偏性,可得
将式(28)~(31)带入式(32)得到关于θ的边缘似然函数为
4.未知参数估计。对于利用二分法对式(32)进行求解得到θ的估计值将估计值分别带入式(28)~(30)可得到模型参数a,b,的估计值参数可由式(31)得到。
5.性能退化模型方差项未知参数估计。将参数估计值分别带入式(22)~(25)可得到模型方差项的参数估计值,进而进行可靠性评估。
步骤五:离子推力器极小子样可靠性评估。具体实施过程如下:
1.可靠度估计。若b>0,离子推力器在给定寿命t的可靠度为
若b<0,离子推力器在给定寿命t的可靠度为
式中,Df为离子推力器中心区域凹槽最大腐蚀深度失效阈值。
2.平均失效时间估计。离子推力器平均失效时间为

Claims (3)

1.一种基于性能退化的离子推力器小子样可靠性评估方法,其特征在于,包括以下步骤:
步骤一:对离子推力器的结构特点和关键失效模式进行分析,选取影响其寿命的关键性能参数;
步骤二:根据关键性能参数y随时间的退化规律,建立离子推力器小子样性能退化模型;
步骤三:确定离子推力器平均失效时间tMTTF和可靠度函数R(t);
步骤四:对第i台离子推力器,记录在时间tij时的测试性能退化值yij,i=1,2,...,m,j=1,2,...,ni,m为投入试验的离子推力器台数,ni为第i台离子推力器性能测试的总次数;
步骤五:估计性能退化模型中的未知参数;
步骤六:利用可靠度函数对离子推力器进行小子样可靠性评估;
在步骤一中所述的选取影响其寿命的关键性能参数,为选取加速栅中心区域凹槽最大腐蚀深度作为影响离子推力器寿命的关键性能参数y;
在步骤二中所述的建立离子推力器小子样性能退化模型,包括以下步骤:
步骤2.1.离子推力器性能退化路径的确定;
离子推力器中心区域凹槽最大腐蚀深度与时间的关系为线性关系,则对第i台离子推力器,有
yi(t)=ai+bit,i=1,2,...,m (1)
步骤2.2.样本方差的确定;
在t时刻,关键性能参数的方差的无偏估计为
其中,则,
因此,离子推力器性能退化过程的方差用二次函数表示;
步骤2.3.建立离子推力器小子样性能退化模型;
考虑到同一台离子推力器不同测试时刻测得的性能参数之间的相关性,因此采用独立增量线性过程性能退化模型,表示为
式中,a,b,d1,d2和d0均为未知参数,其中,d0≥0,d2≥0;令yij=y(tij),i=1,2,...,m,j=1,2,...,ni,则对于测试时间yi1与性能退化增量之间相互独立,ε(t)为退化模型在t时刻的方差项;
其中,
2.根据权利要求1所述的一种基于性能退化的离子推力器小子样可靠性评估方法,其特征在于:在步骤三中所述的确定离子推力器平均失效时间tMTTF和可靠度函数R(t),包括以下步骤:
步骤3.1.假设离子推力器中心区域凹槽最大腐蚀深度达到Df时发生失效,则其平均失效时间公式为
步骤3.2.可靠度函数R(t)由以下公式给出
3.根据权利要求1所述的一种基于性能退化的离子推力器小子样可靠性评估方法,其特征在于:在所述步骤五中,估计性能退化模型中的未知参数的具体实现过程如下:
步骤5.1.重新参数化;令
式中,t1为各样本的第一个测量时间;则性能退化过程y(t)的方差项变为
步骤5.2.定义似然函数;
令Δyij=yij-yi(j-1),Δtij=tij-ti(j-1)则有
根据离子推力器性能退化数据,得到样本的对数似然函数为
步骤5.3.确定边缘似然函数;
固定参数θ,分别对未知参数a,b,求对数,同时考虑到参数估计的无偏性,得到
将式(17)—(20)带入式(16)得到关于θ的边缘似然函数为
步骤5.4.未知参数估计;
对于利用二分法对式(21)进行求解得到θ的估计值然后将带入式(17)~(19)得到模型参数a,b,的估计值参数由式(20)得到;
步骤5.5.性能退化模型方差项未知参数估计;
将参数估计值带入式(9)~(12)得到模型方差项的未知参数估计值,进而进行可靠性评估。
CN201610165279.8A 2016-03-22 2016-03-22 一种基于性能退化的离子推力器小子样可靠性评估方法 Active CN105844095B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610165279.8A CN105844095B (zh) 2016-03-22 2016-03-22 一种基于性能退化的离子推力器小子样可靠性评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610165279.8A CN105844095B (zh) 2016-03-22 2016-03-22 一种基于性能退化的离子推力器小子样可靠性评估方法

Publications (2)

Publication Number Publication Date
CN105844095A CN105844095A (zh) 2016-08-10
CN105844095B true CN105844095B (zh) 2018-08-17

Family

ID=56588078

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610165279.8A Active CN105844095B (zh) 2016-03-22 2016-03-22 一种基于性能退化的离子推力器小子样可靠性评估方法

Country Status (1)

Country Link
CN (1) CN105844095B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107966311B (zh) * 2017-11-24 2019-08-09 中国空间技术研究院 基于加速栅数据的极小子样离子推力器可靠度确定方法
CN108280253B (zh) * 2017-12-20 2020-07-14 中国空间技术研究院 基于栅极腐蚀形貌和电子返流的离子推力器寿命评估方法
CN111199099B (zh) * 2019-12-26 2023-01-31 兰州空间技术物理研究所 一种基于栅极腐蚀的离子推力器运行寿命评估方法
CN114580118B (zh) * 2022-03-25 2023-12-15 兰州空间技术物理研究所 一种离子推力器寿命及可靠性定量评估方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
基于性能参数退化的航天继电器可靠性建模与评估研究;孟彦辰;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20160215;第2016年卷(第02期);C031-86页 *
基于性能退化数据的可靠性评估;邓爱民等;《宇航学报》;20060531;第27卷(第3期);546-552页 *
极小子样加速寿命试验设计和可靠性分析方法;傅惠民;《机械强度》;20131231;第35卷(第1期);22-27页 *
考虑电荷交换的栅极区离子流数值模拟;刘洋等;《真空与低温》;20111231;第17卷(第4期);198-208页 *

Also Published As

Publication number Publication date
CN105844095A (zh) 2016-08-10

Similar Documents

Publication Publication Date Title
CN105844095B (zh) 一种基于性能退化的离子推力器小子样可靠性评估方法
CN103995970B (zh) 一种离子推力器极小子样可靠性评估方法
CN106772080B (zh) 空间锂离子电池加速退化试验时间等效性建模方法
CN111199099B (zh) 一种基于栅极腐蚀的离子推力器运行寿命评估方法
Thiel et al. An improved model for gaseous amplification in the environmental SEM
CN112347610B (zh) 一种基于空心阴极的霍尔推力器的寿命评估方法
CN112329247A (zh) 一种霍尔电推力器放电等离子体和自溅射的仿真方法
CN114388855B (zh) 一种燃料电池寿命预测方法
CN108280253B (zh) 基于栅极腐蚀形貌和电子返流的离子推力器寿命评估方法
Huismann Improving hall thruster plume simulation through refined characterization of near-field plasma properties
Upadhyay et al. Experimentally validated computations of simultaneous ion and fast neutral energy and angular distributions in a capacitively coupled plasma reactor
CN110851939A (zh) 圆柱形阳极层霍尔推力器寿命评估的仿真计算方法
Dolenc et al. Hybrid approach to remaining useful life prediction of solid oxide fuel cell stack
CN113409363B (zh) 一种基于bp-pmbm滤波算法的多目标跟踪方法
CN107995766A (zh) 一种基于空心阴极磁特性的出口电子参数获得方法
Boyd Computation of the plume of an anode-layer Hall thruster
CN115495930B (zh) 一种变孔径栅极离子推力器性能计算方法
Oblapenko et al. Hedging direct simulation Monte Carlo bets via event splitting
CN114580118B (zh) 一种离子推力器寿命及可靠性定量评估方法
US6199029B1 (en) Topography simulation method and system of plasma-assisted etching process
Walton et al. Time-resolved ion flux measurements in pulsed, electron-beam-generated plasmas
Tartz et al. Validated simulation of the ion extraction grid lifetime
WAng et al. Lifetime prediction method for electron multiplier based on accelerated degradation test
Boyd et al. Grid erosion analysis of the T5 ion thruster
CN111460751A (zh) 一种SiC MOSFET器件损伤机理的分析方法、装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210207

Address after: 315100 b-211-47, Kexin building, 655 bachelor Road, Yinzhou District, Ningbo City, Zhejiang Province

Patentee after: NINGBO ZHIZHENG WEIYING INFORMATION TECHNOLOGY Co.,Ltd.

Address before: 100191 No. 37, Haidian District, Beijing, Xueyuan Road

Patentee before: BEIHANG University