Optimal analysis of TSM fitting considering contact interface TECs of meshing gears

The pre-research project aims to propose a method of optimizing tooth surface modification (TSM) fitting under thermoelastic lubrication (TEL) conditions to reveal the most popular and concerned mechanism issues in the field of mechanical engineering. This is an exploratory study, mainly considering the gear comprehensive error and TSM state of TEL contact interfaces are extremely harsh, which complicates simulation analysis and optimal design. TSM simulation has fitted numerically agreement with optimized results obtained experimentally, which means that they can, whether isolated or by using a multiscale coupling method, start to be adopted to revise, and verify meta-models for thermoelastic characteristics (TECs) under TEL problems. This subject involves the TSM fitting of the theoretical teeth surface superimposed structure, and performs 3D and diagonal modification optimization design, obtains the modified surface position and normal vectors, according to load teeth contact analysis (LTCA), a variety of optimized modification models are established and complex curved surfaces are analyzed to fit the actual gear teeth surface thermoelastic contact numerical simulation, which has been further demonstrated and expressed in the type experiment combined with actual key working conditions and multiple influencing parameters, which is of great significance to the development of modern gear transmission system.


Introduction
The TSM fitting optimization is related to the manufacturing and installation errors and the meshing impact caused by the thermoelastic deformation of alternating meshing. 1 It is an important factor affecting the coupling of dynamic thermoelastic characteristics (TECs) of gears. The uniform teeth surface pressure makes the teeth profile deformation run smoothly, reduces the vibration noise and then attenuates the impact, in order to improve the influence of thermoelastic deformation of teeth surface on gear transmission performance and improve the anti-scuffing bearing capacity of teeth surface, which should not be ignored in the analysis of gear vibration noise and description of meshing gears instantaneous contact interface lubrication. [2][3][4] Considering the characteristics of gear transmission, there are many methods of TSM, and the influencing 1 factors are complex. 5 Many scholars have noticed the traditional gear meshing theory to obtain the modified teeth surface by using the modification cutter are too theoretical, and have not considered the comprehensive error of the gear, the precise geometric parameters of the teeth surface, and the uncertain machining error, it is difficult to meet the requirements of high-precision, high-speed, and high-load design. [6][7][8] TEL is the most complex and common reality online state. Figure 1 reveals the mixed lubrication state of the gear pair in the actual operating conditions during meshing time. Many TEL problems of alternate meshing gear teeth point and line contact belong to this research field. [9][10][11] Subject to the influence of pure elastic deformation between the contact surfaces of gear pairs, most previous scholars have revealed the time-varying law of lubricating oil film thickness along the meshing teeth surface along the finite length contact line. 12 However, the relative influences of lubricating oil viscositytemperature -pressure on film thickness have not been highlighted. 13 A key characterization analysis model for Dowson elastic flow pressure oil film was shown in Figure 2.
In the past few decades, a significant progress has been made in numerical methods, analytical solutions and experimental techniques. 14 TSM technology has been the most effective way to improve the gear teeth meshing performance and reduce vibration and noise of the gear transmission system (GTS), the TSM effect mechanism and modification parameter selection principles have become a consensus. [15][16][17] Moreover, most of these research results are based on the assumption of the static state of GTS, and there are few studies on TEL considering the transient dynamic conditions of GTS meshing. There are two common problems in their research. One is that the precise teeth surface geometry conditions under different error conditions or TSM are not considered, and the actual meshing teeth surface and the theoretical conjugate teeth surface have a magnitude error, the other is to rely on modification tool to obtain the modified teeth surface is derived from the traditional meshing theory, which makes the modified teeth surface too theoretical, and its numerical analysis and uncertainty errors in the machining process are difficult to meet the design requirements of high precision, high speed, and high load bearing performance of modern gears.
The empirical simplified formula has been difficult to adapt to the technical requirements of high-index GTSs. Many scholars have carried out a lot of research work in the field of TSM and obtained results for reference. [18][19][20] Related to the key characterization index of gear mesh stiffness density, the load transmission error (LTE) of the gear pair is further realized from twodimensional teeth profile modification to threedimensional TSM and static trajectory to dynamic target parameter process step.
A lot of research on TSM technology based on stiffness, load and stress distributions have been performed, some factors such as load size, working condition parameters, lubricating oil parameters, and cooling measures have been considered, which will affect the TECs (temperature field, thermoelastic deformation, thermal expansion, and thermoelastic stress) of GTS, which further increases the difficulty in accurately determining the modification curve, modification index, and maximum modification amount.
The significant TECs of temperature field of a heavy-duty and high-speed gear bodies have greatly affected the resistance to thermoelastic deformation of the GTS. [21][22][23] The refinement design and optimization process of TSM considering the influence of gear temperature field is quite difficult, which becomes necessary for GTS to analyze the TECs of TSM. 4,6,24 For high-speed and heavy-duty gears, teeth surface thermoelastic anti-scuffing damage is the main failure form of the GTS. 25 The instantaneous contact temperature rise of the meshing teeth surface is caused by transient overload of GTS, which leads to rapid overheating of teeth  surface and thermoelastic anti-scuffing damage. [26][27][28] The thermoelastic anti-scuffing load-bearing performance evaluation is the core content of deep optimization and meshing process verification of high-speed and heavy-duty gears. 29 In view of the relevant literature on gears thermoelastic anti-scuffing load-bearing performance, the research results show that the meshing gears hybrid lubrication characteristics are closely related to gears thermoelastic anti-scuffing load-bearing performance in time domain. In the following investigation, the optimal design of three-dimensional/diagonal modifications and numerical fitting of modified surface are proposed. Adherence to theoretical teeth surface superposition mechanism, a modified surface fitting evolution has been realized and three dimensional/diagonal modification optimal design processes have been implemented. The gears comprehensive error and TSM status and precise teeth surface geometric parameters have been noted, and then the modified surface position and normal vectors have also been deduced, in view of the multiple TSM optimized meta-models, the thermoelastic contact numerical simulation of complex curved surface fitting the actual teeth surface is analyzed to reduce the transmission error and vibration noise in order to improve meshing gear teeth thermoelastic anti-scuffing load-bearing performance.

TSM fitting and position/normal vectors numerical calculation
In addition, the data of teeth profile and teeth profile can be collected uniformly along the teeth profile m 3 n. The parameters of three-dimensional topological and diagonal modification surfaces are used to describe the teeth profile, teeth orientation, three-dimensional topology and diagonal modification surface d ij (x, y). In which i = 1 . . . m, j = 1 . . . n Based on the cubic B-spline, a smooth modification surface is obtained by fitting the mesh node data of the teeth surface. The teeth profile modification surface are cylindrical surfaces, and the three-dimensional topological and diagonal modification surfaces are the complex surfaces of teeth profile and teeth direction superposition Dabc and Ddef , the surface modified by rotational projection is represented by data d(x, y), the ideal modification surface is realized based on the numerical method of fitting optimization design. The modification amount on the mesh node of the modification area is shown in Figure 3. The distance from point a to BC MP and the distance from point d to EF are defined, in the modification area D ABC, the modification amount along the MP direction increases along the parabola curve, and the free form 5-axis linkage CNC machine tool is used to process according to the parabola twice, then M is any point in the modification area DABC, and the modification amount is as follows: where,y 1 and y 3 respectively represents the modification amount (y 1 tooth top modification amount, y 2 tooth root modification amount), y 2 indicates the length of unmodified area, H describes the working teeth height, L is the working teeth length, l ij depicts the distance between grid node and repair boundary area; l pm , l qn respectively expresses the maximum distance between grid node and repair boundary area.

Characterizations of k-order B-spline curves and surfaces
The concept of B-spline was put forward by Schoenberg, and the recursive definition of B-spline was proposed by De door and Cox respectively, which made the calculation of B-spline stable and simple. [30][31][32] Assumed here, B-spline curve of order k has the following characteristics: (1) In the definition domain of k-order B-spline curve, every curve segment is infinitely differentiable, that is, infinite order continuous differentiable. (2) The convex hull characteristics of curve is defined as the union of control points of each curve segment k-orders B-spline curves are placed in the union of convex hull. If the secondary B-spline curve is a straight line segment, then it is a common vertex in order k + 1.
(3) The number of intersections between k-order B-spline curve and control polygon is less than that of any plane and its intersection points, and does not include the whole control polygon plane, namely, attenuation variation characteristics. (4) The larger the distance between k-orders Bspline curves indicates that the larger the control polygon plane is, the higher its order is, but it does not contain collinear vertices, also known as, polishing characteristics. (5) Geometric and affine invariant characteristics.
Here, considering the non-uniform partition of k-order B-spline, its basis function is given where, d i Is the control point, P 1 P 2 Á Á Á P n is the control vertex k À 1 order B-spline basis function of order N i, k (t) by Among them, t i is a node, and if t jÀ1 \t j = t j + 1 = Á Á Á t j + 1Àm \t j + 1 , said t j , t j + 1 , :::, t j + 1Àm , is the node vector T of m duplicate node.
Then B-spline curves and surfaces are characterized as Considering the uniform division of the parameter space (u, v), the cubic uniform B-spline curve and contact interface are represented as where P i, j is (m + 1) 3 (n + 1) points forming the control polygon mesh, u and v is the parameter of two directions in the rectangular field and its times are respectively k and L.

Determination of position/normal vectors of gear modification surface
In the actual modified surface, the analytic control points (or control points) are inversely obtained by fitting and adjusting the space value points for many times, and then the inverse calculation process of the difference between the two surfaces is obtained. N unordered value points are set P 1 , P 2 , Á Á Á , P n , resolve B-spline curve to control polygon vertex V i (0 1 , 1, Á Á Á , n + 1). Then the correlation between the value points and the control points of the space type is obtained.
Set spatial value point P ij (i = 1, 2, Á Á Á , n; j = 1, 2, Á Á Á , m), and analyze the corresponding bicubic spline surface control points Considering the influence of gear manufacturing, installation error, load distribution between teeth and TSM, the alternative meshing teeth surface is non-conjugate. [33][34][35] That is, any meshing position u, when two pairs of teeth are meshed at the same time, the transmission error will be produced d tI , d tII , the normal linear displacement of teeth surface is also called the normal vector of teeth surface in view of the comprehensive error of bearing transmission of gear pair, as shown in Figure 4.
where, r b is the base circle radius of the driven gear, b b is the helix angle of the base circle.
Represented in Figure 4, M 0 expresses the transient meshing point of two teeth surfaces S 1 and S 2 (based on the fixed coordinate system S f ), n f ! is the common normal vector of the two teeth surfaces (at the meshing point M 0 ), M j represents any discrete point (along the long axis of the instantaneous contact ellipse), then the position vector at the meshing point M j in the fixed coordinate system S f .
The parallel line n f ! passing through the meshing point M j intersects the two teeth surfaces at points M j1 and M j2 . M j1 M j2 is the normal vector at the discrete point M j of the two contact interfaces. Then cross the meshing point M j parallel line n f ! equation expression x À x j n x = y À y j n y = z À z j n z ð10Þ Set is the equation of teeth surface S 2 ,where, u 1 , u 1 , u 2 , u 2 describes the two teeth surfaces parameters, u 1 is the rotation angle of gear 1, and u 2 is the rotation angle of gear 2, referring to the fixed coordinate system S f .M j1 (x j1 , y j1 , z j1 ) and M j2 (x j2 , y j2 , z j2 ) are the coordinates of the intersection of the straight line L and the two teeth surfaces, then the normal vector expression of the two teeth surfaces When calculating the actual normal vector of teeth surface, it is necessary to fit the meshing data of teeth surface deviation, especially with large load and high precision. The gear surface deviation should not be more than 50 mm, and the fitting accuracy of gear modification surface should not be greater than 1.0 mm. Based on the cubic B-spline fitting teeth surface meshing deviation and node data, assuming that the meshing time domain is fixed fitting spline surface boundary and collecting spline data, and extending the rectangular grid boundary data to the whole teeth surface boundary, the analytic control matrix equation is inversely solved and the change is little compared with the adjacent boundary surface, so as to ensure the accuracy of teeth surface contact analysis.

Optimization analysis of TSM based on genetic algorithm
The optimization of modified teeth surface is an iterative process of solving nonlinear contact problem, which is also called the process of changing the state of modified teeth surface to improve the contact condition of gear meshing surface. 34 Set the gear pair modification optimization variable as its curve parameter, the LTCA numerical analysis result is the optimization objective function. This topic is based on genetic algorithms (GA) to simulate the genetic and evolutionary processes of organisms in the natural environment, and then searches for adaptive global optimization probabilities, reproduces and optimizes populations and converges to the most suitable environmental individuals. The GA process is shown in Figure 5, and its steps can be described.
Step1. Select the objective function and determine the variable domain and coding accuracy to form the code. Step2. Perform crossover operations on individuals selected to enter the matching pool to form a new population. Step3. Select population individuals with a small probability to perform mutation operations to form a new population. Step4. Calculate individual fitness based on fitness function and select probability to form a new population of individuals. Step5. If the check conditions are met, the GA ends, and the individual with the highest current population fitness is the expected solution.
Considering that the B-spline surface uses nonparametric implicit representation of the smooth connection of multiple surfaces, the relationship between variables in the defined domain is described. If the slope of a point in the point set tends to be parallel to the axis of coordinates, the slope at this point is close to infinity. In view of the fact that curves and surfaces in a non-identical plane are mathematically expressed by parametric equation. 35 In the LTCA equation, the parameter ua, ub, La, Lb, u2, the gear parameter ua, the tool parameter La, the gear parameter ua and the tool parameter La are mapped to the point (x, y) on the rotating projection surface, and then the space (x, y) is mapped to the spline surface parameter space (u, v), the three-coordinate modified surface mesh data are all within the actual teeth surface. With the help of Matlab to solve the nonlinear equations, the iterative analysis gear parameter ua and tool parameter La exceed the boundary of the teeth surface, and the modified surface must be extrapolated.
Based on the data of the rectangular grid, the control vertices are inversely calculated and kept unchanged. The four boundaries of the space (x, y) rectangular grid are naturally extended to the left and right boundaries of the teeth tip, teeth root, and teeth side, and then mapped to the parameter space domain (u, v) of the spline surface. The optimization analysis process of LTCA for TSM is shown in Figure 6.
The gear system transmits power and the meshing teeth bear load. The thermoelastic deformation of the teeth surface causes the distortion of the teeth profile and teeth direction, resulting in impact, vibration, and unbalanced load in the process of alternate meshing. Based on the mechanism of teeth profile modification and teeth alignment modification of involute gear, the time-varying law of thermoelastic deformation of alternate meshing contact interface with meshing position is considered, the distribution coefficient of transient contact dynamic load is shown in Figure 7.
The maximum value of contact dynamic load on the teeth during the alternate meshing process of gear pair is placed in the meshing in and meshing out transients in the single teeth meshing area, which are also shown in the Figure 7C and D, the initial point of the gear teeth is deformed greatly points B 1 and endings B 2 , the reason of the meshing is the interference process of gears C, the point is caused by the maximum transient contact dynamic load B 2 , if the point meshing interference occurs, the top of the driving gear (the position  where the gear teeth are engaged) must be modified, and the maximum modification amount is set to be d B2 , point D maximum transient contact dynamic load leads to point B 1 , in case of meshing interference, the top of driven gear teeth (teeth engagement position) must be modified, and the maximum modification amount is d B1 , considering the influence of gear comprehensive error (base pitch error, teeth profile error, etc.) on the transmission system, the teeth profile modification amount of meshing-in and meshing-out position is taken into account D in and D out where, d C1 is the thermoelastic deformation at point C of the driving gear, d C2 is the thermoelastic deformation of the driven gear at point C, d D1 is the thermoelastic deformation at point D of the driving gear, d D2 is the thermoelastic deformation at point D of the driven gear, Df b is the pitch error caused by the comprehensive error of the gear, which is selected as positive or negative according to its direction. Taking into account the modification of involute teeth profile in the alternate meshing area of double teeth, the ideal modification curve is the change law of the modification amount increasing to the maximum value. The calculation formula of the modification curve is expressed as follows where, D max is the maximum modification amount of gear teeth,x is the distance between the critical point of alternating meshing area of single and double teeth, and is also the distance from any position on the segment B 1 C to point C, L is the length from the critical point of single and double teeth contact area on the alternate meshing line to the starting point or end point of meshing, b is the contact load change index between the pair of meshing teeth, and its value is in the range of 1:0;2:0, in order to reduce the sudden change of contact load and the alternating meshing gear interface impact. Gear parameters and teeth profile shape are affected by load distribution and contact pressure distribution between teeth and relative sliding of contact teeth to meshing teeth surface. Proper amount of addendum trimming along normal direction of driving and driven gear profile improves load contact pressure distribution between teeth, reduces relative sliding of meshing teeth   surface, and optimizes contact temperature of teeth surface. 36 Taking into account the TECs of teeth surface (teeth contact deformation and bending deformation, shear deformation, and teeth root deformation, etc.), and considering the distribution of instantaneous contact dynamic load (teeth meshing stiffness) and the time-varying law of teeth contact dynamic stress and teeth root bending dynamic stress, an estimated empirical formula of the modified curve is as follows where d a is the thermoelastic deformation of the teeth profile, v t is the instantaneous contact dynamic load per unit teeth width (N =mm), v t = F t =b, F t is the tangential force on the dividing circle of the gear (N ), and c g is the meshing stiffness (N=(mm Á mm)).
Teeth profile and teeth alignment modification also affects the load distribution between teeth in meshing state. Several pairs of teeth are engaged in alternate meshing and share the load (contact ratio e a .1). The loaded teeth produce bending, shearing, contact deformation, and thermoelastic deformation, which inevitably leads to uneven load distribution between the teeth in the same meshing state. The load distribution coefficient (z) of spur gear with different modification index (b) is given considering the coupling of thermal and elastic characteristics of teeth surface, as shown in Figure 8.
It can be seen from Figure 8 that the load of unmodified gear changes suddenly at the meshing in end and meshing out end, and the load distribution coefficient of unmodified gear changes from 0 to 1, and the load distribution coefficient is increased from 0 to 1 after modification, the results show that the reduction of meshing impact, dynamic load, and instantaneous temperature rise is suitable for improving the thermoelastic anti-scuffing load-bearing performance of meshing gears. The load distribution coefficient curves of the spur gear are shown in Figure 9.
Considering the thermoelastic deformation of meshing teeth surface, the driven gear lags behind the driving gear in the gear transmission process to produce a rotation angle error. The fluctuation of the rotation angle error in any meshing cycle restricts the transmission stability of the gear system, that is, the smaller the angle error fluctuation; the more stable the gear system transmission. The variation law of rotation angle error (err/ rad) of spur gear with different base pitch error (err) is revealed as shown in Figure 10.
Vibration acceleration RMS response curves of spur/helical gears along meshing line direction are shown in Figure 11, the operating conditions are set to 1500 Nm and 3000 r/mm. The above figures have been revealed and illustrated, the unmodified gear rotation angle error is mutated twice in one meshing cycle and the load fluctuates greatly. After the modification, the angle error does not change suddenly in one meshing cycle and the load fluctuation is small. When the modification index b = 1.0, the dynamic system error of the double-tooth meshing time domain is obviously  convex, the modification index becomes smaller, the curve convex change is obvious, and the dynamic system error fluctuation range becomes larger. When the modification index b . 1.23, the dynamic system error change of double-tooth meshing time domain is obviously concave. The modification index becomes larger, the concave trend is obvious, and the dynamic system error change range becomes smaller. When the modification index b is 1.23, the angle error fluctuation is the minimum and the transmission is the most stable. Modified vibration acceleration RMS of spur/helical gears are reduced from 232.182/107.737 m/s 2 in the unmodified state to 35.376/16.581 m/s 2 . Therefore, the modified system vibration acceleration RMS is greatly reduced and tends to be constant, which can be effectively improved the gear transmission system vibrates and reduces noise.
Gear machining and installation errors lead to errors in determining the optimal modification amount, which is not the ideal modification amount. The variation law of load distribution coefficient (z) of spur gear with different base pitch error (err) is derived in Figure 12. When the base pitch error err .0, the length of single teeth meshing area increases and the length of double teeth meshing area decreases or even disappears; when the base pitch error (err \ 0), the length of single and double teeth meshing area is constant, but the load fluctuates and changes abruptly.
Considering different base pitch errors (err/rad) of spur gears, as shown in Figure 13, when the base pitch error (err . 0), the angular error (err/rad) is upward convex and wave distribution; when the base pitch error (err \ 0), the angular error (err/rad) has mutation and causes vibration to produce dynamic load. Moreover, the fault-tolerant ability of alternate meshing in realtime all working conditions become weak.

Influence analysis of TSM fitting optimization on temperature field
The TSM fitting optimization will cause redistribution of contact load and friction heat flux load on alternate meshing teeth surface, and then affect the temperature distribution of teeth body and instantaneous contact temperature of teeth surface. The structural parameters of single teeth spur gear are proposed in Table 1, and  Table 2 describes the comprehensive parameters of TECs. Limited to top modification of the pinion and large gear when pinion is meshing, the thermoelastic comprehensive deformation is obtained based on the coupled finite element analysis of the TECs, and then the influence of gear body temperature field on its modification and deformation is analyzed. The lubricating oil temperature is set to 55°C.
For the dimensionless linear coordinate system (considering temperature or not), the position of meshing line is represented, and then the modification curve is determined. The equation of spur gear profile modification curve D = D max x=L ð Þ 1:22 , the parameters x and L in the profile modification equation are respectively dimensionless coordinate values and set L = 0:46.
It can be seen from Table 2 that the elastic deformation of driving gear and driven gear is 2.56 and 8.73 mm at point D of meshing position, and 3.36 and 6.52 mm of driven gear and driven gear respectively at point C of B2 meshing position. Considering the influence of temperature, that is, the teeth deformation is the comprehensive deformation of TECs.
The comprehensive deformation of the driving and driven gears is 12.48 and 12.37 mm at point D of B1 meshing position, and 10.76 and 12.23 mm respectively at point C of B2 meshing position.  The influence of temperature on teeth modification curve cannot be ignored, especially for high-speed and heavy-duty gear modification design. In Figure 14, based on the position of gear meshing starting point B1, the theoretical modification amount of driving gear top is D = 24:58 mm, and the actual modification amount is C a1 = 21:0 mm. The position of gear meshing end point B2, the theoretical modification amount of driving gear root is D = 22:73 mm, and the actual modification amount is C eff = 23:64 mm.

Analysis of coupling influence of TSM fitting on its TECs
The gear teeth optimized for TSM fitting can further improve its thermoelastic coupling characteristics, in order to enhance the comprehensive performance of thermoelastic anti-scuffing load-bearing, in order to improve the GTS operating life and its stability, which can also reduce vibration and noise.

Analysis of coupling influence of spur gear modification fitting on its TECs
Two check conditions of Spur optimal modification (smooth meshing time domain) and improvement of loadbearing performance (anti-thermoelastic scuffing) are set to analyze the effect of TECs coupling, see Table 3.
From the data in the above chart, the maximum value of the unmodified spur gear interface temperature and the teeth surface instantaneous temperature (high thermal T.) along the meshing line in time domain are greater than those of the modified spur gear, but the minimum film thickness is less than the modified. The minimum film thickness (minimum oil film thickness) is the main thermoelastic characteristic quantity for evaluating the anti-thermoelastic scuffing load-bearing performance, see Figure 15 and Table 4 above.

Analysis of coupling influence of helical gear modification fitting on its TECs
This section only considers the helical gear TSM as the optimal check condition. Table 5 indicates the parameters of helical gear pairs. According to the ISO standard, the temperature distribution trend of the unmodified helical gear along the meshing line is analyzed and compared, as shown in Figure 16.
The maximum instantaneous temperature of the unmodified helical gear teeth surface is located at the initial point of meshing time domain. The coordinates of the highest instantaneous temperature of the optimal modified helical gear teeth surface are 20.237 (dimensionless). The parameters related to the highest point of the teeth surface instantaneous temperature in the meshing time domain of unmodified helical gears are shown in Table 6.
From the data in Figure 16 and Table 6, it can be seen that the optimal modified helical gear interface temperature along the meshing line in time domain and the maximum instantaneous teeth surface temperature (high thermal T.) are less than the unmodified values, and the film thickness is the minimum. The value is larger than that of the unmodified helical gear.

Type experiment of gear teeth interface temperature monitoring in meshing time domain
Bench test of meshing gear teeth interface temperature of closed power cycle transmission system The meshing gear teeth interface temperature test bench of the closed power cycle transmission system (CPCTS) is constructed on the main test bench of the power flow cycle closed transmission gear pair and the  Figure 17.
Thermoelectric coupling type temperature acquisition/monitoring and control system Taking into account the transient process of adapting to alternate meshing of gear pairs, the dynamic monitoring of instantaneous temperature of the CPCTS meshing gear teeth interface is based on thermoelectric coupling temperature acquisition and monitoring control system.), thermocouple temperature sensor, collector ring (collecting slip ring), collecting signal, monitoring signal and real-time recording, signal online processing, and display, as shown in Figure 18.   Tables 7 to 9.
In view of the collection of real-time test data and subsequent records and post-analysis and processing through Matlab programming software are shown in Figures 19 to 21.
The results can be obtained by analyzing the test data of the above graphic table.
(1) The surface meshing temperature of the gear teeth increases significantly when the load is high and the rotation speed is high. It can also be understood as: the greater the influence of the load-speed parameter interconnection, the higher the power transmission input causes the gear tooth surface temperature to rise rapidly.

Conclusions
For this pre-research topic, considering the gear comprehensive error, modification status, and accurate teeth surface geometric parameters, based on the  theoretical teeth surface superposition, the TSM fitting is proposed, and the 3D and diagonal modification optimization design is revealed, and the position and normal vectors of the modified surface are derived.
Based on LTCA, a variety of optimization modification models are established, and the complex surface is analyzed to fit the thermoelastic contact numerical simulation of actual teeth surface, it is expected to reduce the transmission error and vibration noise of GTS in order to improve the meshing gear teeth loadbearing performance. The results show that, compared with the previous analysis methods, the proposed analysis method is effective and simple.
(1) The results show that there is no edge contact after three-dimensional topology modification, but the teeth tip edge contact cannot be avoided after diagonal modification. The multi load transmission error and load transmission error amplitude are significantly reduced after fitting optimization. (2) The modified teeth profile reduces the impact of meshing, the optimally fitted meshing-in and meshing-out ends do not bear load, and the alternate meshing area of the modified single tooth is increased. (3) In view of the three-dimensional topology modification, the teeth surface lubrication state is improved, and the accumulated amplitude of the instantaneous teeth surface temperature rise at the meshing-in and meshing-out ends is obviously reduced, but the accumulated amplitude of the instantaneous temperature rise of teeth surface near the node element is not easy to fall. (4) After modification, the load transmission error and meshing impact are reduced, and the root mean square of acceleration is also reduced. It shows that the multi-objective optimization analysis of TSM can significantly reduce the transmission error, meshing impact and vibration acceleration RMS, and further improve     the alternate meshing performance of thermoelastic properties of gear teeth interface. (5) Based on the width/narrow and long and small helical teeth tip modification value is not greater than its optimal modification value and the large helical teeth tip modification value is greater than its optimal modification value, it can cause a large helical teeth tip circle radius to reduce and the actual length of the meshing line to decrease. The near area of point A presents a point that does not mesh alternately, and the near area gear temperature of point A is regarded as wide/narrow and large/small helical gear body temperature (teeth temperature). (6) The results show that the length of the addendum is larger than the optimal value of the addendum, when the top modification value of small helical gear is not less than its optimal value and the modification value of large and small helical gear is not greater than its optimal  value, the radius of small helical gear top circle can be reduced and the actual length of meshing line can be shortened. The near field of point E also shows a point which does not participate in alternate meshing, and the temperature of near field gear at point E is regarded as the interface temperature of the gear teeth body.