Next Article in Journal
Prediction of Oil Saturation during Water and Gas Injection Using Controllable Convolutional Long Short-Term Memory
Next Article in Special Issue
Optimization of the Delivery Time within the Distribution Network, Taking into Account Fuel Consumption and the Level of Carbon Dioxide Emissions into the Atmosphere
Previous Article in Journal
Methods of Assessing the Effectiveness of Filter Elements in Power Electronics
Previous Article in Special Issue
Compressed Air Energy Storage Capacity Configuration and Economic Evaluation Considering the Uncertainty of Wind Energy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Light Utility Vehicle Readiness in Military Transportation Systems Using Markov and Semi-Markov Processes

by
Mateusz Oszczypała
,
Jarosław Ziółkowski
* and
Jerzy Małachowski
Institute of Mechanics and Computational Engineering, Faculty of Mechanical Engineering, Military University of Technology, gen. Sylwestra Kaliskiego Street 2, 00-908 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Energies 2022, 15(14), 5062; https://doi.org/10.3390/en15145062
Submission received: 17 June 2022 / Revised: 7 July 2022 / Accepted: 7 July 2022 / Published: 11 July 2022

Abstract

:
This paper presents the issues of modeling the operation process of light utility vehicles operating in military transport systems. The required condition for the effective operation of the system is to maintain the means of transport at the appropriate level of technical readiness. For this purpose, it is necessary to equip the technical system with appropriate resources enabling the efficient implementation of fuel refilling, maintenance and repair processes. Each failure of the means of transport causes a significant reduction in transport capacity, which then results in the inability to perform the planned tasks. Quality control and vehicle operation process management require advanced mathematical methods and tools. Three indicators have been proposed as quantitative characteristics for assessing and optimizing the availability of military vehicles: functional readiness, technical efficiency and airworthiness. To determine their value, a stochastic exploitation model was developed based on the application of the theory of Markov processes. Based on the collected empirical data, a nine-state phase space of the studied process was identified. Operating states were distinguished relating to the implementation of the transport task, refueling, parking in the garage, as well as maintenance and repairs. As part of the considerations for the continuous time, verification of the distributions of time characteristics led to the development of a semi-Markov model. The ergodic probabilities calculated based on the conditional probability matrix of interstate transitions and the expected values of the time spent in the states were used to determine the indicators of functional availability, efficiency and technical suitability. In order to determine the possibility of optimizing the process, a sensitivity analysis was performed. Reducing the amount of time the vehicles must wait for repair by about 50% can improve the values of the indexes from 0.91 to 0.95.

1. Introduction

Military transport systems aim to ensure that transport capacity corresponds to the existing needs. The determinants of the level of transport needs include combat operations, training plans and the current activity of military units. The fleet of reliable vehicles is one of the main factors determining the high quality and timely implementation of processes in modern transport systems. The condition for the effective operation of the system is to maintain the means of transport in a state of technical efficiency and be ready to perform tasks [1,2]. The fulfilment of this condition is possible thanks to the organization of operating systems with appropriate technical resources to carry out processes of diagnosis, servicing and repair of vehicles. Along with the increase in the intensity of the use of means of transport, the demand for fuel and other consumables increases significantly, especially in the case of military vehicles traveling outside the area of public roads [3,4,5].
Military exploitation systems are largely based on a plan-preventive maintenance strategy, aiming to maximize the technical availability of facilities [6,7,8,9]. This strategy assumes the implementation of maintenance with a specific labor intensity, in accordance with the required scope. The guidelines for maintenance activities, time intervals and the size of the service life between maintenance are defined based on the manufacturer’s technical specifications and the knowledge and experience of specialists dealing with planning and standardization of operation at individual management levels. The main disadvantages of this strategy are its cost-intensive nature and low flexibility.
The operation process covers activities related to a technical object from the moment of its production to its liquidation. During this period, there are essentially two overlapping sub-processes usage and maintenance. The rational use of machines and devices allows for extending the intervals between subsequent maintenance services, recreating the technical service life of the facility, and also reducing the current operating costs. The model of the operation process should, on the one hand, reflect the basic technical characteristics of the object consistent with the modeling objective, and on the other hand, it should be used for the rational forecasting of use and maintenance [10,11,12,13].
Markov’s theory has found applications in many fields of science and technology. There are many scientific studies in the literature on the use of Markov processes for reliability modeling [14,15], the operation of objects [16,17,18,19] and technical systems [20,21]. Depending on the case study and the purpose of modeling, the authors constructed models based on a diverse number of operational states in the phase space. The least complicated Markov chain model was presented in [22] and applied to simulate and optimize energy savings for machines operating in production systems. Models with three states were constructed to analyze and assess the technical availability of buses [23], special vehicles [24], operational readiness of wind turbine elements [17], working time of production machines [25], the time interval of preventive maintenance [26] and the reliability of technical facilities [27]. In [28], the authors developed a five-state semi-Markov model with the Weibull distribution of the residence times of four types of buses in the states of the renewal process. This model allowed for profit optimization per unit of time and availability depending on the duration of preventive maintenance.
The comparison of the values obtained by the six-state Markov and semi-Markov models, which the authors of the publication [29] developed for production machines, indicated that the unverified assumption of the exponential distribution of the time the object stays in states may lead to significant errors in the calculation of the readiness indices. In the presented case study, the difference between the results of the semi-Markov model and the erroneous Markov model was as much as 0.40, which is less than half of the actual value of the readiness index.
Markov models with 9- and 16-state phase spaces are presented in [8,30] with much more elaborate models. The increased number of states allows a detailed analysis of the process and factors affecting the technical readiness of the facility. The multi-state models presented in [31] accurately reflected the stochastic nature of the electric vehicle driving cycle during their use in urban areas, outside the city, on the motorway and during road congestion. Table 1 summarizes the literature review containing the latest publications in the field of modeling the exploitation process with the use of Markov theory.
In this publication, the authors addressed the issue related to the operation of light utility vehicles operating in military transport systems. The research sample consisted of 19 Honker vehicles for which detailed data were collected during the three-year research period. The operating system is focused on maintaining the high reliability of vehicles through an appropriately planned and implemented maintenance strategy. The plan-preventive system each time assumes the scope of maintenance works after a specified amount of work (mileage) or time elapsed. Unfortunately, the records of operation are still kept in the form of traditional documentation registered on an ongoing basis by direct users (drivers) and in relation to inspections and repairs by service and repair workshops. Source documents that create departure orders, technical service cards and operation plans were used to prepare detailed databases individually for each facility.
The current review of the literature allows the authors to state there are no studies on the analysis and evaluation of the operation of heavy goods vehicles with the use of the Markov theory. It was a premise for conducting scientific research and developing an exploitation model for the aforementioned group of military vehicles. In addition, Honker vehicles have a significant share in the structure of the military transport fleet of the Polish Armed Forces. Carrying out the modeling of the operation process based on the application of the Markov theory requires a thorough understanding of the examined process and enables the analysis and assessment of the basic operational indicators of the studied object.
This article presents the original methodology for creating a stochastic model of Honker vehicles based on the Markov theory. An algorithm for creating a mathematical model was developed, the practical usefulness of which was verified on the actual operation process of the said sample of vehicles. An event model covering the nine-state phase space of the process was developed. The determined values of ergodic probabilities for individual operational states constituted the basis for calculating the values of functional availability, efficiency and technical suitability indicators. From the point of view of the transport capacity of the entire system as well as economic and technical conditions, the personnel managing the vehicle operation process aim to maximize the presented readiness and/or reliability measures. The proposed methodology makes it possible to indicate possible components influencing the improvement of exploitation indicators. The performed sensitivity analysis of the model allows for the examination of the impact of improving the organization of the repair subsystem, consisting of the shortening/elimination of waiting time for spare parts, on the readiness indicators of the tested sample of vehicles.
The nine-state model adds to the current state of the literature both in terms of the subject of the study and the application of sensitivity analysis to identify opportunities for improving process efficiency. No previous scientific studies have addressed the issue of modeling the operation of light utility vehicles (trucks) using Markov and semi-Markov process theory. In addition, a completely novel way of analyzing the sensitivity of the semi-Markov model in terms of the dependence of the values of ergodic probabilities on the values of expected dwell times was proposed.
The article has been divided into the following main chapters. The introduction reviews the current state of knowledge on the application of Markov theory to the exploitation process. Section 2 describes the methodology of creating event-based models of the operation process using the Markov theory. Section 3 provides a statistical analysis of the source data constituting the basis for the development of the model. In Section 4, the semi-Markov model is described, and the process research results are presented together with the model sensitivity analysis. The values of ergodic probabilities of the semi-Markov model were confronted with the standard Markov model. Finally, the Section 5 includes the conclusions from the conducted research.

2. Methods

The actual operation processes are a composition of deterministic and random sub-processes. Random components are usually interpreted as stochastic processes X(t) reflecting changes in the operational states of the tested object in discrete or continuous time. In the processes of exploitation at a random moment t, the object is only in one of the states identified in the phase space S = X(t). This assumption requires precise determination of all possible operational states in which vehicles may be in the course of the operation process. The stochastic processes fulfilling the Markov property are essential in terms of applicability. According to Markov’s theory, the conditional probabilities of reaching the future states X(tz+1) result only from the current state X(tz) [32]. Mathematically, this property is consistent with Equation (1) [33,34,35]:
P { X ( t z ) = x z | X ( t z 1 ) = x z 1 , X ( t z 2 ) = x z 2 , , X ( t 0 ) = x 0 } = P { X ( t z ) = x z | X ( t z 1 ) = x z 1 } .
The literature is dominated by the division of Markov processes with regard to time and state space, which distinguishes four types of processes, i.e., processes:
  • Discrete in time and discrete in states;
  • Continuous in time and discrete in states;
  • Discrete in time and continuous in states;
  • Continuous in time and continuous in states.
In operation, the most frequently used models are based on discrete processes in states, developed for both discrete and continuous time [8,16,24,28,29].

2.1. Markov Chains

The Markov chain assumes the discretization of time into specific Δt intervals, the width of which depends on the characteristics of the process, measurement technique, and, above all, the adopted modeling objective [36]. With the increase in the dynamics of the process course, the possibility, and, at the same time, the necessity to perform frequent measurements and the focus on creating a very accurate model, the values of time Δt decrease. The dynamics of changes in the operational states identified in the phase space for the operation of light utility vehicles determine the assumption of the duration of moments Δt for the time interval equal to 1 min. Increasing this value could cause the undesirable omission of registering the occurrence of conditions, which are usually short-lived. At the same time, reducing the time intervals is practically impossible due to recording the course of operation processes of vehicles operating in military transport systems.
Constructing the Markov Chain Model based on an empirical process flow requires acquiring data on interstate transitions. For this purpose, it is reasonable to create a matrix of the number of interstate transitions according to the Formula (2):
N = [ n 11 n 12 n 1 ( k 1 ) n 1 k n 21 n 22 n 2 ( k 1 ) n 2 k n ( k 1 ) 1 n ( k 1 ) 2 n ( k 1 ) ( k 1 ) n ( k 1 ) k n k 1 n k 2 n k ( k 1 ) n k k ] .
The values of the N matrix correspond to the total number of observed interstate transitions in the analyzed period of the process implementation, where nij is the transition from the i state to the j state.
For a homogenous Markov chain, the conditional probability pij of transition from state i to state j in one step is the same for every moment t. The homogeneity of the process indicates the invariability of the rules affecting the state changes at any time of its implementation. If the analyzed realizations of the process are included in the same phase of operation, then the course of the process should be homogenous. The probability values of the conditional interstate transitions are presented by means of the stochastic matrix P, according to the Formula (3) [33,37,38]:
P = [ p 11 p 12 p 1 ( k 1 ) p 1 k p 21 p 22 p 2 ( k 1 ) p 2 k p ( k 1 ) 1 p ( k 1 ) 2 p ( k 1 ) ( k 1 ) p ( k 1 ) k p k 1 p k 2 p k ( k 1 ) p k k ] ,
provided that the following formula [37] is fulfilled:
j = 1 k p i j = 1 .
The probabilities of the interstate transitions of the P stochastic matrix can be obtained using the values of the N interstate number matrix by estimation [31,35,39] according to the relationship:
p i j = n i j j = 1 k n i j ,
where the standard error of the conditional probability estimation [40,41] is calculated according to the formula:
S E ( p i j ) = p i j ( 1 p i j ) j = 1 k n i j .
The values of ergodic probabilities πj are calculated by solving the following matrix Equation (7) [37]:
( P T I ) Π = 0 ,
assuming that the normalization condition is met, according to the formula:
j = 1 n π j = 1 .

2.2. Markov and Semi-Markov Processes

Discrete Markov models in states and continuous in time allow for the analysis of the operation process under constant supervision and monitoring of the course of changes in operational states. This approach excludes the influence of the size of the time moments Δt on the values of the instantaneous and ergodic characteristics of the process. In principle, the change of state from Si to Sj can take place at any time during the process. It is also assumed that only one state change can occur at any time t. The Markov process assumes the presence of exponential distributions of the characteristics of individual states describing the analyzed process [42,43,44].
The transition intensity matrix Λ is the quantitative characteristic of the Markov process:
Λ = [ λ 11 λ 12 λ 1 ( k 1 ) λ 1 k λ 21 λ 22 λ 2 ( k 1 ) λ 2 k λ ( k 1 ) 1 λ ( k 1 ) 2 λ ( k 1 ) ( k 1 ) λ ( k 1 ) k λ k 1 λ k 2 λ k ( k 1 ) λ k k ] ,
whose elements meet the following dependencies:
i , j , i j , λ i j = d d t p i j = lim Δ t 0 p i j ( t + Δ t ) p i j ( t ) Δ t ,
i , λ i i = d d t p i i = lim Δ t 0 1 ( p i i ( t + Δ t ) p i i ( t ) ) Δ t .
All elements of the matrix on the main diagonal have non-positive values, while all other elements are non-negative and the sum of all elements for each row of the transition intensity matrix is equal to 0.
In the operation processes of objects, the estimators of the values of the elements λij of the interstate intensity matrix are the reciprocal of the average residence times in the Si state before the transition to the Sj state. On the other hand, the value of λii is taken as the reciprocal of the sum of the remaining elements in the i row. The presented description of the estimation of the intensity matrix elements can be written using the relationship:
λ i j = 1 T i j ,
λ i i = j i λ i j ,
where Tij is the average transition time from state Si to state Sj.
The condition of exponential distributions of the stochastic process characteristics narrows the spectrum of applications of the Markov model in continuous time. Some characteristics of the operation process may not meet it. From the point of view of the reliability of the mapping, this excludes the possibility of applying the Markov model [29,45]. The generalization of the Markov process is the semi-Markov process, which does not require the fulfilment of the condition of exponential distributions of interstate transition times. In the semi-Markov process, the durations of states are independent random variables with any distribution function [27,43,46].
The basic characteristic of the semi-Markov process is the matrix of the renewal kernel Q(t), the elements of which are the products of the probability pij and the transition between the states Si and Sj and the distribution function of the conditional duration distribution of the state Si before the transition to Sj, according to the equation [30,47]:
Q ( t ) = [ 0 Q 12 ( t ) Q 1 ( k 1 ) ( t ) Q 1 k ( t ) Q 21 ( t ) 0 Q 2 ( k 1 ) ( t ) Q 2 k ( t ) Q ( k 1 ) 1 ( t ) Q ( k 1 ) 2 ( t ) 0 Q ( k 1 ) k ( t ) Q k 1 ( t ) Q k 2 ( t ) Q k ( k 1 ) ( t ) 0 ] ,
wherein:
Q i j ( t ) = p i j F i j ( t ) ,
where pij is the probability of transition from the Si state to Sj, and Fij(t) is the distribution function of the residence time in state Si before the transition to state Sj.
An embedded Markov chain is formulated for a semi-Markov process in continuous time, which represents changes in the process state without taking into account the residence times in individual states.
The embedded Markov chain assumes the possibility of transition from the Si state to Sj, assuming that ij. The interstate transition probability matrix for such a chain can have non-zero elements only outside the main diagonal, which can be written by the formula:
P = [ 0 p 12 p 1 ( k 1 ) p 1 k p 21 0 p 2 ( k 1 ) p 2 k p ( k 1 ) 1 p ( k 1 ) 2 0 p ( k 1 ) k p k 1 p k 2 p k ( k 1 ) 0 ] .
If the embedded Markov chain is ergodic and there are expected E(Ti) values of the times in individual states, then the ergodic values of the probabilities πj for the semi-Markov process can be determined using the following dependencies:
π j = p j E ( T j ) i = 1 k p i E ( T i ) ,
E ( T j ) = i = 1 k p i j T i j ,
where pi is the ergodic probability of the embedded Markov chain for state Si, and Tij is the average transition time from state Si to state Sj.
Figure 1 shows a block diagram of creating a model of the exploitation process based on the theory of Markov and semi-Markov. The first stages of stochastic modeling are the identification of the technical object, the selection of a statistical sample, and the collection of empirical data in the form of databases. Then a choice is made between discrete-time and continuous-time models. In the case of discrete-time, a model is constructed based on a Markov chain. When analyzing a process in continuous time, verification of the exponential distribution of time characteristics, a condition for the possibility of using a Markov model, is carried out. Two non-parametric consistency tests are proposed as statistical verification tools: χ2 and Kolmogorov, depending on the size of the research samples. For empirical data containing statistical samples less than 80, the Kolmogorov test is recommended, while otherwise, there is no contraindication to using the χ2 test [48,49]. For processes satisfying the condition of the exponential distribution, Markov models are used, while otherwise, a semi-Markov model is an appropriate solution. The calculation of the values of operating indicators is carried out on the basis of the ergodic probabilities of the corresponding model.

2.3. Functional Readiness, Technical Efficiency and Technical Suitability

Functional readiness is usually understood as the ability of a technical system or object to undertake and perform tasks consistent with its intended use in the required time [50,51]. The availability index Kr reflects the quantitative characteristic of functional availability, which for the Markov model is expressed as the sum of the probabilities of ergodic operational states kSr, in which the object can start the task or is in the process of its implementation. The mathematical notation of this relationship is presented by the formula:
K r = k S r π k ,
where πk is the ergodic probability of the desired (from the standpoint of readiness) set of operational states kSr.
In functional readiness, the vehicle is technically fully operational, i.e., it has an adequate supply of fuel and consumables and is not being serviced. The functional availability indicator shows the share of time that can be allocated to the implementation of tasks by the technical object during its operation period.
The concept of technical efficiency refers to a wider set of operational states than in the case of functional availability. The technically efficient facility has a technical resource to perform the tasks. However, it may require refueling or short-term maintenance before or after use. For the Markov model, the technical efficiency coefficient Ke is the sum of the ergodic probabilities of operational states lSe in which the technical object has a technical service life and does not require periodic maintenance, which is shown in the relationship:
K e = l S e π l ,
where πl is the ergodic probability of the desired set of operational states lSe.
Technical suitability expresses the condition of a technical object, in which it is not damaged, or there is no need to repair it. The technical suitability condition is an extension of the technical efficiency by the time needed to perform periodic maintenance in order to restore its technical life. The technical suitability index Ks for the Markov model is therefore, the sum of the probabilities of ergodic operational states mSs, in which the object is not damaged. This is expressed in the equation:
K s = m S s π m ,
where πm is the ergodic probability of the set of operational states mSs.
There is the following relationship between the sets defining the states of functional readiness, technical efficiency and technical suitability:
S r S e S s .

3. Object of Analysis

3.1. Case Study

The case study of the conducted research is a trial of Honker light utility vehicles, constituting a fleet of vehicles of the military unit transport system. The analyzed technical facilities perform tasks related to the transport of people and small loads weighing up to 1000 kg.
Determining the phase space of the studied process requires the reproduction of a detailed phase trajectory for each object in the three-year research period. In the case of the operation process, one should also take into account the conditions resulting from the operation organization system as well as standards and guidelines for the maintenance subsystem of technical facilities. In the case of military transport systems, the implementation of the vehicle operation process depends on both instructions and procedures developed by the operators of military equipment as part of the adopted operational strategy. As a result of the analysis of the process of light utility vehicles carried out by the authors of this study, a nine-state phase space was identified, for which the possible interstate transitions were specified in Table 2.
The distinguished operating states are mutually disjoint subsets, which means that at any time t the object may be in only one of them. For example, the vehicle is in the S3 garage state while waiting to be used. Before starting the task, it is necessary for the direct user to check the correct operation of the systems and mechanisms that determine the safe use of the vehicle. These activities correspond to state S4. State S1 means completion of the transportation task, upon completion of which the object should be serviced after use (S5). Its purpose is to re-check its operational condition, including the removal of minor defects and cleaning the body. The implementation of the S2 refueling state may occur before, during, or after the task commencement, depending on the identified needs. State S6 corresponds to the performance of periodic vehicle maintenance in accordance with the operational strategy as well as instructions and guidelines. S7, S8 and S9 states are undesirable from the standby point of view, as they symbolize damage to the object and the need for repair.
The database was developed on the basis of operational documentation covering a three-year research period. Its fragment is presented in Table 3. In the following rows of the database, the transitions between the various operational states were identified, along with detailed records in the following system: transition number, date, hour and minute. The accuracy of the measurement in the order of 1 min is determined by the occurrence of short-term states and is necessary to verify the correctness of the time balance over the entire research period.

3.2. Statistical Analysis of Data

For the verification of the empirical distribution with the theoretical distribution, non-parametric consistency tests are commonly used [52,53,54]. The choice of test depends on the type of distribution being verified, the random variable and the sample size. The most common test χ2 requires a sample size of at least 80 realizations of a random variable. For a small number of interstate transitions in the operation process, this condition may not be met. In this case, the verification of the assumptions of the Markov process is performed using the Kolmogorov test.
The condition for the exponential distribution of transition times between states was verified with the use of the Kolmogorov conformance test (Table 4). As the H0 hypothesis, it was assumed that the distribution of Tij times for individual allowed interstate transitions is consistent with the exponential distribution. The alternative hypothesis (H1) contradicts this assumption. The Uy statistic of the Kolmogorov test has the form (23):
U y = m a x 1 i y ( d y , d y + ) ,
wherein:
d y = m a x 1 i y | F ( x i ) i 1 y | ,
d y + = m a x 1 i y | i y F ( x i ) | ,
where F(xi) is the cumulative value of the theoretical distribution according to the H0 hypothesis, and y is the sample size.
For the Uy statistic values within the range of critical values <D(y;α), 1> for the significance level α, the H0 hypothesis is rejected and H1 is accepted. Otherwise, there is no reason to reject the H0 hypothesis. For the significance level α = 0.05, the value D(y;α) can be approximated using the formula:
D ( 0.05 ; y ) = 1.358 y + 0.12 + 0.11 y .
At the adopted significance level of α = 0.05, only two time characteristics, T16 and T71, achieved the values of the Kolmogorov test statistics, which are not included in the critical intervals. This means there are no grounds to reject the null hypothesis assuming the exponential distribution of the times of stay in individual states. Thus, for the remaining time characteristics, the null hypothesis was rejected and an alternative was adopted. The results of the statistical test for all variables are presented in Table 4. Figure 2 presents the frequency graphs for two exemplary times T16 and T34 and the course of the exponential distribution with the parameter estimated on the basis of the reciprocal of the mean values from both statistical samples. The number of histogram intervals was determined using the square root method.
The exploitation process does not fulfill the condition of exponential distributions and, therefore, cannot be considered a Markov process in continuous time. In this case, semi-Markov processes should be used for modeling.

4. Results and Discussions

4.1. Semi-Markov Model (SMM)

The embedded Markov chain in the semi-Markov process, in accordance with the adopted assumptions, does not allow for the possibility of returns (transition from the Si state to the Si state). This assumption adopts that every change of the process state is recorded, while its absence means that the object is in the Si state before transitioning to the next Sj state for a period of time equal to Tij [26,28,29,30]. The matrix (27) shows the empirical numbers of transitions between particular exploitation states as a result of observation of the process. On the other hand, matrix (28) represents the probabilities of transitions between states estimated on the basis of the matrix of the number of transitions.
N = [ 0 798 0 0 3340 32 0 32 0 396 0 366 0 695 0 0 0 0 0 365 0 4063 0 76 0 0 0 3769 294 0 0 0 0 0 0 0 0 0 3988 0 0 0 0 0 70 11 0 84 0 23 0 0 0 0 27 0 67 0 0 10 0 0 44 0 0 0 0 0 0 63 0 39 0 0 0 0 0 0 85 70 0 ] ,
P = [ 0 0.189910 0 0 0.794860 0.007615 0 0.007615 0 0.271791 0 0.251201 0 0.477008 0 0 0 0 0 0.081039 0 0.902087 0 0.016874 0 0 0 0.927640 0.072360 0 0 0 0 0 0 0 0 0 0.982750 0 0 0 0 0 0.017250 0.093220 0 0.711864 0 0.194915 0 0 0 0 0.182432 0 0.452703 0 0 0.067568 0 0 0.297297 0 0 0 0 0 0 0.617647 0 0.382353 0 0 0 0 0 0 0.548387 0.451613 0 ] .
Table 5 summarizes the values of the standard error for the conditional probabilities (28) estimated on the basis of the matrix of the number of interstate transitions (27). According to Formula (6), with the increase in the number of transitions from the Si state, the standard error SE(pij) decreases. For this reason, the values of the standard error are higher for operational states in which the technical object is relatively rare. Nevertheless, for all conditional probabilities, the SE(pij) did not exceed the value of 0.05 [55,56]. The result at this level is considered acceptable.
Figure 3 presents a graph illustrating possible transitions between states during the implementation of the operation process. According to the assumptions made for the embedded Markov chain, the fact that the object remains in the same state is not treated as a SiSi transition. For this reason, the SMM model graph does not have connections coming from the Si state and going directly to the Si state. This means a lack of return to the same state, which is commonly used in modeling the operation processes.
Assuming that at t = 0, the vehicle is technically efficient and awaits the appearance of the task, it is possible to determine the instantaneous probabilities of the embedded Markov chain. This assumption reflects the initiation of the operation process for a vehicle included in the transport system. The instantaneous probabilities are the matrix product of the initial distribution matrix p0 and the n-th power of the conditional probability matrix P, where n corresponds to the number of transitions between states. The development of the dependence of the instantaneous probabilities on the number of interstate transitions allows us to determine the period after which their values stabilize at a certain level and the stochastic process reaches the equilibrium state.
Figure 4 and Figure 5 show changes in the value of the instantaneous probabilities of the embedded Markov chain in the semi-Markov process. In the range from t = 0 to t = 30, there are fluctuations with a large amplitude of changes in the values, which decrease with time and after overcoming about 50 transitions, the probability values remain constant.

4.2. Ergodic Probabilities of SMM

The ergodic probabilities of the semi-Markov process are determined on the basis of the probabilities of the embedded Markov chain and the values of the expected time in individual states of the process. For the embedded Markov chain, the ergodic probabilities are computed using the matrix equation:
( P T I ) P j = [ 1 p 21 0 p 41 0 p 61 p 71 0 0 p 12 1 p 32 p 42 0 0 0 0 0 0 p 23 1 0 p 53 p 63 p 73 0 0 0 0 p 34 1 0 0 0 0 0 p 15 p 25 0 0 1 p 65 0 0 0 p 16 0 p 36 0 0 1 p 76 0 0 0 0 0 0 0 0 1 p 87 p 97 p 18 0 0 0 0 0 0 1 p 98 0 0 0 0 p 59 0 p 79 p 89 1 ] [ p 1 p 2 p 3 p 4 p 5 p 6 p 7 p 8 p 9 ] = [ 0 0 0 0 0 0 0 0 0 ] ,
which can be written as a system of Equation (30) together with the condition (31) of system normalization:
{ p 1 + 0.271791 p 2 + 0.927640 p 4 + 0.093220 p 6 + 0.182432 p 7 = 0 0.189910 p 1 p 2 + 0.081039 p 3 + 0.072360 p 4 = 0 0.251201 p 2 p 3 + 0.982750 p 5 + 0.711864 p 6 + 0.452703 p 7 = 0 0.902087 p 3 p 4 = 0 0.794860 p 1 + 0.477008 p 2 p 5 + 0.194915 p 6 = 0 0.007615 p 1 + 0.016874 p 3 p 6 + 0.067568 p 7 = 0 p 7 + 0.617647 p 8 + 0.548387 p 9 = 0 0.007615 p 1 p 8 + 0.451613 p 9 = 0 0.017250 p 5 + 0.297297 p 7 + 0.382353 p 8 p 9 = 0 ,
j = 1 9 p j = 1 .
After solving the system of equations, the values of ergodic probabilities pj of the inserted Markov chain were obtained. The expected E(Tj) times of staying in individual operating states were determined as average values based on empirical data. The values of pj and E(Tj) were used to calculate the values of ergodic probabilities in the semi-Markov model. The results of the conducted analyses are summarized in Table 6.
The garage state S3 had the highest value of the ergodic probability of over 77%. On the basis of the calculated values of probabilities, it can be concluded that the trucks and passenger cars stay together in other states for approximately 13% of the time during the entire three-year research period. The S9 state also obtained a significant level of ergodic probability, which indicates that the vehicles remain in a state of waiting for repair for almost 9% of their operational time. The remaining operational states obtained probability values below 1% and do not have a significant impact on the availability rates.

4.3. Calculations of Indicators in SMM

In the nine-state exploitation model, state subsets were distinguished corresponding to the functional readiness Sr, technical efficiency Se and technical suitability Ss. The mathematical notation is presented by the Formulas (32)–(34):
S r = { S 1 , S 3 } ,
S e = { S 1 , S 2 , S 3 , S 4 , S 5 } ,
S s = { S 1 , S 2 , S 3 , S 4 , S 5 , S 6 } .
Functional readiness corresponds to the vehicle being in the task (S1) or garage (S3) state. Technical efficiency extends the set of these states with fuel refilling (S2) and the implementation of maintenance before starting the task (S4) and after its completion (S5). On the other hand, technical suitability also takes into account the implementation of periodic maintenance (S6), the purpose of which is to restore the technical service life. The graphical diagram of the division of the operating conditions set into individual subsets is shown in Figure 6.
For such defined subsets, the values of functional availability, technical efficiency and technical suitability indicators were calculated. The results are presented in Table 7. The values of all three indicators are similar, which means a small share of the time of current and periodic maintenance and refueling in the entire test period.
The high probability value of the S9 ergodic state has the greatest impact on the reduction of the indicators of functional availability, technical efficiency and technical suitability of vehicles. This state corresponds to the situation when the vehicle has suffered a breakdown, is out of order and requires repair. Due to the lack of technical possibilities, it is not in the repair state (S7) and is not in the diagnosis phase (S8). The main reasons for such a situation are logistic delays related to the limited availability of spare parts and the lack of qualified technical personnel within the specified time.

4.4. Sensitivity Analysis of SMM

The calculated ergodic probability π9 of the S9 state has the strongest impact on the values of the Kr, Ke and Ks indicators. Due to this fact, the sensitivity analysis of the SMM model was carried out in order to investigate the impact of the presence of objects in the S9 state on the ergodic probabilities and technical readiness rates. The parameter against which the variability analysis was performed is the expected value of the duration of the S9 state. For the change of this parameter, expressed as a percentage and amounting to ΔE(T9), the values of the ergodic probabilities of the semi-Markov process can be determined using the relationship:
π j = p j E ( T j ) i = 1 8 p i E ( T i ) + p 9 E ( T 9 ) ( 1 Δ E ( T 9 ) )   for   j = { 1 , 2 , 3 , , 8 } ,
π 9 = p 9 E ( T 9 ) ( 1 Δ E ( T 9 ) ) i = 1 8 p i E ( T i ) + p 9 E ( T 9 ) ( 1 Δ E ( T 9 ) ) .
The undoubted advantage of the SMM model is the ability to perform a sensitivity analysis for a wide range of changes in the E(T9) parameter without the need to solve many matrix equations. Table 8 shows the results of the analysis for selected ΔE(T9) values. Figure 7 presents a graph of changes in ergodic probabilities for continuous changes in the value of ΔE(T9).
On the basis of the determined ergodic probabilities, the possibilities of improving the values of the readiness ratios were carried out at the assumed levels of reduction of the time spent in the state of incapacity and waiting for repair. The results are presented in Table 9 and Figure 8.
Reducing the waiting time for repairs by 50% would result in an increase in the values of the Kr, Ke and Ks indicators to the level of about 0.95, which at the same time means their improvement by over 0.04. The sensitivity analysis of the SMM model was carried out with regard to the impact of changes in the expected value of the time of staying in the S9 state.

4.5. Markov Model (MM)

Markov models, assuming exponential distributions of time characteristics, often use stochastic models to describe the operation processes of machines, devices and technical systems. For small deviations of empirical distributions of interstate times from the theoretical assumptions of the Markov theory, the results obtained for Markov and semi-Markov models may be similar or even negligibly small.
For the analyzed sample of means of transport, the non-parametric Kolmogorov test showed grounds for rejecting the null hypothesis about exponential distributions. However, the authors of the publication conducted an analysis of the applicability of the Markov model as a simplification of the generalized SMM presented in Section 4.1, Section 4.2, Section 4.3 and Section 4.4. The basic characteristic of the Markov model is the matrix of interstate transition intensity Λ. In the case study under consideration, the value of the Λ matrix presented in the matrix, which were estimated on the basis of the mean values of the transition times between the states according to the Equations (11) and (12).
Λ = [ 0.012136 0.000890 0 0 0.001593 0.007880 0 0.001773 0 0.214402 0.657021 0.229036 0 0.213583 0 0 0 0 0 0.000241 0.000650 0.000257 0 0.000152 0 0 0 0.199841 0.200000 0 0.399841 0 0 0 0 0 0 0 0.122711 0 0.270079 0 0 0 0.147368 0.003216 0 0.003819 0 0.003446 0.010481 0 0 0 0.007004 0 0.005194 0 0 0.010753 0.028343 0 0.005392 0 0 0 0 0 0 0.045652 0.062586 0.016934 0 0 0 0 0 0 0.000068 0.000088 0.000156 ] .
Ergodic probabilities of the Markov process for the entire set of operational states are calculated by solving the matrix Equation (38) together with the system normalization condition (39).
Π T Λ = [ π 1 π 2 π 3 π 4 π 5 π 6 π 7 π 8 π 9 ] T [ λ 11 λ 21 0 λ 41 0 λ 61 λ 71 0 0 λ 12 λ 22 λ 32 λ 42 0 0 0 0 0 0 λ 23 λ 33 0 λ 53 λ 63 λ 73 0 0 0 0 λ 34 λ 44 0 0 0 0 0 λ 15 λ 25 0 0 λ 55 λ 65 0 0 0 λ 16 0 λ 36 0 0 λ 66 λ 76 0 0 0 0 0 0 0 0 λ 77 λ 87 λ 97 λ 18 0 0 0 0 0 0 λ 88 λ 98 0 0 0 0 λ 59 0 λ 79 λ 89 λ 99 ] = [ 0 0 0 0 0 0 0 0 0 ] T ,
j = 1 9 π j = 1 .
Mathematica and MS Excel software were used for the presented calculations. The results for MM are summarized in Table 10 and Figure 9, comparing them with the values of ergodic probabilities obtained for SMM. Percentage differences between the models reached significant values, with the most similar probabilities occurring for the S2 state and differing by over 41.0%. The greatest discrepancies in the results occurred for the S6 state, for which the ergodic probability obtained in MM was as much as 1153.69% higher than in SMM. In addition, the S6, S7, S8 and S9 states, which are a subset of the technical inoperability and failure states, achieved positive differences. The use of MM to evaluate the operation process would result in a significant reduction of the values of the Kr, Ke, and Ks indices, inconsistent with the actual state.
Compared to SMM, the Markov model achieved a mean absolute percentage error (MAPE) of 351.92%. Therefore, the hypothesis about the possibility of using the Markov process to describe state changes by the analyzed technical objects in the considered nine-element set of operational states should be unambiguously rejected.

5. Conclusions

The theory of Markov processes was used to model the operation of light utility vehicles. The paper presents a semi-Markov model that allows for the assessment the technical readiness of means of transport operating in a real system. On the basis of the systems of Chapman–Kolmogorov equations and the values of the expected times of stay in operating states, the values of the ergodic probabilities of the process were determined. The main reason for the reduction of functional availability, technical efficiency and technical suitability indicators was the high probability value of the S9 state (Awaiting repair). This indicates the occurrence of significant delays in the implementation of repairs of damaged vehicles. Therefore, in eliminating the reasons for the presence of technical objects in the S9 state, the operating system does not have sufficient technical resources to carry out repairs without unnecessary time delay. Repair delays are mainly caused by the unavailability of spare parts and a shortage of qualified technical personnel. In order to maximize the technical readiness of vehicles, it is necessary to focus on improving the organization of the spare parts delivery system, which will significantly reduce the time spent by technical facilities in the S9 state.
For the current state of the operation system, the functional availability ratio Kr has reached the value of 0.907334, which should be interpreted as follows: for more than 90% of the duration of the operation process, vehicles in good technical condition await the appearance of a task or are in the process of its implementation. A slight difference between the technical efficiency index Ke and the functional readiness index Kr indicates that the process of refueling and servicing is carried out efficiently. The aforementioned processes constitute a set of activities preparing a technically efficient vehicle to perform the task, as well as control and check after its completion. On the other hand, the technical suitability index Ks amounting to 0.911343 means that for over 91% of the duration of the operation process, the vehicles are fit for use. Its high value indicates a well-thought-out and proper implementation of the exploitation strategy. The small value of the ergodic probability for the S6 state resulted in a slight difference between the technical efficiency index Ke and the technical suitability index Ks. The above-mentioned results show that the capabilities of the technical subsystem, which carries out the periodic maintenance process, are adjusted to the requirements of the plan-preventive strategy used.
The ergodic probability π2 of 0.000292 and the expected residence time in the S2 state of 4.60 (min) indicate an efficiently implemented refueling process. For the current level of intensity of use of vehicles, in which over 13.0% of the operation time is during the implementation of transport tasks, the technical system provides sufficient resources of diesel oil and appropriate distribution equipment.
Despite the assumed high level of readiness, efficiency and suitability indicators of the tested military vehicles, an attempt was made to optimize the process by determining the impact of a potential reduction in the duration of the S9 state, equipped with logistic delays occurring in the operation system. An analysis of the optimization possibilities was carried out on the basis of changes in the expected time of stay in the S9 state based on analytical relationships, which enabled the analysis based on a continuous reduction of changes in E(T9) in the percentage range of 0.0–90.0%. As a result, the reduction of the expected time of the vehicle’s stay in the S9 state by 50.0% resulted in the increase of all indicators by over 0.041, and the reduction by 90% increased the values of the indicators by over 0.077.
The attempt to use the Markov process, despite not meeting the condition of exponential time characteristics, showed significant discrepancies between MM and SMM. The high degree of mismatch between the Markov model and the actual process is evidenced by the value of the MAPE error amounting to 351.92%. Therefore, it is inappropriate to use MM to describe the nine-state process of operation of heavy goods vehicles operating in the military transport system. The performed non-parametric Kolmogorov test resulted in the rejection of the hypothesis concerning the compatibility of empirical distributions of the duration of individual states with the theoretical exponential distribution.
The proposed method allows a detailed analysis of vehicle operation as a stochastic process with a multi-state phase space. The unquestionable advantage of the nine-state semi-Markov model is the possibility of evaluating the operation process using indicators of functional readiness, efficiency, and technical efficiency, calculated based on the basis of ergodic probabilities of the process. Additionally, the model sensitivity analysis makes it possible to determine the impact of reducing the values of expected vehicle dwell times in each state on the efficiency of the operation process. However, a disadvantage of the proposed method is the inability to predict the values of indicators for the increased or decreased intensity of vehicle use expressed by means of average daily mileage. This is due to the condition of determining the characteristics of all states in the same time domain.
The proposed methodology for creating stochastic exploitation models can be applied to a wide range of facilities and technical devices. The developed model can be used to analyze and evaluate the operation process of other vehicles operating in technical systems with an analogous or similar operation strategy.

Author Contributions

Conceptualization, M.O.; methodology, M.O. and J.Z.; software, M.O. and J.Z.; validation, M.O.; formal analysis, M.O.; investigation, M.O.; resources, M.O.; data curation, M.O. and J.Z.; writing—original draft preparation, M.O.; writing—review and editing, M.O., J.Z. and J.M.; visualization, M.O.; supervision, J.Z. and J.M.; project administration, J.Z. and J.M.; funding acquisition, J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Military University of Technology, grant number UGB 22-771.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Oszczypała, M.; Ziółkowski, J.; Małachowski, J. Reliability Analysis of Military Vehicles Based on Censored Failures Data. Appl. Sci. 2022, 12, 2622. [Google Scholar] [CrossRef]
  2. Betkier, I.; Mitkow, S.; Kijek, M. Analysis of Vehicle Stability Loss Due to Strong Crosswind Gusts Using Web Services in the Route Planning Process. Arch. Transp. 2019, 52, 47–56. [Google Scholar] [CrossRef]
  3. Betkier, I.; Zak, J.K.; Mitkow, S. Parking Lots Assignment Algorithm for Vehicles Requiring Specific Parking Conditions in Vehicle Routing Problem. IEEE Access 2021, 9, 161469–161487. [Google Scholar] [CrossRef]
  4. Wróblewski, P.; Lewicki, W. A Method of Analyzing the Residual Values of Low-Emission Vehicles Based on a Selected Expert Method Taking into Account Stochastic Operational Parameters. Energies 2021, 14, 6859. [Google Scholar] [CrossRef]
  5. Dziubak, T. The Effects of Dust Extraction on Multi-Cyclone and Non-Woven Fabric Panel Filter Performance in the Air Filters Used in Special Vehicles. Eksploat. Niezawodn. 2016, 18, 348–357. [Google Scholar] [CrossRef]
  6. Wang, L.; Yang, Q.; Tian, Y. Reliability Analysis of 6-Component Star Markov Repairable System with Spatial Dependence. Math. Probl. Eng. 2017, 2017, 9728019. [Google Scholar] [CrossRef]
  7. Galiev, I.F.; Sabitov, A.E. Analysis of the Reliability and Efficiency of Local Power Supply Systems at Major International Events. Lect. Notes Civ. Eng. 2022, 190, 269–278. [Google Scholar] [CrossRef]
  8. Ziółkowski, J.; Małachowski, J.; Oszczypała, M.; Szkutnik-Rogoż, J.; Lęgas, A. Modelling of the Military Helicopter Operation Process in Terms of Readiness. Def. Sci. J. 2021, 71, 602–611. [Google Scholar] [CrossRef]
  9. Ziółkowski, J.; Żurek, J.; Małachowski, J.; Oszczypała, M.; Szkutnik-Rogoż, J. Method for Calculating the Required Number of Transport Vehicles Supplying Aviation Fuel to Aircraft during Combat Tasks. Sustainability 2022, 14, 1619. [Google Scholar] [CrossRef]
  10. Yahaya, M.; Azis, N.; Mohd Selva, A.; Ab Kadir, M.; Jasni, J.; Kadim, E.; Hairi, M.; Yang Ghazali, Y. A Maintenance Cost Study of Transformers Based on Markov Model Utilizing Frequency of Transition Approach. Energies 2018, 11, 2006. [Google Scholar] [CrossRef] [Green Version]
  11. Dziubak, T.; Dziubak, S.D. A Study on the Effect of Inlet Air Pollution on the Engine Component Wear and Operation. Energies 2022, 15, 1182. [Google Scholar] [CrossRef]
  12. Wróblewski, P.; Drożdż, W.; Lewicki, W.; Miązek, P. Methodology for Assessing the Impact of Aperiodic Phenomena on the Energy Balance of Propulsion Engines in Vehicle Electromobility Systems for given Areas. Energies 2021, 14, 2314. [Google Scholar] [CrossRef]
  13. Paś, J.; Rosiński, A.; Wiśnios, M.; Stawowy, M. Assessing the Operation System of Fire Alarm Systems for Detection Line and Circuit Devices with Various Damage Intensities. Energies 2022, 15, 3066. [Google Scholar] [CrossRef]
  14. Mengistu, T.M.; Che, D.; Alahmadi, A.; Lu, S. Semi-Markov Process Based Reliability and Availability Prediction for Volunteer Cloud Systems. In Proceedings of the 2018 IEEE 11th International Conference on Cloud Computing (CLOUD), San Francisco, CA, USA, 2–7 July 2018; IEEE: Piscataway, NJ, USA; pp. 359–366. [Google Scholar]
  15. Zhang, T.; Liu, T.; Liu, D.; Sun, F. Reliability Evaluation of Modular Multilevel Converter System Based on Semi-Markov Model. In Proceedings of the 2021 IEEE 4th International Conference on Electronics Technology (ICET), Chengdu, China, 7–10 May 2021; pp. 513–517. [Google Scholar]
  16. Papadopoulos, C.T.; Li, J.; O’Kelly, M.E.J. A Classification and Review of Timed Markov Models of Manufacturing Systems. Comput. Ind. Eng. 2019, 128, 219–244. [Google Scholar] [CrossRef]
  17. Wang, J.; Miao, Y. Optimal Preventive Maintenance Policy of the Balanced System under the Semi-Markov Model. Reliab. Eng. Syst. Saf. 2021, 213, 107690. [Google Scholar] [CrossRef]
  18. Zhiyong, G.; Jiwu, L.; Rongxi, W. Prognostics Uncertainty Reduction by Right-Time Prediction of Remaining Useful Life Based on Hidden Markov Model and Proportional Hazard Model. Eksploat. Niezawodn. 2021, 23, 154–165. [Google Scholar] [CrossRef]
  19. Seyr, H.; Muskulus, M. Use of Markov Decision Processes in the Evaluation of Corrective Maintenance Scheduling Policies for Offshore Wind Farms. Energies 2019, 12, 2993. [Google Scholar] [CrossRef] [Green Version]
  20. Li, Y.; Cui, L.; Lin, C. Modeling and Analysis for Multi-State Systems with Discrete-Time Markov Regime-Switching. Reliab. Eng. Syst. Saf. 2017, 166, 41–49. [Google Scholar] [CrossRef]
  21. El-Awady, A.; Ponnambalam, K. Integration of Simulation and Markov Chains to Support Bayesian Networks for Probabilistic Failure Analysis of Complex Systems. Reliab. Eng. Syst. Saf. 2021, 211, 107511. [Google Scholar] [CrossRef]
  22. Wang, J.; Feng, Y.; Fei, Z.; Li, S.; Chang, Q. Markov Chain Based Idle Status Control of Stochastic Machines for Energy Saving Operation. In Proceedings of the 2017 13th IEEE Conference on Automation Science and Engineering (CASE), Xi’an, China, 20–23 August 2017; IEEE: Piscataway, NJ, USA; pp. 1019–1023. [Google Scholar]
  23. Kamlu, S.; Laxmi, V. Condition-Based Maintenance Strategy for Vehicles Using Hidden Markov Models. Adv. Mech. Eng. 2019, 11, 1–13. [Google Scholar] [CrossRef]
  24. Borucka, A.; Niewczas, A.; Hasilova, K. Forecasting the Readiness of Special Vehicles Using the Semi-Markov Model. Eksploat. Niezawodn.-Maint. Reliab. 2019, 21, 662–669. [Google Scholar] [CrossRef]
  25. Yang, W.; Chen, L. Machine Condition Recognition via Hidden Semi-Markov Model. Comput. Ind. Eng. 2021, 158, 107430. [Google Scholar] [CrossRef]
  26. Sánchez-Herguedas, A.; Mena-Nieto, A.; Rodrigo-Muñoz, F. A New Analytical Method to Optimise the Preventive Maintenance Interval by Using a Semi-Markov Process and z-Transform with an Application to Marine Diesel Engines. Reliab. Eng. Syst. Saf. 2021, 207, 107394. [Google Scholar] [CrossRef]
  27. Grabski, F. Semi-Markov Failure Rates Processes. Appl. Math. Comput. 2011, 217, 9956–9965. [Google Scholar] [CrossRef]
  28. Migawa, K.; Borowski, S.; Neubauer, A.; Sołtysiak, A. Semi-Markov Model of the System of Repairs and Preventive Replacements by Age of City Buses. Appl. Sci. 2021, 11, 10411. [Google Scholar] [CrossRef]
  29. Świderski, A.; Borucka, A.; Grzelak, M.; Gil, L. Evaluation of Machinery Readiness Using Semi-Markov Processes. Appl. Sci. 2020, 10, 1541. [Google Scholar] [CrossRef] [Green Version]
  30. Migawa, K. Availability Control for Means of Transport in Decisive Semi-Markov Models of Exploitation Process. Arch. Transp. 2012, 24, 497–508. [Google Scholar] [CrossRef] [Green Version]
  31. Zhao, L.; Li, K.; Zhao, W.; Ke, H.-C.; Wang, Z. A Sticky Sampling and Markov State Transition Matrix Based Driving Cycle Construction Method for EV. Energies 2022, 15, 1057. [Google Scholar] [CrossRef]
  32. Wang, Y.; Xie, B.; E, S. Adaptive Relevance Vector Machine Combined with Markov-Chain-Based Importance Sampling for Reliability Analysis. Reliab. Eng. Syst. Saf. 2022, 220, 108287. [Google Scholar] [CrossRef]
  33. Şahin, İ. Markov Chain Model for Delay Distribution in Train Schedules: Assessing the Effectiveness of Time Allowances. J. Rail Transp. Plan. Manag. 2017, 7, 101–113. [Google Scholar] [CrossRef]
  34. Rajarshi, M.B. Statistical Inference for Discrete Time Stochastic Processes; Springer Briefs in Statistics; Springer: Delhi, India, 2013; ISBN 978-81-322-0762-7. [Google Scholar]
  35. Shepero, M.; Munkhammar, J. Spatial Markov Chain Model for Electric Vehicle Charging in Cities Using Geographical Information System (GIS) Data. Appl. Energy 2018, 231, 1089–1099. [Google Scholar] [CrossRef]
  36. Pourmousavi Kani, S.A.; Ardehali, M.M. Very Short-Term Wind Speed Prediction: A New Artificial Neural Network–Markov Chain Model. Energy Convers. Manag. 2011, 52, 738–745. [Google Scholar] [CrossRef]
  37. Karakaya, E.; Vinel, A.; Smith, A.E. Relocations in Container Depots for Different Handling Equipment Types: Markov Models. Comput. Ind. Eng. 2021, 157, 107311. [Google Scholar] [CrossRef]
  38. González-Domínguez, J.; Sánchez-Barroso, G.; García-Sanz-Calcedo, J. Scheduling of Preventive Maintenance in Healthcare Buildings Using Markov Chain. Appl. Sci. 2020, 10, 5263. [Google Scholar] [CrossRef]
  39. Iversen, E.B.; Morales, J.M.; Madsen, H. Optimal Charging of an Electric Vehicle Using a Markov Decision Process. Appl. Energy 2014, 123, 1–12. [Google Scholar] [CrossRef] [Green Version]
  40. Brayer, E.F. Calculating the Standard Error of a Proportion. Appl. Stat. 1957, 6, 67–68. [Google Scholar] [CrossRef]
  41. Gurland, J.; Tripathi, R.C. A Simple Approximation for Unbiased Estimation of the Standard Deviation. Am. Stat. 1971, 25, 30–32. [Google Scholar] [CrossRef]
  42. van Staden, H.E.; Boute, R.N. The Effect of Multi-Sensor Data on Condition-Based Maintenance Policies. Eur. J. Oper. Res. 2021, 290, 585–600. [Google Scholar] [CrossRef]
  43. Jafari, L.; Makis, V.; Gb, A.K. Optimal Maintenance Policy for a Partially Observable Two-Unit System. World Acad. Sci. Eng. Technol. Int. J. Mech. Aerosp. Ind. Mechatron. Manuf. Eng. 2014, 8, 1612–1619. [Google Scholar]
  44. Agrawal, A.K.; Murthy, V.M.S.R.; Chattopadhyaya, S. Investigations into Reliability, Maintainability and Availability of Tunnel Boring Machine Operating in Mixed Ground Condition Using Markov Chains. Eng. Fail. Anal. 2019, 105, 477–489. [Google Scholar] [CrossRef]
  45. Kampitsis, D.; Panagiotidou, S. A Bayesian Condition-Based Maintenance and Monitoring Policy with Variable Sampling Intervals. Reliab. Eng. Syst. Saf. 2022, 218, 108159. [Google Scholar] [CrossRef]
  46. Girtler, J. Physical Aspects of Application and Usefulness of Semi-Markovian Processes for Modelling the Processes Occurring in Operational Phase of Technical Objects. Pol. Marit. Res. 2004, 11, 25–30. [Google Scholar]
  47. Ivanchenko, O.; Kharchenko, V.; Moroz, B.; Kabak, L.; Smoktii, K. Semi-Markov Availability Model Considering Deliberate Malicious Impacts on an Infrastructure-as-a-Service Cloud. In Proceedings of the 2018 14th International Conference on Advanced Trends in Radioelecrtronics, Telecommunications and Computer Engineering (TCSET), Lviv-Slavske, Ukraine, 20–24 February 2018; pp. 570–573. [Google Scholar]
  48. Conover, W.J. A Kolmogorov Goodness-of-Fit Test for Discontinuous Distributions. J. Am. Stat. Assoc. 1972, 67, 591–596. [Google Scholar] [CrossRef]
  49. Quine, M.P.; Robinson, J. Efficiencies of Chi-Square and Likelihood Ratio Goodness-of-Fit Tests. Ann. Stat. 1985, 13, 727–742. [Google Scholar] [CrossRef]
  50. Shepelev, S.; Shepelev, V.; Cherkasov, Y. Differentiation of the Seasonal Loading of Combine Harvester Depending on Its Technical Readiness. Procedia Eng. 2015, 129, 161–165. [Google Scholar] [CrossRef] [Green Version]
  51. Zhao, J.; Jia, Y.; Shi, L.; Zhao, X.; Kong, L. Research on Technical Readiness Evaluation Method in Model Development. In Proceedings of the 2015 8th International Symposium on Computational Intelligence and Design (ISCID), Hangzhou, China, 12–13 December 2015; Volume 1, pp. 392–396. [Google Scholar]
  52. Selech, J.; Andrzejczak, K. An Aggregate Criterion for Selecting a Distribution for Times to Failure of Components of Rail Vehicles. Eksploat. Niezawodn.-Maint. Reliab. 2019, 22, 102–111. [Google Scholar] [CrossRef]
  53. Fadhil, R.M.; Unami, K. A Multi-State Markov Chain Model to Assess Drought Risks in Rainfed Agriculture: A Case Study in the Nineveh Plains of Northern Iraq. Stoch. Environ. Res. Risk Assess. 2021, 35, 1931–1951. [Google Scholar] [CrossRef]
  54. Huang, X.-W.; Emura, T. Model Diagnostic Procedures for Copula-Based Markov Chain Models for Statistical Process Control. Commun. Stat. Simul. Comput. 2021, 50, 2345–2367. [Google Scholar] [CrossRef]
  55. Roy, V.; Tan, A.; Flegal, J.M. Estimating Standard Errors for Importance Sampling Estimators with Multiple Markov Chains. Stat. Sin. 2018, 28, 1079–1101. [Google Scholar] [CrossRef] [Green Version]
  56. Rosychuk, R.J.; Sheng, X.; Stuber, J.L. Comparison of Variance Estimation Approaches in a Two-State Markov Model for Longitudinal Data with Misclassification. Stat. Med. 2006, 25, 1906–1921. [Google Scholar] [CrossRef]
Figure 1. Flowchart of stochastic modeling.
Figure 1. Flowchart of stochastic modeling.
Energies 15 05062 g001
Figure 2. Fitting exponential distributions to sojourn times.
Figure 2. Fitting exponential distributions to sojourn times.
Energies 15 05062 g002
Figure 3. Transition diagram for nine-state semi-Markov model of light utility vehicle operations process.
Figure 3. Transition diagram for nine-state semi-Markov model of light utility vehicle operations process.
Energies 15 05062 g003
Figure 4. The course of changes in the probabilities of the states of the embedded Markov chain for the initial vector p0 = [0,0,1,0,0,0,0,0,0].
Figure 4. The course of changes in the probabilities of the states of the embedded Markov chain for the initial vector p0 = [0,0,1,0,0,0,0,0,0].
Energies 15 05062 g004
Figure 5. The course of changes in the probabilities of S1S9 states for the initial vector p0 = [0,0,1,0,0,0,0,0,0]. (a) p1(t); (b) p2(t); (c) p3(t); (d) p4(t); (e) p5(t); (f) p6(t); (g) p7(t); (h) p8(t); (i) p9(t).
Figure 5. The course of changes in the probabilities of S1S9 states for the initial vector p0 = [0,0,1,0,0,0,0,0,0]. (a) p1(t); (b) p2(t); (c) p3(t); (d) p4(t); (e) p5(t); (f) p6(t); (g) p7(t); (h) p8(t); (i) p9(t).
Energies 15 05062 g005
Figure 6. Subsets of light utility vehicle exploitation states.
Figure 6. Subsets of light utility vehicle exploitation states.
Energies 15 05062 g006
Figure 7. Ergodic probabilities of SMM obtained in simulation by reduction of the expected value of sojourn time E(T9).
Figure 7. Ergodic probabilities of SMM obtained in simulation by reduction of the expected value of sojourn time E(T9).
Energies 15 05062 g007
Figure 8. Readiness indicators values obtained in sensitivity analysis of SMM.
Figure 8. Readiness indicators values obtained in sensitivity analysis of SMM.
Energies 15 05062 g008
Figure 9. Differences of ergodic probabilities in Markov and semi-Markov models.
Figure 9. Differences of ergodic probabilities in Markov and semi-Markov models.
Energies 15 05062 g009
Table 1. Literature review of Markov and semi-Markov modeling in engineering.
Table 1. Literature review of Markov and semi-Markov modeling in engineering.
PaperCase StudyModelResults and Conclusions
[22]Machines in manufacturing systemTwo-state Markov chainThe model was developed to simulate and optimize energy saving.
[8]Military helicoptersNine-state Markov modelBoundary probabilities were calculated for all states and the functional readiness index reached 0.9223.
[17]Spring shock absorber in wind turbine systemThree-state semi-Markov model of maintenanceErlang distribution of sojourn time. The system operational availability reached 0.6354–0.6497.
[25]Machines in wafer fabrication work centreThree-state hidden semi-Markov modelAccuracy of the machine condition recognition was 95.56% and prediction accuracy of job processing time was 94.91%.
[29]Production machinesSix-state Markov and semi-Markov modelsReadiness index: 0.85 (semi-Markov model) and 0.45 (Markov model). Groundless assumption of exponential distribution lead to incorrect results.
[24]Special vehiclesThree-state semi-Markov modelTechnical readiness factor was 0.95.
[30]Means of transport16-state semi-Markov model of exploitation processGenetic algorithm was proposed for determining the optimal strategy to control availability.
[23]Buses in transportation systemThree-state hidden Markov modelProbability of availability reached values in the range 0.896–0.969.
[28]Four types of city bus renewal processesFive-state semi-Markov model (Weibull distribution) The model allows to optimize the profit per unit time and readiness to carry out transport tasks depending on the time to preventive maintenance.
[27]Reliability of technical objectsThree-state semi-Markov modelReliability function for Poisson and Furry-Yule failure rate processes.
[26]Marine diesel engines Three-state semi-Markov model (Weibull distribution)The optimal preventive maintenance interval was 1095 h.
[31]Driving cycles of electric vehiclesMultistate Markov models of accelerationThe models describe the stochastic nature of driving cycles in four scenarios: rural, highway, urban and congestion.
[19]Offshore wind farmsSix-state Markov modelThe proposed method allows for improving maintenance efficiency of offshore wind farms.
[10]TransformersFive-state Markov chainPrediction of maintenance cost in 20-year forecast horizon.
Table 2. State space of the light utility vehicles operation process.
Table 2. State space of the light utility vehicles operation process.
StateNamePossible Transition fromPossible Transition to
S1Task executionS2, S4, S6, S7S2, S5, S6, S8
S2RefuelingS1, S3, S4S1, S3, S5
S3Parking in garageS2, S5, S6, S7S2, S4, S6
S4Pre-task serviceS3S1, S2
S5Service after taskS1, S2, S6S3, S9
S6Periodic maintenanceS1, S3, S7S1, S3, S5,
S7RepairS8, S9S1, S3, S6, S9
S8DiagnosticsS1, S9S7, S9
S9Awaiting repairS5, S7, S8S7, S8
Table 3. Example of a part of database sheet.
Table 3. Example of a part of database sheet.
No.DateTimeS1S2S3S4S5S6S7S8S9
26730.03.202008:40:00 1
26830.03.202008:45:001
26930.03.202013:31:00 1
27030.03.202013:35:00 1
27130.03.202013:40:00 1
27231.03.202007:00:00 1
27331.03.202007:30:00 1
27431.03.202011:30:00 1
27504.04.202007:20:00 1
27604.04.202007:25:001
27705.04.202008:40:00 1
27805.04.202008:45:00 1
Table 4. Results of Kolmogorov test.
Table 4. Results of Kolmogorov test.
TijyUy StatisticsCritical RangeHypothesis
T127980.2585<0.0479, 1>H1
T1533400.2019<0.0234, 1>H1
T16320.1834<0.2343, 1>H0
T18320.3587<0.2343, 1>H1
T213960.4773<0.0678, 1>H1
T233660.4970<0.0705, 1>H1
T256950.4731<0.0513, 1>H1
T323650.1428<0.0706, 1>H1
T3440630.2867<0.0213, 1>H1
T36760.2907<0.1534, 1>H1
T4137690.6316<0.0221, 1>H1
T422940.6151<0.0786, 1>H1
T5339880.4533<0.0215, 1>H1
T59700.5071<0.1598, 1>H1
T61110.4395<0.3914, 1>H1
T63840.2843<0.1461, 1>H1
T65230.4757<0.2750, 1>H1
T71270.2320<0.2544, 1>H0
T73670.3295<0.1632, 1>H1
T76100.4754<0.4094, 1>H1
T79440.4082<0.2006, 1>H1
T87630.3348<0.1683, 1>H1
T89390.2445<0.2128, 1>H1
T97850.2186<0.1452, 1>H1
T98700.5196<0.1598, 1>H1
Table 5. Standard errors of probabilities estimation.
Table 5. Standard errors of probabilities estimation.
SE(pij)S1S2S3S4S5S6S7S8S9
S100.0061000.00620.001300.00130
S20.011700.011400.01310000
S300.004100.004400.0019000
S40.00410.00410000000
S5000.0020000000.0020
S60.026800.041700.03650000
S70.031700.0409000.0206000.0376
S80000000.048100.0481
S90000000.04000.04000
Table 6. Ergodic probabilities of states in embedded Markov chain and SMM.
Table 6. Ergodic probabilities of states in embedded Markov chain and SMM.
StateS1S2S3S4S5S6S7S8S9
pj0.2235370.0775050.2395750.2161170.2158720.0062670.0077220.0053430.008063
E(Tj) (min)717.524.603956.585.008.13271.95174.6336.1113,240.50
πj0.1313110.0002920.7760240.0008850.0014360.0013950.0011040.0001580.087396
πj (%)13.13110.029277.60240.08850.14360.13950.11040.01588.7396
Table 7. Values of Kr, Ke and Ks indicators in SMM.
Table 7. Values of Kr, Ke and Ks indicators in SMM.
IndicatorValue
Functional readiness Kr0.907334
Technical efficiency Ke0.909947
Technical suitability Ks0.911343
Table 8. Ergodic probabilities obtained in sensitivity analysis of SMM.
Table 8. Ergodic probabilities obtained in sensitivity analysis of SMM.
ΔE(T9) (%)π1π2π3π4π5π6π7π8π9
00.1313110.0002920.7760240.0008850.0014360.0013950.0011040.0001580.087396
100.1324680.0002940.7828650.0008930.0014490.0014070.0011140.0001590.079350
200.1336460.0002970.7898280.0009010.0014620.0014200.0011240.0001610.071162
300.1348460.0003000.7969160.0009090.0014750.0014330.0011340.0001620.062826
400.1360670.0003020.8041330.0009170.0014880.0014460.0011440.0001640.054339
500.1373100.0003050.8114810.0009260.0015020.0014590.0011540.0001650.045698
600.1385770.0003080.8189650.0009340.0015150.0014720.0011650.0001670.036897
700.1398670.0003110.8265880.0009430.0015300.0014860.0011760.0001680.027932
800.1411810.0003140.8343540.0009520.0015440.0015000.0011870.0001700.018799
900.1425200.0003170.8422680.0009610.0015590.0015140.0011980.0001710.009492
Table 9. Indicators values obtained in sensitivity analysis of SMM.
Table 9. Indicators values obtained in sensitivity analysis of SMM.
ΔE(T9) (%)KrΔKrKeΔKeKsΔKs
00.9073340.0000000.9099470.0000000.9113430.000000
100.9153340.0080000.9179700.0080230.9193770.008035
200.9234760.0161420.9261350.0161880.9275550.016213
300.9317640.0244300.9344470.0245000.9358800.024538
400.9402020.0328680.9429100.0329620.9443560.033013
500.9487940.0414600.9515270.0415800.9529860.041643
600.9575450.0502110.9603030.0503560.9617750.050433
700.9664590.0591250.9692430.0592950.9707290.059386
800.9755410.0682060.9783500.0684030.9798500.068508
900.9847940.0774600.9876300.0776830.9891450.077802
Table 10. Comparison of results obtained by MM and SMM.
Table 10. Comparison of results obtained by MM and SMM.
π1π2π3π4π5π6π7π8π9
SMM0.1313110.0002920.7760240.0008850.0014360.0013950.0011040.0001580.087396
MM0.0127840.0001720.2759310.0001770.0004350.0174890.0037800.0013250.687906
difference−0.118527−0.000120−0.500093−0.000708−0.0010010.0160940.0026760.0011670.600510
difference in %−90.26−41.10−64.44−80.00−69.711153.69242.39738.61687.11
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Oszczypała, M.; Ziółkowski, J.; Małachowski, J. Analysis of Light Utility Vehicle Readiness in Military Transportation Systems Using Markov and Semi-Markov Processes. Energies 2022, 15, 5062. https://doi.org/10.3390/en15145062

AMA Style

Oszczypała M, Ziółkowski J, Małachowski J. Analysis of Light Utility Vehicle Readiness in Military Transportation Systems Using Markov and Semi-Markov Processes. Energies. 2022; 15(14):5062. https://doi.org/10.3390/en15145062

Chicago/Turabian Style

Oszczypała, Mateusz, Jarosław Ziółkowski, and Jerzy Małachowski. 2022. "Analysis of Light Utility Vehicle Readiness in Military Transportation Systems Using Markov and Semi-Markov Processes" Energies 15, no. 14: 5062. https://doi.org/10.3390/en15145062

APA Style

Oszczypała, M., Ziółkowski, J., & Małachowski, J. (2022). Analysis of Light Utility Vehicle Readiness in Military Transportation Systems Using Markov and Semi-Markov Processes. Energies, 15(14), 5062. https://doi.org/10.3390/en15145062

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