[發(fā)明專利]大型貨輪整體換件維修多Weibull單元部件備件需求量確定方法

在審

• 申請?zhí)枺?/td>: 201811360545.8
• 申請日: 2018-11-15
• 公開（公告）號: CN109492974A
• 公開（公告）日: 2019-03-19
• 發(fā)明（設(shè)計(jì)）人: 趙中元
• 申請（專利權(quán)）人: 青島航訊科技有限公司
• 主分類號: G06Q10/08
• 分類號: G06Q10/08
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.大型貨輪整體換件維修多Weibull單元部件備件需求量確定方法，包括建立Gamma等效部件步驟、設(shè)置初值步驟、計(jì)算保障概率步驟和判斷步驟，其特征在于：

(1)建立Gamma等效部件步驟，包含下述子步驟：

(1.1)建立單元壽命矩陣：建立L行、N列的單元壽命矩陣M，M中的每個(gè)元素隨機(jī)產(chǎn)生，服從W(a0,b0)分布，a0、b0分別為W(a0,b0)的尺度參數(shù)和形狀參數(shù)，由單元可靠性壽命試驗(yàn)數(shù)據(jù)進(jìn)行Weibull分布擬合計(jì)算得到；L為模擬次數(shù)，1000≤L≤10000；N為單元裝機(jī)數(shù)，由部件特性確定；

(1.2)提取部件首次故障壽命數(shù)組：對單元壽命矩陣M中的每行元素，從小到大排序，得到新壽命矩陣M1；取新矩陣M1中第(N-k+1)列數(shù)據(jù)，令其為T1；其中，k為保證部件正常工作的最少單元數(shù)量，由部件特性確定，1≤k≤N；

(1.3)建立Gamma等效部件：

對T1進(jìn)行Gamma分布擬合計(jì)算，得到G(a,b)，作為該部件的Gamma等效部件，a、b分別為G(a,b)的形狀參數(shù)和尺度參數(shù)；

(2)設(shè)置初值步驟：

設(shè)置保障概率目標(biāo)值P0＝0.5～0.99，置備件需求變量j＝0；

(3)計(jì)算保障概率步驟：

計(jì)算保障概率P：

其中，j為備件需求變量，t為保障任務(wù)時(shí)間，x為時(shí)間變量；

(4)判斷步驟：

判斷是否P≥P0，是則得到備件需求量j×N；否則j＝j(luò)+1，轉(zhuǎn)步驟(3)。

2.如權(quán)利要求1所述的多Weibull單元部件備件需求量確定方法，其特征在于：

所述建立單元壽命矩陣子步驟中，使用Matlab中的wblrnd函數(shù)建立單元壽命矩陣M：M＝wblrnd(a0,b0,L,N)；a0、b0分別為W(a0,b0)的尺度參數(shù)和形狀參數(shù)，由單元可靠性壽命試驗(yàn)數(shù)據(jù)進(jìn)行Weibull分布擬合計(jì)算得到；L為模擬次數(shù)，1000≤L≤10000；N為單元裝機(jī)數(shù)，由部件特性確定。

3.如權(quán)利要求1所述的多Weibull單元部件備件需求量確定方法，其特征在于：

所述建立Gamma等效部件子步驟中，使用Matlab中的gamfit函數(shù)，對T1進(jìn)行Gamma分布擬合計(jì)算，得到變量phat：phat＝gamfit(T1)；從而得到該部件的Gamma等效部件G(a,b)；其中，a＝phat(1)，b＝phat(2)。

4.如權(quán)利要求1所述的多Weibull單元部件備件需求量確定方法，其特征在于：

所述計(jì)算保障概率步驟中，利用Matlab中的gamcdf函數(shù)，計(jì)算保障概率P：P＝1-gamcdf(t,a(j+1),b)。

5.如權(quán)利要求1所述的多Weibull單元部件備件需求量確定方法，其特征在于：

所述計(jì)算保障概率步驟中，可用度作為保障指標(biāo)時(shí)，所述保障概率P修改為：

其中，R(x)為G(a×(j+1),b)在x時(shí)刻的可靠度，j為備件需求變量，t為保障任務(wù)時(shí)間，x、y為時(shí)間變量。

6.如權(quán)利要求5所述的多Weibull單元部件備件需求量確定方法，其特征在于：

所述計(jì)算保障概率步驟中，先利用Matlab中的gamcdf函數(shù)，計(jì)算可靠度R(x)：R(x)＝1-gamcdf(x,a(j+1),b)，再利用Matlab中的quad積分函數(shù)來計(jì)算保障概率P，P＝quad(@R(x),0,t)/t。