Research Article | Open Access
Wenting Liu, Qingliang Zeng, Lirong Wan, Chenglong Wang, "A Comprehensive Method of Apportioning Reliability Goals for New Product of Hydraulic Excavator", Mathematical Problems in Engineering, vol. 2019, Article ID 5030187, 11 pages, 2019. https://doi.org/10.1155/2019/5030187
A Comprehensive Method of Apportioning Reliability Goals for New Product of Hydraulic Excavator
It is important to allocate a reliability goal for the hydraulic excavator in the early design stage of the new system. There are some effective methods for setting reliability target and allocating its constituent subsystems in the field of aerospace, electric, vehicles, railways, or chemical system, but until now there is no effective method for the hydraulic excavator or engineering machinery. In this paper, an approach is proposed which combines with the conventional reliability allocation methods for setting reliability goals and allocating the subsystem and parts useful in the early design stage of the hydraulic excavator newly developed. It includes Weibull analysis method, modified Aeronautical Radio Inc. (ARINC) method, and modified systematic failure mode and effect analysis (FMEA) method. After completing reliability allocation, it is necessary to organize the designers and experts to evaluate the rationality of the reliability target through FEMA analysis considering feasibility of the improvement technically for the part which was new developed or had fault in its predecessor. The proposed approach provides an easy methodology for allocate a practical reliability goal for the hydraulic excavator capturing the real life behavior of the product. It proposes a simple and unique way to capture the improvement of the subsystems or components of the hydraulic excavator. The proposed approach could be extended to consider other construction machinery equipment and have practicality value to research excellent mechanical product.
Reliability goals setting and allocation are important tasks in the early design stage of the hydraulic excavator newly developed; it is essential to allocate scientific reliability target and then apportion it to the subsystems or components of the hydraulic excavator. The advantages of reliability and functions of the hydraulic excavator are key factors for enhancing core competitiveness to a manufacturing enterprise. In recent years, attention to excavator’s system reliability has risen for the increase in sophistication of engineering systems used in high technology industry processes . It is a difficult task to set reasonable and realizable reliability goals for excavator which has parallel system. The historical frequency of failures could be measured quantitatively but it is more difficult to quantify the cost sensitivity to the reliability level .
Reliability allocation is a top-down method to apportion accuracy goals in a system. Several reliability allocation methods are proposed by many researchers, including Aeronautical Radio Inc. (ARINC) method ; Advisory Group of Reliability of Electronic Equipment (AGREE) method ; Feasibility of Objectives (FOO) method which was mentioned in MIL-HDBK-338B ; The Boyd method which combines equal method and ARINC method ; Critical Flow Method (CFM) ; Analytic Critical Flow Method (ACFM) ; Karmiol method ; Bracha method ; Integrated Factors Method (IFM) [11–14]; Analytic IFM method (A-IFM) [15, 16]; The Maximum Entropy Ordered Weighted Averaging (ME-OWA) method [17–19]; and the method using risk priority number (RPN) generated during failure modes and effect analysis (FMEA) to link reliability allocation [20–25], each with its advantages and disadvantages, but none of them is suitable for the reliability allocation of the hydraulic excavators completely. Some conventional reliability allocation methods do not consider the reliability improvement potential except the FEMA method; corrective actions are defined and implemented by identifying potential problems and calculating the risk to eliminate or reduce their occurrence possibility [26–28]. The existing methods as previously mentioned for setting achievable reliability goals to each subsystem or component are characterized by its own advantages and disadvantages. The methods were summarized and introduced in these papers [7, 13, 15] in great detail. Table 1 shows the advantages and disadvantages of the literature methods which could be applied to the hydraulic excavator.
Parallel-series systems hydraulic excavators have hundreds of components; project progress will be prolonged if each one is analyzed with the FMEA method. It is scientific to combine several methods by utilizing their advantages and remedying their disadvantages that not only conform to the actual situation, but also increase the reliability of hydraulic excavators effectively.
2. Analysis of Reliability Allocation Method
If sufficient historical failure information from previous design was available, the Aeronautical Radio Inc. (ARINC) method could be used to estimate the failure rate of subsystems easily .
The Aeronautical Radio Inc. (ARINC) is an apportionment technique method which is based on the failure rates of units or subsystems . It can be used if database or previous experience on similar subsystem, allowing predicting failure rate of the subsystem above, is expressed as
where is the failure rate of unit I and is a reliability allocation weight for unit i, for .
This method seems to be the first reliability allocation methods that consider historical information for obtaining allocation weights but not consider the improvement of the subsystems and the components.
Several papers proposed reliability allocation approaches considering the field failure and failure analysis information generated during historical warranty data. The modified failure mode and effect analysis (FMEA) method which uses the risk priority number (RPN) was defined as the multiplication of the severity ranking and occurrence ranking of a failure modes in a subsystem that makes it good for verifying the feasibility of the reliability target of the subsystems or the components .
FMEA was developed in the 1960s by the aerospace industry. It was a systematic method for identifying different failure modes and their effects occurring in the system. A few recent papers [20–22, 26–30] use RPN information generated during the criticality analysis of FMEA considering the failure effect in reliability allocation. The RPN of failure mode j in subsystem i was given byThen, the reliability allocation weight is given by where and where ’s are severity ranking, an evaluation of how serious the effect would be if a given failure modes did occur, ’s are occurrence ranking, the failure frequency, or the qualitative failure rate, and ’s are detection ranking, an assessment of an ability of the failure being detected.
In a mission critical system, however, the detection ranking is not included in the RPN calculation because the detection of a failure mode is considered when assigning a severity ranking. Thus, the RPN for a mission critical system is given by where
Yadav proposed a reliability allocation method considering failure criticality and functional dependency . The reliability allocation weight for subsystem i is given as where and denote functional dependency and criticality of subsystem i, whereas and represent the importance assigned to dependency and criticality, factors respectively, denote the functions supported by subsystem i, and , indicate the severity and occurrence rankings for failure mode j of function k supported by subsystem i.
Then Kyungmee O. Kim proposed a new approach  that suggests exponential transformation of the original 10 point ordinal severity rating to reduce severe failure effect effectively. Therefore, the reliability weight is given as where is transformed severity and is failure frequency of subsystem and is the number of failure modes in subsystem having maximum severity rating.
Because these ratings are follow linearity assumption in rating scale, a general understanding among design community is easy to think that improving reliability of a poor subsystem than improving reliability of a highly reliable subsystem. And sometimes lower reliability is assigned to a subsystem already having lower failure frequency.
Yadav proposed a more realistic and effective reliability approach  based on nonlinear phenomena in severity rating and failure rate; the severity of subsystem i has multiple failure mode defined as
where , .
The reliability allocation method ensures each subsystem getting improvement target is assigned in proportion of potential for improvement and exposes technological difficulty in achieving improvement goals when negative failure rate target is assigned.
While calculating RPN, although these methods have been considering potential for subsystems or components improvement, technology difficulty and complexity of a given subsystem, to some extent, allocate reliability in accordance with degrees of difficulty in achieving reliability target; however, if the improvement effort is effective in reality, the kinds of measures that could be taken for achieving the reliability targets were not taken into account. So it is imperative to develop practical reliability allocation method for hydraulic excavators considering effective measures in accordance with degree of difficulty for improvement effort.
3. Proposed Comprehensive Methods for Reliability Allocation
In this section, based on the database of operation and maintenance of the existing hydraulic excavators, the failure rate is used for the reliability evaluation indexes; a comprehensive allocation method is proposed to perform the reliability allocation for the hydraulic excavator newly developed which considers effective measures for reliability improvement to allocate reliability. Decision-making procedure used in this study is shown in Figure 1.
3.1. Perform Reliability Goal for New System
Weibull distribution [32–34] and modifications of the Weibull distribution [35, 36] for reliability analysis were discussed by several authors. Many researchers have suggested the lifetime distribution of the hydraulic excavator follows the Weibull distribution [37–40]; they stated above use Weibull distribution as theoretical reliability model, according to the experimental data, analyzed the reliability of the subsystems for the excavator. So the reliability goal of the new systems in this study is gained based on Weibull analysis. Assume the lifetime of the hydraulic excavator follows two-parameter Weibull distribution.
Failure frequency function for Weibull distribution is given as follows:Reliability function for Weibull distribution is given as follows:And unreliability function or cumulative distribution function is given asAnd failure rate function is given aswhere is shape parameter and is scale parameter.
3.2. Reliability Allocation
The failure data of subsystems for existing hydraulic excavators should be utilized in design of new hydraulic excavator. The subsystem or component which has a low failure frequency is also reliable if unchanged in the new design, and it is natural to allocate a high reliability to these subsystems.
Based on the maintenance data of existing hydraulic excavators, the failure frequency factor of the subsystem or components is given as follows:
where represents the subsystem or components failure number and represents the total number of failures.
We allocated the subsystems and the components which are newly designed or have recorded faults by the modified ARINC method as in (14).
3.3. Amend the Reliability Target considering the RPN of the Improvements
Commonly, the design community thought that the subsystems or components with higher failure rate must be assigned lower reliability target and with lower failure rate must be assigned higher reliability target in new corresponding subsystems, but reliability target of the subsystems or components is disproportionate to the failure rate. In the hypothesis that linear or nonlinear relationship between failure rate and improvement effort does not usually conform to the reality, the key is that if the design community can take effective measures to reduce the failure rate or increase Mean Time Between Failure (MTBF), higher subsystems or components reliability target should be assigned; otherwise, it is adjusted to lower reliability level.
The experts are organized to evaluate the allocation reliability target in Section 3.2 term by term and provide effective measures for subsystems or components. Considering the restriction by intricacy, state of the art, maintainability, manufacturing technology, working conditions, and cost, reliability target for subsystems or components should be amended by the modified FEMA method as shown in Table 4 when an effective measures are provided; otherwise, we follow the original reliability target.
Table 2 describes the severity ranking and linguistic terms used to describe the failure effect. Table 3 presents the linguistic terms for the occurrence ranking and the corresponding quantitative failure rate; both of them are calculated through market share and failure rate of one type of particular brand hydraulic excavators and could be revised according to the actual applying.
The allocation weight of the subsystems or components that have effective measures is recalculated bySo the upgrade rate of the subsystems or components is defined asReplace the subsystem or component reliability target with
3.4. Applicability of the Methodology
The proposed approach was applied to the hydraulic excavator system in the early design stage for new product development based on the ARINC method and FEMA method. It is a simple and effective approach which uses failure and maintenance data of the previous generation. The most important is that it provides a realistic reliability allocation method considering the improvement of the parts and components and can be used in different industries and filed.
3.5. Execution Steps
The proposed approach for reliability allocation applied to the hydraulic excavator is structured in the following steps.
Step 1. The following 6 units of the hydraulic excavator for facilitating subsequent improvement are pointed out: Unit 1: hydraulic subsystem. Unit 2: electric subsystem. Unit 3: power subsystem. Unit 4: working device subsystem. Unit 5: chassis subsystem. Unit 6: body and accessories subsystem.
Step 2. Data was collected over a period of 2 or 3 years and sorted out according to the subsystems and parts which have been pointed out in Step 1. Because the hydraulic excavators guarantee period is 2000~3000h, working hours ≤3000h are chosen.
Step 3. The reliability index λ is gained based on the statistic data, as shown in Table 5.
Step 4. The reliability target (failure rate) is allocated to the subsystems in proportion to allocation weights () by the modified ARINC method.
Step 5. Reliability target for subsystems or components will be amended by the modified FEMA method, replace the subsystem or component reliability target with if the target allocated in Step 4 is hard to achieve.
4. Illustrative Example
In this section the proposed approach is applied to the reliability goal setting and allocation for the hydraulic excavator SC5532 in the early design stage to illustrate its feasibility and applicability. Furthermore, an analysis of the proposed methodology is subsequently performed to validate its effectiveness.
The case study will be discussed involving a period of 2-year maintenance data of a type of SC5532 belonging to the company LISHIDE to illustrate the modeling approach proposed in Section 3. All of the failures happening time as shown in Table 6 counted 146 times for 151 excavators. Then the f(t), F(t), R(t), and λ(t) were calculated range from 0 to 3000 hours and the interval time was 300 hours, as shown in Table 7. The total failure rate of the current hydraulic excavator systems at 3000 hours was calculated as λ = 0.01383012. Now suppose we want to improve the failure rate by 20%, so that the system failure rate goal was given by = 0. 0.01106409.
The failures times for failure components of the hydraulic excavator subsystems and the weighs were calculated with (14), as shown in Table 8. The fault-free components with failure rate = 0 were excluded in Table 8.
The reliability target was allocated to the parts in the new system which had fault recording in existing hydraulic excavators.
After reliability allocation, the reliability targets need to be verified considering the improvement of the components by the designers or the experts who were responsible for subsystems or components. In this case the hydraulic cylinders were selected to verify the rationality of the reliability target. The boom cylinders, arm cylinders, and bucket cylinders will be analyzed together for they are working on the same principle and with the same failure modes. The designers and experts were organized to do FEMA analysis and put forward the improving measures and predict , of each failure mode for cylinders as shown in Table 9.
= 0.001515629 was obtained from Table 8, and = 0.001016059; then the reliability target of the cylinder was replaced with . Similar analysis as described above could be used for all the other components of the hydraulic excavators.
4.2. Comparison and Discussion
For further validation of the effectiveness of the approach proposed in this article, the comparative analysis is performed with the conventional FEMA. It must be a tedious process to do reliability goal setting and allocation to the parts or components of the hydraulic excavator which have a large, complex system by the conventional FEMA method. So the modified ARINC method here is used to extract the parts or components which can be improved as the reliability target by statistics analysis of the historical failures and maintenance data from the previous generation. As shown in Table 8, the components’ preliminary reliability target was set, respectively, by the modified ARINC method. But in this step, the same reliability target is allocated to the different components with the same failure times without considering the improvement of them. That is, by using the traditional FMEA, different sets of O, S ratings can yield the same RPN value, which entails that these components should be given the same reliability target. Since the restriction by intricacy, state of the art, maintainability, manufacturing technology, working conditions, and cost, the reliability target could be higher or lower than allocated. Then the modified FEMA method could be used to amend the components’ reliability target by the value considering the improvement measures of them. The proposed approach has some advantages, such as the following:(1)The proposed approach is effective for providing more meaningful information for risk management decision made based on the modified ARNIC method in many practical applications.(2)The proposed approach fuses personal judgments of FEMA team members into group assessments; it is able to relieve the influence of unfair opinions on the risk analysis result.(3)Using the proposed approach we are able to realize the reliability target allocated to the components or the parts taking into account the improvement measures, both subjective and objective weights could be considered in the reliability allocation of the hydraulic excavator systems.
5. Conclusion and Future Work
In this paper, we described an approach based on ARINC method and FEMA method for reliability target goal set and allocated the parts and components of the hydraulic excavator newly developed taking into account the reliability improvement. It is worth mentioning the following:(1)The reliability goal of the hydraulic excavator newly developed have been set utilizing the indicator failure rate = 0.01106409/3.0·103[h] with analytical methods based on the failure data which obeys the Weibull (2P) distribution.(2)For testing the validity of the reliability target considering the improvement of the components, the designers or the experts responsible for subsystem or components were organized to do FEMA analysis of the components. The reliability target should be revised if there are effective measures for improvement.(3)Increasing the reliability of the hydraulic excavator is a full life-cycle activity, including the process of designing, experimental testing, technology route planning, parts purchasing, production, and maintenance. So it is important to control and prevent failures of the components and subsystems at each stage of development to achieve the reliability goal and to improve the reliability of the hydraulic excavator new products after finishing the reliability allocation.(4)The proposed approach in this paper could be used not only to allocate reliability for hydraulic excavators, but also for other construction machinery equipment after a few modifications.
In the future, the research could focus on developing a software tool to make the proposed FEMA approach simple and convenient to help reliability allocation more efficient for enhancing hydraulic excavators’ system reliability and safety.
The data used to support the findings of this study have not been made available because we have used third party data; it belongs to the LISHIDE Company which is mentioned in our manuscript; we do not own the data, and for the privacy and commercial confidentiality of the company, we cannot deposit them in a publicly available data repository. But we make sure all the data supporting the conclusions of the study is real and effective. The readers who want to obtain the source of the data could contact the company technical director: Qingxiang Huan (Email: email@example.com).
Conflicts of Interest
The authors declare that there are no conflicts of interest.
This work was supported by the National Natural Science Foundation of China (Grant no. 51674155) and Innovative Team Development Project of Ministry of Education (Grant no. IRT_16R45).
- O. P. Yadav and X. Zhuang, “A practical reliability allocation method considering modified criticality factors,” Reliability Engineering & System Safety, vol. 129, pp. 57–65, 2014.
- M. Feizabadi and A. E. Jahromi, “A new model for reliability optimization of series-parallel systems with non-homogeneous components,” Reliability Engineering & System Safety, vol. 157, pp. 101–112, 2017.
- W. H. Alven, Reliability Engineering Prepared by ARINC Research Corporation, Prentice Hall, Englewood Cliffs, NJ, USA, 1964.
- Advisory Group of Reliability of Electronic Equipment (AGREE), Reliability of Military Electronic Equipment, Office of the Assistant Secretary of Defense Research and Engineering, Washington, DC, USA, 1957.
- R. T. Anderson, “Reliability Design Handbook,” Defense Technical Information Center, ITT Research Institute, Chicago, Ill, USA, 1976.
- J. Boyd, “Allocation of reliability requirements: a new approach,” in Proceedings of the Reliability and Maintainability Symposium, pp. 5-6, Las Vegas, Nev, USA, January 1992.
- G. Di Bona, A. Forcina, and A. Silvestri, “Critical flow method: a new reliability allocation approach for a thermonuclear system,” Quality and Reliability Engineering International, vol. 32, no. 5, pp. 1677–1691, 2016.
- G. Di Bona and A. Forcina, “Analytic critical flow method (acfm): a reliability allocation method based on analytic hierarchy process,” Journal of Failure Analysis and Prevention (JFAP), vol. 17, no. 6, pp. 1149–1163, 2017.
- E. D. Karmiol, Reliability Apportionment. Preliminary Report EIAM 5, General Electric, Schenectady, NY, USA, 1965.
- V. J. Bracha, “The methods of reliability engineering,” Machine Design, vol. 7, pp. 70–76, 1964.
- D. Falcone, F. De Felice, G. Di Bona, and A. Silvestri, “Analysis in a sintering plant by the employment of a new reliability allocation method,” Modeling and Simulation, vol. 1, pp. 383–388, 2004.
- F. D. Felice, D. Falcone, A. Silvestri, and G. D. Bona, “Proposal of anew reliability allocation methodology: the integrated factors method,” International Journal of Operations and Quantitative Management, vol. 16, no. 1, pp. 67–85, 210.
- G. D. Bona, V. Duraccio, A. Silvestri et al., “Validation and application of a reliability allocation technique (advanced integrated factors method) to an industrial system,” in Proceedings of the IASTED International Conference on Modelling, Identification and Control, (MIC 2014), pp. 75–79, Austria, February 2014.
- G. Di Bona, D. Falcone, A. Silvestri, and A. Forcina, “IFM target 2.0: an innovative method to define reliability target for prototype systems,” The International Journal of Advanced Manufacturing Technology, vol. 95, pp. 3349–3367, 2018.
- G. Di Bona, A. Forcina, A. Petrillo, F. De Felice, and A. Silvestri, “A-IFM reliability allocation model based on multicriteria approach,” International Journal of Quality & Reliability Management, vol. 33, no. 5, pp. 676–698, 2016.
- G. Di Bona, A. Silvestri, A. Forcina, and D. Falcone, “AHP-IFM target: an innovative method to define reliability target in an aerospace prototype based on analytic hierarchy process,” Quality and Reliability Engineering International, vol. 33, no. 8, pp. 1731–1751, 2017.
- Y. Chang, K. Chang, and C. Liaw, “Innovative reliability allocation using the maximal entropy ordered weighted averaging method,” Computers & Industrial Engineering, vol. 57, no. 4, pp. 1274–1281, 2009.
- C.-S. Liaw, Y.-C. Chang, K.-H. Chang, and T.-Y. Chang, “ME-OWA based DEMATEL reliability apportionment method,” Expert Systems with Applications, vol. 38, no. 8, pp. 9713–9723, 2011.
- R. Fullér and P. Majlender, “An analytic approach for obtaining maximal entropy OWA operator weights,” Fuzzy Sets and Systems, vol. 124, no. 1, pp. 53–57, 2001.
- S. Yousefi, A. Alizadeh, J. Hayati, and M. Baghery, “HSE risk prioritization using robust DEA-FMEA approach with undesirable outputs: a study of automotive parts industry in Iran,” Safety Science, vol. 102, pp. 144–158, 2018.
- K. O. Kim, Y. Yang, and M. J. Zuo, “A new reliability weight for reducing the occurrence of severe failure effects,” Reliability Engineering & System Safety, vol. 117, pp. 81–88, 2013.
- K. O. Kim and M. J. Zuo, “General model for the risk priority number in failure mode and effects analysis,” Reliability Engineering & System Safety, vol. 169, pp. 321–329, 2018.
- Y. Wang, R. C. M. Yam, M. J. Zuo, and P. Tse, “A comprehensive reliability allocation method for design of CNC lathes,” Reliability Engineering & System Safety, vol. 72, no. 3, pp. 247–252, 2001.
- O. P. Yadav, N. Singh, and P. S. Goel, “Reliability demonstration test planning: a three dimensional consideration,” Reliability Engineering & System Safety, vol. 91, no. 8, pp. 882–893, 2006.
- R. R. Itabashi-Campbell and O. P. Yadav, “Gauging quality and reliability assurance efforts in product development process,” International Journal of Process Management and Benchmarking, vol. 2, no. 3, pp. 221–233, 2008.
- H. C. Liu, L. Liu, and N. Liu, “Risk evaluation approaches in failure mode and effects analysis: a literature review,” Expert Systems with Applications, vol. 40, no. 2, pp. 828–838, 2013.
- H.-W. Lo and J. J. H. Liou, “A novel multiple-criteria decision-making-based FMEA model for risk assessment,” Applied Soft Computing Journal, vol. 73, pp. 684–696, 2018.
- Y. Du, X. Lu, X. Su, Y. Hu, and Y. Deng, “New failure mode and effects analysis: an evidential downscaling method,” Quality and Reliability Engineering International, vol. 32, no. 2, pp. 737–746, 2016.
- S. Hu, Q. Fang, H. Xia, and Y. Xi, “Formal safety based on relative risks model in ship navigation,” Reliability Engineering & System Safety, vol. 92, no. 3, pp. 369–377, 2007.
- Department of the Army, TM 5-689-4. Failure Modes, Effects and Criticality Analysis (FMECA) for Command, Control, Communications, Computer, Intelligence, Surveillance, and Reconnaissance (C4ISR) Facilities, 2006.
- O. P. Yadav, “System reliability allocation methodology based on three-dimensional analyses,” International Journal of Reliability and Safety, vol. 1, no. 3, pp. 360–375, 2007.
- X. F. Wang, Fundamentals of Automotive Reliability Engineering, Tsinghua University Press, Beijing, China, 2007.
- D. N. P. Murthy, M. Bulmer, and J. A. Eccleston, “Weibull model selection for reliability modelling,” Reliability Engineering & System Safety, vol. 86, no. 3, pp. 257–267, 2004.
- S. J. Almalki and S. Nadarajah, “Modifications of the Weibull distribution: a review,” Reliability Engineering & System Safety, vol. 124, pp. 32–55, 2014.
- J. Yazhou, W. Molin, and J. Zhixin, “Probability distribution of machining center failures,” Reliability Engineering & System Safety, vol. 50, no. 1, pp. 121–125, 1995.
- P. P. Guilani, P. Azimi, S. T. A. Niaki, and S. A. A. Niaki, “Redundancy allocation problem of a system with increasing failure rates of components based on Weibull distribution: A simulation-based optimization approach,” Reliability Engineering & System Safety, vol. 152, pp. 187–196, 2016.
- R. Jiang, M. Zuo, and D. N. P. Murthy, “Two sectional models involving two Weibull distributions,” International Journal of Reliability, vol. 6, pp. 103–122, 1999.
- A. t. Pratama, A. O. Moeis, and Farizal, “Pratama A t, Moeis A O, Farizal. Effort on reliability and availability improvement of excavator HITACHI ex-2500 at open-pit coal mine areas by using importance measure,” in Proceedings of the International Seminar on Industrial Engineering and Management, 2011.
- P. Dašić, M. Živković, and M. Karić, “Reliability analysis of the hydro-system of excavator SchRs 800 using Weibull distribution,” Applied Mechanics and Materials, vol. 806, pp. 173–180, 2015.
- S. J. Back, Y. K. Son, J. H. Kim, and J. C. Lee, “Reliability demonstration of a robot system for an unmanned controlling excavator,” Applied Mechanics and Materials, vol. 563, pp. 215–218, 2014.
Copyright © 2019 Wenting Liu 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.