CN110288099B - Preventive maintenance spare part scheme optimization method and device - Google Patents

Abstract

The invention discloses a preventive maintenance spare part scheme optimization method and device. Wherein the method comprises the steps of categorizing at least one component having the same initial state in the same component category; acquiring at least one spare part scheme in each component category; combining any of the spare part solutions in each of the component categories into at least one overall spare part solution, and evaluating an achievement probability of the overall spare part solution; and selecting the total spare part scheme with the standard reaching probability meeting the preset conditions as the preventive maintenance spare part scheme. The invention can reasonably realize the units which need to be replaced in the components and the number of the units under the condition of not clearly describing the physical indexes of the aging state of the product.

Description

Method and device for optimizing preventive maintenance spare parts plan

technical field

The invention relates to the field of failure prediction and spare parts guarantee, in particular to a method for optimizing a preventive maintenance spare parts plan.

Background technique

Carrying out preventive maintenance work on products generally quantifies the aging state of products by detecting clear physical indicators, so that preventive maintenance can be carried out in time, and products can be replaced before they fail to avoid major accidents.

For products that do not have clear physical indicators, or have clear physical indicators, but the cost of detectors that sense these physical indicators is too high, it is difficult to carry out and promote general preventive maintenance work.

Contents of the invention

The embodiment of the present invention discloses at least one method for optimizing the preventive maintenance spare parts plan, which can reasonably realize the units and the number of units that need to be replaced in the components without clearly describing the physical indicators of the aging state of the product.

The methods include:

S100. Classifying at least one component with the same initial state in the same component category;

S200. Acquire at least one spare parts solution in each component category;

S300. Combining any one of the spare parts schemes in each of the component categories into at least one overall spare parts scheme, and evaluating the compliance probability of the overall spare parts scheme;

S400. Select the overall spare parts plan whose compliance probability meets a preset condition as the preventive maintenance spare part plan.

In some embodiments disclosed in the present invention, S100 is configured with:

S110. Obtain the state of all or part of the units of the component;

S120. Determine the initial state according to the sound state or failure state of all the units, and classify the components with the same initial state into the same component category.

In some embodiments disclosed in the present invention, S110 is configured to describe all or part of the units of the component through an array dyFlag having dyN items,

The i-th variable of dyFlag is dyFlag _i , which is used to describe the sound state or the failure state of the unit in the i-th item in the component;

S120 is configured to describe the initial state and the component category through the matrix bjState, the j-th row vector bjState(j,:) of bjState is [dyFlag nFlag _j ], and nFlag _j is used to describe the component category of the j-th category The number of parts.

In some embodiments disclosed in the present invention, S200 is configured as:

S200 is configured as:

S210, make the category number j=1;

S220, initialization scheme serial number j2=0,

A spare part scheme for components in the component category described in category j is described by an array Ns,

Ns _i is the number of spare parts of the unit described in item i 1≤i≤dyN;

S230. Invoke a component spare part solution evaluation module, and evaluate the maintenance effect of the spare part solution Ns as bjP;

S240,

Take j2=j2+1,

acquisition fee

Update matrix plan, plan(j2,:)=[j2 bjP Mnow Ns];

在

时进入S260，在

时生成下一备件方案；S250, judge

exist

When entering S260, in

generate the next spare parts plan;

S260. Take planAll(1:j2, iy1:iy2)=plan, iy1=(j-1)(3+dyN)+1, iy2=j(3+dyN), take planN _j =j2;

S270. Take j=j+1, return to S220 when j≤N1, and enter S300 when j>N1;

S300 is configured to use planAll to combine a spare part plan for the parts in each category of the category parts into at least one overall spare part plan.

In some embodiments disclosed in the present invention,

在

时进入S260，在

时生成下一备件方案，进入S251；The S250 is configured to judge

exist

When entering S260, in

When generating the next spare parts plan, enter S251;

The S250 is configured with:

S251. Take the candidate sequence number k=1,

S252. Generate a candidate solution Nt,

S253. According to the candidate solution Nt, call the parts and spare parts solution evaluation module, and evaluate the maintenance effect of the candidate solution Nt as bjPt, and the marginal benefit cost ratio

S254, take k=k+1, return to S252 when k≤dyN, enter S255 when k>dyN,

S255. Find the maximum value of p2m _k among all p2m _k , 1≤k≤dyN, record the candidate serial number Im corresponding to the maximum value of p2m _k , update Ns _Im =nFlag _j in Ns, and return to S230.

In some embodiments disclosed in the present invention, S300 is configured with:

S310. Initialize the array planID describing the overall spare parts plan, planID _j = 1, 1≤j≤N1; the planID _j of the jth item is the serial number of the spare parts plan in the component category of the jth class,

Initialize the optimization flag bjF _j =1, 1≤j≤N1 for each type of component category,

Initialize the overall scheme number i2=0;

S320. Invoking the evaluation module of the overall spare parts scheme, obtaining the probability of reaching the standard of the overall spare parts scheme as Pm, and the cost of the spare parts as Mm,

Take i2=i2+1,

Vector result(i2,:)=[Mm Pm planID] in row i2 of matrix result;

The S330 is configured with:

在

时进入S332，在

时进入S400，S331. Judgment

exist

Enter S332 when

enter S400,

S332. Take the serial number j=1 of the component category,

S333. Get IDt = planID, IDt _j = IDt _j + 1,

S334. When IDt _j > planN _j , take bjF _j =0 and p2m _j =0,

When IDt _j ≤ planN _j , the overall spare parts program evaluation module is called, and the attainment probability of the overall spare parts program is denoted as Pt _j , the cost of spare parts is denoted as Mt _j , and the marginal benefit-cost ratio

S335, set j=j+1, enter S333 when j≤N1, enter S336 when j>N1,

S336. Find the serial number corresponding to the maximum value in all p2m _j , 1≤j≤N1, which is Im,

时，取planID_Im＝planID_Im+1，Pm＝Pt_Im，Mm＝Mt_Im，保存结果到矩阵result：取i2＝i2+1，取result的第i2行向量result(i2,:)＝[Mm Pm planID]后返回S331；exist

, take planID _Im ＝planID _Im + 1, Pm＝Pt _Im , Mm＝Mt _Im , save the result to the matrix result: take i2＝i2+1, take the i2th row vector of result result(i2,:)＝[Mm Pm planID] and return to S331;

时，转S400；exist

, turn to S400;

S400 is configured to acquire the overall spare parts plan with the attainment probability Pm=P (n>M) as the preventive maintenance spare part plan according to the matrix result.

In some embodiments disclosed in the present invention, S400 is configured to draw an efficiency curve according to the matrix result;

The X-axis of the cost curve is the cost of the overall spare parts plan, and the Y-axis is the probability of reaching the standard Pm=P (n>M);

The preventive maintenance spare parts plan is obtained by querying the curve to obtain the overall spare parts plan with the attainment probability Pm=P (n>M) and the least spare part cost.

In some embodiments disclosed in the present invention, calling the component spare parts solution evaluation module to evaluate the maintenance effect of the spare parts solution is configured as:

Pick

dyP _i is the reliability of the i-th unit after the preset time, p0 and p1 are the values that the life of the i-th unit obeys the gamma distribution, lognormal distribution, normal distribution or Weibull distribution.

In some embodiments disclosed in the present invention, calling the parts and spare parts plan evaluation module to evaluate the overall spare parts plan is configured as:

S321. Initialize the category number j=1, and spare part cost Mt=0.

S322. Obtain the spare part plan, reliability bjP, and spare part cost Mnow of the spare part plan number j2=planID _j in the component category of the jth category from planAll;

Take Mt=Mt+Mnow;

1≤k≤1+nFlag_j，得到向量p1；S323. Take x1 as an increasing sequence from 0 to nFlag _j , with a step size of 1, and iteratively calculate the probability corresponding to x1 _k

1≤k≤1+nFlag _j , get the vector p1;

S324. Calculate the convolution of p1 and p2 when j>1, and assign the result to p2;

Otherwise, take p2=p1.

S325. Take j=j+1, return to S322 when j≤N1, and enter S326 when j>N1;

总体备件方案planID的备件费用Mm＝Mt。S326. Obtain the probability of reaching the standard of the overall spare parts plan planID

The spare parts cost Mm=Mt of the overall spare parts plan planID.

The embodiment of the present invention discloses at least one device for constructing a preventive maintenance spare parts plan.

The devices include:

A classification module, configured to classify at least one component with the same initial state in the same component category;

a spare parts module for obtaining at least one spare parts solution in each of said component categories;

A combination module, combining any one of the spare parts schemes in each of the component categories into at least one overall spare parts scheme, and evaluating the compliance probability of the overall spare parts scheme;

A selection module, selecting the overall spare parts plan whose compliance probability meets a preset condition as the preventive maintenance spare parts plan.

In view of the above solutions, the present invention will describe the disclosed exemplary embodiments in detail below with reference to the accompanying drawings, and also make other features and advantages of the embodiments of the present invention clear.

Description of drawings

In order to illustrate the technical solutions of the embodiments of the present invention more clearly, the accompanying drawings used in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of the present invention, and thus It should be regarded as a limitation on the scope, and those skilled in the art can also obtain other related drawings based on these drawings without creative work.

Fig. 1 is the flow chart of the method for optimizing the preventive maintenance spare parts scheme provided by the present embodiment;

Fig. 2 is the functional block diagram of the optimized preventive maintenance spare parts scheme provided by the present embodiment;

Fig. 3 is an efficiency-cost curve diagram of the spare parts solution in the example provided in this embodiment.

detailed description

In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments It is a part of embodiments of the present invention, but not all embodiments. The components of the embodiments of the invention generally described and illustrated in the figures herein may be arranged and designed in a variety of different configurations.

As used herein, the term "if" is optionally interpreted to mean "when" or "upon" or "in response to determining" or "in response to detecting," depending on the context. Similarly, the phrases "if it is determined that ..." or "if [the stated condition or event] is detected" are optionally construed to mean "when determining ..." or "in response to determining ..." depending on the context Or "upon detection of [stated condition or event]" or "in response to detection of [stated condition or event]".

In addition, the steps shown in the flow diagrams of the figures may be performed in a computer device, such as a set of computer-executable instructions, and, although a logical order is shown in the flow diagrams, in some cases, may be implemented in The steps shown or described are performed in an order different than here.

The components mentioned in this embodiment generally refer to a product of an intermediate layer structure composed of at least one unit; the units included in the component can generally be described by a reliable connection relationship.

The common reliability connection relationships between single rooms are series relationship, parallel relationship and mixed connection relationship. If all the units in the component are selected according to whether the unit is a critical component, the component can be simplified into a product composed of at least one critical component in series.

From the above, the components referred to in the following embodiments in this embodiment are products in which the reliability relationship between the units is a series relationship; the units are all critical parts that meet the standards.

Specifically, the components are composed of units described in the dyN item, and the life distribution law and unit price dyM _i of each unit are known, 1≤i≤dyN; the current maintenance time is T1; there are currently a total of N components in the same batch; In the way of replacing the faulty unit or aging unit with spare parts, preventive maintenance work is carried out on these N components at time T1; after the maintenance, the probability that the N components work until the time T2 is intact and the number n is greater than M is recorded as P(n> M), the probability is used to evaluate the maintenance effect of the preventive maintenance spare parts plan, which is called the probability of reaching the standard in this embodiment; the lower limit of P (n>M) is recorded as P1, and P1 is the requirement of the probability of reaching the standard for the plan; Then only when P(n>M)≥P1 is satisfied, the preventive maintenance spare parts plan is considered qualified.

This embodiment discloses a method for optimizing a preventive maintenance spare parts plan, which is used to construct a qualified and cost-effective preventive maintenance spare parts plan for the above-mentioned components with dyN units.

Please refer to FIG. 1 , the method of this embodiment has the following steps.

S100 , before the preventive maintenance, perform initial state detection and classification of components according to all units.

S110. Describe the component and all its units through an array dyFlag with dyN items.

Among them, dyFlag _i describes the i-th unit in the component; when the status of the i-th unit in the component is intact, dyFlag _i =1; when the status of the i-th unit in the component is failure, dyFlag _i =0.

S120. Deem the components with the same unit state of dyN items as having the same initial state and classify them into the same component category, and mark different component categories with the category number j.

The matrix bjState is used to describe the classification results of the components according to the initial state of the unit; the j-th row vector bjState(j,:) in bjState is described as [dyFlag nFlag _j ] in the data structure, which is used to describe the initial state check and The results of the classification, namely:

bjState(j,:)=[dyFlag _j1 dyFlag _j2 ... dyFlag _jdyN nFlag _j ], nFlag _j describes the number of jth type components.

In this embodiment, when facing more than one component, the same maintenance scheme is adopted for components of the same category, that is, if a certain unit in the component needs to be replaced, all components in the same category will replace the unit consistently.

S200. Acquire a series of spare parts solutions for preventive maintenance components of the category of the component according to the components in each category of initial state.

S210. Acquire a spare parts solution with category number j=1.

S220. Describe the spare parts scheme of the preventive maintenance component with category number j=1 through the array Ns;

Among them, Ns＝[Ns ₁ Ns ₂ ... Ns _dyN ], Ns _i is the number of spare parts of unit i, 1≤i≤dyN;

The serial number of the initialization scheme j2=0.

S230. For the spare part scheme Ns whose category number j=1, call the component spare part scheme evaluation module to evaluate the maintenance effect of the spare part scheme Ns, and record the evaluation result as bjP.

S240, updating the matrix plan.

j2=j2+1.

cost

The j2th row vector of the configuration matrix, plan(j2,:)=[j2 bjP Mnow Ns].

即Ns_i是否全为1；是则进入S260，否则生成下一备件方案，进入S251。S250, judge

That is, whether Ns _i are all 1s; if yes, go to S260; otherwise, generate the next spare parts plan and go to S251.

S251, Ns _i are not all 1, then take the candidate sequence number k=1.

S252. Generate a candidate solution Nt; wherein,

S253. According to the candidate plan Nt, call the parts and spare parts plan evaluation module to evaluate the maintenance effect of the candidate plan Nt, and record the evaluation result as bjPt, the marginal cost-benefit ratio

S254, k=k+1; if k≤dyN, then return to S252, otherwise go to S255.

S255. After traversing all k less than or equal to dyN, find the maximum value among all p2m _k , 1≤k≤dyN, and record the candidate serial number Im corresponding to the maximum value; then update Ns to make Ns _Im = nFlag _j , return to S230.

S260. Save the plan, and update planN _j and the matrix planAll.

Take planN _j = j2, planN _j actually records the number of preventive maintenance spare parts for the component with state number j;

Place the matrix planAll in the area starting from the first row and the (j-1)(3+dyN)+1 column of the matrix planAll, or as the area intercepted from the matrix planAll, that is: planAll( 1:j2, iy1:iy2) = plan;

Wherein, iy1=(j-1)(3+dyN)+1, iy2=j(3+dyN).

S270. Assuming that there are N1 initial states of different categories, N1 is a positive integer; j=j+1, if j≤N1, go to S220, otherwise go to S300.

S300, using planAll to combine all types of spare parts plans to obtain a series of overall spare parts plans.

S310. Initialize the array planID describing the overall spare parts plan, planID _j = 1, 1≤j≤N1; wherein the jth item planID _j is the component spare part plan serial number of the jth type of initial state component, that is: initialize the overall spare part plan

Initialize various component optimization flags bjF _j =1, 1≤j≤N1.

Initialize the overall scheme number i2=0.

S320. For the overall spare parts plan planID, call the overall spare parts plan evaluation module, mark the obtained intact parts compliance probability as Pm, mark the spare parts cost as Mm, and save the result to the matrix result:

i2=i2+1, the i2th row of result vector result(i2,:)=[Mm Pm planID];

S330. Optimizing and generating the next overall spare parts plan

是则进入S332；否则进入S400。S331. If

If yes, enter S332; otherwise, enter S400.

S332. Get the part category number j=1.

S333. Get IDt=planID, IDt _j =IDt _j +1.

S334. If IDt _j >planN _j , bjF _j =0, p2m _j =0.

If IDt _j ≤ planN _j , then for the overall spare parts plan IDt, the overall spare parts plan evaluation module is called, and the obtained intact parts compliance probability is marked as Pt _j , and the cost of spare parts is marked as Mt _j , and the marginal benefit-cost ratio is taken

S335, j=j+1, if j≤N1; then go to S333, otherwise go to S336.

S336. Find the maximum value among all p2m _j , 1≤j≤N1, and record its serial number as Im;

是则取planID_Im＝planID_Im+1，Pm＝Pt_Im，Mm＝Mt_Im，保存结果到矩阵result：取i2＝i2+1，取result的第i2行向量result(i2,:)＝[Mm Pm planID]后转S331；if

If yes, take planID _Im = planID _Im +1, Pm=Pt _Im , Mm=Mt _Im , save the result to the matrix result: take i2=i2+1, take the i2th row vector of result result(i2,:)=[Mm Pm planID] then transfer to S331;

则进入S400。if

Then enter S400.

S400. Use the result to draw an efficiency curve and display it on the display screen. The X-axis of the efficiency curve is the cost of the overall spare parts solution, and the Y-axis is the compliance probability Pm=P (n>M). Using the curve, the overall spare parts scheme and the corresponding component spare parts scheme that meet the requirement of Pm≥P1 and have the least spare parts cost can be obtained from the query.

The display screen may be used to display information input by or provided to the user and various applications installed in the electronic terminal. The display screen may include a display panel, and optionally, the display panel may be configured in the form of a liquid crystal display (Liquid, Crystal, Display, LCD), an organic light-emitting diode (Organic, Light-Emitting, Diode, OLED), or the like. Further, the touch panel can cover the display panel. When the touch panel detects a touch operation on or near it, it will be sent to the processor to determine the type of the touch event, and then the processor will display the touch event on the display panel according to the type of the touch event. Provide corresponding visual output.

Preferably, the component spare part solution evaluation module mentioned in S230 and S253 is used to evaluate the maintenance effect of the spare part solution Ns and the candidate solution Nt, specifically to evaluate the maintenance effect of the spare part solution Ns, examples are as follows.

式中；

Pick

in the formula;

dyP _i is the reliability of the i-th unit at time T2.

The unit reliability calculation formula for common life distribution types is as follows:

If the lifetime of the i-th unit obeys the Gamma distribution Ga(α,b), then there is

If the lifetime of the i-th unit obeys the lognormal distribution LN(μ,σ ² ), then

If the lifetime of the i-th unit obeys the normal distribution N(μ,σ ² ), then

If the life of the i-th unit obeys the Weibull distribution W(α,b), then there is

且p1＝p0。由于指数分布具有“剩余寿命等于寿命”的特性；那么本实施例的部件备件方案评估模块不对寿命服从指数分布Exp(μ)的第i项单元的预防性维修备件问题进行处理，而是本实施例在遭遇第i项单元的寿命服从指数分布Exp(μ)时，仅直接的更换失效单元或不更换未失效单元。Further, in this embodiment, if the lifetime of the i-th unit obeys the exponential distribution Exp(μ), then there is

And p1=p0. Because the exponential distribution has the characteristic of "the remaining life is equal to the life"; then the parts and spare parts program evaluation module of this embodiment does not process the preventive maintenance spare parts of the i-th unit whose life obeys the exponential distribution Exp(μ), but this implementation For example, when the lifetime of the i-th unit obeys the exponential distribution Exp(μ), only the failed unit is directly replaced or the non-failed unit is not replaced.

Preferably, the overall spare parts program evaluation module mentioned in S320 and S334 in this embodiment is used to obtain the compliance probability and spare parts cost of intact components; the specific steps for obtaining the compliance probability Pm and spare parts cost Mm are as follows.

S321. Initialize the category number j=1, and spare part cost Mt=0.

S322. Take out from planAll the spare parts plan serial number j2=planID _j related information of the spare parts plan of the jth type of initial state parts: part reliability bjP and spare parts cost Mnow;

Take Mt=Mt+Mnow;

1≤k≤1+nFlag_j，得到向量p1；S323. Take x1 as an increasing sequence from 0 to nFlag _j , with a step size of 1, and iteratively calculate the probability corresponding to x1 _k

1≤k≤1+nFlag _j , get the vector p1;

S324. If j>1, calculate the convolution of p1 and p2, and assign the result to p2, that is: p2=p1*p2, "*" in this embodiment is used to indicate convolution calculation;

Otherwise, take p2=p1.

S325. Take j=j+1, if j≤N1 go to S322, otherwise go to S326.

总体备件方案planID的备件费用Mm＝Mt。S326. Obtain the probability of reaching the standard of the overall spare parts plan planID

The spare parts cost Mm=Mt of the overall spare parts plan planID.

Further, the method or device of this embodiment further provides the following simulation model after obtaining the scheme, which is used to simulate the guarantee effect of the spare parts scheme Ns on any component after the preventive maintenance is completed.

S100. Get the unit serial number i=1.

S200. Generate a random number sT, and sT obeys the life distribution law of the unit;

simT_i用于模拟完成预防性维修后该单元的寿命。Pick

simT _i is used to simulate the life of the unit after preventive maintenance has been completed.

S300, i=i+1, if i≤dyN, go to S200, otherwise go to S400.

S400. Find the minimum value from all simT _i , 1≤i≤dyN, and denote it as Tm.

S500. If Tm>T2, the component is in good condition at T2, and record the state value djF=1; otherwise, record djF=0.

After simulating the status values of all components, the cumulative result of all status values is the number n of intact components at time T2 after the preventive maintenance ends.

After a large number of simulations based on the above simulation model, the frequency of n>M is counted, which is the simulation result of the probability of reaching the standard P (n>M).

Please refer to FIG. 2 , this embodiment discloses a device for constructing a preventive maintenance spare parts plan.

The device includes a classification module, a spare parts module, a combination module and a selection module.

Wherein the classification module is used to classify at least one component with the same initial state in the same component category. The spare parts module is used to obtain at least one spare parts solution in each of the parts categories; the combination module combines any one of the spare parts solutions in each of the parts categories into at least one overall spare parts solution, and evaluates the overall spare parts solution Probability of reaching the standard; the selection module selects the overall spare parts plan whose probability of reaching the standard meets the preset condition as the preventive maintenance spare parts plan.

This embodiment uses the following example to illustrate the method for optimizing the preventive maintenance spare parts plan.

A component is composed of 4 units, and the life distribution of each unit is shown in Table 1. There are a total of 18 sets of the same batch of parts in the warehouse. In the first year of T1, the state inspection of all parts is carried out. The inspection results are shown in Table 2, and the preventive maintenance of all parts is completed. After the preventive maintenance is over, the probability that the number n of intact parts is greater than M at year T2 is recorded as P(n>M). It is now known that T1=4, T2=8 and M=14.

Follow the above steps to develop a spare parts plan and draw an efficiency curve. If it is required that P(n>M) is not less than 0.85, find out the spare parts plan that meets the requirements and has the least spare parts expenditure.

Table 1. Life distribution of each unit

Table 2. Result of component status inspection before preventive maintenance

St100, before preventive maintenance of components, the results of initial state inspection and classification of components are shown in Table 2;

St200, for each type of initial state parts, optimize a series of spare parts plans for preventive maintenance of various parts.

For component category 1, after calculation, the results of the component spare parts scheme are shown in Table 3.

Table 3. Component category 1

For component category 2, after calculation, the results of the component spare parts scheme are shown in Table 4.

Table 4. Component category 2

For component category 3, after calculation, the results of the component spare parts scheme are shown in Table 5.

Table 5. Component category 3

St300, using the above results to optimize the combination of spare parts schemes of various components to obtain a series of overall spare parts schemes.

St310. Initialize the overall spare parts plan planID=[1 1 1], initialize various component optimization flags bjF=[11 1], initialize the overall plan number i2=0.

St320. For the overall spare parts plan planID, call the overall spare parts plan evaluation module, obtain the qualified probability of intact parts Pm=0.6177, spare parts cost Mm=607.6, save the result to the matrix result: take result(1,:)=[MmPm planID].

St330, optimize to generate the next overall spare parts plan, after multiple optimizations, the final result is shown in Table 6.

Table 6. Overall Spare Parts Program Results

St400, the calculation is terminated; please refer to Fig. 3, and Fig. 3 also describes the simulation results obtained by simulating the probability of reaching the standard of the overall spare parts scheme through the simulation model provided by this embodiment, and the method or device of this embodiment The evaluation result obtained by evaluating the compliance probability of the overall spare parts scheme. Figure 3 clearly shows that the error of the simulation results and evaluation results of the probability of reaching the standard is small, which meets the engineering requirements. The method or device of this embodiment can realize the optimization of the spare parts solution.

From the above results, it can be seen that the probability of reaching the standard of the No. 6 overall spare parts plan is P(n>M)=0.885, which meets the requirement that P(n>M) is not less than 0.85, and the cost of spare parts is 1429. 3, 2, 3, at this time: for the first category of components, the No. 3 component spare parts plan is adopted, and the number of spare parts of each unit is 7, 0, 7, 0 respectively; for the second category of components, the No. 2 component spare parts plan is adopted, The number of spare parts for each unit is 6, 0, 6, 0 respectively; for the third category of parts, the No. 3 spare parts plan is adopted, and the number of spare parts for each unit is 5, 5, 5, 5 respectively. General Spare Parts Program No. 6 may be used as the final preventive maintenance spare parts program.

The above descriptions are only preferred embodiments of the present invention, and are not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and changes. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention shall be included within the protection scope of the present invention.

Claims (10)

1. A preventive maintenance spare parts program optimization method, characterized in that, The methods include: S100. Classifying at least one component with the same initial state in the same component category; S200. Acquire at least one spare parts solution in each component category; S300. Combining any one of the spare parts schemes in each of the component categories into at least one overall spare parts scheme, and evaluating the compliance probability of the overall spare parts scheme; S400. Select the overall spare parts plan whose compliance probability meets a preset condition as the preventive maintenance spare part plan. 2. the preventive maintenance spare parts scheme optimization method as claimed in claim 1, is characterized in that, The S100 is configured with: S110. Obtain the state of all or part of the units of the component; S120. Determine the initial state according to the sound state or failure state of all the units, and classify the components with the same initial state into the same component category. 3. the preventive maintenance spare parts scheme optimization method as claimed in claim 2, is characterized in that, S110 is configured to describe all or part of the units of the component through an array dyFlag having dyN items, The i-th variable of dyFlag is dyFlag _i , which is used to describe the sound state or the failure state of the unit in the i-th item in the component; S120 is configured to describe the initial state and the component category through the matrix bjState, the j-th row vector bjState(j,:) of bjState is [dyFlag nFlag _j ], and nFlag _j is used to describe the component category of the j-th category The number of parts. 4. the preventive maintenance spare parts scheme optimization method as claimed in claim 3, is characterized in that, S200 is configured as: S210, make the category number j=1; S220, the serial number of the initialization scheme j2=0, A spare part scheme for components in the component category described in category j is described by an array Ns,

Ns _i is the number of spare parts of the unit described in item i 1≤i≤dyN; S230. Invoke a component spare part solution evaluation module, and evaluate the maintenance effect of the spare part solution Ns as bjP; S240, Take j2=j2+1, acquisition fee

Update matrix plan, plan(j2,:)=[j2 bjP Mnow Ns]; S250, judge

exist

When entering S260, in

generate the next spare parts plan; S260. Take planAll(1:j2, iy1:iy2)=plan, iy1=(j-1)(3+dyN)+1, iy2=j(3+dyN), take planN _j =j2; S270. Take j=j+1, return to S220 when j≤N1, and enter S300 when j>N1; S300 is configured to use planAll to combine a spare part plan for the parts in each category of the category parts into at least one overall spare part plan. 5. the preventive maintenance spare parts scheme optimization method as claimed in claim 4, is characterized in that, The S250 is configured to judge

exist

When entering S260, in

When generating the next spare parts plan, enter S251; The S250 is configured with: S251. Take the candidate sequence number k=1, S252. Generate a candidate solution Nt,

S253. According to the candidate solution Nt, call the parts and spare parts solution evaluation module, and evaluate the maintenance effect of the candidate solution Nt as bjPt, and the marginal benefit cost ratio

S254, take k=k+1, return to S252 when k≤dyN, enter S255 when k>dyN, S255. Find the maximum value of p2m _k among all p2m _k , 1≤k≤dyN, record the candidate serial number Im corresponding to the maximum value of p2m _k , update Ns _Im =nFlag _j in Ns, and return to S230. 6. the preventive maintenance spare parts scheme optimization method as claimed in claim 3, is characterized in that, The S300 is configured with: S310. Initialize the array planID describing the overall spare parts plan, planID _j = 1, 1≤j≤N1; the planID _j of the jth item is the serial number of the spare parts plan in the component category of the jth class, Initialize the optimization flag bjF _j =1, 1≤j≤N1 for each type of component category, Initialize the overall scheme number i2=0; S320. Invoking the evaluation module of the overall spare parts scheme, obtaining the probability of reaching the standard of the overall spare parts scheme as Pm, and the cost of the spare parts as Mm, Take i2=i2+1, Vector result(i2,:)=[Mm Pm planID] in row i2 of matrix result; The S330 is configured with: S331. Judgment

exist

Enter S332 when

enter S400, S332. Take the serial number j=1 of the component category, S333. Get IDt = planID, IDt _j = IDt _j + 1, S334. When IDt _j > planN _j , take bjF _j =0 and p2m _j =0, When IDt _j ≤ planN _j , the overall spare parts program evaluation module is called, and the attainment probability of the overall spare parts program is denoted as Pt _j , the cost of spare parts is denoted as Mt _j , and the marginal benefit-cost ratio

S335, set j=j+1, enter S333 when j≤N1, enter S336 when j>N1, S336. Find the serial number corresponding to the maximum value in all p2m _j , 1≤j≤N1, which is Im, exist

, take planID _Im ＝planID _Im + 1, Pm＝Pt _Im , Mm＝Mt _Im , save the result to the matrix result: take i2＝i2+1, take the i2th row vector of result result(i2,:)＝[Mm Pm planID] and return to S331; exist

, turn to S400; S400 is configured to acquire the overall spare parts plan with the attainment probability Pm=P (n>M) as the preventive maintenance spare part plan according to the matrix result. 7. The preventive maintenance spare parts scheme optimization method as claimed in claim 6, is characterized in that, S400 is configured to draw an efficiency curve according to the matrix result; The X-axis of the cost curve is the cost of the overall spare parts plan, and the Y-axis is the probability of reaching the standard Pm=P (n>M); The preventive maintenance spare parts plan is obtained by querying the curve to obtain the overall spare parts plan with the attainment probability Pm=P (n>M) and the least spare part cost. 8. The preventive maintenance spare parts scheme optimization method as claimed in claim 4, is characterized in that, Invoking the component spare parts solution evaluation module to evaluate the maintenance effect of the spare parts solution is configured as: Pick

dyP _i is the reliability of the i-th unit after the preset time, p0 and p1 are the values that the life of the i-th unit obeys the gamma distribution, lognormal distribution, normal distribution or Weibull distribution. 9. The preventive maintenance spare parts scheme optimization method as claimed in claim 6, characterized in that, Invoking the component spare parts scheme evaluation module to evaluate the overall spare parts scheme is configured as: S321. Initialize the category number j=1, spare parts cost Mt=0; S322. Obtain the spare part plan, reliability bjP, and spare part cost Mnow of the spare part plan number j2=planID _j in the component category of the jth category from planAll; Take Mt=Mt+Mnow; S323. Take x1 as an increasing sequence from 0 to nFlag _j , with a step size of 1, and iteratively calculate the probability corresponding to x1 _k

Get the vector p1; S324. Calculate the convolution of p1 and p2 when j>1, and assign the result to p2; Otherwise, take p2=p1; S325. Take j=j+1, return to S322 when j≤N1, and enter S326 when j>N1; S326. Obtain the probability of reaching the standard of the overall spare parts plan planID

The spare parts cost Mm=Mt of the overall spare parts plan planID. 10. A device for constructing a preventive maintenance spare parts plan, characterized in that, The devices include: A classification module, configured to classify at least one component with the same initial state in the same component category; a spare parts module for obtaining at least one spare parts solution in each of said component categories; A combination module, combining any one of the spare parts schemes in each of the component categories into at least one overall spare parts scheme, and evaluating the compliance probability of the overall spare parts scheme; A selection module, selecting the overall spare parts plan whose compliance probability meets a preset condition as the preventive maintenance spare parts plan.

Images