大型货轮长期保障任务多正态单元部件备件需求量确定方法技术

大型货轮长期保障任务多正态单元部件备件需求量确定方法，属于机械设备保障性的设计方法，解决机械设备保障性设计过程中备件需求量计算困难、不准确的问题。本发明专利技术以可用度为保障要求，包括设置初值步骤、建立Gamma等效部件步骤、计算保障概率步骤和判断步骤；本发明专利技术针对正态单元的可靠性特性、立足多单元k/N表决部件的真实可靠性结构，进行快速计算备件数量，解决现有机械设备保障性设计中备件需求量计算不准确的问题。在机械设备的保障性设计阶段和列装运行阶段，能准确确定备件需求量进而能较为准确地估计全寿命周期机械设备相关费用，为评估机械设备费效比状态、实现机械设备“好保障、保障好”提供技术支撑。

【技术实现步骤摘要】
大型货轮长期保障任务多正态单元部件备件需求量确定方法
本专利技术属于机械设备保障性的设计方法，特别涉及一种以可用度为保障要求的大型货轮长期保障任务多正态单元部件备件需求量确定方法。
技术介绍
备件需求量是机械设备保障性中的一种具体属性，备件所对应的费用是全寿命期间机械设备费用的重要组成部分。长期保障任务在本文中指机械设备执行任务的时间与单元平均寿命的比例较大，本文中该比例不小于4。准确计算长期保障任务下的备件需求量，不论是对机械设备的保障性设计还是机械设备有效运行都具有重要的意义。在机械设备的层次结构中：机械设备由多种部件组成，部件由N个单元构成(N≥1)；单元与单元之间以串联或某种可靠性结构形式进行连接，构成部件，该部件称之为单元装机数为N的部件，简称多单元部件。正态单元指该单元的寿命服从正态分布，记为N(a1,b1),a1、b1分别为N(a1,b1)的寿命均值和标准差，其失效密度函数对单元或部件进行可靠性寿命试验，利用寿命试验数据进行正态分布拟合计算，得到a1、b1。正态寿命件主要适用于机械件，如：汇流环、齿轮箱、减速器等。多正态单元部件指由多个正态单元组成的部件。k/N表决部件是由N个单元组成的多单元部件，其中至少有k个单元正常工作，部件才正常工作，1≤k≤N，记为k/N(G)，它的可靠性框图如图1所示；N/N(G)为串联部件，1/N(G)为并联部件。与串联、并联结构相比，表决结构是一种更为一般、普遍的可靠性结构形式。表决结构常见于两种情况：一种在机械设备的可靠性设计中，针对单元可靠度达不到要求的情况，采用表决结构作为一种冗余设计来提高部件的整体可靠度；另一种情况，则是在部件出现劣化或部分故障导致性能下降后，通过定义k来决定所能接受的部件最低性能，小于k则认为部件整体性能已不能满足工作要求。备件需求量与维修模式紧密相关。换件维修主要针对不可修单元，或由于条件限制只能在现场进行换件维修的情况----例如大型货轮在海上运输期间，由于受维修条件、维修能力等限制，对出现故障的单元，船员只能采取换件维修的方法排除故障、恢复机械设备性能。换件维修模式下，单元装机数为N的表决部件，现有确定备件需求量的方法有：(1)卷积法。在数学理论上换件维修等同于卷积计算。备件需求量为m的计算过程意味着一个(m+1)重卷积计算。对多重卷积而言，此时卷积表达式形式极其复杂、难以推导；并且，利用计算机进行多重卷积的数值求解计算，也随着备件需求量m每增加1个,其计算量以10倍以上方式递增，不仅耗时长且难以保证计算精度。(2)模拟方法。该方法以多次模拟表决部件在保障期间工作运行的方式，来统计计算保障概率要求下所需的备件需求数量。与计算备件需求量相比，该方法更适合用于备件数量确定后的保障效果评估。这是由于随机因素的影响，使得备件数量对应的保障概率同时具有数学期望值和标准差两种属性，模拟方法得到的保障概率往往在保障概率期望值附近波动，这种波动在某些情况下(例如较长任务时间)会导致备件数量出现一个较大的波动，而导致模拟方法往往不能得到稳定的备件需求量结果。(3)近似法。由一个正态单元构成的部件，其备件需求量S的近似计算公式：式中：t是保障任务时间，E是平均寿命，σ是标准差，up是正态分布分位数，可从GB/T4086.1从查出。该式不能用于多单元部件的情况。上述方法在实际应用中都存在着结果不准确、计算耗时长的问题，可行性差。
技术实现思路
本专利技术提供一种大型货轮长期保障任务多正态单元部件备件需求量确定方法，以可用度为保障要求，解决机械设备保障性设计过程中备件需求量计算困难、不准确的问题。本专利技术所提供的一种大型货轮长期保障任务多正态单元部件备件需求量确定方法，以可用度为保障要求，包括设置初值步骤、建立Gamma等效部件步骤、计算保障概率步骤和判断步骤，其特征在于：(1)设置初值步骤：设置保障概率目标值P0＝0.5～0.99，置j＝0；(2)建立Gamma等效部件步骤，包含下述子步骤：(2.1)建立单元寿命数组T0：建立包含N个元素的数组，数组中的每个元素表示单元寿命值，随机产生并服从N(a0,b0)分布，a0、b0分别为N(a0,b0)的寿命均值和标准差，由单元可靠性寿命试验数据进行正态分布拟合计算得到；将数组中的所有元素按从小到大顺序排列，得到单元寿命数组T0：T0＝[t01t02…t0N]，其中的任意2个元素满足条件：当i＜j时，t0i≤t0j，部件单元装机数N＞1；(2.2)建立备件单元寿命矩阵：建立L行Q列的备件单元寿命矩阵M：M中的每个元素表示备件单元寿命值，随机产生并服从N(a0,b0)分布，a0、b0分别为N(a0,b0)的寿命均值和标准差，由单元可靠性寿命试验数据进行正态分布拟合计算得到；矩阵表达式中，Q＝N-k+1,k为部件正常工作所需的最小单元数，由部件特性确定，1≤k≤N；L为模拟次数，1000≤L≤10000；(2.3)建立部件寿命数组T1，包括如下过程：(2.3.1)建立包括L个元素的部件寿命数组T1，置序号变量m＝1；(2.3.2)置部件寿命数组T1中的第m个元素T1(m)＝T0(N-k+1)，T0(N-k+1)为单元寿命数组T0中第(N-k+1)个元素；然后更新单元寿命数组T0：(2.3.3)置单元寿命数组T0中的第i个元素T0(i)＝tmi,i≤Q，tmi为备件单元寿命矩阵M中的第m行第i列元素；(2.3.4)将单元寿命数组T0中的元素从小到大排序，得到再次更新的单元寿命数组T0：T0＝[t01t02…t0N]，其中的任意2个元素满足条件：当i＜j时，t0i≤t0j；(2.3.5)判断是否m>L，是则转步骤(2.4)，否则置m＝m+1，转过程(2.3.2)；(2.4)Gamma拟合：对部件寿命数组T1进行Gamma分布拟合计算，得到该部件的Gamma等效部件，记为G(a,b)，a、b分别为G(a,b)的形状参数和尺度参数；(3)计算保障概率步骤：计算保障概率P：其中，R(x)为等效Gamma部件可靠度，t为保障任务时间，j为备件需求变量，Gamma函数x、y为时间变量；(4)判断步骤：判断是否P≥P0，是则得到备件需求量j×(N-k+1)；否则j＝j+1，转步骤(3)。所述的大型货轮长期保障任务多正态单元部件备件需求量确定方法，其特征在于：所述的建立备件单元寿命矩阵步骤中，使用Matlab中normrnd函数建立备件单元寿命矩阵M：M＝normrnd(a0,b0,L,N-k+1)；a0、b0分别为N(a0,b0)的寿命均值和标准差，由单元可靠性寿命试验数据进行正态分布拟合计算得到；L为模拟次数，1000≤L≤10000；N为单元装机数，由部件特性确定。所述的大型货轮长期保障任务多正态单元部件备件需求量确定方法，其特征在于：所述Gamma拟合子步骤中，使用Matlab中的gamfit函数，对T1进行Gamma分布拟合计算，得到变量phat：phat＝gamfit(T1)；从而得到该部件的Gamma等效部件G(a,b)；其中a＝phat(1)，b＝phat(2)。所述的大型货轮长期保障任务多正态单元部件备件需求量确定方法，其特征在于：所述计算保障概率步骤中，先利用Matlab中的gamcdf函数，计算等本文档来自技高网...

【技术保护点】
1.一种大型货轮长期保障任务多正态单元部件备件需求量确定方法，以可用度为保障要求，包括设置初值步骤、建立Gamma等效部件步骤、计算保障概率步骤和判断步骤，其特征在于：(1)设置初值步骤：设置保障概率目标值P0＝0.5～0.99，置j＝0；(2)建立Gamma等效部件步骤，包含下述子步骤：(2.1)建立单元寿命数组T0：建立包含N个元素的数组，数组中的每个元素表示单元寿命值，随机产生并服从N(a0,b0)分布，a0、b0分别为N(a0,b0)的寿命均值和标准差，由单元可靠性寿命试验数据进行正态分布拟合计算得到；将数组中的所有元素按从小到大顺序排列，得到单元寿命数组T0：T0＝[t01 t02 … t0N]，其中的任意2个元素满足条件：当i＜j时，t0i≤t0j，部件单元装机数N＞1；(2.2)建立备件单元寿命矩阵：建立L行Q列的备件单元寿命矩阵M：

【技术特征摘要】
1.一种大型货轮长期保障任务多正态单元部件备件需求量确定方法，以可用度为保障要求，包括设置初值步骤、建立Gamma等效部件步骤、计算保障概率步骤和判断步骤，其特征在于：(1)设置初值步骤：设置保障概率目标值P0＝0.5～0.99，置j＝0；(2)建立Gamma等效部件步骤，包含下述子步骤：(2.1)建立单元寿命数组T0：建立包含N个元素的数组，数组中的每个元素表示单元寿命值，随机产生并服从N(a0,b0)分布，a0、b0分别为N(a0,b0)的寿命均值和标准差，由单元可靠性寿命试验数据进行正态分布拟合计算得到；将数组中的所有元素按从小到大顺序排列，得到单元寿命数组T0：T0＝[t01t02…t0N]，其中的任意2个元素满足条件：当i＜j时，t0i≤t0j，部件单元装机数N＞1；(2.2)建立备件单元寿命矩阵：建立L行Q列的备件单元寿命矩阵M：M中的每个元素表示备件单元寿命值，随机产生并服从N(a0,b0)分布，a0、b0分别为N(a0,b0)的寿命均值和标准差，由单元可靠性寿命试验数据进行正态分布拟合计算得到；矩阵表达式中，Q＝N-k+1,k为部件正常工作所需的最小单元数，由部件特性确定，1≤k≤N；L为模拟次数，1000≤L≤10000；(2.3)建立部件寿命数组T1，包括如下过程：(2.3.1)建立包括L个元素的部件寿命数组T1，置序号变量m＝1；(2.3.2)置部件寿命数组T1中的第m个元素T1(m)＝T0(N-k+1)，T0(N-k+1)为单元寿命数组T0中第(N-k+1)个元素；然后更新单元寿命数组T0：(2.3.3)置单元寿命数组T0中的第i个元素T0(i)＝tmi,i≤Q，tmi为备件单元寿命矩阵M中的第m行第i列元素；(2.3.4)将单元寿命数组T0中的元素从小到大排序，得到再次更新的单元寿命数组T0：T0＝[t01t02…t0N]，其中的任意2个元素满足条件：当i＜j...

【专利技术属性】
技术研发人员：赵中元，
申请(专利权)人：青岛航讯网络技术服务有限公司，
类型：发明
国别省市：山东,37