Condition-Based Maintenance Optimization for Motorized Spindles Integrating Proportional HazardModel with SPC Charts

Reliability of motorized spindles has a great effect on the performance and productivity of computer numerical control (CNC) machine tools for intelligent manufacturing. Condition-based maintenance (CBM) is an efficient method to prevent serious failures, to improve system reliability, and to reduce management costs for motorized spindles. However, owing to various degradation features acquired during condition monitoring, the challenge is to propose an appropriate feature to evaluate the reliability level of motorized spindles and to set up optimal CBM policies. Based on the motivation, a three-stage approach is proposed in this paper. In the first stage, proportional hazard model (PHM) is developed to describe the reliability considering failure events together with multiple degradation features. Next, statistical process control (SPC) charts are constructed for condition monitoring and anomaly detection in order to achieve early detection of potential failures. At last, a CBM schedule is modeled in consideration of maintenance cost minimization; the maintenance plan is optimized by determining the optimal control limits of SPC charts.


Introduction
As an emerging technology in the field of intelligent manufacturing, the motorized spindle of the computer numerical control (CNC) machine tool is regarded as a core component of all manufacturing facilities; its performance could determine the machining accuracy and production efficiency of CNC machine tools, even the whole manufacturing supply chain, to a large extent [1]. us, its reliability improvement and health management cause widespread concern. For the purpose of improving the availability and reducing production expenses of the motorized spindle, condition-based maintenance (CBM) has drawn much attention [2]. It focuses on utilizing the condition monitoring information collected in real time to assess health status and to develop optimal maintenance strategies [3].
As a condition monitoring method, control charts of the statistical process control (SPC) theory can be incorporated with maintenance analysis to realize CBM optimization.
Shrivastava et al. [4] integrated the cumulative sum (CUSUM) control chart with preventive maintenance (PM) planning for joint optimization of maintenance activities and chart parameters. Integrated models combining the X control chart with PM actions were proposed in [5,6] in order to protect the equipment from quality decline. Lampreia et al. [7] modified CUSUM and exponentially weighted moving average (EWMA) control charts for online vibration monitoring, thus reducing false alarms and planning early interventions. e similar modelling idea was also adopted for the economic design of maintenance plans with the Shewhart control chart, chi-square control chart, etc., as discussed in [8][9][10][11][12]. Mehrafrooz and Noorossana [13] developed a model under the assumption of perfect maintenance; the model could express six alternative scenarios of the monitored process to achieve the most economic maintenance schedule. e integrated model was further developed by Yin et al. [14] considering ten scenarios of the delayed monitoring process. From the perspective of reliability improvement, some SPC charts were also established to monitor the feature of time-between-events (TBE) assuming that TBE followed the gamma distribution [15], the normal distribution [16], the exponential distribution [17], the Weibull distribution [18], etc.
However, these aforementioned control charts are all univariate, which are commonly used to monitor individual observations of a certain type of characteristic variable. For the case of the motorized spindle, condition monitoring usually contains multiple sensor signals, such as vibration, bearing preload, acoustic emission, sound, and temperature [19], and the signals are often autocorrelated [20]. erefore, the conventional SPC charts will not be effective to extract useful and reliable information from the monitoring process, leading to larger errors in CBM optimization.
To overcome this problem, the paper develops a threestage approach for the purpose of extracting the fusion feature, monitoring the health degradation, and optimizing the CBM strategy for motorized spindles. e main contributions of the proposed approach are described as follows: (1) Proportional hazard model (PHM) is introduced to extract the health indicator. e benefits can be demonstrated in two aspects: on the one hand, PHM can integrate multivariate information from multiple monitoring signals together with failure data, which results in the more accurate degradation detection for CBM optimization. On the other hand, the hazard rate from PHM is selected as the health indicator of SPC charts. It can reduce random interference for SPC charts when monitoring the autocorrelation process; the detection error, such as the false alarm or undetected failure, can be effectively avoided.
(2) CBM optimization is combined with SPC chart design based upon economic considerations, the health condition of motorized spindles can be constantly updated, and two types of the failure mechanism, including the random failure and the fatigue failure, and four possible scenarios are taken into consideration. e integrated model can avoid unnecessary or untimely maintenance and achieve considerable economic benefit by reducing maintenance costs and raising production efficiency, which provide effective service and guidance for practitioners. e remainder of the paper is organized as follows. Section 2 describes the CBM optimization problem and assumptions in detail, in which the failure mechanism, SPC charts, and maintenance policy are discussed. In Section 3, the integrated model for combined SPC chart design and CBM scheduling is developed. e hazard rate extracted from the PHM is used to characterize the health condition, and the SPC scheme considering four possible scenarios is established to monitor the degradation trend in real time.
en, Section 4 introduces the proposed three-stage approach for optimizing the CBM schedule upon economic considerations. In association with the PHM and the SPC charts, the average long-run cost during a renewal cycle is minimized by optimizing the upper control limit (UCL). Next, a real case of the motorized spindle is studied to verify the feasibility and practicability of the proposed approach in Section 5. Finally, conclusions are presented in the last section.

Problem Statement and Assumptions
For this research, the degradation of the motorized spindle during the long life cycle is due to both external factors, like cumulative shocks, environmental stress, etc., and internal factors, like continuous aging, etc. erefore, the system will suffer from random failures and fatigue failures. Multiple condition features together with failure event data can be used as the health indicator to characterize the effects of degradation.
en, the SPC chart is used as a failure diagnosis tool for monitoring the health indicator and detecting the two types of failures: the random failure, also known as the hard failure, which is self-announcing and could be detected immediately, and the fatigue failure, also known as the soft failure, which is defined to occur only when the degree of degradation reaches the preset control limits. e SPC chart scheme is shown in Figure 1, and its control limits can be usually computed following the formula [21] in Table 1 according to different SPC charts.
When combined with maintenance planning, the control limits of SPC charts can be regarded as the trigger of a maintenance action; the optimization of the CBM strategy is transformed into determining the optimal control limits in order to minimize the maintenance cost from the perspective of economic efficiency. e integrated CBM model is developed under the assumptions as follows: (1) e changing trend of the health indicator may not be obvious over time and increase rapidly after the system enters the wear-out stage. In order to grasp system degradation level and implement maintenance actions in time, inspections are scheduled every Δ time unit to detect the health condition periodically. e inspection time is short so that the inspection activities are considered instantaneous and nondestructive. e cost per inspection C i is predetermined as a constant according to the actual condition.
(2) e two types of failures correspond to two types of maintenance activities, which are the preventive replacement and the corrective replacement, respectively. To be more specific, if a hard failure occurs unexpectedly, the corrective replacement is performed instantly, and the average cost per corrective replacement C c covers all costs with respect to the failure and its associated problems. On the contrary, if the monitoring health indicator is detected to exceed the control limits at the inspection time, which is interpreted as a soft failure, the preventive replacement is carried out. e average cost per preventive replacement C p covers the costs for performance maintenance, improvement, and other related operations. In general, C c > C p since the corrective replacement is relatively more complex than the preventive replacement. (3) After each replacement, the condition of the monitored motorized spindle will be restored to the initial state, representing that the system is under the "as good as new" maintenance strategy. In addition, it is assumed that the time for both preventive and corrective replacements is not considered in the CBM model.

Condition-Based Maintenance (CBM)
Model Development e integrated model for SPC chart design and CBM scheduling is developed from an economic point of view, whose objective is to minimize the average long-run cost considering the costs of all inspections and maintenance activities.

Proportional Hazard Model (PHM). PHM proposed by
Cox [22] is a popular model which is used to describe the relationship between the degradation influence and the hazard rate. e degradation level is characterized by multiple features during lifetime operation, which are considered as covariates. e advantage of using the PHM in CBM is that it could deduce the hazard rate considering all covariates in the current state, such as the failure event and condition monitoring factors. According to the PHM, the hazard rate at time t is defined by where h 0 (t) is the baseline hazard rate at time t, and it will be continually modified by Z(t); let Z(t) � [Z 1 (t), Z 2 (t), · · · Z n (t)] be a vector made up of values of n condition monitoring features at time t, which is a time-dependent stochastic process representing the degradation evolution. And c is a vector made up of regression coefficients, which determines the influence weight of the monitoring features on the failure rate. e time to failure of the electromechanical system commonly follows the Weibull distribution [23]. As for motorized spindles, their baseline hazard rate of the PHM is assumed to be derived from a two-parameter Weibull distribution, which has the expression as where β > 0 represents the shape parameter of the Weibull distribution and η > 0 represents the scale parameter of the Weibull distribution. Substituting the Weibull baseline hazard rate function, Z(t), and c vectors into the PHM, Weibull proportional hazard model (WPHM) is derived as Monitored variable at t μ 0 e mean of in-control data σ 0 e standard deviation of in-control data n Control limit width W Moving average width k e slope of the lower arm D Distance typically set to 10 λ e weighted parameter, 0 < λ ≤ 1 UCL Upper control limit LCL Lower control limit Time Figure 1: e SPC chart scheme of motorized spindles.

Mathematical Problems in Engineering
Given the WPHM, then the conditional reliability can be obtained as where the data of failure time are denoted by T and Z(s) is the actual value of Z(t) at time s. e corresponding probability density function (PDF) is derived as Based on the condition monitoring system, the degradation features are acquired in real time, and the covariates could be predicted and updated, expressed as erefore, the remaining useful life of motorized spindles can be predicted as derived in [24] by According to the monitoring data of multiple features and the failure event data, the parameters of the WPHM can be estimated by the maximum likelihood method. Since the motorized spindle may run to failure or be maintained before failures in the actual CBM process, the failure event data collected involve both time to failure and censored time. Taking these data types into consideration, the likelihood function is given by where t i represents the value of time to failure, t s represents the value of censored time, n is the number of failure times, and m is the number of censored time. e parameter estimate of β, η, and c can be calculated by maximizing equation (7). As a result, the hazard rate in real-time operation can be extracted by the PHM; it is adopted as the health indicator interpreting the current performance state. In the next step, SPC techniques will be applied to monitor the changing trend of the hazard rate and provide CBM decision support.

Integrated Model and Optimization.
In the case of combining SPC charts with CBM scheduling for maintenance strategy optimization of motorized spindles, four possible scenarios of the integrated CBM model can occur during a renewal cycle, as shown in Figure 2. A renewal cycle is identified as the period between start of process monitoring and a maintenance activity [25], and the four different scenarios are further illustrated as follows: Scenario 1 (S 1 ): the performance of the motorized spindle remains in control with no false out-of-control alarm, which is interpreted as h(t, Z(t)) ≤ UCL When a hard failure occurs in the interval (kΔ, (k + 1)Δ), k � 0, 1, · · · , n, a corrective replacement is immediately implemented and restores the spindle to its initial condition. Scenario 2 (S 2 ): the performance of the motorized spindle remains in control until a false out-of-control alarm is given at (k + 1)Δ, k � 0, 1, · · · , n. A preventive replacement is performed and restores the spindle to its initial condition. Scenario 3 (S 3 ): the performance of the motorized spindle shifts to the out-of-control condition in the interval (k 1 Δ, (k 1 + 1)Δ), k 1 � 0, 1, · · · , n, which is interpreted as h(t, Z(t)) ≤ UCL. However, the SPC chart fails to capture an alert signal until the system results in a hard failure in the interval (k 2 Δ, (k 2 + 1)Δ), k 2 � 0, 1, · · · , n. A corrective replacement is then executed and restores the spindle to its initial condition. Scenario 4 (S 4 ): the performance of the motorized spindle degrades to the out-of-control condition in the interval (k 1 Δ, (k 1 + 1)Δ), k 1 � 0, 1, · · · , n − 1, and the SPC chart succeeds to detect the out-of-control signal at k 2 Δ, k 2 � k 1 + 1, · · · , n. A preventive replacement is carried out and restores the spindle to its initial condition.
Considering the four possible scenarios, the average long-run CBM cost c can be transformed to the expected cost over the expected length of a renewal cycle [26], which is given by where . p(S i ) represents the probability of the ith scenario, p(S ij ) represents the probability of the jth state in the ith scenario, and m i represents the number of states of the ith scenario, which are expressed as  e optimization objective is to obtain the minimum value of the average long-run CBM cost by setting up an optimal control limit UCL * as follows: Taking the situation when n � 1for illustration, the CBM elements of each scenario are presented in Tables 2-4.   4 Mathematical Problems in Engineering

Scenario Symbol
Probability of occurrence e threshold of the degradation feature, which has a relationship with UCL as D � 1/cIn(UCL/h 0 (t)) P d e probability of the signal when the process is out of control P u e probability of no signal when the process is out of control P f e probability of the signal when the process is in control P n e probability of no signal when the process is in control

The Three-Stage Approach for CBM Optimization
In this section, the proposed approach developed to monitor the health degradation and optimize the CBM schedule for motorized spindles is described in Figure 3. e approach is composed of three stages involving health indicator extraction, degradation condition monitoring, and CBM optimization. In order to implement this approach, the PHM is constructed, and the hazard rate is estimated, which is then applied in SPC charts as the monitoring statistic. At last, combining SPC with CBM management, an integrated optimization model is given to develop the maintenance plan and minimize the average long-run cost of maintenance by searching for the optimal control limits. e three stages of the proposed approach are outlined as follows: Stage 1 (construction of the PHM): the model is adopted to extract the estimation of the Stage 2 (establishment of the SPC charts): as long as the parameters of the PHM are estimated with the training data, the test data can be collected and plugged into the model to obtain the hazard rate as the degradation feature to be monitored. e samples of the hazard rate are then plotted on the control chart with the corresponding central line (CL) and UCL, which can be computed as Table 1 according to the selected chart type.Stage 2 (establishment of the SPC charts): as long as the parameters of the PHM are estimated with the training data, the test data can be collected and plugged into the model to obtain the hazard rate as the degradation feature to be monitored. e samples of the hazard rate are then plotted on the control chart with the corresponding central line (CL) and UCL, which can be computed as Table 1 according to the selected chart type.hazard rate by the fusion of event data and condition monitoring data. is stage is illustrated in detail as follows: Step 1 (data acquisition). e performance of motorized spindles mainly depends on their key functional components, such as the bearings, the tool clamping system, and the rotary unions [26]. erefore, these components would be focused by installing additional sensors, and different signals are collected comprising vibration, current, force, acoustic emission, temperature, leakage flow, etc. [19], in order to describe the degradation characteristics of the system more accurately.
Step 2 (data preprocessing and feature extraction). Using the collected life cycle signals of the monitored spindle to select the typical features in time, frequency, and time-frequency domains, approaches such as principal component analysis can be adopted for the purpose of reducing the dimension of features while keeping the information with the largest contribution to describe system degradation process.
Step 3 (health indicator estimation). Combining the extracted features from the condition monitoring system with failure event data in the lifetime of motorized spindles, the PHM is constructed, and its parameters are estimated with the training samples. erefore, the hazard rate observation can be obtained as the health indicator for the next stage of degradation condition monitoring.
Stage 2 (establishment of the SPC charts): as long as the parameters of the PHM are estimated with the training data, the test data can be collected and plugged into the model to obtain the hazard rate as the degradation feature to be monitored. e samples of the hazard rate are then plotted on the control chart with the corresponding central line (CL) and UCL, which can be computed as Table 1 according to the selected chart type.
e chart is then adopted as a tool for anomaly detection and a trigger for maintenance activities. In particular, if a sample point plotted on the SPC chart is within the UCL, the system is considered as in-control, and the periodic inspection and daily health management are maintained. Otherwise, the system is defined as out of control, which implies that the degradation level exceeds an acceptable threshold, and a preventive replacement is performed. Meanwhile, a failure of the system could occur during both in-control and out-ofcontrol phases, which will trigger a corrective replacement. Stage 3 (optimization of the CBM): based on the aforementioned maintenance strategy, the monitoring system integrating SPC charts and PHM is further developed for CBM optimization. Four different scenarios that can describe the possibilities in monitoring process are taken into account in order to establish the CBM model. Each scenario is ended by a maintenance activity, and the monitored system is restarted as a new cycle with an "as good as new" condition. Finally, the SPC chart is updated, and the optimal CBM schedule is obtained through minimizing the average long-run cost during a renewal cycle by optimizing the UCL.

Case Study
e proposed three-stage approach is applied to a real case of the motorized spindle in order to illustrate its strong practicability. e developing process of degradation and failure of motorized spindles are mainly blamed on bearing failures, clamping device rubbing, and rotary union loosening, which can be characterized by vibration signals [27]. erefore, the case reported here is based on vibration signal analysis in particular.

Experimental Setup.
e failure events and condition monitoring data sets are generated from the motorized spindle going through run-to-failure experiments under constant rotary speed at 4000 RPM. e experiment rig is shown in Figure 4. e experimental setup is designed considering the structure of the motorized spindle, which is shown in Figure 4(a); the spindle can be subjected to the cutting torque provided by the dynamometer joined to the tool holder with the couplings and the biaxial cutting forces provided by the axial and radial loading systems. e spindle is kept running 12 hours per day except for irregular equipment adjustment and unexpected failures.
Based on experimental requirements, two PCB 356A15 triaxial high sensitivity, ceramic shear ICP accelerometers are installed on the front face and the back face of the motorized spindle, respectively. e specific placement is circled in Figures 4(b) and 4(c), and the vibration signal is collected every 4 hours by ADLINK's USB-24054-CH, 24bit, 128 kS/s dynamic signal acquisition. e test and the condition monitoring system are controlled by LabVIEW program.

Feature Extraction and PHM Construction.
e training data collected include time of inspections, time to failure, censored time, and vibration signal. Each inspection lasts at least 30 seconds, which contains 150,000 vibration signals with the sampling rate of 5 kHz. A total of 246 inspections spanning over 1000-hour operation test on the tested motorized spindle are conducted. During the whole experimental period, 7 hard failures and 5 suspensions are observed. Among them, the 7 failures are considered as unexpected functional failures, and the 5 suspensions are due to adjustments of experimental equipment.
Only the x-axial vibration measurements were used for analysis, and three commonly used time-domain features are extracted from the lifetime vibration data of the tested spindle. e technicians believe that the selected features can contain essential information to interpret the specific failure well [28]. e features as derived in [29] are summarized in Table 5.
PHMs with different features are fitted, and K-S goodness-of-fit tests are performed to obtain a better model. e parameters and the test statistics of hazard rate functions are estimated and listed in Table 6. e traditional Weibull distribution without incorporating any feature is also fitted, and its result is shown in Table 6 as Model 1.
It implies in Table 6 that the significance of RMS in the PHM is greater than kurtosis and crest factor. In addition, the fitting effect of the PHM integrating with condition

SPC Establishment and Degradation
Monitoring. For a monitored renewal cycle of the motorized spindle, a life cycle of 372 hours with 93 inspections in total is adopted to obtain the 93 testing values of RMS as the lifetime features, which are shown in Figure 5(a). en, the set of RMS is plugged into the PHM obtained from the training data set in Section 5.2, and the hazard rate is calculated and plotted in Figure 5(b). e hazard rate set of the testing cycle is then used as monitored health indicators. e SPC charts introduced in Section 2 are constructed as supports for anomaly detection and CBM optimization.
e Shewhart I-chart and the EWMA chart are taken as examples, which are shown in Figures 6(a) and 6(b). In the case of monitoring the hazard rate, only UCL and CL are essential. e red solid lines represent the UCLs of the charts, the green-dotted lines represent the CLs of the charts, the black dots are corresponding statistics derived from the testing data of the hazard rate, and the marked purple circles are degradation states determined to be out of control.
To be more specific, the two charts are established following the rules in Table 1. As for the Shewhart I-chart in Figure 6(a), the estimated process mean μ 0 � 0.0219, the estimated process standard deviation σ 0 � 0.0013, UCL � 0.0258, CL � 0.0219, and the control limit width n � 3, and the process is determined to be out of control at the 62nd inspection. Figure 6(b) depicts the results of the EWMA chart, whose weight parameter is set to be λ �0.25, the UCL is stable at UCL � 0.0234, and the process is determined to be out of control at the 60th inspection.
Since the Shewhart I-chart monitors the performance of the motorized spindle only based on observed data at individual inspections, it is less sensitive in detecting degradation with small and moderate change trends in comparison with the EWMA chart [30]. However, in the case of monitoring the hazard rate of the tested spindle, the results of the EWMA chart show little significant deviation comparing to the Shewhart I-chart, which is mostly the consequence of using the hazard rate from the proposed PHM as the monitoring statistic. e model could compromise multiple information based on time series in order to eliminate autocorrelation in the process and reduce the influence of the random error. As a result, the two types of charts' abilities for gradual degradation detection are satisfactory.

CBM Schedule Optimization.
To build the final CBM model, the established SPC chart is keen to be updated and optimized by finding optimal UCL * , which is used to trigger corresponding maintenance actions so that the inspection and maintenance cost can be minimized.
In the current case, the tested spindle is inspected every △ � 4 hours, the cost of each inspection C i is ￥30. e preventive replacement cost C p is ￥500, and the corrective replacement cost C c is ￥6000. To simplify the analysis, P d and P f are set to be 0, and P u and P n are set to be 1, which reduce the four scenarios to two for the control chart. Finally, the CBM model can be interpreted as follows: Scenario 1 (S 1 ): the motorized spindle remains in control until a hard failure occurs in the interval (4k, 4(k + 1)). A corrective replacement is implemented and restores the spindle to its initial condition. Scenario 2 (S 2 ): the motorized spindle degrades to the out-of-control condition in the interval (4k, 4(k + 1)), and the SPC chart succeeds to detect it at 4(k + 1). A preventive replacement is carried out and restores the spindle to its initial condition.
As a result, the long-run cost rate is obtained through the following formula: Based on the vibration signals at each inspection of the degradation process, RMS is found to follow the normal distribution with parameters μ(t)and σ(t). e changes of the two parameters with time are studied with respect to the 93 inspections, which are expressed as μ(t) � 0.11 exp 1.14 × 10 − 3 t + 8.96 × 10 −16 exp(0.09t).
σ(t) � 0.04 exp 2.65 × 10 − 4 t + 1.87 × 10 −8 exp(0.04t). (12) erefore, the probability that the hazard rate exceeds UCL at time t is derived as rough optimization of equation (11), the optimal CBM schedule is finally obtained with UCL * � 0.0524, and the minimal long-run cost rate is c � 15.1942, which is shown as the red-marked circle in Figure 7(a). According to equation (3), the RMS threshold versus the testing time of the motorized spindle can be obtained and depicted in Figure 7(b). e result also has certain directive significance in practical application because the preventive replacement can be simply triggered when the RMS observation at the inspection exceeds the RMS threshold.

Conclusions
e proposed three-stage approach in this paper is committed to fusing various degradation features, monitoring the degradation process, and developing the optimal CBM plans for motorized spindles. e combined effect of condition features degrading over time together with failure events is interpreted by the PHM, and SPC charts are then introduced to monitor the changing trend of the hazard rate obtained from the PHM. Application of the SPC charts with periodic hazard rate inspections helps in CBM modelling, and four possible scenarios are defined to constitute the maintenance model. At last, the CBM plan with minimum long-run cost is obtained by seeking for the optimal UCL. e feasibility and practicability of the approach are proved by a case study using the actual spindle vibration data. e analysis results indicate that the PHM contributes to describing the degradation states of the motorized spindle more accurately, and three different time-domain features are extracted from the raw vibration signals and are used as covariates in the PHM; RMS is verified to have the best fitting effect of the experimental data. e proposed approach allows not only maintenance optimization but also accurate characterization of health degradation. It can help maintenance technicians and managers reduce breakdown time, delete unnecessary maintenance actions, and minimize maintenance costs in an intuitive and efficient way. e approach in this paper can also be applied for the motorized spindle operating in CNC machine tools for intelligent manufacturing in order to achieve better performance and a longer life of the whole manufacturing supply chain.

Data Availability
No data were used to support this study.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.