Dynamic Torque Coordinated Control Considering Engine Starting Conditions for a Power-Split Plug-In Hybrid Electric Vehicle

: For a power-split plug-in hybrid electric vehicle (PS-PHEV), the control system cannot effectively reduce jerks under some engine-starting conditions due to the pulsating torque of the engine, elastic characteristics of the transmission, battery power limitations, and the limitation of the maximum torque of the motor. To solve this problem, a dynamic torque-coordinated control method considering engine starting conditions (DTCC-ESC) is proposed in this paper. The proposed DTCC-ESC mainly includes three parts: the engine segment active control before and after engine ignition to minimize the impact of the engine ripple torque on the shocks of the powertrain system; the feedforward control of the engine starting conditions and the feedback control of the engine’s optimal target speed; and the active damping feedback compensation control for system resonance to reduce system vibrations. The results show that the proposed DTCC-ESC can effectively improve driving comfort under different engine starting conditions.


Introduction
The plug-in hybrid electric vehicle (PHEV) is one kind of new energy vehicle that has attracted widespread attention from automobile manufacturers [1]. Compared with traditional internal combustion engine (ICE) vehicles, PHEVs have better system efficiency, fuel economy and emission performance because the driver can choose different working modes [2,3]. During the mode-switching process of starting the engine, the vehicle driveline will often produce obvious shock and vibration. For a power-split PHEV, the engine ripple torque (ERT) will be directly transmitted to the powertrain system [4]. In the process of mode transition (MT), the transmitted ERT will cause obvious impacts and vibrations on the power transmission system, especially during engine starting [5,6]. Such transients make driving less comfortable and reduce the service life of the transmission. Some control methods have been applied to the typical engine starting process to reduce vehicle shocks. Huang [7] formulated dynamic and economic shift schedules in an MT process and coordinated the motor torque to compensate for the fluctuation torque. In [8], considering the inertial torque of the motor, a motor compensation coordinated control strategy was designed to reduce unwanted jerks. Taking the speed fluctuation of engine and the fuel consumption rate as the control objective, Zuo [9] designed a modeswitching control strategy based on a logical threshold. Gao [10] proposed a control strategy of torque-dynamic coordination by means of a sliding mode control based on a disturbance compensation. These methods made the coordinated control strategy during MT process reduce vehicle jerks, but they ignored the influence of the ERT and elastic characteristics of the transmission system on the shocks of the powertrain.
Based on the data from an engine start-up test, researchers in [11,12] proposed a lookup table model of the ERT during engine start-up to estimate the engine pulsating torque and then used motors to compensate or offset the torque fluctuation of the powertrain. Canova [13] and Davis [14] set up a reduced-order engine thermodynamic model to estimate the ERT and used a belted starter/alternator to reduce the shocks and vibrations caused by ERT. Zhao [15] analyzed the MT process during the engine start-up and used engine torque estimation based on planetary dynamic equations to reduce vehicle jerks. In the implementation of such schemes, these methods helped estimate ERT to reduce vehicle jerks. Shen [16] designed a coordinated control strategy to determine target engine torque, motor torque, clutch torque and the moment of fuel injection to suppress vehicle jerks and improve riding comfort during engine start. Lin [17] divided the engine start into three phases, analyzed the influences of various parameters on the jerk in each phase, and then proposed a coordinated control strategy between the clutch pressure and motor torque and between the motor torque and engine torque. However, the influence of the elastic (stiffness-damping) characteristics of transmission on the powertrain vibrations was not taken into account.
Tomura [18] used the approximate function of ERT and the proportional feedback control of the speed difference to reduce the fluctuations of the powertrain. Wang [19] proposed an active damping-control strategy during the starting process to reduce the speed fluctuations of the transmission system. Su [20] used a staged torque estimation of the engine and an active damping-compensation control method based on drive shaft torque observations to reduce the jerks of the powertrain. Li [21] proposed a fuzzy control and a two degrees of freedom damping-control algorithm to reduce vehicle jerks. In these control methods, the ERT and the elastic characteristics of the transmission system were considered, but the influence of the limitations of battery-motor characteristics on vehicle jerks under different engine starting conditions was ignored.
In [22][23][24], the dynamic programming algorithm was used to analyze and formulate the optimal brake pressure, engine speed trajectory and control logic of MT process to minimize vehicle jerks. However, these methods cannot be applied to actual vehicle control.
In [25,26], the motor torque change rate was limited using dynamic equations of planetary gear to reduce vehicle jerks. Zeng [27] used feedforward control for one motor and compensation control for the other based on the engine torque estimation to improve driving comfort. Syed [28] designed an active damping-control system at the wheels to suppress system oscillation. Tang [29] designed a motor torque compensation control method and established a vibration path control to change the dynamic characteristics of the engine during starting and stopping. Zhu [30] divided the MT process into four operational stages and adopted the fuzzy gain scheduling proportion integration differentiation (PID) controller to reduce fluctuations of the powertrain system. To some extent, although the above methods can reduce vehicle jerks and vibrations, they also ignore the influence of the limitations of battery-motor characteristics on vehicle jerks under different engine starting conditions.
The above research during engine start-up for a hybrid powertrain system mainly focuses on the engine torque estimation, the limitation of the motor torque change rate, the compensation control based on the active damping and dynamic programming algorithm and so on. However, these methods do not fully consider the influence of factors such as ERT, the elastic characteristics of the transmission system, the battery power limitations or the maximum motor torque. In fact, these factors will diminish the effectiveness of the control system to reduce vehicle jerks under different engine starting conditions.
To solve this problem, this paper proposes a dynamic torque coordinated control method during engine starting conditions (DTCC-ESC). In comparison with previous research, the main contributions of this paper include the following three aspects: (1) In the MT process during engine starting, an engine segment active control method before and after engine ignition is proposed to minimize the impact of ERT on system jerks and vibrations; (2) During starting, the engine demand-angular acceleration is optimized to modify the feedforward control of the engine starting process before ignition; and (3) Based on the zero-pole distribution method and the estimation algorithms of the system torque and speed, an active-damping feedback compensation control for system resonance (ADFC-SR) is modified to restrain the vibrations of the transmission system. By this method, the torque of the engine output shaft fluctuations can be reduced; meanwhile the torque of the two motors to compensate or offset driveline fluctuations can be optimally distributed, so that the compensation torque from MG1 and MG2 will not exceed the limit of its own maximum torque. Therefore, vehicle jerks in some engine-starting conditions can be effectively reduced.
The rest of this article is organized as follows: In Section 2, the model of a compound power-split power transmission system is established, and this model has been verified by some experiments. In Section 3, the dynamic analysis during the engine starting process is carried out by using the equivalent lever method, and the problem to be solved in this paper is further described by using a flow chart. To address this problem, the DTCC-ESC method is proposed in Section 4. Then, in Section 5, the three control methods are simulated and compared under different engine-starting conditions. Finally, the simulation conclusions of the proposed control method are provided in Section 6. Figure 1 shows a compound power-split plug-in hybrid electric vehicle (PS-PHEV), which is studied in this article. The compound PS-PHEV includes an engine, a power battery pack, two motors, a compound power-split transmission based on a Ravigneaux planetary gear train, two brakes (B 1 and B 2 ), a torsional damper spring (TDS) and a main reducer and differential. The Ravigneaux planetary gear train includes the small sun gear S 1 , the large sun gear S 2 , the front short planetary gear P b , the rear long planetary gear P a , the front and rear common planet carrier C, and the front and rear common gear ring R. In the powertrain system of the compound PS-PHEV, the engine is connected to the planet carrier C through the TDS. The planet carrier C can be braked by the brake B 1 . The MG1 connects to the sun gear S 1 , which can be braked by the brake B 2 . The MG2 connects to the large sun gear S 2 . The ring gear R is the output terminal of power, which is directly connected to the final drive.  Figure 1 shows a compound power-split plug-in hybrid electric vehicle (PS-PHEV), which is studied in this article. The compound PS-PHEV includes an engine, a power battery pack, two motors, a compound power-split transmission based on a Ravigneaux planetary gear train, two brakes (B1 and B2), a torsional damper spring (TDS) and a main reducer and differential. The Ravigneaux planetary gear train includes the small sun gear S1, the large sun gear S2, the front short planetary gear Pb, the rear long planetary gear Pa, the front and rear common planet carrier C, and the front and rear common gear ring R. In the powertrain system of the compound PS-PHEV, the engine is connected to the planet carrier C through the TDS. The planet carrier C can be braked by the brake B1. The MG1 connects to the sun gear S1, which can be braked by the brake B2. The MG2 connects to the large sun gear S2. The ring gear R is the output terminal of power, which is directly connected to the final drive.     Figure 2 shows the working principle of a single-cylinder engine. It includes intake and exhaust components, mechanical moving parts, fuel injectors, and spark plugs. The engine dynamics model mainly includes two parts: the engine thermodynamic model in the cylinder and the engine pulsating torque model. A powertrain plant model will be established using MATLAB/Simulink based on formulations and look-up tables as follows: Figure 2 shows the working principle of a single-cylinder engine. It includes intake and exhaust components, mechanical moving parts, fuel injectors, and spark plugs. The engine dynamics model mainly includes two parts: the engine thermodynamic model in the cylinder and the engine pulsating torque model.

Engine Thermodynamic Model in the Cylinder
The general formula of the engine working process was obtained using the first thermodynamic formula [31,32]. In the i-th cylinder (I = 1,2,3,4), the cylinder pressure P cyl,i (ϕ) can be obtained by Formula (1): here, where T temp is engine cylinder temperature; m is the gas mass in the engine cylinder; ϕ is the crankshaft rotation angular displacement; u is the specific thermodynamic energy; V cyl is the engine's instantaneous work volume; λ is the engine excess air ratio; R is a gas constant value; h is the specific enthalpy; V s is the engine's physical work volume; ε c is the engine compression ratio, and λ s is the engine crank-link ratio. Q is the engine exchange heat of the cylinder gas. Subscript s c represents the gas entering the engine cylinder. Subscript e represents the gas entering the engine cylinder through the exhaust; subscript B represents the heat released from gasoline combustion; subscript w is the heat exchange between the wall and the outside; λ s = r c /l, l is the engine-connecting rod length; and r c is engine crankshaft radius.

Engine Pulsating Torque Model
For a gasoline engine with four cylinders, the dynamic torque T ERT of the engine in front of the engine flywheel can be expressed as the sum of the torque of each cylinder with a difference of 180 • , that is, [20,26], where T cyl,i is the dynamic torque of a single cylinder; T P i is the pulsation torque of the cylinder pressure; T I i is the engine inertia torque; and T F is the engine friction torque.

Model of Torsional Damper Spring
According to the torsional stiffness characteristic curve (as shown in Figure 3) and the hysteresis damping moment of the torsional damper, the torsional damper model is established by using the

Motor Dynamic Model
Considering the battery SOC and the characteristics of the two motors, a three-phase permanent magnet synchronous motor (PMSM) with sinusoidal back electromotive force is built in the rotor reference frame (q-d frame). The PMSM is connected to the battery through a DC/DC converter and a three-phase inverter [30], and the battery provides energy to the motor. The basic electromagnetic torque model of the motor is as follows [26,33]: s 1. [ where,

Coupling Mechanism Modeling
For the power-split transmission, the motion equations of the compound planetary gear model are given using the centralized mass method [20,34] as follows: where I is the inertia of the component; r is the base circle radius of gear;  is the meshing angle of a pair of gears;  is the gear rotation angle; T is the input torque; and y F is the meshing force between gears [34]. Furthermore, a and b are the rear and front row

Motor Dynamic Model
Considering the battery SOC and the characteristics of the two motors, a three-phase permanent magnet synchronous motor (PMSM) with sinusoidal back electromotive force is built in the rotor reference frame (q-d frame). The PMSM is connected to the battery through a DC/DC converter and a three-phase inverter [30], and the battery provides energy to the motor. The basic electromagnetic torque model of the motor is as follows [26,33]: where, p s is the number of pole pairs; λ a is the amplitude of induced magnetic flux; L q and L d are the inductance of axis q and axis d, respectively; and i q and i d are the current of axis q and axis d, respectively.

Coupling Mechanism Modeling
For the power-split transmission, the motion equations of the compound planetary gear model are given using the centralized mass method [20,34] as follows: F ys 2 a n · r s 2 (6) . θ a n c = F ya n b n · r a − F ys 2 a n · r a (9) .
where I is the inertia of the component; r is the base circle radius of gear; α is the meshing angle of a pair of gears; θ is the gear rotation angle; T is the input torque; and F y is the meshing force between gears [34]. Furthermore,a and b are the rear and front row planetary gears of the compound planetary gears, and T TI is the equivalent TI shaft torque in Figure 3. N (n) is the number of planetary gears, and I ce = I c + N I a + N I b . Note that the vehicle equivalent inertia to the common gear ring is taken into account in I R .

Brake Dynamic Model
The brake model has three working states: disengagement, sliding, and engagement [26]. The transmission torque T B of the brake is expressed as where . θ c_g is the velocity of the planetary frame; µ sl is the sliding friction coefficient of the friction plate; µ st is the static friction coefficient of the friction plate; N B is the friction disk number; r B is the friction radius of the friction plate; A B is the friction area; P B is the oil pressure; and P 0 is the pre-filled oil pressure.

Vehicle Dynamic Model in Longitudinal Direction
The longitudinal dynamic driving model of the vehicle is expressed as where C TI and k TI are the damping and stiffness at the TI shaft in Figure 8; . θ R is the ring speed; . θ L is the rotation speed at the point B in Figure 8; .. θ R is the ring angular acceleration; .. θ L is the rotation acceleration of point B; i 0 is the speed ratio of the main reducer; r v is the radius of the wheel; g is the acceleration of gravity; m v is the equivalent mass of the vehicle; V x is the speed of the car; f is the rolling resistance coefficient; β is the road slope; T f is the wheel load torque; C d is the coefficient of air resistance; ρ is the density of air; and A is the windward area of the vehicle.

Validation of PS-PHEV Powertrain Model
A PHEV test bench was built to verify the validity of the established model as shown in Figure 4. The engine startup process is tested when the dynamometer is running at 27 rad/s. As shown in Figure 5a, the MT process with the engine startup includes four work stages: (1) pure electric stage with B1 locked, (2) pure electric stage with B1 opened, (3) a motoring stage before engine ignition, and (4) a hybrid driving stage after engine ignition. The engine shaft torque and speed obtained by simulation are basically consistent with the test results. As shown in Figure 5b, the driving shaft torque and speed obtained by simulation are basically consistent with the test results. Figure 5c shows the error between the simulation results and the test results of the drive shaft speed. It can be seen from these results that the maximum error of the drive shaft speed is less than 1.2 rad/s, and the relative error of the driveshaft speed is less than 5%. By comparing and analyzing the simulation results with the test results, the established model can basically reflect response characteristics of the power transmission in the mode switching process. The powertrain system model lays the foundation for the subsequent evaluation of the effectiveness of the proposed dynamic torque-coordinated control method. The main PS-PHEV parameters are listed in Table 1     effectiveness of the proposed dynamic torque-coordinated control method PHEV parameters are listed in Table 1.     Figure 6 shows the lever diagram in the engine startup process. The black solid line in this figure is the lever balance line. The speed on the balance line is defined as zero. The speed below the balance line is negative, and the speed above the balance line is positive. The up arrow = indicates that the torque is positive, and the down arrow indicates that the torque is negative.  Figure 6 shows the lever diagram in the engine startup process. The black solid line in this figure is the lever balance line. The speed on the balance line is defined as zero. The speed below the balance line is negative, and the speed above the balance line is positive. The up arrow = indicates that the torque is positive, and the down arrow indicates that the torque is negative. Figure 6. Lever chart during engine starting process.

Dynamic Analysis During Engine Starting Process
As shown in Figure 6, in the process of engine starting, MG1 torque pluses 1 , to start the engine and bring the engine to its ignition speed. The output power of MG1 is negative, indicating that MG1 is in the power generation state. MG2 adopts a compensation torque 2 , on the basis of 2 to ensure the vehicle runs forward. Ignored were the stiffness-damping elastic effect on the MG1 shaft, MG2 shaft, engine output shaft, drive shaft, as well as the dynamic meshing force and efficiency loss between gears. Based on the lever method [35], the torques of MG1 and MG2 during the engine startup can be expressed as follows:   As shown in Figure 6, in the process of engine starting, MG1 torque pluses T , MG1 to start the engine and bring the engine to its ignition speed. The output power of MG1 is negative, indicating that MG1 is in the power generation state. MG2 adopts a compensation torque T , MG2 on the basis of T MG2 to ensure the vehicle runs forward. Ignored were the stiffness-damping elastic effect on the MG1 shaft, MG2 shaft, engine output shaft, drive shaft, as well as the dynamic meshing force and efficiency loss between gears. Based on the lever method [35], the torques of MG1 and MG2 during the engine startup can be expressed as follows: where T MG1 is the MG1 torque; T MG2 is the MG2 torque; T r_req is the ring torque required; T Eng is the engine torque; I MG1 is the moment of inertia of MG1; and I MG2 is the moment of inertia of MG2. Additionally,I S1 is the moment of inertia of the little sun wheel; I S2 is the moment of inertia of the big sun wheel; i 0 is the speed ratio of the main reducer; .. θ W H is the angular acceleration of the wheel; I Eng is the moment of inertia of the engine. I C is the moment of inertia of the planetary frame; .. θ Eng_des is the engine acceleration; i 1 is the front-row speed ratio; and i 2 is the rear-row speed ratio.

Problem Description during Engine Starting Process
Fuel economy can be effectively improved using the switching working modes for a PHEV [36,37]. However, as shown in Figure 7 for the compound PS-PHEV, during the mode switching process, which includes the engine start-up process, the distribution torque command of the two motors can be obtained by the superposition calculation of the motor distribution torque from Equation (14) and the active damping compensation control (ADCC) torques [19,20]. However, in some engine starting conditions, due to the ERT, the elastic characteristics of the transmission system, limitations of battery power, and limitations of motor maximum torque, the torques of the two motors, are limited, which make them unable to effectively compensate or offset the fluctuating torque of the powertrain. These factors will result in obvious transmission jerks and resonance, which will eventually make driving less comfortable and reduce the service life of the vehicle.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 10 of 2 mode switching process, which includes the engine start-up process, the distributio torque command of the two motors can be obtained by the superposition calculation o the motor distribution torque from Equation (14) and the active damping compensatio control (ADCC) torques [19,20]. However, in some engine starting conditions, due to th ERT, the elastic characteristics of the transmission system, limitations of battery power and limitations of motor maximum torque, the torques of the two motors, are limited which make them unable to effectively compensate or offset the fluctuating torque of th powertrain. These factors will result in obvious transmission jerks and resonance, whic will eventually make driving less comfortable and reduce the service life of the vehicle. Therefore, in view of the above problem, it is necessary to design a dynamic torque coordinated control that considers engine starting conditions to further reduce the system shocks and vibrations.

Design of Dynamic Torque-Coordinated Control That Considers Engine Starting Conditions
To solve the problem of insufficient motor torque compensation under some engin starting conditions, as illustrated in Figure 8, a dynamic torque-coordinated contro method that takes into account engine starting conditions (DTCC-ESC) is proposed in thi paper for the first time. The DTCC-ESC mainly includes three parts: the engine segmen active control, the feedforward and feedback control of engine starting, and the activ damping feedback compensation control for the resonance of the transmission system.
(1) The engine segment active control, also known as the engine transient torque contro before and after the engine ignition, is used to reduce the impact of the ERT on system shocks and vibrations in the engine startup from the perspective of the engine powe source. Regarding this engine segment active control: before ignition, the inputs o control are the engine's intake valve close timing (IVCT) and initial crankshaft angl (ICA); after ignition, the inputs of control are fuel injection and engine ignition timin optimized by a genetic algorithm based on the engine working cycle. The output o the control is the engine output torque. (2) Feedforward and feedback control of engine starting is mainly realized by solvin Therefore, in view of the above problem, it is necessary to design a dynamic torquecoordinated control that considers engine starting conditions to further reduce the system shocks and vibrations.

Design of Dynamic Torque-Coordinated Control That Considers Engine Starting Conditions
To solve the problem of insufficient motor torque compensation under some engine starting conditions, as illustrated in Figure 8, a dynamic torque-coordinated control method that takes into account engine starting conditions (DTCC-ESC) is proposed in this paper for the first time. The DTCC-ESC mainly includes three parts: the engine segment active control, the feedforward and feedback control of engine starting, and the active damping feedback compensation control for the resonance of the transmission system.
(1) The engine segment active control, also known as the engine transient torque control before and after the engine ignition, is used to reduce the impact of the ERT on system shocks and vibrations in the engine startup from the perspective of the engine power source. Regarding this engine segment active control: before ignition, the inputs of control are the engine's intake valve close timing (IVCT) and initial crankshaft angle (ICA); after ignition, the inputs of control are fuel injection and engine ignition timing optimized by a genetic algorithm based on the engine working cycle. The output of the control is the engine output torque. (2) Feedforward and feedback control of engine starting is mainly realized by solving the motor torque distribution algorithm from Equation (14) through different engine demand angular accelerations before and after ignition. For the feedforward and feedback control of the engine startup, the inputs of control are a) the required angular acceleration of the engine based on the engine starting conditions before ignition and b) the required angular acceleration of the engine based on the engine optimal target speed after ignition. The outputs of control are the distributed torques of the two motors based on Equation (14). (3) The active damping feedback compensation control for the system resonance (ADFC-SR) is adopted to restrain the vibration of the transmission system. The inputs of the ADFC-SR are the estimated ring speed and torque, the estimated engine output shaft torque, and the estimated wheel speed. Its outputs are the active damping feedback compensation control torques from Equation (22).     According to the engine thermodynamic Equation (1), since the temperature and gas constant in the cylinder are little changed before ignition, the engine's instantaneous cylinder pressure P cyl,i (ϕ) is mainly related to the instantaneous cylinder volume V(ϕ) and air mass m(ϕ) in the cylinder. The intake valve close timing (IVCT) of the engine can change the initial volume V initial (ϕ) and gas mass m(ϕ) of the closed cylinder in the process of engine intake and compression, and the initial crankshaft angle (ICA) can change the initial volume V initial (ϕ) and gas mass m(ϕ) of the closed cylinder in the compression stage of the first working cycle of the engine startup. Therefore, the IVCT and ICA can directly change the cylinder pressure in the engine compression stage. Combining Equation (3), the relationship between engine ripple torque T ERT and IVCT and ICA can be obtained.
where θ ivct is the angle between IVCT and bottom dead center; θ ica is the ICA between the piston stop position (PSP) and the bottom dead center in the engine compression stage, as shown in Figure 9.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 12 of 2 change the cylinder pressure in the engine compression stage. Combining Equation (3) the relationship between engine ripple torque ERT T and IVCT and ICA can be obtained.
where ivct  is the angle between IVCT and bottom dead center; ica  is the ICA between the piston stop position (PSP) and the bottom dead center in the engine compression stage as shown in Figure 9. As shown in Figure 10, before engine ignition, when ica  is 85° and the cooling wa ter temperature is 25 °C, the engine-output torque decreases significantly with the in crease in IVCT. This is because the actual air volume in the cylinder decreases with the increase in IVCT, resulting in a significant reduction of engine cylinder pressure ripple torque ( ) , i P ivct ica T  . Therefore, the engine output ripple torques can be significantly re duced by delaying IVCT. In Figure 11, it is observed that a different ICA has an obviou impact on engine output ripple torque when ivct  is 150° and the cooling water tempera ture is 25 °C. When the ICA is larger (closer to the top dead center), the engine's outpu shaft ripple torque decreases significantly. This is mainly because the closer the PSP is to top dead center, the smaller the initial transient cylinder pressure of the first working cycl in the compression stage will be. Therefore, according to the influence of the IVCT and ICA on engine output torque the original IVCT and ICA are optimized by ivct  = 150deg and ica  = 165deg to minimiz the engine output ripple torque to reduce the impact of the ERT on the system shocks and vibrations before ignition.  As shown in Figure 10, before engine ignition, when θ ica is 85 • and the cooling water temperature is 25 • C, the engine-output torque decreases significantly with the increase in IVCT. This is because the actual air volume in the cylinder decreases with the increase in IVCT, resulting in a significant reduction of engine cylinder pressure ripple torque T P i (θ ivct , θ ica ). Therefore, the engine output ripple torques can be significantly reduced by delaying IVCT. In Figure 11, it is observed that a different ICA has an obvious impact on engine output ripple torque when θ ivct is 150 • and the cooling water temperature is 25 • C. When the ICA is larger (closer to the top dead center), the engine's output shaft ripple torque decreases significantly. This is mainly because the closer the PSP is to top dead center, the smaller the initial transient cylinder pressure of the first working cycle in the compression stage will be. change the cylinder pressure in the engine compression stage. Combining Equation (3 the relationship between engine ripple torque ERT T and IVCT and ICA can be obtained. where ivct  is the angle between IVCT and bottom dead center; ica  is the ICA betwee the piston stop position (PSP) and the bottom dead center in the engine compression stage as shown in Figure 9.  Figure 9. Intake valve close timing optimization before engine ignition.
As shown in Figure 10, before engine ignition, when ica  is 85° and the cooling wa ter temperature is 25 °C, the engine-output torque decreases significantly with the in crease in IVCT. This is because the actual air volume in the cylinder decreases with th increase in IVCT, resulting in a significant reduction of engine cylinder pressure rippl torque ( ) , i P ivct ica T  . Therefore, the engine output ripple torques can be significantly re duced by delaying IVCT. In Figure 11, it is observed that a different ICA has an obviou impact on engine output ripple torque when ivct  is 150° and the cooling water tempera ture is 25 °C. When the ICA is larger (closer to the top dead center), the engine's outpu shaft ripple torque decreases significantly. This is mainly because the closer the PSP is t top dead center, the smaller the initial transient cylinder pressure of the first working cycl in the compression stage will be. Therefore, according to the influence of the IVCT and ICA on engine output torque the original IVCT and ICA are optimized by ivct  = 150deg and ica  = 165deg to minimiz the engine output ripple torque to reduce the impact of the ERT on the system shocks an vibrations before ignition.

Controls of Engine Fuel Injection and Ignition Timing
According to the bench test results, the equivalent jerk based on an engine working cycle number is obtained as shown in Figure 12. Note that the vehicle jerk is the second derivative of the vehicle speed [38] and the engine crankshaft angle every 720 deg of rotation is an engine working cycle. It can be seen from this figure that the transmission system has obvious jerks during engine starting. For the jerks after ignition, it is found that the equivalent jerk of the first three working cycles (3rd, 4th and 5th) after ignition is greater than that of subsequent engine working cycles. Therefore, this part is mainly focused on reducing the shocks and vibrations of the first three working cycles after ignition. Moreover, as shown in Figure 12, we can see that the second derivative of engine speed, Eng_out  has a similar fluctuation pattern with an equivalent jerk, which shows that the  Since the running time of the first three working cycles after ignition is about 300 ms, it is necessary to realize the transient torque change by means of fast engine response control. Although the engine output torque can be greatly adjusted by the throttle control method, its control response time is longer. Compared to engine throttle control, fuel injection and ignition controls have a faster control effect on engine transient torque [39]. Therefore, controls of fuel injection and ignition timing based on an engine working cycle will be used to rapidly change engine output torque to reduce jerks of the powertrain system during initial engine combustion, especially the first three working cycles after engine ignition in stage 4. As shown in Figure 13, in this paper, according to the engine working cycle and cooling water temperature of engine, the engine fuel injection  Therefore, according to the influence of the IVCT and ICA on engine output torque, the original IVCT and ICA are optimized by θ ivct = 150 deg and θ ica = 165 deg to minimize the engine output ripple torque to reduce the impact of the ERT on the system shocks and vibrations before ignition.

Controls of Engine Fuel Injection and Ignition Timing
According to the bench test results, the equivalent jerk based on an engine working cycle number is obtained as shown in Figure 12. Note that the vehicle jerk is the second derivative of the vehicle speed [38] and the engine crankshaft angle every 720 deg of rotation is an engine working cycle. It can be seen from this figure that the transmission system has obvious jerks during engine starting. For the jerks after ignition, it is found that the equivalent jerk of the first three working cycles (3rd, 4th and 5th) after ignition is greater than that of subsequent engine working cycles. Therefore, this part is mainly focused on reducing the shocks and vibrations of the first three working cycles after ignition. Moreover, as shown in Figure 12, we can see that the second derivative of engine speed, ... θ Eng_out has a similar fluctuation pattern with an equivalent jerk, which shows that the ... θ Eng_out can indirectly reflect the fluctuation phenomenon of a vehicle jerk. This provides a basis for choosing an objective function F GA1 optimized by a genetic algorithm.
Since the running time of the first three working cycles after ignition is about 300 ms, it is necessary to realize the transient torque change by means of fast engine response control. Although the engine output torque can be greatly adjusted by the throttle control method, its control response time is longer. Compared to engine throttle control, fuel injection and ignition controls have a faster control effect on engine transient torque [39]. Therefore, controls of fuel injection and ignition timing based on an engine working cycle will be used to rapidly change engine output torque to reduce jerks of the powertrain system during initial engine combustion, especially the first three working cycles after engine ignition in stage 4.

Controls of Engine Fuel Injection and Ignition Timing
According to the bench test results, the equivalent jerk based on an engine working cycle number is obtained as shown in Figure 12. Note that the vehicle jerk is the second derivative of the vehicle speed [38] and the engine crankshaft angle every 720 deg of rotation is an engine working cycle. It can be seen from this figure that the transmission system has obvious jerks during engine starting. For the jerks after ignition, it is found that the equivalent jerk of the first three working cycles (3rd, 4th and 5th) after ignition is greater than that of subsequent engine working cycles. Therefore, this part is mainly focused on reducing the shocks and vibrations of the first three working cycles after ignition. Moreover, as shown in Figure 12, we can see that the second derivative of engine speed, Eng_out  has a similar fluctuation pattern with an equivalent jerk, which shows that the Since the running time of the first three working cycles after ignition is about 300 ms, it is necessary to realize the transient torque change by means of fast engine response control. Although the engine output torque can be greatly adjusted by the throttle control method, its control response time is longer. Compared to engine throttle control, fuel injection and ignition controls have a faster control effect on engine transient torque [39]. Therefore, controls of fuel injection and ignition timing based on an engine working cycle will be used to rapidly change engine output torque to reduce jerks of the powertrain system during initial engine combustion, especially the first three working cycles after engine ignition in stage 4. As shown in Figure 13, in this paper, according to the engine working cycle and cooling water temperature of engine, the engine fuel injection _ iv cycj F and ignition timing _ it cycj  of the first three working cycles after engine ignition are optimized by using a genetic As shown in Figure 13, in this paper, according to the engine working cycle and cooling water temperature of engine, the engine fuel injection F iv_cycj and ignition timing θ it_cycj of the first three working cycles after engine ignition are optimized by using a genetic algorithm [40]. Since the main purpose of this paper is to reduce shocks in mode switching, which that includes engine starting, the effect of emissions will not be considered as a subject of interest.
The injection volume F iv_cycj of the first three working cycles after ignition is expressed as follows: F iv_cycj = C f rc_cycj * F b f _cycj (16) where C f rc_cycj and F b f _cycj are the fuel-ratio coefficient and basic fuel injection at the jth working cycle (j = 3, 4, 5), respectively. In the initial combustion of an engine startup, a too large or too small ignition angle will lead to a fuel injection combustion misfire [41]. Therefore, in practice, the engine ignition angle and fuel injection need to be limited to a certain range to ensure a smooth start. The combustion misfire constraints of the engine from engine experimental data are as follows.
where the ignition timing at the top dead point is defined as zero. The negative θ it_cycj is the ignition advance angle, and the positive θ it_cycj is the ignition delay angle.
Since the second derivative of engine speed, ... θ Eng_out can reflect the fluctuation of a vehicle jerk, the absolute value of ... θ Eng_out is used as the fitness function F GA1 of the genetic algorithm:  Figure 13. Flow diagram of the fuel injection and ignition timing optimization.
The optimization is achieved in steps illustrated in the Figure13, which are explained as follows: Step 1: In accordance with the established powertrain model, the engine system enters the working cycle j (j = 3, 4, 5) after receiving the ignition command. The optimization is achieved in steps illustrated in the Figure 13, which are explained as follows: Step 1: In accordance with the established powertrain model, the engine system enters the working cycle j (j = 3, 4, 5) after receiving the ignition command.
Step 2: Under the engine coolant temperature and combustion misfire constraints of the engine, the genetic algorithm [40] automatically initializes the fuel ratio coefficient C f rc_cycj and ignition timing θ it_cycj to obtain the initial transient cylinder pressure, and then calculates the engine ripple torque.
Step 3: Using the powertrain simulation model with the obtained ERT, the fitness function F GA1 is then calculated.
Step 4: The termination condition of genetic algorithm iteration is then checked. If it is completed, the optimization of the first three working cycles after ignition will be carried out in sequence according to the above three steps.
When the cooling water temperature is set to 25 • C, results before and after optimization of the first three working cycles after ignition are obtained by the genetic algorithm, as shown in Figures 14 and 15. From Figure 14a, we can see that in the 4th and 5th engine working cycles, the injection volumes are reduced to 40% of the basic fuel injection. In Figure 14b, the ignition timings in these three working cycles are all delayed. As observed in Figure 15, the amplitude of ... θ Eng_out is reduced significantly by the optimization, which can reflect the reduction of the vehicle jerks. In addition, the fuel injection volume of the third working cycle (that is the first working cycle after ignition) in Figure 14a is 20% higher than it was before optimization. This is mainly because the firing reliability can be improved by properly increasing the injection volume of the first working cycle after a coldstart ignition, and enough oil film can be formed in the engine intake port. However, due to the delay of ignition timing, the ... θ Eng_out is still less than the value before optimization in Figure 15 in the 3rd engine working cycle.

Feedforward and Feedback Controls of Engine Starting
According to Equation (14), the required acceleration of the engine Eng_req  is an important control variable in the engine starting process. As shown in Figure 8, feedforward and feedback controls of engine starting are mainly realized by solving the motor torque distribution algorithm, Equation (14), based on different demand-angular acceleration before and after engine ignition.

Feedforward Control before Ignition
As shown in Figure 8, to start the engine reliably and smoothly before ignition, the

Feedforward and Feedback Controls of Engine Starting
According to Equation (14), the required acceleration of the engine Eng_req  is an important control variable in the engine starting process. As shown in Figure 8, feedforward and feedback controls of engine starting are mainly realized by solving the motor torque distribution algorithm, Equation (14), based on different demand-angular acceleration before and after engine ignition.

Feedforward Control before Ignition
As shown in Figure 8, to start the engine reliably and smoothly before ignition, the

Feedforward and Feedback Controls of Engine Starting
According to Equation (14), the required acceleration of the engine .. θ Eng_req is an important control variable in the engine starting process. As shown in Figure 8, feedforward and feedback controls of engine starting are mainly realized by solving the motor torque distribution algorithm, Equation (14), based on different demand-angular acceleration before and after engine ignition.

Feedforward Control before Ignition
As shown in Figure 8, to start the engine reliably and smoothly before ignition, the feedforward control method before ignition was adopted. Firstly, the "S-curve" function f (x, a, b, c) from Equation (19) was taken as the required acceleration of the engine, .. θ Eng_des .  Then, the f (x, a, b, c) was brought into Equation (14) to obtain the two motor-distribution torques, which are required by the feedforward control before ignition.
In the feedforward control, under the given cooling water temperature and different engine starting conditions, taking the absolute value of vehicle jerk as the fitness function F GA2 , the genetic algorithm was used to optimize a, b and c to obtain the optimal "S-curve" function (the optimal engine demand-angular acceleration) under different engine starting conditions before engine ignition.
The S-curve function f (x, a, b, c) is defined as follows. ..
where a, b and c are the parameters of the ..
θ Eng_des . The jerk absolute value was used as the fitness function F GA2 of the genetic algorithm: where j erk is the vehicle jerk.

Feedback Control after Ignition
In order to make the actual engine speed follow the engine optimal target speed in time after ignition, the engine angular acceleration required by feedback control will be obtained by a PID controller (p = 5, I = 1 and d = 0.4). As shown in Figure 8, the inputs of the PID controller were engine optimal speed and the estimated carrier speed. The estimated carrier speedˆ. θ c_axle is expressed as follows: where . θ MG1 and . θ MG2 are the speeds of MG1 and MG2, respectively.

Active Damping Feedback Compensation Torques
The ADFC torques of the two motors (T MG1_Damp and T MG2_Damp ) can be obtained by the equivalent lever method [40] in Figure 6, which compensates for the torque fluctuations transmitted to the wheels. ..
where T R_ad is the feedback damping torque and I R is the ring inertia. In addition,T c_g_est_f was obtained by Kalman filtering of the estimated carrier torqueT c_g_est ; .. θ MG1 is the MG1 acceleration; and .. θ MG2 is the MG2 acceleration. The formula of the estimated carrier torqueT c_g_est (engine output shaft torque) is obtained as follows: To reduce the resonance of the transmission system itself and the torque fluctuations transmitted to the wheel, the feedback damping torque T R_ad is expressed as whereT R_g_est_ f is obtained by Kalman filtering of the estimated ring torqueT R_g_est ; f d is the feedback gain damping coefficient;Ŵ ob is the estimated wheel speed;ˆ. θ R is the estimated speed at ring gear; and T PID (Ŵ ob i 0 ,ˆ. θ R ) is the adjustment of the speed error at the gear ring, which is from a PID controller whose inputs areŴ ob i 0 andˆ. θ R . The estimated ring torqueT R_g_est is determined as follows: The estimated speed at ring gear,ˆ. θ R , can be obtained as follows:

Wheel Speed Estimation
A wheel-speed observer is designed to obtain an accurate, real-time wheel speedŴ ob . In the following formula, X and U are selected as the state variables and control variables, respectively.
The discretization state equation is expressed as follows: with M = I L 0 0 I R , N = −1 0 0 1 , C = −C TI C TI C TI −C TI , and K = −k TI k TI k TI −k TI , whereX is the estimation state variables; I is the unit matrix; and k is the sample number.

Determination of Feedback Gain Damping Coefficient
Based on the simplified spring-mass model between point A and point B in Figure 8, the system transfer function G R (s) from point A to point B can be obtained as follows: where I R is the equivalent ring inertia; I L is the equivalent moment of inertia at point B in Figure 8; and s is the Laplace operator. Figure 16 is the Bode diagram of the system before and after the damping feedback control. It shows the amplitude-frequency and phase-frequency characteristics of the transfer function G R (s). From the solid red line before the damping feedback control, it could be seen that the transmission system resonates at 9.86 Hz. To reduce the resonance of the powertrain system further, the speed error of the output shaft was controlled by the damping feedback method, which is defined as ' ' '  Figure 17 is a zero-pole distribution diagram based on the root locus, which is used to allocate the poles of the system. The transfer function ' ( ) R Gs has three poles before the damping feedback control, one of which is at the origin of Figure 17; the other two are a pair of conjugate complex poles with fd = 0. Under the damping feedback closed-loop control, with the increase of fd, the two poles without origin will gradually approach the real axis and become two double roots at point A in Figure 17. When fd continues to increase, one pole of the system moves to the left and the other pole moves to the right, and the response time of the system will increase. Based on this, in this paper, the two poles are located at the same point A in Figure 17, and the feedback gain damping coefficient fd is 5.6228. In addition, the dotted line after the damping feedback control in Figure 16 is the Bode diagram result when fd = 5.6228. The feedback gain-damping control optimized by the zero-pole distribution method can reduce the resonance peak of the transmission system better. 80 Figure 16. Bode diagram of system before and after damping feedback control.
By introducing the feedback gain damping coefficient f d into Equation (32), the transmission system transfer function from point A to point B under closed-loop control can be obtained again as follows: Figure 17 is a zero-pole distribution diagram based on the root locus, which is used to allocate the poles of the system. The transfer function G R (s) has three poles before the damping feedback control, one of which is at the origin of Figure 17; the other two are a pair of conjugate complex poles with f d = 0. Under the damping feedback closed-loop control, with the increase of f d , the two poles without origin will gradually approach the real axis and become two double roots at point A in Figure 17. When f d continues to increase, one pole of the system moves to the left and the other pole moves to the right, and the response time of the system will increase. Based on this, in this paper, the two poles are located at the same point A in Figure 17, and the feedback gain damping coefficient f d is 5.6228. In addition, the dotted line after the damping feedback control in Figure 16 is the Bode diagram result when f d = 5.6228. The feedback gain-damping control optimized by the zero-pole distribution method can reduce the resonance peak of the transmission system better.
pair of conjugate complex poles with fd = 0. Under the damping feedback closed-loop control, with the increase of fd, the two poles without origin will gradually approach the real axis and become two double roots at point A in Figure 17. When fd continues to increase, one pole of the system moves to the left and the other pole moves to the right, and the response time of the system will increase. Based on this, in this paper, the two poles are located at the same point A in Figure 17, and the feedback gain damping coefficient fd is 5.6228. In addition, the dotted line after the damping feedback control in Figure 16 is the Bode diagram result when fd = 5.6228. The feedback gain-damping control optimized by the zero-pole distribution method can reduce the resonance peak of the transmission system better.

Simulation Conditions
To verify the effectiveness of the proposed DTCC-ESC method more comprehensively, six typical engine starting conditions (Case A, Case B, . . . , Case F), as shown in Figure 18 and Table 2, were set as simulation conditions during mode switching that included engine starting. Figure 18 shows the required ring torque converted from the acceleration pedal position (APP) and vehicle speed in EV mode. Table 2 further shows other engine-starting conditions, which include engine coolant temperature, battery temperature, and battery SOC. Note that the influence of SOC, temperature, and other factors on the dynamic response characteristics of the vehicle during engine starting can be referred to the literature [42]. To verify the effectiveness of the proposed DTCC-ESC method more comprehensively, six typical engine starting conditions (Case A, Case B, …, Case F), as shown in Figure 18 and Table 3, were set as simulation conditions during mode switching that included engine starting. Figure 18 shows the required ring torque converted from the acceleration pedal position (APP) and vehicle speed in EV mode. Table 3 further shows other engine-starting conditions, which include engine coolant temperature, battery temperature, and battery SOC. Note that the influence of SOC, temperature, and other factors on the dynamic response characteristics of the vehicle during engine starting can be referred to the literature [42].

Comparative Analysis of Different Control Methods
In order to verify the effectiveness of the proposed DTCC-ESC method, as shown in Table 3, the active damping compensation control (ADCC) method, the active damping compensation control with engine active control (ADCC-EAC) method, and the DTCC-ESC method proposed in this paper were selected for simulation and comparative analysis on the established compound power-split PHEV model. In addition, taking point D of the engine starting condition as an example, the simulation results were analyzed in detail, and then six typical engine-starting condition points in Figure 18 were taken as examples to analyze the vehicle jerks under different engine starting conditions.  Figure 19 shows the results of engine-demand angular acceleration before and after optimization under engine-starting condition point D. It can be seen that the amplitude of the optimized demand-angle acceleration of the engine is smaller, and the change trend is gentler than that before optimization. This is mainly because the optimized engine demand-angular acceleration reduces the motor-distribution torques of Equation (14), so as to prevent the motor-output torques from being limited by battery charge and discharge power.  Figure 19. Optimized engine demand-angular acceleration before ignition under case D.
The torque results of the power source obtained by the three control methods during stages 3 and 4 under case D are shown in Figure 20, which shows that when the ADCC method started the engine the torque of the engine output shaft fluctuated clearly. The MG1 torque is basically at the maximum peak value. The MG2 torque shows the fluctuation phenomenon of torque compensation, but some of the torque reached the maximum peak value of MG2, which cannot completely compensate or offset the fluctuating torque of the transmission system. From Figure 20b, as the ADCC-EAC method was used to start the engine, the torque fluctuation of the engine output shaft were significantly reduced. This shows that the engine active control can effectively reduce the torque fluctuation of the engine output shaft. However, the torque of MG1 is still at the maximum peak. Some of the MG2 torque still reached the maximum peak value of MG2. It can be seen from Figure 20c that when the proposed DTCC-ESC method started the engine, the MG1 torque gradually increased to its peak value and then decreased. The pulsating torque of the engine output shaft was effectively reduced by the segmented active control method of engine. The MG2 torque did not exceed the limit of its own maximum torque, which means that the proposed DTCC-ESC method can compensate or offset the pulsating torque of the transmission system better. In addition, it was further found that the compensation fluctuating torque of MG2 had a highly symmetrical trend with the corresponding fluctuating torque of the engine output shaft. This showed that the fluctuating torque of the transmission system was mainly compensated for by MG2, while MG1 was mainly used for starting the engine rather than compensating for the fluctuating torque of the system.
The jerking results of the three control methods during the engine-starting process under case D are shown in Figure 21 and Table 4. From Figure 21, when the engine was started by the three control methods, respectively, the vehicle jerk fluctuation generated by the ADCC method was the largest, while that generated by DTCC-ESC method was the smallest. As shown in Table 4, the maximum value of absolute vehicle jerks, J a-max , reached 114.5 m/s 3 using the ADCC method. It is difficult for drivers to accept such a large jerk. The J a-max value produced by the DTCC-ESC method is 34.58 m/s 3 . Compared with the result produced by the ADCC method, the jerk produced by DTCC-ESC improved by 69.8%, which showed that the proposed engine active control method can reduce system shocks. When the DTCC-ESC method was used to start the engine, vehicle jerks were successfully reduced to 9.94 m/s 3 . This result is 91.32% better than that of the ADCC method, and 71.26% better than that of the ADCC-EAC method. This showed that the DTCC-ESC method in this paper can improve the driving comfort during engine starting.

Jerk Analysis under Different Engine Starting Conditions
Since the vehicle jerk produced by the ADCC method was significantly less than 10 m/s 3 , which can be accepted by drivers (as shown in Figure 23) under cases A and B, the demand-angular acceleration of the engine before ignition was optimized for the other four cases (C, D, E and F) as shown in Figure 22. It can be seen from Figure 22 that different engine demand-angular accelerations are presented under different enginestarting conditions. It is mainly the influence of the engine-starting conditions that prevents the motor output torque from being limited by battery power and gives as much maximum torque as possible. Figure 23 shows J a-max values produced by the three control methods under six engine starting conditions. The J a-max values with the ADCC method in cases C, D, E and F are greater than 85 m/s 3 . It is difficult for drivers and passengers to accept such large vehicle jerks. Compared with the ADCC method, although the ADCC-EAC method can effectively reduce vehicle jerks in cases A, B, C and F, the ADCC-EAC method could not effectively reduce jerks under cases D and E. Moreover, the J a-max value in the case E is 60.8m/s 3 , which is also difficult to accept. For the J a-max values generated by the DTCC-ESC method, it could be seen that the proposed DTCC-ESC method can effectively controlled the J a-max values within 10m/s 3 in the six cases of engine starting. Compared with the ADCC-EAC method, the J a-max value of DTCC-ESC method at engine starting point E improved by 97.1%.
Therefore, the DTCC-ESC method proposed in this paper can effectively reduce jerks and improve driving comfort under different engine-starting conditions.  Since the vehicle jerk produced by the ADCC method was significantly less than 10m/s 3 , which can be accepted by drivers (as shown in Figure 23) under cases A and B, the demand-angular acceleration of the engine before ignition was optimized for the other four cases (C, D, E and F) as shown in Figure 22. It can be seen from Figure 22 that different engine demand-angular accelerations are presented under different engine-starting conditions. It is mainly the influence of the engine-starting conditions that prevents the motor output torque from being limited by battery power and gives as much maximum torque as possible. Figure 23 shows Ja-max values produced by the three control methods under six engine starting conditions. The Ja-max values with the ADCC method in cases C, D, E and F are greater than 85 m/s 3 . It is difficult for drivers and passengers to accept such large vehicle jerks. Compared with the ADCC method, although the ADCC-EAC method can effectively reduce vehicle jerks in cases A, B, C and F, the ADCC-EAC method could not effectively reduce jerks under cases D and E. Moreover, the Ja-max value in the case E is 60.8m/s 3 , which is also difficult to accept. For the Ja-max values generated by the DTCC-ESC method, it could be seen that the proposed DTCC-ESC method can effectively controlled the Ja-max values within 10m/s 3 in the six cases of engine starting. Compared with the ADCC-EAC method, the Ja-max value of DTCC-ESC method at engine starting point E improved by 97.1%.
Therefore, the DTCC-ESC method proposed in this paper can effectively reduce jerks and improve driving comfort under different engine-starting conditions.

Conclusions
A compound power-split power transmission system model was first established and verified by bench experiments. The engine starting process was analyzed based on an equivalent lever method. Then, the problem of system jerks and vibrations caused by insufficient motor torque compensation under some engine starting conditions was analyzed. To solve this problem, this paper proposed a dynamic torque coordinated control for engine starting conditions (DTCC-ESC). Finally, three control methods (ADCC, ADCC-EAC, and DTCC-ESC) were selected for simulation and comparative analysis on the established hybrid electric vehicle model. The results showed that the proposed DTCC-ESC strategy can reduce system shocks under different engine start-up conditions during the mode switching process that includes engine start-up. Compared with the other two methods, the DTCC-ESC method had a better control effect.
In future research, the controller signal network delay can be considered to optimize the dynamic coordination control method to further reduce the system shock and vibration.

Conflicts of Interest:
The authors declare no conflict of interest.

ADCC
Active damping compensation control ADCC-EAC Active damping compensation control with engine active control ADFC-SR Active damping feedback compensation control considering the system resonance APP