Friction Coefficient Compensation Control in Synchronizer Synchronization Process for Transmission

In the synchronization process of a synchronizer, the friction coefficient of the friction cone is changed because of friction heat. In this paper, the influence of a changing friction coefficient on shift quality was analyzed. The quantitative relationship between friction coefficient and its influencing factors was studied through a synchronizer synchronization process test bench. Based on the quantitative analysis and an optimal shift force control method, a compensation control strategy for friction coefficient was established. Moreover, the effectiveness of the compensation control was verified through simulation and experiment. The results showed that the friction coefficient was maintained near the expected value of 0.08 after the compensation control, the shifting speed difference the synchronous time was shortened by nearly 0.12 s, and the sliding friction was reduced by 4.64 J under the experimental conditions. The analysis and compensation of the friction coefficient provide a theoretical reference for improving shift quality.


Introduction
As one of the main components of a transmission, the synchronizer generates friction torque during friction synchronization, thereby eliminating the difference in rotational speed between the two ends of the synchronizer, to achieve shifting [1]. During the synchronization of the synchronizer, the friction coefficient of the friction cone is changed because of friction heat, which ultimately affects the shift control accuracy of the transmission [2].
A lot of basic research has been proposed to improve shift control accuracy [3,4]. The temperature changing regularities of the friction surface under different working conditions were analyzed by simulation and experiment [5,6]. The influence of temperature on the synchronizer's working performance was verified in literature [7]. The friction and heat generation process of the synchronizer with different friction coefficients were analyzed in the literature [8,9]. The results showed that the greater the friction coefficient is, within a certain range, the higher the temperature between the friction surfaces. Additionally, Li presented a model of predicting wear in the friction lining of a wet clutch and introduced that the wear mechanism is related to thermal degradation [10]. Marklund developed a friction model which takes temperature, speed, and nominal pressure into account for the wet clutch during boundary lubrication conditions [11]. The friction surface temperature and the friction coefficient of the friction cone are the critical influences on the shift quality, which need further study.
Currently, research on car synchronizers has mainly focused on the influence of the geometry of the synchronizer and the synchronous inertia on the synchronous performance, and considered the friction coefficient between the friction pairs of the synchronizer as a constant. In reference [12], the ADAMS software was used to establish the virtual prototype model of the lock-ring synchronizer, and the influence of each parameter on the shift synchronization time was studied. Huang Hua regarded the friction coefficient as a constant in the process of establishing the synchronizer model [13]. In the synchronous shift phase, the synchronizer realizes the synchronization function by sliding between the friction cones. The length of the synchronization time is directly determined by the friction coefficient between the friction cones, so the calculation method of considering the friction coefficient as a constant is not accurate. Besides this, with the deepening of the research on the synchronizer by domestic and foreign research institutions, research on the friction coefficient between the friction pairs of the synchronizer has deepened. In reference [14], the dynamic model of the double-cone synchronizer was established based on the ADAMS software. On this basis, the effects of spring stiffness, the cone angle of friction ring, and the friction coefficient on shift force and the synchronization process were analyzed. Daniel Häggström analyzed the influence of temperature, speed difference, and shift force on the friction coefficient during synchronizer synchronization [15]. Sun Jun compared the influence of the friction coefficient on synchronizer performance and structure parameters under different inertia and speed difference [16]. Xi Junqiang conducted a bench test on the working process of the synchronizer. The results showed that the friction coefficient was not maintained at a constant value. In order to accurately reflect the real situation of the synchronization process, the friction coefficient must be considered [17].
In this paper, we study the influencing factors of the friction coefficient change of the synchronizer, and explore the quantitative variation law between the friction coefficient and each influencing factor, then develop the friction coefficient compensation control strategy to maintain the friction coefficient near the expected value and shorten the synchronization time. This is important to improve the shift quality.

Evaluation Index of Synchronization Performance
The shifting process should mainly follow the following principles: shorten shift synchronization time and drive interruption time as much as possible; the shifting process should avoid collisions and impacts to achieve smooth shifting; and minimize the sliding-friction work of the synchronization process to extend the synchronizer's service life. Therefore, the shift quality can be evaluated in terms of power interruption time, smoothness, and service life [18]. The locking ring synchronizer of automated mechanical transmission (a two-degrees-of-freedom direct-drive electromagnetic gearshift system) was the research object in this paper, as shown in Figure 1. The schematic diagram of the locking ring synchronizer is shown in Figure 1a, which consists of a synchromesh sliding sleeve, a splined hub, a synchronous ring, a target ring gear, and a locating slider [4,19].
Currently, research on car synchronizers has mainly focused on the influence of the geometry of the synchronizer and the synchronous inertia on the synchronous performance, and considered the friction coefficient between the friction pairs of the synchronizer as a constant. In reference [12], the ADAMS software was used to establish the virtual prototype model of the lock-ring synchronizer, and the influence of each parameter on the shift synchronization time was studied. Huang Hua regarded the friction coefficient as a constant in the process of establishing the synchronizer model [13]. In the synchronous shift phase, the synchronizer realizes the synchronization function by sliding between the friction cones. The length of the synchronization time is directly determined by the friction coefficient between the friction cones, so the calculation method of considering the friction coefficient as a constant is not accurate. Besides this, with the deepening of the research on the synchronizer by domestic and foreign research institutions, research on the friction coefficient between the friction pairs of the synchronizer has deepened. In reference [14], the dynamic model of the double-cone synchronizer was established based on the ADAMS software. On this basis, the effects of spring stiffness, the cone angle of friction ring, and the friction coefficient on shift force and the synchronization process were analyzed. Daniel Häggström analyzed the influence of temperature, speed difference, and shift force on the friction coefficient during synchronizer synchronization [15]. Sun Jun compared the influence of the friction coefficient on synchronizer performance and structure parameters under different inertia and speed difference [16]. Xi Junqiang conducted a bench test on the working process of the synchronizer. The results showed that the friction coefficient was not maintained at a constant value. In order to accurately reflect the real situation of the synchronization process, the friction coefficient must be considered [17].
In this paper, we study the influencing factors of the friction coefficient change of the synchronizer, and explore the quantitative variation law between the friction coefficient and each influencing factor, then develop the friction coefficient compensation control strategy to maintain the friction coefficient near the expected value and shorten the synchronization time. This is important to improve the shift quality.

Evaluation Index of Synchronization Performance
The shifting process should mainly follow the following principles: shorten shift synchronization time and drive interruption time as much as possible; the shifting process should avoid collisions and impacts to achieve smooth shifting; and minimize the sliding-friction work of the synchronization process to extend the synchronizer's service life. Therefore, the shift quality can be evaluated in terms of power interruption time, smoothness, and service life [18]. The locking ring synchronizer of automated mechanical transmission (a two-degrees-of-freedom direct-drive electromagnetic gearshift system) was the research object in this paper, as shown in Figure 1. The schematic diagram of the locking ring synchronizer is shown in Figure 1a, which consists of a synchromesh sliding sleeve, a splined hub, a synchronous ring, a target ring gear, and a locating slider [4,19].

Shift Time
The total time t used for shifting is composed of the time t 1 taken for picking up the gear, the time t 2 used for pre-synchronization, the time t s used for synchronization, and the time t 4 taken for the combination-of these, the proportion of synchronization time t s is larger. The phase represented by the synchronization time has the greatest influence on the sliding-friction work and shock degree [20]. The following mainly studies the synchronization time t s .
Establish the motion equation of the inertial mass at the input end as: The synchronization time: (3) where T f is the friction torque, N·m; µ is the friction coefficient; F s is the axial thrust, N; R is the average radius of the friction cone, m; J r is the rotational inertia of the synchronizer's input, Kg·m 2 ; ω r is the angular speed of the synchronizer's input, r/min; ω c is the angular speed of the synchronizer's output, r/min; ω e is the speed of the motor's input shaft, which is also the initial angular speed of synchronizer's input, r/min; α is half the cone angle of the friction cone, rad.

Shock Degree During Synchronous Phase
The shock degree (j) during the synchronous phase refers to the change rate of the angular speed of the synchronizer's input under the action of synchronous torque. In the synchronization phase, the smoothness of the system is closely related to the shock degree. The greater the shock degree, the more serious the vibration and noise will be. The shock degree during the synchronous phase is an important index to measure shift performance. Its calculation formula is, where T c is the resistance moment of the transmission, j max is the maximum shock degree during synchronous phase. It can be seen from Equation (5) that the change rate of synchronous torque determines the magnitude of the shock degree during synchronous phase. During the shift synchronization, the corresponding control strategy should be added to reduce the shock degree of this process, so as to improve the shift quality.

The Sliding-Friction Work
The synchronizer synchronizes between the synchronizing ring and the target ring gear, so that the components with different rotational speeds are synchronized, and friction heat generation is the most important implementation. The sliding-friction work is the work done by the friction torque on the relative angular displacement between the friction elements during the entire synchronization process. If excessive heat generates in the unit area during the sliding in a short time, it will cause irreversible damage to friction components and reduce the service life of synchronizer. The formula for calculating sliding-friction work (W) is as follows: The evaluation indexes restrict each other during the shift synchronization process. In the process of pursuing an evaluation index improvement, it is necessary to coordinate with other indicators to ensure that the shift quality is comprehensively optimized.

Simulation Analysis of Synchronization Performance
The simulation model of a synchronizer in synchronization phase was built, and the influence of friction coefficient change between the synchronizer's friction pairs on synchronization performance was analyzed. According to the value range of the friction coefficient of the locking ring synchronizer's friction pair selected in this paper, the friction coefficient was set to 0.05, 0.07, and 0.09 in the simulation process, and the rotational inertia of the synchronized partial was 0.04 Kg·m 2 . The influences of speed difference and friction coefficient on synchronization time were analyzed, as shown in Figure 2.
The evaluation indexes restrict each other during the shift synchronization process. In the process of pursuing an evaluation index improvement, it is necessary to coordinate with other indicators to ensure that the shift quality is comprehensively optimized.

Simulation Analysis of Synchronization Performance
The simulation model of a synchronizer in synchronization phase was built, and the influence of friction coefficient change between the synchronizer's friction pairs on synchronization performance was analyzed. According to the value range of the friction coefficient of the locking ring synchronizer's friction pair selected in this paper, the friction coefficient was set to 0.05, 0.07, and 0.09 in the simulation process, and the rotational inertia of the synchronized partial was 0.04 Kg·m 2 . The influences of speed difference and friction coefficient on synchronization time were analyzed, as shown in Figure 2. The larger the friction coefficient was, the shorter the synchronization time that was needed when the initial shifting speed difference was the same. With the increase in shifting speed difference, the effect of friction coefficient on synchronization time became more and more obvious. When the shifting speed difference was 600 r/min and the friction coefficient was 0.09, the synchronization time was 0.58 s. When the friction coefficient was 0.05, the synchronization time was 0.79 s. The difference between them was 0.21 s. According to the simulation, the friction coefficient should be increased in order to shorten the synchronization time.
When the inertia of the synchronized part was 0.04 Kg·m 2 , the influence of speed difference and friction coefficient on sliding-friction work was analyzed, as shown in Figure 3. The variation trend of sliding-friction work was basically the same under different shifting speed differences and friction coefficients. When the shifting speed difference was the same, the friction coefficient had no effect on the final value of sliding-friction work, but only affected the increase rate of sliding-friction work. The greater the friction coefficient was, the shorter the time it took for the sliding-friction work to reach the final value. When the shifting speed difference was 600 r/min, the sliding-friction work consumed under the three friction coefficients was 78.96 J. However, when the friction coefficient was 0.09, the time needed to reach the final value of sliding-friction work was 0.21 s shorter than when the friction coefficient was 0.05. When the final value of sliding-friction work was determined, the shorter the synchronization time was, the more heat generated per unit of time. Moreover, excessive heat generated per unit of time may have caused damage to the friction surface.
When the inertia of the synchronized part was 0.04 Kg·m 2 , the influence of speed difference and friction coefficient on shock degree were analyzed, as shown in Figure 4. Since the maximum shock degree occurred at the initial stage in the synchronization process, the shock degree at the initial stage was intercepted for comparative analysis. It can be seen from the graph that the maximum shock The larger the friction coefficient was, the shorter the synchronization time that was needed when the initial shifting speed difference was the same. With the increase in shifting speed difference, the effect of friction coefficient on synchronization time became more and more obvious. When the shifting speed difference was 600 r/min and the friction coefficient was 0.09, the synchronization time was 0.58 s. When the friction coefficient was 0.05, the synchronization time was 0.79 s. The difference between them was 0.21 s. According to the simulation, the friction coefficient should be increased in order to shorten the synchronization time.
When the inertia of the synchronized part was 0.04 Kg·m 2 , the influence of speed difference and friction coefficient on sliding-friction work was analyzed, as shown in Figure 3. The variation trend of sliding-friction work was basically the same under different shifting speed differences and friction coefficients. When the shifting speed difference was the same, the friction coefficient had no effect on the final value of sliding-friction work, but only affected the increase rate of sliding-friction work. The greater the friction coefficient was, the shorter the time it took for the sliding-friction work to reach the final value. When the shifting speed difference was 600 r/min, the sliding-friction work consumed under the three friction coefficients was 78.96 J. However, when the friction coefficient was 0.09, the time needed to reach the final value of sliding-friction work was 0.21 s shorter than when the friction coefficient was 0.05. When the final value of sliding-friction work was determined, the shorter the synchronization time was, the more heat generated per unit of time. Moreover, excessive heat generated per unit of time may have caused damage to the friction surface.
When the inertia of the synchronized part was 0.04 Kg·m 2 , the influence of speed difference and friction coefficient on shock degree were analyzed, as shown in Figure 4. Since the maximum shock degree occurred at the initial stage in the synchronization process, the shock degree at the initial stage was intercepted for comparative analysis. It can be seen from the graph that the maximum shock degree decreased with the increase of friction coefficient when the shifting speed difference was the same. When the shifting speed difference was 600 r/min, the maximum shock degree was 160.3 rad/s 3 when the friction coefficient was 0.09, and 290.4 rad/s 3 when the friction coefficient was 0.05.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 5 of 15 degree decreased with the increase of friction coefficient when the shifting speed difference was the same. When the shifting speed difference was 600 r/min, the maximum shock degree was 160.3 rad/s 3 when the friction coefficient was 0.09, and 290.4 rad/s 3 when the friction coefficient was 0.05.  degree decreased with the increase of friction coefficient when the shifting speed difference was the same. When the shifting speed difference was 600 r/min, the maximum shock degree was 160.3 rad/s 3 when the friction coefficient was 0.09, and 290.4 rad/s 3 when the friction coefficient was 0.05.  In conclusion, with the increase in the friction coefficient, the synchronization time of the synchronizer decreased. The friction coefficient had no effect on the final value of sliding-friction work in the synchronous process, but it affected the increase rate of the sliding-friction work. At the same time, with the increase in friction coefficient, it took less time for the sliding-friction work to reach the final value, and the maximum shock degree at the initial stage decreased.

Analysis of Friction Coefficient's Influencing Factors
The premise of formulating the friction coefficient compensation control strategy is to obtain the quantitative relationship between the friction coefficient and its influencing factors. Thus, the friction test of the synchronization process must be performed first. A synchronizer test-control platform was built, and the friction coefficient MAP was obtained by collecting and calculating data during the test.

The Synchronizer Test-Control Platform
Based on the study of the synchronization process and mechanical analysis, combined with a modern electromechanical design concept, the synchronizer test-control platform was completed. The designed bench had the functions of variable voltage, variable speed, and constant output. The physical object of the synchronizer test-control platform is shown in Figure 5 [19,21]. The mechanical system of the test-control platform included inertia simulator, synchronizer component, shifting fork, displacement sensor (DS), and the connector that connected the components. The electronic control system included a driving device control system and three-phase asynchronous motor. The electromagnetic linear actuator (ELA) was controlled and driven by a driving device control system. The three-phase asynchronous motor output different speeds according to the driving requirements for the input shaft [22,23]. The inertia simulator on the test-control platform could simulate the range of rotational inertia from 0.04 to 0.09 kg·m 2 . The three-phase asynchronous motor was controlled by a frequency converter, and its maximum driving speed could reach 1200 r/min. The electromagnetic linear actuator was used as a shift driving device, the maximum shift force was 1000 N. The temperature was tested through an infrared camera (Optris PI450) with −20~900 • C measuring range.  In conclusion, with the increase in the friction coefficient, the synchronization time of the synchronizer decreased. The friction coefficient had no effect on the final value of sliding-friction work in the synchronous process, but it affected the increase rate of the sliding-friction work. At the same time, with the increase in friction coefficient, it took less time for the sliding-friction work to reach the final value, and the maximum shock degree at the initial stage decreased.

Analysis of Friction Coefficient's Influencing Factors
The premise of formulating the friction coefficient compensation control strategy is to obtain the quantitative relationship between the friction coefficient and its influencing factors. Thus, the friction test of the synchronization process must be performed first. A synchronizer test-control platform was built, and the friction coefficient MAP was obtained by collecting and calculating data during the test.

The Synchronizer Test-Control Platform
Based on the study of the synchronization process and mechanical analysis, combined with a modern electromechanical design concept, the synchronizer test-control platform was completed. The designed bench had the functions of variable voltage, variable speed, and constant output. The physical object of the synchronizer test-control platform is shown in Figure 5 [19,21]. The mechanical system of the test-control platform included inertia simulator, synchronizer component, shifting fork, displacement sensor (DS), and the connector that connected the components. The electronic control system included a driving device control system and three-phase asynchronous motor. The electromagnetic linear actuator (ELA) was controlled and driven by a driving device control system. The three-phase asynchronous motor output different speeds according to the driving requirements for the input shaft [22,23]. The inertia simulator on the test-control platform could simulate the range of rotational inertia from 0.04 to 0.09 kg·m 2 . The three-phase asynchronous motor was controlled by a frequency converter, and its maximum driving speed could reach 1200 r/min. The electromagnetic linear actuator was used as a shift driving device, the maximum shift force was 1000 N. The temperature was tested through an infrared camera (Optris PI450) with −20~900 °C measuring range.

Quantitative Relationship between the Friction Coefficient and its Influencing Factors
The results show that the friction coefficient between the friction cones of synchronizer is related to temperature, load, and speed difference between the synchronizer ring and target gear ring after the material and lubrication conditions of the synchronizer's friction cones are determined [15,[24][25][26][27]. The premise of formulating the friction coefficient compensation control strategy was to obtain the quantitative relationship between friction coefficient and its influencing factors in the synchronization process. Firstly, the friction test of the synchronization process was carried out to obtain the sampling data of relevant parameters under different working conditions. The collected data was calculated and processed to obtain the control MAP, which was stored in the controller.
When the speed difference and shift force of the synchronizer were constant, the friction test between the synchronization ring and target gear ring was carried out. The controller collected state variables. In the same step, friction tests were carried out under different working conditions. The data was analyzed as typical when the synchronization time was 0.5 s under different working conditions. The operating conditions of the actual synchronization process were arbitrary, but the operating points of the MAP obtained by collecting data were limited and discrete, so the bicubic interpolation optimization method was needed to obtain the parameter values between the known operating conditions and the temperature and friction coefficient in the whole working condition could be obtained by the parameter values of finite working conditions [28]. After interpolation optimization, the friction coefficient MAP of 0.5 s synchronization time under different working conditions was obtained, as shown in Figure 6.
The results show that the friction coefficient between the friction cones of synchronizer is related to temperature, load, and speed difference between the synchronizer ring and target gear ring after the material and lubrication conditions of the synchronizer's friction cones are determined [15,[24][25][26][27]. The premise of formulating the friction coefficient compensation control strategy was to obtain the quantitative relationship between friction coefficient and its influencing factors in the synchronization process. Firstly, the friction test of the synchronization process was carried out to obtain the sampling data of relevant parameters under different working conditions. The collected data was calculated and processed to obtain the control MAP, which was stored in the controller.
When the speed difference and shift force of the synchronizer were constant, the friction test between the synchronization ring and target gear ring was carried out. The controller collected state variables. In the same step, friction tests were carried out under different working conditions. The data was analyzed as typical when the synchronization time was 0.5 s under different working conditions. The operating conditions of the actual synchronization process were arbitrary, but the operating points of the MAP obtained by collecting data were limited and discrete, so the bicubic interpolation optimization method was needed to obtain the parameter values between the known operating conditions and the temperature and friction coefficient in the whole working condition could be obtained by the parameter values of finite working conditions [28]. After interpolation optimization, the friction coefficient MAP of 0.5 s synchronization time under different working conditions was obtained, as shown in Figure 6. When the shifting speed difference was 600 r/min and the shifting force was 600 N, the corresponding relationship of the friction coefficient changes with temperature was obtained, as shown in Figure 7. It can be seen from the graph that the friction coefficient decreased with the increase in temperature, and the friction coefficient decreased slowly when the temperature was low, and rapidly when the temperature was high. Due to the temperature rise, the specific surface energy and lattice resistance of the friction pair friction material changed. The temperature rise had an important influence on the mechanical properties and surface film formation of the friction material's surface layer, which made it easy to produce liquid lubrication film and resulted in the reduction in the friction coefficient.
Above is the corresponding relationship between the friction coefficient and speed difference, shift force and temperature when synchronization time was 0.5 s. In order to obtain a quantitative relationship between the friction coefficient and the influencing factors in the whole working condition, a large number of tests were still needed, which will not be repeated here. The above work laid a foundation for formulating the friction coefficient compensation control strategy. When the shifting speed difference was 600 r/min and the shifting force was 600 N, the corresponding relationship of the friction coefficient changes with temperature was obtained, as shown in Figure 7. It can be seen from the graph that the friction coefficient decreased with the increase in temperature, and the friction coefficient decreased slowly when the temperature was low, and rapidly when the temperature was high. Due to the temperature rise, the specific surface energy and lattice resistance of the friction pair friction material changed. The temperature rise had an important influence on the mechanical properties and surface film formation of the friction material's surface layer, which made it easy to produce liquid lubrication film and resulted in the reduction in the friction coefficient.
Above is the corresponding relationship between the friction coefficient and speed difference, shift force and temperature when synchronization time was 0.5 s. In order to obtain a quantitative relationship between the friction coefficient and the influencing factors in the whole working condition, a large number of tests were still needed, which will not be repeated here. The above work laid a foundation for formulating the friction coefficient compensation control strategy.

Optimal Control
In the synchronization process, the synchronization ring and the target gear ring are compressed by the shifting force, and the speed difference between them is eliminated by friction torque. The synchronization process takes a certain amount of time, the sliding-friction work and synchronization impact are also produced. The length of the synchronization time affects the shift interruption time. The magnitude of sliding-friction work affects the service life of the synchronizer, and the magnitude of shock degree affects the smoothness of synchronization process. The friction coefficient compensation control strategy of the synchronization process studied in this paper was bound to cause an increase in the shifting force. Increasing the shifting force would shorten the synchronization time, but due to the limitation of friction pair material and other two evaluation indexes, a greater shifting force would not be better. When the initial speed difference was fixed, shortening the synchronization time meant that the speed difference decreased faster and the shock degree was higher. However, in order to reduce the shock degree and reduce the shifting force, the synchronization time was prolonged and the synchronizer wear was increased. Therefore, in order to balance the evaluation indexes of each synchronization performance, it was necessary to limit the magnitude of the shifting force acquired by the friction coefficient compensation control strategy. Other evaluation indicators of synchronization performance were also optimized, while reducing synchronization time. According to the principle of minimum value, the optimal control method was adopted, and the evaluation indexes of the synchronization process of the synchronizer were optimized [29,30].
The mathematical model of the shift synchronization process is expressed as: where Jn = [Jc + Jm(ig, in)], Jm is the rotational inertia of the rotating part of the transmission equivalent to the input shaft, Jc is the rotational inertia of clutch driven disk, ig is transmission ratio, in is the final drive ratio, Tm is the external resistance moment of the vehicle. The state vector is taken as:

Optimal Control
In the synchronization process, the synchronization ring and the target gear ring are compressed by the shifting force, and the speed difference between them is eliminated by friction torque. The synchronization process takes a certain amount of time, the sliding-friction work and synchronization impact are also produced. The length of the synchronization time affects the shift interruption time. The magnitude of sliding-friction work affects the service life of the synchronizer, and the magnitude of shock degree affects the smoothness of synchronization process. The friction coefficient compensation control strategy of the synchronization process studied in this paper was bound to cause an increase in the shifting force. Increasing the shifting force would shorten the synchronization time, but due to the limitation of friction pair material and other two evaluation indexes, a greater shifting force would not be better. When the initial speed difference was fixed, shortening the synchronization time meant that the speed difference decreased faster and the shock degree was higher. However, in order to reduce the shock degree and reduce the shifting force, the synchronization time was prolonged and the synchronizer wear was increased. Therefore, in order to balance the evaluation indexes of each synchronization performance, it was necessary to limit the magnitude of the shifting force acquired by the friction coefficient compensation control strategy. Other evaluation indicators of synchronization performance were also optimized, while reducing synchronization time. According to the principle of minimum value, the optimal control method was adopted, and the evaluation indexes of the synchronization process of the synchronizer were optimized [29,30].
The mathematical model of the shift synchronization process is expressed as: where J n = [J c + J m (i g , i n )], J m is the rotational inertia of the rotating part of the transmission equivalent to the input shaft, J c is the rotational inertia of clutch driven disk, i g is transmission ratio, i n is the final drive ratio, T m is the external resistance moment of the vehicle. The state vector is taken as: The control vector is: The state equation is established as: The maximum driving force produced by the electromagnetic linear actuator determined the range of synchronous moment. The maximum driving force of the actuator used in this paper was 1000 N, so the maximum synchronous torque was 13.2 N·m, that 0 ≤ y 1 ≤ [y 1 ] = 13.2 N·m. y 2 stands for the change rate of synchronous torque, while the magnitude of the shock degree in the synchronous phase is determined by the change rate of synchronous torque. Therefore, the shock degree can be well limited by limiting the change rate of synchronous torque in synchronization phase. The constraint condition of calculating y 2 is 0 ≤ y 2 ≤ [y 2 ] = 1596. In order to optimize each evaluation index, the optimal objective function is established as: where J is the scalar function, Q 1 and Q 2 are the weight coefficients of sliding-friction work and shock degree, respectively. Due to the short shifting time, it is generally assumed that the external resistance torque of the vehicle is unchanged during this process, .
T m = 0. Let k = r/i n , the Hamiltonian function introduced based on the minimum principle is: where λ 1 , λ 2 , and λ 3 are co-state variables. It is generally believed that the shifting force increases monotonously during the synchronization process, and the synchronous torque increases monotonously. Therefore, as long as the boundary of the optimal control variable is limited, the relative optimal trajectory and the optimal control result in the control domain range can be obtained. Computation is based on the necessary conditions to be met by minimum principles, and the intermediate steps are omitted. Thus, Finally, the optimal shifting force trajectory is obtained by the optimal control method:

Compensation Control Principle
The optimal shifting force trajectory F* was obtained by the optimal control method. Through the force controller, the voltage input to the electromagnetic linear actuator at different synchronization times was obtained, and the actuator output different shifting forces to the synchronizer according to different voltage values. During the synchronization process of the synchronizer, due to the difference in the rotation speed between the synchronizing ring and the target ring gear, the shifting force acted on the friction pair to generate the friction torque. The friction torque caused the temperature to rise during the synchronous speed difference, and the temperature accumulation caused the friction coefficient between the friction pairs to decrease. The real-time shifting force, the speed difference, and the temperature during the synchronization process were transmitted to the controller through the sensor [31][32][33]. The controller obtained the real-time friction coefficient µ by querying the relevant friction coefficient MAP. The real-time friction coefficient was compared with the expected friction coefficient µ* to obtain the difference. The compensation controller calculated the compensation shifting force F that needed to be output according to the difference. The compensation control force corrected the optimal control force in real time, and changed the force signal input to the force converter to maintain the friction coefficient of the friction pair of the synchronizer near the expected value. The feedback compensation control schematic is shown in Figure 8. Finally, the optimal shifting force trajectory is obtained by the optimal control method:

Compensation Control Principle
The optimal shifting force trajectory F* was obtained by the optimal control method. Through the force controller, the voltage input to the electromagnetic linear actuator at different synchronization times was obtained, and the actuator output different shifting forces to the synchronizer according to different voltage values. During the synchronization process of the synchronizer, due to the difference in the rotation speed between the synchronizing ring and the target ring gear, the shifting force acted on the friction pair to generate the friction torque. The friction torque caused the temperature to rise during the synchronous speed difference, and the temperature accumulation caused the friction coefficient between the friction pairs to decrease. The real-time shifting force, the speed difference, and the temperature during the synchronization process were transmitted to the controller through the sensor [31][32][33]. The controller obtained the real-time friction coefficient μ by querying the relevant friction coefficient MAP. The real-time friction coefficient was compared with the expected friction coefficient μ* to obtain the difference. The compensation controller calculated the compensation shifting force F that needed to be output according to the difference. The compensation control force corrected the optimal control force in real time, and changed the force signal input to the force converter to maintain the friction coefficient of the friction pair of the synchronizer near the expected value. The feedback compensation control schematic is shown in Figure 8.

Results and Analysis
Based on the friction coefficient compensation control principle and the simulation model of the friction coefficient compensation control in the synchronizer synchronization process, the compensation control process was simulated and analyzed, and the synchronization performance index change after friction coefficient compensation was obtained. Because the friction coefficient compensation control system in the synchronization process of the synchronizer was a complex system with strong time-varying and multi-dimensional nonlinearity, the experiment was also carried out to verify the feasibility of the proposed friction coefficient compensation control strategy.

Synchronization Time
According to the literature and the performance of the synchronization loop, the expected friction coefficient of the synchronization process was set to be 0.08, and the moment of inertia of the synchronized part was 0.04 Kg·m 2 . The compensation control simulation and test were carried out

Results and Analysis
Based on the friction coefficient compensation control principle and the simulation model of the friction coefficient compensation control in the synchronizer synchronization process, the compensation control process was simulated and analyzed, and the synchronization performance index change after friction coefficient compensation was obtained. Because the friction coefficient compensation control system in the synchronization process of the synchronizer was a complex system with strong time-varying and multi-dimensional nonlinearity, the experiment was also carried out to verify the feasibility of the proposed friction coefficient compensation control strategy.

Synchronization Time
According to the literature and the performance of the synchronization loop, the expected friction coefficient of the synchronization process was set to be 0.08, and the moment of inertia of the synchronized part was 0.04 Kg·m 2 . The compensation control simulation and test were carried out under different initial shift speeds of the shift. The results of the simulation and test synchronization time with or without compensation were obtained, as shown in Figure 9. When the initial rotational speed difference of the shift was the same, the synchronization time used before the compensation in the test was longer than the compensation, and the regular was similar under different rotational speed differences. The shift speed difference was 600 r/min, the synchronization time used in the pre-compensation test was 0.86 s, and the synchronization time used after compensation was 0.74 s. The compensation time required after compensation was shortened by 0.12 s. The simulation results match the test results. The shift speed difference was 600 r/min. The synchronization time used for the pre-compensation simulation was 0.70 s. The synchronization time used after compensation was 0.62 s. The compensation time required after compensation was shortened by 0.08 s. It can be concluded that the synchronization time can be significantly shortened by compensating the friction coefficient, to ensure its maintenance near the expected value. time with or without compensation were obtained, as shown in Figure 9. When the initial rotational speed difference of the shift was the same, the synchronization time used before the compensation in the test was longer than the compensation, and the regular was similar under different rotational speed differences. The shift speed difference was 600 r/min, the synchronization time used in the precompensation test was 0.86 s, and the synchronization time used after compensation was 0.74 s. The compensation time required after compensation was shortened by 0.12 s. The simulation results match the test results. The shift speed difference was 600 r/min. The synchronization time used for the pre-compensation simulation was 0.70 s. The synchronization time used after compensation was 0.62 s. The compensation time required after compensation was shortened by 0.08 s. It can be concluded that the synchronization time can be significantly shortened by compensating the friction coefficient, to ensure its maintenance near the expected value.

Shock Degree
The initial speed difference was 600 r/min, and the moment of inertia of the synchronized part was 0.04 Kg·m 2 . The bench test of the synchronous process was carried out, and the comparison of the shock degree test before and after compensation was obtained, as shown in Figure 10. The initial synchronization phase had a relatively large shock degree before the compensation control, but the shock degree during the synchronization process was small and relatively stable. After the compensation control, the shock degree of the whole synchronization phase was relatively large, but the shock degree decreased with the continuous synchronization process. It can be seen from the synchronous phase shock degree calculation formula that the synchronous shock degree was related to the synchronous torque change rate during the synchronization of the synchronizer. The synchronous torque at this stage was the friction torque generated between the friction pairs of the synchronizer, and the magnitude of the friction torque was closely related to the shifting force acting on the shift fork. Therefore, the magnitude of the synchronous shock degree depended on the rate of change of the shifting force acting on the shift fork. It can be seen that in the process of compensating the friction coefficient of the synchronous process, the rate of change of the shifting force was gradually reduced, and the compensation control caused the shock degree in the synchronization process to become large.

Shock Degree
The initial speed difference was 600 r/min, and the moment of inertia of the synchronized part was 0.04 Kg·m 2 . The bench test of the synchronous process was carried out, and the comparison of the shock degree test before and after compensation was obtained, as shown in Figure 10. The initial synchronization phase had a relatively large shock degree before the compensation control, but the shock degree during the synchronization process was small and relatively stable. After the compensation control, the shock degree of the whole synchronization phase was relatively large, but the shock degree decreased with the continuous synchronization process. It can be seen from the synchronous phase shock degree calculation formula that the synchronous shock degree was related to the synchronous torque change rate during the synchronization of the synchronizer. The synchronous torque at this stage was the friction torque generated between the friction pairs of the synchronizer, and the magnitude of the friction torque was closely related to the shifting force acting on the shift fork. Therefore, the magnitude of the synchronous shock degree depended on the rate of change of the shifting force acting on the shift fork. It can be seen that in the process of compensating the friction coefficient of the synchronous process, the rate of change of the shifting force was gradually reduced, and the compensation control caused the shock degree in the synchronization process to become large.

Sliding-Friction Work
When the initial speed difference of the shift was 600 r/min, the expected friction coefficient of the synchronization process was 0.08. The simulation results of sliding-friction work before and after compensation were compared, as shown in Figure 11.

Sliding-Friction Work
When the initial speed difference of the shift was 600 r/min, the expected friction coefficient of the synchronization process was 0.08. The simulation results of sliding-friction work before and after compensation were compared, as shown in Figure 11. Figure 10. Comparison of shock degree before and after compensation.

Sliding-Friction Work
When the initial speed difference of the shift was 600 r/min, the expected friction coefficient of the synchronization process was 0.08. The simulation results of sliding-friction work before and after compensation were compared, as shown in Figure 11. The final value sliding-friction work after the synchronization phase of the synchronizer was less than that before the compensation, and the time to reach the final value after compensation was shorter than that before compensation. It can be seen from the trend of the sliding-friction work that after the compensation, the sliding-friction work growth in the early stage was slow, and the growth in the later stage was fast, while that before the compensation sliding-friction work had a fast growth in the early stage and a slow growth in the later stage. The shifting speed difference was 600 r/min, the final value sliding-friction work in the before compensation phase was 77.89 J, and the final value sliding-friction work in the after compensation was 73.25 J. After the compensation, the final slidingfriction work value was reduced by 4.64 J. It can be seen from the comparison of the sliding-friction work before and after compensation that the sliding-friction work can be reduced by compensating the friction coefficient.

Friction Coefficient
When the initial shift speed difference was 600 r/min, the moment of inertia of the synchronized part was 0.04 Kg·m 2 . After compensation, the test and simulated friction coefficient results were compared, as shown in Figure 12. The final value sliding-friction work after the synchronization phase of the synchronizer was less than that before the compensation, and the time to reach the final value after compensation was shorter than that before compensation. It can be seen from the trend of the sliding-friction work that after the compensation, the sliding-friction work growth in the early stage was slow, and the growth in the later stage was fast, while that before the compensation sliding-friction work had a fast growth in the early stage and a slow growth in the later stage. The shifting speed difference was 600 r/min, the final value sliding-friction work in the before compensation phase was 77.89 J, and the final value sliding-friction work in the after compensation was 73.25 J. After the compensation, the final sliding-friction work value was reduced by 4.64 J. It can be seen from the comparison of the sliding-friction work before and after compensation that the sliding-friction work can be reduced by compensating the friction coefficient.

Friction Coefficient
When the initial shift speed difference was 600 r/min, the moment of inertia of the synchronized part was 0.04 Kg·m 2 . After compensation, the test and simulated friction coefficient results were compared, as shown in Figure 12. In the experiment, the temperature increased first and then decreased. At the initial stage of friction, the energy of heat flux density input was greater than the energy of heat conduction and convective heat transfer taken away, and the temperature increased speedily. With the decrease in speed difference at both ends of the synchronizer, the heat flux and convective heat transfer In the experiment, the temperature increased first and then decreased. At the initial stage of friction, the energy of heat flux density input was greater than the energy of heat conduction and convective heat transfer taken away, and the temperature increased speedily. With the decrease in speed difference at both ends of the synchronizer, the heat flux and convective heat transfer coefficient of the friction cone decreased gradually. When the energy of heat flux density input was equal to the energy taken away, the temperature of the synchronizer ring friction cone reached the maximum. As the speed difference at both ends of the synchronizer decreased continuously, the temperature of the synchronizer ring friction surface reduced gradually until the end of synchronization process. It can be seen from Figure 12 that after the compensation control, the friction coefficient simulation and test values during the synchronization of the synchronizer had been fluctuating around the expected value, which satisfied the expected control target. However, the fluctuation coefficient of the friction coefficient test value during the synchronization process was larger than expected-the expected friction coefficient was 0.08, and the friction coefficient obtained by the fitting was slightly lower than 0.08. Therefore, it was necessary to carry out a large number of tests in the later stage, and a better friction coefficient compensation control of various MAP, and consider various errors in the compensation control process to obtain a better control performance.
By comparing the synchronization time, shock degree, sliding-friction work, and friction coefficient before and after the synchronization of the synchronizer, it can be concluded that the compensation control strategy can maintain the friction coefficient of the synchronization near the expected friction coefficient. The synchronization time was shortened, and the sliding-friction work consumed in the synchronization phase was reduced. It is only the shock degree in the synchronization phase that was slightly larger, but the maximum value of the shock degree in the synchronous phase was smaller than the shock degree in the uncompensated phase. It can be seen from the performance evaluation indicators of these synchronization stages that the friction coefficient compensation control strategy can have a good effect. From the perspective of shortening the synchronization time, the reduction in the synchronization time can reduce the shift interruption time, which is crucial for shift quality.

Conclusions
The quantitative relationship between the friction coefficient and its influencing factors was studied. Based on the quantitative analysis and an optimal shift force control method, a compensation control strategy for the friction coefficient was established and verified. The main results can be summarized as follows:

1.
After the friction coefficient compensation control, the friction coefficient was maintained near the expected value of 0.08, which verified the effectiveness of the friction coefficient compensation; 2.
Under the compensation control strategy based on an optimal shift force control method, the shifting speed difference of the synchronous time was shortened by nearly 0.12 s and the sliding friction was reduced by 4.64 J under the experimental conditions. The shift quality was improved significantly.