In the rapidly evolving automotive industry, electric cars have emerged as a pivotal force driving global economic growth and sustainable development. The propulsion systems of electric cars, particularly the gear transmission systems, are subjected to a myriad of nonlinear excitations due to high-speed motor operations, which can significantly impact their dynamic performance and stability. Unlike traditional internal combustion engines, electric cars often employ motors that can reach speeds exceeding 10,000 r/min, imposing unique challenges on the gear transmission systems. This study delves into the nonlinear dynamics of gear transmission systems in electric cars, focusing on the influences of motor speed excitation frequencies and internal system parameters. By constructing a comprehensive nonlinear dynamic model that accounts for time-varying meshing stiffness, time-varying support stiffness, tooth-side clearance, and bearing clearance, we aim to unravel the intricate behaviors of these systems. Through numerical simulations using the Runge-Kutta method and analytical approaches like the multi-scale method, we explore the dynamic responses, bifurcations, and resonance conditions. The integration of time-delay control parameters further enhances our understanding of stability regimes. This work not only provides insights into optimizing the design and operation of electric car gear transmissions but also underscores the importance of nonlinear dynamics in advancing the reliability and efficiency of electric vehicles.
The gear transmission system in an electric car serves as a critical component that transmits power from the high-speed motor to the wheels. Its dynamic behavior is inherently nonlinear due to factors such as varying stiffness, clearances, and external excitations. In electric cars, the absence of a traditional engine means that the motor directly couples with the gear system, leading to distinct vibrational characteristics. Understanding these nonlinear dynamics is essential for mitigating noise, vibration, and harshness (NVH) issues, improving power density, and ensuring long-term durability. This article presents a detailed analysis from a first-person perspective, where we, as researchers, investigate the system’s response under different operational conditions. We emphasize the role of key parameters like motor speed, meshing damping, and stiffness variations, all within the context of electric car applications. The use of tables and formulas facilitates a concise summary of complex relationships, while the inclusion of nonlinear dynamics principles highlights the sophistication required in modern electric car engineering.
To begin, we establish the mathematical foundation for the gear transmission system in an electric car. The system comprises multiple gears, bearings, and shafts, modeled with 30 degrees of freedom to capture translational and rotational motions. The time-varying meshing stiffness, a primary source of nonlinearity, is derived using the slice integration method for helical gears. For a gear pair, the meshing stiffness \( K_m(t) \) varies with the number of engaged teeth and can be expressed as:
$$ K_m(t) = \left[ \frac{1}{K_h(t)} + \sum_{c=1,2} \left( \frac{1}{K_{a,c}(t)} + \frac{1}{K_{b,c}(t)} + \frac{1}{K_{f,c}(t)} + \frac{1}{K_{s,c}(t)} \right) \right]^{-1}, $$
where \( K_h \) is the Hertzian contact stiffness, \( K_{a,c} \), \( K_{b,c} \), \( K_{f,c} \), and \( K_{s,c} \) represent axial compression, bending, foundation deformation, and shear stiffnesses, respectively. For helical gears in electric cars, the overlapping contact leads to a fluctuating stiffness that excites the system at the meshing frequency. The time-varying support stiffness from bearings is also considered, modeled using Hertzian contact theory for rolling elements. The bearing force components are given by:
$$ F_{bx} = \sum_{i=1}^{Z_b} \cos \gamma’_i K_b \delta^n H(\delta) \cos \theta_i, $$
$$ F_{by} = \sum_{i=1}^{Z_b} \cos \gamma’_i K_b \delta^n H(\delta) \sin \theta_i, $$
$$ F_{bz} = \sum_{i=1}^{Z_b} \sin \gamma’_i K_b \delta^n H(\delta), $$
with \( n = 2/3 \), \( Z_b \) as the number of rolling elements, \( \gamma’ \) as the contact angle, \( \delta \) as the deformation, and \( H(\delta) \) as the Heaviside function. These equations encapsulate the nonlinear bearing dynamics prevalent in electric car transmissions.
The gear meshing force incorporates dynamic backlash, which is crucial for accurate modeling. The composite meshing error \( e(t) \) is represented as a harmonic function: \( e(t) = e_a + e_b \sin(\omega_m t + \phi_0) \), where \( \omega_m \) is the meshing frequency. The nonlinear function for backlash \( f(x_n) \) defines the contact condition:
$$ f(x_n) = \begin{cases}
x_n – b_m & \text{if } x_n \geq b_m, \\
0 & \text{if } |x_n| < b_m, \\
x_n + b_m & \text{if } x_n < -b_m,
\end{cases} $$
with \( b_m \) as half the backlash. The meshing force \( F_n \) and its components are then:
$$ F_n = K(t) f(x_n) + C_m \frac{dx_n}{dt}, $$
$$ F_{mx} = F_n \cos \alpha_n \cos \beta_i, \quad F_{my} = F_n \sin \alpha_n, \quad F_{mz} = F_n \cos \alpha_n \sin \beta_i, $$
where \( C_m \) is the meshing damping, \( \alpha_n \) is the pressure angle, and \( \beta_i \) is the helix angle. These forces drive the vibrational responses in electric car gear systems.
The dynamic model of the electric car gear transmission system is derived from Newton’s second law, resulting in a set of dimensionless differential equations. For the first gear pair, the equation along the meshing line is:
$$ \frac{d^2 \bar{x}_{n1}}{d\tau^2} + 2\zeta_1 \frac{d\bar{x}_{n1}}{d\tau} \cos \alpha_n + \kappa_1(\tau) f(\bar{x}_{n1}) \cos \alpha_n = \bar{f}_0 + \bar{f} \cos(\omega \tau), $$
where \( \bar{x}_{n1} \) is the dimensionless displacement, \( \zeta_1 \) is the equivalent meshing damping ratio, \( \kappa_1(\tau) \) is the dimensionless time-varying meshing stiffness, \( \bar{f}_0 \) is the static load, \( \bar{f} \) is the dynamic load amplitude, and \( \omega \) is the excitation frequency ratio. Similar equations apply to the second gear pair and bearing supports. The dimensionless parameters are summarized in Table 1, which provides key values used in our simulations for electric car applications.
| Parameter | Symbol | Value | Description |
|---|---|---|---|
| Equivalent mass | \( m_e \) | 1.0 kg | Base mass for normalization |
| Natural frequency | \( \omega_n \) | \( \sqrt{K_p / m_e} \) rad/s | Reference frequency |
| Meshing damping ratio | \( \zeta_1, \zeta_2 \) | 0.05, 0.03 | For first and second gear pairs |
| Stiffness fluctuation amplitude | \( \kappa \) | 0.0437 (first), 0.1623 (second) | Dimensionless time-varying stiffness |
| Static load | \( \bar{f}_0 \) | 0.06 (first), 0.04 (second) | Dimensionless external torque |
| Backlash half-width | \( b_m \) | 1.0 (normalized) | Reference for displacement |
| Pressure angle | \( \alpha_n \) | 20° | Gear geometry parameter |
| Helix angle | \( \beta_i \) | 22° (first), 12.5° (second) | For helical gears in electric cars |
To solve these equations, we employ the fourth-order Runge-Kutta method, a numerical technique suitable for handling nonlinearities. The dynamic responses are analyzed through bifurcation diagrams, phase portraits, Poincaré sections, and frequency spectra. These tools help identify periodic, chaotic, and bifurcation behaviors in electric car gear systems. For instance, the excitation frequency ratio \( \omega_e = \omega_m / \omega_n \) is varied to study its impact. Results show that the system exhibits rich nonlinear dynamics, with stability regions interspersed with chaotic zones. This is critical for electric cars, where motor speed variations can induce undesirable vibrations.
The influence of motor speed excitation frequency on the electric car gear system is profound. As \( \omega_e \) increases, the system transitions through different dynamic states. For the first gear pair, when \( 0.20 < \omega_e < 0.23 \) and \( 0.26 < \omega_e < 0.32 \), the motion is periodic, indicating stable operation. However, for \( 0.23 < \omega_e < 0.25 \) and \( 0.32 < \omega_e < 0.48 \), chaotic behavior emerges, leading to instability. Similarly, the second gear pair shows periodic motion for \( 0.40 < \omega_e < 0.52 \), \( 0.70 < \omega_e < 1.02 \), and \( \omega_e > 1.08 \), but chaos appears in intervals like \( 0.53 < \omega_e < 0.70 \). These findings imply that electric car motors should avoid operating at specific speed ranges to maintain system stability. Table 2 summarizes the stability intervals based on \( \omega_e \), which can guide the design of control strategies for electric vehicles.
| Gear Pair | Stable Intervals (\( \omega_e \)) | Unstable Intervals (\( \omega_e \)) | Dynamic Behavior |
|---|---|---|---|
| First Pair | 0.20–0.23, 0.26–0.32, \( \omega_e > 0.48 \) | 0.23–0.25, 0.32–0.48 | Periodic to chaotic transitions |
| Second Pair | 0.40–0.52, 0.70–1.02, \( \omega_e > 1.08 \) | 0.53–0.70, 1.02–1.08 | Chaos and bifurcations observed |
Meshing damping is another vital parameter for electric car gear systems. We analyze the effect of the damping ratio \( \zeta_1 \) on the first gear pair. As \( \zeta_1 \) increases from 0.03 to 0.13, the system evolves from chaotic to periodic motion, enhancing stability. This suggests that higher damping can suppress nonlinear oscillations in electric car transmissions. For the second gear pair, variations in \( \zeta_2 \) lead to complex behaviors: chaos for \( 0.020 < \zeta_2 < 0.024 \), bifurcation at \( \zeta_2 = 0.024 \), and a return to chaos before stabilizing at higher values. These trends underscore the need for optimal damping design in electric cars to mitigate vibrations. The relationship between damping and stability is quantified through Lyapunov exponents, but for brevity, we focus on the qualitative outcomes relevant to electric vehicle applications.
Time-varying meshing stiffness also plays a crucial role. In electric cars, the stiffness fluctuations arise from gear tooth interactions and manufacturing errors. We model \( \kappa(\tau) \) as a periodic function with amplitude \( \kappa \). The system’s response to stiffness variations is studied through parametric analysis. For instance, increasing \( \kappa \) amplifies the nonlinear effects, potentially leading to resonance conditions. This is particularly important for electric cars, where lightweight materials may alter stiffness properties. Table 3 lists the stiffness parameters and their effects, derived from our simulations.
| Stiffness Parameter | Value Range | Impact on System | Recommendation for Electric Cars |
|---|---|---|---|
| Fluctuation amplitude \( \kappa \) | 0.01–0.20 | Higher values increase vibration amplitudes | Keep below 0.15 for stability |
| Mean stiffness \( K_p \) | 2.5 × 10⁷ N/m | Affects natural frequency and resonance | Optimize for operating speed ranges |
| Stiffness phase | 0–2π rad | Influences synchronization with excitations | Align to minimize dynamic loads |
To delve deeper into resonance phenomena, we apply the multi-scale method to the system with time-delay control parameters. This approach is essential for electric cars, where active control systems can introduce delays. The governing equation with time-delay feedback is:
$$ \frac{d^2 \bar{x}_{n1}}{d\tau^2} + \left[ 2\zeta_1 \frac{d\bar{x}_{n1}}{d\tau} + \left( 1 + \varepsilon \kappa \cos(\omega \tau) \right) f(\bar{x}_{n1}) \right] \cos \alpha_n = \bar{f}_0 + \bar{f} \cos(\omega \tau) + g_d \bar{x}_{n1}(\tau – \tau_d) + g_v \frac{d\bar{x}_{n1}}{d\tau}(\tau – \tau_v), $$
where \( g_d \) and \( g_v \) are displacement and velocity control gains, and \( \tau_d \), \( \tau_v \) are time delays. Using the multi-scale expansion with \( T_0 = \tau \), \( T_1 = \varepsilon \tau \), we derive the amplitude-frequency response equation for primary resonance. Assuming \( f(\bar{x}_{n1}) \approx \rho_1 \bar{x}_{n1} + \rho_2 \bar{x}_{n1}^3 \) with \( \rho_1 = \omega_n^2 \), the steady-state amplitude \( a \) satisfies:
$$ \left( a \zeta_1 \cos \alpha_n – W_1 \right)^2 + \left( a \sigma – W_2 \right)^2 = W_3^2, $$
with \( \sigma = \omega – \omega_0 \) as the frequency detuning, and:
$$ W_1 = \frac{g_v a \omega_0 \cos(\omega_0 \tau_v) – g_d a \sin(\omega_0 \tau_d)}{2\omega_0}, $$
$$ W_2 = \frac{3\rho_2 a^3 \cos \alpha_n + 12\rho_2 a \bar{f}_0^2 \cos \alpha_n}{8\omega_0} + \frac{g_d a \cos(\omega_0 \tau_d) – g_v a \omega_0 \sin(\omega_0 \tau_v)}{2\omega_0}, $$
$$ W_3 = \frac{\bar{f} – \omega_0^2 \kappa \bar{f}_0 \cos \alpha_n}{2\omega_0}. $$
The stability condition is derived from the Jacobian of the averaged equations:
$$ \left[ \frac{g_v \omega_0 \cos(\omega_0 \tau_v) – g_d \sin(\omega_0 \tau_d)}{2\omega_0} – \zeta_1 \cos \alpha_n \right]^2 + \left[ \sigma – \frac{3\rho_2 a^2 \cos \alpha_n + 12\rho_2 \bar{f}_0^2 \cos \alpha_n}{8\omega_0} \right] \left[ \sigma – \frac{12\rho_2 a^2 \cos \alpha_n + 12\rho_2 \bar{f}_0^2 \cos \alpha_n}{8\omega_0} \right] > 0. $$
This inequality defines the stable operating regions for electric car gear systems under time-delay control. We analyze the effects of control parameters on primary resonance. For example, increasing the displacement gain \( g_d \) initially amplifies the response but can reduce instability at higher values. Similarly, the velocity gain \( g_v \) influences the amplitude growth, with optimal ranges enhancing stability. The external load fluctuation \( \bar{f} \) directly affects resonance peaks; thus, in electric cars, controlling torque ripples from the motor is vital. These insights are summarized in Table 4, which links control parameters to system performance in electric vehicle contexts.
| Parameter | Typical Range | Effect on Amplitude \( a \) | Stability Implication for Electric Cars |
|---|---|---|---|
| Displacement gain \( g_d \) | 0–0.5 | Increases then decreases with optimal near 0.3 | Moderate gains improve vibration suppression |
| Velocity gain \( g_v \) | 0–0.5 | Monotonically increases amplitude | Keep low to avoid excessive oscillations |
| Time delay \( \tau_d, \tau_v \) | \( T/9 \) (used) | Phase shifts can stabilize or destabilize | Design delays to match system dynamics |
| Load fluctuation \( \bar{f} \) | 0–0.1 | Directly proportional to resonance peak | Minimize through motor control in electric cars |
Experimental validation is conducted using a parallel-shaft gearbox test rig, simulating electric car transmission conditions. The setup includes a drive motor, torque sensor, gearbox, and brake. Measurements at 1000 r/min and 1 N·m torque show good agreement with simulation results in both time and frequency domains. The dominant frequencies align, confirming the model’s accuracy for electric car applications. This validation underscores the practicality of our nonlinear dynamics approach in real-world electric vehicle systems.
In conclusion, the gear transmission system in electric cars exhibits complex nonlinear dynamics due to factors like time-varying stiffness, backlash, and bearing clearances. Our analysis reveals that motor speed excitation frequencies can induce chaotic behaviors in specific intervals, necessitating careful speed management in electric vehicles. Meshing damping and stiffness parameters significantly influence stability, with optimal designs enhancing performance. The incorporation of time-delay control parameters offers avenues for active vibration suppression, crucial for the refinement of electric cars. By avoiding unstable speed ranges and tailoring internal parameters, the reliability and efficiency of electric car gear transmissions can be substantially improved. Future work may explore advanced materials and real-time control algorithms to further optimize these systems for the burgeoning electric car industry.
Throughout this study, the term “electric car” has been emphasized to highlight the application context. The nonlinear dynamics principles discussed are universally applicable but gain particular relevance in the high-speed, high-performance environment of electric vehicles. As the adoption of electric cars accelerates, understanding and mitigating vibrational issues will remain a key research frontier, driving innovations in transmission design and control strategies.