一种卫星平台可靠性增长评估方法和系统与流程

1.本发明属于卫星总体设计技术领域，尤其涉及一种卫星平台可靠性增长评估方法和系统。


背景技术：

2.对卫星系统的可靠性评估一直都存在评估对象样本不足的问题。除了大型星座系统，一般很少有多颗卫星是完全一样的，即使是同一型号，装备星也会随着质量问题归零和技术升级在首发星的基础上做出适当的改进。样本量小和技术状态持续变化，是卫星系统的普遍状态。
3.虽然大部分卫星载荷差别很大，但对作为卫星系统中技术状态较为固定的卫星平台来说，基于同一平台的卫星，其平台部分技术状态基本都较为相近。
4.经典的可靠性评估方法和贝叶斯方法都需要假定被评估样本技术状态稳定不变，故，对于技术状态处于不断改进中的样本来说，经典的可靠性评估方法和贝叶斯方法都不能适用。


技术实现要素：

5.本发明的技术解决问题：克服现有技术的不足，提供一种卫星平台可靠性增长评估方法和系统，将卫星平台的不断改进看作是一个无计划的可靠性增长过程，使用可靠性增长模型来评估卫星平台的可靠性，给出卫星平台当前的可靠性水平，解决了传统可靠性评估方法无法应对样本技术状态不断变化的情况。
6.为了解决上述技术问题，本发明公开了一种卫星平台可靠性增长评估方法，包括：
7.收集卫星平台上所有已发射卫星的在轨可靠性信息，确定各已发射卫星的累积运行时间和累积关联故障次数；
8.按照各已发射卫星的发射时间的先后顺序，对卫星平台上所有已发射卫星进行排序并重新编号：1,2,...,i,...,k；其中，k表示卫星平台上所有已发射卫星的总数量；
9.根据确定的各已发射卫星的累积运行时间和累积关联故障次数，计算得到卫星平台截止到第i颗已发射卫星的平台累积运行时间ti、平台累积关联故障次数ni和平台累积平均故障间隔时间mtbf
∑i
；
10.对ti和mtbf
∑i
取自然对数后进行线性拟合，得到相关系数ρ，并确定第一经验公式和第二经验公式；
11.根据第一经验公式和第二经验公式，解算得到卫星平台截至到当前时刻的平台瞬时平均故障间隔时间；
12.根据解算得到的卫星平台截至到当前时刻的平台瞬时平均故障间隔时间，采用指数分布假设，评估得到所述待分析卫星平台在典型寿命期的可靠性水平。
13.在上述卫星平台可靠性增长评估方法中，卫星平台截止到第i颗已发射卫星的平台累积运行时间ti为：
[0014][0015]
其中，tj表示第j颗已发射卫星的累积运行时间，j＝1,2,...,i。
[0016]
在上述卫星平台可靠性增长评估方法中，卫星平台截止到第i颗已发射卫星的平台累积关联故障次数ni为：
[0017][0018]
其中，nj表示第j颗已发射卫星的累积关联故障次数。
[0019]
在上述卫星平台可靠性增长评估方法中，卫星平台截止到第i颗已发射卫星的平台累积平均故障间隔时间mtbf
∑i
为：
[0020][0021]
在上述卫星平台可靠性增长评估方法中，相关系数ρ的解算公式如下：
[0022][0023]
在上述卫星平台可靠性增长评估方法中，
[0024]
第一经验公式表示如下：
[0025][0026]
其中，表示卫星平台截至到当前时刻的平台累积平均故障间隔时间，a和m为待求解参数；
[0027]
第二经验公式表示如下：
[0028][0029]
其中，表示卫星平台截至到当前时刻的平台瞬时平均故障间隔时间。
[0030]
在上述卫星平台可靠性增长评估方法中，确定第一经验公式和第二经验公式，包括：
[0031]
对ti和mtbf
∑i
的自然对数ln(ti)和ln(mtbf
∑i
)进行多点最小二乘法拟合，得到待求解参数a和m的值；
[0032]
将得到的待求解参数a和m的值分别代入公式(2)和(3)，即可得到第一经验公式和第二经验公式。
[0033]
在上述卫星平台可靠性增长评估方法中，待求解参数a和m的解算公式如下：
[0034][0035][0036]
在上述卫星平台可靠性增长评估方法中，在轨可靠性信息，包括：发射时间、在轨关联故障的发生次数和每次在轨关联故障发生时间；其中，经证实是未按规定的条件使用而引起的故障，或已经证实仅属某项将不采用的设计所引起的故障统称为在轨非关联故障，否则就是在轨关联故障。
[0037]
相应的，本发明还公开了一种卫星平台可靠性增长评估系统，包括：
[0038]
数据获取模块，用于收集卫星平台上所有已发射卫星的在轨可靠性信息，确定各已发射卫星的累积运行时间和累积关联故障次数；
[0039]
排序模块，用于按照各已发射卫星的发射时间的先后顺序，对卫星平台上所有已发射卫星进行排序并重新编号：1,2,...,i,...,k；其中，k表示卫星平台上所有已发射卫星的总数量；
[0040]
第一计算模块，用于根据确定的各已发射卫星的累积运行时间和累积关联故障次数，计算得到卫星平台截止到第i颗已发射卫星的平台累积运行时间ti、平台累积关联故障次数ni和平台累积平均故障间隔时间mtbf
∑i
；
[0041]
第二计算模块，用于对ti和mtbf
∑i
取自然对数后进行线性拟合，得到相关系数ρ，并确定第一经验公式和第二经验公式；
[0042]
第三计算模块，根据第一经验公式和第二经验公式，解算得到卫星平台截至到当前时刻的平台瞬时平均故障间隔时间；
[0043]
评估模块，根据解算得到的卫星平台截至到当前时刻的平台瞬时平均故障间隔时间，采用指数分布假设，评估得到所述待分析卫星平台在典型寿命期的可靠性水平。
[0044]
本发明具有以下优点：
[0045]
本发明公开了一种卫星平台可靠性增长评估方案，将卫星平台的不断改进看作是一个无计划的可靠性增长过程，使用可靠性增长模型来评估卫星平台的可靠性，给出卫星平台当前的可靠性水平，解决了传统可靠性评估方法无法应对样本技术状态不断变化的情况。其次，通过本发明评估出的可靠性增长率可用来评估卫星平台的可靠性改进提升情况。
附图说明
[0046]
图1是本发明实施例中一种卫星平台可靠性增长评估方法的步骤流程图；
[0047]
图2是本发明实施例中一种某卫星平台ln(ti)与ln(mtbf
∑i
)的线性拟合图。
具体实施方式
[0048]
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明公开的实施方式作进一步详细描述。
[0049]
如图1，在本实施例中，该卫星平台可靠性增长评估方法，包括：
[0050]
步骤101，收集卫星平台上所有已发射卫星的在轨可靠性信息，确定各已发射卫星的累积运行时间和累积关联故障次数。
[0051]
在本实施例中，轨可靠性信息包括但不仅限于：发射时间、在轨关联故障的发生次数和每次在轨关联故障发生时间；其中，经证实是未按规定的条件使用而引起的故障，或已经证实仅属某项将不采用的设计所引起的故障统称为在轨非关联故障，否则就是在轨关联故障。
[0052]
步骤102，按照各已发射卫星的发射时间的先后顺序，对卫星平台上所有已发射卫星进行排序并重新编号：1,2,...,i,...,k。
[0053]
在本实施例中，k表示卫星平台上所有已发射卫星的总数量。例如，排序并重新编号的已发射卫星表示如下：
[0054]
序号卫星关联故障数已运行时间1卫星169年8个月2卫星207年7个月3卫星307年1个月4卫星405年8个月5卫星505年8个月
…………
[0055]
表1，已发射卫星排序表
[0056]
步骤103，根据确定的各已发射卫星的累积运行时间和累积关联故障次数，计算得到卫星平台截止到第i颗已发射卫星的平台累积运行时间ti、平台累积关联故障次数ni和平台累积平均故障间隔时间mtbf
∑i
。
[0057]
在本实施例中，按照各已发射卫星的发射时间的先后顺序，卫星平台截止到第i颗已发射卫星的平台累积运行时间ti、平台累积关联故障次数ni和平台累积平均故障间隔时间mtbf
∑i
的解算公式分别如下：
[0058][0059][0060][0061]
其中，tj表示第j颗已发射卫星的累积运行时间，nj表示第j颗已发射卫星的累积关联故障次数，j＝1,2,...,i。
[0062]
优选的，基于表1可知，tj和nj与ti、ni和mtbf
∑i
的关系如表2所示：
[0063]
卫星tj(h)nj(次)ti(h)ni(次)mtbf
∑i
卫星185032685032614172卫星215146462364961219708卫星321384064503361825018.6
卫星426388067142162429759卫星5313920610281363034271.2
………ꢀꢀ…
[0064]
表2，tj和nj与ti、ni和mtbf
∑i
的关系示意表
[0065]
步骤104，对ti和mtbf
∑i
取自然对数后进行线性拟合，得到相关系数ρ，并确定第一经验公式和第二经验公式。
[0066]
在本实施例中，相关系数ρ的解算公式如下：
[0067][0068]
第一经验公式表示如下：
[0069][0070]
其中，表示卫星平台截至到当前时刻的平台累积平均故障间隔时间；a和m为待求解参数，用于表征卫星平台的初始mtbf水平，m表示卫星平台的mtbf增长率。
[0071]
第二经验公式表示如下：
[0072][0073]
其中，表示卫星平台截至到当前时刻的平台瞬时平均故障间隔时间。
[0074]
对ti和mtbf
∑i
的自然对数ln(ti)和ln(mtbf
∑i
)进行多点最小二乘法拟合，得到待求解参数a和m的值：
[0075][0076][0077]
将通过公式(4)和(5)解算得到的待求解参数a和m的值分别代入公式(2)和(3)，即可得到第一经验公式和第二经验公式。
[0078]
步骤105，根据第一经验公式和第二经验公式，解算得到卫星平台截至到当前时刻的平台瞬时平均故障间隔时间。
[0079]
步骤106，根据解算得到的卫星平台截至到当前时刻的平台瞬时平均故障间隔时间，采用指数分布假设，评估得到所述待分析卫星平台在典型寿命期的可靠性水平。
[0080]
在上述实施例的基础上，下面以一个具体实例进行说明。
[0081]
(1)收集卫星平台上所有已发射卫星的在轨可靠性信息。
[0082]
某卫星平台在2004年～2014年的10年间累积发射入轨卫星20颗，截至2014年5月23日，该平台卫星累积运行75.42星*年。该卫星平台经过不断的改进和完善，目前技术成熟度和在轨可靠性都非常高。
[0083]
(2)按照各已发射卫星的发射时间的先后顺序，对卫星平台上所有已发射卫星进行排序并重新编号：1,2,...,i,...,k。
[0084]
按照各已发射卫星的发射时间的先后顺序，对2004年～2014年(截至2014-05-23)该卫星平台发射入轨的20颗卫星进行排序，得到下表3：
[0085]
卫星tj(h)nj(次)卫星1850326卫星2664320卫星3623760卫星4500400卫星5500400卫星6488640卫星7396721卫星8369600卫星9369600卫星10369601卫星11318000卫星12252482卫星13246961卫星14130560卫星15130560卫星16130560卫星1774880卫星1863120卫星1963120卫星2063120
[0086]
表3
[0087]
(3)计算得到卫星平台截止到第i颗已发射卫星的平台累积运行时间ti、平台累积关联故障次数ni和平台累积平均故障间隔时间mtbf
∑i
，如表4所示：
[0088][0089][0090]
表4
[0091]
(4)对ti和mtbf
∑i
取自然对数后进行线性拟合，得到相关系数ρ，并确定第一经验公式和第二经验公式。
[0092]
对ti和mtbf
∑i
取自然对数后进行线性拟合，如图2所示，得到相关系数ρ＝0.0940。
[0093]
对表4中的数据ln(ti)和ln(mtbf
∑i
)进行多点最小二乘法拟合，得到m＝0.6126，a＝0.0548。
[0094]
第一经验公式表示如下：
[0095][0096]
第二经验公式表示如下：
[0097][0098]
将m＝0.6126，a＝0.0548代入式(2)和式(3)可得：
[0099]
[0100][0101]
(5)根据第一经验公式和第二经验公式，解算得到卫星平台截至到当前时刻的平台瞬时平均故障间隔时间。
[0102]
通过计算，卫星平台截至到2014-5-23的平台瞬时平均故障间隔时间为173126.2h，记作：mtbf
now
＝173126.2h。
[0103]
(6)根据解算得到的卫星平台截至到当前时刻的平台瞬时平均故障间隔时间，采用指数分布假设，评估得到所述待分析卫星平台在典型寿命期的可靠性水平。
[0104]
对该卫星平台的可靠性增长水平进行分析，如下表5所示：
[0105][0106]
表5
[0107]
在开展有计划的可靠性增长试验时，一般取m＝0.3～0.7。当m＝0.l时，说明在增长过程中基本没有采取改进措施；当m＝0.6～0.7时，说明在增长过程中采取了强有力的故障分析和改进措施，取得了较大的增长效果。虽然本方法只是将卫星平台的不断改进看作是一个无计划的可靠性增长过程，但具体实施例的分析结果说明了该卫星平台的可靠性增长显著。
[0108]
在上述实施例的基础上，本发明还公开了一种卫星平台可靠性增长评估系统，包括：数据获取模块，用于收集卫星平台上所有已发射卫星的在轨可靠性信息，确定各已发射卫星的累积运行时间和累积关联故障次数；排序模块，用于按照各已发射卫星的发射时间的先后顺序，对卫星平台上所有已发射卫星进行排序并重新编号：1,2,...,i,...,k；第一计算模块，用于根据确定的各已发射卫星的累积运行时间和累积关联故障次数，计算得到卫星平台截止到第i颗已发射卫星的平台累积运行时间ti、平台累积关联故障次数ni和平台累积平均故障间隔时间mtbf
∑i
；第二计算模块，用于对ti和mtbf
∑i
取自然对数后进行线性拟合，得到相关系数ρ，并确定第一经验公式和第二经验公式；第三计算模块，根据第一经验公式和第二经验公式，解算得到卫星平台截至到当前时刻的平台瞬时平均故障间隔时间；评估模块，根据解算得到的卫星平台截至到当前时刻的平台瞬时平均故障间隔时间，采用指数分布假设，评估得到所述待分析卫星平台在典型寿命期的可靠性水平。
[0109]
对于系统实施例而言，由于其与方法实施例相对应，所以描述的比较简单，相关之处参见方法实施例部分的说明即可。
[0110]
本发明虽然已以较佳实施例公开如上，但其并不是用来限定本发明，任何本领域技术人员在不脱离本发明的精神和范围内，都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改，因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰，均属于本发明技术方案的保护范围。
[0111]
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。