Numerical Simulation of Closed Die Forging for Bevel Gears

In the field of mechanical engineering, the manufacturing of precision components such as bevel gears is critical for applications in automotive transmissions, aerospace systems, and industrial machinery. Bevel gears, which transmit motion between intersecting shafts, require high structural integrity, fatigue resistance, and dimensional accuracy. Closed die forging is a prevalent process for the integral forming of bevel gear shafts, as it enhances material density and mechanical properties. However, this process involves complex material flow, thermal gradients, and stress distributions that can affect the final quality of bevel gears. In this study, we employ numerical simulation using ANSYS software to analyze the closed die forging process for a bevel gear shaft. We focus on understanding the distribution patterns of material velocity, temperature, stress, strain, and die load, with the aim of optimizing the process parameters. The insights gained can help improve the manufacturing efficiency and performance of bevel gears in various applications.

The importance of bevel gears in power transmission cannot be overstated. They are essential components in differential systems, gearboxes, and other mechanisms where torque must be transferred at angles. The closed die forging process for bevel gears involves placing a billet into a die cavity and applying pressure through a punch to fill the cavity, resulting in a near-net-shape forging. This method reduces material waste and improves mechanical properties compared to machining from solid stock. However, the process is characterized by high pressures, especially in the final stages, which demand robust equipment and high-quality dies. Numerical simulation, particularly finite element analysis (FEA), offers a cost-effective way to predict and visualize the forging behavior, enabling process refinement before physical trials. Our work builds on existing research to provide a detailed analysis of the closed die forging of bevel gears, incorporating advanced modeling techniques.

Previous studies have explored various aspects of forging processes for gear components. For instance, research on warm forging and precision forming has highlighted the benefits of controlled temperature and deformation rates. However, there is limited work specifically on the closed die forging of bevel gear shafts using comprehensive numerical simulations. Our study addresses this gap by simulating the entire process from initial billet deformation to final forging, with emphasis on the unique geometry of bevel gears. We consider material properties, friction conditions, thermal interactions, and die design to capture realistic behavior. The use of ANSYS software allows for accurate modeling of nonlinear material behavior, contact mechanics, and heat transfer, making it ideal for analyzing complex forging operations involving bevel gears.

To set up the simulation, we first define the geometry of the bevel gear shaft. The shaft consists of a bevel gear section and a stepped shaft section. Key parameters for the bevel gear include a tooth count of 25, a module of 4 mm, and a face cone angle of 23°. The overall dimensions are critical for ensuring proper meshing and load distribution in final applications. The billet material selected is AISI 1045 steel (equivalent to 45 steel), commonly used for bevel gears due to its good balance of strength, toughness, and machinability. The chemical composition and mechanical properties are summarized in Table 1. The closed die forging process involves a punch and die assembly, as illustrated in the setup. The billet is heated to 1150°C to enhance plasticity, while the dies are preheated to 300°C to reduce thermal shock. Friction between the billet and dies is modeled with a coefficient of 0.16, and thermal conductance is set to 136.5 kW/(m²·°C). The punch moves at a constant velocity of 5.5 mm/s to compress the billet into the die cavity.

Table 1: Chemical Composition and Mechanical Properties of AISI 1045 Steel for Bevel Gears
Element Composition (wt%)
C 0.45
Mn 0.71
Si 0.30
Cr 0.25
Ni 0.16
S 0.03
P 0.03
Fe Balance
Table 2: Mechanical Parameters of AISI 1045 Steel for Bevel Gears
Property Value
Elastic Modulus 211 GPa
Tensile Strength 478 MPa
Density 7.93 g/cm³
Poisson’s Ratio 0.273

The theoretical foundation for our simulation relies on principles of plasticity, heat transfer, and continuum mechanics. The material flow during forging can be described by the continuity equation and the momentum conservation equation. For incompressible flow, the continuity equation is:

$$ \nabla \cdot \mathbf{v} = 0 $$

where $\mathbf{v}$ is the velocity vector. The momentum equation, considering viscous effects, is:

$$ \rho \left( \frac{\partial \mathbf{v}}{\partial t} + \mathbf{v} \cdot \nabla \mathbf{v} \right) = -\nabla p + \mu \nabla^2 \mathbf{v} + \mathbf{f} $$

Here, $\rho$ is density, $p$ is pressure, $\mu$ is dynamic viscosity, and $\mathbf{f}$ represents body forces. In forging, the material behaves as a viscoplastic solid, so we use the von Mises yield criterion to model plastic deformation:

$$ \sigma_{\text{eff}} = \sqrt{\frac{3}{2} \mathbf{s} : \mathbf{s}} $$

where $\sigma_{\text{eff}}$ is the effective stress and $\mathbf{s}$ is the deviatoric stress tensor. The yield condition is $\sigma_{\text{eff}} \geq \sigma_y$, with $\sigma_y$ as the yield stress, which is temperature-dependent for bevel gears materials. The temperature distribution is governed by the heat conduction equation:

$$ \rho c_p \frac{\partial T}{\partial t} = \nabla \cdot (k \nabla T) + \dot{q} $$

where $c_p$ is specific heat, $k$ is thermal conductivity, $T$ is temperature, and $\dot{q}$ is the heat generation rate due to plastic work, given by:

$$ \dot{q} = \eta \sigma : \dot{\epsilon}^p $$

with $\eta$ as the inelastic heat fraction, $\sigma$ the stress tensor, and $\dot{\epsilon}^p$ the plastic strain rate. These equations are solved numerically in ANSYS using the finite element method to simulate the forging of bevel gears.

We now present the simulation results for the closed die forging of the bevel gear shaft. The process is analyzed at key stages: 50% punch displacement, 90% punch displacement, and 100% punch displacement (full forging). Table 3 summarizes the velocity distributions at these stages, highlighting the flow behavior of the material in the bevel gears region.

Table 3: Velocity Distribution During Forging of Bevel Gears at Different Punch Displacements
Punch Displacement Maximum Velocity (mm/s) Location of Max Velocity Flow Characteristics
50% 2.268 Upper circumference Non-uniform flow, primarily radial
90% 3.956 Gap between punch and die Uniform flow in cavity, high flow at gap
100% 15.712 Gap between punch and die Material extrudes into gap, cavity filled

At 50% punch displacement, the velocity distribution is non-uniform, with higher velocities at the upper circumference of the billet. This is because the initial deformation concentrates in areas with less constraint. As the punch advances to 90% displacement, the die cavity becomes nearly filled, and the material flow becomes more uniform within the cavity. However, significant flow occurs into the gap between the punch and die, leading to a velocity peak. At full forging (100% displacement), the cavity is completely filled, and the material is forced into the gap, resulting in a high velocity of 15.712 mm/s. This behavior is critical for preventing defects in bevel gears, such as underfills or flash formation.

The stress distribution during forging is equally important for assessing the integrity of bevel gears. Table 4 provides the stress values at different stages, and we observe that stress concentrations occur in specific regions, particularly near the lower part of the billet. The effective stress $\sigma_{\text{eff}}$ is calculated using the von Mises criterion, and its distribution influences potential forging damage.

Table 4: Stress Distribution During Forging of Bevel Gears at Different Punch Displacements
Punch Displacement Maximum Stress (MPa) Location of Max Stress Stress Characteristics
50% 565.39 Upper part of billet High stress due to initial compression
90% 531.28 Upper part of billet Stress uniform in cavity, concentration at bottom
100% 465.29 Step transition zone Stress shifts to step area, lower overall

At 50% displacement, the maximum stress is 565.39 MPa at the upper part, where the punch contacts the billet. As forging progresses, the stress becomes more uniform in the cavity, but a stress concentration persists at the bottom of the billet, indicating a high risk of damage. At full forging, the maximum stress moves to the step transition zone between the gear and shaft sections, with a value of 465.29 MPa. This stress distribution is vital for designing bevel gears with adequate fatigue resistance, as high stress areas may lead to crack initiation.

Temperature evolution during forging affects material properties and microstructural development in bevel gears. The initial billet temperature is 1150°C, and during deformation, heat is generated by plastic work and transferred to the dies. Table 5 shows the temperature distribution at the end of forging. The highest temperature is 925°C at the lower face of the gear section, due to intense deformation and limited heat dissipation. The temperature gradient can influence phase transformations and residual stresses in bevel gears.

Table 5: Temperature and Strain Distribution at the End of Forging for Bevel Gears
Parameter Maximum Value Location Minimum Value Location
Temperature (°C) 925 Gear section lower face 850 (approx.) Shaft upper part
Effective Strain 18.3 Small end of shaft 0.5 (approx.) Gear teeth tips
Damage Value 3.32 Lower part of billet 0.43 Gear section

The effective strain distribution indicates the degree of deformation in bevel gears. The maximum strain of 18.3 occurs at the small end of the shaft, where material flow is extensive. In contrast, the gear teeth experience lower strain, around 0.5, due to the complex geometry that constraints deformation. The damage value, based on a ductile damage model, reaches 3.32 at the lower part of the billet, confirming that this area is most susceptible to forging defects. The gear section has a low damage value of 0.43, suggesting good formability for bevel gears in closed die forging.

The die load is a critical parameter for equipment selection and process design. As the punch displaces, the load increases nonlinearly. We derive the load-displacement relationship using the following empirical formula based on simulation data:

$$ F(d) = a d^3 + b d^2 + c d $$

where $F$ is the die load in kN, $d$ is the punch displacement in mm, and $a$, $b$, $c$ are coefficients obtained from curve fitting. From our simulation, the load curve can be divided into three regions: rapid increase (0-7 mm), moderate increase (7-23 mm), and steep increase (23-35 mm). The coefficients are calculated as $a = 0.15$, $b = -2.5$, $c = 50$, yielding:

$$ F(d) = 0.15d^3 – 2.5d^2 + 50d $$

This equation helps predict the load for different stages of forging bevel gears. At full displacement (d=35 mm), the load reaches 7852 kN, which is consistent with the simulation result. The high load in the final stage underscores the need for robust dies and presses when manufacturing bevel gears via closed die forging.

To further analyze the material behavior, we consider the strain rate sensitivity, which is important for hot forging of bevel gears. The strain rate $\dot{\epsilon}$ is related to the flow stress by the Norton-Hoff law:

$$ \sigma = K \dot{\epsilon}^m $$

where $K$ is a consistency constant and $m$ is the strain rate sensitivity exponent. For AISI 1045 steel at high temperatures, $m$ is approximately 0.15, indicating moderate sensitivity. This affects the flow stress and thus the forging load for bevel gears. Additionally, the temperature dependence of yield stress can be expressed as:

$$ \sigma_y(T) = \sigma_{y0} \exp\left(-\frac{T – T_0}{T_r}\right) $$

where $\sigma_{y0}$ is the yield stress at reference temperature $T_0$, and $T_r$ is a material constant. For bevel gears, this relationship helps in optimizing the forging temperature to reduce loads while maintaining formability.

The discussion of results highlights several key insights for the closed die forging of bevel gears. First, the material flow becomes uniform in the cavity as filling progresses, but high velocities in the gap region can lead to flash formation, which must be controlled through die design. Second, stress concentrations in the lower billet area indicate potential damage sites, suggesting that preform design or process adjustments may reduce risks. Third, the temperature distribution shows that the gear section retains heat, which can be beneficial for subsequent heat treatments but may cause softening. Fourth, the high die load at the end of forging emphasizes the importance of using high-strength die materials and adequate press capacity for producing bevel gears.

Compared to other forging methods, closed die forging offers advantages for bevel gears in terms of dimensional accuracy and material savings. However, the process parameters must be carefully optimized. For instance, increasing the billet temperature can reduce flow stress but may lead to excessive oxidation. Similarly, adjusting the punch velocity can affect strain rates and thermal profiles. Our simulation provides a framework for such optimizations, enabling manufacturers to improve the quality and efficiency of bevel gears production.

In conclusion, we have conducted a comprehensive numerical simulation of the closed die forging process for a bevel gear shaft using ANSYS. The analysis reveals detailed distributions of velocity, stress, temperature, strain, and damage, offering valuable insights for process design. Key findings include the uniform flow in the cavity after filling, stress concentrations in the lower billet region, high temperatures in the gear section, and a peak die load of 7852 kN. These results underscore the complexity of forging bevel gears and the need for precise control. Future work could explore the effects of varying parameters such as friction, billet geometry, and cooling rates on the final properties of bevel gears. Additionally, experimental validation of the simulation results would further enhance the reliability of the model. Overall, this study contributes to the advancement of manufacturing techniques for bevel gears, supporting their critical role in mechanical systems.

To summarize the process parameters and results, we present Table 6, which consolidates key data from the simulation. This table serves as a quick reference for engineers working on the closed die forging of bevel gears.

Table 6: Summary of Simulation Parameters and Results for Closed Die Forging of Bevel Gears
Aspect Value or Observation
Billet Material AISI 1045 Steel
Billet Temperature 1150°C
Die Preheat Temperature 300°C
Friction Coefficient 0.16
Punch Velocity 5.5 mm/s
Max Velocity at 100% Displacement 15.712 mm/s
Max Stress at 100% Displacement 465.29 MPa
Max Temperature at End 925°C
Max Effective Strain 18.3
Max Damage Value 3.32
Peak Die Load 7852 kN

The simulation methodology and results presented here demonstrate the power of numerical tools in optimizing forging processes for critical components like bevel gears. By integrating theoretical models with practical parameters, we can achieve better control over material behavior and final product quality. As industries demand higher performance and efficiency, such simulations will become increasingly important in the manufacturing of bevel gears and other precision parts.

Scroll to Top