Research on High-Precision Form Grinding Technology of Gear Based on Ambient Temperature Adaptability

The thermal error of high-precision computer numerical control (CNC) form grinding machine has become the critical factor affecting its locating precision. Because the thermal error is more complicated to be measured directly, most of the measurement tests are usually aimed at no-load conditions. In this article, a method to evaluate the machine tool thermal error based on the principle of tooth profile error through the gear grinding precision is presented. Based on high-precision CNC form grinding machine need to maintain isoperibol, this paper proposed the ambient temperature as a variable, for different ambient temperature to measure the thermal error of the machine tool. According to the measurement results of the thermal error, the method of partial least squares neural network is used to structure the thermal error prediction model of the machine tool. The experimental data show that the grade of the gear precision can reach four levels at different temperatures. This approach shows a promising prospect in the application of high-precision CNC gear form grinding machine in the future.


Introduction
New research shows that the thermal error of CNC accounts for 70 percent of total error.Moreover, there are many means to reduce the thermal error of the CNC, for instance, to design the thermal symmetry structural of the machine tool and to separate the heat source from the main machine body.Nevertheless, it must result in much material cost and higher manufacture cost, due to change the structure design.Error compensation is an economical and effective compensation method [1][2][3].
The vibration and noise of the gear transmission system are one of the hotspots concerned by the researchers of the gear transmission system.Due to the competition in the automotive market and many countries' requests to the vehicle noise has been more and more strict.German Volkswagen put forward some requests accordingly on one tooth and three cross sections for gear form grinding process in 2007, as shown in Figure 1.There are few manufacturers of gear grinding machines that can satisfy some performance specification requirements.With the popularity of new energy vehicles, it is necessary to set higher requirements for the accuracy of gear manufacturing.It will lead to abnormal gear meshing, due to gear manufacturing error, installation error, and inadequate support system stiffness, gear modification becomes an important method to solve this problem.In order to obtain the best transmission performance, Wang et al. discussed the form grinding method orient gear meshing performance, while the high-performance gear manufacture cannot separate from the high-precision gear form grinding technology [4].At the same time, in order to get highperformance gear, grinding wheel accuracy is very important.Ding et al. made further research in the related fields [5][6][7][8][9].
In order to obtain high-precision gear, universities and other institutions of scientific research institutions spare no effort in relate fields and achieved fruitful results [10,11].The model of thermal error in NC machine tool is usually divided into empirical thermal error model and theoretical thermal error model.The empirical thermal error model is Mathematical Problems in Engineering based on the mathematical model of the temperature rise of the sensitive point, the machine tool speed, and the feed rate, and the thermal drift is the output.The establishment of the theoretical thermal error model requires the three forms of heat transfer as the basis for calculation.After solving the temperature field of a certain part or machine of the machine tool, the corresponding displacement field is calculated.Because of the large space of factory, it will require tremendous amounts of energy, if the whole space keeps at a constant temperature.Therefore, it is necessary to improve the precision grinding technology to a new level through technological innovation.
In the thermal error compensation of machine tool, in order to improve the efficiency, some scholars consider that combining the measured error with the comprehensive model [12][13][14][15][16][17].For example, Zhang et al. have presented a method of obtaining 21 error elements in space by 22 straight lines in 1988 [12] and Chen et al. on the base of the previous study that decrease from 22 lines to 15 lines [13].Domestic researchers also proposed 14-line and 9-line measurement method.Efficiency has also increased with the line gradually reducing, but the actual measurement is more complicated and more time-consuming.The less measurement lines mean less workload of measurement but cannot fully reflect the overall thermal deformation of the measurement system, and sometimes through a variety of measurements, it may not be able to find the best measurement points.So it cannot be good for measuring the thermal error of machine tools to provide better measurement ideas.
At the present stage, in order to make it easier to measure the thermal error, usually use no-load way to preheat the machine to obtain the correlation data about thermal error.In the actual production process, the machine tool cannot always be in no-load condition.The no-load and normal condition belongs to two vastly different types of states.By the no-load condition, we arrive at some conclusions which are insufficient guiding significance to production practice.Based on the above shortcomings, according to the multiple of thermal errors, taking gear form grinding process for research background, based on the formation principle of spiral bevel gears surface errors, this manuscript studied the mechanism of thermal error compensation about the CNC machine tool under nonisothermal conditions.This paper presents a method to measure the CNC thermal error indirectly, and then through the modeling and prediction of related parameters, we get the method of the error compensation which can meet the actual processing requirements.

Error Analysis Hypothesis
In this paper, we take the high-precision NC forming grinding machine developed independently as study object and centers on its related problem.The machining region is shown in Figure 2. In order to grasp the principle contradictions about thermal error for NC machine tool, some assumptions are made as follows.
ree Sections ree Sections  First, due to the continuity of the production process, it will not study the compensation before the machine reaches thermal equilibrium.
Second, because the machine tool is in a steady state of operation, small machining allowance, and small mutual interaction between grinding wheel and gear, so we take no account of the dynamic change of grinding force and other factors.
Third, the wear of grinding wheel throughout the production process is not considered.
Fourth, the products produced by the factory in a period of time are relatively simple and the process is relatively fixed.
Fifth, the deformation of the parts of the machine tool has a single-directional.
Sixth, after dressing grinding wheel, the wheel profile has met the requirements.

Error Analysis of Moving Parts in
High-Precision Computer Numerical Control Form Grinding Machine   machine is working; the error can be divided into translation error and angular error [18].
If the part moves along the X-coordinate axis, the six degree errors of the moving part are as follows: the straightness error of the X direction is   (, ), the straightness error of the Y direction is   (, ), the straightness error of the Z direction is   (, ), the rolling error around X direction is   (, ), the rolling error around Y direction is   (, ), and the rolling error around Z direction is   (, ), as shown in Figure 3.There are 2 main prismatic pairs in the highprecision computer numerical control forming gear grinding machine.The major form of error is similar to this, so no longer elaborated.The error of X direction and Y direction is shown in Table 1.

Analysis of Motion Error of Revolute Pair.
Rotating parts of the high-precision computer numerical control forming grinding machine also have six degree errors.Assuming that the origin position and the coordinate axis of the motion coordinate system A and the relative reference coordinate system R are exactly the same before starting the motion, when the rotation angle around the A axis's Z direction is 6 degrees of freedom error, rotation angle error is   (, ), rotation angle error is   (, ), rotation angle error is   (, ), translation error is   (, ), translation error is   (, ), and translation error is   (, ), as shown in Figure 4.There are 2 main revolute pairs in the highprecision computer numerical control forming gear grinding machine.

Analysis of Vertical
Error.In the movement process, besides the motion errors and revolute pair errors, the high accuracy computer numerical control forming gear grinding machine will produce perpendicularity error because the axes are not perpendicular to each other.
Because of the high-precision computer numerical control forming grinder only has two coordinate axes X and Y, only one perpendicularity is produced during the movement, that is the perpendicularity error   between the shaft X and the shaft Y.In this paper, it is converted into motion error or revolution error by projection method to be calculated.

The Comprehensive Error Model of High-Precision Computer Numerical Control Gear Grinding Machine
There are two transmission chains in the high-precision computer numerical control forming grinding machine: one is through the bed of machine tool to the gear, and the other one comes from the grinding machine to the grinding wheel.The interaction between the two transmission chains is established by the contact between the grinding wheel and the gear, and the relative position error between the wheel and the gear directly affects the accuracy of the machining.In this paper, the error synthesis model of a high-precision computer numerical control gear grinding machine is discussed by use of homogeneous coordinate transformation method.The thermal and geometric errors are compensated by this error synthesis model.

Homogeneous Coordinate Transformation.
In the threedimensional coordinate, the calculation methods of homogeneous coordinate transformation are as follows, which is shown in Figure 5. Hypothetical space has a coordinate system         , the coordinates of a point  in the coordinate system are         , when coordinate system         moves  along the   axis to get a new coordinate system         , the relationship between the coordinates of point  in the new coordinate system (  ,   ,   ) and the coordinates in the original coordinate system (  ,   ,   ) can be expressed as follows: Three-dimensional translation homogeneous coordinate transformation matrix of () expressed among   axis to translate  is Similarly, the three-dimensional homogeneous coordinate transformation matrix that translates  along   axis and  along   axis, respectively, is When the coordinate system         is rotated around the   axis   to obtain a new coordinate system         , the relationship between the coordinates of the point  in the new coordinate system (  ,   ,   ) and the coordinates in the original coordinate system (  ,   ,   ) can be expressed as where (  ) represents a three-dimensional homogeneous coordinate transformation matrix rotated   about the   axis: Similarly, three-dimensional homogeneous coordinate transformation matrix around the   axis rotation   and   axis rotation   , respectively, is When the coordinate system         , respectively, transforms , ,  along the   ,   ,   axis and then revolves   ,   ,   , respectively, around   ,   ,   to obtain a new coordinate system         , then the three-dimensional homogeneous coordinate transformation matrix of coordinate system         to the coordinate system         is Take ( 2), ( 3), ( 4), ( 6), (7), and ( 8) into ( 9): When the rotation angle   ,   ,   is very small, there is sin   ≈   , sin   ≈   , sin   ≈   , cos   ≈ 1, cos   ≈ 1, cos   ≈ 1.When the translation displacement is very small, use   ,   ,   to separate table , , and , respectively.Ignoring the second-coordinate transformation matrix is expressed as 4.2.Coordinate System Setting.In this paper, the error synthesis model of a high-precision computer numerical control gear grinding machine is discussed by use of homogeneous coordinate transformation method.The model contains geometric error and thermal error, which provides a theoretical basis for carrying out error compensation.According to the above-described high-precision computer numerical control gear grinding machine basic movement process, we establish a series of coordinate system as shown in Figure 6: machine tool coordinate system O; sliding table coordinate system X; engine base coordinate system Y; wheel frame disk coordinate system B; spindle of grinding wheel coordinate system S; grinding wheel coordinate system W; gear spindle coordinate system A; and gear coordinate system G.
For the convenience of calculations, we placed the origin of base coordinates and local coordinates in the zero point , and the same axis direction remains same.The two branches of the coordinate system drive chain are shown in Figure 7.
Branch One.Machine tool, gear spindle, and gear.
Branch Two.Machine tool, slide, machine base, wheel frame disk,, wheel spindle, and wheel.
In the ideal situation, the corresponding positions on the grinding wheel and gear are the same in space.Assume that the coordinate transformation matrix    of the gear coordinate system  to the grinding wheel coordinate system  is Based on the formation principle of different kinds of errors, the coordinate transformation matrix of the gear coordinate system  to the grinding wheel coordinate system  should be Erratic movements are indicated by superscripts .In reality, the coordinate transformation matrix of the gear coordinate system  to the grinding wheel coordinate system  is equal to the multiplication of the ideal motion matrix and an error matrix   , that is, Based on the premise of small errors, the error motion conversion matrix of the gear coordinate system  to the grinding wheel coordinate system  can be set as ,   ,   represent the relative error between the actual machining position of the grinding wheel and the ideal machining position;   ,   ,   represent the relative error between the actual machining direction of the grinding wheel and the ideal machining direction.
Based on the small error principle, after calculating the theoretical and practical transformation matrices of each sliding pair and the revolute pair, the three matrix equations ( 12)∼( 14) can be used to get the error motion matrix   :    represents the coordinate conversion matrix of the wheel spindle coordinate system  to the wheel coordinate system .The grinding wheel is installed in the spindle, moving together with the spindle.Therefore, the coordinate transformation matrix between them is unit matrix ,    = .Due to space constraints, there is seldom detailed description for the error solution procedure.Solving the matrix equation by MATLAB and neglecting the highorder infinitesimal small quantities, we finally obtain the comprehensive error model.Wheel angle deviations in three directions are   ,   ,   .Deviations in three directions are   ,   ,   .

Error Analysis of Gear Forming Grinding.
In this paper, the involute spur gear is taken as the research object, through measuring the thermal error, and then studying the thermal error compensation.In general, the error factors that affect the machining accuracy of the gear come from six directions, as shown in Figure 8.The involute spur gear is studied in this paper, so the error in Y direction has no influence on the gear machining.After dressing grinding wheel, the wheel profile has met the requirements and does not consider the deviation of revolute pair of Y direction.Because the grinding wheel is revolving and machined around the axis Z, this article does not study the above three errors; the main study is the displacement deviation along the direction X, the rotation deviation of direction X, and the displacement deviation of direction Z on the gear grinding accuracy.The above three errors, respectively, are referred to as radial error, rotating error, and vertical error.

The Impact of Machine Error on the Tooth Profile Accuracy
While using grinding wheel to grind spur gear with grinding method, the precision of tooth curve is guaranteed by the accuracy of the slide rail and fixture.With many measurements, the accuracy of the tooth curve is very high, reaching 0 or 1, as shown in Figure 9. Therefore, the main consideration of this manuscript is the influence of thermal error on the accuracy of tooth profile.The position deviation of grinding wheel refers to the distance of wheel along the axial or vertical direction deviating from the theoretical position, which can be expressed by Δ and Δ.As Figure 10 shows,  is the theoretical involute,  is the actual tooth profile,   is the radius of rolling circle,  0 is the theoretical tooth space half-angle, and  is the tooth profile rolled angle of the point  on the theoretical tooth profile.According to Figure 10, this article deduces the accurate calculation formula of the tooth profile error caused by the Δ, Δ in order to accurately analyze the grinding tooth error.
When the position deviation of grinding wheel is Δ, Δ, the actual tooth profile after manufacture is From any point ( 1 ,  1 ) of the actual tooth profile tangent for the rolling circle,  as the tangent point, the angle between the straight line  and the axis  is The angle between the line  and  is Let the intersection of the straight line  and the theoretical tooth profile be (, ), and it can be seen from Figure 11 that Tooth space half-angle  0 is known, you can get the point  of the tooth profile rolled angle that The curvature radius of the point  is The length of the line segment  is  We can get the tooth profile error at any point (, ) on the theoretical tooth profile, that is, Based on the derivation and discussion of the equation, the tooth profile deviation is deeply studied in this article.The basic parameters of the gears are shown in Table 2, and the rest are the standard parameters.
Radial error caused by the tooth profile error curve is shown in Figure 12.Vertical error caused by the tooth profile error curve is as shown in Figure 13.

Tooth Profile Error Caused by the Inclination Deviation of the Grinding Wheel
In the actual design, as shown in Figure 14, the cross section of gear is commonly used to be described by polar coordinates.When discussing involute spur gear, making the helical angle to 0 that can meet the requirements.
From Figure 14, the polar form of each point on the section of the transverse can be transformed into In the formula,   is the polar angle of the point  on the section of the gear transverse;   =   (  ) is the vector length of the point  on the section of the gear transverse, which is a function of the polar angle   .
The polar equation of helix surface is represents the rotation angle of the curve around the axis , and  represents the parameter variable corresponding to the position of each point on the surface.Taking ( 26) into (27) to get the gear helical equation is The normal vector on the helical surface of the gear is projected on each axis as is the angle between the vector direction of the point  on the section of the gear transverse and the tangent direction of the section of the gear transverse.  ,   can be obtained by the following two formulas: The formula for calculating the contact line is Substituting ( 28) and (29) into the contact line condition (32), after transformation and simplification, the contact line condition can be obtained as In the formula, represents the center distance of gear.In order to meet the needs of engineering appliance, the error estimate of solution of matrix equation is less than 0.0001.According to the above formula, let the spiral angle be 0; then the deviation of the grinding wheel can be obtained.
When the wheel deviation is 0.1 ∘ , the tooth profile deviation is shown in Figure 15.
The error that occurs when the wheel is deflected from 0 degrees to 1 degree is shown in Figure 16.

Thermal Error Compensation Model for Machine Tool
The developed model must have the ability to deal with nonlinear problems and the ability to accurately predict thermal errors, which is the key to thermal error compensation.At present, neural networks and grey theory are used to model the thermal error of machine tools.
In this paper, a model based on partial least squares neural network is presented, which takes full advantage of partial least squares regression in solving the multiple linear problem.At the same time, this method is improved by importing the neural network algorithm, the parameters of partial least squares regression are trained by this learning algorithm, which makes the new model have the ability of dealing with nonlinear problems.

Partial Least Squares Neural Network Modeling Method.
Partial least squares regression (PLSR) was first presented by Wold to deal with the multiple linear problem of variables in data analysis.This approach is a bit similar to principal component analysis (PCA), but the advantage of PLSR is that both independent variable and dependent variables are used in the process, not just for independent variables.The specific theory of partial least squares regression can refer to literature [19][20][21][22][23][24][25][26][27][28].
Partial least squares regression, although it can deal with the problem of multiple linear between independent variables, is a linear algorithm.For the nonlinearity of machine tool thermal error, this method needs to be improved.We introduce the learning rules of neural network to establish the algorithm of partial least squares neural network.
The structure of partial least squares neural network can be shown in Figure 17; different stages of calculation correspond to different PLS factors ℎ and are determined by the optimal number of factors ℎ ∘ .
By Partial Least Squares Neural Network (PLSNN) method, the Hebbian learning plan in neural networks is used to calculate ŵℎ , Vℎ , and bℎ .The error function and learning algorithm are as follows:    (36) Among them, coefficient of learning rate   ,   ,  V satisfies 0 <   < (∑  =1  2  ) −1 and   ≥   ≥  V .  ,   is positive symmetric weight matrix,   is positive rational number, and (•) is derivative of Nonlinear Weight Function.

Experimental Results
From this we know that we can judge the tooth surface accuracy by measuring the gear.The gear is grinded and measured under different ambient temperatures, and the thermal error of the machine could be obtained indirectly.At a certain temperature, the tooth profile deviation of the gear is shown in Figure 18.
In this paper, the temperature distribution of machine tool is measured by a noncontact thermal infrared imager.When the machine tool is on a working state, the grinding test of gear every 2 ∘ C, according to the change of the ambient temperature and the precision of the tooth profile, is measured.The method of partial least squares neural Although the model of partial least squares neural network has good predictive ability for modeling data, the ability of predicting thermal error is the key which can be applied in thermal error compensation.Therefore, in order to verify the prediction ability of the model for thermal error, a set of thermal error data was measured and used as validation data of thermal error data.As space is limited, only the validation data for radial deviations is shown in Figure 22.
It can be obtained from the diagram that the maximum value of the predicted residual is 0.85m, and the model has good prediction ability for the thermal error of the machine tool.According to the experimental results and fitting curve data, the thermal error of machine tools is compensated for different ambient temperature.After compensation, the gear grinding results are shown in Figure 23, and the accuracy level of gear teeth reaches 4.
In order to make further study on the performance of the model, taking the radial deviation as the comparison object,    3.
The partial least squares neural network modeling method presented in this paper is superior to the commonly used method in fitting ability, and it also has good prediction accuracy and robustness, which is shown in Table 3.

Conclusions
In this paper, the thermal error of the high-precision CNC forming grinding machine was tested.The thermal error of the machine tool was analyzed by the formation of tooth profile deviation.The partial least squares neural network method was presented to model and predict the thermal error.Field test data show that, at different ambient temperatures, after thermal error compensation, the accuracy of the gear can reach four levels.To solve the problem of thermal error modeling and prediction in high-precision CNC forming grinding machine, we can reach the following conclusions.
(1)By measuring the thermal error indirectly, the thermal error of the machine tool under load can be measured.A simple mathematical model can achieve high prediction accuracy.
(2) In this paper, the relationship between the machine tool thermal error and the ambient temperature is studied by taking the ambient temperature as a variable.This method provides a new method and train of thought for the thermal error compensation of the machine tool and provides reference for the factory production.
(3) After compensating by this method, the precision of the gears processed under different ambient temperatures reaches four levels, which provides a new method for highprecision gear grinding.

Figure 1 :
Figure 1: Diagrammatic sketch of the three sections of one tooth.

Figure 3 :
Figure 3: Six pairs of motion errors.

Figure 5 :
Figure 5: The diagram of homogeneous coordinate transformation.

Figure 6 :
Figure 6: The coordinate system of high-precision CNC forming grinding machine setting.

Figure 7 :
Figure 7: The transformation of coordinates system about highprecision CNC forming grinding machine.

Figure 11 :
Figure 11: Wheel along the y direction deviation of the error (radial error).

Figure 12 :
Figure 12: Radial error caused by the tooth profile error curve.

Figure 13 :Figure 14 :
Figure 13: Vertical error caused by the tooth profile error curve.

Figure 15 :
Figure 15: The tooth profile error when the wheel deviation is 0.1 ∘ .

Figure 16 :
Figure 16: The influence of the angle error of the grinding wheel on the tooth profile error.

Figure 19 :
Figure 19: Thermal error (radial deviation) of measurement and model fitting.

Figure 20 :
Figure 20: Thermal error (vertical deviation) of measurement and model fitting.
and using other common methods to model and predict thermal error.The three models are partial least squares regression model (PLSRM), multiple regression model (MRM), and neural network model (NNM), and the performance of each model is compared as shown in Table

Figure 21 :Figure 22 :
Figure 21: Thermal error of measurement and model fitting (rotating error).

Table 1 :
The error of sliding pair.

Table 3 :
The performance comparison of different models.