CN107220216B - 一种利用特征数的威布尔型备件需求量的近似计算方法 - Google Patents
一种利用特征数的威布尔型备件需求量的近似计算方法 Download PDFInfo
- Publication number
- CN107220216B CN107220216B CN201710343600.1A CN201710343600A CN107220216B CN 107220216 B CN107220216 B CN 107220216B CN 201710343600 A CN201710343600 A CN 201710343600A CN 107220216 B CN107220216 B CN 107220216B
- Authority
- CN
- China
- Prior art keywords
- distribution
- skewness
- weibull
- kurtosis
- parameters
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Evolutionary Biology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种利用特征数的威布尔型备件需求量的近似计算方法,该近似计算方法主要包括如下步骤:(1)利用所述威布尔分布参数α、b计算伽玛分布的参数αg、λ,计算正态分布的参数μ、σ;(2)计算偏度和峰度,依据所述伽玛分布的参数αg、λ及所述正态分布的参数μ、σ,按下式计算三种分布的所述特征数中的偏度和峰度;(3)比较所述特征数中的偏度和峰度,若所得伽玛分布与威布尔分布的偏度绝对差值较小,则按照一种方法计算备件保障概率;否则,按照另外的方法计算备件保障概率Ps。按照本发明实现的备件需求量的近似计算方法,能够简化计算过程以及提高近似计算的精度。
Description
技术领域
本发明属于备件需求量计算领域,特别是涉及一种利用特征数的威布尔型备件需求量的近似计算方法。
背景技术
威布尔分布常用来描述因逐渐老化导致故障的元器件寿命,具有这种威布尔分布的元器件为威布尔型单元。威布尔型单元主要适用于机电件,如:滚珠轴承、继电器、开关、断路器、某些电容器、电子管、磁控管、电位计、陀螺、电动机、航空发电机、蓄电池、液压泵、空气涡轮发动机、齿轮、活门、材料疲劳件等。
在上述单元使用于各类系统中时,需要对其备件的需求量进行预先评估计算,备件是在考虑备件寿命的情况下保障装备可持续工作的物质条件,在理论上,备件需求量计算涉及多重卷积。由于威布尔分布的多重卷积形式极为复杂,以致难以获得其多重卷积的数值积分结果。因此,在工程上,一般都采用近似方法来计算威布尔型备件需求量(例如指数近似、正态近似),但目前在工程中使用的近似方法误差较大,其中指数近似是在威布尔形状参数接近于1的时候计算效果好,正态近似只在威布尔形状参数大于3,并且还要保证近似计算方法合理才能达到较好的计算效果,上述计算方式不仅计算过程复杂,并且不能有效覆盖形状参数的可能取值范围的所有情况,使得不能执行有效的备件需求量的计算。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种利用特征数的威布尔型备件需求量的近似计算方法,利用威布尔型单元的寿命分布参数来计算伽玛分布和正态分布各自的参数,再分别计算这三种分布各自的偏度和峰度,从中选择和威布尔分布的偏度和峰度更接近的一种分布(伽玛或正态),用于近似描述该威布尔型单元寿命,并以此计算备件需求量。
为实现上述目的,按照本发明,提供一种利用特征数的威布尔型备件需求量的近似计算方法,所述威布尔型备件的寿命服从威布尔分布W(α,b),α、b为威布尔分布参数,α为尺度参数,b为形状参数;所述特征数为均值、方差、偏度和峰度,其特征在于,该计算方法包括如下步骤:
步骤一:利用所述威布尔分布参数α、b计算伽玛分布的参数αg、λ,
计算正态分布的参数μ、σ,
当正态分布的参数为μ、σ时,其均值为μ,方差为σ2;按照所求正态分布的均值和方差,与威布尔分布的均值和方差相等的原则,计算出μ、σ,
步骤二:计算偏度和峰度,
依据所述伽玛分布的参数αg、λ及所述正态分布的参数μ、σ,按下式计算三种分布的所述特征数中的偏度和峰度:
正态分布,偏度为0,峰度为0;
步骤三:按如下规则比较所述特征数中的偏度和峰度,
判断所述步骤二中的伽玛分布与正态分布中偏度和峰度与所述威布尔分布的偏度和峰度的绝对差值情况;
步骤三:计算备件需求量,
以所述偏度特征数的绝对差值为首要比较条件,相比正态分布,若所述步骤三所得伽玛分布与所述威布尔分布的偏度绝对差值较小,则按照下式计算备件保障概率:
否则,按照下式计算备件保障概率Ps:
其中,Tw为保障任务时间,所述保障任务时间为所述备件完成任务的预期累积工作时间;
设置所述备件保障概率阈值,令j从0开始逐一递增,使得所述保障概率Ps大于或等于所述概率阈值的j值即为计算出的备件需求量。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,具有以下有益效果:
(1)现有技术中也有近似方法,但是由于存在近似程度无法预估的问题导致算法的精度无法达到有效的目标,本发明率先提出了采用两种分布来执行近似计算,并且计算出偏度和峰度,依据两个指标上挑选近似方案,并且尤其是偏度上来进行判断,这样可以比较出近似程度,从而选取相应的近似分布执行备件需求量的计算;
(2)在现有技术的结论上,在威布尔的形状参数取极端时,一个极端趋向于伽玛,另外一个极端趋向于正态,在这种情况下,利用偏度和峰度这两个指标来选择趋向的近似,能够显著地提高计算的精确程度;
(3)按照本发明的近似方法,在近似成伽玛和正态近似分布的时候,提出了参数转化公式,按照上述的参数转化公式直接执行计算,无需迭代过程,从而节约了复杂的迭代过程;
(4)本发明还提出了参数转化公式的近似原则,即不管是伽玛或者是正态近似,务必满足在均值和方差这两个指标上与威布尔分布相等,对任何一个分布,都可以使用均值、方差、偏度和峰度来描述,本发明的计算方法中采用均值和方差相等,利用偏度和峰度作为近似计算指标的约定,从而简化了计算过程并且提高了近似的精度。
附图说明
图1为威布尔分布的形状参数为1.3情况下的三种分布的概率密度曲线对比情况;
图2为威布尔分布的形状参数为1.9情况下的三种分布的概率密度曲线对比情况;
图3为威布尔分布的形状参数为2.5情况下的三种分布的概率密度曲线对比情况;
图4为威布尔分布的形状参数为3.5情况下的三种分布的概率密度曲线对比情况。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
记随机变量X服从威布尔分布W(α,b),威布尔分布密度函数如式(1),
其中α>0为尺度参数,b>0为形状参数。
记随机变量X服从伽玛分布Ga(αg,λ),其中αg>0为形状参数,λ>0为尺度参数,伽玛分布密度函数如式(2)。
记随机变量X服从正态分布N(μ,σ2),其中μ为位置参数,μ的物理含义为寿命均值;σ为尺度参数,σ2的物理含义为寿命的方差。正态分布密度函数如式(3)。
对于寿命服从威布尔分布W(α,b)的单元,本发明计算备件需求量的步骤如下:
1)计算伽玛分布的参数αg、λ
按照所求伽玛分布的均值和方差,与威布尔分布的均值和方差相等的原则,按式(4)计算αg、λ,
2)计算正态分布的参数μ、σ
按照所求正态分布的均值和方差,与威布尔分布的均值和方差相等的原则,按式(5)计算μ、σ
3)计算偏度和峰度
已知伽玛分布的参数αg、λ,正态分布的参数μ、σ和威布尔分布α、b参数,按表1计算三种分布的偏度和峰度。
表1偏度和峰度的计算式
4)计算备件需求量
4.1)在偏度和峰度这两个特征数上,相比正态分布,若步骤3)所得伽玛分布更接近威布尔分布的偏度和峰度(以偏度的误差绝对值指标为主),则按照式(6)计算备件保障概率;否则,按照式(7)计算备件保障概率Ps
式(6)、(7)中,Tw为保障任务时间。
4.2)设置所述单元保障概率阈值,令j从0开始逐一递增,使得所述保障概率Ps大于或等于所述概率阈值的j值即为所计算出的备件需求量。
描述随机变量的常见工具除了分布函数、概率密度函数外,也可用均值、方差、偏度和峰度这4种特征数来进行描述。这4种特征数由该分布的1~4阶矩决定。在大部分情况下,只要知道1阶矩到4阶矩就已经足够描述分布的情况。因此,通过比较两种分布在均值、方差、偏度和峰度这4种特征数的差异程度,可以了解二者的相似程度。
由于本发明的方法采用“令伽玛/正态分布的均值和方差,与原威布尔分布的均值和方差都相等”的原则,因此首先保证了所求得的伽玛/正态分布与原威布尔分布具有一定的相似性,再通过进一步比较偏度和峰度,从中选择更为相似的近似分布结果。表3列出了威布尔形状参数1.1~4.1范围内,按照“令伽玛/正态分布的均值和方差,与原威布尔分布的均值和方差都相等”原则,计算伽玛和正态分布参数后,这三种分布的偏度和峰度情况。
表3偏度和峰度结果
从表2可以看出,当威布尔分布的形状参数在2.1以内时,伽玛分布比正态分布更接近威布尔分布。图1~图4展示了威布尔分布的形状参数4种典型取值时,三种分布的概率密度曲线对比情况,与上述结论相符。
为了解释上述算法的准确性,本实施例1执行上述方法来进行备件需求量的计算,并且利用以下备件保障仿真模型开展仿真验证。
对于某个不可修单元,配置n个备件,该类单元的寿命服从威布尔分布W(α,b),保障任务时间记为Tw,则模拟一次备件保障的过程如下:
(1)产生1+n个随机数ti(1≤i≤1+n),随机数ti服从威布尔分布W(α,b);
(3)当simT≥Tw时,保障任务成功,输出结果flag=1;否则保障任务失败,输出结果flag=0。
多次重复运行上述备件保障仿真模型,对所有模拟结果flag进行统计,flag均值即备件保障概率。
实施例1:某单元寿命服从威布尔分布W(200,1.3),保障任务时间为1000h,要求备件保障概率不小于0.8,计算备件需求量。
1)计算伽玛分布的参数αg、λ
按照所求伽玛分布的均值和方差,与威布尔分布的均值和方差相等的原则,计算αg、λ
2)计算正态分布的参数μ、σ
按照所求正态分布的均值和方差,与威布尔分布的均值和方差相等的原则,计算μ、σ
3)计算偏度和峰度
伽玛分布Ga(1.66,0.009)的偏度、峰度为1.55、3.61;
正态分布N(184.7,143.32)的偏度、峰度为0、0;
威布尔分布W(200,1.3)的偏度、峰度为1.35、2.43。
4)计算备件需求量
4.1)经比较三种分布的偏度和峰度,相比正态分布,伽玛分布Ga(1.66,0.009)与威布尔分布W(200,1.3)的相似程度更高,因此按下式计算备件保障概率:
4.2)令j从0开始逐一递增,使得所述保障概率Ps大于或等于所述概率阈值的j值即为所计算出的备件需求量。计算过程中的结果如表4.
表4计算过程的结果
从表4可知,算例1的备件需求量为7。
实施例2:某单元寿命服从威布尔分布W(200,2.9),保障任务时间为1000h,要求备件保障概率不小于0.8,计算备件需求量。
1)计算伽玛分布的参数αg、λ
按照所求伽玛分布的均值和方差,与威布尔分布的均值和方差相等的原则,计算αg、λ
2)计算正态分布的参数μ、σ
按照所求正态分布的均值和方差,与威布尔分布的均值和方差相等的原则,计算μ、σ
3)计算偏度和峰度
伽玛分布Ga(7.12,0.04)的偏度、峰度为0.75、0.84;
正态分布N(178.3,66.82)的偏度、峰度为0、0;
威布尔分布W(200,2.9)的偏度、峰度为0.20、-0.26。
4)计算备件需求量
4.1)经比较三种分布的偏度和峰度,相比伽玛分布,正态分布N(178.3,66.82)与威布尔分布W(200,2.9)的相似程度更高,因此按下式计算备件保障概率:
4.2)令j从0开始逐一递增,使得所述保障概率Ps大于或等于所述概率阈值的j值即为所计算出的备件需求量。计算过程中的结果如表5.
表5计算过程的结果
从表5可知,算例2的备件需求量为6。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种利用特征数的威布尔型备件需求量的近似计算方法,所述威布尔型备件的寿命服从威布尔分布W(α,b),α、b为威布尔分布参数,α为尺度参数,b为形状参数;所述特征数为均值、方差、偏度和峰度,其特征在于,该计算方法包括如下步骤:
步骤一:利用所述威布尔分布参数α、b计算伽玛分布的参数αg、λ,
计算正态分布的参数μ、σ,
当正态分布的参数为μ、σ时,其均值为μ,方差为σ2;按照所求正态分布的均值和方差,与威布尔分布的均值和方差相等的原则,计算出μ、σ,
步骤二:计算偏度和峰度,
依据所述伽玛分布的参数αg、λ及所述正态分布的参数μ、σ,按下式计算三种分布的所述特征数中的偏度和峰度:
正态分布,偏度为0,峰度为0;
步骤三:按如下规则比较所述特征数中的偏度和峰度,
判断所述步骤二中的伽玛分布与正态分布中偏度和峰度与所述威布尔分布的偏度和峰度的绝对差值情况;
步骤三:计算备件需求量,
以所述偏度特征数的绝对差值为首要比较条件,相比正态分布,若所述步骤三所得伽玛分布与所述威布尔分布的偏度绝对差值较小,则按照式(1)计算备件保障概率;否则,按照式(2)计算备件保障概率Ps:
式(1)、(2)中,Tw为保障任务时间,所述保障任务时间为所述备件完成任务的预期累积工作时间;
设置所述备件保障概率阈值,令j从0开始逐一递增,使得所述保障概率Ps大于或等于所述概率阈值的j值即为计算出的备件需求量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710343600.1A CN107220216B (zh) | 2017-05-16 | 2017-05-16 | 一种利用特征数的威布尔型备件需求量的近似计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710343600.1A CN107220216B (zh) | 2017-05-16 | 2017-05-16 | 一种利用特征数的威布尔型备件需求量的近似计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107220216A CN107220216A (zh) | 2017-09-29 |
CN107220216B true CN107220216B (zh) | 2020-07-28 |
Family
ID=59945251
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710343600.1A Active CN107220216B (zh) | 2017-05-16 | 2017-05-16 | 一种利用特征数的威布尔型备件需求量的近似计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107220216B (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108509390B (zh) * | 2018-03-25 | 2021-11-02 | 中国人民解放军海军工程大学 | 一种伽玛型串联部件的备件需求量计算方法 |
CN108509389B (zh) * | 2018-03-25 | 2022-01-14 | 中国人民解放军海军工程大学 | 一种威布尔型串联部件的备件需求量计算方法 |
CN108470250B (zh) * | 2018-03-25 | 2021-11-02 | 中国人民解放军海军工程大学 | 一种正态型串联部件的备件需求量计算方法 |
CN108460495B (zh) * | 2018-03-28 | 2020-12-18 | 联想(北京)有限公司 | 一种数据处理方法及装置 |
CN108564270B (zh) * | 2018-04-09 | 2021-11-02 | 中国人民解放军海军工程大学 | 一种贮存失效风险下的伽玛型单元备件需求量计算方法 |
CN108334720B (zh) * | 2018-04-09 | 2021-08-31 | 中国人民解放军海军工程大学 | 一种贮存失效风险下的正态型单元备件需求量计算方法 |
CN108520145B (zh) * | 2018-04-09 | 2022-04-05 | 中国人民解放军海军工程大学 | 一种贮存失效风险下的威布尔型单元备件需求量计算方法 |
CN108932358B (zh) * | 2018-04-16 | 2021-08-13 | 中国人民解放军海军工程大学 | 一种确定威布尔型单元后续备件需求量的方法 |
CN109145502B (zh) * | 2018-09-17 | 2023-05-12 | 中国人民解放军海军工程大学 | 一种威布尔型单元寿命分布参数估计方法 |
CN109543276B (zh) * | 2018-11-15 | 2022-11-01 | 青岛航讯网络技术服务有限公司 | 大型货轮长期保障任务多Gamma单元备件需求量确定方法 |
CN110737865B (zh) * | 2019-09-09 | 2023-05-19 | 中国人民解放军海军工程大学 | 一种基于分段函数的威布尔型备件需求量确定方法及系统 |
CN110704884B (zh) * | 2019-09-30 | 2023-05-12 | 中国人民解放军海军工程大学 | 一种对数正态型通用件备件需求量计算方法及装置 |
CN110717265B (zh) * | 2019-09-30 | 2022-06-14 | 中国人民解放军海军工程大学 | 一种对数正态型单元的备件需求量近似计算方法及装置 |
EP3916495A1 (de) | 2020-05-28 | 2021-12-01 | Carl Zeiss Industrielle Messtechnik GmbH | Verfahren und anordnung zum auswerten oder zum vorbereiten einer auswertung einer stichprobe von messdaten aus einer vermessung einer vielzahl von werkstücken |
CN114529018B (zh) * | 2022-01-14 | 2024-05-31 | 中国人民解放军海军工程大学 | 一种基于Gamma分布的舰船备件需求近似计算方法 |
CN116579573B (zh) * | 2023-05-23 | 2024-02-09 | 中国人民解放军海军工程大学 | 一种满足维修需求的机电件备件保障分析方法 |
CN116611654B (zh) * | 2023-05-23 | 2024-02-09 | 中国人民解放军海军工程大学 | 机电通用件在维修工作下的备件需求量估计方法及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1363211A3 (en) * | 2002-05-13 | 2004-06-23 | Honda Giken Kogyo Kabushiki Kaisha | System for predicting a demand for repair parts |
-
2017
- 2017-05-16 CN CN201710343600.1A patent/CN107220216B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1363211A3 (en) * | 2002-05-13 | 2004-06-23 | Honda Giken Kogyo Kabushiki Kaisha | System for predicting a demand for repair parts |
Non-Patent Citations (3)
Title |
---|
可靠性试验数据分布类型的综合判别法;白和丰 等;《武汉工业大学学报》;19920330;第14卷(第1期);第96-97页 * |
备件保障度评估与备件需求量模型研究;张建军 等;《电子产品可靠性与环境试验》;20041230(第6期);第18-22页 * |
威布尔型可修备件需求量的解析算法研究;刘天华 等;《系统工程与电子技术》;20120530;第34卷(第5期);第966-972页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107220216A (zh) | 2017-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107220216B (zh) | 一种利用特征数的威布尔型备件需求量的近似计算方法 | |
CN106971048B (zh) | 一种伽马型有寿件备件需求量计算方法 | |
CN106777819B (zh) | 一种正态型有寿件更换周期的计算方法 | |
CN108509389A (zh) | 一种威布尔型串联部件的备件需求量计算方法 | |
CN109165425B (zh) | 一种齿轮接触疲劳可靠性分析方法 | |
CN108520145B (zh) | 一种贮存失效风险下的威布尔型单元备件需求量计算方法 | |
CN102682348B (zh) | 复杂装备部件维修级别优化系统及其建立方法 | |
CN112884237A (zh) | 配电网预测辅助状态估计方法及系统 | |
Hong et al. | A tool for evaluating time-varying-stress accelerated life test plans with log-location-scale distributions | |
CN107392352B (zh) | 一种基于融合极限学习机的电池未来温度预测方法及系统 | |
CN109992880B (zh) | 一种燃气轮机性能适配方法、设备以及系统 | |
CN109492974A (zh) | 大型货轮整体换件维修多Weibull单元部件备件需求量确定方法 | |
Shayanfar et al. | An adaptive importance sampling-based algorithm using the first-order method for structural reliability | |
CN111931395A (zh) | 一种降低应变场重构误差的传感器测点优化方法 | |
CN108932358B (zh) | 一种确定威布尔型单元后续备件需求量的方法 | |
US7369974B2 (en) | Polynomial generation method for circuit modeling | |
CN117132144A (zh) | 一种考虑维修耗时的机械单元使用可用度估计方法和系统 | |
CN117093836A (zh) | 一种基于vmd-ewoa-lstm的末级再热器汽温预警方法及系统 | |
CN109583779B (zh) | 一种电-气互联系统的薄弱元件确定方法 | |
CN111350638A (zh) | 计算风力发电机组的发电量损失的方法及其装置 | |
CN112507516B (zh) | 一种基于可靠性的电气设备预防性维修优化方法及装置 | |
CN110826234B (zh) | 一种基于仿真的多应力加速寿命试验方案优化方法 | |
CN109145502B (zh) | 一种威布尔型单元寿命分布参数估计方法 | |
CN108038317B (zh) | 精密仪器性能参数保持期预测方法和系统 | |
CN108595767B (zh) | 一种基于可靠性的海洋立管viv疲劳安全系数的确定方法 |
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 |