Eccentricity Parameters Identification for a Motorized Spindle System Based on Improved Maximum Likelihood Method

As a kind of rotor system, the electric spindle system is the core component of the precision grinding machine. +e vibration caused by the mass imbalance is the main factor that causes the vibration of the grinding machine. Identifying the eccentricity parameters in an electric spindle system is a key issue in eliminating mass imbalances. It is difficult for engineers to understand the approximate range of eccentricity by experience; that is, it is difficult to obtain a priori information about eccentricity. At the same time, due to the geometric characteristics of the electrospindle system, the material factors and the randomness of the measurement response, these uncertain factors, even in a small case, are likely to cause large deviations in the eccentricity recognition results. +e search algorithm used in the maximum likelihood method to identify the eccentricity parameters of the electrospindle system is computationally intensive, and the sensitivity in the iterative process brings some numerical problems. +is paper introduces an Advance-Retreat Method (ARM) of the search interval to the maximum likelihood method, the unknown parameter increment obtained by the maximum likelihood method is used as the step size in the iteration, and the Advance-Retreat Method of the search interval is used to adjust the next design point so that the objective function value is gradually decreasing. +e recognition results under the three kinds of measurement errors show that the improved maximum likelihood method improves the recognition effect of the maximum likelihood method and can reduce the influence of uncertainty factors on the recognition results, and the robustness is satisfactory.


Introduction
e high-speed motorized spindle is supported by bearings, which integrates the functions of the machine tool spindle and the high-speed motor in structure, and realizes the "zero transmission" between the variable frequency motor and the machine tool spindle. It is a product that combines the motor with the spindle. e rotor of the motor is the rotating part of the spindle. In theory, the motorized spindle can be regarded as a high-speed motor. It focuses on many key technologies, such as the design, manufacture, assembly and control of high-speed motor, and high-performance machine tool spindle, and has the characteristics of high-speed and zero-transmission chain. e motorized spindle is the core component of optical grinding machine tools, and its performance has a significant effect on the machining accuracy because of the combination of spindle and tool. In recent years, with the increase of grinding machine speed, machine tool failures frequently occur, and the research on safe operation and fault diagnosis of high-speed grinding machines has become an important research content in the development of machine tool technology [1,2]. Unbalance is one of the most common faults in high-speed grinding machine, which causes equipment vibration, coupling deflection, bearing wear and tear, and other hazards [3,4]. e fault will cause the electric spindle system to bear the eccentricity parameter, which will cause changes in the dynamic performance of the spindle system, such as eigenvalues, thus aggravating the operation fault. Motorized spindles widely used for high-speed precision machine tools are very sensitive to the mass unbalance of rotors; thus, their balancing problem is always a research hotspot. Although many significant studies were done regarding the theory and application of various rotor balancing technologies for motorized spindles. Dynamic balance technology including influence coefficient method, cross-correlation method, and mode method [5,6] is commonly used to balance machine tool spindle faults at present. e basic process involves installing the acceleration and speed sensors, collecting the spindle vibration signal, using a suitable data processing method to the extract amplitude and phase of the signal, and outputting the correction masses and mounting positions. However, the particularity of motorized spindles is not considered carefully in the existing balancing approaches. When the rotor unbalance of a motorized spindle occurs in operation, it is subject to both the mass unbalance-induced inertia force and air gap unbalance-induced electromagnetic force, which is an important feature that distinguishes the motorized spindle from a mechanical spindle [7]. In addition, dynamic balance technology of machine tool depends on whether the unbalance parameters of the machine tool are obtained accurately or not. e eccentricity of grinding wheel and motor is the main factor of unbalanced fault in high-speed motorized spindle system. e eccentricity parameters should be accurately obtained before using dynamic balance technology. Using parameter identification technology to identify the eccentricity parameter of motorized spindle system is a necessary work in diagnosing unbalance fault of rotor system. However, for the electric spindle system, in engineering practice, the eccentricity parameters are mainly caused by the disk parts of the motor and the grinding wheel. It is affected by many factors, and it is difficult for the engineer to understand the approximate range of the eccentricity parameters based on experience; that is, it is difficult to obtain prior information of the eccentricity parameter. At the same time, due to the randomness of rotor geometry, material factors, and measurement response, these uncertainties, even in a relatively small case, may also lead to large deviations in the identification results of eccentricity parameter. How to reduce the influence of uncertainties such as randomness and identify quickly eccentricity parameter is a core of researches on engineering inverse problem under uncertainty.
On engineering inverse problem under uncertainty, when the sample information of uncertain parameters is sufficient, the probability density can be used to describe the uncertainty of parameters. When the prior information of the parameters identified is unknown, the maximum likelihood method is often used to identify structural parameters or load parameters. In the maximum likelihood method, the parameter value when the maximum value is taken by the likelihood function is used as the identification result, and the corresponding confidence interval is calculated [8][9][10]. Liu et al. [11] proposed the maximum likelihood method based on the sensitivity matrix method (SMM), which is suitable for dealing with the uncertainty inverse problems with insufficiency and imprecision in the input and/or output parameters.
is proposed method can transform some known parameters into explicit solvable form and identify unknown parameters by iteration method. However, due to the unknown prior information of the parameters identified, that is, the interval of the parameters identified is unknown, and the sensitivity will bring some numerical problems in the iterative process, which affects the accuracy and reliability of the parameter identification result, and sometimes causes the identification process to diverge and the effective identification result is not obtained [12]. In this paper, Advance-Retreat Method (ARM) of the search interval is introduced to the SMM, the strict requirements of the traditional iterative numerical method for the search space are avoided, and the SMM is improved to make the iteration suitable for complex engineering optimization problems. It is possible to improve the identification effect of the SMM.

e Basic Idea of SMM to Identify Eccentricity Parameters.
e motion equation (1) of the motorized spindle system with unbalance is expressed as follows [13]: where M, C, K, and G are the mass, damping, stiffness, and gyroscope matrices which can be obtained by the method used in [13]; a, _ a, and € a are the unbalanced displacement, velocity, and acceleration vectors of the motorized spindle system, respectively; and F e is an unbalanced force.
In the paper, the unbalanced force of the motorized spindle system is caused by the eccentricity parameters of the motor and grinding wheel. e eccentricity parameters are unknown; the unbalanced response measured by experiment is known. e identification problem is to determine the unknown eccentricity parameter from the observed unbalanced response. However, the input part parameters, such as motor length and motor radius, and output parameters such as unbalanced response are all partly known; that refers to the measurement errors. e identification problem belongs to the uncertainty inverse problem with both insufficiency and randomness in the input and output of structural systems. e unbalanced response of equation (1) is expressed as a forward solver as follows: where X is vector of input parameters, Y is vector of output parameters, and T is the forward solver representing the translation process from input to output. e subscript u denotes unknown and the subscript k denotes known. e input vector X consists of two parts; one part is the known parameter vector X k , which obeys the probability density distribution, such as the length and radius of the motor. e other part is the unknown parameter vector X u , such as the eccentricity parameters of the motorized spindle system. e output vector Y also contains two parts: the unknown Y u and the known Y k . Y k is a known vector derived from experimental observations and subject to a certain probability density distribution, such as the unbalanced displacement. e maximum likelihood method is a reverse method for finding the mean and confidence interval of an unknown parameter X u , Y u by the probability density distribution of the known parameters X k , Y k in the input and/or output of the system. Introducing the sensitive matrix S [14], such as equation (3), the incomplete information of input and output known parameters in equation (2) can be transformed into explicit solvable form of equation (4): is the change in input parameters, and m and n are the number of output and input parameters, respectively.
where D � According to the maximum likelihood principle, equation (4) is transformed into a similar least squares form, as in where σ 2 1 and σ 2 2 are variances, A w � WA, D w � W D,W � diag (1/σ 1 , t1/σ 2 ), σ xi k and σ yi k are deviations of the ith parameter in the input and output known vector e iteration method is used to solve the increment B of unknown parameters. According to probability theory, the least square of equation (5) is the maximum likelihood solution: Considering the nonlinearity of the structural system, the confidence interval is usually calculated by the Markov Chain Monte Carlo method (MCMC). 95% confidence interval of unknown parameters [15] is where cov (B) is the covariance of unknown parameters, H � B T − (B c ) T is an N rows and m columns matrix, representing the distance between the calculated value of parameter B and the average value B c . N is the number of calculations.

ARM for Determining Search
Interval. e basic idea of sensitive matrix method and maximum likelihood method (SMM) is to transform the inverse problem of stochastic uncertainty into an explicit iterative optimization problem similar to the least squares form. e unknown parameters are identified by minimizing the deviation between the experimental and calculated values. e mean and confidence intervals of unknown parameters are gained by using input and output known parameters with normal distribution characteristics. at is, the unknown eccentricity parameters are identified by solving the optimization problem that minimizes the error between the unbalanced displacement calculated and measured at the corresponding measuring point by adjusting repeatedly the eccentricity parameters of the motorized spindle system, as in (8) where X i u and X i+1 u are the input unknowns parameters corresponding to current step i and next step i + 1, respectively. Y i+1 k and Y i k are the unbalanced displacement numerically and accurately calculated based on equation (2), Y s k is the output known parameters obtained via measurement. ε 1 and ε 2 are allowable errors of iteration stopping.
Optimization methods such as Gradient Descendent [16] require step size and search direction as iterative guidance in the search interval. e SMM can obtain the increment B of unknown parameters in each iteration process, but the prior information of unknown parameters is not considered, so the search interval in the iteration process is unknown. It will lead to blindness of search and affect the accuracy and reliability of parameter identification results. In this paper, the maximum likelihood function is improved by introducing the ARM for determining search interval [17,18], the trend of "high -low-high" of the value of the objective function formed by the design points is explored, and the value of the objective function is towed toward the direction of gradual descent. Equation (8) shows that the SMM should continuously reduce the error between the calculated results and the measured results in the iteration process.
is results in the comparison of the objective function values of the next step with those of the previous step, as well as the comparison of the calculated objective function values with the test results. In this way, the value of the objective function will be pulled down gradually during the iteration process, and the error between the calculated value of the objective function and the test results will be reduced or increased. Considering the unknown search interval and the problems arising from the iteration process, the ARM is introduced to improve equation (8).
e unknown parameter increment B obtained by the SMM is used as the step size in the iteration process and the ARM of the search interval to adjust the next design point, so that the objective function value is towed in the direction of gradual descent. When the error between the objective function value and the test result becomes larger, the design point of the previous step is selected as the initial design point and the search is restarted. e selection of design points is shown infd9 Shock and Vibration

Solution Procedure
To sum up the above arguments, the flow chart of the present method is shown in Figure 1. e main steps are as follows: Step 1. Set the error constant ε and the probability density distribution of known parameters, such as normal distribution , to test and obtain the unbalanced displacement Y s k at the sensitive points of sensors that are easy to arrange Step 2. Assume the unknown parameter value X i u , combining them together with the known parameters X k into the forward solver to calculate the corresponding output 〈Y i u , Y i k 〉 Step 3. Calculate the sensitivity matrix S i centering on the present parameters 〈X i u , X k , Y i u , Y i k 〉 based on equation (3), and the increment B of unknown parameters in equation (6) is obtained by SMM. Update Step 4. If all these errors are within the predefined toleration, this set of parameters X i+1 u is considered to be the total solution of the problem, and the solution procedure ends. Otherwise, the unknown parameter X i+1 u is updated by the ARM equation (9), and return to Step 2. Namely, Step 5. Use equation (7) to calculate the confidence interval of unknown parameters.

Description of the Numerical Example.
e motorized spindle system that referred to an inverse problem of identification for the eccentricity parameters is given in the section to validate the proposed method. e model parameters of the motorized spindle system are shown in Figure 2. e diameters of the grinding wheel, the motor, and shaft are 1.4 m, 2.4 m, and 0.6 m, respectively, their lengths are 0.04 m, 0.5 m, and 1.3 m, respectively, and the front and rear bearings are installed at a distance of 0.5 m and 0.2 m from both ends of the shaft. e performance parameters of the motorized spindle system are listed in Table 1.
e average air gap length is 8 mm without eccentric, the air permeability coefficient is 1.256 * 10 − 6 , the air gap fundamental wave magnetomotive force coefficient is 5.2, and the rotating frequency is 100 Hz.
According to the solution procedure described in Figure 2, the unknown input parameters are eccentricity parameters (e 1 , e 2 ); the known input parameters are motor length and diameter (L, D). e output parameters are the radial displacement (d 1 , d 2 , d 3 , d 4 ).
If the radial displacement of this measured point is very sensitive to eccentric parameters, but it is not convenient to arrange sensors at this point, the radial displacement of this measured point is attributed to unknown parameters. Considering the influence of random errors in measurement results, the known parameters in input and output are random variables subject to a normal distribution. Taking 3% measurement deviation as an example, the known parameters are expressed as L ∝ N(0.5, 0.015 2 ) and D ∝ N(2.4, 0.072 2 ).

Sensitivity Analysis.
e optimization method implementation criterion requires that the unbalanced response of the measuring point is sensitive to the eccentricity parameter. It is necessary to evaluate the influence of the eccentricity parameter on the response data. is is accomplished using a sensitivity analysis combined with forward solver  based on Finite Element Simulation Analysis to spindle system. e results of the sensitivity analysis in the x-direction and the y-direction of each point on the spindle are shown in Figure 3. It can be seen that the eccentricity is sensitive to the unbalanced response. It is feasible to identify eccentric parameters based on optimization method through minimizing the error squared of the unbalance response between the experiment results and the computational ones. e experiment unbalance response is measured by displacement sensors. Due to the electromagnetic effects of eddy current sensor, it should be installed in a suitable position. In this paper, there are two unknown parameters and two known parameters in the input parameters. In order to avoid the morbidity of the sensitive matrix in formula (3), the number of known parameters and unknown parameters in the output parameters should be the same. Figure 3 shows that the unbalanced response in x-direction at both ends of the spindle is large, and the sensor can be easily arranged, so it can be used as a measurable point. e y-direction unbalanced response of the first and second bearings is of great concern. It is not easy to arrange sensors and treat them as unmeasurable points. Table 2 gives the corresponding solution for each iteration. According to the identification process, the parameter vector is written as follows:  Table 2. en, an initial increment B 0 is obtained from this set of initial parameters    Table 2). And subsequently, substituting the originals of eccentricity parameters with the newly computed values of 0.0639 and 0.1676, and combining them with the known X k into the forward solver FEM once again, a new set of radial displacement is calculated to be 0.1055, 0.035, 0.0366, 0.0459 { } μm, respectively, as shown in the second row of Table 2. en, the error defined in equation (8) is examined to decide if the solution procedure ends. When the stopping criterion in equation (8) is not satisfied, the unknown parameters are updated by the ARM from equation (9), the updated eccentricity parameters e 1 1 , e 1 2 are 0.1916 μm and 0.5027 μm, and the solving process continues to the next iteration. After 981 times of such iterations, the maximal error of the sought parameters with   Table 2 gives the corresponding solutions at each of the iterations. Figure 4 shows that eccentricity parameter identification iterative process converges; the displacement error is decreasing at every iteration, which makes the eccentricity parameter to be identified approach the most advantageous. e eccentric parameters are substituted into improved maximum likelihood method from zero. With the increase of iterative steps, the simulated unbalanced response at the grinding wheel and shaft end approaches the experimental displacement measured by the displacement sensor step by step, and the displacement error decreases, about 200 steps. e change of eccentric parameters is small. e given convergence criterion is fulfilled after 981 iterations of training. e iterative evolution process of parameter identification in Figure 4 shows that the improvement measures made in this paper to the maximum likelihood method are effective.

Identification and Discussion for Different Measurement
Deviations.
e influence of the three measured deviations of 1%, 3%, and 5% of the known parameters on the identification results of eccentricity parameters is considered at the motor and grinding wheel. e identification results are    given in Table 3. Comparing to the different measurement deviations, the improved maximum likelihood method can reduce the influence of uncertainty factors on the recognition results, and the robustness is satisfactory. In each iteration step, equation (7) is used to calculate 95% confidence interval. Comparisons of the confidence interval of the eccentricity parameters calculated in each iteration step are shown in Figure 5. Figure 5 shows the 95% confidence interval iterative process for parameters under three errors. e eccentric parameters change a little at about 500 steps and float near the mean. It can be seen that the confidence interval increases with the increase of the measurement error.

Conclusions
e unbalanced fault of the motorized spindle system of the grinder is caused by the coupling of several mass eccentricity parameters, so it is difficult for engineers to understand the approximate range of mass eccentricity parameters based on experience. Due to the existence of uncertain parameters such as the structure parameters and test conditions of the grinder motorized spindle system and the unknown prior information of the mass eccentricity parameters, the search calculation used in the maximum likelihood method to identify the mass eccentricity parameters of rotor unbalanced faults is faced with a large amount of calculation, and the sensitivity in the iterative process will bring some numerical problems.
e Advance-Retreat Method of the search interval is introduced into the sensitive matrix method and the maximum likelihood method to improve the maximum likelihood method. e iteration in the improved maximum likelihood method is suitable for dealing with complex engineering optimization problems and avoids the stringent requirements of the traditional iterative numerical method on the search space. In this paper, the sensitive matrix method and maximum likelihood method are developed based on Advance-Retreat Method to identify the eccentricity parameters in a motorized spindle system with input size and output test response with random measurement error. e feasibility and validity of developed algorithm has been tested with numerical example. ree measurement deviations are considered in the numerical application. e comparison of parameter identification means and 95% confidence interval indicates that it is robust to the insufficient observations.
Data Availability e authors confirm that the data supporting the conclusion of the article are shown in the relevant figures and tables in the article.

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