Study of the Impact of E-Machine in Hybrid Dual Clutch Transmission Powertrain

: Recent studies show that hybrid powertrains with dual clutch transmission are becoming more common. This is mainly inﬂuenced by the fact that they have higher efﬁciency, allow power-on shifting, and are more compact. However, the electric motor attached to the primary shaft of the transmission might represent some critical issues in statics and dynamics due to additional torque, inertia, and working irregularities. This paper aims to investigate the inﬂuence of the electric motor integration in the transmission on the reliability of the overall system. The focus is on the investigation of potential over stresses of the ICE bearings due to the bending of the gearbox primary shaft. To achieve the aim, the static and the rotor dynamic analysis of the inner and outer primary shafts is performed. The numerical analyses show the potential critical issues that can arise from the excessive manufacturing and assembling tolerances, the interaction of the e-machine torque irregularities and inertia with the bending dynamics of the primary shafts.


Introduction
New ambitious targets on the amount of CO 2 emissions set by European Commission [1,2] require implementation of even more efficient technologies from car manufacturers. Since the last decade of the 20th century, a wide variety of different technologies and policies aiming the reduction of vehicle fuel consumption and emissions at different vehicle and infrastructure level are widely implemented [3,4]. Most regulations in the past were focused on the reduction of the pollutants. The real motivation to the hybridization is the more stringent limit to the emissions. This reduction is not achieved by simple ICE tuning and pushes towards hybrids [5].
Hybridized powertrains can be generally of a dedicated or an add-on type [6,7]. Dedicated Hybrid Transmissions (DHT) use custom designed transmissions that can be mostly of a planetary gear train type [6,7]. Add-on types are obtained by adding an Electric Motor (EM) to conventional powertrain at different positions. Based on the position of the electric motor they can be of P0 to P5 types [7,8]. Add-on type hybrid transmissions are gaining great interest of the manufacturers as they can be developed as hybridization kits [9]. In this regard, most hybridization kits are obtained by integrated design of the Dual Clutch Transmissions (DCT) and the electric machine to obtain a Hybrid Dual Clutch Transmissions (HDT). The choice of the DCT as transmission is largely driven by the fact that they have high efficiency, simple construction, lower inertia in case of the adoption of wet clutches and smooth gear shifting [9,10]. Furthermore, DCT can allow power-on gear shifting. All this makes it the architecture of choice in many hybrid solutions [11][12][13]. The use of DCT as a part of HEV (i.e., HDT) allows the Internal Combustion Engine (ICE) The aim of the present paper is to study the influence of the above-mentioned factors on the potential damage of the ICE and DMF bearings and bushings. The novelty of the work lies in the study of the influence of the electromagnetic and inertia torque of the electric machine on the flexural deflection of the gear shafts and on the modification of the natural frequencies in the rotor dynamic case. The coupling between torsional and flexural behavior highlights the potential static and dynamic issues that might arise during the integration of ICE with HDT module. The emerging trend in automotive industry is integration of HDT from external suppliers with the in-house ICE. Therefore, the practical utility of the study is to understand the possible operational issues of the integrated HDT and ICE module.
The HEV powertrain configuration analyzed in this study is shown in Figure 1. The power sources are an ICE, which is connected to a DMF to reduce the engine torque's cyclic irregularity and an EM, which is connected to the gearset with even gear numbers. Both are combined with a DCT consisting of two separate friction clutches C1 and C2 as well as two sub-gearsets. its misalignment based on numerical analysis. Additionally, there is no analysis on bending-torsional coupling of the HDT from static and dynamic point of view. The aim of the present paper is to study the influence of the above-mentioned factors on the potential damage of the ICE and DMF bearings and bushings. The novelty of the work lies in the study of the influence of the electromagnetic and inertia torque of the electric machine on the flexural deflection of the gear shafts and on the modification of the natural frequencies in the rotor dynamic case. The coupling between torsional and flexural behavior highlights the potential static and dynamic issues that might arise during the integration of ICE with HDT module. The emerging trend in automotive industry is integration of HDT from external suppliers with the in-house ICE. Therefore, the practical utility of the study is to understand the possible operational issues of the integrated HDT and ICE module.
The HEV powertrain configuration analyzed in this study is shown in Figure 1. The power sources are an ICE, which is connected to a DMF to reduce the engine torque's cyclic irregularity and an EM, which is connected to the gearset with even gear numbers. Both are combined with a DCT consisting of two separate friction clutches C1 and C2 as well as two sub-gearsets. The first part of the paper deals with static analysis where different operating modes of the hybrid DCT transmission are considered and the static deflection due to the torque of the electric motor in different operating modes is calculated using a FEM.
In the second part the rotor dynamic analysis is performed using DYNROT, a FEM code to study the dynamics of the rotating parts (Campbell diagram and mode shapes) including the effect of EM and ICE torque irregularities on the displacement at the input splines of the inner primary shaft. Finally, the results are discussed, and the values of deflections are compared with those indicated in the service manuals for normal operation of the DMF to highlight potential critical issues.

Methodology
The present section describes the methodology used to carry out the static and rotor dynamic analysis on the HDT.
The general layout of the HDT that will be used for the static and rotor dynamic analysis is shown in Figure 2. The ICE is connected to the DMF, whose output is coupled to the two input splines of dual clutch system C1 and C2. The two clutches transmit the torque to the inner (IPS) and outer primary shaft (OPS) of the drivetrain including even and odd gears, respectively. The EM is installed above the primary shaft, its pinion (G1) transmits torque to the outer primary shaft by means of an idler (G2). The inner and outer primary shafts are supported by ball bearings BB1 and BB2 while the inter shaft needle bearings NB1 and NB2 are installed between the two primary shafts. The first part of the paper deals with static analysis where different operating modes of the hybrid DCT transmission are considered and the static deflection due to the torque of the electric motor in different operating modes is calculated using a FEM.
In the second part the rotor dynamic analysis is performed using DYNROT, a FEM code to study the dynamics of the rotating parts (Campbell diagram and mode shapes) including the effect of EM and ICE torque irregularities on the displacement at the input splines of the inner primary shaft. Finally, the results are discussed, and the values of deflections are compared with those indicated in the service manuals for normal operation of the DMF to highlight potential critical issues.

Methodology
The present section describes the methodology used to carry out the static and rotor dynamic analysis on the HDT.
The general layout of the HDT that will be used for the static and rotor dynamic analysis is shown in Figure 2. The ICE is connected to the DMF, whose output is coupled to the two input splines of dual clutch system C1 and C2. The two clutches transmit the torque to the inner (IPS) and outer primary shaft (OPS) of the drivetrain including even and odd gears, respectively. The EM is installed above the primary shaft, its pinion (G1) transmits torque to the outer primary shaft by means of an idler (G2). The inner and outer primary shafts are supported by ball bearings BB1 and BB2 while the inter shaft needle bearings NB1 and NB2 are installed between the two primary shafts. The reverse gear as well as the final gear connected to the differential are not represented in Figure 2.

Static Analysis
The aim of the present section is to evaluate the deflections of the HDT primary shaft due to the transmitted torque, including the contribution of the EM. Figure 3 shows the layout of the FEM model used for the static analysis. Starting from the top, the EM is connected to the pinion gear (G1) that delivers the torque to the idler (G2). This torque is then transmitted to the gear (G3) on the outer primary shaft (OPS) and finally to the secondary shafts (not shown). The outer and inner primary shafts are supported using the ball bearings BB1 and BB2. While needle bearings NB1 and NB2 are installed between the inner and outer primary shafts. As the main possible criticality could arise on the primary shaft due to the coupling between the HDT and ICE, therefore, this shaft is considered in the study. Furthermore, it is the one that exchanges the radial load between the two modules. The reverse gear as well as the final gear connected to the differential are not represented in Figure 2.

Static Analysis
The aim of the present section is to evaluate the deflections of the HDT primary shaft due to the transmitted torque, including the contribution of the EM. Figure 3 shows the layout of the FEM model used for the static analysis. Starting from the top, the EM is connected to the pinion gear (G1) that delivers the torque to the idler (G2). This torque is then transmitted to the gear (G3) on the outer primary shaft (OPS) and finally to the secondary shafts (not shown). The outer and inner primary shafts are supported using the ball bearings BB1 and BB2. While needle bearings NB1 and NB2 are installed between the inner and outer primary shafts. As the main possible criticality could arise on the primary shaft due to the coupling between the HDT and ICE, therefore, this shaft is considered in the study. Furthermore, it is the one that exchanges the radial load between the two modules. The reverse gear as well as the final gear connected to the differential are not represented in Figure 2.

Static Analysis
The aim of the present section is to evaluate the deflections of the HDT primary shaft due to the transmitted torque, including the contribution of the EM. Figure 3 shows the layout of the FEM model used for the static analysis. Starting from the top, the EM is connected to the pinion gear (G1) that delivers the torque to the idler (G2). This torque is then transmitted to the gear (G3) on the outer primary shaft (OPS) and finally to the secondary shafts (not shown). The outer and inner primary shafts are supported using the ball bearings BB1 and BB2. While needle bearings NB1 and NB2 are installed between the inner and outer primary shafts. As the main possible criticality could arise on the primary shaft due to the coupling between the HDT and ICE, therefore, this shaft is considered in the study. Furthermore, it is the one that exchanges the radial load between the two modules.  The FEM modeling has been carried out in SOLIDWORKS simulation software. The shafts have been modelled as solid parts with material properties of AISI 4340. The ball bearings are modeled as stiffness element reacting just in radial and axial direction. The needle bearings are modeled as inter-shaft stiffness element acting just in the radial direction. The bearings stiffnesses are summarized in Table 1. The radial stiffness of the ball bearing was estimated using the Hertz theory [33,34] while the needle bearings stiffness using the Palmgren theory [35]. With reference to Figure 4, the operating mode of the EM (b-boost, c-recuperation) defines the directions of the forces acting on the driving and driven helical gear. In Figure 4, the EM provides a torque of 55 Nm rotating in clockwise direction seen from left hand side. This is then delivered to the idler (G2) which engages to the helical gear (G3) on the outer primary shaft. The forces exchanged between the gears are calculated [36] as follows: The forces on the gear G3 used for the FEM simulation are calculated as follows: In boost case: In recuperation case: While the axial forces used for FEM simulation remain the ones depicted in Figure  4a.

Scenario 2: 2nd Gear in Boost and Recuperation
Similar to the 1st scenario, the EM works in both boost and recuperation mode even if the torque path is different, as the EM delivers the torque to the OPS using the gears G1-G2-G3 ( Figure 3) and it is finally delivered to the 2nd gear on the secondary shaft 2 using the helical gear (G4). Therefore, there are 2 helical gears on the OPS transmitting torque as depicted in Figure 5a. In Figure 5b,c the force contribution coming from the G3 and G4 gears is shown in boost and recuperation mode, respectively. G3 gear forces work in B-B plane and are shown with solid lines, while the G4 gear forces work in A-A plane and are shown with dashed lines. The tangential force (F t1 ) at the driven gear (G3) is given by: where, D G3 indicates the pitch diameter of the helical gear G3, T EM is EM Torque (56 Nm taken as 70% of the Maximum torque, that is equal to 80 Nm [9]) and τ G13 is gear ratio between the EM pinion gear (G1) and the outer primary shaft gear (G3). The radial (F r1 ) and axial (F a1 ) forces at the driven gear (G3) are given by: where, β is gear helix angle (25 • ) and φ is helical gear normal pressure angle (20 • ).
In Figure 4 only the forces acting on the outer primary shaft have been depicted for simplification purposes. These forces are the ones generated due to the EM torque and will cause the deflection in the shaft. In Figure 4b is considered the EM in boost mode where the torque is being transferred from the G2 (driver) gear to the G3 (driven) gear, resulting in the tangential, radial and axial forces depicted with red arrows (A-A plane). This torque is then delivered from G3 gear to the 4th gear on the secondary shaft 1 resulting in the forces F t2 and F r2 . The resultant forces acting on the G3 gear on the outer primary shaft are shown with black arrows. Similarly, the EM in recuperation mode, as shown in Figure 4c, will receive the torque from G2 gear that is delivered from the 4th gear on the secondary shaft 1 to the gear G3. Hence, the direction of tangential and axial forces will be opposite to the boost case.
The forces on the gear G3 used for the FEM simulation are calculated as follows: In boost case: In recuperation case: While the axial forces used for FEM simulation remain the ones depicted in Figure 4a.

Scenario 2: 2nd Gear in Boost and Recuperation
Similar to the 1st scenario, the EM works in both boost and recuperation mode even if the torque path is different, as the EM delivers the torque to the OPS using the gears G1-G2-G3 ( Figure 3) and it is finally delivered to the 2nd gear on the secondary shaft 2 using the helical gear (G4). Therefore, there are 2 helical gears on the OPS transmitting torque as depicted in Figure 5a. In Figure 5b,c the force contribution coming from the G3 and G4 gears is shown in boost and recuperation mode, respectively. G3 gear forces work in B-B plane and are shown with solid lines, while the G4 gear forces work in A-A plane and are shown with dashed lines.
The forces on the gear G3 and G4 that will be used for the FEM simulation are calculated as follows: In boost case: In recuperation case: While the axial forces used for FEM simulation remain the ones depicted in Figure 5.
While the axial forces used for FEM simulation remain the ones depicted in Figure 5.

Scenario 3: ICE + EM in 4th Gear Boost
In this scenario torque from both the power sources (ICE and EM) is considered, with the EM operating in boost mode only. In Figure 6a the forces due to the torque from ICE and EM are shown. The EM torque is introducing the forces on gear G3 similar to the 4th gear boost mode, in addition there are forces due to torque of ICE that will act on the gear G3 depicted using the light blue lines. As a result, the forces experienced by gear G3 will be higher compared to the other two scenarios and is considered as the scenario with maximum torque that can be transferred by the transmission. The total torque that can be transmitted by the HDT is taken as 80% of the maximum torque mentioned in the manufactures catalogue [9]. In Figure 6b the force acting in the A-A plane are shown.
The forces on the gear G3 that will be used for the FEM simulation are calculated as follows: In boost case: For both cases, the axial forces used for FEM simulation remain the ones depicted in Figure 6a. In this scenario torque from both the power sources (ICE and EM) is considered, with the EM operating in boost mode only. In Figure 6a the forces due to the torque from ICE and EM are shown. The EM torque is introducing the forces on gear G3 similar to the 4th gear boost mode, in addition there are forces due to torque of ICE that will act on the gear G3 depicted using the light blue lines. As a result, the forces experienced by gear G3 will be higher compared to the other two scenarios and is considered as the scenario with maximum torque that can be transferred by the transmission. The total torque that can be transmitted by the HDT is taken as 80% of the maximum torque mentioned in the manufactures catalogue [9]. In Figure 6b the force acting in the A-A plane are shown.
The forces on the gear G3 that will be used for the FEM simulation are calculated as follows: In boost case: For both cases, the axial forces used for FEM simulation remain the ones depicted in Figure 6a.

Numerical Simulation Parameters
The mesh used for the FEM analysis is curvature-based tetrahedron mesh as shown in Figure 7a,b. The mesh density is chosen considering the computational cost and the convergence of the deflection results as shown in convergence plot in Figure 7c. It can be seen that the displacement results saturate with the increase of the mesh density. Hence, mesh with maximum element size of 1.5 mm is selected. The quality of the mesh elements is also an important aspect and could influence the results. Therefore, the FE mesh elements were run through an element quality check based on the aspect and Jacobian ratio. 99.8% of the elements fall under the aspect ratio of 3, with the maximum aspect ratio reaching the value of 6. Based on the Jacobian ratio results, no distorted elements were detected.

Numerical Simulation Parameters
The mesh used for the FEM analysis is curvature-based tetrahedron mesh as shown in Figure 7a,b. The mesh density is chosen considering the computational cost and the convergence of the deflection results as shown in convergence plot in Figure 7c. It can be seen that the displacement results saturate with the increase of the mesh density. Hence, mesh with maximum element size of 1.5 mm is selected. The quality of the mesh elements is also an important aspect and could influence the results. Therefore, the FE mesh elements were run through an element quality check based on the aspect and Jacobian ratio. 99.8% of the elements fall under the aspect ratio of 3, with the maximum aspect ratio reaching the value of 6. Based on the Jacobian ratio results, no distorted elements were detected.
For the FEM analysis, the axial length of the primary shaft is considered equal to 325 mm, while the average inner and outer diameters corresponding to the inner shaft are 12 and 28 mm and for the outer shaft are 31 and 40 mm respectively, typical for a DCT. For the rest, the geometrical and inertial properties of a 7 speed DCT (7DCT300) are used in the simulation for the calculations. The details can be inferred from the Appendix A.

Static Analysis Results
In this section, the results obtained from the Solidworks FEM simulation for the static analysis are presented and discussed. In Table 2 the deflection at the end of the IPS is summarized for all the 3 scenarios presented before. In 2nd gear boost and recuperation mode the displacement reaches a maximum value of 46 µm (Figure 8) compared to the 4th gear, where it reaches a value of 79 µm (Figure 9). The deflection in 2nd gear is less compared to other scenarios due to presence of the BB1 next to G4, where the forces are applied. The maximum displacement is seen when the torque is transmitted from both the power sources that can reach up to 120 µm as shown in Figure 10. mesh with maximum element size of 1.5 mm is selected. The quality of the mesh elements is also an important aspect and could influence the results. Therefore, the FE mesh elements were run through an element quality check based on the aspect and Jacobian ratio. 99.8% of the elements fall under the aspect ratio of 3, with the maximum aspect ratio reaching the value of 6. Based on the Jacobian ratio results, no distorted elements were detected. For the FEM analysis, the axial length of the primary shaft is considered equal to 325 mm, while the average inner and outer diameters corresponding to the inner shaft are 12 and 28 mm and for the outer shaft are 31 and 40 mm respectively, typical for a DCT. For the rest, the geometrical and inertial properties of a 7 speed DCT (7DCT300) are used in the simulation for the calculations. The details can be inferred from the Appendix A.

Static Analysis Results
In this section, the results obtained from the Solidworks FEM simulation for the static analysis are presented and discussed. In Table 2 the deflection at the end of the IPS is summarized for all the 3 scenarios presented before. In 2nd gear boost and recuperation mode the displacement reaches a maximum value of 46 μm (Figure 8) compared to the 4th gear, where it reaches a value of 79 μm (Figure 9). The deflection in 2nd gear is less compared to other scenarios due to presence of the BB1 next to G4, where the forces are applied. The maximum displacement is seen when the torque is transmitted from both the power sources that can reach up to 120 μm as shown in Figure 10.

Rotor Dynamic Analysis
In this section are presented the modeling and results of the rotor dynamic analysis carried out using DYNROT FEM code [37] developed in MATLAB environment at our University. The main objective is to compute the deformations of the HDT primary shaft at the DMF side caused by the irregularities of the ICE and EM.
The model is shown in Figure 11. The inner and outer primary shafts have been modeled using Timoshenko beam elements (63 in total). The material properties for the beam elements are defined as AISI 4340. The ball bearings BB1 and BB2 are modeled as spring elements with the radial stiffness reported in Table 1. Each of the two needle bearings NB1 and NB2 is modeled as a combination of three spring elements working in parallel with stiffness equal to one third of the radial stiffness of the corresponding needle bearings.

Rotor Dynamic Analysis
In this section are presented the modeling and results of the rotor dynamic analysis carried out using DYNROT FEM code [37] developed in MATLAB environment at our University. The main objective is to compute the deformations of the HDT primary shaft at the DMF side caused by the irregularities of the ICE and EM.
The model is shown in Figure 11. The inner and outer primary shafts have been modeled using Timoshenko beam elements (63 in total). The material properties for the beam elements are defined as AISI 4340. The ball bearings BB1 and BB2 are modeled as spring elements with the radial stiffness reported in Table 1. Each of the two needle bearings NB1 and NB2 is modeled as a combination of three spring elements working in parallel with stiffness equal to one third of the radial stiffness of the corresponding needle bearings.

Rotor Dynamic Analysis
In this section are presented the modeling and results of the rotor dynamic analysis carried out using DYNROT FEM code [37] developed in MATLAB environment at our University. The main objective is to compute the deformations of the HDT primary shaft at the DMF side caused by the irregularities of the ICE and EM.
The model is shown in Figure 11. The inner and outer primary shafts have been modeled using Timoshenko beam elements (63 in total). The material properties for the beam elements are defined as AISI 4340. The ball bearings BB1 and BB2 are modeled as spring elements with the radial stiffness reported in Table 1. Each of the two needle bearings NB1 and NB2 is modeled as a combination of three spring elements working in parallel with stiffness equal to one third of the radial stiffness of the corresponding needle bearings.

Rotor Dynamic Analysis
In this section are presented the modeling and results of the rotor dynamic analysis carried out using DYNROT FEM code [37] developed in MATLAB environment at our University. The main objective is to compute the deformations of the HDT primary shaft at the DMF side caused by the irregularities of the ICE and EM.
The model is shown in Figure 11. The inner and outer primary shafts have been modeled using Timoshenko beam elements (63 in total). The material properties for the beam elements are defined as AISI 4340. The ball bearings BB1 and BB2 are modeled as spring elements with the radial stiffness reported in Table 1. Each of the two needle bearings NB1 and NB2 is modeled as a combination of three spring elements working in parallel with stiffness equal to one third of the radial stiffness of the corresponding needle bearings. The rotating parts are assumed to be axially symmetrical so that linearized equation of motion for the undamped system in complex coordinate frame can be written as follows [38]:   The rotating parts are assumed to be axially symmetrical so that linearized equation of motion for the undamped system in complex coordinate frame can be written as follows [38]: [M] ..
where, M is mass matrix, G is gyroscopic matrix, K is stiffness matrix, F is vector of unbalances, q = u ψ is vector of complex coordinates and ω is rotation speed.
The model includes 55 nodes in total with 2 complex degrees of freedom each, resulting in 110 degrees of freedom. The nodal displacements and rotations in complex coordinates involved in the flexural behavior are u = u x + iu y , ψ = ψ y − iψ x where u x , u y , ψ x and ψ y are the nodal displacements and rotations. The axial displacement u z and the torsional rotation ψ z are 2 additional real degrees of freedom that take the axial and torsional behavior into account, respectively. Each node has thus a total of 2 complex plus two real degrees of freedom. For this analysis only the flexural and torsional behavior are studied.
The model is used to determine the rotor Campbell diagram and critical speeds with the corresponding mode shapes. The effect of EM is evaluated considering two kind of dynamic effects:

1.
Effect of EM inertia torque during shafts angular accelerations.

2.
Effect of EM torque irregularities.
Both effects produce a torque that translates in contact forces acting on gears G1-G2-G3. Therefore, the problem has been studied by computing at varying speeds the transfer functions between the radial displacement at node 1 and the forces acting on node 54 as shown in Figure 11.
For transformation of the matrices from complex coordinates to real coordinates, the formulation in [33] has been used. Hence, the q = u ψ in complex coordinates is in real coordinates.
The equation of motion to study the flexural behavior in X-Y coordinates is written as in Equation (19). The description of variables is given in Table 3.

M xy
.. q xy + L xy + ωG xy . q xy + K ωoxy + ω 2 K ω2xy + ωH xy q xy (19) = f sxy + ω 2 f unb sin(ωt) cos(ωt) + S ixy f xy Each node in real coordinates has four degrees of freedom compared to two degrees of freedom in complex coordinates. So, there are 220 degrees of freedom in total to describe the flexural behavior of the rotating system. Dynamic Equation (19) can be reported, with explicit spin speed, in the state space representation in the following way, for X-Y behavior. The matrices are described in Table 4. . (20) y(t) = Cx(t) + Du(t) (21)  Output selection matrix Table 4. Description of the matrices.

Name Description Description
Finally, the transfer function between the displacement at node 1 and force at node 54 is computed from the state space that is used to calculate the displacement in the scenarios discussed later. The transfer function for different values of ω is depicted below in Figures 12 and 13. It can be seen that there are two peaks corresponding to the natural frequencies at 586 and 763 Hz. In Equation (19), it can be seen that the gyroscopic matrix (G xy ) is speed dependent and as a result the gyroscopic effect increases with the speed. Hence, in Figure 13, the absolute value of the transfer function increases with increase in speed.
The block shown in Figure 14 is used to calculate the displacement at the node 1, using the transfer function values shown in Figure 12. The inputs are ω, f xy (forces generated due to ICE or EM torque irregularities which are described in the next subsection) and the frequency of the applied forces (2nd order harmonics for the ICE irregularities and 3rd order harmonics for the EM torque irregularities) and the output is displacement at node 1. ( ) is speed dependent and as a result the gyroscopic effect increases with the speed. Hence, in Figure 13, the absolute value of the transfer function increases with increase in speed. The block shown in Figure 14 is used to calculate the displacement at the node 1, using the transfer function values shown in Figure 12. The inputs are , (forces generated due to ICE or EM torque irregularities which are described in the next subsection) and the frequency of the applied forces (2nd order harmonics for the ICE irregularities and 3rd order harmonics for the EM torque irregularities) and the output is displacement at node 1.   The block shown in Figure 14 is used to calculate the displacement at the node 1, using the transfer function values shown in Figure 12. The inputs are , (forces generated due to ICE or EM torque irregularities which are described in the next subsection) and the frequency of the applied forces (2nd order harmonics for the ICE irregularities and 3rd order harmonics for the EM torque irregularities) and the output is displacement at node 1.  .

Scenario 1: Effect of EM Inertia Torque
Torque irregularities coming from the ICE (2nd order) plays a crucial role in the vibrations and deflections in the transmission. In this scenario, as depicted in Figure 15, the effect of the angular acceleration acting at the input of the transmission on the generation of the EM inertia torque are analyzed.
The speed of the ICE can be written as: where, ω 0 is the average ICE speed and δω · cos(2ω 0 t) represents the rotational speed irregularity that is present due to the 2nd order harmonics of the ICE.
where, is the average ICE speed and ⋅ (2 ) represents the rotation irregularity that is present due to the 2nd order harmonics of the ICE. The angular acceleration can be obtained by differentiating the Equation (22 The EM inertia torque is then calculated as follows: = where, is the inertia of the EM, is the angular acceleration due to the irregularities and is the gear ratio between primary shaft and EM. Inertia torque ( _ ) acts on gear G3 and generates tangential ( _ dial ( _ ) and axial ( _ ) forces with a frequency that corresponds to order harmonics of the ICE, which produces the displacement at the front end of Figure 15. Setup for the EM inertia torque calculations.

Scenario 2: Effect of EM Torque Irregularities
Similar to ICE torque irregularities, the EM also have torque irregularities. T important is represented by the 3rd order frequencies of the EM speed (3 ). EM irregularities ( _ ) generate tangential ( _ ) , ( _ ) and axial ( _ ) forces that acts on the gear G3 as show ure 16 and produce the displacements at the input of the inner primary shaft. In nario is considered as 1% of the total torque being transmitted by the EM, function of the speed of the EM. .
The EM inertia torque is then calculated as follows: where, J EM is the inertia of the EM, . ω ICE is the angular acceleration due to the torque irregularities and τ G13 is the gear ratio between primary shaft and EM.
Inertia torque (T EM_inertia ) acts on gear G3 and generates tangential (F t_inertia ), radial (F r_inertia ) and axial (F a_inertia ) forces with a frequency that corresponds to the 2nd order harmonics of the ICE, which produces the displacement at the front end of IPS.

Scenario 2: Effect of EM Torque Irregularities
Similar to ICE torque irregularities, the EM also have torque irregularities. The most important is represented by the 3rd order frequencies of the EM speed (3ω EM ). EM torque irregularities (T EM_irregularities ) generate tangential F t_irregularities , radial (F r_irregularities ) and axial (F a_irregularities ) forces that acts on the gear G3 as shown in Figure 16 and produce the displacements at the input of the inner primary shaft. In this scenario T EM is considered as 1% of the total torque being transmitted by the EM, that is a function of the speed of the EM.

Rotor Dynamic Analysis Results
In this section, the results obtained from the rotor dynamic analysis carried out with DYNROT FEM code are presented and discussed. Campbell diagram for the model in Figure 11 is shown in the Figure 17. From the Campbell diagram it can be seen that there are two natural frequencies of the system at 586 and 763 Hz calculated at 0 rpm that can reach up to 590 and 764 Hz at 6000 rpm due to variation of the stiffness matrix ( ) in Equation (19) with speed. As far as the critical speeds are concerned, it can be seen that in the working range of the ICE from 800 to 6000 rpm, no critical speeds are present for synchronous speed. Moreover the 2nd and 4th order ICE harmonics do not show any critical speeds either. The only critical speed exists for the 3rd order harmonics of the EM, which is at 4700 rpm. Figure 17. Campbell diagram.
The mode shape for the first critical speed is reported below in Figure 18. It can be

Rotor Dynamic Analysis Results
In this section, the results obtained from the rotor dynamic analysis carried out with DYNROT FEM code are presented and discussed. Campbell diagram for the model in Figure 11 is shown in the Figure 17. From the Campbell diagram it can be seen that there are two natural frequencies of the system at 586 and 763 Hz calculated at 0 rpm that can reach up to 590 and 764 Hz at 6000 rpm due to variation of the stiffness matrix (K ω2xy ) in Equation (19) with speed. As far as the critical speeds are concerned, it can be seen that in the working range of the ICE from 800 to 6000 rpm, no critical speeds are present for synchronous speed. Moreover the 2nd and 4th order ICE harmonics do not show any critical speeds either. The only critical speed exists for the 3rd order harmonics of the EM, which is at 4700 rpm.

Rotor Dynamic Analysis Results
In this section, the results obtained from the rotor dynamic analysis carried out with DYNROT FEM code are presented and discussed. Campbell diagram for the model in Figure 11 is shown in the Figure 17. From the Campbell diagram it can be seen that there are two natural frequencies of the system at 586 and 763 Hz calculated at 0 rpm that can reach up to 590 and 764 Hz at 6000 rpm due to variation of the stiffness matrix ( ) in Equation (19) with speed. As far as the critical speeds are concerned, it can be seen that in the working range of the ICE from 800 to 6000 rpm, no critical speeds are present for synchronous speed. Moreover the 2nd and 4th order ICE harmonics do not show any critical speeds either. The only critical speed exists for the 3rd order harmonics of the EM, which is at 4700 rpm. The mode shape for the first critical speed is reported below in Figure 18. It can be seen that the IPS and OPS both have a deflection at the ends. The mode shape for the first critical speed is reported below in Figure 18. It can be seen that the IPS and OPS both have a deflection at the ends. Figure 19 shows the displacement due to the EM inertia torque. It can be seen that the displacement increases as the ω increase. This can be attributed to the fact that the torque irregularities coming from the ICE increases with increase in the speed of ICE up to a certain rpm (2500 rpm). Above this speed the angular acceleration is considered constant and the increase in displacement at higher rpm is due to the high value of the transfer function that are dependent on speed.  Figure 19 shows the displacement due to the EM inertia torque. It can be seen that the displacement increases as the increase. This can be attributed to the fact that the torque irregularities coming from the ICE increases with increase in the speed of ICE up to a certain rpm (2500 rpm). Above this speed the angular acceleration is considered constant and the increase in displacement at higher rpm is due to the high value of the transfer function that are dependent on speed.  Figure 20 shows the displacement due to the EM torque irregularities. It can be seen that the displacement reaches a maximum value of 92 μm at = 11750 . This is due to the fact that at this speed the 3rd order harmonics of the EM intersects with the first natural frequency as seen in the Campbell diagram ( Figure 17).  Figure 19 shows the displacement due to the EM inertia torque. It can be seen that the displacement increases as the increase. This can be attributed to the fact that the torque irregularities coming from the ICE increases with increase in the speed of ICE up to a certain rpm (2500 rpm). Above this speed the angular acceleration is considered constant and the increase in displacement at higher rpm is due to the high value of the transfer function that are dependent on speed.  Figure 20 shows the displacement due to the EM torque irregularities. It can be seen that the displacement reaches a maximum value of 92 μm at = 11750 . This is due to the fact that at this speed the 3rd order harmonics of the EM intersects with the first natural frequency as seen in the Campbell diagram ( Figure 17).  Figure 20 shows the displacement due to the EM torque irregularities. It can be seen that the displacement reaches a maximum value of 92 µm at ω ICE = 11750 rpm. This is due to the fact that at this speed the 3rd order harmonics of the EM intersects with the first natural frequency as seen in the Campbell diagram ( Figure 17).

Conclusions
Multiple sources of excitation and frequently changing operating modes may introduce potential issues in terms of reliability for P2.5 hybrid powertrain with dual clutch transmission. The paper presents static and rotor dynamic analysis of the HDT in terms of radial displacements that arises at the DMF level due to additional torque transmission from the EM. The responsible of such radial displacement is the gear train transmission that links the EM to the primary shaft. The torque is transmitted by a tangential/radial/axial force causing bending on the primary shafts. Static analysis show that the maximum displacement of the shaft end can reach 120 μm in the case of torque transmission from both the power sources. On the other hand, rotor dynamic analysis shows that the displacements due to EM inertia torque are not high and can reach up to 33 μm at 6000 rpm, while they can reach up to 92 μm due to the EM torque irregularities. The combined static and dynamic effects can lead up to a maximum displacement of 242 μm, which is below the critical values (1.6 mm) indicated in service manual of the DMF manufactures [27]. However, this radial displacement can be amplified due to the manufacturing and assembling tolerances of the HDT and may reach close to the critical values of the radial gap in the DMF. Therefore, the manufacturing and assembling tolerances should be kept below the critical values of the DMF clearance, including the deflection coming from the EM torque transmission.

Conclusions
Multiple sources of excitation and frequently changing operating modes may introduce potential issues in terms of reliability for P2.5 hybrid powertrain with dual clutch transmission. The paper presents static and rotor dynamic analysis of the HDT in terms of radial displacements that arises at the DMF level due to additional torque transmission from the EM. The responsible of such radial displacement is the gear train transmission that links the EM to the primary shaft. The torque is transmitted by a tangential/radial/axial force causing bending on the primary shafts. Static analysis show that the maximum displacement of the shaft end can reach 120 µm in the case of torque transmission from both the power sources. On the other hand, rotor dynamic analysis shows that the displacements due to EM inertia torque are not high and can reach up to 33 µm at 6000 rpm, while they can reach up to 92 µm due to the EM torque irregularities. The combined static and dynamic effects can lead up to a maximum displacement of 242 µm, which is below the critical values (1.6 mm) indicated in service manual of the DMF manufactures [27]. However, this radial displacement can be amplified due to the manufacturing and assembling tolerances of the HDT and may reach close to the critical values of the radial gap in the DMF. Therefore, the manufacturing and assembling tolerances should be kept below the critical values of the DMF clearance, including the deflection coming from the EM torque transmission.

Appendix A
For this study, the DCT data used is from the 7-speed transmission from Magna (7DCT300) [39]. Gear ratios of the DCT used in this work are given below in the Table A1: Detailed CAD data regarding the DCT can be assessed from [40].