Research Article  Open Access
Chunping Li, Huibing Hao, Fang Xu, Guotao Zhao, "Reliability Evaluation of Multiple DCFP System subject to Shifting Threshold", Shock and Vibration, vol. 2020, Article ID 9206239, 11 pages, 2020. https://doi.org/10.1155/2020/9206239
Reliability Evaluation of Multiple DCFP System subject to Shifting Threshold
Abstract
This paper focuses on system reliability analysis with dependent competing failure process due to soft failure and hard failure. Some new probabilistic methods based on cumulative shock model and nonlinear Wiener process under different shifting thresholds situation are obtained. Considering that nonlinearity exists extensively in practice, the continuous soft failure process is governed by random effected nonlinear Wiener process. Firstly, reliability evaluation models for hard failure and soft failure are obtained under the cumulative shock, respectively. Furthermore, some system reliability models under different shifting thresholds situation are studied, in which failure threshold will decrease after a certain number of shocks. A real numerical example about fatigue crack growth dataset is carried out to demonstrate the proposed procedure. Numerical results indicate that both random shocks and shifting threshold have significant effect on system reliability. Finally, some sensitivity analysis are also been given.
1. Introduction
In actual operation environment, most of products are often exposed to different types of risk. Those risks may come from products themselves or external environments, and the same risk can also result in different failure categories. Specifically, we consider two failure processes: soft failure process and hard failure process, where the soft failure is caused by a gradual degradation and additional abrupt degradation damages due to a shock process, and hard failure is caused by an abrupt stress from the same shock process. As we know that these two failure processes exist simultaneously for most products (or systems). Therefore, the failure of product is usually the result of the competition between them, and we call it competing failure.
In recent years, theoretical studies for competing failure have been extensively explored in reliability analysis. Li and Pham [1] conducted system reliability analysis subject to degradation and shock process. Lehman [2] studied the competing failure model due to degradation and random shock. Klutke and Yang [3] studied the availability of maintained systems affected by the degradation process and random shock process. Zhu et al. [4] established an independent competing failure model based on sudden failure and degradation process. In [5], a continuous degradation system was proposed by considering soft failure and hard failure. In other studies, several random shock models are constructed, such as cumulative shock model [6, 7], extreme shock model [8], mixed model [9], and δshock model [10, 11]. Most previous studies assumed that the two failure processes were independent. However, in actual operation circumstances, these two failure processes may be dependent because the effects from the same shock process contribute to both failure processes. When they are dependent, it creates a new problem to analyze the system reliability.
Recently, the reliability analysis of system with multiple dependent competing failure processes (DCFP) has been analyzed by many researchers. Peng et al. [12] introduced reliability modeling for systems with multiple DCFP. Wang et al. [13] proposed the dependent reliability assessment model under the degradation and shock process. Qiu and Cui [14] proposed the dependent twostage failure process model with competing failures. Song et al. [15] considered multicomponent DCFP. Raflee et al. [16] proposed a reliability model for devices subject to DCFP with the changing degradation rate. An and Sun [17] proposed the reliability model for system experiencing DCFP with shock loads above a certain level. Hao and Yang [18] proposed a reliability analysis for DCFP with changing degradation rate and hard failure threshold levels. Fan et al. [19] proposed Bayesian approach for remaining useful life prediction of DCFP. Kong et al. [20] developed reliability models with sensorbased calibrations with multiple DCFP. Gao et al. [21] considered reliability models subject to DCFP with degradationshock dependences. In all of these DCFP models, the failure thresholds are considered to be fixed constants, implying that the component resistance to failure does not deteriorate as the system ages. However, when the system sustains a series of shocks, the product may become more vulnerable and its resistance to failure is weakening. Therefore, the failure threshold decreases after multiple shocks, but only few papers considered the changing failure threshold [22, 23].
Considering that the stochastic process model can provide flexibility to describe failure generating mechanisms and operating environment characteristics, many articles used stochastic process approach to model degradation path [24–26], such as the Wiener process and Gamma process. In addition, most of the above papers assumed that the continuous degradation path was represented by a linear function of time t, but the nonlinearity exists extensively in practice [27]. It is worth pointing out that little research has been carried out to consider nonlinearity stochastic process and changing failure threshold for system experiencing DCFP and those have become a critical issue for many complex systems. In this paper, some new reliability evaluation models are conducted on the basis of DCFP. Firstly, some reliability models with hard failure and soft failure are established, where the soft failure process is governed by the random effect nonlinear Wiener process. On that basis, the DCFP models under the different shifting threshold are considered. Numerical example indicates that the shocks and the shifting threshold have significant effects on system reliability.
The remainder of this paper is organized as follows. Failure process including degradation and random shock is described in Section 2. Some reliability models about hard failure and soft failure are established in Section 3. In Section 4, some new reliability models are obtained for system experiencing DCFP with different shifting thresholds. A numerical example is presented in Section 5, and the sensitivity analysis is also given. In Section 6, some main conclusions are obtained.
2. Dependent Competing Failure Process
Generally speaking, product may fail due to two different failure modes: soft failure and hard failure. These two failure processes are dependent because both of them are subject to the same shock. The product will fail when either of them exceeds its corresponding threshold value. In this section, the cumulative shock model is applied.
As shown in Figure 1(a), let be the degradation volume of the continuous degradation until time t. Shock loads will cause additional abrupt damage, Y_{i} (i = 1, 2, 3, …,), which contribute to acceleration of the degradation process. When the overall degradation is beyond critical strength level H_{1}, the soft failure occurs.
(a)
(b)
Considering that so many external shocks can affect the degradation performance of the product, it is a reasonable assumption that the product may become more vulnerable after multiple shocks. Therefore, in Figure 1(b), we suppose that the critical value of degradation failure decreases from H_{1} to H_{2} after the arrival of m shocks and m is the number of shocks’ occurrences.
In addition, the same shock loads may also result in hard failure, as shown in Figure 2(a). Let denote the magnitude of ith (i = 1, 2, 3,…,) shock, and hard failure occurs when cumulative shock loads magnitude exceeds threshold value D_{1}. Similarly, the hard threshold value D_{1} can be affected by the random shock process. In Figure 2(b), the critical strength value decreases from D_{1} to D_{2} after the arrival of m shocks, and m is the number of shocks’ occurrences.
(a)
(b)
3. Reliability Modeling of Hard Failure and Soft Failure
3.1. Reliability Modeling of Hard Failure due to Random Shock
During normal use stage, suppose that the product is affected by random shock, and random shocks reach according to the Poisson process with rate λ. Let denote the total number of shocks in given time interval (0, t]. Then, the probability function of N (t) can be obtained as
Let denote the magnitude of the ith shock, and assume that is an independent and identically distributed random variable following a normal distribution as W_{i}∼N () and = 0, and then we can obtain
Let S_{i} denote the total damage of the ith shock magnitude W_{i} (i = 1, 2, 3, …,); then, we can obtain
Via the property that the sum from independent normal distribution is still normal distribution, by using equations (2) and (3), the probability function of the total damage S_{i} can be obtained as
Let D (t) denote the total damage in given time interval (0, t] and N (t) denote the total number of shocks until time t. Then, we can obtain
3.1.1. Case I: Reliability Analysis under a Fixed Hard Threshold Value
As shown in Figure 2(a), let T be the hard failure time and D_{1} be threshold value of the hard failure. By using the total probability formula and the compound Poisson process (5), the reliability of hard failure can be obtained as
3.1.2. Case II: Reliability Analysis under a Shifting Hard Threshold Value
As shown in Figure 2(b), after m shocks, the hard failure critical value decreases from D_{1} to D_{2}. Similarly, in such a shifting hard threshold situation, the reliability function of hard failure can be obtained aswhere m is the number of shocks’ occurrences when failure threshold decreases.
3.2. Reliability Modeling of Soft Failure due to Degradation and Shock
Considering that the nonlinearity exists extensively in practice and the linear model cannot trace the dynamics, a nonlinearity stochastic process is used to describe the degradation path. A welladopted form for the nonlinearity Wiener process can be expressed aswhere a is the drift degradation rate, σ is the diffusion coefficient, is the transformed time scale, and is the standard Brownian motion representing a time correlated structure.
Considering that product possibly experiences different sources of variations during its operation, for a degradation model to be realistic, it is more appropriate to incorporate unit to unit variability in the degradation process. To capture the unit to unit differences, the random effect model is used in this paper. The nonlinearity degradation model with random effect can be expressed as
For simplicity, we assume that random variables a and are independent. From the properties of the Wiener process, we can obtain that the nonlinearity degradation performance X(t) is normally distributed as
3.2.1. Case III: Reliability Analysis under Fixed Degradation Threshold Value
Suppose that the random shock can cause an abrupt damage. The abrupt damages are measured by the shock damage sizes as {Y_{1}, Y_{2}, … ,}. Let S (t) be the cumulative damage size due to random shocks until time t; then, we can obtain
As shown in Figure 1(a), the total damage of system is the sum of the cumulative effect of continuous degradation damage and abrupt degradation damages due to shocks until time t. The soft failure occurs when the overall degradation is beyond the threshold value H_{1}. From Figure 1(a), we can obtain
In this article, we assume that the magnitude of gradual degradation X (t) follows normal distribution as N (μ (t), σ^{2} (t)). And the shock damage sizes Y_{i} are independently distributed random variables taking the form as Y_{i}∼N (). Considering the independence of the two random variables, follows the Gaussian distribution as
In most engineering problems, the mean and standard deviation of a degradation process increase in terms of time. The Gaussian process with increasing mean and standard deviation could be used to describe these engineering problems. Furthermore, when is a Gaussian process with increasing mean and standard deviation, is a monotone decreasing function with respect to time.
As shown in Figure 1(a), let T be the soft failure time and H_{1} be the corresponding threshold value. Then, considering the monotone increasing features, the reliability function of the soft failure process can be obtained as
3.2.2. Case IV: Reliability Analysis under a Shifting Degradation Threshold Value
As shown in Figure 1(b), after m shocks, the soft failure critical value decreases from H_{1} to H_{2}. In such a shifting threshold situation, the reliability function can be obtained as
4. Reliability Modeling and Analysis for DCFP
From the above analysis, the reliability functions under hard failure and soft failure are established, respectively. In the following research, considering that the system subjects to both degradation and random shocks, some new dependent competing reliability modeling and assessment models are established.
4.1. Reliability Analysis due to a Fixed Threshold Value
As shown in Figures 1(a) and 2(a), we can find that the system includes two types of dependent competing failure processes. Thus, the DCFP system reliability at time t is the probability that the total degradation is less than the soft threshold level () and the cumulative magnitude is less than the hard threshold value (). That is to say, the system reliability function can be obtained as
4.2. Reliability Analysis due to a Shifting Hard Threshold Value
As shown in Figure 2(b), after the shocks happened m times, the hard failure threshold decreases from D_{1} to D_{2}. In such a shifting hard threshold situation, we can get the DCFP system reliability at time t as
4.3. Reliability Analysis due to a Shifting Soft Threshold Value
Similar to Section 4.2, in Figure 1(b), after the shocks happened m times, the soft failure threshold decreases from H_{1} to H_{2}. In such a shifting soft failure threshold situation, we can get DCFP system reliability at time t as
5. Numerical Example: A Fatigue Crack Growth Data Application
In this section, a case study is given to illustrate the proposed models and methods. This example is based on the fatigue crack growth data of an alloy in paper from Lu and Meeker [28]. The datasets were collected from 13 test samples, and all samples have the same initial crack length of 0.90 inches. The loading cycles of the datasets are considered to be the loading time. Figure 3 shows part of the data, and time unit is in 10^{4} cycles.
Similarly to [13], we suppose sizes of random shock loads, (i = 1, 2, 3, … ,), are assumed to follow normal distribution, W_{i}∼N (2, 0.5), which are measured in units of component life. The shock damage sizes Y_{i}∼N (0.02, 0.01), the degradation failure threshold value H = 2.0 inches, catastrophic failure threshold value D_{1} = 35 units, and the parameter λ = 0.5. In addition, we assume that the number m = 3, the lower level of degradation failure threshold value H_{2} = 1.8 mm, and the lower level of catastrophic failure threshold value D_{2} = 25 units.
5.1. Estimation of Unknown Parameters
Firstly, we use the random effect nonlinear Wiener process model to fit the fatigue crack growth data. According to the suggestion of Pan et al. [29], let be a reasonable choice for the degradation modeling about the fatigue crack growth data.
Here, we discuss the estimation of parameters required to implement the reliability function in equation (9). The unknown parameters are
By using MCMC method and WinBUGS [30], the parameters posterior mean, standard error, Monte Carlo error, and 95% HPD interval are obtained in Table 1.

5.2. Reliability Analysis Subject to Only Hard Failure or Soft Failure
According to equations (6), (7), (14), and (15), under the fixed and shifting failure threshold, the reliability curves are plotted in Figures 4 and 5, respectively. From Figures 4 and 5, we can obtain that the shifting threshold have significant effect on the failure process, both hard failure and soft failure.
5.3. System Reliability Analysis Subject to DCFP
According to equations (16)–(18), the reliability curve of different cases in Section 4 are plotted in Figure 6, respectively. From Figure 6, we can obtain the following conclusions:(a)Compared with the degradation process without shocks, the reliability of system with shocks is significantly lower. Conclusion indicates that the shocks have significant effect on the reliability of system.(b)When system sustains a series of shocks, it is easier to fail. Conclusion indicates that the shock can decrease the failure threshold, both the soft failure and hard failure threshold. Thus, when the shifting thresholds are considered, the reliability of system is also significantly lower than the fixed thresholds.(c)For different cases, we can find that the reliability almost keeps 1 when t < 10^{5} cycles because the effect of random shocks is not significant during that time, and the gradual degradation amount is not large enough to cause any failure. In the next time period, with the gradual degradation and random shocks, reliability is decreasing quickly; the random shocks become a significant factor for the system reliability. Eventually, when t > 2.5 × 10^{5} cycles, reliability drops to a quite low level.
Reliability curve based on degradation without shocks is also provided in Figure 6. According to the results, it can be concluded that random shocks will accelerate the failure process which validates the effectiveness of models constructed in Sections 3 and 4. Moreover, when hard failure threshold value decreases, failure will occur sooner.
5.4. Sensitivity Analysis of Different Parameters
This section only discusses the sensitivity analysis of the fixed threshold case, and different shifting threshold cases can be analyzed similarly.
Sensitivity analysis under the threshold H_{1} of soft failure, the threshold D_{1} of hard failure, and the random shocks rate λ are investigated with regard to reliability distribution, and the corresponding curve are plotted in Figures 7–9, respectively.
From Figures 7 and 8, we can find that both soft failure threshold and hard failure threshold have significant effect on the reliability performance. Reliability curve shifts to right as the threshold value increases, which implies a better reliability performance with a larger threshold value. However, the changing thresholds by manmodified can improve reliability, but that does not mean system becomes more reliable. In engineering practice, many factors should be considered comprehensively in order to scientifically determine the reliability threshold, such as the system characteristics factor, operating conditions factor, and working environment factor.
In addition, by comparing Figures 7 and 8, it can be found that the influence of soft failure threshold on system reliability exists in the whole life cycle; however, the hard failure threshold has different influences on system reliability in different time periods, and the closer it is to the end of the life cycle, the greater the influence will be. This is because the number of initial shocks is small, and the cumulative effect of shocks is not obvious. After a period of time, the cumulative impact effect has been very significant, and if hard failure threshold drops, the impact on system reliability will become more obvious.
Figure 9 shows that the shock arrival rate, λ, has a significant effect with system reliability. Obviously, comparing with the effect of the failure threshold, the reliability is more sensitive for the change of impact arrival rate λ. It is reasonable to get such results because intensively frequent random shocks, larger sizes of shock loads, and shock damages will resultantly make a component more vulnerable. Therefore, when the number of shocks is large enough, the reliability of the system will drop to a very low level; the higher the impact arrival rate, the earlier the system enters the unreliable state.
Therefore, by setting higher failure threshold H_{1}, H_{2}, D_{1}, and D_{2}, the system seems to be more reliable. However, it is very difficult to make such assumption in reality, since reliability threshold depends on many conditions, including manufacturing processes, importance, work environment, and maintaining measures.
6. Conclusions
In this paper, some new reliability evaluation models are established by considering dependent competing failure process due to soft failure and hard failure under different shifting thresholds. A real numerical example about fatigue crack growth dataset is carried out to demonstrate the proposed procedure. Compared with the degradation process without shocks, models established in this paper indicate that random shocks have significant effect on the failure process. When considering the shiftingthreshold situation, reliability decreases even faster. From the sensitivity analysis, the shifting threshold and the shock arrival rate have a significant effect with the system reliability.
The proposed models in this paper only consider one degradation path. In real situation, components may have multiple degradation measures, and this will be the focus of future research. In addition, this paper only considers the cumulative shock model and wiener process model; in the future, some other works can be carried out, such as extreme shock model, δshock model, Gamma process model, and inverse Gaussian process model.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this article.
Authors’ Contributions
The authors contributed equally to the presented mathematical framework and the writing of the paper. All authors have read and approved the final manuscript.
Acknowledgments
This work was supported by the Humanity and Social Science Foundation of Ministry of Education of China (no. 19YJAZH039), Key Project of Hubei Provincial Education Department (no.D20172701), Technology Creative Project of Excellent Middle & Young Team of Hubei Province (no. T201920), National Bureau of Statistics of China (no.2017LY73), and Humanity and Social Science Foundation of Ministry of Education of China (no. 20YJAZH035).
References
 W. Li and H. Pham, “Reliability modeling of multistate degraded systems with multicompeting failures and random shocks,” IEEE Transactions on Reliability, vol. 54, no. 2, pp. 297–303, 2005. View at: Publisher Site  Google Scholar
 A. Lehmann, “Joint modeling of degradation and failure time data,” Journal of Statistical Planning and Inference, vol. 139, no. 5, pp. 1693–1706, 2009. View at: Publisher Site  Google Scholar
 G.A. Klutke and Y. Yang, “The availability of inspected systems subject to shocks and graceful degradation,” IEEE Transactions on Reliability, vol. 51, no. 3, pp. 371–374, 2002. View at: Publisher Site  Google Scholar
 Y. Zhu, E. A. Elsayed, H. Liao, and L. Y. Chan, “Availability optimization of systems subject to competing risk,” European Journal of Operational Research, vol. 202, no. 3, pp. 781–788, 2010. View at: Publisher Site  Google Scholar
 E. Deloux, B. Castanier, and C. Bérenguer, “Predictive maintenance policy for a gradually deteriorating system subject to stress,” Reliability Engineering & System Safety, vol. 94, no. 2, pp. 418–431, 2009. View at: Publisher Site  Google Scholar
 M. Kijima and T. Nakagawa, “A cumulative damage shock model with imperfect preventive maintenance,” Naval Research Logistics, vol. 38, no. 2, pp. 145–156, 1991. View at: Publisher Site  Google Scholar
 Y. Wang and H. Pham, “Imperfect preventive maintenance policies for twoprocess cumulative damage model of degradation and random shocks,” International Journal of System Assurance Engineering and Management, vol. 2, no. 1, pp. 66–77, 2011. View at: Publisher Site  Google Scholar
 J. Chen and Z. Li, “An extended extreme shock maintenance model for a deteriorating system,” Reliability Engineering & System Safety, vol. 93, no. 8, pp. 1123–1129, 2008. View at: Publisher Site  Google Scholar
 G. Allan, “Mixed shock models,” Bernoulli, vol. 7, no. 3, pp. 541–555, 2001. View at: Publisher Site  Google Scholar
 J.M. Bai, Z.H. Li, and X.B. Kong, “Generalized shock models based on a cluster point process,” IEEE Transactions on Reliability, vol. 55, no. 3, pp. 542–550, 2006. View at: Publisher Site  Google Scholar
 Y. Lam, “A geometric process deltashock maintenance model,” IEEE Transactions on Reliability, vol. 58, no. 2, pp. 389–396, 2009. View at: Publisher Site  Google Scholar
 H. Peng, Q. Feng, and D. W. Coit, “Reliability and maintenance modeling for systems subject to multiple dependent competing failure processes,” IIE Transactions, vol. 43, no. 1, pp. 12–22, 2010. View at: Publisher Site  Google Scholar
 Z. Wang, H.Z. Huang, Y. Li, and N.C. Xiao, “An approach to reliability assessment under degradation and shock process,” IEEE Transactions on Reliability, vol. 60, no. 4, pp. 852–863, 2011. View at: Publisher Site  Google Scholar
 Q. Qiu and L. Cui, “Reliability evaluation based on a dependent twostage failure process with competing failures,” Applied Mathematical Modelling, vol. 64, pp. 699–712, 2018. View at: Publisher Site  Google Scholar
 S. Song, D. W. Coit, Q. Feng, and H. Peng, “Reliability analysis for multicomponent systems subject to multiple dependent competing failure processes,” IEEE Transactions on Reliability, vol. 63, no. 1, pp. 331–345, 2014. View at: Publisher Site  Google Scholar
 K. Rafiee, Q. Feng, and D. W. Coit, “Reliability modeling for dependent competing failure processes with changing degradation rate,” IIE Transactions, vol. 46, no. 5, pp. 483–496, 2014. View at: Publisher Site  Google Scholar
 Z. An and D. Sun, “Reliability modeling for systems subject to multiple dependent competing failure processes with shock loads above a certain level,” Reliability Engineering & System Safety, vol. 157, pp. 129–138, 2017. View at: Publisher Site  Google Scholar
 S. Hao and J. Yang, “Reliability analysis for dependent competing failure processes with changing degradation rate and hard failure threshold levels,” Computers & Industrial Engineering, vol. 118, pp. 340–351, 2018. View at: Publisher Site  Google Scholar
 M. Fan, Z. Zeng, E. Zio, R. Kang, and Y. Chen, “A sequential bayesian approach for remaining useful life prediction of dependent competing failure processes,” IEEE Transactions on Reliability, vol. 68, no. 1, pp. 317–329, 2018. View at: Google Scholar
 D. Kong, C. Qin, Y. He, and L. Cui, “Sensorbased calibrations to improve reliability of systems subject to multiple dependent competing failure processes,” Reliability Engineering & System Safety, vol. 160, pp. 101–113, 2017. View at: Publisher Site  Google Scholar
 H. D. Gao, L. R. Cui, and Q. G. Qiu, “Reliability modeling for degradationshock dependence systems with multiple species of shocks,” Reliability Engineering and System Safety, vol. 185, pp. 133–143, 2019. View at: Publisher Site  Google Scholar
 L. Jiang, Q. Feng, and D. W. Coit, “Reliability and maintenance modeling for dependent competing failure processes with shifting failure thresholds,” IEEE Transactions on Reliability, vol. 61, no. 4, pp. 932–948, 2012. View at: Publisher Site  Google Scholar
 L. Jiang, Q. Feng, and D. W. Coit, “Reliability analysis for dependent failure processes and dependent failure threshold,” in Proceedings of the 2011 International Conference on Quality, Reliability, Risk, Maintenance, and Safety Engineering, pp. 30–34, Xi’an, China, June 2011. View at: Publisher Site  Google Scholar
 Q. Dong, L. Cui, and S. Si, “Reliability and availability analysis of stochastic degradation systems based on bivariate wiener processes,” Applied Mathematical Modelling, vol. 79, pp. 414–433, 2020. View at: Publisher Site  Google Scholar
 D. Kong and L. Cui, “Bayesian inference of multistage reliability for degradation systems with calibrations,” Proceedings of the Institution of Mechanical Engineers, Part O: Journal of Risk and Reliability, vol. 230, no. 1, pp. 18–33, 2016. View at: Publisher Site  Google Scholar
 Q. A. Qiu and L. R. Cui, “Gamma process based optimal mission abort policy,” Reliability Engineering and System Safety, vol. 190, Article ID 106496, 2019. View at: Publisher Site  Google Scholar
 X.S. Si, W. Wang, C.H. Hu, D.H. Zhou, and M. G. Pecht, “Remaining useful life estimation based on a nonlinear diffusion degradation process,” IEEE Transactions on Reliability, vol. 61, no. 1, pp. 50–67, 2012. View at: Publisher Site  Google Scholar
 C. J. Lu and W. O. Meeker, “Using degradation measures to estimate a timetofailure distribution,” Technometrics, vol. 35, no. 2, pp. 161–174, 1993. View at: Publisher Site  Google Scholar
 Z. Pan, N. Balakrishnan, Q. Sun, and J. Zhou, “Bivariate degradation analysis of products based on wiener processes and copulas,” Journal of Statistical Computation and Simulation, vol. 83, no. 7, pp. 1316–1329, 2013. View at: Publisher Site  Google Scholar
 I. Ntzoufras, Bayesian Modeling Using WinBUGS, John Wiley & Sons, New York City, NY, USA, 2009.
Copyright
Copyright © 2020 Chunping Li 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.