CN106682273B - 一种串并联混合航天装备系统寿命重要度确定方法 - Google Patents

一种串并联混合航天装备系统寿命重要度确定方法 Download PDF

Info

Publication number
CN106682273B
CN106682273B CN201611103084.7A CN201611103084A CN106682273B CN 106682273 B CN106682273 B CN 106682273B CN 201611103084 A CN201611103084 A CN 201611103084A CN 106682273 B CN106682273 B CN 106682273B
Authority
CN
China
Prior art keywords
life
subsystem
service life
series
members
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.)
Expired - Fee Related
Application number
CN201611103084.7A
Other languages
English (en)
Other versions
CN106682273A (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.)
CHINA AEROSPACE STANDARDIZATION INSTITUTE
Original Assignee
CHINA AEROSPACE STANDARDIZATION INSTITUTE
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 CHINA AEROSPACE STANDARDIZATION INSTITUTE filed Critical CHINA AEROSPACE STANDARDIZATION INSTITUTE
Priority to CN201611103084.7A priority Critical patent/CN106682273B/zh
Publication of CN106682273A publication Critical patent/CN106682273A/zh
Application granted granted Critical
Publication of CN106682273B publication Critical patent/CN106682273B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Abstract

本发明公开了一种串并联混合系统寿命重要度确定方法,可以建立串并联模型下设备寿命和系统寿命之间的关系,定量化设备寿命的变化对系统寿命影响,度量不同位置不同寿命水平的设备对实现系统寿命指标的重要度。为串并联混合系统中不同设备重要度提供分析方法,可以用于卫星、导弹、空间站等大型航天装备系统开展不同寿命设计方案的比较和优化,并可以用于确定寿命薄弱环节,从而可以支持在有限的时间和经费条件下实现寿命提升效益的最大化。

Description

一种串并联混合航天装备系统寿命重要度确定方法
技术领域
本发明涉及可靠性技术领域,具体涉及一种串并联混合系统寿命重要度确定方法。
背景技术
随着航天技术的发展,对大型复杂系统寿命的要求越来越高,例如一般要求高轨卫星寿命15年,低轨卫星寿命8年,导弹贮存寿命20年。如何有效提高复杂系统的寿命是航天可靠性技术领域的一项难题。大型复杂系统寿命是以设备的寿命为基础的,当系统寿命不满足要求时,需要改进设备的寿命来实现系统寿命的提升。但是,航天型号构成复杂,构成设备达数百种,而且不同设备之间的冗余逻辑关系也非常复杂,很难直接判断不同设备影响系统寿命的程度,造成系统寿命工作难以开展。
重要度可以表征反映系统中各单元的重要程度。目前只有系统寿命建模方法和可靠性重要度计算的相关方法,没有针对寿命指标的重要度计算方法,无法支持复杂系统寿命设计及寿命提升工作。另外,当前的系统寿命模型中,都是假定不同设备寿命为独立同分布的变量,复杂系统中的设备多种多样,这个假定不符合工程实际情况。
发明内容
航天装备中最常见的设备之间的逻辑关系就是串联和并联,有鉴于此,本发明提供了一种串并联混合系统寿命重要度确定方法。
一种串并联混合航天装备系统寿命重要度确定方法,包括如下步骤:
步骤1、根据系统中设备之间的串并联关系,将系统从上到下进行分解,得到具有多个层级的逻辑关系树;所述逻辑关系树的每一层级中,均为独立的串联或并联关系的子系统,子系统中包括至少2个互为串联或并联关系的成员;逻辑关系树中除最高层级的任意层级X的子系统为上一层级X-1中一个子系统的成员;
步骤2、针对逻辑关系树的各个层级中子系统,基于子系统与成员之间的串联或并联关系,根据子系统中各个成员的寿命,确定子系统的寿命,由此得到子系统寿命与成员寿命之间的模型;
步骤3、首先,从所述逻辑关系树的最底层开始:
A、根据步骤2获得的子系统寿命与成员寿命之间的模型,结合成员的寿命分布函数,基于子系统与成员之间的串联或并联关系,获得用成员寿命分布函数表征的各个子系统的寿命分布函数;
B、再根据子系统的寿命分布函数获得各个子系统的寿命均值;
然后,遍历最底层的各个子系统,得到子系统的寿命均值后,将子系统的寿命均值作为上一层级中所在子系统成员的寿命,采用步骤A和步骤B的方法,获得该层级子系统的寿命均值;依次类推,直到得到最高层级的寿命均值,即系统的寿命均值;
步骤4、用系统寿命均值对任意成员寿命求偏导数,即为该成员的寿命重要度。
较佳的,所述步骤2中,当子系统的成员之间为并联关系时,子系统寿命为成员中寿命最大值;当子系统的成员之间为串联关系时,子系统寿命为成员中寿命最小值。
较佳的,所述步骤3中,组成串联或并联关系子系统的单个成员的寿命分布函数服从指数分布。
较佳的,所述步骤4中,计算设备级成员的寿命重要度时,先求系统寿命对子系统偏导数,再乘上子系统寿命对成员寿命的偏导数。
本发明具有如下有益效果:
本发明可以建立串并联模型下设备寿命和系统寿命之间的关系,定量化设备寿命的变化对系统寿命影响,度量不同位置不同寿命水平的设备对实现系统寿命指标的重要度。为串并联混合系统中不同设备重要度提供分析方法,可以用于卫星、导弹、空间站等大型航天装备系统开展不同寿命设计方案的比较和优化,并可以用于确定寿命薄弱环节,从而可以支持在有限的时间和经费条件下实现寿命提升效益的最大化。
附图说明
图1为本发明的方法流程图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了一种串并联混合系统寿命重要度排序方法,如图1所示,包括如下步骤:
步骤1、将串并联混合系统分解为单一的串联子系统或者并联子系统。根据系统中设备之间的串并联关系,将系统从上到下进行分解,得到具有多个层级的逻辑关系树;所述逻辑关系树的每一层级中,均为独立的串联或并联关系的子系统,子系统中包括至少2个互为串联或并联关系的成员;逻辑关系树中除最高层级的任意层级X的子系统为上一层级X-1中一个子系统的成员。例如,首先将系统进行分解,假设可分解为S个子系统,记为A1,...,As。A1,...,As如果不是单一串联子系统或并联子系统,还可以再往下分解,例如将Ar,r=1,...,s分解为
Figure GDA0001228169910000031
hr为第r个子系统可以再往下分解的子系统的数量。直到所有的子系统都是单一串联子系统或并联子系统或者单个设备。
步骤2、建立系统寿命与设备寿命之间的模型。系统寿命为L,系统由n个设备构成,设备的寿命为L1,L2,...,Ln
当子系统的成员之间为串联关系时,子系统寿命为成员中寿命最小值,即:
L=min(L1,L2,...,Ln) (1)
当子系统的成员之间为并联关系时,子系统寿命为成员中寿命最大值,即:
L=max(L1,L2,...,Ln) (2)
统一记为L=f(L1,L2,...,Ln),对于串联系统,f为min;对于并联系统,f为max。
对串并联混合系统,子系统A1,...,As的寿命对应为L(1),L(2),...,L(s)。单一串联或并联子系统寿命模型可由构成设备寿命按公式(1)或公式(2)获得,则
L=f(L(1),L(2),...,L(s)) (3)
步骤3、从所述逻辑关系树的最底层开始向上计算,得到系统寿命均值,具体为:由于复杂系统中设备寿命一般服从指数分布,则可以假定设备寿命分布为:
Figure GDA0001228169910000043
即第i个设备的寿命均值为θi。下面首先分别给出串联系统和并联系统的计算方法,再推广到串并联混合系统,给出串并联混合系统的一般性的计算方法。
3.1并联系统
根据子系统寿命与成员寿命之间的模型,结合成员的寿命分布函数,基于子系统与成员之间的并联关系,获得用成员寿命分布函数表征的各个子系统的寿命分布函数;
Figure GDA0001228169910000041
m1,...,mn都可以取1,2,...,n中的任一个值,这组下标用于表示k(k=1,2,...,n)个不同设备。
然后得到子系统寿命均值:
Figure GDA0001228169910000042
对于相同产品并联,有以下更简单的计算方法:
n个产品独立同分布,分布函数为
Figure GDA0001228169910000051
Figure GDA0001228169910000052
Figure GDA0001228169910000053
Figure GDA0001228169910000054
所以
Figure GDA0001228169910000055
3.2串联系统
根据子系统寿命与成员寿命之间的模型,结合成员的寿命分布函数,基于子系统与成员之间的串联关系,获得用成员寿命分布函数表征的各个子系统的寿命分布函数;
Figure GDA0001228169910000056
则子系统寿命均值为:
Figure GDA0001228169910000057
3.3串并联混合系统
遍历最底层的各个子系统,得到子系统的寿命均值后,将子系统的寿命均值作为上一层次中所在子系统成员的寿命,获得该层次子系统的寿命均值;依次类推,直到得到最高层次的寿命均值,即系统的寿命均值,具体过程为:
单一串联子系统或并联子系统寿命均值可按公式(5)或公式(10)计算,统一记为θ=g(θ1,...,θn)。子系统A1,...,As寿命均值记为θ(1)(2),...,θ(s),子系统寿命均值可由再下一层的子系统或者设备寿命均值按照公式(5)或(10)计算获得。则串并联混合系统寿命均值为:
θ=g(θ(1)(2),...,θ(s)) (11)
步骤4、计算寿命重要度。
寿命重要度为设备寿命变化对系统寿命的影响,即系统寿命对第i个设备寿命的偏导数:
Figure GDA0001228169910000061
对于复杂系统,直接按公式(12)计算难度比较大,可以先求系统寿命对子系统偏导数,再乘上子系统寿命对设备寿命的偏导数。若第i个设备属于第r个子系统,则
Figure GDA0001228169910000062
根据实际需求,根据公式(12)或(13)计算的结果对不同设备寿命重要度进行排序,得到
I(1)≤...≤I(n) (14)
I(i)是对设备重要度按从小到大排序后的第i个重要度。寿命重要度最高的设备可作为系统寿命改进的关键设备。
下面举实施例来进一步说明本发明所提供的方法。
某系统由四个设备构成,设备之间的逻辑关系为设备1与设备2并联,然后与设备3串联,再与设备4并联。四个设备的寿命分别为L1,L2,L3,L4,均服从指数分布,寿命均值θ1为4年,θ2为6年,θ3为5年,θ4为7年。
步骤1、对系统进行分解
首先将系统进行分解,可以分解为子系统A1,A2,A1由设备1,2,3构成,A2为设备4。A1不是单一串联系统或并联系统或单个设备,还需往下分解。A2为单个设备,不需要再往下分解。A1可以分解为A11,A12,A11由设备1,2构成,A12为设备3。A11为单一并联系统,A12为单个设备,都不需要再往下分解。
步骤2、建立系统寿命与设备寿命之间的模型。
子系统A11,A12的寿命记为L(1,1),L(1,2)。由公式(2),L(1,1)=min(L1,L2),L(1,2)=L3。同样,由公式(1)可得到A1的寿命L(1)=max(L(1,1),L(1,2))。A2的寿命L(2)=L4。系统寿命模型为:
L=f(L(1),L(2))=min(L(1),L4)=min(max(L(1,1),L(1,2)),L4)
=min(max(min(L1,L2),L3),L4)
步骤3、计算系统寿命均值。
子系统A11,A12的寿命均值记为θ(1,1)(1,2)。由公式(5)可得
Figure GDA0001228169910000071
θ(1,2)=θ3=5
再由公式(10),可得到子系统A1的寿命均值
Figure GDA0001228169910000072
则系统均值为:
Figure GDA0001228169910000073
步骤4、计算寿命重要度。
由公式(13),可计算各设备寿命重要度:
设备1属于子系统A11,则设备1的重要度为:
Figure GDA0001228169910000074
设备2属于子系统A11,则设备2的重要度为:
Figure GDA0001228169910000081
设备3属于子系统A12,则设备3的重要度为:
Figure GDA0001228169910000082
设备4的重要度为:
Figure GDA0001228169910000083
所以
I4>I3>I2>I1,系统中设备4的寿命重要度最高。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (3)

1.一种串并联混合航天装备系统寿命重要度确定方法,其特征在于,包括如下步骤:
步骤1、根据系统中设备之间的串并联关系,将系统从上到下进行分解,得到具有多个层级的逻辑关系树;所述逻辑关系树的每一层级中,均为独立的串联或并联关系的子系统,子系统中包括至少2个互为串联或并联关系的成员;逻辑关系树中除最高层级的任意层级X的子系统为上一层级X-1中一个子系统的成员;
步骤2、针对逻辑关系树的各个层级中子系统,基于子系统与成员之间的串联或并联关系,根据子系统中各个成员的寿命,确定子系统的寿命,由此得到子系统寿命与成员寿命之间的模型;
步骤3、首先,从所述逻辑关系树的最底层开始:
A、根据步骤2获得的子系统寿命与成员寿命之间的模型,结合成员的寿命分布函数基于子系统与成员之间的串联或并联关系,获得用成员寿命分布函数表征的各个子系统的寿命分布函数;
B、再根据子系统的寿命分布函数获得各个子系统的寿命均值;
然后,遍历最底层的各个子系统,得到子系统的寿命均值后,将子系统的寿命均值作为上一层级中所在子系统成员的寿命,采用步骤A和步骤B的方法,获得该层级子系统的寿命均值;依次类推,直到得到最高层级的寿命均值,即系统的寿命均值;
步骤4、用系统寿命均值对任意成员寿命求偏导数,即为该成员的寿命重要度;
其中,所述步骤3中,根据子系统寿命与成员寿命之间的模型,结合成员的寿命分布函数,基于子系统与成员之间的并联关系,获得用成员寿命分布函数表征的各个子系统的寿命分布函数如下:
Figure FDA0002632523380000011
子系统寿命为L,系统由n个成员构成,成员的寿命为L1,L2,...,Ln;m1,...,mn都可以取1,2,…,n中的任一个值,这组下标用于表示k,k=1,2,…,n个不同设备;θi表示第i个设备的寿命均值;
Figure FDA0002632523380000012
表示第m1个设备的寿命均值;
Figure FDA0002632523380000013
表示第mk个设备的寿命均值;
然后得到子系统寿命均值:
Figure FDA0002632523380000021
所述步骤3中,根据子系统寿命与成员寿命之间的模型,结合成员的寿命分布函数,基于子系统与成员之间的串联关系,获得用成员寿命分布函数表征的各个子系统的寿命分布函数如下:
Figure FDA0002632523380000022
则子系统寿命均值为:
Figure FDA0002632523380000023
遍历最底层的各个子系统,得到子系统的寿命均值后,将子系统的寿命均值作为上一层级中所在子系统成员的寿命,获得该层级子系统的寿命均值;依次类推,直到得到最高层级的寿命均值,即系统的寿命均值,具体过程为:
单一串联子系统或并联子系统寿命均值按公式(5)或公式(10)计算,统一记为θ=g(θ1,...,θn);子系统A1,…,As寿命均值记为θ(1),θ(2),...,θ(s),子系统寿命均值可由再下一层的子系统或者设备寿命均值按照公式(5)或(10)计算获得;则串并联混合系统寿命均值为:
θ=g(θ(1),θ(2),...,θ(s)) (11)
所述步骤4中,计算设备级成员的寿命重要度时,先求系统寿命对子系统偏导数,再乘上子系统寿命对成员寿命的偏导数。
2.如权利要求1所述的一种串并联混合航天装备系统寿命重要度确定方法,其特征在于,所述步骤2中,当子系统的成员之间为并联关系时,子系统寿命为成员中寿命最大值;当子系统的成员之间为串联关系时,子系统寿命为成员中寿命最小值。
3.如权利要求1所述的一种串并联混合航天装备系统寿命重要度确定方法,其特征在于,所述步骤3中,组成串联或并联关系子系统的单个成员的寿命分布函数服从指数分布。
CN201611103084.7A 2016-12-05 2016-12-05 一种串并联混合航天装备系统寿命重要度确定方法 Expired - Fee Related CN106682273B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611103084.7A CN106682273B (zh) 2016-12-05 2016-12-05 一种串并联混合航天装备系统寿命重要度确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611103084.7A CN106682273B (zh) 2016-12-05 2016-12-05 一种串并联混合航天装备系统寿命重要度确定方法

Publications (2)

Publication Number Publication Date
CN106682273A CN106682273A (zh) 2017-05-17
CN106682273B true CN106682273B (zh) 2021-04-13

Family

ID=58867523

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611103084.7A Expired - Fee Related CN106682273B (zh) 2016-12-05 2016-12-05 一种串并联混合航天装备系统寿命重要度确定方法

Country Status (1)

Country Link
CN (1) CN106682273B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109726492B (zh) * 2019-01-03 2023-04-18 中国人民解放军海军航空大学 一种无失效数据的机载多余度电气线路互联系统可靠性分析方法
CN111553059B (zh) * 2020-04-16 2022-11-01 北京航空航天大学 一种任务健康和风险健康耦合的卫星健康度评估方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102722620A (zh) * 2012-06-06 2012-10-10 中国矿业大学 一种实现寿命均衡的产品可靠性分配方法
CN104361169A (zh) * 2014-11-12 2015-02-18 武汉科技大学 一种基于分解法建模的可靠性监测方法
CN104899646A (zh) * 2015-05-14 2015-09-09 电子科技大学 一种多设备混联系统的预测性维护方法
CN106021838A (zh) * 2015-12-30 2016-10-12 西北工业大学 一种复杂电子系统剩余寿命预测方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120178066A1 (en) * 2011-01-11 2012-07-12 Ben Richard Drum Accelerated Reading Loop Bible, Reading Method and System to Encourage Bible Reading & Familiarity And Book Edge Coloring Technique
CN104217104A (zh) * 2014-08-19 2014-12-17 上海交通大学 一种基于风险评估的电力变压器寿命分析方法及系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102722620A (zh) * 2012-06-06 2012-10-10 中国矿业大学 一种实现寿命均衡的产品可靠性分配方法
CN104361169A (zh) * 2014-11-12 2015-02-18 武汉科技大学 一种基于分解法建模的可靠性监测方法
CN104899646A (zh) * 2015-05-14 2015-09-09 电子科技大学 一种多设备混联系统的预测性维护方法
CN106021838A (zh) * 2015-12-30 2016-10-12 西北工业大学 一种复杂电子系统剩余寿命预测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"复杂系统可靠性工程相关理论及技术研究";孔伟;《中国优秀博硕士学位论文全文数据库(硕士) 工程科技II辑》;20030615(第2期);正文第10-11页第2节,第18-25页第3节,第34-37页第4节 *
"面向维修过程的多态混联系统综合重要度计算方法";涂慧玲等;《自动化学报》;20140115;第40卷(第1期);正文第126-134页第3节 *

Also Published As

Publication number Publication date
CN106682273A (zh) 2017-05-17

Similar Documents

Publication Publication Date Title
CN112186766B (zh) 一种基于Benders分解的多阶段分解方法
CN107977744B (zh) 一种基于传统Benders分解法的电力系统日前鲁棒调度方法
CN107451686B (zh) 考虑随机预测误差的遗传算法的微电网能源优化方法
CN106682273B (zh) 一种串并联混合航天装备系统寿命重要度确定方法
CN109193807B (zh) 一种电力系统经济调度方法及系统
CN110533304B (zh) 一种电力系统负荷不确定性分析方法
CN110442143B (zh) 一种基于组合多目标鸽群优化的无人机态势数据聚类方法
CN109242290B (zh) 一种无人机群行动方案自动生成方法
CN112488371B (zh) 一种基于大数据的管道智能预警方法和系统
CN108388745B (zh) 基于分布并行局部优化参数的最小二乘支持向量机软测量方法
CN108108883B (zh) 一种基于聚类算法的车辆调度网络弹性简化方法
CN110619148A (zh) 基于区间灰数的装备adc效能评估方法
CN105976127B (zh) 综合电力推进系统风险评估方法及装置
CN112380763A (zh) 一种基于数据挖掘的堆内构件可靠性分析系统及方法
Chan et al. An approach to optimally allocate energy storage in naval electric ships
CN111027760A (zh) 一种基于最小二乘向量机的电力负荷预测方法
CN116187801A (zh) 无人体系中表征任务可持续性的体系区间可用度确定方法
Hu et al. Availability equivalence analysis of a repairable series-parallel system
CN111080023A (zh) 一种基于头脑风暴优化的舰载飞机出动调度规划方法
CN105469644B (zh) 飞行冲突解脱方法及设备
MAndelli et al. Mining nuclear transient data through symbolic conversion
CN108268982B (zh) 一种大规模有源配电网分解策略评价方法和装置
CN108898890B (zh) 空管运行效率等级评估方法及其装置
CN109617123B (zh) 基于状态空间组合和聚类简化的风火系统的可靠性灵敏度分析方法
CN109726873B (zh) 一种考虑信息物理耦合的通信节点检修优先级排序方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210413

Termination date: 20211205

CF01 Termination of patent right due to non-payment of annual fee