DOI: 10.15632/jtam-pl.53.3.723
GAS TURBINE RELIABILITY MODEL BASED ON TANGENT HYPERBOLIC RELIABILITY FUNCTION
Ahmed Zohair Djeddi, Ahmed Hafaifa
University of Djelfa, Faculty of Science and Technology, Djelfa, Algeria e-mail: a z djeddi@univ-djelfa.dz; hafaifa.ahmed.dz@ieee.org
Abudura Salam
University of M´ed´ea, Faculty of Science and Technology, M´ed´ea, Algeria e-mail: abudura.salam@univ-medea.dz
The present work deals with the exploration of a new model proposed for the reliability analysis of industrial production systems. This proposed model is mainly based on the tan-gent hyperbolic function, where the survival function is determined and used in the lifetime distribution modeling taking into account of estimation the parameters of the proposed func-tion. On the other side, tests validation is performed using the real data of a gas turbine installation. The obtained results allow the modeling of damage effects, hence the prediction of the performance of the examined gas turbine using the proposed model gives good results in terms of validity.
Keywords: reliability estimation, reliability algorithms, lifetime distribution, Weibull
distri-bution, availability
1. Introduction
Actually, the complexity of industrial plants and their equipment behavior led practically to very complex maintenance strategies that are containing a number of different tasks. Moreover, the random nature of degradation and failure that may occur makes the determination of the required strategy to fulfill the best decisions regarding their maintenance a very difficult practical task (Costa et al., 2014; Lai, 1994; Guemana et al., 2011; Rao et al., 2005; Hasumi et al., 2009). The main aim of the proposed model is to analyze and measure industrial system reliability parameters to achieve the best time determination for making decision on the maintenance actions. Whereas, in the present work, this study is based on a model of a gas turbine which is considered as the main important equipment of the compressor stations, gas pumps, oil pumps and petrol production. The reliability model proposed in this paper is based essentially on the elementary function of tangent hyperbolic which is proved in this paper to be theoretically and practically similar to the well know Weibull distribution which was introduced by the Swedish mathematician in 1951 (Weibull, 1951). The Weibull distribution is considered the first distribution used in reliability analysis due to its very big capacity to adapt a very large number of data sets. However, many lifetime distributions that have a bathtub-shaped hazard rate function have been introduced after Weibull distribution, in which the main aim was to improve the systems reliability. A survey presenting the state-of-the-art on the class of such distributions was presented by Moeini et al. (2013), Scott (1979), Lai et al. (1998, 2001, 2003). This work proposes an exploration of a new model which will contribute to the enhancement of the industrial systems reliability. It is obvious that the hyperbolic tangent function has a curve similar to the cumulative function of the probability distribution. Hence, a new model based on a tangent hyperbolic function defined on ℜ+ can be proposed. It is very clear that
its form is similar to the form of Weibull’s model. In this model, the survival function will be studied based on it.
For validation of the proposed model, a practical application has been made with real data. These data have been collected from a gas turbine operating in a natural gas transportation system which is used in the south of Algeria. The reliability approach developed in this paper allows the modeling of the effects of damage to predict the performance of the examined gas turbine operation and give good results in terms of validity compared with the Weibull approach.
2. Reliability model based on the tangent hyperbolic function
The availability of the industrial equipment control allows the industries to act positively on the conformity of production, the exploitation costs and the production successful competitiveness. Indeed, the life managing of the equipment that is being used in oil and gas facilities is based on taking into account the aspects of aging, the economic and regulatory factors in order to optimize the operation, maintenance and lifespan of the system structures and components. The main goal is to maintain the level of safety and reliability as well as to maximize the investment return on the overall lifespan. The reliability of industrial systems is always essential and attracts much attention among scientific researches. It is especially important for components whose failure can cause major problems in terms of maintainability, availability and security (Halimi et al., 2014; Hafaifa et al., 2013a,b; Sturges, 1926; Ruji, 1990; Trofimov et al., 1978). As it is aforementioned, a new model based on the tangent hyperbolic function is proposed to be used for the analysis and the study of the survival function in the reliability modeling. The tangent hyperbolic function considered in this paper is defined as a function on ℜ+ (Fig. 1). It is expressed as follows
F(x) = (
0 for x∈ ℜ+
tanh(x) for x 0 (2.1)
Fig. 1. Tangent hyperbolic function
To assess the impact of the aging effects on the equipment by annual frequency and to identify the sensitive equipment and the requested strategies for service priority to manage the risks that are associated with the age of the facilities, a cumulative distribution function is used. In this case it is the limit of F (x) equal to 1 when x tends to +infinity and is equal to 0 when x tends to 0. The proposed function is a non-decreasing and a right-continuous function, so the proposed tangent hyperbolic function is fulfils these requirements and it presents a good choice as a candidate cumulative distribution function (Costa et al., 2014; Guaily and Epstein, 2013)
lim
The distribution lifespan in the proposed tangent hyperbolic distribution form is deduced from the Weibull distribution. The survival function is represented as follows
R(t) = 1 − tanh(λt)β (2.3)
where β is the shape parameter and λ is the scale parameter. The cumulative distribution function is expressed as follows
F(t) = 1 − R(t) = tanh(λt)β (2.4)
Based on the derivative of the last function, the probability density function is obtained as follows:.
f(t) = dF(t) dt =
(λt)β−1λβ
cosh2(λt)β (2.5)
The hazard function is determined in the flowing equations Pr(t < T ¬ t + ∆t T > t) = f(t) R(t)∆t= h(t)∆t (2.6) where h(t) = f(t) R(t) = (λt)β−1 λβ cosh2(λt)β 1 − tanh(λt)β
Using a variable converting the reliability R(t) replaced by the hyperbolic function, h(t) can be rewritten using exponential form as follows
h(t) = 2λβ(λt)
β−1
1 + e−2(λt)β (2.7)
2.1. Hazard and density functions
The hazard function in the proposed tangent hyperbolic model is expressed by equation (2.7) and is shown in Fig. 2. It is obvious that the curve of this function depends only on the parameter β. Therefore, the curve of this function can be discussed based on the value of β:
• In the case when β = 1, the hazard function tends quickly to 2λ lim
t→∞h(t) = 2λ (2.8)
• In the case when β = 2, the hazard function becomes h(t) = 4tλ
2
1 + e−2(λt)2 → 4tλ
2 (2.9)
The hazard function h(t) is a straight line with slope equal to 4λ2. • In the case when 0 < β < 1, the hazard function becomes
h(t) → 2λβ(λt)β−1 (2.10)
The hazard function curve follows tp form with −1 < p < 0.
• In the case when β > 2, the hazard function curve follows tp with p > 1 (convex form).
• In the case when 1 < β < 2, the hazard function curve follows tp with 0 > p > 1 (concave
form).
It can be concluded that the curves of the hazard function have their Weibull equivalents cu-rves. On the other side, the density function of the proposed tangent hyperbolic model presented in equation (2.5) is presented in Fig. 3 for different values of β.
Fig. 2. Hazard function determined using the tangent hyperbolic algorithm
Fig. 3. Density function found using the hyperbolic model 3. Application results
The data used in this paper are taken from the data history of a gas turbine (PGT10) used for natural gas transportation and installed in a gas plant located in the south of Algeria. Those data are presented in Table 1.
Table 1.Examined gas turbine data
N 1 2 3 4 5 6 7 8 9 10 TBF∗ 240 264 384 552 624 648 696 720 720 768 N 11 12 13 14 15 16 17 18 19 20 TBF 768 792 960 1272 1344 1464 1632 1680 1776 1944 N 21 22 23 24 25 26 27 28 29 30 TBF 1968 1992 2064 2136 2208 2376 2448 2448 2472 2664 N 31 32 33 34 35 36 37 38 39 40 TBF 2832 2880 3000 3600 3672 3720 4272 4656 5592 5856 ∗ TBF is the Time Between Failures given in hours
3.1. Esytimation of the model parameters
Using equation (2.4), the inverse of the tangent hyperbolic function can be obtained as follows tanh−1
The tangent hyperbolic can be eliminated from this equation based on the following well known expression tanh−1 y= 1 2log 1 + y 1 − y (3.2)
Furthermore, if the log is used to the two members of equation (3.1), the obtained expression is as follows log1 2log 1 + y 1 − y = β log λ + β log t (3.3)
which can be transformed into a linear expression
Y = A + BX (3.4) where Y = log1 2log 1 + y 1 − y X = log t A= β log λ B= β
To estimate A and B, any graphical or analytical methods can be used, a possible approach is the simple regression analysis using (3.2). The estimation of the parameters based on the least--squares fit is shown in Fig. 4. The parameters of the examined gas turbine using the proposed approach are presented in Table 2. The obtained estimated parameters are λ = 0.00034 and β = 1.24568.
Table 2.The obtained parameters for the examined gas turbine Estimate Std. Error t value Pr(> |t|) β 1.24569 0.03435 36.27 <2E-16 βlog λ −7.98420 0.02801 −285.07 <2E-16
Fig. 4. The distribution function of the examined gas turbine determined using the hyperbolic algorithm The plots of Y versus X which are related to the cumulative distribution function and time for the examined gas turbine determined using the proposed tangent hyperbolic algorithm are shown in Fig. 4, where the plotted data are scattered close to the fitted straight line. It can be said that this model is acceptable (Murthy et al., 2004; Yang and Scott, 2013).
3.1.1. Cumulative distribution function
The cumulative distribution function of the examined gas turbine found using the tangent hyperbolic model is shown in Fig. 5. The plotted data are scattered close to the fitted cu-rve. Following the cumulative distribution function expressed in equation (2.4), the parameters λ= 0.00034 and β = 1.24568.
Fig. 5. The cumulative failure rate of the examined gas turbine using the hyperbolic model
3.1.2. Density and hazard functions
The density function and the hazard function of the examined gas turbine determined by using the tangent hyperbolic function model are shown respectively in Figs. 6 and 7. These functions are obtained by using equation (2.5) and taking into account the estimated parameter values λ = 0.00034 and β = 1.24568.
Fig. 6. The distribution function of the model of the examined gas turbine using the hyperbolic model
Fig. 7. The hazard function of the examined gas turbine using the hyperbolic model
3.1.3. Mean time between failures
The Mean Time Between Failures (MTBF) is the predicted elapsed time between inherent failures of the system during exploitation (Yang and Scott, 2013) which is expressed as follows
MTBF = +∞ Z 0 tf(t) dt = +∞ Z 0 t(λt) β−1λβ cosh2(λt)β dt (3.5)
It is worth noting that there is an alternative way for computing the expected value (Murthy et al., 2004) MTBF = +∞ Z 0 S(t) dt = +∞ Z 0 1 − tanh(λt)β dt (3.6) Table 3
Model Weibull Hyperbolic function MTBF 2088.623 h 2075.40 h
4. Conclusion
The goal achieved in this paper is concerned with industrial reliability, which requires the spe-cialists to be more competitive and more responsive to variable market conditions. It is obvious that to achieve better management of the facilities performance, it is necessary to improve the reliability. It has been shown that with the application of the proposed method based on the tangent hyperbolic function good results can be obtained. The main objective is to analyze and to measure the reliability parameters of the examined system in order to determine the best time to take maintenance action based on reliability analysis. The developed in this paper reliability approach allows the modeling of the damage effects to predict the performance of the examined gas turbine operation and to lead to good results in terms of validity compared with the Weibull approach. This model can be improved to have a bathtub-shaped hazard form.
Refernces
1. Costa F.M.P., Rocha A.M.A.C., Fernandes E.M.G.P., 2014, An artificial fish swarm al-gorithm based hyperbolic augmented Lagrangian method, Journal of Computational and Applied
Mathematics, 259, Part B, 868-876
2. Guaily A.G., Epstein M., 2013, Boundary conditions for hyperbolic systems of partial differen-tials equations, Journal of Advanced Research, 4, 4, 321-329
3. Guemana M., Aissani S., Hafaifa A., 2011, Use a new calibration method for gas pipelines: an advanced method improves calibrating orifice flowmeters while reducing maintenance costs,
Hydrocarbon Processing Journal, 90, 8, 63-68
4. Hafaifa A., Belhadef R., Guemana M., 2013a, Reliability model exploitation in industrial system maintainability using expert system evaluation, Proceedings of the 4th International
Con-ference on Integrity, Reliability and Failure, IRF2013, Funchal, Madeira, Portugal, 387-388
5. Hafaifa A., Guemana M., Daoudi A., 2013b, Vibrations supervision in gas turbine ba-sed on parity space approach to increasing efficiency, Journal of Vibration and Control, doi: 10.1177/1077546313499927
6. Halimi D., Hafaifa A., Bouali E., 2014, Maintenance actions planning in industrial centrifugal compressor based on failure analysis, Journal of Maintenance and Reliability, 16, 1, 17-21 7. Hasumi T., Akimoto T., Aizawa Y., 2009, The Weibull-log Weibull distribution for
interoccur-rence times of earthquakes, Physica A: Statistical Mechanics and its Applications, 388, 4, 491-498 8. Lai C.D., 1994, Tests of univariate and bivariate stochastic ageing, IEEE Transactions on
9. Lai C.D., Moore T., Xie M., 1998, The beta integrated model, Prococeding of the International
Workshop on Reliability Modeling and Analysis: from Theory to Practice, 153-159
10. Lai C.D., Xie M., Murthy D.N.P., 2001, Bathtub shaped failure rate distributions, Handbook
in Reliability, 20, 69-104
11. Lai C.D., Xie M., Murthy D.N.P., 2003, A modified Weibull distribution, IEEE Transactions
on Reliability, 52, 1, 33-37
12. Moeini A., Jenab K., Mohammadi M., Foumani M., 2013, Fitting the three-parameter Weibull distribution with Cross Entropy, Applied Mathematical Modelling, 37, 9, 6354-6363
13. Murthy P.D.N., Xie M., Jiang R., 2004, Weibull Models, John Wiley & Sons
14. Rao C.R., Wegman E.J., Solka J.L., 2005, Handbook of Statistics: Data Mining and Data
Visualization, Volume 24 de Handbook of Statistics, Elsevier Science
15. Ruji H., 1990, State space tree method and exact decomposition algorithm for finding network overall reliability, Journal of Electronics (China), 7, 4, 296-305
16. Scott D.W., 1979, On optimal and data-based histograms, Biometrika, 66, 3, 605-610
17. Sturges H.A., 1926, The choice of a class interval, Journal of the American Statistical Association,
21, 153, 65-66
18. Trofimov N.G., Kravchenko B.A., Kramarovskii B.I., Baturov V.B., Kostina G.N., 1978, Increasing the strength and reliability of turbine blades by thermoplastic hardening methods,
Strength of Materials, 10, 8, 990-996
19. Weibull W., 1951, A statistical distribution function of wide applicability, Journal of Applied
Mechanics, 1, 18, 293-297
20. Yang J., Scott D.W., 2013, Robust fitting of a Weibull model with optional censoring,
Compu-tational Statistics and Data Analysis, 67, 1, 149-161