A unified modeling and control design for precision transmission system with friction and backlash

Transcription

1 Special Issue Article A unified modeling and control design for precision transmission system with friction and backlash Advances in Mechanical Engineering 016, Vol. (5) 1 11 Ó The Author(s) 016 DOI: / aime.sagepub.com Xiulan Bao 1, Jincheng Mao and Xin Luo Abstract The structural flexibility, nonlinear friction, and backlash are the major factors limiting the control performance of precision transmission systems. If uncompensated, these factors compromise the positioning and tracking accuracy of precision transmission systems and even cause limit cycles and oscillation. In this article, a framework for integrated design from dynamic modeling to controller design is proposed. A multi-state dynamic model is presented, which can unify the modeling for a multi-state, discontinuous system including the motor state, the motion state, the mechanical contact state, and the friction state. Then, a control design method related to the dynamic modeling using perturbation separation of the model parameters is presented. Using the proposed modeling method, a continuous dynamic model is established to include all different partition models. The model comprehensively describes the mechanical and electrical characteristics of the precision transmission system. A robust controller is designed using the proposed control method. Experimental results demonstrate that the proposed modeling method is accurate and the proposed control method significantly improves accuracy and robustness of the controller compared to traditional control methods. Keywords Precision transmission system, dynamic model, robust control, friction, backlash Date received: 9 November 015; accepted: 1 April 016 Academic Editor: Nasim Ullah Introduction Precision transmission systems have been widely used in industrial manufacturing. Performance standards of precision transmission systems are becoming stricter in accordance with improved machining accuracy. 1, To achieve the increased requirements of steady state, tracking, and transient performances, the dynamic change of the precision transmission system cannot be ignored when designing the controller. The dynamic behavior is usually caused by coupled factors including non-smooth nonlinearities such as friction and backlash, mechanical system dynamic performance variables such as varying transmission ratios and the changing payload, and the motor dynamic characteristics. The motion accuracy is greatly affected by these factors. Friction is a nonlinear phenomenon present in the motion of bodies in contact, which can cause control problems, such as static errors, limit cycles, and stickslip effect, limiting the performance of the system. In precision transmission systems, it is essential to assure 1 College of Engineering, Huazhong Agricultural University, Wuhan, China School of Mechanical and Electrical Engineering, Wuhan Institute of Technology, Wuhan, China State Key Laboratory of Digital Manufacturing Equipment and Technology, Huazhong University of Science and Technology, Wuhan, China Corresponding author: Jincheng Mao, School of Mechanical and Electrical Engineering, Wuhan Institute of Technology, Wuhan, Hubei 007, China. mjincheng1999@16.com Creative Commons CC-BY: This article is distributed under the terms of the Creative Commons Attribution.0 License ( which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages ( open-access-at-sage).

2 Advances in Mechanical Engineering good dynamic characteristics by compensating friction. A precise friction model is necessary in the design of friction compensators. Many friction models were developed to describe the phenomenon of friction. According to whether differential equations are used, friction models can be broadly divided into two categories: static friction models and dynamic friction models. In the static friction models, friction is described as a function of the relative velocity. Static friction models include Karnopp model, Armstrong model, and classical models which are different combinations of Coulomb friction, viscous friction, and Stribeck effect., In dynamic friction models, friction is described as a function of the relative velocity and displacement. Dynamic friction models include Dahl model, Bliman Sorine model, Bristle model, Reset integrator model, LuGre model, Leuven model, and generalized Maxwell slip. Dynamic models can describe the nonlinear behavior of friction more accurately. However, there is a serious drawback: model parameter identification is difficult in practice. Hence, compensation techniques based on static friction models are still needed in practical applications. 5 It is essential in controller design to reduce the effects of friction on a closed loop. 6 In previous works on friction compensation, there exist three major trends: compensation using soft computing methods, observer-based friction compensation, and model-oriented friction compensation techniques. 7, The nonlinear part of the friction model is compensated using a fuzzy system in Garagic and Srinivasah 9 and Mostefai et al. 10 Coulomb friction compensation with a reduce-order observer has been proposed in Barabanov and Ortega. 11 A sliding mode adaptive control with observer is studied for servo actuators. 1 Based on nonlinear parameterizable friction, a Lyapunov-based tracking control is proposed in Makkar et al. 1 A novel continuously differentiable nonlinear friction model is first derived by modifying the traditional piecewise continuous LuGre model. 1 Backlash exists in every transmission system where a driving member (motor) is not directly connected with the driven member (load). 15 Backlash is a common phenomenon in transmission systems. 16 Compensation of backlash has attracted research efforts for several decades and there still exist many theoretical and practical problems. 15 The first approach to backlash compensation is identification of backlash and implementation of backlash inverse. 17 This approach usually requires the location of the backlash block to satisfy certain conditions, since a direct cancelation requires the backlash to be at the plant input. The second approach is to design a linear robust controller, which is a stabilizing controller, robust enough to cope with the backlash effect. 1 In the previous works, the controllers were designed to compensate friction or backlash individually. In practice, however, friction and backlash often coexist. Few studies deal with the control of systems by compensating both backlash and nonlinear friction. 5,19 A controller accounts for not only the structured uncertainties but also the unstructured uncertainties is proposed. 0 Furthermore, in almost all of the studies, the motor dynamic characteristics and the dynamic performance variables of the mechanical system are ignored. In order to improve the control precision and provide the guidance of dynamic design, a framework for integrated design from dynamic modeling to the controller design is proposed in this article. A multi-state dynamic modeling is proposed, which describes friction, backlash, motor dynamic characteristics, and structural flexibility simultaneously and can be easily introduced in robust controller design. A robust control method of perturbation separation of the parameters introduced in the model is also proposed. The rest of the article is organized as follows: in section System descriptions and modeling, a multi-state dynamic modeling method is proposed to describe the precision transmission system. In section Parameter identification of the dynamic model, the parameters of the derived model are identified, and the accuracy of the model is validated by experimental results. In section Robust controller syntheses, a new highperformance robust controller is designed. The effectiveness and robustness of the controller is also validated by experimental results. Finally, section Conclusion concludes this study. System descriptions and modeling Figure 1 shows the setup of a typical precision transmission system. In the transmission part of the system, there are a motor, a two-stage reduction gear, a screw, a nut, and a fork-lever which are typically machine parts. The two-stage reduction gear is composed of a double cylindrical gear and two cylindrical gears. The mechanical transmission ratio is modified by the screw and the nut. The friction and backlash which can cause system instability exist between the gear, screw, nut, and fork-lever. Hence, a typical precision transmission system can be divided into four distinct entities: the mechanical system, the disturbance generating system, the controller, and the spec generating system. The controlled system, which consists of the mechanical system and the disturbance generating system, is manipulated to reach the desired specs by the controlling system including the controller and the spec generating system. In the typical transmission system studied in this article, the mechanical system is used to convert highspeed rotary motion to low-speed rotation. The control performance is affected by the undesirable effects from the disturbance generating system. In a typical

3 Bao et al. Figure 1. Structure of a typical precision transmission system. transmission system, the disturbances are divided into three parts: friction and backlash caused by the gear, screw, nut, and fork-lever; the load change and transmission ratio change; and the motor dynamics characteristics. These disturbances must be taken into account in modeling. The spec generating system generates a number of desired trajectories for the mechanical system. The controller manipulates the motor to force the mechanical system to the desired trajectories according to certain rules, which are the focus of the design. Friction and backlash exist between gears, screw, nut, and fork-lever. The backlash effects between the nut and fork-lever are particularly significant. The friction effects in the fork-lever also deserve most attention. Hence, ignoring the dynamic behavior of nonmovement direction, the dynamic model of the mechanical system is sketched in Figure. In order to describe the mechanical system and the disturbances accurately and sufficiently, the mechanical system is split into four rigid bodies in which the quality of transmission shafts is ignored. The first rigid body contains the motor. The second rigid body contains the first cylindrical gear. The third rigid body contains the double cylindrical gear, the second cylindrical gear, the screw, and the nut. The fourth rigid body contains the fork-lever. The interconnection of the first and second body is a linear torsion spring and damper. The Figure. Dynamic model of mechanical system. interconnection of the second and third body is a gear pair, which represents total transmission ratio of mechanical system. The backlash, together with the linear torsion spring and damper, exists between the third and fourth body. Define u m, u p, u g, u l indicating the rotation angles of four rigid bodies. u p and u g are mutually dependent. A classical dead zone model is adopted to analyze the mechanical system with backlash. Due to backlash, the

4 Advances in Mechanical Engineering dynamic model can be divided into two regimes: the contact mode, when the torque can be transmitted from the third rigid body to the fourth rigid body, and the backlash mode, when the direction of motion changes and there is no contact between the third rigid body and the fourth rigid body. According to the above description, a dynamic model with degrees of freedom is given as follows: Contact mode (if u g u l b) J m u m + C p (_u m h_u g )+K p (u m hu g )=T m >< (J p h + J g ) u g C p _u m +(C p h + C g )_u g C g _u l K p u m +(K p h + K g )u g K g u l K g b sign(u g u l )=0 u l C g (_u g _u l ) K g (u g u l ) >: + K g b sign(u g u l )= t f T l ð1þ Backlash mode (if u g u l \b) < J m u m + C p (_u m h_u g )+K p (u m hu g )=T m (J p h + J g ) u g C p _u m + C p h_u g K p u m + K p hu g = 0 : u l = t f T l ðþ where J m, J p, J g, and denote the inertia of the four rigid bodies; _u m, _u p, _u g, and _u l denote the angular velocities of the four rigid bodies; u m, u p, u g, and u l denote the angular accelerations of the four rigid bodies; h denotes the total transmission ratio of the mechanical system; T m and T l denote the input and external load torque, respectively; K p and C p denote the spring stiffness and the damper coefficient between the first and second body, respectively; K g and C g denote the spring stiffness and damper coefficient between the third and fourth body, respectively; t f denotes the total friction which is exerted on the last rigid body; and b denotes the backlash gap size. The input torque T m is given by the input motor voltage. Considering the dynamic model, the motor dynamic characteristic can be modeled as follows U a = K e _u m + I a R a + L a _I a ðþ T m = K T I a where U a is the armature voltage, I a is the armature current, R a is the total resistance of the armature circuit, L a is the total inductance of the armature circuit, K e is the back electro motive force coefficient, and K T is the torque coefficient. In order to accurately describe friction of the mechanical system, a friction model combining stiction, Coulomb friction, viscous friction, and Stribeck effect is adopted. Taking into account the difference of friction in positive and negative directions, the friction model can written as f ns, _u l = 0 and f e f ns >< f ps, _u l = 0 and f e f ps t f = f e, _u l = 0 and f ns \f e \f ps f pc +(f ps f pc )e j_u l j= jq ps j + fpv _u l, _u l.0 >: f nc ( f ns f nc )e j_u l j= jq ns j + f nv _u l, _u l \0 ðþ where f ps and f ns are the static friction coefficients in the positive and negative direction, respectively; f pc and f nc are the Coulomb friction coefficients in the positive and negative direction, respectively; f pv and f nv are the viscous friction coefficients in the positive and negative direction, respectively; q ps and q ns are the Stribeck velocities in the positive and negative direction, respectively; and f e denotes the generalized tangential drive force exerted on the fourth body by the third body. The nonlinear friction model in equation () is difficult to be applied in controller design when backlash is also present in the controlled mechanical system. Assume that the fork-lever moves in the ½ _u lmn, _u lswp Š velocity domain, where _u lswp, _u lmn are the positive maximum velocity and negative maximum velocity of the fork-lever, respectively. A linear approximate friction model can be used for the design of the controller a p1 + b p1 _u l, if _u l (0, _u lswp Š >< a t linear = p + b p _u l, if _u l (_u lswp, _u lmp Š a n1 + b n1 _u l, if _u l ½ _u lswn, 0Š >: a n + b n _u l, if _u l ½ _u lmn, _u lswn ) ð5þ where, a p1 = f ps, a p = f pc, a n1 = f ns, a n = f nc, _u lswp =(a p1 a p )=(b p b p1 ), _u lswn =(a n1 a n )= (b n b n1 ), b p1 = f pv ( f ps f pc )=q ps, b p = f pv ( f ps f pc )e _u lmp =q ps =q ps, b n1 = f nv ( f ns f nc )=q ns, and b n = f nv ( f ns f nc )e _u lmn =q ns =q ns. Based on the above expressions, it can be easily inferred that the friction behaves are different in the low-velocity and high-velocity regimes. If the absolute value of the velocity is smaller than _u lswn, the friction may induce unstable behaviors and a limit cycle causing control problems. By approximating the friction t f with the linearized friction model t linear and considering backlash, the mechanical system with friction, backlash, and motor dynamic characteristic can be modeled as follows: If the system is in backlash mode L_I a + RI a + K e u m = U >< a J m u m + C p (_u m h_u g )+K p (u m hu g )=K T I a (J p h + J g ) u g C p _u m + C p h_u g K p u m + K p hu g = 0 >: u l = t linear T l ð6þ If the system is in contact mode

5 Bao et al. 5 L_I a + RI a + K e u m = U a J m u m + C p (_u m h_u g )+K p (u m hu g )=K T I a >< (J p h + J g ) u g + C p h_u g _u m + Cg (_u g _u l ) + K p hu g u m + Kg (u g u l b sign(u g u l )) = 0 >: u l C g (_u g _u l ) K g (u g u l ) + K g b sign(u g u l )= t linear T l ð7þ In the above model, the linear approximate friction varies with the fork-lever angular velocity. Note that transmission ratio h and payload T l vary with the fork-lever angle position, that is, h ½h 1, h Š and T l ½ T l max, T l max Š, where h 1 and h are the maximum and minimum transmission ratios, respectively, and T l max is the maximum payload. Therefore, variable payload and variable transmission are also included in the model. According to the definition above, the following equations are obtained h = h + h 1 + h h 1 d, T l = T l max d 5 t linear = a f d + b f max + b f min _u l + b f max b f min _u l d ðþ Obviously, d characterizes the dynamic change caused by varying transmission ratios; d and d characterize the dynamic change caused by nonlinear friction; and d 5 characterizes the dynamic change caused by the varying payload. Clearly, d, d, d, d 5 ½ 1, 1Š. In addition, let d 1 represent the model change caused by backlash, where d 1 ½ 1, 1Š. Then, a multi-state dynamic model is derived, which includes the electric state and different dynamic states of all cases ðj d + DJ d Þ q + ðc d + DC d Þ_q + ðk d + DK d Þq = F ð9þ T, where q = I a, u m, u g, u l Jd = diagf0, J m, h J + h 1 T, p + J g, g, F = U a, 0, 0, T l max d 5 a f d h and DJ d = diag 0, 0, J h 1 p d, 0. The matrices of stiffness and damp can be expressed as follows K d = K d11 K d1 K d1 K d R K T K p h + h 1 K p 0 = h + h 1 0 K p K p + K g K g K 5 g K g ð10þ C d = C d11 C d1 C d1 C d L K e C p h + h 1 C p 0 = h + h 1 0 C p C p + C g C g C g C g + b 7 f max + b 5 f min ð11þ h h 1 6 K p d DK d = h h 1 K p d + K g 0 0 K g d 1 d 1 K g d 1 7 K 5 g d h h 1 C p d 0 DC d = 0 0 h h 1 C p d + C g 6 d 1 Cg d Cg d C g 1 d 1 + b f max b f min ð1þ d 7 5 ð1þ Parameter identification of the dynamic model Parameter identification The imperfection of the transmission model leads to undesirable behaviors in the system. Therefore, proper parameter identification of the system dynamic model is important for realizing good performance. There are five types of parameters to be determined: the inherent property parameters of the mechanical system; the electrical parameters; and the corresponding coefficients of the friction model, the backlash gap size, and the maximum payload of the mechanical system. The inherent property parameters of the mechanical system include inertia of each rigid body, stiffness and damping coefficients of each spring, and transmission ratio. The rigid body inertias are calculated based on the shape and size of the corresponding components. The spring stiffness and damping coefficients are obtained through the finite element analysis of corresponding components using the ANSYS FEM software package. The transmission ratio, which varies with the fork-lever angle position, can be expressed as h = h0 cos ½h 1, h Š ð1þ u l where h 0 is the transmission ratio when the nut is located at the center of the screw. The value of h 0 can

6 6 Advances in Mechanical Engineering velocity is assumed to be equal to the friction torque. The coefficients of friction model, as shown in Table 1, are identified by fitting the corresponding Stribeck function equation () to the experimental data that minimizes the quadratic cost function. By fixing the motor input shaft and measuring the moveable angle of the fork-lever, the backlash gap size is obtained. The maximum payload of the mechanical system is the resistance generated when the fork-lever is in the maximum deflection angle under the working condition. The specific values of all parameters are shown in Table 1. Figure. Friction curve: (*) measurement data and ( ) fitted curve. be easily calculated through assembly relationship among the components of the mechanical system. Therefore, it is easy to calculate h 1 and h, where h 1 = h 0. The electrical parameters can be obtained from the motor performance specifications. The friction curve, which is depicted in Figure, is obtained from a closed-loop experiment where the precision transmission system is controlled with a proportional integral derivative (PID) controller at several constant velocities and the input torque at a constant Dynamic model verification In order to validate the derived dynamic model, the time-domain experiment and simulation are completed. Due to safety requirements and operating conditions of the precision potentiometer, the time-domain experiment is accomplished under closed-loop conditions. A PID controller is used in the experiment to form a closed-loop system. The experimental setup is equipped with a precision potentiometer to measure the angle position of the fork-lever on discrete time samples via a digital signal processor. Additional external reference signals can be applied to the feedback connection of the positioning mechanism to provide sufficient excitation while gathering data for verification. In order to explicitly verify the various effects of the model, especially the friction effects, sinusoidal signals are used as the external reference signal. Table 1. Parameter values of the transmission system. Mechanical inherent property parameters Type Parameter Value Inertia J m (kg m ) l J p (kg m ). 10 J h (kg m ) J 1 (kg m ) Spring stiffness K p (N m/rad).7 10 K g (N m/rad) Damper coefficient C p (N m s/rad) 1 10 C 6 (N m s/rad). 10 Transmission ratio h h. 10 Electrical parameters Total inductance of armature circuit La (H) Total resistance of the armature circuit Ra (O) 1.55 Back electro motive force coefficient K c (V s/rad).9 10 Torque coefficient K T (N m/a).5 10 Friction parameters Positive static friction coefficient f ps (N m) Negative static friction coefficient f ns (N m) Positive Coulomb friction coefficient f pc (N m) Negative Coulomb friction coefficient f nc (N m) Positive viscous friction coefficient F pv (N m s/rad).1 10 Negative viscous friction coefficient f nv (N m s/rad) Positive Stribeck velocity q ps (rad/s) 10 Negative Stribeck velocity q ns (rad/s) 10 Backlash gap size b (m) Maximum payload T Imax (N m) 10 1

7 Bao et al. 7 Figure. Tracking results of dynamic model time-domain simulation and experiment. In order to match the experimental conditions, closed-loop simulation is performed in which the controlled object is the above dynamic model, and the controller is the same as in the experiments. Figure plots the tracking results for a command specifying a - and 1-Hz sinusoidal trace in 1 s. In Figure, the tracking amplitude errors of 0.1 in the time intervals ½0: s, 0:5 sš and ½0:7 s, 0:5 sš occurred when the speed of the fork-lever was near zero. The tracking errors, which are caused by friction and backlash, are accurately captured by our dynamic model. Good agreements between the simulation results and experimental results can be seen, which validates our dynamic model. Robust controller syntheses Therefore, design of the advanced robust controller can be based on the proposed accurate and simple dynamic model to minimize negative factors caused by dynamic changes, which improve the tracking performance. Perturbation separation of the model parameters The system perturbations of the unified multi-state dynamic model (9) caused by dynamic changes need to be separated for controller design. Using the descriptor form, a model is described as follows in which parameter perturbations and nominal system are separated < E + P d i E i _x = A + P d i A i x + B + P d i B i u + Hd : i = 1 i = 1 i = 1 y = Cx + Qd ð15þ where x = ½I a, _u m, _u g, _u l, u m, u g, u l Š T are the state variables; d 1, d, and d are considered to be the model parameter perturbations; and d and d 5 are considered to be the hypothetical external disturbance perturbations, that is, d = ½ d d 5 Š T. u = U a is the input signal, and y is the fork-lever angular position. The corresponding matrices of the nominal system and the parameter uncertainty can be expressed as follows K d11 C d1 0 A = K d1 C d K d 5 ð16þ 0 M 0 M 1 M 5 M 5 M 5 M 5 M 5 M C g C g 0 K g K g A 1 = C g M C g K g 0 K 6 g 7 5 ð17þ M M M M M M M 0 M 0 0 h h 1 h h 1 M C p M K p 0 A = M h h 1 C p M h 6 h 1 7 K p 0 5 M 7 M 6 M 7 M 6 M 6 ð1þ M 1 = 0 R, M = 0 R 1, ð19þ M = 0 R, M = 0 R 1 M 5 = 0 R 1, M 6 = 0 R 1, M 7 = 0 R ð0þ A = diag 0, 0, 0, b f max b f min, 0, 0, 0, ð1þ Q = 0 R 1 " H = # a f T T, M = diagf1, 1, 1g l max ðþ h + h 1 E = diag L, J m, J p + J g,, 1, 1, 1, ðþ E 1 = E = 0 R 7 7 h h 1 E = diag 0, 0, J p, 0, 0, 0, 0, ðþ C = ½ Š B 1 = B = B = 0 R 7 1, B = ½ Š T ð5þ Robust controller design The control goal is to design a feedback controller that stabilizes the closed-loop system and keeps the servo errors and the feedback control actions within required bounds. The dynamics of the precision transmission system strongly depend upon the motion status at any

8 Advances in Mechanical Engineering min W del (I + G M) 1 G 1 W del (I + G M) 1 G ðþ The system performance weighting function W del is set according to the performance requirements of the precision transmission system. Through the above process, the controller can be designed with the identified parameters, which can be expressed as follows Figure 5. Robust synthesis block diagram. M = 0:009s 5 + 0:9s + :1s + 175s + 70s :00001s 6 + 0:009s 5 + 0:51s + 0:7s + 109s + 110s ð9þ moment. Thus, the controller must be robust against these dynamic changes. As shown in Figure 5, the motion control system performing the precision transmission servo-control mainly consists of the nominal dynamic model (G p ), the feedback controller (M), the model parameter uncertainty (D 1 ), the hypothetical external disturbance uncertainty (D ), and the system performance weighting function W del. D 1 = diagfd 1, d, d g and D = diagfd, d 5 g. The input and output of the parameter uncertainty model are denoted as w 1 and z, respectively. The disturbance signal and desired trajectory are denoted by w and r, respectively. The servo error and system output are denoted by e and y, respectively. The feedback control action is denoted by u. The sensors are assumed to be perfect, so that the sensor noise can be neglected. In Figure 5, the model in the dotted-line box can be easily transformed by model (15). The nominal dynamic model G p of the precision transmission system can be expressed as follows A _ S 1 H _ B _ G p = S S A _ S S 1 S H _ S S B _ 5 ð6þ C 0 Q where A _ = E 1 A, B _ = E 1 B, H _ = E 1 H, S i = ½S i1 S i S i Š, and E 1 ½ E i A i B i Š = S 1i ½ S i S i S i Š, i = 1,,. According to the dimension of the input and output of the precision transmission system, G p can be partitioned as follows G p = G 11 G 1 G 1 ð7þ G 1 G G Note that the feedback connection in Figure 5 is internally stable only if kg 11 k 1 and kg 1 k 1. Under the stability conditions of the system, a controller can be designed so that the infinity norm of the transfer function from fv 1, v g to e is minimized, which can be expressed as Control system evaluation To verify the effectiveness of the controller, we carried out experiments comparing the performance of the proposed controller in this article with a friction compensation controller. The structure of experiment device is shown in Figure 6. As shown in Figure 7, the friction compensation is based on the identified friction model in section Parameter identification ; the controller in closedloop is PID controller. In the experiments, payload disturbance exists which proportional to the angle of the fork-lever; the maximum payload is present in Table 1. In the experiments, the sampling frequency is 1000 Hz. The sinusoidal wave signal is used as excitation signal. Figure plots the experimental results for the friction compensation controller and the proposed robust controller. The friction of the systems caused the flat-topped phenomenon in control performance. In Figure, flattopped phenomenon is clearly shown by amplitude errors of 0.1. The backlash causes tracking delay phenomenon, which is shown by blue line in Figure (b). In Figure (a), a sinusoidal input signal with Hz frequency and 10 amplitude is given; good tracking performance can be achieved with both controllers. The tracking errors caused by friction and backlash, as shown in Figure, are well compensated by the two controllers. The trajectory of the friction compensation controller, shown by red line in the figure, has amplitude overshoot at the peaks and troughs of sinusoidal signal, especially in the troughs. The trajectory of our proposed robust controller in this article have no amplitude overshoot at the peaks and troughs. And the phase lag of the proposed controller is 1.76 which is smaller than the 5.67 phase lag of the friction compensation controller. In Figure (b), a sinusoidal input signal with 0 Hz frequency and 10 amplitude is given. The decrease in amplitude of the friction compensation controller is approximately 50%, and the phase lag is There is a large tracking error using friction compensation

9 Bao et al. 9 Figure 6. Experimental schematic representation of the precision servo transmission system. Figure 7. Schematic representation of the friction compensation controller. Figure. Experimental tracking performance: (a) at low frequencies ( Hz) and (b) at high frequencies (0 Hz). controller for high-frequency signals. This is because the friction compensation controller is only for a fixed friction model. Our proposed controller has good agreement with the input signal in amplitude, and the phase lag is only.5. The frequency characteristics of precision transmission system are shown in Figure 9; one adopt the friction compensation controller and another adopt the proposed robust controller. The two blue lines and two black lines are the frequency domain characteristics

10 10 Advances in Mechanical Engineering when using the proposed robust controller. The onefrequency domain characteristics are obtained using 1 amplitude sinusoidal excitation signal series, and the another frequency domain characteristics are obtained using 10 amplitude sinusoidal excitation signal series. This phenomenon of change in frequency characteristics also implies that the precision servo transmission system is a nonlinear system, which is different from the linear system. The two black lines are the frequency domain characteristics when using friction compensation controller, also using 1 and 10 amplitude sinusoidal excitation signals. The db bandwidth is extended from 16.6 rad/s by adopting the friction compensation controller to 61 rad/s by adopting the proposed robust controller, using 1 amplitude sinusoidal excitation signal series under the same conditions. And the bandwidth is extended from rad/s by adopting the friction compensation controller to 5.1 rad/s by adopting the proposed robust controller, using 10 amplitude sinusoidal excitation signal series under the same conditions. The phase lags both have been greatly improved. From Figure 9, the proposed robust controller is better than the friction compensation controller. As mentioned above, because of the compensation for friction, the friction compensation controller has good performance in low frequency, but cannot resist the change of gap and payload in high frequency. Apparently, the proposed controller is much better than the friction compensation controller for the precision transmission system, not only can resist the friction effects at low frequencies but also can resist the change of gap and payload in high frequencies. The robustness and effectiveness of the proposed controller can be clearly seen. Conclusion In this article, a multi-state dynamic model which is used for control design is presented, and a controller design method using perturbation separation of the dynamic model parameters is introduced. The effects of the dynamic performance variables of the mechanical system, the motor dynamic properties, the nonlinear characteristics of friction, and backlash on the accuracy of motion were discussed. A dynamic model for a typical precision transmission system was established. Identification of model parameters based on experimental analyses was presented. Experimental results demonstrated that the proposed dynamic model was accurate and effective. A robust controller was also designed for the dynamic model. Closed-loop control experimental using the proposed robust controller compared with friction compensation controller was carried out. Experimental results show that our proposed robust controller can be achieved by the good tracking performance and greatly improves the tracking accuracy in time domain, especially for high-frequency signal tracking. Compared with friction compensation controller frequency response experimental results, the bandwidth and phase lags both have been greatly improved in frequency domain by adopting the proposed robust controller. The proposed framework for integrated design from dynamic modeling to controller design is verified to be feasible and correct. In this article, the dynamic modeling and controller design are synthesized in a unified form. The proposed method can be modeled for discontinuous vary system. In the process of modeling, the motor and motion state, the mechanical contact state, and the friction state are considered. Through the proposed integrated design process, every model parameters that affect the control performance were analyzed clearly. The proposed method can be used to complex electromechanical system with friction and backlash. Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article. Figure 9. Frequency characteristics of precision transmission systems. Funding The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was supported by the Scientific Research Funds of the Educational Commission of Hubei Province (grant no. Q ), the Natural Science Foundation of Hubei Province, China (grant no. 01CFB), and the Fundamental Research Funds for the Central Universities of China (grant no. 016PY017, 01QC007).

11 Bao et al. 11 References 1. Chen Q, Zhao YC and Yin JB. Automated mechanical transmission technology and development trend. Appl Mech Mater 01; 697: 0.. Zhang H, Zhong Z and Chen X. Systematic synthesis method and prototyping of fixed-axle vehicular electrified mechanical transmission. Int J Automot Techn 015; 16: Liu LL, Liu HZ and Wu ZY. An overview of friction models in mechanical systems. Adv Mech 00; : Pascal B. Reduced nonlinear observer for bounded estimation of the static friction model with the Stribeck effect. Syst Control Lett 009; 5: Lee W, Lee CY and Jeong YH. Friction compensation controller for load varying machine tool feed drive. Int J Mach Tool Manu 015; 96: Lotfi M, Mouloud D and Yoichi H. Robust tracking controller design with uncertain friction compensation based on local modeling approach. IEEE/ASME T Mech 010; 15: Marton L and Lantos B. Modeling, identification, and compensation of stick-slip friction. IEEE T Ind Electron 007; 5: Lee TH, Tan KK and Huang S. Adaptive friction compensation with a dynamical friction model. IEEE/ASME TMech011; 16: Garagic D and Srinivasah K. Adaptive friction compensation for precision machine tool drive. Control Eng Pract 00; 1: Mostefai L, Denai M and Sehoon O. Optimal control design for robust fuzzy friction compensation in a robot joint. IEEE T Ind Electron 009; 56: Barabanov N and Ortega R. Necessary and sufficient conditions for passivity of the LuGre friction model. IEEE T Automat Contr 000; 5: Xie WF. Sliding-mode-observer-based adaptive control for servo actuator with friction. IEEE T Ind Electron 007; 5: Makkar C, Sawyer WG and Hu G. Lyapunov-based tracking control in the presence of uncertain nonlinear parameterizable friction. IEEE T Automat Contr 007; 5: Yao JY, Deng WX and Jiao ZX. Adaptive control of hydraulic actuators with LuGre model-based friction compensation. IEEE T Ind Electron 015; 6: Nordin M and Gutman PO. Controlling mechanical systems with backlash a survey. Automatica 00; : Voros J. Modeling and identification of systems with backlash. Automatica 010; 6: Selmic RR and Lewis FL. Backlash compensation in nonlinear systems using dynamic inversion by neural networks. Asian J Control 000; : Barreiro A and Banos A. Input-output stability of systems with backlash. Automatica 006; : Chen SL, Tan KK and Huang S. Modeling and compensation of ripples and friction in permanent-magnet linear motor using a hysteretic relay. IEEE/ASME T Mech 010; 15: Yao JY, Jiao ZX, Ma DW, et al. High-accuracy tracking control of hydraulic rotary actuators with modeling uncertainties. IEEE/ASME T Mech 01; 19: 6 61.