System reliability simulation calculation method based on conduction effect

Abstract

The invention discloses a system reliability simulation calculation method based on a conduction effect. The method comprises the following steps: step 1, constructing structural matter elements of different components in a system; 2, calculating a conduction effect between the structural matter elements constructed in the step 1, and constructing a conduction chain; 3, according to current historical data or industrial standards, counting the failure quantity of different characteristics of the component, and calculating a basic failure rate mu; and step 4, calculating the reliability of the system based on the failure rates [mu] i and [mu] j and the conduction effect rij of different structural matter elements in the conduction chain. According to the system reliability simulation calculation method based on the conduction effect, the reliability operation rule of the operation parameters can be dynamically calculated through setting of the first step, the second step, the third step and the fourth step, and system reliability estimation is facilitated.

Description

Technical Field

[0001] The present invention relates to a calculation method, and more particularly to a system reliability simulation calculation method based on conduction effect. Background Art

[0002] Due to the large number of components and complex relationships of complex systems, the reliability of complex systems will gradually degrade over time, hiding huge safety risks. Traditional reliability calculation methods mainly calculate the basic failure rate of each component, establish the composition relationship between different product structures through reliability block diagrams, and then generate the reliability of the entire system. This method is simply and intuitively applied in many mechanical transmission structures, but it does not take into account the system operation process and the conduction effect between different parameters, that is, when a certain reliability variable changes, it will affect the change of related reliability variables, and it is impossible to measure the change of reliability.

[0003] As the scale of complex systems becomes larger and larger, many researchers have begun to study the risk assessment of different faults in recent years. For example, the invention patent: A reliability analysis method and system for high-end complex equipment (CN115829200A) discloses a reliability analysis method and system for high-end complex equipment. Through the weight vector of risk factors and the risk assessment matrix, the risk-benefit value of each failure mode is calculated, which reduces the human factors in the reliability analysis process, reduces the uncertainty in the reliability analysis process, and effectively improves the objectivity and accuracy of the reliability analysis of complex equipment. The invention patent: A reliability simulation evaluation method and system for a complex hybrid system (CN202111417023.9) decomposes the reliability model of a complex hybrid system into a combination of reliability models of each single system, performs logical operations on the reliability index of each device in each single system virtual unit, obtains the reliability index of each single system virtual unit, calculates the reliability index of all single system virtual units, obtains the reliability index of the reliability model in any simulation task, repeats the simulation task several times, and obtains the reliability index set of the reliability model in all simulation tasks.

[0004] Although the above methods take into account the impact of reliability parameter changes on system reliability, they still treat the reliability system as a static system and do not take into account the conduction effect between different reliability variables. In engineering, many complex systems transmit actions through large-scale servo motors. When a motor fails, it often affects servo motors with a farther transmission relationship through the transmission chain, rather than affecting adjacent motors. In this case, it is impossible to calculate using traditional reliability block diagrams. Summary of the invention

[0005] In view of the deficiencies in the prior art, the object of the present invention is to provide a method for dynamically calculating system reliability by comprehensively considering the conduction relationship between all operating variables within the system.

[0006] To achieve the above object, the present invention provides the following technical solution: a system reliability simulation calculation method based on conduction effect, characterized in that it includes the following steps:

[0007] Step 1: Construct the structural elements of different components within the system;

[0008] Step 2: Calculate the conduction effect between the structural elements constructed in step 1 and construct a conduction chain;

[0009] Step 3: According to the current historical data or industry standards, count the number of failures of different characteristics of components and calculate the basic failure rate μ;

[0010] Step 4: Based on the failure rate μ of different structural elements in the conduction chain i , μ j , conduction effect ij Calculate the reliability of the system.

[0011] As a further improvement of the present invention, the specific steps of constructing the structural elements of different components in step 1 are as follows:

[0012] The structural element M is used to represent the smallest motion unit in the system:

[0013]

[0014] The structural object element includes an ordered triple consisting of a structure name N, a feature name c, and a value v of N with respect to c. N represents the name of the servo motor, c represents the operating characteristics of the servo motor object element, including speed C1, torque C2, vibration C3, temperature C4, and slip C5, v represents the value before the operating characteristic changes, and v' represents the value after the operating characteristic changes.

[0015] As a further improvement of the present invention, the specific steps of calculating the conduction effect in step 2 are as follows:

[0016] Through the structural element M i Features of c i Calculate the corresponding value change Δv i :

[0017] Δv i =v' i -v i

[0018] By triggering the structural element M j Features of c j Calculate the value change Δvj :

[0019] Δv j =v' j -v j

[0020] When Δv i and Δv j When the transformation characteristics of are the same, that is, c i =c j , the conduction effect r ij Defined as:

[0021]

[0022] Among them, the correlation of the conduction effect is as follows:

[0023] When Δv1 and Δv2 have the same sign, it means that the change of v1 is positively correlated with v2, that is, the change of v1 will accelerate the change of v2; when Δv1 and Δv2 have opposite signs, it means that the change of v1 is negatively correlated with v2, that is, the change of v1 will weaken the change of v2; when Δv2=0, it means that the change of v1 is not correlated with v2.

[0024] As a further improvement of the present invention, the specific steps of constructing the conductive chain in step 2 are as follows:

[0025] Calculate the structural element M i and M j The conduction effect r of all quantities under this feature ij , r with the same sign ij The corresponding structural element M i and M j It constitutes a conduction unit, expressed as: r m (v i -v j );

[0026] Determine the direction of the conduction unit r m (v i -v j ): When the conduction effect r ij is positive, forming a forward conduction unit; the conduction effect r ij is negative, forming a reverse conduction unit; the conduction effect r ij It is zero and does not constitute a conduction unit;

[0027] When several conduction units have the same direction, they form a conduction chain, which is expressed as: m (v i -...-v j ).

[0028] As a further improvement of the present invention, the calculation base failure rate μ in step 3 is calculated by the following formula:

[0029]

[0030] Among them, n ci For feature c i The number of failed parts under c For the total parts data.

[0031] As a further improvement of the present invention, the specific steps of calculating the system reliability in step 4 are as follows: the failure rate μ of different structural elements in the conduction chain is i , μ j , conduction effect ij Substituting into the following formula, we get the characteristic c k Reliability changes:

[0032]

[0033] The reliability change of all characteristics c of the system is calculated according to the following formula:

[0034]

[0035] Among them, when ΔR is greater than 0, it means that under this operating condition, the system reliability tends to increase, and the operating state is good; when ΔR is less than 0, it means that the system reliability tends to decrease, and the system has potential reliability risks; when ΔR is equal to 0, it means that the system reliability is in a stable operation stage.

[0036] The beneficial effects of the present invention are as follows: the system reliability simulation calculation method based on the conduction effect of the present invention constructs structural elements of different components inside the system, and uses the structural element M to represent the smallest motion unit inside the system. When the value change Δvi corresponding to the characteristic ci of the structural element Mi triggers the value change Δvj of the characteristic cj of the structural element Mj, the conduction effect between the structural elements is calculated. According to the positive and negative values ​​of the conduction effect, the conduction unit and the conduction chain are generated. According to historical data or industry standards, the number of failures of different characteristics of the components is counted, and the basic failure rate is calculated. The failure rate and conduction effect of different structural elements in the conduction chain are substituted into the system reliability calculation formula to find the change in system reliability. This method can dynamically calculate the reliability operation law of the operating parameters, which is helpful for system reliability estimation. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] Figure 1 It is a schematic diagram of the calculation flow of the system reliability simulation calculation method based on the conduction effect of the present invention. DETAILED DESCRIPTION

[0038] The present invention will be further described below in detail with reference to the embodiments shown in the accompanying drawings.

[0039] Reference Figure 1 As shown, a system reliability simulation calculation method based on conduction effect in this embodiment includes the following steps:

[0040] Step 1: Construct the structural elements of different components within the system;

[0041] Step 2: Calculate the conduction effect between the structural elements constructed in step 1 and construct a conduction chain;

[0042] Step 3: According to the current historical data or industry standards, count the number of failures of different characteristics of components and calculate the basic failure rate μ;

[0043] Step 4: Based on the failure rate μ of different structural elements in the conduction chain i , μ j , conduction effect ij The reliability of the calculation system can be calculated by setting up steps one to four, and then the conduction chain can be constructed by using structural elements, and then the reliability calculation can be effectively realized through the constructed conduction chain.

[0044] As an improved specific implementation method, the specific steps of constructing the structural elements of different components in step 1 are as follows:

[0045] The structural element M is used to represent the smallest motion unit in the system:

[0046]

[0047] The structural object element includes an ordered triple consisting of a structure name N, a feature name c, and N's value v with respect to c. N represents the name of the servo motor, c represents the operating characteristics of the servo motor object element, including speed C1, torque C2, vibration C3, temperature C4, and slip C5. v represents the value before the operating characteristic changes, and v' represents the value after the operating characteristic changes. Through the above, the structural object element can be simply and effectively constructed, for example:

[0048]

[0049] As an improved specific implementation method, the specific steps of calculating the conduction effect in step 2 are as follows: i Features of c i Calculate the corresponding value change Δv i :

[0050] Δv i =v' i -v i

[0051] By triggering the structural element M j Features of c j Calculate the value change Δv j :

[0052] Δv j =v' j -v j

[0053] When Δv i and Δv j When the transformation characteristics of are the same, that is, c i =c j , the conduction effect r ij Defined as:

[0054]

[0055] Among them, the correlation of the conduction effect is as follows:

[0056] When Δv1 and Δv2 have the same sign, it means that the change of v1 is positively correlated with v2, that is, the change of v1 will accelerate the change of v2; when Δv1 and Δv2 have opposite signs, it means that the change of v1 is negatively correlated with v2, that is, the change of v1 will weaken the change of v2; when Δv2=0, it means that the change of v1 is irrelevant to v2. The above steps can effectively calculate the conduction effect. Take gear wear as an example: if the variable wear hardness of the driving wheel increases by 10Hv (Rockwell hardness), causing the wear hardness of the driven wheel to decrease by 25Hv (Rockwell hardness), then Δv1=10Hv, Δv2=-25Hv, and its conduction effect r=-2.5.

[0057] As an improved specific implementation method, the specific steps of constructing the conductive chain in step 2 are as follows: Calculate the structural element M i and M j The conduction effect r of all quantities under this feature ij , r with the same sign ij The corresponding structural element M i and M j It constitutes a conduction unit, expressed as: r m (v i -v j );

[0058] Determine the direction of the conduction unit r m (v i -v j ): When the conduction effect r ij is positive, forming a forward conduction unit; the conduction effect r ij is negative, forming a reverse conduction unit; the conduction effect r ij It is zero and does not constitute a conduction unit;

[0059] When several conduction units have the same direction, they form a conduction chain, which is expressed as: m (v i -...-v j ), through the setting of the above steps, the conduction chain can be effectively constructed.

[0060] As a specific implementation of the improvement, the calculation base failure rate μ in step 3 is calculated by the following formula:

[0061]

[0062] Among them, n ci For feature c i The number of failed parts under c For the total parts data, the basic failure rate μ can be simply and effectively calculated using the above formula.

[0063] As an improved specific implementation method, the specific steps of calculating the system reliability in step 4 are as follows: the failure rate μ of different structural elements in the conduction chain is i , μ j , conduction effect ij Substituting into the following formula, we get the characteristic c k Reliability changes:

[0064]

[0065] The reliability change of all characteristics c of the system is calculated according to the following formula:

[0066]

[0067] Among them, when ΔR is greater than 0, it means that under this operating condition, the system reliability tends to increase, and the operating state is good; when ΔR is less than 0, it means that the system reliability tends to decrease, and the system has potential reliability risks; when ΔR is equal to 0, it means that the system reliability is in a stable operation stage. Through the above setting, the reliability calculation can be simply and effectively realized.

[0068] This embodiment takes a servo motor system as an example to specifically illustrate the method of this embodiment.

[0069] Step 1: Constructing a set of structural elements

[0070] For example, the structural element M of four servo motors at time 1 is:

[0071]

[0072] The structural element M at time 2:

[0073]

[0074]

[0075] Step 2: Calculate the value changes of different features at time 1 and time 2.

[0076] Since the characteristics before and after the operation are the same, the first conductivity calculation method is used to obtain the impact on other characteristic values ​​when the first characteristic value changes.

[0077] N <![CDATA[Δv1]]> <![CDATA[Δv2]]> <![CDATA[Δv3]]> <![CDATA[Δv4]]> <![CDATA[Δv5]]> 1 10 -1.8 1.57 -9 0.02 2 -130 -6.5 -0.48 13.2 -0.03 3 10 22.1 -1.65 -26 -0.06 4 10 10.3 -1.42 -28.9 -0.04

[0078] Taking the r1 conduction chain as an example, r1(v1-v2)=10, which means that the change in the speed of servo motor 1 is in the same direction as the change in the speed of servo motor 2, which means that when the speed of servo motor 1 increases, the speed of servo motor 2 will also increase accordingly.

[0079] According to rule 2.3, the forward conduction loop: r1(v1-v3-v4), r1(v1-v4); the reverse conduction loop: r1(v1-v2-v3-v4) According to formula (1), the conduction effect r1 of the four servo motors in feature 1 is calculated

[0080]

[0081] Step 4: According to the historical data of the quality department, the number of failures of different characteristics of the servo motor is counted, and the basic failure rate is calculated according to formula (2)

[0082] N μ1 μ2 μ3 μ4 μ5 1 0.09 0.08 0.03 0.07 0.07 2 0.07 0.07 0.02 0.06 0.09 3 0.08 0.08 0.01 0.05 0.05 4 0.09 0.06 0.06 0.08 0.07

[0083] Step 5: Reliability calculation

[0084] According to formula (3), the reliability changes of the four servo motors under different characteristics are calculated respectively.

[0085] Description Feature 1: System reliability increased by 1.62% during stable operation

[0086] Step 6: Similarly, calculate the reliability change of the four servo motors in features 2, 3, 4, and 5 under feature 2:

[0087] <![CDATA[r2]]> <![CDATA[Δv1]]> <![CDATA[Δv2]]> <![CDATA[Δv3]]> <![CDATA[Δv4]]> <![CDATA[Δv1]]> 1 3.6111 -12.2778 -5.7222 <![CDATA[Δv2]]> 1 -3.4000 -1.5846 <![CDATA[Δv3]]> 1 0.4661 <![CDATA[Δv4]]> 1

[0088]

[0089] Reliability change calculation under feature 3:

[0090] <![CDATA[r3]]> <![CDATA[Δv1]]> <![CDATA[Δv2]]> <![CDATA[Δv3]]> <![CDATA[Δv4]]> <![CDATA[Δv1]]> 1 -0.3057 -1.0510 -0.9045 <![CDATA[Δv2]]> 1 3.4375 2.9583 <![CDATA[Δv3]]> 1 0.8606 <![CDATA[Δv4]]> 1

[0091]

[0092] Reliability change calculation under feature 4:

[0093] <![CDATA[r4]]> <![CDATA[Δv1]]> <![CDATA[Δv2]]> <![CDATA[Δv3]]> <![CDATA[Δv4]]> <![CDATA[Δv1]]> 1 -1.4667 2.8889 3.2111 <![CDATA[Δv2]]> 1 -1.9697 -2.1894 <![CDATA[Δv3]]> 1 1.1115 <![CDATA[Δv4]]> 1

[0094]

[0095] Reliability change calculation under feature 5:

[0096] <![CDATA[r5]]> <![CDATA[Δv1]]> <![CDATA[Δv2]]> <![CDATA[Δv3]]> <![CDATA[Δv4]]> <![CDATA[Δv1]]> 1 -1.5000 -3.0000 -2.0000 <![CDATA[Δv2]]> 1 2.0000 1.3333 <![CDATA[Δv3]]> 1 0.6667 <![CDATA[Δv4]]> 1

[0097]

[0098] Step 6: Calculate the reliability change under all features according to formula (4)

[0099] In summary, after considering the conduction effect between different servo motors, the system reliability will increase by 0.16%, indicating that the system reliability has an increasing trend and the system is in good operating condition.

[0100] The above is only a preferred embodiment of the present invention, and the protection scope of the present invention is not limited to the above embodiments. All technical solutions under the concept of the present invention belong to the protection scope of the present invention. It should be pointed out that for ordinary technicians in this technical field, some improvements and modifications without departing from the principle of the present invention should also be regarded as the protection scope of the present invention.

Claims

1. A system reliability simulation calculation method based on conduction effect, characterized in that: The steps include: Step 1: Construct the structural elements of different components within the system; Step 2: Calculate the conduction effect between the structural elements constructed in step 1 and construct a conduction chain; Step 3: According to the current historical data or industry standards, count the number of failures of different characteristics of components and calculate the basic failure rate μ; Step 4: Based on the failure rate μ of different structural elements in the conduction chain i , μ j , conduction effect ij Calculate the reliability of the system.

2. The system reliability simulation calculation method based on conduction effect according to claim 1 is characterized in that: The specific steps of constructing the structural elements of different components in step 1 are as follows: The structural element M is used to represent the smallest motion unit in the system: The structural object element includes an ordered triple consisting of a structure name N, a feature name c, and a value v of N with respect to c. N represents the name of the servo motor, c represents the operating characteristics of the servo motor object element, including speed C1, torque C2, vibration C3, temperature C4, and slip C5, v represents the value before the operating characteristic changes, and v' represents the value after the operating characteristic changes.

3. The system reliability simulation calculation method based on conduction effect according to claim 2 is characterized in that: The specific steps for calculating the conduction effect in step 2 are as follows: Through the structural element M i Features of c i Calculate the corresponding value change Δv i : Δv i =v' i -v i By triggering the structural element M j Features of c j Calculate the value change Δv j : Δv j =v' j -v j When Δv i and Δv j When the transformation characteristics of are the same, that is, c i =c j , the conduction effect r ij Defined as: Among them, the correlation of the conduction effect is as follows: When Δv1 and Δv2 have the same sign, it means that the change of v1 is positively correlated with v2, that is, the change of v1 will accelerate the change of v2; when Δv1 and Δv2 have opposite signs, it means that the change of v1 is negatively correlated with v2, that is, the change of v1 will weaken the change of v2; when Δv2=0, it means that the change of v1 is not correlated with v2.

4. The system reliability simulation calculation method based on conduction effect according to claim 3 is characterized in that: The specific steps of constructing the conduction chain in step 2 are as follows: Calculate the structural element M i and M j The conduction effect r of all quantities under this feature ij , r with the same sign ij The corresponding structural element M i and M j It constitutes a conduction unit, expressed as: r m (v i -v j ); Determine the direction of the conduction unit r m (v i -v j ): When the conduction effect r ij is positive, forming a forward conduction unit; the conduction effect r ij is negative, forming a reverse conduction unit; the conduction effect r ij It is zero and does not constitute a conduction unit; When several conduction units have the same direction, they form a conduction chain, which is expressed as: m (v i -...-v j ).

5. The system reliability simulation calculation method based on conduction effect according to claim 4 is characterized in that: The calculation base failure rate μ in step 3 is calculated by the following formula: Among them, n ci For feature c i The number of failed parts under c For the total parts data.

6. The system reliability simulation calculation method based on conduction effect according to claim 5 is characterized in that: The specific steps of calculating the system reliability in step 4 are as follows: The failure rate μ of different structural elements in the conduction chain i , μ j , conduction effect ij Substituting into the following formula, we get the characteristic c k Reliability changes: The reliability change of all characteristics c of the system is calculated according to the following formula: Among them, when ΔR is greater than 0, it means that under this operating condition, the system reliability tends to increase, and the operating state is good; when ΔR is less than 0, it means that the system reliability tends to decrease, and the system has potential reliability risks; when ΔR is equal to 0, it means that the system reliability is in a stable operation stage.

Images