Design and Performance Simulation Analysis of a Novel Inter-Wheel Differential Using Spur and Pinion Gear Planetary Mechanism

In the pursuit of enhancing automotive efficiency and handling performance, reducing the size and weight of drivetrain components is paramount. Traditional differentials, while effective, often impose constraints on packaging and mass optimization. This paper introduces a novel inter-wheel differential design based on a spur and pinion gear planetary arrangement. The primary innovation lies in its compact axial configuration, achieved through strategic gear positioning and modification, which significantly reduces volume and weight compared to conventional designs. The spur and pinion gear system is central to this advancement, offering simplicity and reliability. I will detail the structural principles, derive the necessary kinematic and dynamic equations, present the design parameters, and validate the design through comprehensive virtual prototyping and simulation using ADAMS. The results confirm that this spur and pinion gear-based differential operates correctly, provides equitable torque distribution, and is ideally suited for applications where installation space is limited or lightweighting is critical.

The fundamental role of a differential is to allow wheels on the same axle to rotate at different speeds during cornering, while transmitting torque from the drivetrain. Conventional inter-wheel differentials often use bevel gears or planetary gearsets with helical gears, which can be bulky. The proposed design employs a spur and pinion gear planetary mechanism. Spur gears are characterized by straight teeth parallel to the gear axis, offering ease of manufacturing and assembly. Pinion gears here refer to the smaller planetary gears. The key structural insight is to arrange the sun gears and planet gears such that the meshing zones of the left and right sides are axially overlapped, thereby minimizing the overall axial footprint. This is made possible through careful application of gear profile shifting (modification).

The core assembly consists of two sun gears (Sun Gear I and Sun Gear II), two planet gears (Planet Gear I and Planet Gear II), and a differential carrier (case). The sun gears are connected to the left and right half-shafts, respectively. The planet gears are mounted on a common planet pinion shaft within the carrier and mesh with each other. Crucially, Sun Gear I meshes with Planet Gear I, and Sun Gear II meshes with Planet Gear II. The meshing between the two planet gears is the central spur and pinion gear interaction that facilitates the differential action. To prevent physical interference between the components of the left and right branches, the gear pairs are designed with specific profile shift coefficients. The condition for interference-free operation and balanced torque transmission in straight-line motion is that the pitch circle radii satisfy the following relation derived from geometric and force balance considerations:

$$ \frac{r_a}{r_b} = \frac{r_2}{r_1} $$

Here, \( r_a \) and \( r_b \) are the pitch radii of Sun Gear I and Sun Gear II, respectively, while \( r_1 \) and \( r_2 \) are the pitch radii of Planet Gear I and Planet Gear II. This relationship ensures that the tangential forces at the meshes are balanced when the differential is not differentiating.

Kinematic analysis describes the rotational speed relationships among the components. Let \( \omega_0 \) be the angular velocity of the differential carrier. Let \( \omega_a \) and \( \omega_b \) be the angular velocities of Sun Gear I and Sun Gear II, respectively. Let \( \omega_1 \) and \( \omega_2 \) be the angular velocities of Planet Gear I and Planet Gear II relative to the carrier (their absolute angular velocities are \( \omega_0 + \omega_1 \) and \( \omega_0 + \omega_2 \), but for analysis, we consider the relative spin). The fundamental kinematic equations are derived from the rolling contacts at the gear meshes. For the mesh between Sun Gear I and Planet Gear I:

$$ (\omega_a – \omega_0) r_a = -\omega_1 r_1 $$

For the mesh between Sun Gear II and Planet Gear II:

$$ (\omega_b – \omega_0) r_b = \omega_2 r_2 $$

For the mesh between Planet Gear I and Planet Gear II:

$$ \omega_1 r_1 = -\omega_2 r_2 $$

Combining these equations, we can eliminate the planet gear speeds and derive the classic differential kinematic equation for this specific geometry. Adding the first two equations and substituting the third yields:

$$ \omega_a r_a + \omega_b r_b = 2 \omega_0 \frac{r_a r_b}{r_a + r_b} \left( \frac{r_1 + r_2}{r_1 r_2} \right) $$

However, under the design condition \( r_a / r_b = r_2 / r_1 \), this simplifies significantly. A more practical form for simulation and analysis is obtained by solving for the half-shaft speeds. The general motion characteristic equation is:

$$ \omega_a = \omega_0 – \frac{r_1}{r_a} \omega_1 $$
$$ \omega_b = \omega_0 + \frac{r_2}{r_b} \omega_2 $$
$$ \text{and } \omega_1 r_1 = -\omega_2 r_2 $$

Therefore, the differential action is governed by the relative spin of the planet gears. When the vehicle travels straight, there is no relative spin (\( \omega_1 = \omega_2 = 0 \)), and both sun gears rotate at the carrier speed: \( \omega_a = \omega_b = \omega_0 \). During a turn, the planet gears spin, causing one sun gear to accelerate and the other to decelerate relative to the carrier, satisfying \( \omega_a + \omega_b = 2\omega_0 \) only when \( r_a = r_b \), which is not strictly required here due to the design ratio, but the principle holds that the average speed is proportional to the carrier speed.

Dynamic analysis aims to determine the torque distribution between the half-shafts and the internal forces. The input torque \( T_0 \) is applied to the differential carrier. This torque is transmitted to the planet gears via the carrier pins. Let \( F_{o1} \) and \( F_{o2} \) be the forces exerted by the carrier on Planet Gear I and Planet Gear II, respectively. Assuming six planet gear pairs (a common design for load sharing), the force from the carrier on one planet pinion is:

$$ F_{o1} = \frac{T_0}{6 (r_a + r_1)} $$
$$ F_{o2} = \frac{T_0}{6 (r_b + r_2)} $$

Given the design condition \( r_a / r_b = r_2 / r_1 \), it follows that \( F_{o1} r_1 = F_{o2} r_2 \), ensuring moment balance on the planet pinion shaft. When the differential is operating (e.g., during a turn), friction at the planet gear pinion shaft bearings generates an internal locking torque. The friction torque \( T_f \) arises from the radial forces on the pinion shafts. A detailed force balance on the planet gears yields expressions for the output torques \( T_a \) and \( T_b \) on the sun gears. The torque difference between the half-shafts is:

$$ \Delta T = T_b – T_a = \frac{\mu r_d T_0}{6} \left( \frac{1}{r_a + r_1} + \frac{1}{r_b + r_2} \right) $$

where \( \mu \) is the sliding friction coefficient between the planet pinion shaft and its bushing, and \( r_d \) is the radius of the pinion shaft’s thrust collar. The locking coefficient \( K \), defined as the ratio of the torque difference to the input torque, is:

$$ K = \frac{\Delta T}{T_0} = \frac{\mu r_d}{6} \left( \frac{1}{r_a + r_1} + \frac{1}{r_b + r_2} \right) $$

The tooth mesh forces are critical for stress analysis. The force between Sun Gear I and Planet Gear I, \( F_{a1} \), and between Sun Gear II and Planet Gear II, \( F_{b2} \), can be expressed as:

$$ F_{a1} = \frac{T_a / 6 + T_0 / 12}{r_a \cos \alpha_1} $$
$$ F_{b2} = \frac{T_b / 6 – T_0 / 12}{r_b \cos \alpha_2} $$

where \( \alpha_1 \) and \( \alpha_2 \) are the pressure angles of the respective sun-planet meshes. The force between the two planet gears, \( F_{12} \), is:

$$ F_{12} = \frac{T_0}{6 r_0 \cos \alpha_0} $$

Here, \( r_0 \) is the pitch radius of the planet-planet mesh and \( \alpha_0 \) is its pressure angle. These formulas are essential for evaluating gear strength and designing the spur and pinion gear teeth.

The design process followed automotive engineering standards. The primary parameters for the spur and pinion gears were selected to meet torque capacity, durability, and packaging constraints. The gear geometry was defined using profile shifting to achieve the necessary center distances and clearance. The key design parameters are summarized in the table below.

Gear Sun Gear I Planet Gear I (Pinion) Sun Gear II Planet Gear II (Pinion)
Number of Teeth (z) 36 13 36 13
Module (m) [mm] 3
Pressure Angle (α) [°] 20
Center Distance [mm] 71 76
Face Width [mm] 22 22 22 40
Profile Shift Coefficient (x) -1.0099 0.3 1.2296 -0.3
Pitch Radius (r) [mm] 52.16 18.84 55.84 20.16
Addendum Circle Diameter [mm] 107.20 46.059 120.80 42.662
Dedendum Circle Diameter [mm] 94.441 33.3 107.878 29.7
Base Circle Diameter [mm] 101.487 36.648 101.487 36.648

These parameters ensure the kinematic condition \( r_a / r_b = r_2 / r_1 \) is satisfied (52.16/55.84 ≈ 18.84/20.16 ≈ 0.934). The spur and pinion gear teeth are standard involute profiles with modification to avoid undercutting and to provide adequate backlash and contact ratio.

A virtual prototype was constructed using Pro/ENGINEER (Creo Parametric) based on the parameters above. The three-dimensional solid model included all differential components: the carrier, sun gears, planet gears, pinion shaft, and housing features. The model was then imported into ADAMS/View using the Parasolid format for multibody dynamics simulation. In ADAMS, materials were assigned (steel for gears, density 7.85e-6 kg/mm³), and constraints were applied: revolute joints between the carrier and ground (for input rotation), between the carrier and each planet gear (allowing rotation of the planet about its pinion shaft), and between each sun gear and ground (allowing rotation for output). The gear meshes were modeled using contact forces rather than kinematic joints to capture the dynamics of tooth engagement, including potential separation and impact.

The contact model in ADAMS used a Hertzian impact-based formulation. For the spur and pinion gear contacts, the normal force was computed as:

$$ F_n = K \delta^e + C \dot{\delta} $$

where \( \delta \) is the penetration depth, \( e \) is the force exponent (set to 1.5), \( C \) is the damping coefficient, and \( K \) is the stiffness coefficient. The stiffness \( K \) for two contacting cylinders is given by:

$$ K = \frac{4}{3} R^{1/2} E^* $$

with \( \frac{1}{R} = \frac{1}{R_1} + \frac{1}{R_2} \) (effective radius) and \( \frac{1}{E^*} = \frac{1-\nu_1^2}{E_1} + \frac{1-\nu_2^2}{E_2} \) (effective modulus). For the gear material (40CrNi2Mo), \( E = 2.06 \times 10^5 \) N/mm², \( \nu = 0.28 \). The calculated stiffness values for the three meshes were approximately \( K_1 = 3.9 \times 10^5 \) N/mm³/², \( K_2 = 3.31 \times 10^5 \) N/mm³/², and \( K_3 = 3.63 \times 10^5 \) N/mm³/². Damping was set to 50 N·s/mm, and the penetration depth was 0.1 mm. Friction was modeled with a static coefficient of 0.05 and a dynamic coefficient of 0.03 at the sliding velocity corresponding to the operating speed.

The simulation was conducted for two primary operating conditions: straight-line driving and cornering. The input carrier speed was set to 1.389 rad/s (approximately 13.27 rpm) to represent a low-speed, high-torque scenario. The simulation time was 0.3 seconds with a step size of 0.001 seconds to capture transient dynamics.

Straight-Line Driving Simulation: For straight-line motion, equal resisting torques were applied to both sun gears, each set to -1,050,000 N·mm. The kinematic results showed that the angular velocities of both sun gears converged rapidly to the carrier speed, confirming the expected motion: \( \omega_a = \omega_b = \omega_0 \). The velocity profiles exhibited minor oscillations due to tooth contact impacts but stabilized around the mean. The torque transmission was analyzed next. The input torque on the carrier fluctuated around a mean of 2,092,899 N·mm, closely matching the theoretical input of 2,100,000 N·mm. The output torques on the sun gears had mean values of 1,043,981.5 N·mm and 1,032,872 N·mm, which are nearly equal and sum to the input torque (allowing for small simulation losses). The small discrepancy (within 5%) is attributable to the dynamic contact model and friction losses. The following table summarizes the torque results.

Torque Component Maximum Value (N·mm) Minimum Value (N·mm) Mean Value (N·mm) Theoretical Value (N·mm) Error vs. Theory
Input (Carrier) 2,169,432 2,016,367 2,092,899 2,100,000 -0.34%
Output (Sun Gear a) 1,078,636 1,009,327 1,043,982 1,050,000 -0.57%
Output (Sun Gear b) 1,067,360 998,384 1,032,872 1,050,000 -1.63%

The gear mesh forces were also extracted. The force between Sun Gear I and Planet Gear I oscillated around a mean of 6,595 N. The force between Sun Gear II and Planet Gear II averaged 6,568 N. The force between the two planet gears (the central spur and pinion gear pair) averaged 7,216 N. These align well with theoretical calculations using the derived formulas, which gave approximately 6,898 N, 6,896 N, and 6,899 N, respectively. The slight variations are due to dynamic effects and the discrete nature of the contact simulation. The consistency validates the force distribution within the spur and pinion gear system.

Cornering Simulation: To simulate a turning condition, unequal resisting torques were applied to the sun gears, reflecting the different road loads on the inner and outer wheels. Based on the calculated locking coefficient (using \( \mu = 0.26 \), \( r_d = 12 \) mm), the torque difference was approximately 60,201 N·mm. Therefore, resisting torques of -1,110,201 N·mm and -989,799 N·mm were applied to Sun Gear a and Sun Gear b, respectively. The kinematic results clearly showed differential action: the sun gear with the lower resisting torque (Sun Gear b) accelerated relative to the carrier, while the other decelerated. The average speeds obtained were \( \omega_a = 0.1204 \) rad/s and \( \omega_b = 10.1524 \) rad/s (values are illustrative; actual simulation values in rad/s). The relationship \( \omega_a + \omega_b \approx 2\omega_0 \) held approximately, confirming the basic differential function.

The dynamic results for cornering are summarized in the table below. The input carrier torque mean was 2,100,459 N·mm. The output torques averaged 976,950 N·mm and 1,108,587.5 N·mm, showing the expected bias due to the internal friction. The torque difference was about 131,637.5 N·mm, which is higher than the simple theoretical calculation due to the dynamic engagement conditions and the modeled friction characteristics. Nevertheless, the trend is correct.

Torque Component Maximum Value (N·mm) Minimum Value (N·mm) Mean Value (N·mm) Theoretical Value (N·mm) Error vs. Theory
Input (Carrier) 2,199,033 2,001,885 2,100,459 2,100,000 +0.02%
Output (Sun Gear a) 1,038,356 915,544 976,950 989,799 -1.30%
Output (Sun Gear b) 1,185,159 1,032,016 1,108,588 1,110,201 -0.15%

The mesh forces during cornering were slightly asymmetric but remained within safe limits. The planet-planet mesh force averaged 6,960 N, while the sun-planet mesh forces averaged 6,914 N and 6,824 N. These are close to the theoretical predictions, confirming that the spur and pinion gear teeth are not overloaded during differential action. The simulation also allowed observation of the contact patterns and force transients, which showed no signs of abnormal shock or sustained vibration, indicating a robust design.

In conclusion, the novel differential design utilizing a spur and pinion gear planetary mechanism has been successfully developed and validated. The key achievement is a significant reduction in axial length compared to conventional designs, achieved through intelligent layout and gear profile modification. The kinematic and dynamic analyses provided a theoretical foundation, which was corroborated by detailed multibody dynamics simulations. The spur and pinion gear assembly performed reliably under both straight-line and cornering conditions, distributing torque appropriately and demonstrating a small, controllable locking effect due to friction. The design parameters, including the profile shift coefficients, were crucial in achieving an interference-free, compact package. This differential is particularly advantageous for electric vehicles, lightweight sports cars, and applications where drivetrain space is at a premium. Future work could involve physical prototyping and testing to validate durability and refine the gear micro-geometry for noise and efficiency optimization. The success of this project underscores the potential of rethinking traditional automotive components through fundamental mechanical design principles and advanced simulation tools.

The entire process from concept to simulation reinforces the versatility and effectiveness of spur and pinion gear systems in transmitting power and enabling controlled differential motion. The use of profile-shifted spur gears allowed for the precise alignment and spacing needed for this compact layout. The simulation results confirm that the meshing behavior of the spur and pinion gear pairs is stable and predictable, even under dynamic loading conditions. This work contributes a viable alternative to existing differential architectures, one that prioritizes compactness and weight savings without compromising functional performance. The spur and pinion gear approach thus stands as a testament to the enduring value of well-executed gear design in automotive engineering.

Scroll to Top