In the evolving landscape of engineering technology, hydraulic systems have become indispensable across numerous industrial applications. Within these systems, the gear pump stands as a critical component due to its simplicity, cost-effectiveness, compact design, and reliable operation. Among the various types, the external gear pump is widely utilized. Its performance and longevity are heavily influenced by the integrity of its core rotating element—the gear shaft. As the gear shaft transmits torque and withstands complex dynamic loads, its structural reliability is paramount. This study focuses on conducting a detailed finite element analysis (FEA) of the gear shaft from a CB-B125 external gear pump. By leveraging SolidWorks for precise three-dimensional modeling and ANSYS for advanced computational analysis, I aim to evaluate the stress distribution and deformation under operational conditions, thereby providing a theoretical foundation for optimized design practices. The gear shaft is the central focus, as its failure can lead to catastrophic pump malfunction.
The process begins with creating an accurate digital twin of the gear shaft. I employ SolidWorks, a powerful CAD software, to develop a parametric three-dimensional model. The initial step involves defining the gear geometry. Instead of modeling the complex involute profile directly within SolidWorks, I use specialized gear design software, CAXA, to generate the precise tooth profile curve based on key parameters: module, pressure angle, and number of teeth. This profile is exported in a DWG format, which is universally compatible. Importing this DWG file into SolidWorks, I use it as a sketch to create the gear’s base feature via the “Extrude Boss/Base” command. For the CB-B125 pump’s gear shaft, the extrusion depth is set to 35 mm to match the gear width. Subsequently, the shaft sections on both ends are modeled by sketching concentric circles on the gear’s end faces and extruding them to the required lengths. Critical features like the keyway, necessary for transmitting torque from the motor, are then cut using the “Extruded Cut” command. Finally, fillets and chamfers are added to edges to reflect the actual manufacturing process and reduce stress concentrations, though for subsequent FEA, minor fillets may be simplified. The final three-dimensional model of the driving gear shaft is a crucial digital asset for analysis.
With the gear shaft model complete, the next phase is to prepare it for structural analysis. The transition from CAD to CAE is facilitated by the seamless data exchange between SolidWorks and ANSYS. To ensure compatibility, I save the SolidWorks assembly containing the gear shaft in the Parasolid (*.x_t) format, a neutral geometry kernel format that preserves geometric data accurately. It is essential to use non-Chinese, alphanumeric filenames to avoid import errors in ANSYS. Launching ANSYS Workbench or the classic ANSYS interface, I use the import utility (File > Import > PARA…) to load the Parasolid file. This process successfully transfers the detailed geometry of the gear shaft into the ANSYS environment, where it can be discretized into finite elements. The integrity of this imported gear shaft geometry is foundational for an accurate simulation.
The core of this investigation is the static structural analysis of the gear shaft. The objective is to determine the stress and strain fields when the gear shaft is subjected to its operational loads. The gear shaft material is specified as 40Cr alloy steel, a common choice for such components due to its good strength and toughness. Its material properties, essential for the FEA, are summarized in the table below.
| Property | Symbol | Value | Unit |
|---|---|---|---|
| Elastic Modulus | E | 2.1e5 | MPa |
| Poisson’s Ratio | ν | 0.3 | – |
| Density | ρ | 7.8e3 | kg/m³ |
| Tensile Strength | σB | 700 | MPa |
| Yield Strength (Approx.) | σs | 550 | MPa |
Before meshing, a careful analysis of the forces acting on the gear shaft is conducted. The gear shaft is supported by bearings at two locations (points B and D in a simplified schematic). It receives input torque from an electric motor via the keyed connection (point A). At the gear mesh region (point C), it experiences the reaction forces from engaging with the driven gear. These include the tangential force (circumferential force) and the radial force. A free-body diagram is essential, and the governing equations for these forces are derived from power transmission fundamentals. The input torque \( M_A \) is related to the motor power \( P \) and angular speed \( \omega \):
$$ M_A = \frac{P}{\omega} $$
For the CB-B125 gear pump operating at typical conditions, if the power is \( P \) (in Watts) and speed is \( n \) (in rpm), then \( \omega = \frac{2 \pi n}{60} \) rad/s. The tangential force \( F_{YC} \) at the pitch circle of the gear is:
$$ F_{YC} = \frac{2 M_A}{d} $$
where \( d \) is the pitch diameter of the gear. The radial force \( F_{ZC} \) is calculated using the pressure angle \( \alpha \):
$$ F_{ZC} = F_{YC} \cdot \tan(\alpha) $$
These calculated forces are applied in the FEA model. The support reactions at the bearing locations (\( F_{XB}, F_{YB}, F_{XD}, F_{YD} \)) are determined by the solver based on equilibrium conditions.
Within ANSYS, the geometry of the gear shaft is prepared for meshing. Small cosmetic features like very small fillets and chamfers that have negligible impact on global stress results are suppressed to simplify the mesh and reduce computational cost without sacrificing accuracy for the gear shaft’s primary stress zones. The element type selected is SOLID185, an 8-node brick element suitable for 3-D modeling of solid structures. It has plasticity, creep, stress stiffening, and large deflection capabilities, though for this linear static analysis, its basic formulation suffices. A global element size is defined, and a mapped meshing approach is used in regular sections of the gear shaft, while a free mesh is applied to the complex gear tooth and keyway regions. The resulting finite element model is detailed below.
| Aspect | Detail |
|---|---|
| Element Type | SOLID185 (3-D 8-Node Structural Solid) |
| Meshing Method | Combination of Mapped & Free |
| Number of Nodes | 16,283 |
| Number of Elements | 9,141 |
| Geometry Features Handled | Gear Teeth, Keyway, Shaft Steps |
The application of boundary conditions and loads is a critical step that defines the simulation scenario. Constraints are applied to represent the physical supports. On the cylindrical surfaces of the gear shaft that interface with bearings or pump cover bores, displacement constraints are applied in the transverse directions (Y and Z axes), allowing only rotation around the shaft axis (X-axis). On the side faces of the gear (not the teeth), displacement in the X-direction is constrained to represent axial location. The most significant load application is the torque. In FEA, torque is often applied as a tangential pressure or force on a surface. For the keyway, which is the torque input location, a distributed surface load is applied on one of the vertical faces of the keyway. The magnitude of this pressure \( p \) is derived from the torque \( M_A \) and the effective radius \( r_{key} \) from the shaft center to the centroid of the keyway face:
$$ p = \frac{F}{A_{face}} = \frac{M_A / r_{key}}{A_{face}} $$
where \( A_{face} \) is the area of the keyway face. Alternatively, a pure moment can be applied via remote forces. Simultaneously, to balance this input torque, the calculated tangential force \( F_{YC} \) is applied in the opposite direction along the line of action on the gear teeth flank in the mesh zone. The radial force \( F_{ZC} \) is applied perpendicular to this direction. This setup creates a state of static equilibrium for the gear shaft.
With the model fully defined, the ANSYS solver is executed for a linear static analysis. The solution provides comprehensive data on nodal displacements, stresses, and strains. Post-processing reveals the critical results. The von Mises equivalent stress is a key indicator for ductile materials like 40Cr steel, as it predicts yielding under multiaxial stress states. The contour plot (not explicitly referenced by number but described) shows a clear stress distribution across the gear shaft. The maximum von Mises stress is found to be concentrated in a specific region. Detailed quantitative results are extracted and summarized.
| Result Parameter | Maximum Value | Location | Unit |
|---|---|---|---|
| Von Mises Equivalent Stress | 206.12 | Root fillet region adjacent to the keyway | MPa |
| Total Deformation | 0.015 | At the free end of the shaft extension | mm |
| Maximum Principal Stress | 198.5 | Similar to max von Mises location | MPa |
| Safety Factor (based on yield) | Calculated as per below | Minimum at max stress location | – |
The analysis conclusively identifies the region near the keyway as the most critically stressed part of the gear shaft. This is a common location for stress concentration due to the sudden change in geometry—the keyway acts as a notch. The maximum stress value of 206.12 MPa must be evaluated against the material’s strength. For 40Cr steel with a yield strength \( \sigma_s \approx 550 \) MPa and applying a safety factor \( n_s = 1.3 \) as per common design practice for dynamic loads, the allowable stress \( [\sigma] \) is:
$$ [\sigma] = \frac{\sigma_s}{n_s} = \frac{550}{1.3} \approx 423.08 \text{ MPa} $$
Since the calculated maximum stress (206.12 MPa) is significantly less than the allowable stress (423.08 MPa), the gear shaft design is deemed safe under the specified operating conditions. The strain distribution follows a predictable pattern, with minimal deformation observed, confirming the gear shaft’s adequate stiffness.
To further generalize the findings, it is insightful to consider the stress concentration factor (\( K_t \)) for the keyway. While the FEA provides the direct result, theoretical estimates can complement it. The stress at the keyway corner can be expressed as:
$$ \sigma_{max} = K_t \cdot \sigma_{nom} $$
where \( \sigma_{nom} \) is the nominal torsional shear stress in the shaft, calculated as \( \tau_{nom} = \frac{16 M_A}{\pi d^3} \) for a solid circular shaft. For a shaft with a keyway, \( K_t \) values typically range from 2 to 3 for torsion. This analytical approach reinforces why the FEA shows elevated stress specifically at that geometric discontinuity on the gear shaft.
The integration of SolidWorks and ANSYS presents a robust workflow for the virtual prototyping and analysis of critical components like the gear shaft. The three-dimensional modeling capability of SolidWorks allows for accurate representation of complex geometries, including the involute gear teeth and the keyway on the gear shaft. The associative nature of the model means design changes can be quickly updated and re-analyzed. The direct import into ANSYS eliminates manual geometry recreation, reducing errors and saving time. The finite element analysis provides a deep insight into the performance of the gear shaft that theoretical calculations alone cannot offer, especially in identifying localized stress concentrations. This process is not limited to the gear shaft alone; it can be extended to other pump components like the housing or bearings for a comprehensive system analysis.
In conclusion, this detailed finite element analysis of the gear pump gear shaft, facilitated by SolidWorks and ANSYS, successfully maps the stress and deformation behavior under working loads. The gear shaft, while generally safe within the elastic limit, exhibits a predictable vulnerability near the keyway region due to stress concentration. This finding is critical for designers, prompting considerations for design enhancements such as optimizing keyway geometry, using larger fillet radii, or even considering alternative torque transmission methods for highly loaded applications. The methodology demonstrated—from precise CAD modeling to rigorous CAE simulation—constitutes a powerful toolset for advancing the reliability and performance of gear pumps and similar mechanical systems. Future work could involve dynamic analysis, fatigue life prediction, or multiphysics simulations including thermal effects on the gear shaft.