RELIABILITY ANALYSES OF AN INDUSTRIAL SYSTEM BASED ON HERMITE POLYNOMIAL AND 2PARAMETER WEIBULL Текст научной статьи по специальности «Медицинские технологии»

Anas Sani Maihulla, Micheal Khoo Boon Chong and Ibrahim Yusuf

RELIABILITY ANALYSES OF AN INDUSTRIAL SYSTEM BASED ON RT&A, No 3 (74) HERMITE POLYNOMIAL AND 2-PARAMETER WEIBULL_Volume 18, September 2023

RELIABILITY ANALYSES OF AN INDUSTRIAL SYSTEM BASED ON HERMITE POLYNOMIAL AND 2-PARAMETER WEIBULL

Anas Sani Maihulla

Department of Mathematics, Sokoto State University, Sokoto- Nigeria anas.maihulla@ssu.edu.ng

Michael Khoo Boon Chong

School of Mathematical Sciences, University Science Malaysia mkbc@usm.my

Ibrahim Yusuf

Department of Mathematical sciences, Bayero University, Kano- Nigeria

iyusuf.mth@buk.edu.ng

Abstract

The current study aims to assess the structure's reliability using stochastic Hermite surface methodology. This approach uses series expansion of standard normal random variables to model uncertainty. (i.e., polynomial chaos expansion). The coefficients of the polynomial chaos expansion are found through stochastic collocation, which only requires a few performance function evaluations. After determining the order of the polynomial and its coefficients, first order methods calculates the reliability index. To demonstrate the applicability of the suggested based reliability analysis, failure rates were used and the duration for the evaluated were set to be 360 days. Numerical result are provided on monthly basics. On the other hand, to achieve our goal, we proposed the 2-parameter modified Weibull distribution. The simulation was performed using Maple software. The evaluation for each subsystem was displayed in the result and analyses section. The conclusion, however, draws a broad conclusion about the study.

Keywords:Reliability, Variables, Polynomials, Stochastic, Simulation, Distributions

I. Introduction

Derivatives are challenging to assess. For reliability analysis in this situation, Bucher and Bourgand [1] developed the response surface method (RSM). Near the failure region, a multidimensional quadratic polynomial in RSM is used to approximate the unknown limit surface. This method has recently been widely used by engineers and researchers for a variety of applications, including performance evaluation, crash simulation, and reliability-based design optimization. However, since the failure region is close to the polynomial approximation of the original limit state, it is frequently challenging to determine the ideal separation for limit states with multiple design points.Additionally, because this is a deterministic representation, it is unable to accurately depict

Anas Sani Maihulla, Micheal Khoo Boon Chong and Ibrahim Yusuf

RELIABILITY ANALYSES OF AN INDUSTRIAL SYSTEM BASED ON RT&A, No 3 (74)

HERMITE POLYNOMIAL AND 2-PARAMETER WEIBULL_Volume 18, September 2023

the stochastic features of the initial limit state. The Askey-Wilson polynomials, which share the characteristics of more conventional orthogonal polynomials like the Legendre polynomials or the spherical harmonics, are the most general orthogonal polynomials [2].

A new dynamic reliability assessment method based on the decomposition method and the Hermite polynomials approximation is presented for the dynamic reliability analysis of stochastic structures under stationary random excitation. In this method, a multi-dimensional dynamic reliability response function is additively broken down into a one-dimensional function, and the one-dimensional dynamic reliability response function is then approximated using Hermite polynomials. Last but not least, the Maple software simulation yields the unconditional reliability of the explicit response function, and two expressions show the logic of the suggested approach. A systematic account of analytical calculus for distributions can be found in Potthoff-Yan [3]. Ordinary differentiation and integration can be thought of as being generalized to arbitrary order in the calculus of fractional order. The uncertainty of the rational order of the derivation gave rise to G.W [4]. Leibniz's question, which gave rise to the fractional calculus, in 1695 [5]. Recent years have seen a significant increase in interest in both theory and application for Jacobi polynomials. For the purpose of solving the multi-term fractional differential equations, the shifted Jacobi operational matrix of fractional derivatives and the spectral tau approach are combined. [6]. To calculate multi-term FDEs, a new explicit solution is created that is targeted for shifted Chebyshev polynomials with flexible degree and fractional order. The work in [7] can be used to solve the same linear problem. Laguerre polynomials have been used in a number of attempts to solve FDEs using various spectral methods within the realm of numerical methods [8]. The research done in [9] led to the creation of a novel time-dependent problem-targeting algorithm built on spectral Laguerre approximations. A novel tau method was proposed in a recent paper that examined the modified Laguerre functions [10]. Comparative analyses of Weibull distribution with Reliability, Availability, Maintainability, and Dependability analyses of an industrial system was conducted [14]. The reliability analyses of filtration system using Copula approach was conducted [15]. Complex reverse osmosis system was studied by [16] and [17] using RAMD analyses. And that of Photovoltaic using the same RAMD analyses was conducted by [18]. There is no comparison of the industrial system's performance using the Hermite polynomial and the two-parameter Weibull distribution in the literature currently available, this prompt our research. Additionally, notice how an industrial process innovation is developing.

II. Methods

Initial value conditions of multi-term fractional differential equations serve as a driving force behind many practical issues. The Hermite tau method is modified in this section to include the operational matrix in order to solve fractional differential equations. The following lists every step in the entire process.

2.2. The Properties of Hermite Polynomials

Let a = (-x>,x>) and y(l) = e-'2be the reliability function on a. Hermite polynomials of degree m are defined in their analytical form [11]

Rm(l) = ld=0 d](m-2d)] (1)

where R0(l) = 1 and R1(l) = 21

This recurrence relation is satisfied by hermite polynomials. Rm+i(l) = 2lRm(l) + 2mRm-i(l) (2)

Anas Sani Maihulla, Micheal Khoo Boon Chong and Ibrahim Yusuf

RELIABILITY ANALYSES OF AN INDUSTRIAL SYSTEM BASED ON RT&A, No 3 (74)

HERMITE POLYNOMIAL AND 2-PARAMETER WEIBULL_Volume 18, September 2023

An orthogonal system is one in which the set of Hermite polynomials is an orthogonal polynomial Rrn(l)Rj (IMl)dx = hjSmj (3)

where Smj = 0 V j and Smj = m V m = j describes the role of Kronecker and hj = 2'j\ Jn

2.3. Fractional integration hermite operational reliability

We aim to derive a modified reliability equation of an industrial system for Hermite polynomials in this section.

Let A(l) £ L2(a), subsequently, A(l) can be defined as follows using Hermite polynomials

W) = ZT=0ajRj(l) (4)

Then, coefficient a.j can be written as

aj=2-^i-:mX(l)Rj(l)Y(l)dl, j = 0,1..........(5)

The initial (N + l) the only consideration is given to terms of Hermite polynomials, such that the modified failure rate of the industrial system is in the equation below:

¿N(i) = zUajRj(l) = ATp(l) (6)

where

AT = [«o % ... aN] and <p(l) = [R0(l) R±(l) ... RN(l)]T (7)

Integration of Hermite vector p(x) by Jqp(x) the system reliability in table 1 is obtained from the following equation.

Ji<P(l) = pWp(l) (8)

where q indicates a fined integer value and P(q) represents the actual operational reliability of integrated p(/).

2.4. For the reliability analyses using 2-parameter Weibull distribution

It is clearly stated by Quek and Ang [12]. If the lifetimes follow the Weibull distribution, the p.d.f.

f(t) = /3a-PtP-1e-(a)P (9)

Where a and ft are the scale and shape parameter of the Weibull distribution respectively. The size of the units in which the random variable, t, is measured is reflected by the scale parameter, a. The distribution's form changes depending on the shape parameter,

ft. We can create a diverse set of curves that reflect real lifetime failure distributions by modifying the value of ft.

From (1), the Cumulative Distribution Function is given by:

ft\B

F(t) = 1- e-(a)P (10)

From the relation,

R(t) = 1-F(t) (11)

We can substitute (2) into (3), and we have:

R(t) = (12)

Anas Sani Maihulla, Micheal Khoo Boon Chong and Ibrahim Yusuf

RELIABILITY ANALYSES OF AN INDUSTRIAL SYSTEM BASED ON RT&A, No 3 (74)

HERMITE POLYNOMIAL AND 2-PARAMETER WEIBULL_Volume 18, September 2023

Where R (t) is the reliability or survival function.

The failure rate function or the hazard function can therefore be derived from the following relation:

h(t)=l§) (13) Substituting (9) and (12) into (13) we have:

h(t) = --(14)

e {a'

h(t) = pa-BtB-1 (15)

If we consider M0 to be the number of hours in a year (the size of the population). Out of which Ms units (the number of hours that the system is upstate) survive the test. While Mf fail, then reliability function R (t) is given by:

R(t)=Ml = M-Ml (16)

v y m0 M0 v '

Differentiating both sides of (16) and taking M0 fixed, the following equation will result

dR(t) _ 1 dMf

dt M0 dt ( )

The rate at which component fails can therefore be defined as:

—f = —M0 (18)

dt 0 dt K '

Dividing both sides of the above equation byMs, we obtain the instantaneous probability g(t) of failure, this is:

g{t)=±djML=-MLdm (19)

' MS dt MS dt v '

Using equation (16) into (19) we get:

"fi—ssr^ (20)

Integrating both sides of the equation (20) we have:

f g(t)dt = -logR(t) (21)

From the above equation, the R(t) will be:

R(t) = e-f09(t)dt (22)

Where x is variable.

The function g( ) is called the hazard function or failure rate. Equation (22) can be considered as a generic expression of failure as it is applicable to both exponential and non-exponential failure distribution.

For our modified Weibull distribution, we compared (22) with (12)

(t)B = f0x g(t)dt (23)

For,

g(t) = —2, and the below facts:

a R(t) dt

The derivative of a definite integral of a function is the function itself only when the lower limit of the integral is a constant and the upper limit is the variable with respect to which we are differentiating. To summarize Ulrich et al. [13]:

• The derivative of an indefinite integral of a function is the function itself. i.e., -^ff(x)dx = f(x)

• The derivative of a definite integral with constant limits is 0. i.e., ^ f^ f(x)dx = 0

• The derivative of a definite integral where the lower limit is a constant and the upper limit is a variable is a function itself in terms of the given variable (upper bound).

i.e., ~^faf(x)dx = fix) where 'a' is a constant and 'x' is a variable.

M0

Ms

The NMW proposed for the present research follows;

& = -?;(«(<)) (24)

R(t) = -M0& (25)

Anas Sani Maihulla, Micheal Khoo Boon Chong and Ibrahim Yusuf RELIABILITY ANALYSES OF AN INDUSTRIAL SYSTEM BASED ON HERMITE POLYNOMIAL AND 2-PARAMETER WEIBULL

III. Results

The numerical values for the 2-parameter Weibull distribution and Hermite polynomial derived equations were presented in this section. For the two approaches, comparative analyses were conducted. 360 days were allotted for the analyses using the failure rate range of 0.75 to 0.95.

Table 1: Hermite Polynomials reliability with time

Tim el X^ 0.75 0.85 0.95 Exp(-m)

30 0.9483 0.9385 0.8383 0.7653

60 0.9245 0.9134 0.8134 0.7446

90 0.9188 0.9076 0.8075 0.7146

120 0.9176 0.8962 0.8056 0.7087

150 0.9171 0.8759 0.7946 0.6801

180 0.9148 0.8451 0.7676 0.6479

210 0.9117 0.8153 0.7373 0.6178

240 0.9092 0.8009 0.7187 0.5929

270 0.8977 0.7852 0.7089 0.5708

300 0.8832 0.7657 0.6806 0.5479

330 0.8633 0.7352 0.6735 0.5188

360 0.8558 0.7143 0.6566 0.5070

Table 2: 2-Parameter Weibull distribution's reliability with time

Timel a^ 0.75 0.85 0.95 Exp(-m)

30 0.9954 0.7854 0.5855 0.3667

60 0.9857 0.7789 0.5734 0.3584

90 0.9787 0.7607 0.5677 0.3499

120 0.9644 0.7525 0.5537 0.3312

150 0.9567 0.7457 0.5344 0.3290

180 0.9439 0.7346 0.5270 0.3043

210 0.9387 0.7277 0.5145 0.2932

240 0.9298 0.7190 0.5044 0.2823

270 0.9124 0.6932 0.4944 0.2732

300 0.9097 0.6854 0.4797 0.2617

330 0.8934 0.6724 0.4653 0.2476

360 0.8862 0.6656 0.4522 0.2332

Anas Sani Maihulla, Micheal Khoo Boon Chong and Ibrahim Yusuf RELIABILITY ANALYSES OF AN INDUSTRIAL SYSTEM BASED ON HERMITE POLYNOMIAL AND 2-PARAMETER WEIBULL

Figure 1: Reliability analyses For failure against repair rate of a system using Hermite Polynomial

Figure 2: Reliability analyses For failure against repair rate of a system using 2-parameter Weibull distribution

IV. Discussion

Using the Hermite polynomial method, the result shown in table 1 that corresponds to Figure 1 demonstrated how the reliability of an industrial system is impacted by the rise in failure rate. The reliability of an industrial system was then shown to be affected as the failure rate was changed from 0.75 to 0.95 in table 2 that corresponds to Figure 2. When comparing the results in figures 1 and 2, it is obvious that the 2-parameter Weibull yields a more accurate result than the Hermite polynomials in terms of reliability estimation. Therefore, it is advised that future studies modify the Hermite polynomial to increase the dependability of industrial systems.

Acknowledgement

I recommended the effort by Tertiary Education Trust Fund (TETFund) for their kindness to sponsorship for Ph.D research as well as the benchwork here in Malaysia.

Anas Sani Maihulla, Micheal Khoo Boon Chong and Ibrahim Yusuf RELIABILITY ANALYSES OF AN INDUSTRIAL SYSTEM BASED ON HERMITE POLYNOMIAL AND 2-PARAMETER WEIBULL

References

[1] Bucher, C.G. Bourgund, U. (1990). A fast and efficient response surface approach for structural reliability problems. Struct Saf 7:57-66

[2] Szego, G. (1975). Orthogonal polynomials. 4th ed., Amer. Math. Soc, Providence, RI

[3] Yan, A. (1994). Linear stochastic differential equations and Wick products, Probability. Theory Rei. Fields, Vol. 99, 501-526

[4] Abramowitz, M. and I. S. Stegun, I. S. (1965). Handbook of Mathematical Functions." Dover Publications, Inc., New York, NY.

[5] Carpinteri, A. and Mainardi, F. (1997). Fractals and Fractional Calculus in Continium Mechanics. Springer-Verlag Wien New-York.

[6] Doha, E.H. Bhrawy, A. H. and Ezz-Eldien, S.S. (2013). A New Jacobi Operational Matrix: An application for Solving Fractional Differential Equations. Appl. Math. Modell, 36:4931-4943. https://doi.org/10.10167j.apm.2011.12.031.

[7] Bhrawy A. H. and Alofi, A. S. (2013). The Operational Matrix of Fractional Integration for Shifted Chebyshev Polynomials. Applied Mathematics Letters, 26:25-31. https://doi.org/10.10167j.aml.2012.01.027.

[8] Guo, B. Y. and Zhang, X. Y. (2005). A New Generalized Laugerre Approximation and Its Applications. J. Of Comput. App. Math, 181:342-363. https://doi.org/10.1016/j.cam.2004.12.008.

[9] Mikhailenko, B. G. (1999). Spectral Laugerre Method for the Approximate Solution of time dependent problems. Applied Mathematics Letters, 12:105-110. https://doi.org/10.1016/S0893-9659(99)00043-9.

[10] Bhrawy, A. H. Alghamdi M. M. and Taha, M. T. (2012). A New Modified Generalized Laugerre Operational Matrix of Fractional Integration for Solving Fractional Differential Equations on the Half Line Vol. 179. https://doi.org/10.1186/1687-1847- 2012-179.

[11] A.D Poularikas, , (1999) "The Handbook of Formulas and Tables for Signal Processing, Hermite Polynomials."

[12] Quek, S. T. and Ang A. H. S. (1986). Structural System Reliability by The Method of Stable Configuration. National Science Foundation, University of Illinois , ISSN: 0069-4274, NO. 529, Washington.

[13] Ulrich, L. Rohde, G. C. Jain, A. K. Poddar, G. (2011). Introduction to Differential Calculus. Systematic Studies with Engineering Applications for Beginners ISBN:9781118117750 I D0I:10.1002/9781118130155

[14] Anas, S. M. Ibrahim Y. Saminu I. B. (2023). WEIBULL COMPARISON BASED ON RELIABILITY, AVAILABILITY, MAINTAINABILITY, AND DEPENDABILITY (RAMD) ANALYSIS. Reliability: Theory & Applications. March 1(72):120-132. https://doi.org/10.24412/1932-2321- 2023-172-120-132

[15] Maihulla, A. S. and Yusuf, I. (2022). RELIABILITY ANALYSIS OF REVERSE OSMOSIS FILTRATION SYSTEM USING COPULA. Reliability: Theory & Applications, 17(2):

163-177. Retrieved from http://www.gnedenko.net/RTA/index.php/rta/article/view/890.

[16] Maihulla, A. S. Yusuf, I. and Bala, S. I. (2021). Performance Evaluation of a Complex Reverse Osmosis Machine System in Water Purification using Reliability, Availability, Maintainability and Dependability Analysis. Reliability: Theory & Applications, 16(3):115-131. Retrieved from http://www.gnedenko.net/RTA/index.php/rta/article/view/764

[17] Maihulla, A. Yusuf, I. (2022). Reliability modeling and performance analysis of reverse osmosis machine in water purification using Gumbel-Hougaard family copula," Life Cycle Reliab Saf Eng . https://doi.org/10.1007/s41872-022-00214-2

[18] Maihulla, A. S. and Yusuf, I. (2021). Performance Analysis of Photovoltaic Systems Using (RAMD) Analysis. Journal of the Nigerian Society of Physical Sciences, 3(3): 172-180. https://doi.org/10.46481/jnsps.2021.194