CC BY
72
RELIABILITY, RISK AND AVAILABILITY ANLYSIS OF A CONTAINER GANTRY
CRANE
Joanna Soszynska-Budny.
Gdynia Maritime University, Gdynia, Poland joannas@am.gdynia.pl
ABSTRACT
The joint model of the system operation process and the system multi-state reliability is applied to the reliability, risk and availability evaluation of the container gantry crane. The container gantry crane is described and the mean values of the container gantry crane operation process unconditional sojourn times in particular operation states are found and applied to determining this process transient probabilities in these states. The container gantry crane different reliability structures in various its operation states are fixed and their conditional reliability functions on the basis of data coming from experts are approximately determined. Finally, after applying earlier estimated transient probabilities and system conditional reliability functions in particular operation states the unconditional reliability function, the mean values and standard deviations of the container gantry crane lifetimes in particular reliability states, risk function and the moment when the risk exceeds a critical value are found. Next the renewal and availability characteristics for the considered gantry crane are determined.
1 INTRODUCTION
Most real technical systems are very complex and it is difficult to analyze their reliability, availability and safety. Large numbers of components and subsystems and their operating complexity cause that the identification, evaluation, prediction and optimization of their reliability, availability and safety are complicated. The complexity of the systems' operation processes and their influence on changing in time the systems' structures and their components' reliability characteristics are very often met in real practice.
Taking into account the importance of the safety and operating process effectiveness of such systems it seems reasonable to expand the two-state approach to multistate approach (Aven 1985, Kolowrocki 2004, Kolowrocki, Soszynska-Budny 2011) in their reliability analysis. The assumption that the systems are composed of multistate components with reliability states degrading in time gives the possibility for more precise analysis of their reliability and operation processes' effectiveness. This assumption allows us to distinguish a system reliability critical state (Kolowrocki 2004, Kolowrocki, Soszynska 2010, Kolowrocki, Soszynska-Budny 2011, Soszynska 2010) to exceed which is either dangerous for the environment or does not assure the necessary level of its operation process effectiveness. Then, an important system reliability characteristic is the time to the moment of exceeding the system reliability critical state and its distribution, which is called the system risk function. This distribution is strictly related to the system multistate reliability function that is basic characteristics of the multi-state system.
The complexity of the systems' operation processes and these processes influence on changing in time the systems' structures and their components' reliability characteristics (Kolowrocki, Soszynska 2010, Kolowrocki, Soszynska-Budny 2010, Kolowrocki, Soszynska-Budny 2011, Soszynska 2010) is often very difficult to fix and to analyse. A convenient tool for solving this problem is semi-Markov (Grabski 2002, Limnios 2001) modelling of the systems operation processes which is proposed in the paper. Using the joint model of the system multi-state reliability and the system semi-Markov operation process (Kolowrocki, Soszynska 2010, Kolowrocki, Soszynska-Budny 2011, Soszynska 2010) it is possible to point out the variability of system
components reliability characteristics by introducing the components' conditional multi-state reliability functions determined by the system operation states. Consequently it is possible to find the system conditional reliability function dependent on the operation states and next it is possible to find its unconditional reliability function and renewal and availability characteristics. This way the obtained results concerned with multi-state systems in its varying in time operation states can be applied to the reliability, risk and availability evaluation of the container gantry crane.
2 CONTAINER GANTRY CRANE SYSTEM ANALYSIS
We analyse the reliability of the container gantry crane that is operating at the container terminal placed at the seashore (Kolowrocki, Soszynska-Budny 2010, Kolowrocki, Soszynska-Budny 2011). The considered container terminal is engaged in trans-shipment of containers. The loading of containers is carried out by using the gantry cranes called Ship-To-Shore (STS). We consider the STS container gantry crane that is composed of 5 basic subsystems S1, S2, S3,
S4 and S5 having an essential influence on its reliability. Those subsystems are as follows:
Sj - the power supply subsystem,
52 - the control and monitoring subsystem,
53 - the arm getting up and getting down subsystem,
54 - the transferring subsystem,
55 - the loading and unloading subsystem.
The gantry crane power supply subsystem S1 consists of:
- a high voltage cable delivering the energy from the substation to the gantry crane E^1-1,
- a drum allowing the cable unreeling during the crane transferring E 2(1),
- an inner crane power supply cable E3(1) ,
- a device transmitting the energy from the high voltage cable to the inner crane cable E 4(1),
- main and supporting voltage transformers E5(1),
- a low voltage power supply cable E6(1),
- relaying and protective electrical components E^ .
The gantry crane control and monitoring subsystem S2 consists of:
- a crane software controller precisely analyzing the situation and takes suitable actions in order to assure correct work of the crane E1(2),
- a measuring and diagnostic device sending signals about the crane state to the software controller
E 22),
- a transmitter of signals from the controller to elements executing the set of commands E3(2),
- devices carrying out the controller's orders (a permission to work, a blockade of work, etc.) E42),
- control panels (an engine room, an operator's cabin, a crane arm cabin) E5(2),
- control and steering cables' connections E62).
The gantry crane arm getting up and getting down subsystem S3 consists of:
- a propulsion unit (an engine, a rope drum, a transmission gear, a clutch, breaks, a rope) E1(3),
- a set of rollers and multi-wheels E 23),
- a crane arm (joints, hooks fastening the arm) E3(3). The gantry crane transferring subsystem S4consists of:
- a driving unit (an engine, a clutch, breaks, a transmission gear, gantry crane wheels) Ej-4-1.
The gantry crane loading and unloading subsystem S5 consists of the winch unitEj(5) composed of:
- a propulsion unit (an engine, a clutch, breaks, a transmission gear, ropes),
- a winch head (which a container grab is connected to),
- a container's grab,
- a container's grab stabilizing unit and the cart unit E 25) composed of:
- a propulsion unit (an engine, a clutch, breaks, a transmission gear, cart wheels, ropes),
- rails which cart is moving on during the operation,
- a crane cart.
The subsystems S1, S2, S3, presented in Figure 1.
S4, S5 are forming a general series gantry crane reliability structure
Si S2
S,
Figure 1. General scheme of gantry crane reliability structure
3 CONTAINER GANTRY CRANE OPERATION PROCESS CHARACTERISTICS PREDICTION
The container gantry crane reliability structure and the subsystems and components reliability depend on its changing in time operation states.
Taking into account expert opinions on the varying in time operation process of the considered container gantry crane we fix the number of the system operation process states v = 6 and we distinguish the following as its six operation states:
• an operation state z1 - the crane standby with the power supply on and the control system off,
• an operation state z2 - the crane prepared either to starting or finishing the work with the crane arm angle position of 90o,
• an operation state z3 - the crane prepared either to starting or finishing the work with the crane arm angle position of 0o,
• an operation state z4 - the crane transferring either to or from the loading and unloading area with the crane arm angle position of 90o,
• an operation state z5 - the crane transferring either to or from the loading and unloading area with the crane arm angle position of 0o,
• an operation state z6 - the containers' loading and unloading with the crane arm angle position of 0o.
Moreover, we fix that there are possible the transitions between all system operation states. To identify all parameters of the container gantry crane operation process the statistical data about this process was needed. All statistical data are collected in (Kolowrocki, Soszynska-Budny 2010). On the basis of statistical data from (Kolowrocki, Soszynska-Budny 2010) the following matrix
[Pbl ] =
0 0.648 0.336 0.008 0 0.008
0.525 0 0.373 0.093 0 0.009
0.105 0.111 0 0 0.118 0.666
0.417 0.583 0 0 0 0
0.005 0 0.220 0 0 0.775
0.012 0 0.628 0 0.360 0
of the probabilities pbl, b,l = 1,2,...,6, of the container gantry crane operation process transitions from the operation state zb to the operation state zl has been fixed.
The mean values Mbl = E[dbl ], b,l = 1,2,...,6, b ^ l, of the container gantry crane operation process
conditional sojourn times at the particular operation states according to (2.12) (Kolowrocki, Soszynska-Budny 2011) are as follows:
M12 = 456.98, M13 = 36.86, M14 = 50, M16 = 3, M21 = 7.89, M23 = 9.12, M24 = 1.55, M26 = 16,
M31 = 5.50, M32 = 4.34, M35 = 6.82, M36 = 7.86, M41 = 2, M42 = 2.14, M51 = 10, M53 = 2.90,
M56 = 24.68, M61 = 22.60. M63 = 23.12, M65 = 20.51.
(1)
After considering the results (1) and applying the formula (2.21) from (Kolowrocki, Soszynska-Budny 2011), the unconditional mean sojourn times of the container gantry crane operation process at the particular operation states are given by:
M = E[01] = 0.648 • 456.98 + 0.336 • 36.86 + 0.008 • 50 + 0.008 • 3 = 308.93,
M2 = E[02] = 0.525 • 7.89 + 0.373 • 9.12 + 0.093 • 1.55 + 0.009 • 16 = 7.83,
M3 = E[63] = 0.105 • 5.50 + 0.111 • 4.34 + 0.118 • 6.82 + 0.666 • 7.86 = 7.09,
M4 = E[64] = 0.417 • 2 + 0.583 • 2.14 = 2.08,
M5 = E[65 ] = 0.005 • 10 + 0.220 • 2.90 + 0.775 • 24.68 = 19.82,
M6 = E[66] = 0.012• 22.60 + 0.628 • 23.12 + 0.360 • 20.51 = 22.17.
Since, according to (2.23) (Kolowrocki, Soszynska-Budny 2011), from the system of equations
= Is,6
I nx + n2 + n3 + n4 +n5 + n6 = 1,
we get
Kx = 0.0951, TT2 = 0.1020, tT3 = 0.3100, TT4 = 0.0102, = 0.1547, TT6 = 0.3280.
Then, the limit values of the transient probabilities pb (t) of the gantry crane operation process at the operation states zb, according to (2.22) in (Kolowrocki, Soszynska-Budny 2011), are given by
p1 = 0.6874, p2 = 0.0187, p3 = 0.0515, p4 = 0.0005, p5 = 0.0717, p6 = 0.1702. (2)
4 CONTAINER GANTRY CRANE IN VARIABLE OPERATION CONDITION RELIABILITY AND RISK EVALUATION
After discussion with experts, taking into account the effectiveness of the operation of the container gantry crane, we fix that the system and its components have four reliability states 0, 1, 2, 3, i.e. z = 3. And consequently, at all operation states zb, b = 1,2,...,6, we distinguish the following reliability states of the system and its components:
• a reliability state 3 - the gantry operation is fully effective,
• a reliability state 2 - the gantry operation is less effective because of ageing,
• a reliability state 1 - the gantry operation is less effective because of ageing and more dangerous,
• a reliability state 0 - the gantry is destroyed.
We assume that there are possible the transitions between the components reliability states only from better to worse ones and we fix that the system and components critical reliability state is r = 2.
Consequently, we assume that the gantry crane subsystems Sv, v = 1,2,...,5 are composed of four-state components, i.e. z =3, with the multi-state reliability functions
[R(v)(t,-)](b)= [1, [RV (t,1)](b),[RV(t,2)](b),[R^-1 (t,3)](b) ], b = 1,2,...,6,
with exponential co-ordinates [R(v)(t,1)](b), [R(v)(t,2)](b), [R(v)(t,3)](b) different in various operation states zb, b = 1,2,...,6.
In (Kolowrocki, Soszynska-Budny 2010), on the basis of expert opinions, the reliability functions of the gantry crane components in different operation states are approximately determined. We will use them in our further system reliability analysis and evaluation.
At the system operation state z1, the container gantry crane is composed of the subsystem S1 which is a series system composed of n = 7 components E(\ i = 1,2,..,7 (subsystems) with the structure showed in Figure 2.
S1
Figure 2. The scheme of the container gantry crane at operation state z1
Thus, at the system operation state z1, the container gantry crane is identical with subsystem S1, that is a four-state series system with its structure shape parameter n = 7 and according to (1.22)-(1.23) (Kolowrocki, Soszynska-Budny 2011), its four-state reliability function is given by the vector
[R(t, -)](1) = [1,[R(t,1)](1), [R(t,2)](1), [R(t,3)](1)], t > 0,
with the coordinates
[ R(t,1)](1)= exp[-0.020t]exp[-0.040t]exp[-0.040t]exp[-0.020t]exp[-0.020t]
exp[-0.018t]exp[-0.033t] = exp[-0.191t], (3)
[ R(t,2)](1)= exp[-0.033t]exp[-0.050t]exp[-0.050t]exp[-0.033t]exp[-0.030t]
exp[-0.028t]exp[-0.040t] = exp[-0.264t], (4)
[R(t,3)](1) = exp[-0.050t]exp[-0.066t]exp[-0.066t] exp[-0.050t]exp[-0.040t]
exp[-0.040t] exp[-0.050t] = exp[-0.362t]. (5)
The expected values of the container gantry crane conditional lifetimes in the reliability state subsets {1,2,3}, {2,3}, {3} at the operation state z1, calculated from the results given by (3)-(5), according to (3.8) (Kolowrocki, Soszynska-Budny 2011), respectively are:
^(1) = 5.24, ^j(2) = 3.79, ^j(3) = 2.76 years.
(6)
At the system operation states z2 and z3, the container gantry crane is composed of the subsystems S1, S2 and S3 forming a series structure shown in Figure 3. The subsystem Sx is a series system composed of n = 7 components E(\ i= 1,2,..,7, the subsystem S2 is a series system composed of n = 6 components Ei(2), i = 1,2,..,6, and the subsystem S3 is a series system composed of n = 3 components Ei(3), i = 1,2,3 .
Figure 3. The scheme of the container gantry crane at operation states z2 and z3
Thus, at the system operation state z2, the container gantry crane is composed of the subsystems S j, S 2 and S3 forming a series structure.
At this operation state, the subsystem S1 is a four-state series system with its structure shape parameter n = 7 and according to (1.22)-(1.23) (Kolowrocki, Soszynska-Budny 2011), its four-state reliability function is given by the vector
[R(1)(t,-}](2) = [1,[R(1)(t,1}](2),[R(1)(t,2}](2), [R(1)(t,3}](2)], t > 0,
with the coordinates
[R(1)(t,1)](2) =exp[-0.020t] exp[-0.040t] exp[-0.040t] exp[-0.020t]
exp[-0.022t] exp[-0.018t] exp[-0.033t] = exp[-0.193t], (7)
[R(1) (t,2)](2) = exp[-0.033t] exp[-0.050t] exp[-0.050t] exp[-0.033t]
exp[-0.027t] exp[-0.028t] exp[-0.040t] = exp[-0.261t], (8)
[R(1)(t,3)](2)= exp[-0.050t] exp[-0.066t] exp[-0.066t] exp[-0.050t]
exp[-0.048t] exp[-0.040t] exp[-0.050t] = exp[-0.370t]. (9)
The subsystem S2 at the operation state z 2, is a four-state series system with its structure shape parameter n = 6 and according to (1.22)-(1.23) (Kolowrocki, Soszynska-Budny 2011), its four-state reliability function is given by the vector
[R(2)(t,-)](2) = [1, [R(2)(t,1)](2), [R(2)(t,2)](2), [R(2)(t,3)](2)], t > 0,
with the coordinates
[ R (2)(t ,1)](2) = exp[-0.053t]exp[-0.048t]exp[-0.048t] exp[-0.048t]
exp[-0.020t]exp[-0.018t] = exp[-0.235t], (10)
[R(2) (t,2)](2) = exp[-0.059t] exp[-0.053t] exp[-0.053t] exp[-0.053t]
exp[-0.025t] exp[-0.029t] = exp[-0.272t], (11)
[R(2)(t,3)](2)= exp[-0.066t] exp[-0.059t] exp[-0.059t] exp[-0.059t]
exp[-0.033t]exp[-0.040t] = exp[-0.316t]. (12)
The subsystem S3 at the operation state z2, is a four-state series system with its structure shape parameter n = 3 and according to (1.22)-(1.23) (Kolowrocki, Soszynska-Budny 2011), its four-state reliability function is given by the vector
[R(3)(t,-)](2) = [1, [R(3)(t,1)](2), [R(3)(t,2)](2), [R(3){t,3}](2)], t > 0,
with the coordinates
[R(3)(t,1)](2) = exp[-0.025t] exp[-0.033t] exp[-0.033t] = exp[-0.091t], (13)
[R(3)(t,2)](2) = exp[-0.040t] exp[-0.040t] exp[-0.040t] = exp[-0.120t], (14)
[R(3)(t,3)](2)= exp[-0.066t] exp[-0.066t] exp[-0.066t] = exp[-0.198t]. (15)
Considering that the container gantry crane at the operation state z2 is a four-state series system composed of subsystems S1 , S2 and S3, after applying (1.22)-(1.23) (Kolowrocki, Soszynska-Budny 2011), its conditional four-state reliability function is given by the vector
[R(t,-)](2) = [1, [R(t,1)](2), [R(t,2)](2), [R(t,3)](2)], t > 0,
with the coordinates
[R(t,1)](2)= exp[-0.193t] exp[-0.235t] exp[-0.091t] = exp[-0.519t], (16)
[R(t, 2)](2)= exp[-0.261t] exp[-0.272t] exp[-0.120t] = exp[-0.653t], (17)
[R(t,3)](2)= exp[-0.370t] exp[-0.316t] exp[-0.198t] = exp[-0.884t]. (18)
The expected values of the container gantry crane conditional lifetimes in the reliability state subsets {1,2,3}, {2,3},{3} at the operation state z2, calculated from the results given by (16)-(18), according to (3.8) (Kolowrocki, Soszynska-Budny 2011), respectively are:
^(1) = 193, ^(2) = 153, ^2(3) = 1.13 year. (19)
After proceeding in the analogous way in the system reliability analysis and evaluation at the remaining operation states z3, z4, z5 and z6, we may determine the system conditional reliability function that are presented below.
„ill ■ET1'
-
Figure 4. The scheme of the container gantry crane at operation states z4 and z5
Figure 5. The scheme of the container gantry crane at operation state z6
At the operation state z3, the container gantry crane conditional reliability function of the system is given by the vector
[R(t, -)](3) = [1, [R(t,1)](3), [R(t,2)](3), [R(t,3)](3) ], t > 0,
with the coordinates
[tf(t,1)](3)= exp[-0.196t]exp[-0.235t] exp[-0.091t] = exp[-0.522t], (20)
[^{t,2}](3)= exp[-0.264t] exp[-0.272t]exp[-0.120t] = exp[-0.656t], (21)
[R(t,3)](3)= exp[-0.375t] exp[-0.316t]exp[-0.198t] = exp[-0.889t]. (22)
The expected values of the container gantry crane conditional lifetimes in the reliability state subsets {1,2,3}, {2,3}, {3} at the operation state z3, calculated from the results given by (20)-(22), according to (3.8) (Kolowrocki, Soszynska-Budny 2011), respectively are:
^(1) = 1.91, ^(2) = 1.52, ^(3) = 1.12 year. (23)
At the system operation states z4 and z5, the container gantry crane is composed of the subsystems S1, S2, S3 and S 4 forming a series structure shown in Figure 4. The subsystem S1 is a series system composed of n = 7 components Ei(1), i = 1,2,..,7, the subsystem S2 is a series system composed of n = 6 components Ei(2), i = 1,2,..,6, the subsystem S3 is a series system composed of n = 3 components Ei(3), i = 1,2,3, and the subsystem S4 consists of a component E1(4). Thus, at the operation state z4 , the container gantry crane conditional reliability function of the system is given by the vector
[R(t, -)](4) = [1, [R(t,1)](4), [R(t,2)](4), [R(t,3)](4)], t > 0,
with the coordinates
[R(t,1)](4) = exp[-0.216t] exp[-0.241t] exp[-0.061t] exp[-0.029t] = exp[-0.547t], (24)
[R(t, 2)](4) = exp[-0.289t] exp[-0.278t] exp[-0.091t] exp[-0.04t] = exp[-0.698t], (25)
[R(t,3)](4) = exp[-0.428t] exp[-0.328t] exp[-0.133t] exp[-0.066t]= exp[-0.955t]. (26)
The expected values of the container gantry crane conditional lifetimes in the reliability state subsets {1,2,3}, {2,3},{3} at the operation state z4, calculated from the results given by (24)-(26), according to (3.8) (Kolowrocki, Soszynska-Budny 2011), respectively are:
^(1) = 1.83, ^(2) = 1.43, ^(3) = 1.05 year. (27)
At the operation state z5, the container gantry crane conditional reliability function of the system is given by the vector
[R(t,-)](5) = [1,[R(t,1}](5), [R(t,2)](5), [R(t,3)](5)], t > 0,
with the coordinates
[ R(t,1)](5) =exp[-0.216t]exp[-0.241t]exp[-0.061t]exp[-0.025t] = exp[-0.543t], (28)
[ R(t,2)](5)=exp[-0.289t]exp[-0.278t] exp[-0.091t]exp[-0.029t] = exp[-0.687t], (29)
[ R(t ,3)](5)=exp[-0.428t]exp[-0.328t]exp[-0.133t]exp[-0.050t] = exp[-0.939t]. (30)
The expected values of the container gantry crane conditional lifetimes in the reliability state subsets {1,2,3}, {2,3},{3} at the operation state z5, calculated from the results given by (28)-(30), according to (3.8) (Kolowrocki, Soszynska-Budny 2011), respectively are:
U5(1) = 1.84, u5(2) = 1.46, ¡u5(3) = 1.06 year.
(31)
At the system operation state z 6, the container gantry crane is composed of the subsystems Sx, S2, S3 and S5 forming a series structure shown in Figure 5. The subsystem Sx is a series system composed of n = 7 components E(\ i = 1,2,..,7, the subsystem S2 is a series system composed of n = 6 components Ei(2), i = 1,2,..,6, the subsystem S3 is a series system composed of n = 3 components Ei(3), i = 1,2,3 and the subsystem S5 is a series system composed of n = 2 components E?\ i = 1,2.
Thus, at the operation state z6, the container gantry crane conditional reliability function of the system is given by the vector
[R(t, -)](6) = [1,[R(t,1)](6), [R(t,2)](6), [R(t,3)](6)], t > 0,
with the coordinates
[ R(t,1)](6)=exp[-0.201t]exp[-0.250t]exp[-0.087t]exp[-0.080t] =exp[-0.618t], (32)
[R(t, 2)](6)= exp[-0.273t] exp[-0.29t] exp[-0.12t] exp[-0.1t]= exp[-0.783t], (33)
[R(t,3)](6)= exp[-0.396t] exp[-0.337t]exp[-0.16t] exp[-0.132t]= exp[-1.025t]. (34)
The expected values of the container gantry crane conditional lifetimes in the reliability state subsets {1,2,3}, {2,3},{3} at the operation state z6, calculated from the results (32)-(34), according to (3.8) (Kolowrocki, Soszynska-Budny 2011), respectively are:
U6(1) = 1.62, u6(2) = 1.28, u6(3) = 0.98 year.
(35)
In the case when the operation time is large enough the unconditional four-state reliability function of the container gantry crane is given by the vector
R(t, •) = [1, R(t,1), R(t,2), R(t,3)], t >0, (36)
where according to (3.5)-(3.6) (Kolowrocki, Soszynska-Budny 2011) and considering (2), the vector coordinates are given respectively by
R(t,1) = 0.6874 • [R(t,1}](j) + 0.0187 • [R(t,1}](2) + 0.0515 • [R(t,1}](3) + 0.0005 • [R(t,1}](4) + 0.0717 • [R(t,1}](5) + 0.1702 • [R(t,1}](6) for t > 0,
( 3 7 )
R(t,2) = 0.6874 • [R(t,2}](1) + 0.0187 • [R(t,2}](2) + 0.0515 • [R(t,2}](3) + 0.0005 • [R(t,2}](4)
+ 0.0717 • [R(t,2}](5) + 0.1702 • [R(t,2}](6) for t > 0, (38)
R(t,3) = 0.6874 • [R(t,3}](1) + 0.0187 • [R(t,3}](2) + 0.0515 • [R(t,3}](3) + 0.0005 • [R(t,3}](4)
+ 0.0717 • [R(t,3}](5) + 0.1702 • [R(t,3}](6) for t > 0, (39)
and the coordinates [R(t,u}](b), b = 1,2,...,6, u = 1,2,3, are given by (3)-(5), (16)-(18), (20)-(22), (24)-(26), (28)-(30), (32)-(34). The graphs of the coordinates of the container gantry crane reliability function are presented in Figure 6.
0 • V ...............................................r 7 7 = 5= = i
0 2 A 6 8 10 t 12
Figure 6. The graph of the container gantry crane reliability function R(t, •) coordinates
The expected values and standard deviations of the container gantry crane unconditional lifetimes in the reliability state subsets {1,2,3}, {2,3},{3) calculated from the results given by (37)-(39), according to (3.7)-(3.9) (Kolowrocki, Soszynska-Budny 2011) and considering (2), (6), (19), (23), (27), (31), (35), respectively are:
H(1) = 0.6874 • 5.24 + 0.0187 • 1.93 + 0.0515 • 1.91 + 0.0005 • 1.83 + 0.0717 • 1.84
+ 0.1702 • 1.62 = 4.14 years, (40)
cr (1) = 4.71years,
ju(2) = 0.6874 • 3.79 + 0.0187 • 1.53 + 0.0515 • 1.52 + 0.0005 • 1.43 + 0.0717 • 1.46
+ 0.1702 • 1.28 = 3.04 years, (41)
cr(2) = 3.43 years, (42) ju(3) = 0.6874 • 2.76 + 0.0187 • 1.13 + 0.0515 • 1.12 + 0.0005 • 1.05 + 0.0717 • 1.06
+ 0.1702 • 0.98 = 2.22 years, (43) cr(3) = 2.50 years,
Further, considering (3.10) from (Kolowrocki, Soszynska-Budny 2011) and (40), (41), (43), the mean values of the unconditional lifetimes in the particular reliability states 1, 2, 3 respectively are:
£(1) = M1) - M2) = 1.10, £(2) = ¡d(2) - n(3) = 0.82, £(3) = ^(3) = 2.22 years. (44)
Since the critical reliability state is r = 2, then the system risk function, according to (3.11) (Kolowrocki, Soszynska-Budny 2011), is given by
r(t) = 1 -R(t, 2), (45)
where R(t ,2) is given by (38).
Hence, the moment when the system risk function exceeds a permitted level, for instance 5 = 0.05, from (3.12) (Kolowrocki, Soszynska-Budny 2011), is
t = r_1(5) = 0.126 year.
(46)
The graph of the risk function r (t) of the container gantry crane operating at the variable conditions is given in Figure 7.
i
0.8 0.6 -0.4 0.2 0
0 2 4 6 8 10 t 12
Figure 7. The graph of the container gantry crane risk function r (t)
5 CONTAINER GANTRY CRANE AVAILABILITY PREDICTION
Using the results of the container gantry crane reliability prediction given by (41)-(42) and the results of the classical renew theory presented in (Kolowrocki, Soszynska-Budny 2011), we may predict the renewal and availability characteristics of this system in the case when it is repairable and its time of renovation is either ignored or non-ignored.
First, assuming that the container gantry crane is repaired after the exceeding its reliability critical state r = 2 and that the time of the system renovation is ignored we obtain the following results:
a) the time SN (2} until the Nth exceeding by the system the reliability critical state r = 2, for sufficiently large N, has approximately normal distribution N (3.04 N,3.43VN), i.e.,
t - 3 04N
F(N)(t,2} = P(Sn (2} < t) = Fn№1)( 3 ^ ), t E (-x, x);
b) the expected value and the variance of the time SN (2} until the Nth exceeding by the system the reliability critical state r = 2 are respectively given by
E[Sn (2}] = 3.04N, D[SN (2}] = 11.76N;
c) the number N (t,2) of exceeding by the system the reliability critical state r = 2 up to the moment t, t > 0, for sufficiently large t, approximately has the distribution of the form
^ ^ ,3.04(N +1)-1 „ 3.04N-tN ir _
P(N(t'2) = N) s F'(",j)(-jkjГ-) - Fn(»,1)(^}, n =0'1''.;
d) the expected value and the variance of the number N (t,2) of exceeding by the system the reliability critical state r = 2 up to the moment t, t > 0, for sufficiently large t, approximately are respectively given by
H(t,2) = 0.3289?, D(t,2) = 0.419t.
Further, assuming that the container gantry crane is repaired after the exceeding its reliability critical state r = 2 and that the time of the system renovation is not ignored and it has the mean value yU0(2) = 0.0027 and the standard deviation cr0(2} = 0.0014, we obtain the following results:
a) the time SN (2) until the Nth exceeding by the system the reliability critical state r = 2, for sufficiently large N, has approximately normal distribution
N(3.04N + 0.0027(N -1)^11.76N - 0.0000019(N -1)), i.e.,
F (N) (t,2) = P(Sn (2) < t) . Fn(01), >- 3.0427N + qoo27 ), ( E (-x, x),
N V11.760002N - 0.0000019
b) the expected value and the variance of the time SN (2) until the Nth exceeding by the system the reliability critical state r = 2, for sufficiently large N, are respectively given by
E[Sn (2)] = 3.04N + 0.0027(N -1), D[SN (2)] = 11.76N + 0.0000019(N -1) ;
c) the number N (t,2) of exceeding by the system the reliability critical state r = 2 up to the moment t, t > 0, for sufficiently large t, has approximately distribution of the form
, , ^ ,3.0427( N +1) -1 - 0.0027 N P( N (t,2) = N) = Fn (0 j) (-^ ■ '-)
1.966V t + 0.0027
3 0427N -1 - 0 0027 -F (3 42 0 0027) n = 01 •
1 N(0,1) V I--JV "jV-J
( 1.96Wt + 0.005
d) the expected value and the variance of the number N (t,2) of exceeding by the system the reliability critical state r = 2 up to the moment t, t > 0, for sufficiently large t, are respectively given by
H(t,2) = 0.329(t + 0.0027), D(t,2) = 0.417(t + 0.0027);
e) the time SN (2) until the Nth system's renovation, for sufficiently large N, has approximately normal distribution N(3.0427N,3.429a/N), i.e.,
F(N) (^,2) = P(SN (2) < t) 3 Fn(01) (t 3340427N)), 1 G
f) the expected value and the variance of the time SN (2) until the Nth system's renovation, for sufficiently large N, are respectively given by
E[Sn (2)] = 3.0427N, D[Sn (2)] = 11.76002N;
g) the number N(t,2) of the system's renovations up to the moment t, t > 0, for sufficiently large t, has approximately distribution of the form
P(N(t 2) = N) 3 F (3 0427(N +1 31 ) 3f (3 0427N31 ) N = 01 • P(N (t,2) N) 3 Fn(01)( 1.96^/F j FN(01)( 1.966V1 N 0
h) the expected value and the variance of the number N(t,2) of system's renovations up to the moment t, t > 0, for sufficiently large t, are respectively given by
H (t,2) = 0.3286t, D (t,2) = 0.417t;
i) the steady availability coefficient of the system at the moment t, t > 0, for sufficiently large t, is given by
A(t,2) = 0.9989, t > 0;
j) the steady availability coefficient of the system in the time interval < t, t + r),r > 0, for sufficiently large t, is given by
x
A(t,T,2) = 0.3291 R(t,2)dt, t > 0, r > 0,
where R(t, 2) is given by (7.97).
6 CONCLUSION
In the paper the multi-state approach to the analysis and evaluation of systems' reliability and risk has been practically applied. The container gantry crane has been considered at varying in time operation conditions. The system reliability structure and its components reliability functions were changing at variable operation conditions. The paper proposed an approach to the solution of practically very important problem of linking the systems' reliability and their operation processes. To involve the interactions between the systems' operation processes and their varying in time reliability structures a semi-markov model of the systems' operation processes and the multi-state system reliability functions were applied. This approach gives practically important in everyday usage tool for reliability evaluation of the systems with changing reliability structures and components reliability characteristics during their operation processes what exemplary was illustrated in its application to the gantry crane.
The characteristics of the gantry crane operation process are of high quality because of the very good statistical data necessary for their estimation. Unfortunately, the reliability characteristics of the gantry crane components are evaluated on non sufficiently exact data coming from experts and concerned with the mean values of the components lifetimes only that because of the complete lack of statistical data about their failures are strongly lowered. Also, the system and its components reliability states are defined on a high level of generality and should be described more precisely. All these inaccuracies causes that the evaluation of the gantry crane reliability, risk and availability characteristics should be consider as an illustration of the proposed approach application.
Acknowledgements
The paper describes part of the work in the Poland-Singapore Joint Research Project titled "Safety and Reliability of Complex Industrial Systems and Processes" supported by grants from the Poland's Ministry of Science and Higher Education (MSHE grant No. 63/N-Singapore/2007/0) and the Agency for Science, Technology and Research of Singapore (A*STAR SERC grant No. 072 1340050).
7 REFERENCES
Aven, T. 1985. Reliability evaluation of multi-state systems with multi-state components. IEEE
Transactions on reliability 34, pp 473-479. Aven, T., Jensen, U. 1999. Stochastic Models in Reliability. Springer-Verlag, New York. Grabski, F. 2002. Semi-Markov Models of Systems Reliability and Operations. Warsaw: Systems
Research Institute, Polish Academy of Science. Kolowrocki, K. 2004. Reliability of Large Systems. Elsevier, Amsterdam - Boston - Heidelberg -London - New York - Oxford - Paris - San Diego - San Francisco - Singapore - Sydney -Tokyo.
Kolowrocki, K., Soszynska, J. 2010. Integrated Safety and Reliability Decision Support System -IS&RDSS. Tasks 10.0-10.15 in WP10: Safety and Reliability Decision Support Systems for Various Maritime and Coastal Transport Sectors. Poland-Singapore Joint Research Project. MSHE Decision No. 63/N-Singapore/2007/0. Gdynia Maritime University.
Kolowrocki, K., Soszynska-Budny, J. 2010. The container gantry crane operation, reliability, risk, availability identification, prediction and optimization - Testing IS&RDSS. Task 9.5 in WP9: Applications and Testing of Packages of Tools in Complex Maritime Transportation Systems and Processes Safety and Reliability Optimization. Poland-Singapore Joint Research Project. MSHE Decision No. 63/N-Singapore/2007/0.Gdynia Maritime University.
Kolowrocki, K., Soszynska-Budny, J. 2011. Reliability and safety of complex technical systems and processes. Modeling-Identification-Prediction-Optimization. Springer.
Limnios, N., Oprisan, G. 2001. Semi-Markov Processes and Reliability. Birkhauser, Boston.
Soszynska, J. 2010. Reliability and risk evaluation of a port oil pipeline transportation system in variable operation conditions. International Journal of Pressure Vessels and Piping Vol. 87, No. 2-3, 81-87.
