Research Article  Open Access
A Fuzzy Reliability Model of Blades to Avoid Resonance and Its Convergence Analysis
Abstract
The inherent frequencies of blades, especially with friction structure of shroud, are related to many uncertain factors, which include not only random factors but also fuzzy or other uncertain factors. In this paper, the fuzzy reliability model of blades to avoid resonance is investigated. By regarding the inherent frequency as a fuzzy random variable with trapezoidal or flatnormal membership function, a fuzzy reliability model of blades to avoid resonance is proposed based on the fuzzy cutset theory, and the corresponding numerical solution is given. Further, a theorem is proposed and proved to indicate the deficiency of the previous cutset distributions on convergence, and, then, a new cutset distribution is given to cover this shortage and ensure that the model is of good convergence. Finally, the proposed model is applied to evaluate the reliability of one certain blade. Some simulations are carried out to compare the new cutset distribution with the previous ones and study the influence of fuzzy degree on the reliability.
1. Introduction
Blade is one of the important components of aeroengine, whose failure rate is always higher than other components due to its complex service environment. The failure of blades is often related to vibration, in which resonance failure accounts for 30%~40% [1]. Therefore, it is of great significance to study the reliabilities of blades to avoid resonance.
The essence of the conventional random reliability theory is to regard the basic variables as random variables. However, in the process of vibration, the inherent frequencies of blades have not only randomness but also great fuzziness. The fuzziness is mainly due to the following. (1) The inherent frequencies of blades (especially with damping structures, such as shrouded blades) are affected by many factors, some of which have complex influence on inherent frequencies (e.g., positive pressure on contact surface has a high nonlinear effect [2]), some factors (such as circumferential position excursion of shroud to firtree serration and installation tightness between tenon and mortise) affect it irregularly, and even some unknown factors also might affect the inherent frequencies. Because of the multiple and composite influences, the inherent frequencies of blades with damping structures are of great fuzziness. (2) In the acquisition (test or simulation) process of the inherent frequencies, many steps depend on the experiences and subjective consciousness of engineers, such as the design of experiment scheme, experiment equipment, algorithm, and program. Experience and consciousness are both fuzzy concepts. These indicate that the fuzziness of inherent frequencies is inevitable and cannot be neglected. Just as the inherent frequencies of blades, the variables which include both fuzziness and randomness are called fuzzy random variables and should be handled by the fuzzy reliability theory in the reliability analysis.
The fuzzy reliability theory is based on the fuzzy set theory, which was proposed by Zadeh in 1965 [3] and had been applied in many fields [4–9]. The fuzziness in reliability problems mainly involves the fuzziness of failure state and the fuzziness of basic variables. According to the type of fuzziness in the model, the fuzzy reliability models could be classified into three groups [10–12]: the model which is only based on the fuzziness of failure state is called Profust model, the model which is only based on the fuzziness of basic variables is called Posbist model, and the model which is based on the fuzziness of both failure state and basic variables is called Posfust model. At present, the studies on the fuzzy reliability to avoid resonance are mostly the Profust models, which ignore the fuzziness of basic variable (inherent frequency) [13–16]. Therefore, this paper is aimed to study the Posbist reliability model of blades to avoid resonance.
To deal with fuzzy random variables, some methods in the existing studies (such as “stressstrength” interference model with the fuzziness of stress or strength [17–20]) are worth learning. The basic idea is to translate the fuzzy random variables into pure random variables by some defuzzification methods, and then the mature random theories could be used. The cutset method is one of the common defuzzification methods, in which it is very important to choose an appropriate cutset distribution due to its significant effects on the analysis results. At present, the common cutset distributions include uniform distribution, linear distribution, and truncated normal distribution [21–23]. However, when the cutset is longer, these distributions become unreasonable and cannot reflect the random characteristics of variables. This deficiency will bring trouble to the convergence of the model, and thus a new cutset distribution will be proposed in this paper to cover the shortage of the previous cutset distributions.
This paper is organized as follow. In Section 2, two kinds of membership functions were adopted to describe the fuzziness of the inherent frequency, and their fuzzy degrees were discussed. In Section 3, a Posbist reliability model of blades to avoid resonance was proposed based on the fuzzy cutset theory, and a numerical solution was given. In Section 4, the convergence of the Posbist model was analyzed. A convergence theorem was given and proved to indicate the deficiency of the previous cutset distributions, and a new distribution was proposed to ensure the convergence of the Posbist model. And in Section 5, simulations were carried out in order to verify the theoretical results and compare the new cutset distribution with the previous ones, and the influences of fuzzy degree on the reliability were investigated. Some conclusions were made in Section 6.
2. Description and Measurement of the Fuzzy Random Inherent Frequency
The fuzziness of fuzzy random variable can be described by a membership function , where is a set in the domain and = {the possible values of this variable}. For some , the value represents the degree of event and is called membership degree. Therefore, represents , represents , and represents whether belongs to is uncertain, and the possibility is . These show that the elements in are fuzzy and uncertain, so is not an ordinary set but a fuzzy set.
In engineering, the membership function of a fuzzy random variable generally satisfies the following format: where is called flat zone, is an increasing left continuous function, and ; is a decreasing right continuous function and . The linear or normal function is usually chosen as and in engineering.
2.1. Membership Function of Inherent Frequency
The inherent frequency is a fuzzy random variable, so the set the possible values of the inherent frequency} is a fuzzy set and could be described by the membership function . can be determined through the expression (1). Generally speaking, the inherent frequency of one blade is designed to be some value, but the dispersion is inevitable due to the manufacture dispersion. A certain frequency tolerance is allowed in engineering. In the predetermined tolerance range, the inherent frequency can be deemed to be clear. So the tolerance interval is just the flat zone, where the membership degree is the constant 1. Practically, the ultratolerance phenomenon sometimes may occur; that is to say, the inherent frequency is fuzzy within a certain range on both sides of the flat zone, where the membership degree is determined by functions and . If and are linear functions, the membership function is trapezoidal (expression (2)), and if they are normal functions, the membership function is flatnormal (expression (3)).(I)Trapezoidal membership function is as follows:
where and are fuzzy zones and and are fuzzyinformation parameters and are used to characterize the length of fuzzy zones.(II)Flatnormal membership function is as follows: where and are fuzzy zones and and are fuzzyinformation parameters and are used to adjust the shape of membership functions.
If there is no additional information, it should be taken that (). In practice, for a specific blade, the membership function of its inherent frequency can be chosen as the specific one from (2) and (3) by some methods, such as statistics, function fitting, or the experience of experts.
2.2. Fuzzy Degree of the Inherent Frequency
Fuzzy degree is used to measure the ambiguity of a fuzzy random variable through its membership function. In continuous domain, fuzzy degree can be calculated according to the following expression [24]: where . The fuzzy degrees of trapezoidal and flatnormal fuzzy random variables can be calculated according to the expression (4), and the results are shown as follows.(I)The fuzzy degree of trapezoidal membership function is as follows: (II)The fuzzy degree of Flatnormal membership function is as follows:
In summary, the fuzzy degree of a fuzzy random variable only relates to the fuzzyinformation parameters of its membership function and has nothing to do with its precise information, that is, the flat zone .
Particularly, when (or ), the membership function (2) (or (3)) will degenerate into the constant function: In this case, the set degenerates into ordinary set and the fuzzy degree of the inherent frequency will be 0 according to the expression (4), which indicates that the inherent frequency is no longer a fuzzy random variable but a pure random variable.
3. Posbist Reliability Model of Blades to Avoid Resonance
Suppose that the inherent frequency of one blade is a fuzzy random variable, whose membership function is . The excitation frequency is a random variable and follows normal distribution , whose density function is .
3.1. The Posbist Reliability Model
The reliability of this blade to avoid resonance is shown as follows:
Define and ; thus there are and . Suppose that the density function and survival function of are and , respectively and those of are and . Then,
The reliability to avoid resonance is studied based on the fuzzy cutset theory. A cut level is given first, and, then, the cutset of can be calculated as , which is an ordinary set. By introducing a random density function into the cutset , the cutset reliability can be obtained based on the conventional random reliability theory, where Because is the random density function and defined in the cutset , it is called cutset distribution. Taking the cut level as the integral variable and integrating into the interval , the Posbist fuzzy reliability model of blades to avoid resonance can be proposed as follows:
3.2. Calculation of the Model
According to the formulas in (10), after a group of are simulated randomly in the cutset , can be calculated by the two formulas below:
The cut level is regarded as a random variable which follows uniform distribution in , and the formula (11) can be written as . Therefore, a group of cut levels can be simulated randomly, and the corresponding cutset reliabilities can be calculated according the formulas in (12). Then, the fuzzy reliability to avoid resonance can be calculated as
The calculation results are related to the sample sizes and , and thus an iterative method can be used. The flowchart of the numerical simulation method is shown in Figure 1.
4. Convergence of the Posbist Model and a New CutSet Distribution
According to the model in (11), the cutset distribution introduced into the cutset is of significant effect on the reliability, and, thus, it is very important to choose an ideal cutset distribution. At present, the common cutset distributions are uniform distribution (UD), linear distribution (LD), and truncated normal distribution (TND). In order to point out the deficiency of these previous cutset distributions, the convergence of the Posbist model should be analyzed first.
4.1. Convergence of the Posbist Model
As is well known, the pure random variable is a special fuzzy random variable, which no longer has fuzziness. When the fuzzy random inherent frequency gradually converges to a pure random variable, the Posbist model should converge to the corresponding random model. Particularly, when the membership function of the inherent frequency degenerates into the constant function , , the inherent frequency should degenerate into a random variable which obeys the probability function in . In this case, according to the Posbist model, it can be seen that the cutsets corresponding to arbitrary cut levels are all . Thus, has nothing to do with and can be denoted by . Therefore, the cutset reliabilities to arbitrary cut levels are equal and have no relationship to and can be denoted by . Through the expression (11), it is obtained that . Because is the random reliability model where the inherent frequency follows the random density function in , the Posbist model degenerates into a random model.
In fact, should reflect the random characteristic of the inherent frequency. Therefore, it should just be the probability density function of random inherent frequency in conventional random model. In this way, the conventional random model is a special Posbist model, while the Posbist method is the development and extension of the random reliability method. Because is the limit function of when the cutset tends to , an ideal cutset distribution should well reflect the random characteristic of the inherent frequency and converge to the probability density function of the inherent frequency in the conventional random model.
4.2. Deficiency of the Previous CutSet Distributions
At present, the common is UD, LD, or TND. These cutset distributions can reflect the random characteristic of the inherent frequency when the cutset is short. But if the cutset is longer and tends, especially, to , LD and TND tend to UD, which means that the inherent frequency takes equal probability everywhere in (Figure 2). This is unreasonable and brings a deficiency to the convergence of the Posbist model, which is shown in the following lemma and theorem.
Lemma 1. If the cutset distribution is UD, LD, or TND, when the cutset of the inherent frequency converges to real field, that is, , the cutset reliability to avoid resonance will converge to the fixed value 1 and has nothing to do with the specific issues.
Proof. When , the cutset distributions LD and TND will converge to UD. Thus, it is just needed to prove UD.
For arbitrary fixed, when is long enough and tends to , there is always . Therefore, according to the formulas in (10), and can be written in the equivalent forms below:
When the cutset distribution is UD, substitute it into the above formulas and the following expression can be obtained:
In specific engineering problems, different blades will lead different parameters in the Posbist model. However, the above analysis shows that no matter what values of these parameters are, the cutset reliability will always converge to 1 when . In other words, it has nothing to do with the specific issues.
According the lemma above, the following theorem can be derived.
Theorem 2 (convergence theorem of the Posbist model under the cutset distribution UD, LD, or TND). If the cutset distribution is UD, LD, or TND, when the membership function of the inherent frequency gradually converges to , , the reliability assessment to avoid resonance using the Posbist model will converge to 1 and has nothing to do with the specific issues.
Proof. When the membership function converges to , , for an arbitrary fixed , the corresponding cutset always converges to . According to the lemma above, if the cutset distribution is UD, LD, or TND, the corresponding cutset reliability . Thus, based on the formula (11), it can be obtained that and will have nothing to do with the specific problem.
The analysis in Section 4.1 shows that an ideal cutset distribution should ensure that the Posbist model converges to the conventional random model when , . However, the above theorem shows that UD, LD, and TND will make the Posbist reliability assessments converge to 1 invariably, which is obviously unreasonable. The essence for this is because UD, LD, and TND cannot reflect the random characteristic of the inherent frequency when the cutset is long enough, and, thus, a new cutset distribution should be proposed.
4.3. A New CutSet Distribution
As the analysis in Section 4.1, should converge to when . is just the probability density function of the inherent frequency in the conventional random method, and it always is assumed as a normal distribution. Therefore, a new cutset distribution called modified truncated normal distribution (MTND) is proposed in this paper, whose expression is shown as follows: where , , and . is the standard deviation coefficient. is given by the expertise and generally takes the standard deviation of the inherent frequency in the random reliability model.
According to (16), MTND will converge to normal distribution when (Figure 3), and, therefore, the reliability assessment of model (11) will converge to the result of the random model.
5. Simulation Study
The first bend frequency of a known shroud blade is designed to be about 180 Hz. The manufacture tolerance is ±5%. While the ultratolerance phenomenon may appear in practice and the acceptable ultratolerance range is ±8% in engineering. When the blade is shocked by the excitation frequency which follows random distribution , the resonance may occur. Now the reliability to avoid resonance will be estimated.
From the known conditions, the median and length of the flat zone of the first bend frequency are 180 and 18, respectively; that is, . The trapezoidal and flatnormal membership functions are considered. For the former, the lengths of left and right fuzzy zones are both . So the value can be taken as . For the latter, [25]. The two membership functions of the inherent frequency are shown in Figure 4. Reliabilities to avoid resonance can be calculated as 0.9926 and 0.9929, respectively, for these two membership functions by the Posbist model proposed.
5.1. Comparison of MTND with the Previous CutSet Distributions (UD, LD, and TND)
In order to study the convergence of the Posbist model under different cutset distributions, the reliability variations with the fuzzyinformation parameters and are calculated and shown in Tables 1 and 2 and Figures 5 and 6. It is shown in Tables 1 and 2 that when the membership function of the inherent frequency converges to , ( or ), the reliability assessments under UD, LD, and TND always converge to 1, which is consistent with the theorem in Section 4.2. Meanwhile, it is shown in Figures 5 and 6 that the reliability assessments under MTND converge to 0.9908, which is just the result of the random reliability model.


In addition, it is also shown in Figures 5 and 6 that the reliability assessments under MTND decreases with the increasing of and the decreasing of . According to the formulas (5) and (6), when the fuzzy information flank is symmetrical, the higher the fuzzy degree of the inherent frequency is, the smaller the reliability will be and vice versa. This is consistent with the objective reality. However, as shown in Tables 1 and 2, with the increasing of the fuzzy degree, the reliability assessments under UD, LD and TND all decrease first and then increase instead, which is another evidence to show the deficiency of UD, LD, and TND on convergence property.
Therefore, the proposed MTND should have better performance than UD, LD, and TND.
5.2. Regularity of Fuzzy Reliability with the Variation of FuzzyInformation Parameters
In practice, if the engineers only know that the inherent frequency is about 180 Hz and have no more information, the fuzzy information flank of the flat zone is symmetrical; that is, (). Conversely, if there is additional information to ascertain that it is more possible that “the inherent frequency is smaller (or bigger) than 180 Hz,” the fuzzy information flank is asymmetrical; that is, (). For the symmetrical case, it is shown in Figures 5 and 6 that the higher the fuzzy degree of the inherent frequency is, the smaller the fuzzy reliability will be and vice versa. Thus, numerical simulations will be carried out with the asymmetrical case in Section 5.2.
The variation surface of fuzzy reliabilities with the fuzzyinformation parameters and is shown in Figure 7 and that with and is shown in Figure 8. The variation curve families of fuzzy reliabilities with each parameter are shown in Figure 9–12.
It is shown in Figures 7, 9, and 10 that the fuzzy reliability increases with the increasing of , for given , and decreases with the increasing of for given . Similarly, Figures 8, 11, and 12 show that the fuzzy reliability decreases with the increasing of , for given , and increases with the increasing of for given . In a word, for the given rightside fuzzy information, the fuzzy reliability will increase with the increasing of the fuzzy degree of leftside fuzzy information and vice versa. While, for the given leftside fuzzy information, the fuzzy reliability will decrease with the increasing of the fuzzy degree of rightside fuzzy information and vice versa. This is because of the fact that when the leftside ambiguity increases, the possibility of “the inherent frequency in the left fuzzy zone” will increase. While the mean of the excitation frequency is on right side of the flat zone , thus, the possibility that “the inherent frequency away from the excitation frequency” increases; that is to say, the reliability increases and vice versa. Similarly, the regularity will be contrary when the mean of the excitation frequency is on the left side of the flat zone.
5.3. Influences of Other Factors on the Fuzzy Reliability
Fuzzyinformation parameters and fuzzy degree can only reflect the influence of fuzzy information on reliability. But, besides the fuzzy information, the fuzzy reliability is also related to precise information, that is, the flat zone. Therefore, the conclusions in Section 5.2 are made based on the premise that the flat zone is fixed. Eight experiments are carried out to show the influence of the flat zone on reliability in Table 3.

It is shown in Table 3 that the change of the position and length of the flat zone will affect the fuzzy reliability even for the same fuzzyinformation parameters (or fuzzy degree): (1) comparing ② with ③ (or ⑥ with ⑦), it can be seen that although the fuzzy degree increases, the reliability does not decrease. Instead, it increases due to the shortening of the flat zone. (2) Comparing ③ with ⑤ (or ④ with ⑥), it can be seen that although the fuzzy degree decreases, the reliability does not increase. Instead, it decreases due to the position change of the flat zone.
6. Conclusions
The inherent frequencies of engine blades with damping structures often have a lot of fuzziness. In this paper, the inherent frequencies of blades are regarded as fuzzy random variables. The Posbist fuzzy reliability model of blades to avoid resonance is proposed based on the fuzzy cutset theory. The cutset distribution is of a significant effect on the Posbist reliability assessment result. A convergence theorem is proposed and proved to show the deficiency of the previous commonly used distributions UD, LD, and TND, and a new distribution MTND is further presented. Simulation results show that MTND is of better performance to ensure a good convergence of the Posbist model. When the fuzzy information flank is symmetrical, the bigger the fuzzy degree of the inherent frequency is, the smaller the fuzzy reliability will be and vice versa, while when the fuzzy information flank is asymmetrical, the bigger the fuzzy degree on the side near the mean of the excitation frequency is, the smaller the fuzzy reliability will be, and the bigger the fuzzy degree on the side away from the mean of the excitation frequency is, the bigger the fuzzy reliability will be.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work has been supported by the National Natural Science Foundation of China (nos. 70471057 and 71171164), and the Natural Science Foundation of the Education Department of Shaanxi Province (no. 03JK065). The authors express their sincere appreciation to Dr. Yang Yang and all the experts for their useful comments and help and are very grateful to the reviewers and the editor for their helpful comments and suggestions which have helped them in improving the quality of this paper.
References
 X. Jin, D. Gao, X. Cai, and L. Wu, “Analysis of integrated centrifugal impeller blade vibration reliability,” Journal of Aerospace Power, vol. 19, no. 5, pp. 610–613, 2004. View at: Google Scholar
 H. Y. Zhang, The Shroud Fan Blade Design for Reducing the Vibration, Chinese Aeronautical Establishment, Beijing, China, 2011.
 L. A. Zadeh, “Fuzzy sets,” Information and Computation, vol. 8, no. 4, pp. 338–353, 1965. View at: Google Scholar  MathSciNet
 V. Ravi, P. J. Reddy, and H. Zimmermann, “Fuzzy global optimization of complex system reliability,” IEEE Transactions on Fuzzy Systems, vol. 8, no. 3, pp. 241–248, 2000. View at: Publisher Site  Google Scholar
 F. Biondini, F. Bontempi, and P. G. Malerba, “Fuzzy reliability analysis of concrete structures,” Computers and Structures, vol. 82, no. 1314, pp. 1033–1052, 2004. View at: Publisher Site  Google Scholar
 H. Wu, “Bayesian system reliability assessment under fuzzy environments,” Reliability Engineering and System Safety, vol. 83, no. 3, pp. 277–286, 2004. View at: Publisher Site  Google Scholar
 S. S. Halilčević, F. Gubina, and A. F. Gubina, “The composite fuzzy reliability index of power systems,” Engineering Applications of Artificial Intelligence, vol. 24, no. 6, pp. 1026–1034, 2011. View at: Publisher Site  Google Scholar
 M. Tanrioven, Q. H. Wu, D. R. Turner, C. Kocatepe, and J. Wang, “A new approach to realtime reliability analysis of transmission system using fuzzy Markov model,” International Journal of Electrical Power and Energy System, vol. 26, no. 10, pp. 821–832, 2004. View at: Publisher Site  Google Scholar
 H. Z. Huang, M. J. Zuo, and Z. Q. Sun, “Bayesian reliability analysis for fuzzy lifetime data,” Fuzzy Sets and Systems, vol. 157, no. 12, pp. 1674–1686, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 K. Y. Cai, C. Y. Wen, and M. L. Zhang, “Survivability index for CCNs: a measure of fuzzy reliability,” Reliability Engineering and System Safety, vol. 33, no. 1, pp. 71–99, 1991. View at: Publisher Site  Google Scholar
 K. Y. Cai, C. Y. Wen, and M. L. Zhang, “Fuzzy variables as a basis for a theory of fuzzy reliability in the possibility context,” Fuzzy Sets and Systems, vol. 42, no. 2, pp. 145–172, 1991. View at: Publisher Site  Google Scholar  MathSciNet
 K. Y. Cai, C. Y. Wen, and M. L. Zhang, “Posbist reliability behavior of typical systems with two types of failure,” Fuzzy Sets and Systems, vol. 43, no. 1, pp. 17–32, 1991. View at: Publisher Site  Google Scholar  MathSciNet
 K. J. Xu, L. P. Jiang, and Y. S. Sui, “Applying fuzzy method to assessment vibration reliability of structural components,” Journal of Mechanical Strength, vol. 25, no. 5, pp. 523–526, 2003. View at: Google Scholar
 L. Ma, H. X. Zhang, B. Xu et al., “Fuzzy estimation of aeroengine blades operational vibration reliability,” Mechanical Science and Technology, vol. 11, no. 25, pp. 1365–1368, 2006. View at: Google Scholar
 L. Ma, J. Xie, B. Xu et al., “Study on blade vibration reliability under working environment,” Turbine Technology, vol. 49, no. 1, pp. 40–42, 2007. View at: Google Scholar
 F. D. Pi, Study on the Structure Vibration Fuzzy Reliability Analysis Method and Its Application, Northwestern Polytechnical University, Xi'an, China, 2006.
 J. Sun, Study on the Fuzzy Random Reliability Analysis Methods, Northwestern Polytechnical University, Xi'an, China, 2005.
 H. Z. Huang, “Structural reliability analysis using fuzzy sets theory,” Maintenance and Reliability, vol. 14, no. 4, pp. 284–294, 2012. View at: Google Scholar
 H. Huang, Z. Sun, D. Guo, L. Li, F. Zhou, and Z. Tian, “Fuzzy reliability analysis in the case of random stress and fuzzy strength,” Journal of Mechanical Strength, vol. 23, no. 3, pp. 305–307, 2001. View at: Google Scholar
 Y. G. Dong, “Fuzzy reliability design with random variable and fuzzy variable,” Chinese Journal of Mechanical Engineering, vol. 36, no. 6, pp. 25–29, 2000. View at: Google Scholar
 Y. G. Dong, W. Y. Zhu, and X. Z. Chen, “Study on a calculating method of machine fuzzy reliability,” Journal of Systems Engineering, vol. 15, no. 1, pp. 7–12, 2000. View at: Google Scholar
 Q. M. Jiang and C. H. Chen, “A numerical algorithm of fuzzy reliability,” Reliability Engineering and System Safety, vol. 80, no. 3, pp. 299–307, 2003. View at: Publisher Site  Google Scholar
 Z. Lu, J. Sun, and Y. Xu, “Numerical algorithm for fuzzy reliability analysis of mechanical structure system,” Chinese Journal of Mechanical Engineering, vol. 41, no. 9, pp. 19–23, 2005. View at: Google Scholar
 S. L. Chen, J. G. Li, and X. G. Wang, Fuzzy Set Theory and Application, Science Press, Beijing, China, 1st edition, 2005.
 H. He and Z. Lu, “Reliability analysis method for structure with fuzzy variable,” Journal of Mechanical Strength, vol. 31, no. 4, pp. 609–614, 2009. View at: Google Scholar
Copyright
Copyright © 2014 Meng Zhang 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.