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

Authors: Xinrong Liu; Zhonghou Wang
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

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.


---

## Body

## 1. 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–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 Figure1. 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 high-performance gear, grinding wheel accuracy is very important. Ding et al. made further research in the related fields [5–9].Figure 1
Diagrammatic sketch of the three sections of one tooth.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 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–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.


## 2. 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 Figure2. In order to grasp the principle contradictions about thermal error for NC machine tool, some assumptions are made as follows.Figure 2
Machine tool structure.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.


## 3. Error Analysis of Moving Parts in High-Precision Computer Numerical Control Form Grinding Machine

### 3.1. The Error Analysis about the Slide Pair Kinematic Parameters

It will have the corresponding error when the high-precision computer numerical control forming grinding 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δxx,t, the straightness error of the Y direction is δyx,t, the straightness error of the Z direction is δzx,t, the rolling error around X direction is εxx,t, the rolling error around Y direction is εyx,t, and the rolling error around Z direction is εzx,t, as shown in Figure 3. There are 2 main prismatic pairs in the high-precision 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.Table 1
The error of sliding pair.
Direction X Direction Y Straightness error of the X direction δ x x , t δ x y , t Straightness error of the Y direction δ y x , t δ y y , t Straightness error of the Z direction δ z x , t δ z y , t Rolling error around X direction ε x x , t ε x y , t Rolling error around Y direction ε y x , t ε y y , t Rolling error around Z direction ε z x , t ε z y , tFigure 3
Six pairs of motion errors.


### 3.2. 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εxφ,t, rotation angle error is εyφ,t, rotation angle error is εzφ,t, translation error is δxφ,t, translation error is δyφ,t, and translation error is δzφ,t, as shown in Figure 4. There are 2 main revolute pairs in the high-precision computer numerical control forming gear grinding machine.Figure 4
Six motion errors of revolute pairs.


### 3.3. 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 errorSxy 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.


## 3.1. The Error Analysis about the Slide Pair Kinematic Parameters

It will have the corresponding error when the high-precision computer numerical control forming grinding 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δxx,t, the straightness error of the Y direction is δyx,t, the straightness error of the Z direction is δzx,t, the rolling error around X direction is εxx,t, the rolling error around Y direction is εyx,t, and the rolling error around Z direction is εzx,t, as shown in Figure 3. There are 2 main prismatic pairs in the high-precision 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.Table 1
The error of sliding pair.
Direction X Direction Y Straightness error of the X direction δ x x , t δ x y , t Straightness error of the Y direction δ y x , t δ y y , t Straightness error of the Z direction δ z x , t δ z y , t Rolling error around X direction ε x x , t ε x y , t Rolling error around Y direction ε y x , t ε y y , t Rolling error around Z direction ε z x , t ε z y , tFigure 3
Six pairs of motion errors.


## 3.2. 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εxφ,t, rotation angle error is εyφ,t, rotation angle error is εzφ,t, translation error is δxφ,t, translation error is δyφ,t, and translation error is δzφ,t, as shown in Figure 4. There are 2 main revolute pairs in the high-precision computer numerical control forming gear grinding machine.Figure 4
Six motion errors of revolute pairs.


## 3.3. 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 errorSxy 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.


## 4. 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.


### 4.1. Homogeneous Coordinate Transformation

In the three-dimensional coordinate, the calculation methods of homogeneous coordinate transformation are as follows, which is shown in Figure5. Hypothetical space has a coordinate system OAXAYAZA, the coordinates of a point P in the coordinate system are OAXAYAZA, when coordinate system OAXAYAZA moves x along the XA axis to get a new coordinate system OBXBYBZB, the relationship between the coordinates of point P in the new coordinate system (xB,yB,zB) and the coordinates in the original coordinate system (xA,yA,zA) can be expressed as follows:Figure 5
The diagram of homogeneous coordinate transformation.Three-dimensional translation homogeneous coordinate transformation matrix ofTrans(x) expressed among XA axis to translate x is(1)xAyAzA1=100x010000100001xByBzB1=TransxxByBzB1(2)Transx=100x010000100001Similarly, the three-dimensional homogeneous coordinate transformation matrix that translatesy along YA axis and z along ZA axis, respectively, is(3)Transy=1000010y00100001(4)Transz=10000100001z0001When the coordinate systemOAXAYAZA is rotated around the XA axis εx to obtain a new coordinate system OBXBYBZB, the relationship between the coordinates of the point P in the new coordinate system (xB,yB,zB) and the coordinates in the original coordinate system (xA,yA,zA) can be expressed as(5)xAyAzA1=10000cos⁡εx-sin⁡εx00sin⁡εxcos⁡εx00001xByBzB1=RotεxxByBzB1whereRot(εx) represents a three-dimensional homogeneous coordinate transformation matrix rotated εx about the XA axis:(6)Rotεx=10000cos⁡εx-sin⁡εx00sin⁡εxcos⁡εx00001Similarly, three-dimensional homogeneous coordinate transformation matrix around theYA axis rotation εy and ZA axis rotation εz, respectively, is(7)Rotεy=cos⁡εy0sin⁡εy00100-sin⁡εy0cos⁡εy00001(8)Rotεz=cos⁡εz-sin⁡εz00sin⁡εzcos⁡εz0000100001When the coordinate systemOAXAYAZA, respectively, transforms x, y, z along the XA, YA, ZA axis and then revolves εx, εy, εz, respectively, around XA, YA, ZA to obtain a new coordinate system OBXBYBZB, then the three-dimensional homogeneous coordinate transformation matrix of coordinate system OAXAYAZA to the coordinate system OBXBYBZB is(9)TBA=Transx×Transy×Transz×Rotεx×Rotεy×RotεzTake (2), (3), (4), (6), (7), and (8) into (9):(10)TBA=cos⁡εycos⁡εz-cos⁡εysin⁡εzsin⁡εyxsin⁡εxsin⁡εycos⁡εz+cos⁡εxsin⁡εz-sin⁡εxsin⁡εysin⁡εz+cos⁡εxcos⁡εz-sin⁡εxcos⁡εyy-cos⁡εxsin⁡εycos⁡εz+sin⁡εxsin⁡εzcos⁡εxsin⁡εysin⁡εz+sin⁡εxcos⁡εzcos⁡εxcos⁡εyz0001When the rotation angleεx, εy, εz is very small, there is sin⁡εx≈εx, sin⁡εy≈εy, sin⁡εz≈εz, cos⁡εx≈1, cos⁡εy≈1, cos⁡εz≈1. When the translation displacement is very small, use δx, δy, δz to separate table x, y, and z, respectively. Ignoring the second-coordinate transformation matrix is expressed as(11)TBA=1-εzεyδxεz1-εxδy-εyεx1δz0001


### 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 Figure6: 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.Figure 6
The coordinate system of high-precision CNC forming grinding machine setting.For the convenience of calculations, we placed the origin of base coordinates and local coordinates in the zero pointO, and the same axis direction remains same. The two branches of the coordinate system drive chain are shown in Figure 7.Figure 7
The transformation of coordinates system about high-precision CNC forming grinding machine.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 matrixTWG of the gear coordinate system G to the grinding wheel coordinate system W is(12)TWG=TAG·TOA·TXO·TYX·TBY·TSB·TWSBased on the formation principle of different kinds of errors, the coordinate transformation matrix of the gear coordinate systemG to the grinding wheel coordinate system W should be(13)TeWG=TeAG·TeOA·TeXO·TeYX·TeBY·TeSB·TeWSErratic movements are indicated by superscriptse. In reality, the coordinate transformation matrix of the gear coordinate system G to the grinding wheel coordinate system W is equal to the multiplication of the ideal motion matrix and an error matrix EWG, that is,(14)TeWG=EWG·TWGBased on the premise of small errors, the error motion conversion matrix of the gear coordinate systemG to the grinding wheel coordinate system W can be set as(15)EWG=1-εzεyδxεz1-εxδy-εyεx1δz0001δ x, δy, δz represent the relative error between the actual machining position of the grinding wheel and the ideal machining position; εx, εy, εz 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 EGW:(16)EWG=TeAG·TeOA·TeXO·TeYX·TeBY·TeSB·TeWS·TWS-1·TSB-1·TBY-1·TYX-1·TXO-1·TOA-1·TAG-1T  W S represents the coordinate conversion matrix of the wheel spindle coordinate system S to the wheel coordinate system W. The grinding wheel is installed in the spindle, moving together with the spindle. Therefore, the coordinate transformation matrix between them is unit matrix I, TWS=I.Due to space constraints, there is seldom detailed description for the error solution procedure. Solving the matrix equation by MATLAB and neglecting the high-order infinitesimal small quantities, we finally obtain the comprehensive error model. Wheel angle deviations in three directions areεx,εy,εz.Deviations in three directions are δx,δy,δz.


### 4.3. 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 Figure8. 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.Figure 8
Wheel thermal error direction diagram.


## 4.1. Homogeneous Coordinate Transformation

In the three-dimensional coordinate, the calculation methods of homogeneous coordinate transformation are as follows, which is shown in Figure5. Hypothetical space has a coordinate system OAXAYAZA, the coordinates of a point P in the coordinate system are OAXAYAZA, when coordinate system OAXAYAZA moves x along the XA axis to get a new coordinate system OBXBYBZB, the relationship between the coordinates of point P in the new coordinate system (xB,yB,zB) and the coordinates in the original coordinate system (xA,yA,zA) can be expressed as follows:Figure 5
The diagram of homogeneous coordinate transformation.Three-dimensional translation homogeneous coordinate transformation matrix ofTrans(x) expressed among XA axis to translate x is(1)xAyAzA1=100x010000100001xByBzB1=TransxxByBzB1(2)Transx=100x010000100001Similarly, the three-dimensional homogeneous coordinate transformation matrix that translatesy along YA axis and z along ZA axis, respectively, is(3)Transy=1000010y00100001(4)Transz=10000100001z0001When the coordinate systemOAXAYAZA is rotated around the XA axis εx to obtain a new coordinate system OBXBYBZB, the relationship between the coordinates of the point P in the new coordinate system (xB,yB,zB) and the coordinates in the original coordinate system (xA,yA,zA) can be expressed as(5)xAyAzA1=10000cos⁡εx-sin⁡εx00sin⁡εxcos⁡εx00001xByBzB1=RotεxxByBzB1whereRot(εx) represents a three-dimensional homogeneous coordinate transformation matrix rotated εx about the XA axis:(6)Rotεx=10000cos⁡εx-sin⁡εx00sin⁡εxcos⁡εx00001Similarly, three-dimensional homogeneous coordinate transformation matrix around theYA axis rotation εy and ZA axis rotation εz, respectively, is(7)Rotεy=cos⁡εy0sin⁡εy00100-sin⁡εy0cos⁡εy00001(8)Rotεz=cos⁡εz-sin⁡εz00sin⁡εzcos⁡εz0000100001When the coordinate systemOAXAYAZA, respectively, transforms x, y, z along the XA, YA, ZA axis and then revolves εx, εy, εz, respectively, around XA, YA, ZA to obtain a new coordinate system OBXBYBZB, then the three-dimensional homogeneous coordinate transformation matrix of coordinate system OAXAYAZA to the coordinate system OBXBYBZB is(9)TBA=Transx×Transy×Transz×Rotεx×Rotεy×RotεzTake (2), (3), (4), (6), (7), and (8) into (9):(10)TBA=cos⁡εycos⁡εz-cos⁡εysin⁡εzsin⁡εyxsin⁡εxsin⁡εycos⁡εz+cos⁡εxsin⁡εz-sin⁡εxsin⁡εysin⁡εz+cos⁡εxcos⁡εz-sin⁡εxcos⁡εyy-cos⁡εxsin⁡εycos⁡εz+sin⁡εxsin⁡εzcos⁡εxsin⁡εysin⁡εz+sin⁡εxcos⁡εzcos⁡εxcos⁡εyz0001When the rotation angleεx, εy, εz is very small, there is sin⁡εx≈εx, sin⁡εy≈εy, sin⁡εz≈εz, cos⁡εx≈1, cos⁡εy≈1, cos⁡εz≈1. When the translation displacement is very small, use δx, δy, δz to separate table x, y, and z, respectively. Ignoring the second-coordinate transformation matrix is expressed as(11)TBA=1-εzεyδxεz1-εxδy-εyεx1δz0001


## 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 Figure6: 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.Figure 6
The coordinate system of high-precision CNC forming grinding machine setting.For the convenience of calculations, we placed the origin of base coordinates and local coordinates in the zero pointO, and the same axis direction remains same. The two branches of the coordinate system drive chain are shown in Figure 7.Figure 7
The transformation of coordinates system about high-precision CNC forming grinding machine.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 matrixTWG of the gear coordinate system G to the grinding wheel coordinate system W is(12)TWG=TAG·TOA·TXO·TYX·TBY·TSB·TWSBased on the formation principle of different kinds of errors, the coordinate transformation matrix of the gear coordinate systemG to the grinding wheel coordinate system W should be(13)TeWG=TeAG·TeOA·TeXO·TeYX·TeBY·TeSB·TeWSErratic movements are indicated by superscriptse. In reality, the coordinate transformation matrix of the gear coordinate system G to the grinding wheel coordinate system W is equal to the multiplication of the ideal motion matrix and an error matrix EWG, that is,(14)TeWG=EWG·TWGBased on the premise of small errors, the error motion conversion matrix of the gear coordinate systemG to the grinding wheel coordinate system W can be set as(15)EWG=1-εzεyδxεz1-εxδy-εyεx1δz0001δ x, δy, δz represent the relative error between the actual machining position of the grinding wheel and the ideal machining position; εx, εy, εz 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 EGW:(16)EWG=TeAG·TeOA·TeXO·TeYX·TeBY·TeSB·TeWS·TWS-1·TSB-1·TBY-1·TYX-1·TXO-1·TOA-1·TAG-1T  W S represents the coordinate conversion matrix of the wheel spindle coordinate system S to the wheel coordinate system W. The grinding wheel is installed in the spindle, moving together with the spindle. Therefore, the coordinate transformation matrix between them is unit matrix I, TWS=I.Due to space constraints, there is seldom detailed description for the error solution procedure. Solving the matrix equation by MATLAB and neglecting the high-order infinitesimal small quantities, we finally obtain the comprehensive error model. Wheel angle deviations in three directions areεx,εy,εz.Deviations in three directions are δx,δy,δz.


## 4.3. 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 Figure8. 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.Figure 8
Wheel thermal error direction diagram.


## 5. 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 Figure9. Therefore, the main consideration of this manuscript is the influence of thermal error on the accuracy of tooth profile.Figure 9
Gear tooth direction measurement results.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Δx and Δy. As Figure 10 shows, ab¯ is the theoretical involute, cd¯ is the actual tooth profile, rb is the radius of rolling circle, θ0 is the theoretical tooth space half-angle, and θ is the tooth profile rolled angle of the point k on the theoretical tooth profile. According to Figure 10, this article deduces the accurate calculation formula of the tooth profile error caused by the Δx, Δy in order to accurately analyze the grinding tooth error.Figure 10
Error along the x direction of the wheel (vertical error).Set the equation of theoretical involute of the right tooth profile is(17)x=rbsin⁡θ+θ0-θcos⁡θ+θ0y=rbcos⁡θ+θ0+θsin⁡θ+θ0When the position deviation of grinding wheel isΔx, Δy, the actual tooth profile after manufacture is(18)x1=rbsin⁡θ+θ0-θcos⁡θ+θ0+Δxy1=rbcos⁡θ+θ0+θsin⁡θ+θ0+ΔyFrom any pointMx1,y1 of the actual tooth profile tangent for the rolling circle, N as the tangent point, the angle between the straight line OM and the axis y is(19)θ1=tan-1⁡x1y1The angle between the lineOM and ON is (20)θ2=cos-1⁡rbx12+y12Let the intersection of the straight lineOM and the theoretical tooth profile be k(x,y), and it can be seen from Figure 11 that(21)θ0+θ=θ1+θ2Figure 11
Wheel along the y direction deviation of the error (radial error).Tooth space half-angleθ0 is known, you can get the point k of the tooth profile rolled angle that(22)θ=θ1+θ2-θ0The curvature radius of the pointk is (23)NK=rbθThe length of the line segmentNM is(24)NM=x12+y12-rb2We can get the tooth profile error at any pointkx,y on the theoretical tooth profile, that is,(25)ΔJ=NM-NK=x12+y12-rb2-rbθ=x12+y12-rb2-rbtan-1⁡x1y1+cos-1⁡rbx12+y12-θ0Based 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 Table2, and the rest are the standard parameters.Table 2
Gear basic parameters.
Parameter Value Number of teeth 35 Modulus 2.5mm pressure angle 20° Tooth width 20mmRadial error caused by the tooth profile error curve is shown in Figure12. Vertical error caused by the tooth profile error curve is as shown in Figure 13.Figure 12
Radial error caused by the tooth profile error curve.Figure 13
Vertical error caused by the tooth profile error curve.


## 6. Tooth Profile Error Caused by the Inclination Deviation of the Grinding Wheel

In the actual design, as shown in Figure14, 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.Figure 14
The cross section of gear.From Figure14, the polar form of each point on the section of the transverse can be transformed into(26)z1u=rycos⁡δyz2u=rysin⁡δyz3u=0In the formula,δy is the polar angle of the point M on the section of the gear transverse; ry=ryδy is the vector length of the point M on the section of the gear transverse, which is a function of the polar angle δy.The polar equation of helix surface is(27)x=x0ucos⁡θ-y0usin⁡θy=x0usin⁡θ+y0ucos⁡θz=pθθ represents the rotation angle of the curve around the axis Z, and u 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(28)x1=rycos⁡δy+θx2=rysin⁡δy+θx2=pθThe normal vector on the helical surface of the gear is projected on each axis as(29)nx1=psin⁡τcos⁡ξynx2=-pcos⁡τcos⁡ξynx3=ryξ y is the angle between the vector direction of the point M on the section of the gear transverse and the tangent direction of the section of the gear transverse. δy, ξy can be obtained by the following two formulas:(30)ξy=arctan⁡rydry/dδy=arctan⁡rydδydry(31)τ=δy+ξy+θThe formula for calculating the contact line is(32)zgnx+angcotΣ+a-xg+pcotΣnz=0Substituting (28) and (29) into the contact line condition (32), after transformation and simplification, the contact line condition can be obtained as(33)n2cos⁡τ-n1sin⁡τ-τ-n3=0In the formula,(34)k1=a+pcotΣp2k2=acotΣpm=rycosξyv=rysinξyn1=k1mn2=k2+m2p2n3=-mvp2-δy-ξya 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 Figure15.Figure 15
The tooth profile error when the wheel deviation is 0.1°.The error that occurs when the wheel is deflected from 0 degrees to 1 degree is shown in Figure16.Figure 16
The influence of the angle error of the grinding wheel on the tooth profile error.


## 7. 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.


### 7.1. 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–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 Figure17; different stages of calculation correspond to different PLS factors h and are determined by the optimal number of factors h∘.Figure 17
Partial least squares neural network structure.By Partial Least Squares Neural Network (PLSNN) method, the Hebbian learning plan in neural networks is used to calculatew^h, v^h, and b^h. The error function and learning algorithm are as follows:(35)ewai=ai-ciwhebai=ai-t^hiwhevci=ci-vht^hi(36)w^hk+1=w^hk+μwSwcigai-ciw^hkb^hk+1=b^hk+μbSbt^higai-t^hib^hkv^hk+1=v^hk+μvSvt^higai-t^hiv^hkAmong them, coefficient of learning rateμw,μb,μv satisfies 0<μw<∑i=1mci2-1 and μw≥μb≥μv. Sw,Sb is positive symmetric weight matrix, Sy is positive rational number, and g(•) is derivative of Nonlinear Weight Function.


## 7.1. 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–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 Figure17; different stages of calculation correspond to different PLS factors h and are determined by the optimal number of factors h∘.Figure 17
Partial least squares neural network structure.By Partial Least Squares Neural Network (PLSNN) method, the Hebbian learning plan in neural networks is used to calculatew^h, v^h, and b^h. The error function and learning algorithm are as follows:(35)ewai=ai-ciwhebai=ai-t^hiwhevci=ci-vht^hi(36)w^hk+1=w^hk+μwSwcigai-ciw^hkb^hk+1=b^hk+μbSbt^higai-t^hib^hkv^hk+1=v^hk+μvSvt^higai-t^hiv^hkAmong them, coefficient of learning rateμw,μb,μv satisfies 0<μw<∑i=1mci2-1 and μw≥μb≥μv. Sw,Sb is positive symmetric weight matrix, Sy is positive rational number, and g(•) is derivative of Nonlinear Weight Function.


## 8. 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 Figure18.Figure 18
Profile deviation.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 network was used to model the errors in three directions. The measured and fitted values of the radial error are shown in Figure19. The maximum residual error can be calculated as 0.51μm from the residual curve in the figure. The measured and fitted values of the vertical error are shown in Figure 20. The residual error curve in the figure shows the maximum residual error of 0.42μm. The measured value and the fitted value of the angular error caused by the angular deviation are shown in Figure 21. The residual error curve in the figure shows the maximum residual error is 0.48μm, and the model has good fitting ability to the thermal error data used in modeling.Figure 19
Thermal error (radial deviation) of measurement and model fitting.Figure 20
Thermal error (vertical deviation) of measurement and model fitting.Figure 21
Thermal error of measurement and model fitting (rotating error).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 Figure22.Figure 22
Model error prediction capability verification (radial error).It can be obtained from the diagram that the maximum value of the predicted residual is 0.85μm, 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.Figure 23
Compensation gear measurement results.In order to make further study on the performance of the model, taking the radial deviation as the comparison object, 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 Table3.Table 3
The performance comparison of different models.
Model Modeling 1(μm) Modeling 2(μm) PLSRM 0.91 1.42 MRM 0.56 1.68 NNM 0.74 1.13 PLSNN 0.51 0.85The 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 Table3.


## 9. 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 high-precision gear grinding.



---

*Source: 4619624-2018-07-15.xml*