Next Article in Journal
Chirped-Pulse Amplification in an Echo-Enabled Harmonic-Generation Free-Electron Laser
Next Article in Special Issue
Quantitative Analysis of the Complex Time Evolution of a Camphor Boat
Previous Article in Journal
Study on Pavement Performance of Recycled Asphalt Mixture Modified by Carbon Nanotubes and Waste Engine Oil
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reliability by Using Weibull Distribution Based on Vibration Fatigue Damage

by
Jesús M. Barraza-Contreras
,
Manuel R. Piña-Monarrez
* and
Roberto C. Torres-Villaseñor
Industrial and Manufacturing Department, Engineering and Technological Institute, Universidad Autónoma de Ciudad Juárez, Ciudad Juárez 32310, Chihuahua, Mexico
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(18), 10291; https://doi.org/10.3390/app131810291
Submission received: 7 August 2023 / Revised: 8 September 2023 / Accepted: 11 September 2023 / Published: 14 September 2023

Abstract

:
In this paper, a Weibull probabilistic methodology is proposed with an approach to model vibration fatigue damage accumulation using two parameters: Weibull distribution and a nonlinear fatigue damage accumulation model. The damage is cumulated based on the application of a vibration stress profile and is used to determine both the Weibull β and η parameters, and the corresponding component reliability R(t). The vibration fatigue damage is analyzed to accumulate the damage as a stress function for a fatigue life exponent derived with the assistance of the acceleration’s force response. The steps to determine the Weibull β and η parameters are estimated based only on the principal vibration stresses σ 1 and σ 2 that allow the reproduction of the vibration fatigue damage. The method’s efficiency is based on the probabilistic approach by using the vibration fatigue damage as the Yi vector that covers the arithmetic mean as well as the β parameter. Finally, the procedure proposed is applied in a practical case where a mechanical component is used as a support for telecommunication connections and is submitted to vibration stress. The results show that using the damage accumulated as the Yi vector to estimate the parameters allows for the analysis of dynamic and individual applications.

1. Introduction

Mechanical components and structures subjected to vibration are affected by dynamic loads, which induce fatigue damage due to the cycling loading application [1,2]. The generated fatigue damage is primarily related to vibration history loading, geometry, and material properties [3]. The generated fatigue damage directly determines the component’s reliability, replacement policy, and warranty costs. Thus, the reliability index characteristic has an important role that is determined mathematically and is used to describe the fatigue damage behavior of mechanical components [4]. That description can be performed by using a probabilistic approach [5]. Then, for mechanical components or systems, according to their application and data complexity, an accurate method must be selected to determine an effective level of service life reliability [6]. Now, since data from fatigue damage are affected by significant scatter, and damage is provoked as a response to random forces, gathering fatigue damage data is generally a difficult activity [7]. One of the more commonly applied methods to work with fatigue damage is the Miner’s rule [8,9]. To consider the random nature of the generated damage, here, we use a probabilistic time-dependent approach [7]. In the fatigue damage accumulation models, the principal factors involved are load sequence, type of load, overloads, plasticization, and type of material. Consequently, for the damage accumulation analysis, we require a probabilistic concept, or a physical quantity related to the probability of occurrence [10]. Thus, the measurement of damage helps us to calculate the probabilities of failure. Since the accumulation of random vibration fatigue damage entails increasing deterioration, an increasing hazard function is required, and the most recommendable cumulative distribution function (cdf) used to estimate the fatigue damage is the Weibull cdf [11,12,13]. In a vibration profile, each row has its own stationary frequency and amplitude which allow for the performance of a vibration analysis with a statistical approach. Since it is possible to estimate the stress from the acceleration responses, here, the Weibull distribution is used. Thus, based on the increasing behavior of the cumulated damage, in this paper, the increasing random damage behavior is used in the Weibull cdf to determine the failure percentile that the observed cumulated damage represents in the used vibration profile. Then, once the damage percentiles are determined based on their corresponding cycles of the S–N curve, the Weibull scale parameter is determined. Similarly, the Weibull shape parameter is determined directly from the cumulated response stress of the used profile (see Equations (9) and (10) and Section 3.2). Thus, in Section 2.3, a probabilistic methodology to characterize the fatigue damage induced by random vibration is developed by using the Weibull distribution, which uses a relation of the scale and shape parameters with the mechanical vibration fatigue damage. The methodology includes probabilistic estimations based on the Weibull scale and shape parameters that are governed by a new way of analyzing the fatigue damage accumulation presented in Section 2.1. Then, the vibration analysis is based on the fatigue damage, which is determined directly from the principal stresses σ 1 and σ 2 ; therefore, the methodology is efficient because it uses the damage as the platform to project and represent its random behavior and consequently the component’s fatigue life. With the purpose to assess it, the methodology is applied in Section 3 to a probabilistic failure analysis of a panel support made of cold drawn steel AISI 1025. The mechanical component is submitted to a random vibration loading profile of 10–55 Hertz with an amplitude of 1.5 mm peak to peak during a 2 h period per block. The testing was performed by using an electrodynamic vibration system in which the vibration profile loading was applied 29 times. The component’s physical damage results are shown in Section 3.1. The purpose of the method proposed lies in the use of the accumulated fatigue damage Di instead of the median rank operation. This is illustrated by Equations (12) and (13) which allow the use of the resultant vector Y in the estimation of Weibull parameters that completely reproduces the principal vibration stress values.
The paper is organized as follows. Section 2 includes the generalities of the vibration fatigue damage accumulation and the proposed Weibull fatigue damage analysis method. In Section 3, a numerical application is presented. Section 4 is related to the median rank method comparison. Finally, in Section 5, the conclusions are given.

2. Fatigue Damage

2.1. Fatigue Damage Accumulation

Fatigue damage can be described as a failure mechanism that is manifested when a material tends to fail or break under repeated deflections [14,15,16]. Thus, a nonlinear model to accumulate the random vibration fatigue damage has been proposed [17] with the purpose of evaluating the fatigue damage of different dynamic loads in mechanical components and structural elements. The acceleration response of the analyzed vibration system is determined by stress as in Equation (1). The applied vibration cycles are determined by the rainflow method [18]. From the S–N material’s curve, the corresponding life cycles are determined by using the Basquin Equation [19] as is in Equation (4).
σ v i b i = σ d y n a m i c × A r e s
σ d y n a m i c = K m e L ^ C I A
A r e s = 2 π 2 F 2 D 2 G
N i × σ v i b i b = a b
In Equation (1), σ d y n a m i c and A r e s are the dynamic load factor and the acceleration response, respectively. In Equation (2) [20], K is the stress concentration factor in the mechanical component, m e is the effective mass, C is the distance to the neutral axis, L ^ is the distance from the fixed point of the component to the point of application load, A is the constant of gravity, and I is the moment of inertia. In Equation (3), F is the frequency applied by the vibration power spectral density (PSD) and G is the gravity constant. In Equation (4), N i   is the maximum number of cycles that the material’s component can sustain at a vibration load with stress amplitude σ v i b i and the parameters a and b are constant variables that represent the intercept and the slope of the S–N curve, respectively.
As shown by Equations (1)–(4), because a vibration has a nonlinear behavior [21], the generated fatigue damage also presents a nonlinear behavior. Consequently, the fatigue damage is determined by using Equation (5), where the damage is described by a curve that represents the effect under two-level loading conditions [22], where n i represents the applied vibration cycles at the stress level σ v i b i .
D = i = 1 2 D 2 = n 2 N 2 , f N 2 , f N 1 , f σ ( v i b ) 1 σ ( v i b ) 2
Now that the random fatigue damage generated by the vibration environment is determined, a Weibull formulation is presented that let us use the generated damage in the Weibull Y vector (see Equation (13)) to determine the reliability of the analyzed element.

2.2. Weibull Analysis

A two-parameter Weibull distribution is used to statistically analyze fatigue behaviors [23,24,25]. It allowed us to perform accurate fatigue failure analysis [26,27]. The probability density function f(t) and cumulative distribution function F(t) are described by Equations (6) and (7), respectively.
f t = β η t η β 1 e x p t η β
F t = 1 e x p t η β
where, β is the shape parameter, η is the scale parameter, and t is the selected random variable (damage or fatigue life). The corresponding reliability function R(t) is given as
R t = e x p t η β
From [28], the Weibull fatigue damage β D and η D parameters are determined as
β D = 4 μ y 0.995 × l n σ 1 σ 2
η D = e x p ( µ x )  
where μ y represents the mean of the Y vector (see Equation (13)) determined by using Equation (5). µ x represents the log-mean of the failure-time data, which is determined here directly from the addressed maximum σ 1 and minimum σ 2 stress values of Section 2.1. Thus, µ x is determined as
µ x = l n σ 1 σ 2 1 2
Here, notice that the efficiency of the Weibull parameters β D and η D only depends on the accuracy with which the σ 1 and σ 2 values are determined by Equation (1). In this paper, they were determined from an electrodynamic shaker acceleration response and by the method performed in [17], Section 3.1.

2.3. Weibull Fatigue Damage Analysis

The expected behavior of the σ 1 and σ 2 parameters is determined by the following steps:
Step 1. By using the fatigue damage accumulation results from the component’s operation, the times that the vibration profile (PSD) loading were applied to the mechanical component are taken as n.
Step 2. The accumulation fatigue damage D i determined by Equation (5) in the end of each one of the n vibration loadings is used by Equation (12) as the corresponding cumulated failure percentile R D i , as
R ( D i ) = 1 D i
Step 3. By using the R D i   elements in the linearized form of the reliability function given in Equation (8), the corresponding Y i elements are determined as in Equation (13). Then, its corresponding arithmetic mean value is computed as in Equation (14).
Y i = L N ( L N ( 1 D i ) )
μ y = i = 1 n Y i n
Step 4. By plugging the μ y value and the σ 1 and σ 2 values into Equation (9), the corresponding Weibull shape β D parameter is determined. Similarly, by plugging the σ 1 and σ 2 values into Equation (11), the corresponding μ x value is determined. Then, by using this corresponding μ x value in Equation (10), the corresponding Weibull scale η D parameter is determined. These β D and η D parameters represent the Weibull fatigue damage family that is used to model the random behavior of the estimated σ 1 and σ 2 principal stress values.
Note 1. Here, notice the random behavior of the σ 1 and σ 2 values. In the proposed Weibull analysis, let us use the σ 1 i values as the minimum required strength that the component’s material must present in order to ensure that the reliability of the component will meet at least (as a minimum) the desired R(t) index.
From the Weibull analysis, by using the β D and η D parameters, the minimum strength σ 1 i values are determined by using the t 0 i value that correspond to each Y i element as
t 0 i = e x p { Y i / β D }
Thus, the σ 2 i value is determined as
σ 2 i = η D × t 0 i
and the σ 1 i value is determined as
σ 1 i = η D / t 0 i
Additionally, from Equation (18), by using the known σ 1 value, the t 01 element that belongs to the σ 1 and σ 2 values determined in Section 2.1 is determined as
t 01 = η D / σ 1
Now, the t 01 and the β D values are used to determine the corresponding Y 1 value as follows,
Y 1 = ln t 01 × β D
Finally, the reliability index that corresponds to the Y 1 value is determined as
R t = exp exp Y 1
Note 2. Here, observe that the R(t) index determined in Equation (20) by using the σ 1 value according to our proposed method corresponds to a component with strength equivalent to the σ 1 value. Thus, if we define the material S y parameter as the actual strength of the component, then by using this Sy value in Equation (18) and the corresponding Y S y value of Equation (19) in Equation (20), the minimum expected reliability of a component that presents a strength of S y is determined. Please also notice from the proposed Weibull analysis that any desired strength value can be used to determine its corresponding reliability. The diagram with the steps required to determine the probability of failure F t and the reliability R t based on the vibration fatigue damage D is shown in Figure 1.
Next, a numerical application is presented where the proposed method is applied to a mechanical component that is submitted to random vibration stress due to its field application.

3. Numerical Application

3.1. Fatigue Damage Accumulation

The accumulation fatigue damage study case is performed by analyzing a mechanical support component that is used to install a fiber optic panel into a frame. It is shown in Figure 2.
The support is made of cold drawn steel AISI 1025, with a modulus of elasticity E = 200 GPa, Poisson’s ratio ɤ = 0.29, yield strength S y   = 430 MPa, ultimate tensile strength Sut = 510 MPa, endurance limit Se = 255 MPa, density ρ   = 7.9 g/cm3, length L = 51 mm, width W = 200 mm, and a wall thickness t = 3 mm. During its function, the component supports a static load of 80 N. It is submitted to an operating random vibration with an input PSD consisting of frequencies ranging from 10 to 55 Hz at an amplitude of 1.5 mm peak to peak for a period of 2 h. The testing is carried out physically by using a vibration system, and the results are as follows. By using Equations (1) and (2), the dynamic factor σ d y n a m i c and the vibration stresses σ v i b i , are calculated, respectively. The acceleration responses A r e s i are obtained from the vibration system but can also be determined by using Equation (3). Then, the vibration cycles applied n i are determined by the rainflow method and the total cycles N i are determined by using Equation (4). Finally, the fatigue damage accumulation is obtained by Equation (5). Table 1 shows the vibration stress and cycle results and Table 2 shows the vibration fatigue damage accumulation results, where the failure is presented when D = 1 [29].
Figure 3 shows the mechanical support areas where the fatigue damage accumulated was presented. The analysis and estimation of the damage were determined from the acceleration responses that were the base to calculate the principal vibration stresses, σ 1 and σ 2 , that are employed in the Weibull probabilistic analysis.
Now, in the following section, the probabilistic approach is applied to the vibration fatigue damage results obtained in Table 2.

3.2. Weibull Fatigue Damage Analysis

Here, it is remarked that if a different sample with the same features was submitted to the same vibration damage accumulation experiment, due to random behavior, the results would be different. For the damage data given in Table 2, where 29 blocks were tested, the probabilistic analysis is as follows. For this purpose, the vector Y that includes the fatigue damage as described by Equation (13) is used, instead of using the median rank approach, to determine the Weibull parameters as follows [30,31].
By selecting the fatigue damage accumulated Di value of each individual block of Table 2, and using it in Equation (13), the n = 29 Y i elements are determined. Thus, from Equation (14), the mean is, μ y = 0.6672 . Next, from the results of Table 1, the principal vibration stress values σ 1 = 304.76 MPa and σ 2 = 15.99 MPa are obtained. With these data, we can proceed to use the Weibull distribution in order to obtain its failure probability, reliability, and random behavior. From Equation (9), the Weibull shape parameter β D value is
β D = ( 4 ) ( 0.6672 ) 0.995 × l n ( 304.76 15.99 ) = 0.9102
and from Equation (11), the logarithm average value is μ x = ln 304.76 × 15.99 = 4.2457 . Thus, from Equation (10), the Weibull parameter ƞ D is ƞ D = e x p { 4.2457 } = 69.8077 . Consequently, the Weibull damage family is W(0.9102, 69.8077 MPa). This Weibull family completely represents the observed principal vibration stresses σ 1 and σ 2 values.
Now, by using the Weibull family results, the random behavior strength can be determined. Since the Weibull parameters only depend on the principal stress values caused by the random vibration σ 1 and σ 2 values, the random behavior can be obtained by performing the Weibull analysis using the following steps.
Since the determination of the fatigue damage is based on the random behavior of the σ2 stress value, here, the random behavior of σ 1   and σ 2 is determined by Equations (16) and (17). Then, by using the β D and the Y i values in Equation (15), the basic Weibull elements [28] t 0 i for each Y i are obtained. Whereas by using the η D and σ 1 values in Equation (18), the Weibull t 01   value from the σ 1 and σ 2 stress values are reproduced. This value is calculated as
t 01 = 69.8077 304.76 = 0.2291
and, by using the β D value in Equation (19), the Y 1 value that belongs to the t 01 value is determined as
Y 1 = ln 0.2291 × 0.9102 = 1.3414
Next, by substituting the Y 1 value in Equation (20), the reliability R t   that belongs to the t 01 element is
R t = exp exp 1.3414 = 0.7699
The previous results shown in this section are included in Table 3.
Here, is important to mention that the reliability obtained R t = 0.7699 corresponds to a designed component or structure with minimal strength of S y = σ 1 = 304.76 MPa. In relation to the mechanical support, it has a S y = 430 MPa, then, its reliability is R t = 0.8260 , which is determined as follows. The minimal reliability of the component is obtained when the Y 1 value that belongs to the S y value is used in Equation (20). The steps to determine the reliability of the design component when the S y value is used as σ 1 in Equation (18) are, the t S y element that belongs to the S y value is t s y = 69.8077 430 = 0.1623 . From Equation (19), the corresponding Y s y value is Y s y = ln 0.1623 × 0.9102 = 1.6548 . The reliability index for the Y s y value is calculated by using Equation (20), R t = exp exp 1.6548 = 0.8260 . Thus, we conclude that the reliability of the design component is R t = 0.8260 . Additionally, regarding the material’s component or structure, it is noticed about that the higher the strength Sy value, the higher the reliability R t will be [28].
Now, as a comparison for the probabilistic cumulative density, the Weibull distribution is used, but in this case, by using the median rank method to estimate the vector Y , which is the base to determine the reliability and the probability of failure.

4. Median Rank Approach

In this section, the median rank method [30] is applied to the data given in Section 3. By using this method, the corresponding cumulated failure percentile F t i   that previously was determined by Equation (12) is now determined as
F t i = i 0.3 n + 0.4
Then, by using the F t i elements in the linearized form of the reliability function given in Equation (8), the corresponding Y i elements are determined as in Equation (22),
Y i = L N ( L N ( 1 ( ( i 0.3 ) / ( n + 0.4 ) ) ) )
By selecting the fatigue damage accumulated Di value of each one block of Table 2 and using it in Equation (22), the n = 29 Y i elements are determined. Thus, from Equation (14), the mean is, μ y = 0.5525 . Next, from the results of Table 1, the principal vibration stress values σ 1 = 304.76 MPa and σ 2 = 15.99 MPa are obtained. With these data, we can now proceed to use the Weibull distribution in order to obtain its failure probability, reliability, and random behavior. From Equation (9), the Weibull shape parameter β D value is
β D = ( 4 ) ( 0.5525 ) 0.995 × l n ( 304.76 15.99 ) = 0.7538
Since the principal vibration stress values σ1 = 304.76 MPa and σ2 = 15.99 MPa are maintained, the Weibull parameter value η D = 69.8077 remains.
Next, the remaining steps followed in Section 3 are applied and the results are included in Table 4.
A comparison of the principal results between both the proposed method and the median rank method by following the formulation given in Section 2 is shown in Table 5.
In Table 5, notice that although we are using the same σ 1 = 304.76, σ 2   = 15.99 and S y values in both methods, the component’s reliability R(t) is different. This difference occurs because in both methods, the Weibull β parameter shape have different values. This is because the arithmetic mean of both methods is different. Thus, because in the proposed method the addressed damage completely represents the analyzed component, we conclude that the real reliability is the one given by the proposed method. Notic that because the damage is random, the proposed method is dynamic, and its efficiency depends only on the accuracy with which σ 1 and σ 2   are determined. Finally, observe from the last row of Table 5 that if R(t) = 0.95 is required, then from the proposed method, the maximum allowed stress is σ 1 = 1824.3190 MPa. At the same time, using the median rank method, it is σ 1 = 3591.0886 MPa, implying that by using the median rank approach, we can predict that the component will become overstressed, lowering its life.

5. Discussion

The paper presents an alternative to fit a probabilistic vibration fatigue analysis based on the use of the damage accumulated for two parameters provided by the Weibull distribution model. The main contribution of the proposed method is the probabilistic approach that the Weibull distribution function has in the mechanical industrial field. A considerable number of industrial standards and guidelines employ the Weibull model for their fatigue analysis; thus, reliability engineers are especially familiar with this model.
This paper presents a probabilistic aspect that involves fatigue experiments and the proper use of vibration fatigue damage during the investigation and design of mechanical components. The use of vibration fatigue damage as the cumulative failure percentile in the process allows the transfer to obtain an estimation of the reliability and a prediction of failure, in spite of the significant variability involved in the vibration fatigue.
Regarding the results obtained for the probabilistic prediction of reliability and cumulated failure, based on the principal vibration stress values σ 1 = 304.76 and σ 2 = 15.99, the proposed method predicts a reliability of R(t) = 0.826 and a cumulative failure of F(t) = 0.174. Here, it is important to notice that if we use the standard median rank approach, then R(t) = 0.776 and F(t) = 0.224. By their comparison, we have an R(t) variation of 5%. Moreover, it can be observed that because the used cumulated damage depends, among other factors, on the vibration load distribution applied, the material’s strength, and the geometry of the tested mechanical component, then because each analysis is different, the proposed method is completely dynamic, and can be used in any application where the accumulated damage can be measured. Here, we highlight that the only restriction to use the cumulated damage is that the failure must be defined when the damage is equal to one. Therefore, among other applications, the proposed model can be used in the automotive and telecommunications industries, where mechanical supports are used to hold electrical and electronic devices that are subjected to environmental vibration due to their field application.

6. Conclusions

  • A probabilistic alternative of the Weibull distribution to vibration fatigue analysis is developed, which allows it to define reliability and probability of failure.
  • Contrary to other models that use the median rank method as the cumulated failure percentile, the proposed methodology considers the accumulated fatigue damage for a platform in which the component’s probabilistic life, the probability of failure, and the reliability index are estimated.
  • A methodology is developed based on the model presented to permit a probabilistic approach to vibration fatigue damage accumulation which allows the probabilistic failure and reliability estimation for mechanical components and structures subjected to variable amplitude loading, specifically random vibration.
  • The model and methodology included in this work are applied to mechanical components used in the telecommunication industry to assist in real and practical structural fatigue analysis.
  • An application case is selected to illustrate the proposed Weibull model and its parameter estimation methodology based on the cumulated vibration damage included in this paper.

Author Contributions

Conceptualization, J.M.B.-C., M.R.P.-M. and R.C.T.-V.; methodology, J.M.B.-C. and M.R.P.-M.; data analysis, J.M.B.-C. and R.C.T.-V.; writing—original draft preparation, J.M.B.-C. and M.R.P.-M.; writing—review and editing, J.M.B.-C., M.R.P.-M. and R.C.T.-V.; supervision, M.R.P.-M.; funding acquisition, J.M.B.-C. and R.C.T.-V. All authors have read and agreed to the published version of the manuscript.

Funding

The research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sosnovskiy, L. Tribo Fatigue; Springer: New York, NY, USA, 2005. [Google Scholar]
  2. Piña-Monarrez, M.R. Weibull analysis for normal/accelerated and fatigue random vibration test. Qual. Reliab. Eng. Int. 2019, 35, 2408–2428. [Google Scholar] [CrossRef]
  3. Mršnik, M.; Slavič, J.; Boltežar, M. Frequency-domain methods for a vibration-fatigue-life estimation—Application to real data. Int. J. Fatigue 2013, 47, 8–17. [Google Scholar] [CrossRef]
  4. Zaitseva, E.; Levashenko, V.; Rabcan, J. A new method for analysis of Multi-State systems based on Multi-valued decision diagram under epistemic uncertainty. Reliab. Eng. Syst. Saf. 2023, 229, 108868. [Google Scholar] [CrossRef]
  5. Isametova, M.E.; Nussipali, R.; Martyushev, N.V.; Malozyomov, B.V.; Efremenkov, E.A.; Isametov, A. Mathematical Modeling of the Reliability of Polymer Composite Materials. Mathematics 2022, 10, 3978. [Google Scholar] [CrossRef]
  6. Martyushev, N.V.; Malozyomov, B.V.; Sorokova, S.N.; Efremenkov, E.A.; Valuev, D.V.; Qi, M. Review Models and Methods for Determining and Predicting the Reliability of Technical Systems and Transport. Mathematics 2023, 11, 3317. [Google Scholar] [CrossRef]
  7. Benaroya, H.; Han, S.M. Probability Models in Engineering and Science; Taylor & Francis Group, LLC: Boca Raton, FL, USA, 2005. [Google Scholar]
  8. Shimokawa, T.; Tanaka, S. A statistical consideration of Miner’s rule. Int. J. Fatigue 1980, 2, 165–170. [Google Scholar] [CrossRef]
  9. Sun, Q.; Dui, H.N.; Fan, X.L. A statistically consistent fatigue damage model based on Miner’s rule. Int. J. Fatigue 2014, 69, 16–21. [Google Scholar] [CrossRef]
  10. Castillo, E.; Fernandez-Canteli, A. A Unified Statistical Methodology for Modeling Fatigue Damage; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  11. Harlow, D.G. Statistical Modeling for Fatigue Life. Metals 2023, 13, 1419. [Google Scholar] [CrossRef]
  12. Pinto, J.M.A.; Pujol, J.C.F.; Cimini, C.A. Probabilistic cumulative damage model to estimate fatigue life. Fatigue Fract. Eng. Mater. Struct. 2014, 37, 85–94. [Google Scholar] [CrossRef]
  13. Rinne, H. The Weibull Distribution: A Handbook; Taylor & Francis Group, LLC: Boca Raton, FL, USA, 2009. [Google Scholar]
  14. Batsoulas, N.D. Cumulative Fatigue Damage: CDM-Based Engineering Rule and Life Prediction Aspect. Steel Res. Int. 2016, 87, 1670–1677. [Google Scholar] [CrossRef]
  15. Abdullah, S.; Beden, S.M.; Ariffin, A.K.; Rahman, M.M. Fatigue Life Assessment for Metallic Structure: A Case Study of Shell Structure under Variable Amplitude Loading. J. Appl. Sci. 2008, 8, 1622–1631. [Google Scholar] [CrossRef]
  16. Burhan, I.; Kim, H. S-N Curve Models for Composite Materials Characterization: An Evaluative Review. J. Compos. Sci. 2018, 2, 38. [Google Scholar] [CrossRef]
  17. Barraza-Contreras, J.M.; Piña-Monarrez, M.R.; Molina, A.; Torres-Villaseñor, R.C. Random Vibration Fatigue Analysis Using a Nonlinear Cumulative Damage Model. Appl. Sci. 2022, 12, 4310. [Google Scholar] [CrossRef]
  18. ASTM E1049-85; ASTM International Standard Practices for Cycle Counting in Fatigue Analysis. ASTM International Standard: West Conshohocken, PA, USA, 2017.
  19. Usabiaga, H.; Muniz-Calvente, M.; Ramalle, M.; Urresti, I.; Fernández Canteli, A. Improving with probabilistic and scale features the Basquin linear and bi-linear fatigue models. Eng. Fail. Anal. 2020, 116, 104728. [Google Scholar] [CrossRef]
  20. Irvine, T. A Fatigue Damage Spectrum Method for Comparing Power Spectral Density Base Input Specifications. Vibrationdata. 2014. Available online: https://vibrationdata.wordpress.com/ (accessed on 28 July 2023).
  21. Harris, C.M.; Piersol, A.G. Shock and Vibration Handbook, 5th ed.; McGraw Hill: New York, NY, USA, 2002; Volume 15, No. 9. [Google Scholar] [CrossRef]
  22. Manson, S.S.; Halford, G.R. Practical implementation of the double linear damage rule and damage curve approach for treating cumulative fatigue damage. Int. J. Fract. 1981, 17, 169–192. [Google Scholar] [CrossRef]
  23. Weibull, W. A Statistical Theory of the Strength of Materials; Generalstabens Litografiska Anstalts Förlag: Stockholm, Sweden, 1939. [Google Scholar]
  24. Harikrishnan, R.; Mohite, P.M.; Upadhyay, C.S. Generalized Weibull model-based statistical tensile strength of carbon fibres. Arch. Appl. Mech. 2018, 88, 1617–1636. [Google Scholar] [CrossRef]
  25. Melikayeva, I.; Riesch-Oppermann, H.; Govorukha, V.B.; Kraft, O. A generalized Weibull approach to interface failure in bi-material ceramic joints. Arch. Appl. Mech. 2011, 81, 1585–1596. [Google Scholar] [CrossRef]
  26. Choi, S.-K.; Grandhi, R.V.; Canfield, R.A. Reliability-Based Structural Design; Springer-Verlag: London, UK, 2006. [Google Scholar]
  27. Barraza-Contreras, J.M.; Piña-Monarrez, M.R.; Molina, A. Fatigue-life prediction of mechanical element by using the Weibull distribution. Appl. Sci. 2020, 10, 6384. [Google Scholar] [CrossRef]
  28. Piña-Monarrez, M.R. Weibull stress distribution for static mechanical stress and its stress/strength analysis. Qual. Reliab. Eng. Int. 2018, 34, 229–244. [Google Scholar] [CrossRef]
  29. Lee, Y.-L.; Pan, J.; Hathaway, R.B.; Barkey, M.E. Fatigue Testing and Analysis: Theory and Practice; Elsevier Butter-worth-Heinemann: Burlington, MA, USA, 2005. [Google Scholar]
  30. Mischke, C.R. Predicting Fatigue Endurance Strengths for High- and Low-Cycle Fatigue. J. Mech. Des. 1982, 104, 653–660. [Google Scholar] [CrossRef]
  31. Dias, J.P.; Ekwaro-Osire, S.; Cunha, A., Jr.; Dabetwar, S.; Nispel, A.; Alemayehu, F.M.; Endeshaw, H.B. Parametric probabilistic approach for cumulative fatigue damage using double linear damage rule considering limited data. Int. J. Fatigue 2019, 127, 246–258. [Google Scholar] [CrossRef]
Figure 1. Diagram of the Weibull Fatigue Damage Analysis.
Figure 1. Diagram of the Weibull Fatigue Damage Analysis.
Applsci 13 10291 g001
Figure 2. Mechanical support.
Figure 2. Mechanical support.
Applsci 13 10291 g002
Figure 3. Fatigue cracks on the support after vibration testing.
Figure 3. Fatigue cracks on the support after vibration testing.
Applsci 13 10291 g003
Table 1. Vibration Stress and Cycles Results.
Table 1. Vibration Stress and Cycles Results.
Frequency
(Hz)
Accel.
Response (G)
Dynamic Factor
σ d y n a m i c
Equation (2)
Vibration Stress
σ v i b i
Equation (1)
Applied
Vibration Cycles ( n i )
Total Cycles
( N i )
Equation (4)
100.7222.2215.9970,3841.36 × 1020
202.6558.86140,1953.04 × 1013
305.62124.8492,6194.42 × 109
409.17203.6910,8071.40 × 107
5013.72304.7629211.23 × 105
5512.36274.557624.20 × 105
Table 2. Vibration Fatigue Damage Accumulation.
Table 2. Vibration Fatigue Damage Accumulation.
10 Hz20 Hz30 Hz40 Hz50 Hz55 Hz
Block No.D1D1+2D1+2+3D1+2+3+4D1+2+3+4+5D1+2+3+4+5+6
1 (2 h)5.18 × 10−165.18 × 10−165.20 × 10−165.73 × 10−162.37 × 10−22.00 × 10−2
2 (2 h)2.00 × 10−22.42 × 10−22.42 × 10−22.42 × 10−24.84 × 10−25.00 × 10−2
3 (2 h)5.00 × 10−24.91 × 10−24.93 × 10−25.01 × 10−27.38 × 10−27.00 × 10−2
4 (2 h)7.00 × 10−27.49 × 10−27.51 × 10−27.63 × 10−21.00 × 10−11.00 × 10−1
5 (2 h)1.00 × 10−11.01 × 10−11.02 × 10−11.03 × 10−11.27 × 10−11.30 × 10−1
6 (2 h)1.30 × 10−11.29 × 10−11.29 × 10−11.31 × 10−11.55 × 10−11.60 × 10−1
7 (2 h)1.60 × 10−11.56 × 10−11.57 × 10−11.59 × 10−11.83 × 10−11.90 × 10−1
8 (2 h)1.90 × 10−11.85 × 10−11.86 × 10−11.89 × 10−12.12 × 10−12.10 × 10−1
9 (2 h)2.10 × 10−12.14 × 10−12.15 × 10−12.18 × 10−12.42 × 10−12.40 × 10−1
10 (2 h)2.40 × 10−12.45 × 10−12.45 × 10−12.49 × 10−12.73 × 10−12.80 × 10−1
11 (2 h)2.80 × 10−12.76 × 10−12.76 × 10−12.81 × 10−13.04 × 10−13.10 × 10−1
12 (2 h)3.10 × 10−13.07 × 10−13.08 × 10−13.13 × 10−13.37 × 10−13.40 × 10−1
13 (2 h)3.40 × 10−13.40 × 10−13.41 × 10−13.46 × 10−13.70 × 10−13.70 × 10−1
14 (2 h)3.70 × 10−13.73 × 10−13.74 × 10−13.80 × 10−14.03 × 10−14.10 × 10−1
15 (2 h)4.10 × 10−14.07 × 10−14.08 × 10−14.14 × 10−14.38 × 10−14.40 × 10−1
16 (2 h)4.40 × 10−14.42 × 10−14.43 × 10−14.50 × 10−14.73 × 10−14.80 × 10−1
17 (2 h)4.80 × 10−14.77 × 10−14.79 × 10−14.86 × 10−15.10 × 10−15.10 × 10−1
18 (2 h)5.10 × 10−15.14 × 10−15.15 × 10−15.23 × 10−15.47 × 10−15.50 × 10−1
19 (2 h)5.50 × 10−15.51 × 10−15.53 × 10−15.61 × 10−15.85 × 10−15.90 × 10−1
20 (2 h)5.90 × 10−15.89 × 10−15.91 × 10−16.00 × 10−16.24 × 10−16.30 × 10−1
21 (2 h)6.30 × 10−16.29 × 10−16.30 × 10−16.40 × 10−16.63 × 10−16.70 × 10−1
22 (2 h)6.70 × 10−16.68 × 10−16.70 × 10−16.80 × 10−17.04 × 10−17.10 × 10−1
23 (2 h)7.10 × 10−17.09 × 10−17.11 × 10−17.22 × 10−17.45 × 10−17.50 × 10−1
24 (2 h)7.50 × 10−17.51 × 10−17.53 × 10−17.64 × 10−17.88 × 10−17.90 × 10−1
25 (2 h)7.90 × 10−17.94 × 10−17.96 × 10−18.08 × 10−18.31 × 10−18.40 × 10−1
26 (2 h)8.40 × 10−18.37 × 10−18.40 × 10−18.52 × 10−18.76 × 10−18.80 × 10−1
27 (2 h)8.80 × 10−18.82 × 10−18.85 × 10−18.97 × 10−19.21 × 10−19.30 × 10−1
28 (2 h)9.30 × 10−19.28 × 10−19.30 × 10−19.44 × 10−19.67 × 10−19.70 × 10−1
29 (2 h)9.70 × 10−19.74 × 10−19.77 × 10−19.91 × 10−11.01 × 1001.02 × 100
Table 3. Weibull vibration fatigue damage statistics analysis for the numerical application data.
Table 3. Weibull vibration fatigue damage statistics analysis for the numerical application data.
niDamage (Di)
Equation (5)
R (Di)
Equation (12)
Yi
Equation (13)
µy
Equation (14)
R(t)
Equation (20)
toi
Equation (15)
σ2i
Equation (16)
σ1i
Equation (17)
F(t)
Equation (7)
10.02420.9758−3.7106−0.12800.97580.01701.18434114.91150.0242
20.04910.9509−2.9881−0.10300.95090.03752.61921860.54200.0491
0.05000.9500−2.9702−0.10240.95000.03832.67121824.31900.0500
30.07490.9251−2.5535−0.08810.92510.06054.22191154.23630.0749
40.10130.8987−2.2366−0.07710.89870.08575.9806814.82660.1013
50.12850.8715−1.9838−0.06840.87150.11317.8951617.22920.1285
0.17400.8260−1.6548−1.81800.82600.162311.3328430.00000.1740
60.15640.8436−1.7713−0.06110.84360.14289.9713488.71230.1564
70.18510.8149−1.5863−0.05470.81490.175012.2184398.83350.1851
80.21450.7855−1.4212−0.04900.78550.209814.6488332.66370.2145
0.23010.7699−1.3414−1.47380.76990.229115.9900304.76000.2301
90.24460.7554−1.2709−0.04380.75540.247517.2778282.04440.2446
100.27560.7244−1.1321−0.03900.72440.288320.1243242.15050.2756
110.30720.6928−1.0023−0.03460.69280.332523.2107209.95070.3072
120.33970.6603−0.8794−0.03030.66030.380526.5641183.44710.3397
130.37290.6271−0.7621−0.02630.62710.432930.2170161.27040.3729
140.40690.5931−0.6492−0.02240.59310.490034.2089142.45150.4069
150.44180.5582−0.5396−0.01860.55820.552838.5881126.28540.4418
160.47740.5226−0.4323−0.01490.52260.621943.4144112.24650.4774
170.51390.4861−0.3266−0.01130.48610.698548.763099.93460.5139
180.55130.4487−0.2215−0.00760.44870.784054.730189.03890.5513
190.58950.4105−0.1162−0.00400.41050.880261.441479.31320.5895
200.62850.3715−0.0097−0.00030.37150.989469.064970.55850.6285
0.63210.36790.00000.00000.36791.000069.806569.80650.6321
210.66850.33150.09900.00340.33151.115077.832762.61010.6685
220.70940.29060.21160.00730.29061.261788.078155.32720.7094
230.75110.24890.32990.01140.24891.4368100.303248.58380.7511
240.79380.20620.45690.01580.20621.6519115.316142.25880.7938
250.83750.16250.59720.02060.16251.9273134.542636.21990.8375
260.88210.11790.75990.02620.11792.3045160.871130.29200.8821
270.92770.07230.96590.03330.07232.8898201.726924.15700.9277
280.97430.02571.29790.04480.02574.1616290.514516.77410.9743
0.97820.02181.34141.47380.02184.3657304.760015.99000.9782
290.99000.01001.52720.05270.01005.3540373.749813.03840.9900
β = 0.9102ƞ = 69.8065µy = −0.6672σ1 = 304.7600σ2 = 15.9900
BOLD: The principal vibration stresses above and below the ƞ parameter, the R(t) index of 95%, and the R(t) index of the Sy 430 MPa.
Table 4. Weibull vibration fatigue damage statistics analysis by using median rank method.
Table 4. Weibull vibration fatigue damage statistics analysis by using median rank method.
niYi
Equation (21)
µy
Equation (14)
R(t)
Equation (20)
toi
Equation (15)
σ2i
Equation (16)
σ1i
Equation (17)
F(t)
Equation (7)
1−3.7256−0.12850.97620.00710.49819783.40300.0238
−2.9702−0.10240.95000.01941.35703591.08860.0500
2−2.8207−0.09730.94220.02371.65462945.17600.0578
3−2.3400−0.08070.90820.04493.13111556.35560.0918
4−2.0062−0.06920.87410.06984.8755999.50450.1259
5−1.7476−0.06030.84010.09846.8706709.26680.1599
6−1.5347−0.05290.80610.13059.1130534.74380.1939
−1.3704−1.81800.77570.162311.3328430.00000.2243
7−1.3524−0.04660.77210.166311.6070419.84080.2279
8−1.1918−0.04110.73810.205814.3630339.28230.2619
−1.1109−1.47380.71940.229115.9900304.75000.2806
9−1.0474−0.03610.70410.249217.3959280.12930.2959
10−0.9154−0.03160.67010.296920.7257235.12380.3299
11−0.7930−0.02730.63610.349224.3773199.90340.3639
12−0.6784−0.02340.60200.406628.3815171.70050.3980
13−0.5699−0.01970.56800.469532.7756148.68110.4320
14−0.4663−0.01610.53400.538737.6055129.58520.4660
15−0.3665−0.01260.50000.614942.9271113.52070.5000
16−0.2697−0.00930.46600.699248.809599.83950.5340
17−0.1751−0.00600.43200.792755.338988.05950.5680
18−0.0819−0.00280.39800.897162.624177.81530.6020
0.00000.00000.36791.000069.806569.80650.6321
190.01070.00040.36391.014370.805168.82430.6361
200.10330.00360.32991.146980.065360.86420.6701
210.19690.00680.29591.298690.651253.75670.7041
220.29250.01010.26191.4741102.904047.35590.7381
230.39130.01350.22791.6805117.314341.53900.7721
240.49500.01710.19391.9285134.621936.19850.8061
250.60620.02090.15992.2349156.015831.23470.8401
260.72880.02510.12592.6298183.582726.54450.8741
270.87030.03000.09183.1730221.496922.00080.9082
281.04740.03610.05784.0133280.159917.39400.9422
1.11091.47380.04804.3656304.750015.99000.9520
291.31850.04550.02385.7498401.379512.14090.9762
BOLD: The principal vibration stresses above and below the ƞ parameter, the R(t) index of 95%, and the R(t) index of the Sy 430 MPa.
Table 5. The proposed method and median rank method comparison results.
Table 5. The proposed method and median rank method comparison results.
FeatureProposed MethodMedian Rank Method
Weibull Shape Parameter
Equation (9)
β = 0.9102β = 0.7538
Weibull Scale Parameter
Equation (10)
ƞ = 69.8077ƞ = 69.8077
Principal Stresses,
Equations (17) and (16)
σ1 = 304.7600 MPa
σ2 = 15.9900 MPa
Material strength,
Sy = 430 MPa
R(t) = 0.8260R(t) = 0.7757
R(t) = 0.95
Equation (8)
σ1 = 1824.3190 MPa
σ2 = 2.6712 MPa
σ1 = 3591.0886 MPa
σ2 = 1.3570 MPa
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Barraza-Contreras, J.M.; Piña-Monarrez, M.R.; Torres-Villaseñor, R.C. Reliability by Using Weibull Distribution Based on Vibration Fatigue Damage. Appl. Sci. 2023, 13, 10291. https://doi.org/10.3390/app131810291

AMA Style

Barraza-Contreras JM, Piña-Monarrez MR, Torres-Villaseñor RC. Reliability by Using Weibull Distribution Based on Vibration Fatigue Damage. Applied Sciences. 2023; 13(18):10291. https://doi.org/10.3390/app131810291

Chicago/Turabian Style

Barraza-Contreras, Jesús M., Manuel R. Piña-Monarrez, and Roberto C. Torres-Villaseñor. 2023. "Reliability by Using Weibull Distribution Based on Vibration Fatigue Damage" Applied Sciences 13, no. 18: 10291. https://doi.org/10.3390/app131810291

APA Style

Barraza-Contreras, J. M., Piña-Monarrez, M. R., & Torres-Villaseñor, R. C. (2023). Reliability by Using Weibull Distribution Based on Vibration Fatigue Damage. Applied Sciences, 13(18), 10291. https://doi.org/10.3390/app131810291

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