Next Article in Journal
Short-Term Photovoltaic Power Prediction Based on a Digital Twin Model
Previous Article in Journal
Underwater Image Enhancement Based on Light Field-Guided Rendering Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulation Model of Hydraulic System States for Ship Cranes

Maritime Department, University of Dubrovnik, Ćira Carića 4, 20000 Dubrovnik, Croatia
J. Mar. Sci. Eng. 2024, 12(7), 1218; https://doi.org/10.3390/jmse12071218
Submission received: 16 May 2024 / Revised: 5 July 2024 / Accepted: 7 July 2024 / Published: 19 July 2024
(This article belongs to the Section Ocean Engineering)

Abstract

:
The aim of this research is to devise a continuous simulation model for predicting ship crane failures to increase their reliability and reduce unplanned downtime during cargo loading and unloading operations. To predict the condition of the hydraulic system, a database from the GALIOT software package was used for carrying out maintenance on cranes at m/v “O” over a period of 120,000 working hours. In the research, fault tree analysis (FTA) was used to identify causal relationships between system failures and basic events, while the Markov mathematical model was used to model the system state and predict transitions between different failure states. A system dynamics simulation model was developed to simulate the behavior of a system using POWERSIM PowerSim Constructor 2.5.d (4002), and regression analysis was performed to analyze the simulation results and understand the relationships between dependent and independent variables. The results show that a model for predicting failures in the hydraulic motors and pumps of ship cranes was developed, and the Markov model makes it possible to estimate the frequency of transitions between states under the condition that the sum of reliability equals one. The simulation model shows high reliability of the cranes and a constant frequency of failures throughout the 120,000 operating hours, while the regression analysis confirms the validity of the simulation model and shows a strong correlation between the analyzed variables. These models are used to improve the planning of ship crane maintenance, reduce unplanned downtime, and predict and promptly detect failures, which overall minimizes maintenance costs and failures.

1. Introduction

Designing a bulk carrier requires taking into consideration the type and number of deck cranes to be installed onboard. The cranes should be reliable and have a minimal possibility of failure, while the ship crew and superintendent are expected to predict possible failures and prevent them.
Papers [1,2,3,4] are typical papers for this area of research. Special attention is given to researching failures occurring during the exploitation of ship cranes; hence, the paper [1] offers recommendations to avoid unforeseen failures, damages, and breakdowns. Suggestions for proper and timely detection of possible failures and maintenance are also given. The main focus is on safety during cargo operations, maintenance, and the repairing of ship cranes. The literature review reveals additional relevant topics, as shown in the paper [2], dealing with training, reliability, maintenance, handling, control systems, mechanical systems, structural damage, and the inspection of ship cranes.
The Markov chain theory or the reliability model in the paper [3] is used to evaluate associated risks and expenses resulting from delays due to system failures. The Markov model not only evaluates economic consequences (costs) due to crane failures but also helps draw attention to critical aspects during the exploitation of cranes. Some methods that incorporate efficient solutions are a combination of mathematical models obtained from Markov chains and system dynamics, which have been proven to be excellent for monitoring and predicting the system status, as shown in the paper [4].
One of the maintenance decision-making models is the Integrated Maintenance Decision Making Model (IMDMM) for the exploitation of cranes and is based on the concept of the Integrated Decision Making Model [5]. This model aims to improve operational efficiency and crane operation, as well as reduce the risk of gantry crane inefficiency (GCI) on the basis of the digital twins (DT) concept for carrying out maintenance. This paper [5] further evaluates the risk of crane operation based on the Monte Carlo simulation model elaborated using Markov chains, while the IMDMM was used for optimization and supported by optimization algorithms. When endeavoring to meet planning requirements for fuel and lubricating oil systems in terms of maintenance schedules for container cranes (CC) as described in the paper [6], where maintenance planning is of extreme importance and should be aligned with the CC work schedule, with due attention paid to avoiding maintenance activities that interfere with cargo loading and unloading schedules. This paper utilizes the system dynamics modeling method to model system performance and determine the minimum reliability and availability index values, which stem from the schedule for proper container crane maintenance. The aim of many papers addressing this area of research, among which is the paper [7] presenting the crane maintenance model aiming to avoid delays, is maximum safety of the crew who handle ship cranes and tailoring of maintenance to the crew who carry out maintenance works. The primary aim of the paper [7] is based on crane maintenance in ports performed by ship crews, with great reliance on crew education and training, which in turn enables them to solve issues promptly and efficiently.
This working concept (Figure 1) ensures that causes are detected and dealt with while preventing or minimizing the consequences. The primary goal is to elaborate a continuous simulation model with feedback within a given period of time, which would anticipate all system states in a given future time, using the methods and analysis shown in Figure 1.
Unlike the previous solutions, this paper includes more methods and analyses (Figure 1), offering a more acceptable solution to defining, monitoring, analyzing, and optimizing preventive maintenance of ship cranes. Besides optimization, education, and training, the derived model also ensures the safety of the crew working on the cranes.
This paper defines and elaborates a system dynamics simulation model capable of anticipating expected ship crane states in the future, thus, preventing major failures or breakdowns of ship cranes. Meeting this goal requires identifying the advantages and quality obtained by implementing the simulation model using the Markov mathematical model, which, in tandem with system dynamics, enables monitoring of ship crane states at any time.
This paper adopted the FTA method utilizing data from the ship crane maintenance software GALIOT to anticipate possible future states of ship crane systems, i.e., states causing breakdowns [8,9]. The FTA method was used as the basis for working out the Markov model, which, in turn, was used to construe the system dynamics simulation model. The results obtained were analyzed using simulation and regression analysis. The sequence of methods and analyses is shown in Figure 1.
Given that ship cranes operate in harsh environments, there is definitely a need to make timely predictions as to when a failure may occur, thus, ensuring smooth crane operation during loading and unloading operations. A combination of all methods to be used in the paper assisted in providing a quality model and provided better possibilities for further application and upgrades.

2. Methods Analyzing Ship Crane States

Ensuring the best analyses of data collected and defining the states as precisely as possible relied on various methods for deriving a quality simulation model.

2.1. Crane State (Failure) Analysis Using the FTA Method

The FTA method uses deductive analysis, which studies failures and identifies relationships between system failures and top events (i.e., detecting causes of events). It is a system state analysis technique (i.e., unwanted system states), and when applying the FTA method, defining the top event and its cause is of utmost importance.
The majority of failures on ship cranes can be anticipated before they occur and, therefore, prevented. Some failures, however, cannot be foreseen because they are accidental failures and happen during the exploitation of ship cranes.
Fault tree analysis (FTA) is a method for diagnosing and analyzing the safety and reliability of complex ship systems. When using the FTA method, such systems are presented in a diagram to evaluate the reliability of ship cranes while operating.
This paper utilizes the FTA method to construct a diagram, i.e., a fault tree, that models the states causing a top fault, i.e., unwanted breakdown or failure. The fault tree provides a graphic illustration of relationships between events and their impact on the top event. This analysis incorporates a deductive method, which is necessary for ship system reliability models as it monitors the frequency of failures due to top events, i.e., the sequence of events to the very bottom.
Constructing the fault tree requires the following steps:
  • Defining and describing the ship system;
  • Becoming familiar with the limitations of the system and goals;
  • Defining the top event;
  • Utilizing the ship system database, which usually includes a history of the system state from the maintenance program package used to construct the fault tree.
The fault tree consists of an event symbol, logic gates, and symbols for transferring the diagram structure so as to further develop the tree at another location or page. The logic symbols are used to define the interdependence and relationship between events. An example would be the OR gate used if an input event occurs and the AND gate if all the input events occur. If the fault tree is complex, it is divided into several smaller trees that describe the subsystems of the complex ship system, making analysis and interpretation easier.
The paper analyses the target event on the part of a ship crane, involving:
  • failure of the hydraulic engine, and
  • failure of the hydraulic pump.
These two failures in the hydraulic system are analyzed due to their complexity, sensitivity, and key importance in the operation of the ship crane. Every complex ship system, including the hydraulic system, due to the sensitive nature of the components and exposure to the sea environment, is prone to damage and may subsequently cause delays in the ship crane performing unloading and loading operations, including financial losses. A solution to prevent or minimize these losses is given in this paper.
Having analyzed the system and its states and using the FTA method, it becomes an ideal situation to construct the Markov model.

2.2. Markov Mathematical Model

The Markov model is often used to analyze system states as, at any moment, the system can remain in the same state or transit into another state, which is called transition. When using the Markov model, the possibilities of the model are solely dependent on the present state and not on the past states, which is why they are commonly called models “without memory” [4]. The greatest advantage and best use of the Markov model is the possibility of creating a large variety of possible system states.

2.3. System Dynamics Simulation Model

Nowadays, system dynamics has become widespread in science as an efficient simulation tool when combined with simulation software such as PowerSim Constructor 2.5.d (4002) POWERSIM [10,11,12,13,14,15,16]. The actions used to elaborate the system dynamics model are given in the following order:
  • Defining the problem;
  • Conceptualization of the system;
  • Constructing the model in the POWERSIM PowerSim Constructor 2.5.d (4002) simulation program;
  • The behavior of the model when changing the initial conditions;
  • Testing the validity of the model;
  • Analysing the obtained results.
The model is described using quantitative and qualitative simulation modeling, leading to the optimization of the dynamic system and the entire process, which is completely based on a mental-verbal model.
System dynamics are utilized in this paper to define the following models for further research:
  • Mental-verbal model;
  • Structural model;
  • Mathematical model;
  • Simulation model using the POWERSIM PowerSim Constructor 2.5.d (4002) computer program.
Modeling the system dynamics was achieved using the visual modeling tool POWERSIM to explore the system further.

2.4. Regression Analysis of the Results Obtained from the System Dynamics Simulation Model

Regression analysis is used to follow up and deduce why the dependent variable is affected by changes in other independent variables or factors [11,12,17]. Regression analysis estimates the relationship between a dependent variable and one or more independent variables, i.e., the regression function. Simple, multiple, linear, and non-linear regressions are used in regression analysis, and in this paper, in particular, a linear regression is used where the regression function is defined with a final number of parameters. The correlation between the two variables is described as:
  • Linear correlation: y = ax + b.
  • Curved line correlation: y = aebx.
In regression analysis, the dispersion diagram is used as a goal, and on the basis of the dispersion, a conclusion about the form is made (defines the mathematical function representing the deterministic part of the regression model), the direction (indicates whether the connection is positive or negative), and the strength of the connection (analyses random variables of the regression model).
In this paper, regression analysis was used to explore correlations between variables obtained from the simulation model, i.e., it analyses the relationship between the dependent variable and one or more independent variables, thus, giving the mathematical equation. A linear and exponential correlation was used to determine the relationship between the dependent and independent variables obtained from simulation model data analysis as it best fits the simple regression model.

3. Results and Discussion

3.1. FTA Method for Analyzing Ship Crane Failures

Data from the maintenance program package installed on board was used to describe complex ship systems, including measurement values and system statuses updated daily. For example, the ship crane on board m/v “O” was used, where its top was “SHIP’S CRANE OUT OF ORDER”, while a fault tree was constructed for major failures that put the crane out of order.
The ship crane fault tree is divided into two main systems where failures most frequently occur and put the crane out of order. These failures are (1) faults on the hydraulic engine and (2) faults on the hydraulic pump, as shown in Figure 2.
Frequent failures putting the crane out of order relate to the hydraulic system or the hydraulic motor (Figure 3) and the hydraulic pump (Figure 4), whereas failures rarely occur on the actual crane structure.
On these systems, as shown in Figure 3, Figure 4 and Figure 5, the faults can be total, when the crane cannot fulfill its required purpose, or partial, when crane operation is compromised.
Figure 3 shows the total faults causing faults on the hydraulic motor and consequently putting the crane out of order. Total faults include:
  • Damage to the hydraulic motor housing;
  • Blocked control valve;
  • Broken hydraulic motor vanes;
  • Blocked safety valve;
  • Damaged hydraulic pump gears.
Figure 4 shows only a particular total fault, i.e., damaged gears, which halted the operation of the entire hydraulic system and the crane itself.
The states mentioned in Figure 2, Figure 3 and Figure 4 are used to make a fault tree incorporating total faults affecting the top event, i.e., failure of the ship crane (Figure 5).
To work out the mathematical model using Markov chains, the fault tree with a total fault on the ship crane will be used (Figure 5).

3.2. Mathematical Model Described Using the Markov Model

The Markov model was elaborated by defining and analyzing data obtained while monitoring state transition dynamics (from one state to another) of the selected systems on ship cranes over a given period of time.
Marks (0, 1, 2, 3, 4, 5, 6) on the Markov model in Figure 6 have the following meanings:
  • Ship crane state: “in operation”;
  • Out of order state: “control valve blocked”;
  • Out of order state: “vanes on hydraulic motor broken”;
  • Out of order state: “safety valve blocked”;
  • Out of order state: “hydraulic motor housing damaged”;
  • Out of order state: “hydraulic pump gears damaged”;
  • Out of order state: “ship crane out of order”.
Figure 6. Markov model for ship cranes on m/v “O”.
Figure 6. Markov model for ship cranes on m/v “O”.
Jmse 12 01218 g006
The frequency of transition from one state to another state until reaching a total fault in the ship crane interrupts unloading operations as well as returning to the state “in operation” (repair) [4,10,11,12,14].
Frequency of transition from state “0” to “out of operation” (Equation (1)):
λ0n = Lam_0n     n = 1, 2, 3, 4, 5
Frequency of transition from “out of order” to “in operation” “0” (Equation (2)):
μn0 = Mu_n0     n = 1, 2, 3, 4, 5, 6
Frequency of transition from “out of order” to “6”, indicating the ship crane is out of order (Equation (3)):
αn6 = Alpha_n6     n = 1, 2, 3, 4, 5
Reliability of the technical system (Equation (4)) [4,8,9]:
Rn(t) = e−λ0nt     n = 0, 1, 2, 3, 4, 5, 6.    where n = state number
A ship crane is a complex technical system consisting of systems and subsystems, each directly affecting the operation of the crane. Should any of the systems or subsystems fail, the crane is immediately out of order. The paper describes the two most important systems, i.e., the hydraulic motor and hydraulic pump.
Besides the hydraulic motor and pump, the hydraulic system also consists of additional components playing a smaller part in eventual crane failure because they can easily be replaced, or the crane can operate at a reduced capacity. Each “out of operation” state (1–5) directly affects the crane and leads to the crane becoming out of order (6). After repairing the hydraulic system, the ship crane moves from the “out of order” state and returns to the 0 or “in operation” state.
Figure 6 shows the Markov mathematical model with states of ship crane systems and subsystems during exploitation. The Markov mathematical model shows system states in the form of linear differential Equations (4)–(11).
State 0:
d P 0 ( t ) d t = ( λ 01 + λ 02 + λ 03 + λ 04 + λ 05 ) P 0 ( t ) + + μ 10 P 1 ( t ) + μ 20 P 2 ( t ) + μ 30 P 3 ( t ) + μ 40 P 4 ( t ) + μ 50 P 5 ( t ) + μ 60 P 6 ( t )
State 1:
d P 1 ( t ) d t = λ 01 P 0 ( t ) ( μ 10 + α 16 ) P 1 ( t )
State 2:
d P 2 ( t ) d t = λ 02 P 0 ( t ) ( μ 20 + α 26 ) P 2 ( t )
State 3:
d P 3 ( t ) d t = λ 03 P 0 ( t ) ( μ 30 + α 36 ) P 3 ( t )
State 4:
d P 4 ( t ) d t = λ 04 P 0 ( t ) ( μ 40 + α 46 ) P 4 ( t )
State 5:
d P 5 ( t ) d t = λ 05 P 0 ( t ) ( μ 50 + α 56 ) P 5 ( t )
State 6:
d P 6 ( t ) d t = α 16 P 1 ( t ) + α 26 P 2 ( t ) + α 36 P 3 ( t ) + α 46 P 4 ( t ) + α 56 P 5 ( t ) μ 60 P 6 ( t )
Initial conditions determined at the time t = 0 are given in Equation (12):
P 0 ( 0 ) = 1 ; P 1 ( 0 ) = 0 ; P 2 ( 0 ) = 0 ; P 3 ( 0 ) = 0 ; P 4 ( 0 ) = 0 ; P 5 ( 0 ) = 0 ; P 6 ( 0 ) = 0
The condition defined by the identity equation (Equation (13)) must be fulfilled at all times:
P u k = P 0 + P 1 + P 2 + P 3 + P 4 + P 5 + P 6 = 1
The stationary solution to the linear differential equations for the system and relating to the Markov model [4,10,11,13] is required (Equation (14)):
d P n ( t ) d t = 0 ; n = 0 , 1 , 2 , 3 , 4 , 5 , 6
The solution to the Markov model for a stationary process state, representing a system of linear equations with unknown factors, is used to determine probabilities for stationary ship cargo crane states, as shown in Equations (15)–(17).
P 0 ( t ) = μ 60 i = 1 5 λ 0 i μ 60 + α i 6 μ i 0 + α i 6 i = 0 , 1 , 2 , 3 , 4 , 5
P 6 ( t ) = 1 μ 60 α 16 P 1 ( t ) + α 26 P 2 ( t ) + α 36 P 3 ( t ) + α 46 P 4 ( t ) + α 56 P 5 ( t )
P j ( t ) = λ 0 n P 0 ( t ) μ n 0 + α n 6 ; j = 0 , 1 , 2 , 3 , 4 , 5 ; n = 0 , 1 , 2 , 3 , 4 , 5
When using the program package GALIOT [8,9] for the maintenance of a ship diesel engine, the model parameters are calculated, as well as the number of transitions from state into state and time in individual states of a ship crane.
The obtained data are used to calculate mean times:
T ¯ i j —mean time in state i until transitioning to state j.
To calculate the frequency of transitions from state to state, the table containing calculated mean times for individual states of a ship crane is used. The mean time is calculated using the frequency of failures, i.e., the frequency of transitions from state to state. Based on this exponential distribution, the frequencies of failures λ(t) or transitions are constant.
The frequency of failures λ(t) is constant and is calculated using Equation (18), where the frequency of transitions from the “in operation” state to the “out of order” state (λ0n for n = 1, 2, 3, 4, 5, and 6) is shown (Equation (18)) [4].
λ ( t ) = 1 T ¯ u r = λ 0 n
The frequency of transitions from any “out of order” state into a “crane out of order” state, given that the failures of these systems halt crane operation when any system is “out of order” α(t) and does not return into the “in operation” state, is calculated using Equation (19) [4,14], where: αn6 = 1, 2, 3, 4, and 5.
The frequency of transitions from any “out of order” state into a “crane out of order” state, given that failures in these systems halt crane operations when any system is “out of order” and does not return to the “in operation” state, is calculated using Equation (19), where:
α ( t ) = 1 T ¯ s t = α n 6
Frequency of repair μ(t) is constant and is calculated using Equation (20) [4]:
μ ( t ) = 1 T ¯ u k = μ n 0 n = 6
Tsrur_n, or mean time in operation until the system enters state n (number of states n = 1, 2, 3, 4, 5, and 6), is defined based on system reliability (Equation (4)) [4,10,11,12] and the frequency of system failures λ0n:
T s r u r _ n = T ¯ u r _ n = 0 R n ( t ) d t 0 e λ 0 n t d t , T s r u r _ n or T ¯ u r _ n
Table 1 shows the frequency of transitions from state to state in case of a failure λ0n, transition to the state “in operation” μ0n, total failure, crane state “out of order” αn6, and probabilities for each state P on the ship cranes. The values were obtained from the GALIOT program package database for the m/v “O” over a period of 80.000 working hours. Consequently, the scenario accounting for 120.000 working hours for the simulation was given to predict the reliability of the ship cranes in the respective period.
All calculations shown in Table 1 were done using time in operation t0n, time out of operation tn0, and time of transition into ship crane total failure state tn6, where n = 1, 2, 3, 4, 5, and 6. When describing the transition from state to state (out of order, in operation, and total failure of ship crane), failures must be random and independent. One of the conditions is that values of time until system failure and time duration of repairs are known. Monitoring over a prolonged period of time indicates that the time values are actually mean times (Equation (22)) to calculate the frequency. This approach is applied in each period, given that components can be renewed [4]:
T ¯ 0 n = t 0 n n 0 n
The same is obtained for the mean time of the “out of order” state (which is the mean time for repairs in the given interval) (Equation (23)) [4]:
T ¯ n 6 = t n 6 n n 6  
Parameters for frequencies of transitions from one state to another in the Markov model are calculated using data from Table 1 and using the given (for frequency of transitions from state to state) Equations (10)–(14).
The values for individual states are calculated using the Equations (15)–(19).

3.3. System Dynamics Simulation Model for Ship Crane Reliability

By applying the mathematical model described in Markov chains for reliability analysis, a dynamic structural model was devised and presented in Figure 7a. It shows how the reliability of a hydraulic pump and hydraulic motor on crane systems was monitored, given that these two systems are the most crucial systems affecting the reliability of ship cranes, and they are liable to failure caused by the aggressive environment in which they operate [15,16].
Figure 7b shows the structural system for the dynamic ship crane reliability model on board the m/v “O” in POWERSIM symbols for states: in operation, out of order, control valve blocked, hydraulic motor vanes broken, safety valve blocked, hydraulic motor housing damaged, hydraulic pump gears damaged, and ship cranes out of operation.

3.4. Analysis of Obtained Results

For this research paper, the initial parameters were taken from the program package GALIOT [8] database for cranes on board the m/v “O”, for a period of 10 years and inserted into the mathematical model construed by using the Markov chains. The initial parameters obtained were used to construct the system dynamics simulation model in the POWERSIM simulation program. For simulation purposes, the time scenario of the system state during 120.000 working hours was used.
Figure 8a shows the reliability functions R_1, R_2, R_3, R_4, and R_5 for each state, unreliability functions F1, F2, F3, F4, and F5, and the condition to be satisfied at all times during monitoring R_1uk, R_2uk, R_3uk, R_4uk, and R_5uk, the sum of probable states must be 1 at all times.
Numerical results under Figure 8b show values for reliability R_1, R_2, R_3, R_4, and R_5 and unreliability F1, F2, F3, F4, and F5 of ship cranes over a time period of 120,000 working hours, obtained from the simulation.
The frequency of failures follows the bathtub curve, which is a base curve for the frequency of failures from which the length of initial failures can be determined, i.e., failures due to the system running in. To confirm that the reliability follows an exponential distribution, the frequency of failure must be constant, as shown in Figure 9.
The curve in Figure 9a shows failure intensity obtained empirically (calculated) and by simulation of the ship crane state over a time period of 120.000 working hours, from which it can be observed that in the first period of operation, i.e., 0–50.000 working hours, the failure intensity drops, and in the period from 50.000 to 120.000 working hours, it remains constant. Therefore, it can be concluded that fault intensities tend to be constant, are random, and cannot be avoided. They also appear with a certain probability, as can be seen from the simulation results and monitoring realities (failure history from the maintenance program package GALIOT software package database). The conclusion is that reliability follows an exponential distribution, as shown in Figure 9a. Failure frequency is directly related to the mean time of operation in individual states, as shown in Figure 10a,b.
From Figure 10a,b and numerical values obtained by simulation, where curves Tsr_4, Tsr_5, and Tsr_2 indicate that the system transitions into failure in states 4, 5, and 2 least frequently, i.e., state 4 “damage to hydraulic motor housing”, state 5 “damage to hydraulic pump gear” and state 2 “broken vanes on hydro motor” rarely occur. Tsr_1 and Tsr_3 show the frequency of transition to state 1 “control valve blocked” and state 3 “safety valve blocked”.
An increase in failure intensity (Figure 9) and mean time in operation (Figure 10) are affected by the tear and wear of components, aggressive environment—corrosion, fatigue, and the lifespan of materials.

3.5. Regression Analysis of Results Obtained from Simulation of a Given Scenario

To verify (confirm) the validity of the simulated mathematical model of ship crane states on board the m/v “O”, this paper used regression analysis as one of the most common statistical methods [4,17,18,19,20,21,22]. Values obtained from the simulation for reliability R_n, failure frequency Lam_n, and mean time Tsr_n depending on the given period of time (120.000 working hours) were statistically analyzed using regression analysis to confirm whether a strong correlation between the monitored parameters and times exists (see Figure 11a–c).
Besides the dependence of reliability, the mean time in operation, and failure frequency on time, it is necessary to show the dependence of the same reliabilities and mean times in operation with respect t to failure frequency for the parameters obtained by simulation (see Figure 12).
Figure 12 shows the dependence of two variables (mean time in operation Tsr_n in relation to failure frequency Lam_n), which is best described using an exponential dependence in the form of a curve, while the two other dependent variables (reliability of system R_n in relation to failure frequency Lam_n) are best shown using linear dependence. As failure frequency Lam_n increases, mean time in operation Tsr_n and system reliability R_n decrease, as shown in Figure 12. It indicates a strong correlation between the data obtained by simulation and shown by the determination coefficient R2, ranging from 0.966 to 0.99, and also indicates a 96–99% deviation described and interpreted by this regression model. It further indicates a large percentage of the representation. The conclusion is that the simulation model corresponds to real situations, and any state can be simulated with great accuracy.

4. Conclusions

The main objective of this research has been to develop a continuous simulation model for predicting ship crane failures in order to increase their reliability and reduce unplanned downtime during cargo loading and unloading operations. The key findings indicate that a model has been developed to predict failures in hydraulic motors and pumps of ship cranes. The Markov model achieves a precise estimation of the frequency of transitions between states, and the simulation model shows high reliability of cranes and a constant frequency of failures over a period of 120.000 operating hours. Regression analysis confirms the validity of the model and shows a strong correlation between the analyzed variables. These findings are significant as they provide insight into the possible development of a simulation model that can effectively predict ship crane failures. The results of the simulation and regression analysis show that the model can accurately estimate the transitions between failure and maintenance states, which supports the initial hypothesis. The practical implications are that the model can be used to plan the maintenance of ship cranes, reduce unplanned downtime, and promptly detect failures. The theoretical implications confirm the validity of the combined approach using FTA, the Markov model, system dynamics, and regression analysis to predict the reliability of complex technical systems. The research had some limitations, including the use of data specific to only a particular type of ship crane and a time period of 120.000 operating hours, which can generalize the results. The recommendation is to research other types of cranes under different operating conditions and environmental factors to determine when the frequency of failures transitions from exponential to a normal distribution, which would deepen our understanding and prediction of ship crane failures.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The author declares no conflicts of interest.

References

  1. Japan P&I Club. Cranes, Their Operation and Reasons for Failures; P&I Loss Prevention Bulletin; The Japan Ship Owners’ Mutual Protection & Indemnity Association Loss Prevention and Ship Inspection Department: Tokyo, Japan, 2015; Volume 36. [Google Scholar]
  2. Walton, P. Ship Cranes/Wires Problems and Their Prevention. Available online: https://www.nautinsthk.com/pluginAppObj/pluginAppObj_437_02/Ships-Cranes-Wires-Problems-Their-Prevention-2016.pdf (accessed on 15 April 2024).
  3. Mennis, E.; Agapios, P.; Ioannis, L.; Nikitas, N. Improving port container terminal efficiency with the use of Markov theory. Marit. Econ. Logist. 2008, 10, 243–257. [Google Scholar] [CrossRef]
  4. Jurjević, M. Simulacijsko Modeliranje Pouzdanosti Tehničkog Sustava Primjenom Sustavne Dinamike. Doctoral Dissertation, Faculty of Electrical Engineering, Mechanical Engineering and Naval Architecture, Split, Croatia, 2009. [Google Scholar]
  5. Szpytko, J.; Salgado Duarte, Y.A. Digital twins concept model for integrated maintenance: A case study for crane operation. J. Intell. Manuf. 2021, 32, 1863–1881. [Google Scholar] [CrossRef]
  6. Widarto, H.; Handani, D.W. Maintenance scheduling for diesel generator support system of container crane using system dynamic modeling. In Proceedings of the 2015 International Conference on Advanced Mechatronics, Intelligent Manufacture, and Industrial Automation (ICAMIMIA), Surabaya, Indonesia, 15–17 October 2015; pp. 188–193. [Google Scholar]
  7. Bello, A.E. Gantry Crane Maintenance: Proactive or Reactive; Port Technology International; Port of Miami Crane Management, Inc.: Miami, FL, USA, 2005; p. 27. [Google Scholar]
  8. A. P. GALIOT-Serving the World’s Shipping, Database for marine propulsion system MAN B&W 6L60MC on the m/v Exp. Galiot Maritime System is dedicated software package special designed for large and small Ship Management Companies, Split, Ship Inspection Report AP, Dubrovnik 2008, Croatia.
  9. Pudar, N. Use of amos in the ship maintenance optimisation process. In Proceedings of the Student International Scientific Conference, Batumi, Georgia, 17 May 2022. [Google Scholar]
  10. Bilić, B.; Jurjević, M.; Barle, J. Reliability Estimation of Technical System Using Markov Models and System Dynamics. Strojarstvo 2010, 52, 271–281. [Google Scholar]
  11. Jurjević, M.; Koboević, N.; Bošnjak, M. Modelling of dynamic reliability stages of a ship Modelling the Stages of Turbocharger Dynamic Reliability by Application of Exploitation Experience. Tech. Gaz. 2018, 25, 792–797. [Google Scholar]
  12. Lui, J.; Zio, E. System dynamic reliability assessment and failure prognostics. In Reliability Engineering & System Safety; Elsevier: Amsterdam, The Netherlands, 2017; Volume 160, pp. 21–36. [Google Scholar]
  13. Grabski, F. Semi-Markov Processes: Applications in System Reliability and Maintenance; Elsevier: Oxford, UK, 2014. [Google Scholar]
  14. Bazovsky, I. Reliability Theory and Practice; Manufactured in the United States by RR Donnelly, Dover Publications, INC: Mineola, NY, USA, 2015. [Google Scholar]
  15. Munitić, A.; Ristov, P. Sistemska Dinamika; Pomorski Fakultet, Sveučilište u Splitu: Split, Croatia, 2009. [Google Scholar]
  16. Baliwangi, L.; Arima, H.; Ishida, K.; Artana, K.B. Simulation on System Operation and Maintenance Using System Dynamics. J. JIME 2007, 42, 129–133. [Google Scholar] [CrossRef] [PubMed]
  17. Štambuk, A.; Biljan-August, M. Regresijska i Korelacijska Analiza; Sveučilišta u Rijeci, Ekonomski Fakultet u Rijeci: Rijeka, Croatia, 2009; ISBN 978-953-6148-86-8. [Google Scholar]
  18. Ebeling, C.E. An Introduction to Reliability and Maintainability Engineering; Waveland Press: Dayton, OH, USA, 2010. [Google Scholar]
  19. Lovrić, J. Osnove Brodske Terotehnologije; Pomorski Fakultet u Dubrovniku, Sveučilište u Splitu: Split, Croatia, 1989. [Google Scholar]
  20. Jakšić, S.; Erjavec, N.; Čeh Časni, A. Metode Primijenjene Matematičke i Statističke Analize; Metoda Primijenjene Matematičke i Statističke Analize, Sveučilište u Zagrebu-Ekonomski Fakultet: Zagreb, Croatia, 2020. [Google Scholar]
  21. Stazic, L.; Komar, I.; Račić, N. Evaluation Methodology for Ship’s Planned Maintenance System Database. Trans. Marit. Sci. 2017, 6, 109–116. [Google Scholar] [CrossRef]
  22. Lindqvist, B.H.; Doksum, A.K. Mathematical and Statistical Methods in Reliability; World Scientific: London, UK, 2003; Volume 7. [Google Scholar]
Figure 1. The sequence of methods and analyses used in the paper.
Figure 1. The sequence of methods and analyses used in the paper.
Jmse 12 01218 g001
Figure 2. Fault tree for the ship crane.
Figure 2. Fault tree for the ship crane.
Jmse 12 01218 g002
Figure 3. Fault tree for the ship crane hydraulic motor.
Figure 3. Fault tree for the ship crane hydraulic motor.
Jmse 12 01218 g003
Figure 4. Fault tree for the ship crane hydraulic pump.
Figure 4. Fault tree for the ship crane hydraulic pump.
Jmse 12 01218 g004
Figure 5. The fault tree with total faults that put the crane out of order.
Figure 5. The fault tree with total faults that put the crane out of order.
Jmse 12 01218 g005
Figure 7. (a) System dynamics qualitative structural model of ship crane reliability on m/v “O”; (b) structural system dynamics ship crane reliability model on m/v “O”.
Figure 7. (a) System dynamics qualitative structural model of ship crane reliability on m/v “O”; (b) structural system dynamics ship crane reliability model on m/v “O”.
Jmse 12 01218 g007
Figure 8. (a) Calculated values using the program; (b) numerical values for reliability.
Figure 8. (a) Calculated values using the program; (b) numerical values for reliability.
Jmse 12 01218 g008
Figure 9. (a) Failure intensity functions over a time period of 120.000 working hours; (b) numerical values for failure intensity.
Figure 9. (a) Failure intensity functions over a time period of 120.000 working hours; (b) numerical values for failure intensity.
Jmse 12 01218 g009
Figure 10. (a) The mean time in operation functions until reaching state 1, 2, 3, 4, and 5; (b) numerical values for mean time in operation.
Figure 10. (a) The mean time in operation functions until reaching state 1, 2, 3, 4, and 5; (b) numerical values for mean time in operation.
Jmse 12 01218 g010
Figure 11. Regression analysis of reliability and failure frequency on ship cranes on board the m/v “O”, obtained by simulation; (a) change of system reliability R_n (n = 1, 2, 3, 4, 5) as a function of time t; (b) change of transition frequency Lam_n (n = 1, 2, 3, 4, 5) as a function of time t; (c) mean time in operation for “n” states Tsr_n (n = 1, 2, 3, 4, 5) as a function of time t.
Figure 11. Regression analysis of reliability and failure frequency on ship cranes on board the m/v “O”, obtained by simulation; (a) change of system reliability R_n (n = 1, 2, 3, 4, 5) as a function of time t; (b) change of transition frequency Lam_n (n = 1, 2, 3, 4, 5) as a function of time t; (c) mean time in operation for “n” states Tsr_n (n = 1, 2, 3, 4, 5) as a function of time t.
Jmse 12 01218 g011
Figure 12. Correlation between reliability R_n and the mean time in operation Tsr_n (n = 1, 2, 3, 4, 5) shown on the y-axis in relation to the failure frequency shown on the x-axis.
Figure 12. Correlation between reliability R_n and the mean time in operation Tsr_n (n = 1, 2, 3, 4, 5) shown on the y-axis in relation to the failure frequency shown on the x-axis.
Jmse 12 01218 g012
Table 1. Frequency of transitions from state to state in case of failure.
Table 1. Frequency of transitions from state to state in case of failure.
λ01μ10λ02μ20λ03μ30λ04μ40
0.00007530.01470590.00002510.00666670.00005020.0160.00001250.02
λ05μ50α16α26α36α46α56μ60
0.00001250.003330.01470590.00666670.0160.020.03330.000009
P0P1P2P3P4P5P6Puk
0.10.00020.00010.000160.000030.000190.91.0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jurjević, M. Simulation Model of Hydraulic System States for Ship Cranes. J. Mar. Sci. Eng. 2024, 12, 1218. https://doi.org/10.3390/jmse12071218

AMA Style

Jurjević M. Simulation Model of Hydraulic System States for Ship Cranes. Journal of Marine Science and Engineering. 2024; 12(7):1218. https://doi.org/10.3390/jmse12071218

Chicago/Turabian Style

Jurjević, Mate. 2024. "Simulation Model of Hydraulic System States for Ship Cranes" Journal of Marine Science and Engineering 12, no. 7: 1218. https://doi.org/10.3390/jmse12071218

APA Style

Jurjević, M. (2024). Simulation Model of Hydraulic System States for Ship Cranes. Journal of Marine Science and Engineering, 12(7), 1218. https://doi.org/10.3390/jmse12071218

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop