Parameter Estimation and Applications of the Weibull Distribution for Strength Data of Glass Fiber

Wu, Yuxuan; Xie, Hanyang; Chiang, Jyun-You; Peng, Gang; Qin, Yan

doi:https://doi.org/10.1155/2021/9175170

Mathematical Problems in Engineering

On this page

Abstract Introduction Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 9175170 | https://doi.org/10.1155/2021/9175170

Parameter Estimation and Applications of the Weibull Distribution for Strength Data of Glass Fiber

Yuxuan Wu,¹Hanyang Xie,¹Jyun-You Chiang,¹Gang Peng,¹and Yan Qin²

Academic Editor: Francesco Zammori

Received18 May 2021

Revised05 Aug 2021

Accepted07 Sept 2021

Published22 Sept 2021

Abstract

Glass fiber is a good substitute for metal materials. However, in the process of manufacturing, it is necessary to carry out sampling inspection on its tensile strength to infer its quality. According to previous literatures, the strength data can be well fitted by the Weibull distribution, while the poor parameter estimation method cannot obtain reliable analysis results. Therefore, a new parameter estimation method is proposed. Based on the simulation results, it is found that the proposed parameter estimation method outperforms the other competitors to obtain reliable estimates of the Weibull parameters. Finally, the proposed parameter estimation method is applied to two real data sets of glass fiber strength for illustration. The results of data analysis show that our proposed parameter estimation method is more suitable for these data sets than other estimation methods.

1. Introduction

Glass fiber is a kind of high-performance composite fiber. Compared with carbon fiber, it has lower cost and wider application in the field of construction, electronics, chemical industry, and so on. In the manufacturing process of glass fiber, it is very important to determine the failure probability distribution of several nominally identical glass fiber samples by destructive test and then statistical analysis of their strength data. Figure 1 shows the typical fracture pattern for low strength fiber, which is from Feih et al. [1]. According to previous literatures, the strength data of fiber is usually described by two-parameter Weibull distribution. For example, Sakin and Ay [2] studied the bending fatigue properties of glass fiber reinforced polyester composite plates and fitted the data with two-parameter Weibull distribution. Acitas et al. [3] assumed that the strength data of glass fiber followed the Weibull distribution and estimated the parameters with maximum likelihood estimation (MLE). Zhu et al. [4] developed the technology of making high strength refractory ceramic fibers using fly ash, and the mechanical properties of a series of fly ash fibers were assessed by the Weibull distribution. Ferreira et al. [5] studied the Weibull analysis of tensile tested piassava fibers with different diameters.

Because of the good properties of the Weibull distribution, it is currently applied in various fields. Many researches have proposed the related distribution, neutrosophic statistics, or extended distribution and discussed its properties and applications, such as Alkarni et al. [6], Almongy et al. [7], Almetwally et al. [8], Ferretti et al. [9], Zhao et al. [10], Aslam et al. [11], Fernández et al. [12], Arif and Aslam [13], Aslam et al. [14], Aslam [15], Aslam et al. [16], Aslam and Arif [17], and Zhang et al. [18].

Before the reliability and failure rate of glass fiber are inferred, sample-based parameters’ estimation is required for the analysis process. There are many methods for estimating the unknown parameters and reliability analysis of the Weibull distribution. Naresh et al. [19] provided reliability analysis of tensile strength of glass/epoxy and carbon/epoxy composites based on the Weibull distribution. Malik et al. [20] evaluated reliability of a non-series-parallel system of seven components, and the Weibull distribution is used for the random variables associated with failure and operating times of the components. Jia [21] studied the reliability for the Weibull distribution with homogeneous heavily censored data. Johnson et al. [22] introduced many commonly used estimation methods for the Weibull distribution in the book, such as MLE, moment estimation (MM), percentile method (PM), method of logarithmic moment (MLM), and Bayesian estimation. Teimouri et al. [23] proposed an L-moment (LM) estimator for the Weibull distribution. Kantar [24] proposed generalized least squares and weighted least squares estimation methods for distributional parameters. Datsiou and Overend [25] investigated four estimation methods of the Weibull distribution for glass strength data. Hidekazu et al. [26] proposed an improved estimation method for MLE, which improved the estimation of Weibull shape parameter, scale parameter, and quantile. Acitas et al. [3] used particle swarm optimization to solve the MLE of the Weibull distribution and applied the results to the strength data of glass fiber. Although the above literatures on parameter estimation and reliability analysis of the Weibull distribution are valuable, there is a lack of research on the parameter estimation method of strength data fitting the Weibull distribution (Datsiou and Overend [25]). The most commonly used estimation method is MLE, while it does not perform well for the Weibull distribution, especially in small samples.

The aim of this article is to propose a new parameter estimation method for the Weibull distribution. The performance of proposed method outperforms other common estimation methods for small sample size. In addition, the proposed estimation method is applied to estimate the reliability and failure rate of products to evaluate quality or reliability.

The rest of this paper is organized as follows. In Section 2, some existing parameter estimation methods for the Weibull distribution are addressed. In Section 3, we propose a new parameter estimation method for the Weibull distribution. Based on the proposed parameter estimation method, it is applied to the estimation of reliability and failure rate. In Section 4, a Monte Carlo simulation is conducted to evaluate the performances of the proposed parameter estimation method. Then, two real strength data sets of glass fiber are used to demonstrate the proposed parameter estimation method and reliability analysis in Section 5, and some concluding remarks are provided in Section 6.

2. Review of Parameter Estimation Methods

Suppose the random variable T has a Weibull distribution with shape parameter and scale parameter , denoted as . The probability density function (PDF), cumulative distribution function (CDF), reliability function (RF), and failure rate function (FRF) are, respectively, given by

We assume is a random sample from the two-parameter Weibull distribution. Five common parameter estimation methods are reviewed as follows.

2.1. Method of Moments (MM)

The MM is a traditional parameter estimation method, which is widely used for parameter estimation of various distributions. This estimation method assumes that the first k population moments are equivalent to the first k sample moments and estimates the unknown parameters by solving the equations derived from the equivalence relation. Therefore, it can be used for two-parameter Weibull distribution, whose parameter estimation satisfies the following two equations:where and are the estimates of and , respectively, is gamma function, and and are the first and second noncentral sample moments, respectively.

2.2. Method of L-Moment (LM)

According to Teimouri et al. [23], they first proposed a L-moment estimator for two-parameter Weibull distribution. The LM estimator of parameters can be obtained by the following equations:where denotes the order statistic in the sample size n, and is the average of sample , i = 1, ..., n.

2.3. Maximum Likelihood Estimator (MLE)

The MLE method is one of the most commonly used parameter estimation methods. The MLE of two-parameter Weibull distribution can be obtained by the following equations:

Since the close-form solution of the MLE is not available, some numerical optimization methods have been proposed to overcome this issue.

2.4. Method of Percentile (PM)

According to Johnson et al. [22], the PM estimates of the two-parameter Weibull distribution arewhere and are the estimators of and percentile, respectively. Former studies show that PM has a good estimation effect when .

3. Proposed Parameter Estimation Method

According to the distribution property, if Y = , then Y becomes a uniform distribution from 0 to 1. Moreover, the mean and variance of Y are 1/2 and 1/12, respectively, which can be expressed as

Suppose is a random sample with sample size n from . Theoretically, , also follows a uniform distribution from 0 to 1. However, the sample has inevitable random error. Therefore, the idea of the proposed parameter estimation method is to find the estimators of and , which can make approximate the uniform distribution as much as possible. In other words, the estimators of and can satisfy as much as possible that the mean and variance of are 1/2 and 1/12, respectively. Because the sample mean and sample variance are unbiased estimators of population mean and population variance, respectively, the proposed parameter estimators can be obtained by solving the following equations:

Equation (15) can be rewritten as follows:

Therefore, the estimators of and can be obtained from equation (16) and denoted as and , respectively. Since the proposed parameter estimation method is based on CDF, we name this method CDFM (cumulative distribution function method). Then, and in equations (3) and (4) are replaced by and , respectively, so the RF and FRF based on CDFM can be obtained and expressed as

Since the close-form solution of the proposed estimation is not available, the estimators can be obtained by using numerical optimization methods such as the quasi-Newton algorithm [27]. Figure 2 shows the flow chart of the parameter estimation process. When we obtain a sample data T, the estimators of the parameters, and , can be obtained according to the proposed parameter estimation method. Then, the goodness of fit test is used to test whether the data comes from the Weibull distribution with parameters and . Common goodness of fit tests include Kolmogorov-Smirnov test (K-S test) and Anderson-Darling test. When the test results show that the data does not follow the Weibull distribution, we must find the exact distribution of the data before inference. When the test results show that the data may come from the Weibull distribution, we can use equations (17) and (18) to estimate RF and FRF.

4. Monte Carlo Simulation Study

In this section, all simulation results are obtained by R software. The R codes can be obtained from the authors upon request. In this simulation study, we compare proposed estimation method CDFM with MLE, LM, PM, and MM. For convenience, we denoted the results estimated by these methods as , , , , and , respectively. Considering the sample size n = 10, 20, 30, 40, 50, 60, 70, and data set are generated from the , where , 6 and . Ten thousand iteration runs are used to evaluate the bias and MSEs of five estimation methods. The bias and MSE can be evaluated based on equations (19) and (20), respectively:where can be or and the subscript i denotes the five estimation methods. All obtained biases and MSEs are reported in Tables 1 and 2. In view of Tables 1 and 2, the following results are found:(1)With the increase of sample size, the performance of bias and MSE of all methods are getting better and better.(2)For the bias of shape parameter , the performance of outperforms other methods in most cases, followed by .(3)For the MSE of , the performance of outperforms other methods in most cases, followed by and .(4)For the bias of scale parameter , the performance of outperforms other methods for n = 10, except . And the performance of outperforms other methods in most cases.(5)For the MSE of , the performance of outperforms other methods, while the performance of is the worst.(6)With the increase of the population shape parameter, the bias of is larger and larger, and the bias of is smaller and smaller for all methods.(7)When = 1, the performance of each estimation method is similar to that of = 2.

According to the above findings, the bias of outperforms other methods. Although the bias of is worse than , the difference is not significant. Therefore, using the proposed estimation method for inference can get reliable results.

In addition, in order to compare the performance of the RF and FRF based on different parameter estimation methods, we consider that the sample size n = 10, 20, 30, 40, 50, 60, 70, , 6, , and strength t = 2. All the simulation results are shown in Figures 3 and 4. In view of Figures 3 and 4, the following results are found:(1)The performance of biases of the RF and FRF estimated by CDFM method outperforms other methods in most cases.(2)With the increase of or sample size n, the performance of biases of each estimation method is better for RF and FRF.(3)Among all estimation methods, PM method performs the worst for biases of RF and FRF.(4)Among all estimation methods, PM method also performs the worst for MSEs of RF and FRF.(5)With the increase of or sample size n, the performance of MSEs of each estimation method is better for RF and FRF.(6)Except PM method, the performance of MSEs of other estimation methods is very close for RF and FRF.

5. Real Data Analysis

5.1. Example 1

The tensile strength data of glass fiber, reported by Smith and Naylor [28] and further studied by Acitas et al. [3], Morais and Barreto-Souza [29], Alizadeh et al. [30], and Jones and Faddy [31], is used to illustrate the methodologies developed in this paper. In the article of Smith and Naylor [28], two experimental strength data sets of glass fiber of two lengths, 1.5 cm and 15 cm, are provided, originally obtained by workers at the UK National Physical Laboratory. Unfortunately, the units of measurement are not given in the paper. In this example, the data of 1.5 cm length is selected for analysis. The strength data is shown in Table 3.

Two-parameter Weibull distribution is used to characterize the data in Table 3. Moreover, the CDFM, MLE, LM, PM, and MM estimation methods are used to obtain the estimates of the Weibull distribution parameters. All the estimation results are reported in Table 4. The estimators of LM are consistent with those of Acitas et al. [3]. The p value of the K–S test in Table 4 indicates that the Weibull distribution can well fit this data set and we found the CDFM estimation method is the best estimation method to obtain the estimates of the Weibull distribution parameters. In addition, Figure 5 presents the histogram of sample and the fitted PDFs of the Weibull distributions by using the parameter estimation results in Table 4. Figure 6 presents empirical distribution function of sample and the fitted CDFs of the Weibull distributions by using the parameter estimation results in Table 4. As shown in Figure 5, the PDF corresponding to CDFM is closest to the sample histogram, which also shows that the estimation result of CDFM is more suitable for this data. From Figure 6, we can also see that the curve corresponding to CDFM is closest to the empirical distribution function of the sample. The results in Figure 6 also show that the CDFM is more suitable for this data.

Suppose we want to estimate the reliability and failure rate of glass fiber, we can substitute equations (17) and (18) according to the parameter estimator results in Table 4, and the results are shown in Figure 7. From Figure 7, the results show that when t < 1.6, the RF based on CDFM is higher than other methods except PM method; when t > 1.6, the RF based on CDFM is lower than other methods except PM method. In addition, the FRF based on CDFM is higher than other methods except PM method.

5.2. Example 2

In this example, the data of 15 cm length reported by Smith and Naylor [28] is selected for analysis. The strength data is shown in Table 5. Two-parameter Weibull distribution is also used to characterize the data in Table 5. Moreover, the CDFM, MLE, LM, PM, and MM estimation methods are used to obtain the estimates of the Weibull distribution parameters. All the estimation results are reported in Table 6. The p value of the K-S test in Table 6 indicates that the Weibull distribution can well fit this data set and we found the CDFM estimation method is the best estimation method to obtain the estimates of the Weibull distribution parameters. In addition, Figure 8 presents the histogram of sample and the fitted PDFs of the Weibull distributions by using the parameter estimation results in Table 6. Figure 9 presents empirical distribution function of sample and the fitted CDFs of the Weibull distributions by using the parameter estimation results in Table 6. The results in Figures 8 and 9 show that the fitting performance of these estimation methods for this sample is similar.

Suppose we want to estimate the reliability and failure rate of glass fiber, we can substitute equations (17) and (18) according to the parameter estimator results in Table 6, and the results are shown in Figure 10. It can be seen from Figure 10 that the RF and FRF results estimated by each method are almost the same, because the parameter estimator results of this data set for each estimation method are almost the same.

6. Conclusions

This study aims to solve the problem of parameter estimation for the Weibull distribution and then infer the reliability of products.

6.1. Unique Contributions of the Study

We propose a new parameter estimation method (CDFM) and use this method to estimate the reliability of products. According to Monte Carlo simulation results, it can be found that the CDFM estimation method outperforms the other competitors to obtain reliable estimates of two-parameter Weibull parameters. Therefore, two strength data sets of glass fiber provided by Smith and Naylor [28] are estimated according to the proposed estimation method. The analysis results show that the estimates of CDFM are more suitable for these data sets. And according to the estimates of CDFM, the reliability and failure rate of glass fiber are also estimated. From the analysis results, it can be seen that the reliability and failure rate estimated by different parameter estimates will be different, so it is very important to use a reliable parameter estimate for inference. In addition, the estimation method proposed in this paper can be applied to other fibers, materials, or tests as long as the data follow the Weibull distribution.

6.2. Limitations

In particular, the proposed parameter estimation method can be used as long as the data follows the Weibull distribution, but it cannot guarantee the same performance when the data follows the other distributions.

6.3. Future Research

In the current work, we only consider the parameter estimation and application of two-parameter Weibull distribution without inferring the theoretical properties and limit distribution of the estimators. In future research, in addition to applying this estimation method to other distributions, we can also infer the limit distribution and statistical properties of these estimators. We are currently exploring possible extensions in this direction.

Data Availability

Two strength data sets of glass fiber were provided by Smith and Naylor (1987).

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the Fundamental Research Funds for the Central Universities of China, Program no. JBK2001024.

References

S. Feih, A. Thraner, and H. Lilholt, “Tensile strength and fracture surface characterisation of sized and unsized glass fibers,” Journal of Materials Science, vol. 40, no. 7, pp. 1615–1623, 2005.
View at: Publisher Site | Google Scholar
R. Sakin and İ. Ay, “Statistical analysis of bending fatigue life data using weibull distribution in glass-fiber reinforced polyester composites,” Materials & Design, vol. 29, no. 6, pp. 1170–1181, 2008.
View at: Publisher Site | Google Scholar
S. Acitas, C. H. Aladag, and B. Senoglu, “A new approach for estimating the parameters of Weibull distribution via particle swarm optimization: an application to the strengths of glass fibre data,” Reliability Engineering & System Safety, vol. 183, 2019.
View at: Publisher Site | Google Scholar
Z. Zhu, C. Zhang, S. Meng, and Z. Shi, “A statistical damage constitutive model based on the weibull distribution for alkali-resistant glass fiber reinforced concrete,” Materials, vol. 12, no. 12, 2019.
View at: Publisher Site | Google Scholar
D. C. D. O. N. Ferreira, A. D. S. Ferreira, and S. N. Monteiro, “Weibull analysis of tensile tested piassava fibers with different diameters,” Matèria. Revista Internacional d’Art, vol. 23, no. 4, 11 pages, 2018.
View at: Publisher Site | Google Scholar
S. Alkarni, A. Z. Afify, and I. Elbatal, “The extended inverse weibull distribution: properties and applications,” Complexity, vol. 2020, Article ID 3297693, 12 pages, 2020.
View at: Publisher Site | Google Scholar
H. M. Almongy, E. M. Almetwally, and R. Alharbi, “The weibull generalized exponential distribution with censored sample: estimation and application on real data,” Complexity, vol. 2021, Article ID 6653534, 15 pages, 2021.
View at: Publisher Site | Google Scholar
E. M. Almetwally, M. A. H. Sabry, and R. Alharbi, “Marshall–Olkin alpha power weibull distribution: different methods of estimation based on type-I and type-II censoring,” Complexity, vol. 2021, Article ID 6653534, 2021.
View at: Publisher Site | Google Scholar
C. Ferretti, P. Ganugi, and G. Pisano, “Abraded glass strength: an ad hoc fitting protocol based on the change of variable theorem,” Mathematical Problems in Engineering, vol. 2020, Article ID 8519426, 12 pages, 2020.
View at: Publisher Site | Google Scholar
J. Zhao, Z. Ahmad, and E. Mahmoudi, “A new class of heavy-tailed distributions: modeling and simulating actuarial measures,” Complexity, vol. 2021, Article ID 6653534, 18 pages, 2021.
View at: Publisher Site | Google Scholar
M. Aslam, S. Balamurali, C.-H. Jun, and M. Ahmad, “A two-plan sampling system for life testing under weibull distribution,” Industrial Engineering and Management Systems, vol. 9, no. 1, pp. 54–59, 2010.
View at: Publisher Site | Google Scholar
A. J. Fernández, C. J. Pérez-González, M. Aslam, and C.-H. Jun, “Design of progressively censored group sampling plans for weibull distributions: an optimization problem,” European Journal of Operational Research, vol. 211, no. 3, pp. 525–532, 2011.
View at: Publisher Site | Google Scholar
O. H. Arif and M. Aslam, “A new sudden death chart for the weibull distribution under complexity,” Complex & Intelligent Systems, vol. 32, pp. 1–9, 2021.
View at: Google Scholar
M. Aslam, G. Srinivasa Rao, and N. Khan, “Single-stage and two-stage total failure-based group-sampling plans for the Weibull distribution under neutrosophic statistics,” Complex & Intelligent Systems, vol. 7, no. 2, pp. 891–900, 2021.
View at: Publisher Site | Google Scholar
M. Aslam, “A new failure-censored reliability test using neutrosophic statistical interval method,” International Journal of Fuzzy Systems, vol. 21, no. 4, pp. 1214–1220, 2019.
View at: Publisher Site | Google Scholar
M. Aslam, P. Jeyadurga, and S. Balamurali, “Time-truncated group plan under a weibull distribution based on neutrosophic statistics,” Mathematics, vol. 7, no. 10, 2019.
View at: Google Scholar
M. Aslam and O. Arif, “Testing of grouped product for the weibull distribution using neutrosophic statistics,” Symmetry, vol. 10, no. 9, p. 403, 2018.
View at: Publisher Site | Google Scholar
C. Zhang, D. Y. Li, and X. P. Kang, “Neutrosophic fusion of rough set theory: an overview,” Computers in Industry, vol. 115, Article ID 103117, 2020.
View at: Google Scholar
K. Naresh, K. Shankar, and R. Velmurugan, “Reliability analysis of tensile strengths using weibull distribution in glass/epoxy and carbon/epoxy composites,” Composites Part B: Engineering, vol. 133, no. 15, pp. 129–144, 2018.
View at: Publisher Site | Google Scholar
S. C. Malik, S. K. Chauhan, and N. Ahlawat, ““Reliability analysis of a non series–parallel system of seven components with weibull failure laws,” International Journal of Systems Assurance Engineering and Management, vol. 11, no. 1, pp. 577–582, 2020.
View at: Publisher Site | Google Scholar
X. Jia, “Reliability analysis for weibull distribution with homogeneous heavily censored data based on bayesian and least-squares methods,” Applied Mathematical Modelling, vol. 83, pp. 169–188, 2020.
View at: Publisher Site | Google Scholar
N. L. Johnson, S. Kotz, and N. Balakrishnan, Continuous Univariate Distributions, Wiley, NewYork, NY, USA, 2nd edition, 1994.
M. Teimouri, S. M. Hoseini, and S. Nadarajah, “Comparison of estimation methods for the Weibull distribution,” Statistics, vol. 47, no. 1, pp. 93–109, 2013.
View at: Publisher Site | Google Scholar
Y. M. Kantar, “Generalized least squares and weighted least squares estimation methods for distributional parameters,” Revstat-Statistical Journal, vol. 13, no. 3, pp. 263–282, 2015.
View at: Google Scholar
K. C. Datsiou and M. Overend, “Weibull parameter estimation and goodness-of-fit for glass strength data,” Structural Safety, vol. 73, pp. 29–41, 2018.
View at: Publisher Site | Google Scholar
T. Hidekazu, P. Nabendu, and K. L. Wooi, “On improved estimation under Weibull model,” Journal of Statistical Theory and Practice, vol. 12, no. 1, pp. 48–65, 2018.
View at: Google Scholar
R. H. Byrd, P. Lu, J. Nocedal, and C. Zhu, “A limited memory algorithm for bound constrained optimization,” SIAM Journal on Scientific Computing, vol. 16, no. 5, pp. 1190–1208, 1995.
View at: Publisher Site | Google Scholar
R. L. Smith and J. C. Naylor, “A comparison of maximum likelihood and bayesian estimators for the three- parameter weibull distribution,” Applied Statistics, vol. 36, no. 3, pp. 358–369, 1987.
View at: Publisher Site | Google Scholar
A. L. Morais and W. Barreto-Souza, “A compound class of Weibull and power series distributions,” Computational Statistics & Data Analysis, vol. 55, no. 3, pp. 1410–1425, 2011.
View at: Publisher Site | Google Scholar
M. Alizadeh, S. Rezaei, and S. F. Bagheri, “On the estimation for the weibull distribution,” Annals of Data Science, vol. 2, no. 4, pp. 373–390, 2015.
View at: Publisher Site | Google Scholar
M. C. Jones and M. J. Faddy, “A skew extension of the t -distribution, with applications,” Journal of the Royal Statistical Society: Series B, vol. 65, no. 1, pp. 159–174, 2003.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Yuxuan Wu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies