Reliability block diagram of series-parallel system based on discrete event simulation

Abstract

The present invention belongs to the technical field related to reliability assessment, and discloses a reliability block diagram and calculation method of a series-parallel system based on discrete event simulation. The reliability block diagram includes multiple reliability components and multiple synchronization components. The reliability components correspond one-to-one to the components in the series-parallel system to be processed, and the connection relationship between the reliability components is connected according to the connection relationship between the components; the starting point and the end point of the model are provided with a starting node and an end point node, and the reliability components are directly connected through connecting lines to form a series relationship, and the front and rear ends of the modules formed by the parallel reliability components are respectively inserted with a node and the reliability components are connected to the nodes through connecting lines, and one node corresponds to one synchronization component. The present invention also discloses a calculation method for the reliability of the above reliability block diagram. Through the present invention, the problem that the reliability model and reliability value calculation of the series-parallel system cannot be established through discrete event simulation is solved.

Description

Technical Field

[0001] The present invention belongs to the technical field related to reliability assessment, and more specifically, relates to a reliability block diagram and a calculation method of a series-parallel system based on discrete event simulation. Background Art

[0002] Reliability refers to the ability of a product to complete specified functions under specified conditions and within a specified time. The product here can be any component, device, equipment, system, etc. A series-parallel system is a system composed of multiple components or subsystems connected in series and parallel. In a series-parallel system, if a component or subsystem fails, the reliability of the entire system will be affected.

[0003] The reliability block diagram is a logic diagram that studies the relationship between the system and components from the perspective of reliability. It is a common tool for modeling and calculating system reliability. The reliability value of the system is calculated based on the reliability block diagram. The main methods include Monte Carlo simulation method, fault tree, etc. The Monte Carlo simulation method simulates and randomly samples the system multiple times, obtains the system state with the simulation results, and estimates the system reliability approximately by the frequency of occurrence of the system state. The Monte Carlo simulation method cannot obtain an accurate solution to the system reliability, and in complex systems or when the reliability of components is high, the number of simulations required to obtain the system reliability solution increases sharply, and the calculation time is greatly increased. The fault tree analysis method analyzes the reliability problem of the system based on graphical deduction and Boolean operations. The traditional static fault tree method cannot describe the random process of fault occurrence, and the dynamic fault tree method also has difficulty in solving when dealing with complex systems.

[0004] The process of discrete event system simulation is as follows: at a discrete time point, a simulated event occurs, changes the state of the system, triggers the activity of the object, and the activity ends after a random period of time. After the end, the simulation generates a new event, changes the state of the system again, and triggers a new activity, and repeats until the simulation ends. The key to discrete event system simulation lies in the simulation of events, the simulation of the start and end of object activities, and the simulation of object state transitions. Establishing a reliability block diagram model in a discrete event simulation model and calculating system reliability through discrete event simulation is an innovative idea, and there is no relevant literature or patent. Therefore, there is an urgent need for a method and system that can use discrete event simulation to calculate system reliability. Summary of the invention

[0005] In view of the above defects or improvement needs of the prior art, the present invention provides a reliability block diagram and calculation method of a series-parallel system based on discrete event simulation, which solves the problem that the reliability model and reliability value calculation of the series-parallel system cannot be established through discrete event simulation.

[0006] To achieve the above-mentioned purpose, according to one aspect of the present invention, a reliability block diagram of a series-parallel system based on discrete event simulation is provided, wherein the reliability block diagram includes multiple reliability components and multiple synchronization components, the reliability components correspond one-to-one to the components in the series-parallel system to be processed, and the connection relationship between the reliability components is connected according to the series or parallel connection relationship between the components in the series-parallel system to be processed; the starting point and the end point of the reliability block diagram are provided with a starting node and an end node, the reliability components in series are directly connected by connecting lines, and a node is inserted into the front end and the rear end of the module formed by the reliability components in parallel respectively, and the reliability components are connected to the nodes by connecting lines, and one node corresponds to one synchronization component.

[0007] Further preferably, the number of predecessors and the number of successors of the reliability component are both 1, and the number of predecessors and the number of successors of the synchronization component object are both non-negative integers.

[0008] According to another aspect of the present invention, there is provided a method for calculating the reliability of a series-parallel system reliability block diagram based on discrete event simulation, the method comprising the following steps:

[0009] S1 sets the first sampling time, total sampling times, sampling interval, and the reliability calculation function of each reliability component with respect to time. n =0;

[0010] S2 calculates the sampling time corresponding to the current sampling number, calculates the reliability of each reliability component at the current sampling time, assigns a token object to the starting node and initializes the reliability value list of the token object to [1.0], and the token object is passed backward from the starting node;

[0011] S3 updates the token object in accordance with the preset updating methods of the first module and the second module in sequence;

[0012] S4: Determine whether there is only one token object in the entire reliability block diagram and the token object is at the terminal node. If not, return to step S3; if yes, the first element of the reliability list corresponding to the current token object is the reliability of the series-parallel system to be processed under the current sampling number, and delete the token object at the terminal node;

[0013] S5 Order n = n+ 1. Determine the current sampling times n Is it greater than the maximum sampling number? If not, return to step S2; if yes, output the reliability of the series-parallel system to be processed corresponding to each sampling number and draw a curve of the reliability changing with the sampling number.

[0014] Further preferably, the first module updating method in S31 includes the following steps:

[0015] S311 i =1;

[0016] S312 determines whether a token object exists on the current reliability component;

[0017] S313: When there is a token object on the current reliability component, the reliability list is updated according to the preset first relational expression, and the token object on the current reliability component is moved backward, and the process goes to step S314;

[0018] When there is no token object on the current reliability component, go to step S314;

[0019] S314 Order i = i +1, return to step S312 until i Greater than the total number of reliability components.

[0020] Further preferably, in step S312, the first relationship is performed according to the following relationship:

[0021] LP 1 = LP 1 * R it

[0022] in, LP 1 Represents the reliability list of token objects on the current reliability component LP The first element of R it Indicates the current reliability component at sampling time t reliability value.

[0023] Further preferably, the second module updating method comprises the following steps:

[0024] S321 Order j =1;

[0025] S322 determines whether the number of tokens on the current synchronization component is equal to the number of predecessors;

[0026] If yes, update the reliability list and the token object on the transfer synchronization component in the manner of the pre-built third module, and proceed to step S323;

[0027] If not, go directly to step S323;

[0028] Order S323 j =j +1, determine the number of traversed synchronization components j Is it greater than the total number of synchronization components? If so, go to step S4; if not, return to step S322.

[0029] Further preferably, the third module is implemented in the following manner:

[0030] S3221 Calculate the preset value according to the preset second relational expression R t * , determine the reliability list in the first token object of the current synchronization component LC Is the number of elements in less than 2?

[0031] If yes, set the default value R t * Assigned to the first element of the reliability list LC 1 ;

[0032] If not, assign a value to the second element of the reliability list according to the preset third relational expression, and delete the reliability list LC The first element in ;

[0033] Get an updated reliability list LC , delete all token objects in the current synchronization component;

[0034] S3222 Calculates the number of successors of the current synchronization component and selects one of the following methods based on the number of successors:

[0035] When the number of successors is 0, a new token object is created on the current synchronization component, and the reliability list corresponding to the token object is the current reliability list LC ;

[0036] When the number of successors is 1, a new token object is created on the current synchronization component, and the reliability list corresponding to the token object is the current reliability list LC , and drives the token to move to the next position;

[0037] When the number of successors is greater than 1, in the current reliability list LC Insert a new element 1.0 before the first element, thereby updating the reliability list LC At the same time, multiple tokens are generated according to the number of successors of the current synchronization component and issued to each successor object respectively. The reliability list corresponding to each token is copied to the current reliability list. LC .

[0038] Further preferably, in step S3221, the preset second relationship is calculated according to the following relationship:

[0039] R t * =(1-∏ P p=1 (1 - L p,1 ))

[0040] in, P is the number of token objects contained in the current synchronization component, L p,1 It is p Reliability list of synchronization components L p The first data of .

[0041] Further preferably, in step S3221, the preset third relationship is performed according to the following relationship:

[0042] LC 2 = LC 2 * R t *

[0043] in, LC 2 Is a reliability list LC The second element in .

[0044] Further preferably, in step S2, the reliability of each reliability component at the current sampling time is calculated according to the following relationship:

[0045] ;

[0046] in, i is the index of the reliability component, t represents the sampling time, R it Reliability Component i At sampling time t The reliability value when Is a reliability component i Reliability over time t The calculation function of .

[0047] In general, compared with the prior art, the above technical scheme conceived by the present invention provides a reliability block diagram and calculation method of a series-parallel system based on discrete event simulation. The above method, device, computer equipment and storage medium can be used to establish a reliability block diagram of a series-parallel system based on a discrete event simulation model, and the discrete event simulation can quickly calculate the reliability. The reliability value is the reliability value corresponding to the system reliability block diagram. The whole process is simple and can be automatically implemented by a computer, avoiding the tedious steps of using traditional mathematical formulas for calculation, avoiding the shortcomings of low efficiency of manual calculation and unsuitability for complex systems, and overcoming the problem that the Monte Carlo simulation method cannot accurately calculate the system reliability value, providing a new idea for analyzing the reliability of complex systems and improving efficiency. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] Figure 1 It is a flow chart of a reliability calculation method of a series-parallel system reliability block diagram based on discrete event simulation constructed according to a preferred embodiment of the present invention.

[0049] Figure 2 It is a flow chart of a first module updating method constructed according to a preferred embodiment of the present invention.

[0050] Figure 3 It is a flow chart of the second module updating method constructed according to the preferred embodiment of the present invention.

[0051] Figure 4 It is a flow chart of the third module updating method constructed according to the preferred embodiment of the present invention.

[0052] Figure 5 It is a schematic diagram of reliability samplers, reliability components, synchronization components, directed connecting lines, series-parallel relationships and complex parallel relationships in a preferred discrete event simulation model of the present invention.

[0053] Figure 6 It is a schematic diagram of each component in the reliability block diagram of the series-parallel system constructed according to the preferred embodiment of the present invention.

[0054] Figure 7 It is a schematic diagram of various components in a series-parallel system reliability block diagram based on discrete event simulation constructed according to a preferred embodiment of the present invention.

[0055] Figure 8 It is a curve diagram showing the reliability varying with the sampling times constructed according to the preferred embodiment of the present invention. DETAILED DESCRIPTION

[0056] In order to make the purpose, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention. In addition, the technical features involved in the various embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

[0057] A reliability block diagram of a serial-parallel system based on discrete events. The modeling objects of the reliability block diagram include: a reliability sampler, a reliability component, a synchronization component, a directed connection line and a token.

[0058] like Figure 5 As shown, the components in the series-parallel system to be processed are set as reliability components, different components correspond to different reliability components, and the reliability components are connected through directed connection line objects according to the series or parallel connection relationship between the components, so that the series-parallel system to be processed is simplified into an initial reliability block diagram; the starting node and the end node are inserted at the starting end and the end end of the initial reliability block diagram, and a node is inserted into the front end and the back end of the parallel module formed by the parallel components respectively, so as to form a final reliability block diagram, wherein each node is defined as a synchronization component.

[0059] Reliability components are permanent objects in discrete event simulation, representing units or functions with reliability values, corresponding to the "box" elements in the reliability block diagram.

[0060] The synchronization component is a permanent object in discrete event simulation, which represents the starting point or end point of the system function or parallel relationship, corresponding to the "node" element in the reliability block diagram.

[0061] Directed connecting lines are permanent objects in discrete event simulation, corresponding to the "connection" elements in the reliability block diagram, reflecting the direction of system functions or processes.

[0062] For components JP (Reliability Sampler, Reliability Component, Synchronization Component) and Components JS They are connected by directed lines, with the direction being JP point to JS ,say JP for JS The forerunner of JS for JP The successor of the reliability sampler is specified as 1 and the predecessor as 0; the predecessor and successor numbers of any reliability component are both 1, and the predecessor and successor numbers of any synchronization component object are non-negative integers. The connection line between two reliability components represents a reliability series relationship, and the structure between two synchronization components represents a reliability parallel relationship.

[0063] The reliability sampler is a permanent object in discrete event simulation, which can generate token objects at preset time intervals and set the maximum number of sampling times. A complete reliability block diagram of a serial-parallel system based on discrete events contains only one reliability sampler.

[0064] Tokens are temporary objects in discrete event simulation and can be generated by reliability samplers. Tokens can flow in sequence between reliability components and synchronization components along directed connections, and perform creation and deletion operations. The movement events of token objects are driven by the discrete event simulation engine, and the time to pass through each object is 0. The token object contains a reliability value list, which records the calculation process data of the system reliability.

[0065] like Figures 1 to 4 As shown, a reliability block diagram modeling method of a series-parallel system based on discrete event simulation, the modeling method steps include:

[0066] Step 1: Based on the schematic diagram (or functional diagram) of the series-parallel system, establish the system reliability block diagram.

[0067] Step 2: Based on the reliability block diagram, establish a reliability model based on discrete event simulation.

[0068] Step 2.1: For the "boxes" in the reliability block diagram, create corresponding reliability components in the discrete event simulation model , set it to change with the simulation time t Changing reliability function ,in, M is the total number of reliability components.

[0069] Step 2.2: For the "nodes" in the reliability diagram (including starting points, intermediate nodes, and end points), create corresponding synchronization components in the discrete event simulation model. ,in, N is the total number of synchronized components.

[0070] Step 2.3: Create a directed connection line object for the "connection" and logical relationship in the reliability block diagram, and connect the reliability components and synchronization components in the corresponding simulation model.

[0071] Step 3: Create a reliability sampler and connect the reliability sampler to the synchronization component corresponding to the system starting point through a directed connection line. Set the first sampling time of the reliability sampler to t 1, the sampling interval is Δt , the maximum number of sampling times is max .

[0072] A method for calculating a reliability block diagram of a series-parallel system based on discrete event simulation, the steps of which include:

[0073] Step 1: Order n = 0 , running the discrete event based series-parallel system reliability block diagram simulation model.

[0074] Step 2: According to the sampling interval preset by the reliability sampler Δt , the simulation engine advances the simulation clock to t = t 1 + n * Δt time. According to the reliability component RB i The reliability function is used to calculate the reliability The reliability sampler generates a token object and initializes the reliability value list to [1.0]; the discrete event simulation engine drives the token object to move along the directed connection line to the successor object.

[0075] Step 3: Go through all reliability components RB i :if RB i If the token object does not exist on the RB i Token object exists on TP , TP The reliability values ​​are listed as LP ,calculate LP First element LP 1= LP 1* R it , the discrete event simulation engine moves the token object to the successor object according to the directed connection line.

[0076] Step 4: Set up synchronization components SB j The number of predecessors is Pred j , the number of successors is Succ j ; Remember to synchronize components when executing the current step SB j The number of token objects on P , traverse the synchronous components: for synchronous components SB j ,if P = 0 ,or Pred j >0 and Pred j > P , then skip the current synchronization component, otherwise, SBj Follow these steps:

[0077] Step 4.1: Calculate the reliability value and update the reliability value list of the token. Specific steps:

[0078] 1) Calculation R t * = (1-∏ P p=1 (1 - L p,1 )),in, L p,1 Token object List of reliability values L p The first element of .

[0079] 2) Copy the token T 1 List of reliability values L 1 , denoted as LC , Representation list LC No. k elements, K For List LC If K ≥ 2, let LC 2 = LC 2 * R t * , and delete the first element of the list; if K = 1, then let LC 1 = R t * .

[0080] 3) Retention LC , delete all token objects on the current synchronization component.

[0081] Step 4.2: Create a new token object and move it to the successor object respectively. Specific steps:

[0082] 1) If the number of successors of the current synchronization component Succ j =0, a new token object is created and the reliability list is set as LC .

[0083] 2) If the number of successors of the current synchronization component Succ j=1, a new token object is created and the reliability list is set as LC ; The discrete event simulation engine drives the token object to the successor object.

[0084] 3) If the number of successors of the current synchronization component Succ j >1, then LC Insert data before the first element: 1.0, then create a new Succ j Token Object T s (1 ≤s≤Succ j ),set up T s Reliability list L s = LC ;Discrete event simulation engine drives token object T s Move to the corresponding s A successor object.

[0085] Step 5: When executing the current step, if there is one and only one token object in the simulation model and it exists on the synchronization component corresponding to the system endpoint (the successor number is 0), execute step 6; otherwise, execute step 3.

[0086] Step 6: Order n = n+ 1. The first element of the reliability list of the token object is R t ,but t The instantaneous reliability of the series-parallel system at time R t , delete the token object; if n ≤ max , then return to step 2; otherwise, go to step 7.

[0087] Step 7: According to R t The sampling calculation data is summarized to obtain Dynamic reliability change curve of series-parallel system within the time period, ends.

[0088] The present invention will be further described below with reference to specific embodiments.

[0089] like Figure 4 As shown, the reliability block diagram of a ship power system and a steering system is used as an example to further illustrate the modeling method of the series-parallel system reliability block diagram based on discrete event simulation of the present invention.

[0090] Step 1: Based on the schematic diagram (or functional diagram) of the series-parallel system, establish the system reliability block diagram.

[0091] Specifically, Figure 6 The figure shows an example of a reliability block diagram of a ship's power system and steering system.

[0092] Step 2: Based on the reliability block diagram, establish a reliability model based on discrete event simulation.

[0093] Step 2.1: For the "boxes" in the reliability block diagram, create corresponding reliability components in the discrete event simulation model , set it to change with the simulation time t Changing reliability function .

[0094] Preferably, this example is implemented using discrete event simulation software FactorySimulation.

[0095] Specifically, Figure 6 As shown, the "box" represents the module or function of the system, and the "box" includes: 1 engine system, 1 control system, 2 transmission devices, 2 auxiliary systems, 3 propellers, 1 steering hydraulic system pump station, 1 whole ship hydraulic system pump station, 2 hydraulic air systems, 2 hydraulic presses, 1 emergency conversion valve, 1 stop valve group, 1 electro-hydraulic controller, 1 emergency manual valve group, 1 isolation hydraulic valve group, and 1 driver.

[0096] Specifically, a corresponding reliability component is newly created in the discrete event simulation model, such as Figure 7As shown, the “engine system box” corresponds to the reliability component “engine system”, the “control system box” corresponds to the reliability component “control system”, the two “transmission boxes” correspond to the reliability components “transmission 1” and “transmission 2” respectively, the two “auxiliary system boxes” correspond to the reliability components “auxiliary system 1” and “auxiliary system 2” respectively, the three propellers correspond to the reliability components “propeller 1”, “propeller 2” and “propeller 3” respectively, the “steering hydraulic system pump station box” corresponds to the reliability component “steering hydraulic system pump station”, the “whole ship hydraulic system pump station box” corresponds to the reliability component “whole ship hydraulic system pump station”, the two “hydraulic air system boxes” correspond to the reliability components “hydraulic air system 1” and “hydraulic air system 2”, the two “hydraulic press boxes” correspond to the reliability components “hydraulic press 1” and “hydraulic press 2”, the “emergency conversion valve box” corresponds to the reliability component “emergency conversion valve”, the “stop valve group” box corresponds to the reliability component “stop valve group”, the “electro-hydraulic controller box” corresponds to the reliability component “electro-hydraulic controller”, The “Emergency manual valve group box” corresponds to the reliability component “Emergency manual valve group”, the “Isolation hydraulic valve group box” corresponds to the reliability component “Isolation hydraulic valve group”, and the “Driver box” corresponds to the reliability component “Driver”.

[0097] Specifically, this example includes a power system and a steering system. The reliability function of the module in the power system is: hour, ,when hour, . Reliability function of the module in the steering system: hour, ,when hour, ;in, The unit is hours.

[0098] Specifically, the reliability functions of the reliability components “engine system”, “control system”, “transmission 1”, “transmission 2”, “auxiliary system 1”, “auxiliary system 2”, “propeller 1”, “propeller 2” and “propeller 3” are set as f 1 ( t ), set the reliability functions of reliability components “steering hydraulic system pump station”, “whole ship hydraulic system pump station”, “hydraulic air system 1”, “hydraulic air system 2”, “hydraulic press 1”, “hydraulic press 2”, “emergency conversion valve”, “stop valve group”, “electro-hydraulic controller”, “emergency manual valve group”, “isolation hydraulic valve group” and “driver” as follows: f 2 ( t ) .

[0099] Step 2.2: For the "nodes" in the reliability diagram (including starting points, intermediate nodes, and end points), create corresponding synchronization components in the discrete event simulation model. .

[0100] Specifically, for the "nodes" in the reliability block diagram, corresponding synchronization components are newly created in the discrete event simulation model, such as Figure 6 and Figure 7 As shown, “starting point node”, “end point node”, “node 1”, “node 2”, “node 3”, “node 4” “node 5”, “node 6”, “node 7”, “node 8” and “node 9” respectively correspond to the synchronization components “starting point”, “end point”, “node 1”, “node 2”, “node 3”, “node 4” “node 5”, “node 6”, “node 7”, “node 8” and “node 9” in the simulation model.

[0101] Step 2.3: Create a directed connection line object for the "connection" and logical relationship in the reliability block diagram, and connect the reliability components and synchronization components in the corresponding simulation model.

[0102] Specifically, Figure 6 and Figure 7 As shown in the figure, the "engine system box" and the "control system box" are in a series relationship, so the reliability component "engine system" is connected to the "control system" through a directed connection line; the two "transmissions" are in a parallel relationship, so the synchronization component "node 1" is connected to the reliability components "transmission 1" and "transmission 2" through a directed connection line, and then the reliability components "transmission 1" and "transmission 2" are connected to the synchronization component "node 2"; the parallel subsystem composed of the "control system box" and the two "transmissions" is in a series relationship, so the reliability component "control system" is connected to the synchronization component "node 1" through a directed connection line. By analogy, other series and parallel relationships in the system are connected using directed connections.

[0103] Step 3: Create a reliability sampler and connect the reliability sampler to the synchronization component corresponding to the system starting point through a directed connection line. Set the first sampling time of the reliability sampler to t 1, the sampling interval is Δt , the maximum number of sampling times is max .

[0104] Specifically, Figure 7 As shown in the figure, create a new reliability sampler "Reliability Sampler", connect the "Reliability Sampler" with the synchronization component "Starting Point" corresponding to the system starting point through a directed connection line, and set the first sampling time of the sampler to t 1 =0, the sampling interval is Δt =1 (unit: hour), the maximum number of sampling times is max =1000.

[0105] Taking the reliability block diagram of a ship power system and a steering system as an example, the reliability block diagram calculation method of the series-parallel system based on discrete event simulation of the present invention is further explained:

[0106] Execute step 1 and set n =0, run the discrete event based serial-parallel system reliability block diagram simulation model.

[0107] Execute step 2 to further simulate the simulation process. n = 10, according to the sampling interval Δt =1 (unit: hour), first sampling time t 1=0, the simulation engine advances the simulation clock to t =0+10*1=10 moments.

[0108] The reliability is calculated based on the reliability function of the reliability component. f 1 (10) = 1-0.001*10=0.99, f 2 (10) =1-0.000001*10 2 =0.99, that is, the reliability components "engine system", "control system", "transmission 1", "transmission 2", "auxiliary system 1", "auxiliary system 2", "propeller 1", "propeller 2" and "propeller 3" are t =10, the instantaneous reliability function is 0.99, that is, the reliability components "steering hydraulic system pump station", "whole ship hydraulic system pump station", "hydraulic air system 1", "hydraulic air system 2", "hydraulic press 1", "hydraulic press 2", "emergency conversion valve", "stop valve group", "electro-hydraulic controller", "emergency manual valve group", "isolation hydraulic valve group", "driver" are all in the state of reliability. t =10, the instantaneous reliability function is 0.99.

[0109] The discrete event simulation engine drives the movement of the token object, which will move along the directed connection line to the successor object: the synchronization component "starting point".

[0110] Execute step 3 and traverse all reliability components. If no token object exists, skip it.

[0111] Execute step 4 and traverse all synchronization components. In the current model, there is only one token object at the "starting point", which is recorded as the token object. T 1 , the number of predecessors of the synchronization component "starting point" is 1, which meets the execution conditions of steps 4.1 and 4.2.

[0112] Execute step 4.1 for the synchronization component "starting point", calculate the reliability value, and update the reliability value list of the token. Specific steps:

[0113] 1) Calculation R t * =(1-∏ P p=1 (1 - L p,1 ))=1.0;

[0114] 2) Copy T 1 List of reliability values LC = L 1 =[1.0], let LC 1 = R t * = 1.0, we get LC =[1.0];

[0115] 3) Retention LC , delete all token objects at the "starting point".

[0116] Execute step 4.2 for the synchronization component "starting point": create a new token object and move it to the successor object respectively. Specific steps:

[0117] The number of successors is equal to 1, so a new token object is created and the reliability list is set as LC =[1.0]; The discrete event simulation engine drives the token object to the successor object "Turbine System".

[0118] Execute step 5 and return to step 3. Traverse all reliability components. In the current model, there is only one token object on the "turbine system", which is recorded as the token object TP , TP Reliability list LP =[1.0], reliability value of “engine system” R it = 0.99, calculated LP 1 = LP 1 * R it =1.0*0.99= 0.99, list updated LP =[0.99]; The discrete event simulation engine moves the token object to the successor object "control system" according to the directed connection line.

[0119] Execute step 4, traverse all synchronization components, and skip.

[0120] Execute step 5 and return to step 3. Traverse all reliability components. There is 1 token object on the "control system" in the current model, which is recorded as a token pair TP , TP Reliability list LP =[0.99], "Control System" R it = 0.99, calculated LP 1 = LP 1 * R it =0.99*0.99=0.9801, list updated LP =[0.9801]; the discrete event simulation engine moves the token object to the successor object "Node 1" according to the directed connection line.

[0121] Execute step 4 and traverse the synchronization components. In the current model, there is only one token object on "node 1", which is recorded as the token object. T 1. The number of predecessors of the synchronization component "Node 1" is 1, which meets the execution conditions of steps 4.1 and 4.2.

[0122] Execute step 4.1 for the synchronization component "Node 1", calculate the reliability value, and update the reliability value list of the token. Specific steps:

[0123] 1) Calculation R t * =(1-∏ P p=1 (1 - L p,1 ))=0.9801;

[0124] 2) Copy T Reliability value list for 1 LC = L 1=[0.9801], let LC 1= R t * =0.9801, we get LC =[0.9801];

[0125] 3) Retention LC , delete the "node1" token object.

[0126] Execute step 4.2 for the synchronization component "Node 1": Create a new token object and move it to the successor object respectively. Specific steps:

[0127] The number of successors is 2, which is greater than 1, so in LC Insert data before the first element: 1.0, LC =[1.0,0.9801], then create 2 new token objects T 1 ,T 2. Settings T 1 ,T 2 reliability list L 1= L 2= LC =[1.0,0.9801]. Discrete event simulation engine driver token object T 1 ,T 2 are moved to the corresponding successor objects "Transmission 1" and "Transmission 2" respectively.

[0128] Execute step 5 and return to step 3. Traverse all reliability components. There are two token objects in the current model, one in "Transmission 1" and one in "Transmission 2", respectively recorded as T 1 ,T 2. The reliability list is L 1= L 2=[1.0,0.9801].

[0129] Reliability value of "Transmission 1" R it = 0.99, calculated L 1,1 = L 1,1 * R it =1.0*0.99= 0.99, list updated L 1=[0.99,0.9801]; the discrete event simulation engine moves the token object according to the directed connection line T 1 to the successor object "Node 2".

[0130] Reliability value of "Transmission 2" R it = 0.99, calculated L 2,1 = L 2,1 * R it =1.0*0.99= 0.99, list updated L 2=[0.99,0.9801]; The discrete event simulation engine moves the token object according to the directed connection line T 2 to the successor object "Node 2".

[0131] Execute step 4 and traverse the synchronization components. In the current model, there are only two token objects on "Node 2", which are recorded as token objects. T 1 ,T 2. The number of predecessors of the synchronization component "Node 1" is 2, which meets the execution conditions of steps 4.1 and 4.2.

[0132] Execute step 4.1 for the synchronization component "Node 2": Calculate the reliability value and update the reliability value list of the token. Specific steps:

[0133] 1) Calculation R t * =(1-∏ P p=1 (1 - L p,1 ))=0.9999;

[0134] 2) Copy T 1 List of reliability values LC = L 1=[0.99,0.9801], at this time K =2, so let LC 2= LC 2* R t * =0.9801*0.9999=0.98000199, delete the first element of the list, now LC =[0.98000199];

[0135] 3) Retention LC , delete all token objects on "Node 1".

[0136] Execute step 4.2 for the synchronization component "Node 2": Create a new token object and move it to the successor object respectively. Specific steps:

[0137] The number of successors is 2, so in LC Insert data before the first element: 1.0, LC =[1.0,0.98000199]. Then create 2 new token objects T 1 ,T 2. Settings T 1 ,T 2 reliability list L 1= L 2= LC =[1.0,0.98000199]. Discrete event simulation engine driver token object T 1 ,T2 are moved to the corresponding successor objects "Auxiliary System 1" and "Auxiliary System 2" respectively.

[0138] Ignore the subsequent steps until you reach step 6.

[0139] Step 6: Order n = n +1=11, the first element of the reliability list of the token object R t ,but t The instantaneous reliability of the series-parallel system at time R t =0.9975997735, delete the token object;

[0140] at this time, n ≤ max , return to step 2; otherwise, go to step 7.

[0141] Ignore the subsequent steps until you reach step 7.

[0142] Execute step 7, according to R t The sampling calculation data is summarized to obtain The dynamic reliability change curve of a ship power system within a certain period of time is as follows: Figure 8 As shown, end.

[0143] It will be easily understood by those skilled in the art that the above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims

1. A reliability block diagram of a series-parallel system based on discrete event simulation, characterized in that: The reliability block diagram includes a plurality of reliability components and a plurality of synchronization components, the reliability components correspond one to one with the components in the series-parallel system to be processed, and the connection relationship between the reliability components is connected according to the series or parallel connection relationship between the components in the series-parallel system to be processed; the starting point and the end point of the reliability block diagram are provided with a starting node and an end point node, the series reliability components are directly connected by connecting lines, the front end and the rear end of the module formed by the parallel reliability components are respectively inserted with a node and the reliability components are connected to the nodes by connecting lines, and one node corresponds to one synchronization component; S1 sets the first sampling time, the total number of samplings, the sampling interval, the calculation function of the reliability of each reliability component with respect to time, and sets the number of samplings n = 0; S2 calculates the sampling time corresponding to the current sampling number, calculates the reliability of each reliability component at the current sampling time, assigns a token object to the starting node and initializes the reliability value list of the token object to [1.0], and the token object is passed backward from the starting node; S3 updates the token object in accordance with the preset updating methods of the first module and the second module in sequence; S4 determines whether there is only one token object in the entire reliability block diagram and the token object is on the terminal node, if not, returns to step S3; If yes, the first element of the reliability list corresponding to the current token object is the reliability of the series-parallel system to be processed under the current sampling number, and the token object at the terminal node is deleted; S5 sets n=n+1, determines whether the current sampling number n is greater than the maximum sampling number, if not, returns to step S2; If yes, then output the reliability of the series-parallel system to be processed corresponding to each sampling number and draw a curve of the reliability varying with the sampling number; In step S2, the reliability of each reliability component at the current sampling time is calculated according to the following relationship: R it =f(i,t) Where i is the index of the reliability component, t is the sampling time, and R it is the reliability value of reliability component i at sampling time t, and f(i,t) is the calculation function of the reliability of reliability component i with respect to time t.

2. The reliability block diagram according to claim 1, characterized in that: The number of predecessors and the number of successors of the reliability component are both 1, and the number of predecessors and the number of successors of the synchronization component object are both non-negative integers.

3. The reliability block diagram according to claim 1, characterized in that: In step S3, the first module update method includes the following steps: S311 i=1; S312 determines whether there is a token object on the current reliability component; S313: When there is a token object on the current reliability component, the reliability list is updated according to the preset first relational expression, and the token object on the current reliability component is moved backward, and the process goes to step S314; When there is no token object on the current reliability component, go to step S314; S314 sets i=i+1 and returns to step S312 until the number i of reliability components traversed is greater than the total number of reliability components.

4. The reliability block diagram according to claim 3, characterized in that: In step S313, the first relational expression is performed according to the following relational expression: LP1=LP1*R it Among them, LP1 represents the first element of the reliability list LP of the token object TP on the current reliability component, R it Represents the reliability value of the current reliability component at sampling time t.

5. The reliability block diagram according to claim 1, characterized in that: In step S3, the second module update method includes the following steps: S321 sets j=1; S322 determines whether the number of tokens on the current synchronization component is equal to the number of predecessors; If yes, update the reliability list and the token object on the transfer synchronization component in the manner of the pre-built third module, and proceed to step S323; If not, proceed to step S323; S323 sets j=j+1, and determines whether the number j of traversed synchronization components is greater than the total number of synchronization components. If yes, proceed to step S4; if not, return to step S322.

6. The reliability block diagram according to claim 5, characterized in that: The third module is carried out in the following manner: S3221 calculates the preset value R according to the preset second relational expression t * , determine whether the number of elements in the reliability list LC in the first token object of the current synchronization component is less than 2; If yes, set the preset value R t * Assigning a value to the first element LC1 of the reliability list LC; If not, assign a value to the second element of the reliability list LC according to a preset third relational expression, and delete the first element in the reliability list LC; In this way, an updated reliability list LC is obtained, and all token objects in the current synchronization component are deleted; S3222 calculates the number of successors of the current synchronization component and selects one of the following methods according to the number of successors: When the number of successors is 0, a new token object is created on the current synchronization component, and the reliability list corresponding to the token object is the current reliability list LC; When the number of successors is 1, a new token object is created on the current synchronization component. The reliability list corresponding to the token object is the current reliability list LC, and the token is driven to move to the next position; When the number of successors is greater than 1, a new element 1.0 is inserted before the first element of the current reliability list LC to update the reliability list LC. At the same time, multiple tokens are generated according to the number of successors of the current synchronization component and issued to each successor object respectively. The reliability list corresponding to each token is copied to the current reliability list LC.

7. The reliability block diagram according to claim 6, characterized in that: In step S3221, the preset second relationship is calculated according to the following relationship: R t * =(1-∏ P p＝1 (1-L p,1 )) Where P is the number of token objects contained in the current synchronization component, L p,1 is the reliability list L of the pth synchronization component p The first element of .

8. The reliability block diagram according to claim 7, characterized in that: In step S3221, the preset third relationship is performed according to the following relationship: LC2=LC2*R t * Among them, LC2 is the second element in the reliability list LC.

Images