发明名称 一种装备研制过程备件需求量预测方法
摘要 本发明公开了一种装备研制过程备件需求量预测方法,包括以下几个步骤:步骤一:直接获取装备的可达可用度指标Aa、平均故障间隔时间MTBF、平均修复时间MTTR和平均预防性维修时间MPMT;步骤二:获取系统级备件平均等待时间WT;步骤三:对系统级备件平均等待时间WT进行分配,获取分系统等待时间WTi;步骤四:采用Dephi分配法获取组成分系统的各LRU(外场可更换单元)级产品等待时间WTij;步骤五:获取LRU备件缺货期望值EBOij、库存水平sij;步骤六:根据库存水平sij获取备件满足率EFRij;本发明建立备件等待时间与装备可用度的量化关系,始终围绕用户需求和研制要求进行备件规划工作;本发明采用等分配法与Dephi法相结合的方法,具有较大的灵活性,且能适用工程研制各个阶段。
申请公布号 CN101853448A 申请公布日期 2010.10.06
申请号 CN201010190151.X 申请日期 2010.05.25
申请人 北京航空航天大学 发明人 郭霖瀚;陶小创;康锐;王乃超
分类号 G06Q10/00(2006.01)I;G06Q50/00(2006.01)I 主分类号 G06Q10/00(2006.01)I
代理机构 北京永创新实专利事务所 11121 代理人 赵文利
主权项 1.一种装备研制过程备件需求量预测方法,其特征在于,包括以下几个步骤:步骤一:直接获取装备的可达可用度指标A<sub>a</sub>、平均故障间隔时间MTBF、平均修复时间MTTR和平均预防性维修时间MPMT;根据用户需求和装备研制要求中规定的可达可用度指标A<sub>a</sub>,从可靠性维修性设计数据中直接获取平均故障间隔时间MTBF、平均修复时间MTTR和平均预防性维修时间MPMT;步骤二:获取系统级备件平均等待时间WT;根据公式(1),获取系统级备件平均等待时间WT;<maths num="0001"><![CDATA[<math><mrow><msub><mi>A</mi><mi>a</mi></msub><mo>=</mo><mfrac><mi>MTBF</mi><mrow><mi>MTBF</mi><mo>+</mo><mi>MTTR</mi><mo>+</mo><mi>MPMT</mi><mo>+</mo><mi>WT</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤三:对系统级备件平均等待时间WT进行分配,获取分系统等待时间WT<sub>i</sub>;如果能够获取分系统平均修复时间、平均供应时间以及分系统故障对任务完成影响程度信息,则采用Dephi分配法,否则采用等分配法;等分配法具体为:根据式(2)得到各分系统等待时间WT<sub>i</sub>:<maths num="0002"><![CDATA[<math><mrow><msub><mi>d</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>WT</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>m</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>d</mi><mi>i</mi></msub><mo>&CenterDot;</mo><mi>WT</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中:d<sub>i</sub>为第i个分系统需求率,i=1,2,…,m,m为装备中分系统个数,具体为:d<sub>i</sub>=λ<sub>i</sub>+γ<sub>i</sub>                                                      (3)式中,λ<sub>i</sub>为第i个分系统的故障率,γ<sub>i</sub>为第i个分系统的更换率,在装备研制阶段,λ<sub>i</sub>和γ<sub>i</sub>从装备的可靠性和维修性设计数据直接获得;Dephi分配法具体为:根据式(4)将系统级备件平均等待时间WT分配到分系统级:<maths num="0003"><![CDATA[<math><mrow><msub><mi>WT</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>k</mi><mi>i</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>d</mi><mi>i</mi></msub></mrow><mrow><msub><mi>d</mi><mi>i</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>k</mi><mi>i</mi></msub></mrow></mfrac><mi>WT</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,k<sub>i</sub>为第i个分系统的期望因子值,第i个分系统中可修产品的比例为B<sub>i</sub>,分配时考虑的因子为任务影响因子k<sub>i1</sub>、供应时间因子k<sub>i2</sub>和维修时间因子k<sub>i3</sub>,不可修产品的比例为1-B<sub>i</sub>,分系统的相关因子取值具体为:任务影响因子k<sub>i1</sub>为:k<sub>i1</sub>与分系统级产品是否存在功能余度和该分系统故障后对系统任务完成的影响程度有关,当分系统功能余度越多因子取值越大,对任务完成影响程度越小,因子取值越大;供应时间因子k<sub>i2</sub>为:k<sub>i2</sub>取值与该分系统的平均供应时间和分系统中最大平均供应时间有关:<maths num="0004"><![CDATA[<math><mrow><msub><mi>k</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>=</mo><mfrac><msub><mi>MTTS</mi><mi>i</mi></msub><mrow><munder><mi>max</mi><mi>i</mi></munder><msub><mi>MTTS</mi><mi>i</mi></msub></mrow></mfrac><mo>&times;</mo><mn>10</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中:MTTS<sub>i</sub>为第i个分系统的平均供应时间,<img file="FSA00000143024800015.GIF" wi="232" he="67" />是分系统中的最大平均供应时间;维修时间因子k<sub>i3</sub>为:k<sub>i3</sub>取值与该分系统的平均修复时间和分系统中最大平均修复时间有关:<maths num="0005"><![CDATA[<math><mrow><msub><mi>k</mi><mrow><mi>i</mi><mn>3</mn></mrow></msub><mo>=</mo><mfrac><msub><mi>MTTR</mi><mi>i</mi></msub><mrow><munder><mi>max</mi><mi>i</mi></munder><msub><mi>MTTR</mi><mi>i</mi></msub></mrow></mfrac><mo>&times;</mo><mn>10</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中:MTTR<sub>i</sub>为第i个分系统的平均修复时间,<img file="FSA00000143024800022.GIF" wi="235" he="69" />是分系统中的最大平均修复时间;k<sub>i</sub>=B<sub>i</sub>(k<sub>i1</sub>+k<sub>i2</sub>+k<sub>i3</sub>)+(1-B<sub>i</sub>)(k<sub>i1</sub>+k<sub>i2</sub>)                              (7)步骤四:采用Dephi分配法获取组成分系统的各LRU级产品等待时间WT<sub>ij</sub>;LRU为外场可更换单元,将第i个分系统的平均等待时间WT<sub>i</sub>分配到LRU级别,具体为:①若第i个分系统的第j个LRU为可修件,则第j个LRU的平均等待时间为:<maths num="0006"><![CDATA[<math><mrow><msub><mi>WT</mi><mi>ij</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>L</mi><mi>ij</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>d</mi><mi>ij</mi></msub></mrow><mrow><msub><mi>d</mi><mi>ij</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>L</mi><mi>ij</mi></msub></mrow></mfrac><mrow><mo>(</mo><msub><mi>WT</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>B</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中:第i个分系统中包括j个LRU级产品,第j个LRU的因子L<sub>ij</sub>=L<sub>ij1</sub>+L<sub>ij2</sub>+L<sub>ij3</sub>,对于可修产品LRU相关因子为任务影响因子L<sub>ij1</sub>、供应时间因子L<sub>ij2</sub>和维修时间因子L<sub>ij3</sub>,LRU相关因子取值具体为:任务影响因子L<sub>ij1</sub>为:L<sub>ij1</sub>与LRU级产品是否存在功能余度和该LRU故障后对分系统任务完成的影响程度有关,当LRU功能余度越多因子取值越大,对任务完成影响程度越小因子取值越大;供应时间因子L<sub>ij2</sub>为:L<sub>ij2</sub>取值与该LRU的平均供应时间和LRU的最大平均供应时间有关:<maths num="0007"><![CDATA[<math><mrow><msub><mi>L</mi><mrow><mi>ij</mi><mn>2</mn></mrow></msub><mo>=</mo><mfrac><msub><mi>MTTS</mi><mi>ij</mi></msub><mrow><munder><mi>max</mi><mi>j</mi></munder><msub><mi>MTTS</mi><mi>ij</mi></msub></mrow></mfrac><mo>&times;</mo><mn>10</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中:MTTS<sub>ij</sub>为第i个LRU的平均供应时间,<img file="FSA00000143024800025.GIF" wi="240" he="72" />是LRU的最大平均供应时间;维修时间因子L<sub>ij3</sub>为:L<sub>ij3</sub>取值与该LRU的平均修复时间和LRU的最大平均修复时间有关:<maths num="0008"><![CDATA[<math><mrow><msub><mi>L</mi><mrow><mi>ij</mi><mn>3</mn></mrow></msub><mo>=</mo><mfrac><msub><mi>MTTR</mi><mi>ij</mi></msub><mrow><munder><mi>max</mi><mi>j</mi></munder><msub><mi>MTTR</mi><mi>ij</mi></msub></mrow></mfrac><mo>&times;</mo><mn>10</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中:MTTR<sub>ij</sub>为第i个LRU的平均修复时间,<img file="FSA00000143024800027.GIF" wi="243" he="73" />是LRU的最大平均修复时间;WT<sub>i</sub>·B<sub>i</sub>为第i个分系统中可修产品的等待时间;d<sub>ij</sub>为第i个分系统中第j个LRU的需求率,具体为:d<sub>ij</sub>=λ<sub>ij</sub>+γ<sub>ij</sub>                                                   (11)式中:λ<sub>ij</sub>为第i个分系统中第j个LRU的故障率,γ<sub>ij</sub>为第i个分系统中第j个LRU的更换率;②若第i个分系统的第j个LRU为不可修件,则它的平均等待时间为:<maths num="0009"><![CDATA[<math><mrow><msub><mi>WT</mi><mi>ij</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>L</mi><mi>ij</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>d</mi><mi>ij</mi></msub></mrow><mrow><msub><mi>d</mi><mi>ij</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>L</mi><mi>ij</mi></msub></mrow></mfrac><mrow><mo>(</mo><msub><mi>WT</mi><mi>i</mi></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>B</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中:第j个LRU的因子L<sub>ij</sub>=L<sub>ij1</sub>+L<sub>ij2</sub>,对于不可修产品相关因子为任务影响因子L<sub>ij1</sub>和供应时间因子L<sub>ij2</sub>,相关因子取值见步骤四;WT<sub>i</sub>·(1-B<sub>i</sub>)为第i个分系统中不可修产品的等待时间;步骤五:获取LRU备件缺货期望值EBO<sub>ij</sub>、库存水平s<sub>ij</sub>;第i个分系统的第j个LRU备件缺货期望值EBO<sub>ij</sub>通过式(13)获取:EBO<sub>ij</sub>=WT<sub>ij</sub>×d<sub>ij</sub>                                                 (13)第j个LRU库存水平s<sub>ij</sub>满足式(14)中的数量关系:s<sub>ij</sub>=OH<sub>ij</sub>+DI<sub>ij</sub>-BO<sub>ij</sub>                                              (14)式中:OH<sub>ij</sub>为第i个分系统的第j个LRU现有备件库存量,DI<sub>ij</sub>为为第i个分系统的第j个LRU来自修复和补给的备件量,BO<sub>ij</sub>为为第i个分系统的第j个LRU备件缺货量,OH<sub>ij</sub>、DI<sub>ij</sub>和BO<sub>ij</sub>均为非负的随机变量;由式(14)可知,当DI<sub>ij</sub>≤s<sub>ij</sub>-1时,OH<sub>ij</sub>≥0,则此时还存在备件库存,此时备件需求能够得到满足,第i个分系统的第j个LRU备件的缺货期望值EBO<sub>ij</sub>为:<maths num="0010"><![CDATA[<math><mrow><msub><mi>EBO</mi><mi>ij</mi></msub><mo>=</mo><mi>P</mi><mo>{</mo><msub><mi>DI</mi><mi>ij</mi></msub><mo>=</mo><msub><mi>s</mi><mi>ij</mi></msub><mo>+</mo><mn>1</mn><mo>}</mo><mo>+</mo><mn>2</mn><mi>P</mi><mo>{</mo><msub><mi>DI</mi><mi>ij</mi></msub><mo>=</mo><msub><mi>s</mi><mi>ij</mi></msub><mo>+</mo><mn>2</mn><mo>}</mo><mo>+</mo><mn>3</mn><mi>P</mi><mo>{</mo><msub><mi>DI</mi><mi>ij</mi></msub><mo>=</mo><msub><mi>s</mi><mi>ij</mi></msub><mo>+</mo><mn>3</mn><mo>}</mo><mo>+</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><msub><mi>s</mi><mi>ij</mi></msub><mo>+</mo><mn>1</mn></mrow><mo>&infin;</mo></munderover><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>s</mi><mi>ij</mi></msub><mo>)</mo></mrow><mi>P</mi><mo>{</mo><msub><mi>DI</mi><mi>ij</mi></msub><mo>=</mo><mi>x</mi><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中P表示概率,根据式(13)计算得出的EBO<sub>ij</sub>,将EBO<sub>ij</sub>代入式(15),反解出库存水平s<sub>ij</sub>;步骤六:根据库存水平s<sub>ij</sub>获取备件满足率EFR<sub>ij</sub>;将步骤五计算出的库存水平s<sub>ij</sub>代入式(16)求出第i个分系统的第j个LRU备件的满足率;EFR<sub>ij</sub>=P{DI<sub>ij</sub>=0}+P{DI<sub>ij</sub>=1}+…+P{DI<sub>ij</sub>=s-1}=P{DI<sub>ij</sub>≤s-1}       (16)备件库存水平s<sub>ij</sub>为备件需求量,EFR<sub>ij</sub>用于表征LRU备件库存水平s<sub>ij</sub>满足装备使用要求的程度。
地址 100191 北京市海淀区学院路37号