In the field of mechanical engineering, the straight spur gear is one of the most widely used transmission components due to its high efficiency, constant transmission ratio, long service life, and compact structure. However, when the gear precision is low, the transmission noise and vibration become significant, affecting the working accuracy and stability of the gear system. Accurate modal analysis of the straight spur gear is essential to understand its vibration characteristics during the design stage, thereby eliminating excessive vibration or noise and improving the ride comfort, handling stability, and fuel economy of vehicles. Although many researchers have performed finite element modal analysis on gears to obtain natural frequencies and mode shapes, the directly established finite element model inevitably involves simplifications in geometry and loading conditions, leading to discrepancies from the real dynamic behavior. Therefore, it is necessary to update the finite element model using experimental modal data. The finite element model updating method adjusts certain initial parameters (such as physical parameters like density, elastic modulus, Poisson’s ratio, geometric dimensions, and boundary conditions) within a certain range to make the dynamic characteristics calculated by the theoretical model as consistent as possible with the actual structure.
This paper presents a finite element modal model updating approach for the straight spur gear based on physical parameter modification. Starting from the finite element modal analysis of the straight spur gear, we take the natural frequencies obtained from experimental modal analysis as the target, and use the ANSYS optimization design module as a tool to adjust the physical parameters of the gear. The work provides a reliable basis for further dynamic simulation analysis of the straight spur gear structure.
Finite Element Modal Analysis of the Straight Spur Gear
Theoretical Background
According to vibration theory and the finite element method, the equation of motion for a straight spur gear can be written as:
$$
[M]\{\ddot{X}\} + [C]\{\dot{X}\} + [K]\{X\} = \{f(t)\}
$$
where $\{X\} = \{x_1, x_2, \ldots, x_n\}^T$; $[M]$, $[C]$, $[K]$ are the global mass, damping, and stiffness matrices, respectively; $\{\ddot{X}\}$, $\{\dot{X}\}$, $\{X\}$ are the acceleration, velocity, and displacement vectors of the structure. In the absence of external forces, i.e., $f(t)=0$, the system performs free vibration. The corresponding undamped vibration equation is:
$$
[M]\{\ddot{X}\} + [K]\{X\} = \{0\}
$$
The associated eigenvalue problem is:
$$
([K] – \omega^2[M])\{\Psi\} = \{0\}
$$
Solving the homogeneous equation yields $n$ eigenvalues $\omega$ (natural frequencies) and the corresponding eigenvectors $\{\Psi\}$ (mode shapes). Each natural frequency and mode shape represents the free vibration of a single-degree-of-freedom system.
Finite Element Modeling of the Straight Spur Gear
The involute profile of the straight spur gear is generated using the parametric equation:
$$
r_k = r_b / \cos \alpha_k
$$
$$
\theta_k = \text{inv}\alpha_k = \tan\alpha_k – \alpha_k
$$
where $r_k$ is the radius vector, $r_b$ is the base circle radius, $\alpha_k$ is the pressure angle, and $\theta_k$ is the involute angle. Using the APDL (ANSYS Parametric Design Language) and a loop to generate multiple consecutive points on the involute curve, we fit a spline curve to quickly and accurately create the end face involute of the straight spur gear. The basic parameters of the gear are: number of teeth $z=16$, module $m=6$ mm, pressure angle $\alpha=20^\circ$, profile shift coefficient $x=0.38$, initial elastic modulus $E=210$ GPa, density $\rho=7800$ kg/m³, and Poisson’s ratio $\mu=0.3$.
Considering the cyclic symmetry of the gear, we model only a single tooth sector. The mesh uses Solid45 elements. The finite element mesh of the single sector is shown conceptually. The complete model of the straight spur gear is obtained using cyclic symmetry analysis. To simplify the modeling, the root transition curve is approximated by a tangent arc, and small fillets are neglected.
Modal Parameters from Finite Element Analysis
Modal analysis of the straight spur gear is performed to obtain the modal parameters under free-free conditions. No constraints or loads are applied. The Block Lanczos method is used for eigenvalue extraction. Since lower-order natural frequencies and mode shapes dominate the structural vibration, we extract the first six natural frequencies. The results are listed in the table below.
| Mode | Natural Frequency (Hz) | Mode Shape Description |
|---|---|---|
| 1 | 5909.2 | 1st radial mode |
| 2 | 5909.2 | 1st radial mode |
| 3 | 7161.6 | 1st bending mode |
| 4 | 7161.6 | 1st bending mode |
| 5 | 14574 | 2nd radial mode |
| 6 | 14574 | 2nd radial mode |
The first six mode shapes are characterized as follows: Modes 1 and 2 are primarily radial stretching in the first order with no axial vibration (1st radial mode). Modes 3 and 4 are the first-order bending vibration of the teeth (1st bending mode). Modes 5 and 6 are the second-order radial vibration (2nd radial mode). These mode shapes are typical for a straight spur gear.
Experimental Modal Analysis of the Straight Spur Gear
Modal Test Setup
For the experimental modal analysis, the straight spur gear is suspended by highly elastic rubber ropes to approximate free-free boundary conditions. This ensures that the first flexible mode frequency is more than 20 times higher than the highest rigid body mode frequency, thus the measured modes represent the gear’s elastic modes. The frequency resolution of the frequency response function (FRF) measurement is set to 2 Hz to improve accuracy. The roving hammer impact method is used for excitation.
Because the straight spur gear is symmetric in both radial and axial directions, there are many coupled modes on the two end faces and the cylindrical face. To accurately decouple them, two accelerometers are placed on each face during separate tests. The FRF matrix rows are measured to identify repeated roots. The reciprocity check between points 4 and 48 on the cylindrical face shows excellent agreement in the range 0-15 kHz, confirming that the measured FRFs are of high quality and satisfy Maxwell’s reciprocity.
Experimental Modal Parameter Identification
The least-squares complex frequency domain method combined with the complex modal indicator function and stabilization diagram is used to identify the closely spaced repeated modes. The stabilization diagrams for the cylindrical face and end face clearly indicate two coupled modes around 6000 Hz, 7000 Hz, and 14000 Hz. The identified natural frequencies and damping ratios are given in the table below.
| Mode | Experimental Frequency (Hz) | Damping (%) | Relative Error (%) (compared to initial FE) | Mode Shape Description |
|---|---|---|---|---|
| 1 | 5745.5 | 0.054 | 2.85 | 1st radial mode |
| 2 | 5766.1 | 0.090 | 2.48 | 1st radial mode |
| 3 | 7104.9 | 0.438 | 0.80 | 1st bending mode |
| 4 | 7105.6 | 0.445 | 0.79 | 1st bending mode |
| 5 | 14151.0 | 0.099 | 2.99 | 2nd radial mode |
| 6 | 14202.2 | 0.074 | 2.62 | 2nd radial mode |
The relative errors between the initial finite element model and the experimental data are small (maximum 2.99%), indicating that the initial model is basically acceptable. However, to achieve higher accuracy for dynamic predictions, model updating is necessary.
Model Updating of the Straight Spur Gear Finite Element Model
Formulation of the Updating Problem
The inaccuracies in the finite element model of the straight spur gear arise from two main sources: structural modeling errors and parameter modeling errors. Since the gear geometry is built according to design drawings, the structural model is relatively accurate. Therefore, we focus on adjusting the material parameters to improve the model. The model updating problem is essentially a constrained structural optimization problem. The objective is to minimize the sum of absolute relative errors between the calculated and experimental natural frequencies by adjusting the design variables within specified bounds.
Mathematically, the optimization can be expressed as:
$$
\min f(\mathbf{v}) = \sum_{i=1}^{6} \frac{|F_i – f_i|}{F_i}
$$
subject to:
$$
\mathbf{v}_{\min} \le \mathbf{v} \le \mathbf{v}_{\max}
$$
$$
w_{\min} \le w \le w_{\max}
$$
where the design variable vector $\mathbf{v} = (E, \rho, \mu)$ includes the elastic modulus, density, and Poisson’s ratio. The state variable vector $w = (\phi_1, \ldots, \phi_6)$ with $\phi_i = |F_i – f_i|/F_i$ is constrained to be less than 2% for each mode. $F_i$ is the experimental natural frequency, $f_i$ is the finite element computed natural frequency, and $f(\mathbf{v})$ is the objective function representing the sum of the six relative errors.
Optimization Using ANSYS
The optimization is performed within the ANSYS optimization module using the zero-order method. After 13 iterations, the optimal design variables and the corresponding objective function value are obtained. The initial and optimized values are compared in the table below.
| Parameter | Initial Value | Optimized Value |
|---|---|---|
| Elastic Modulus $E$ (Pa) | $2.10 \times 10^{11}$ | $2.07 \times 10^{11}$ |
| Density $\rho$ (kg/m³) | 7800.0 | 7941.8 |
| Poisson’s Ratio $\mu$ | 0.3000 | 0.2998 |
| Objective $f(\mathbf{v})$ (%) | 12.53 | 5.74 |
Updated Finite Element Modal Results
Using the optimized material parameters, a new finite element modal analysis of the straight spur gear is performed. The natural frequencies are recalculated and compared with the experimental values. The results are summarized in the following table.
| Mode | Experimental Frequency (Hz) | Initial FE Frequency (Hz) | Updated FE Frequency (Hz) | Initial Relative Error (%) | Updated Relative Error (%) |
|---|---|---|---|---|---|
| 1 | 5745.5 | 5909.2 | 5809 | 2.85 | 1.11 |
| 2 | 5766.1 | 5909.2 | 5809 | 2.48 | 0.74 |
| 3 | 7104.9 | 7161.6 | 7039 | 0.80 | -0.93 |
| 4 | 7105.6 | 7161.6 | 7039 | 0.79 | -0.94 |
| 5 | 14151.0 | 14574 | 14327 | 2.99 | 1.24 |
| 6 | 14202.2 | 14574 | 14327 | 2.62 | 0.88 |
As shown in Table 4, the maximum relative error among the first six modes is reduced from 2.99% to 1.24% after updating. The sum of absolute relative errors decreases from 12.53% to 5.74%, demonstrating a significant improvement in the accuracy of the finite element modal model for the straight spur gear. Moreover, the updated model now provides a consistent set of damping ratios (obtained from the experimental data) that can be incorporated into further dynamic analyses.
Discussion
The successful updating of the straight spur gear finite element model confirms that the approach of using APDL parametric modeling combined with material property optimization is both efficient and effective. The initial model, built using standard involute equations and mesh generation, inherently contains simplifications such as the straight root fillet and omission of small chamfers. Instead of refining the geometry (which could be time-consuming), we compensate for these errors by adjusting the effective material properties within a physically plausible range. The optimized density (7941.8 kg/m³) is slightly higher than the nominal steel density, which may account for the added mass from the neglected fillets and other details. The elastic modulus (207 GPa) and Poisson’s ratio (0.2998) remain close to typical values for gear steel, indicating that the adjustment is reasonable.
The reduction of the maximum relative error to 1.24% ensures that the updated finite element model can reliably predict the dynamic behavior of the straight spur gear under various operating conditions. This is crucial for subsequent analyses such as forced response, noise and vibration prediction, and fatigue life estimation. Additionally, the experimental modal analysis provides damping ratios, which are essential for transient and harmonic response simulations. The updated model thus serves as a robust baseline for future design iterations and optimization of the straight spur gear.
Conclusions
In this work, we have presented a systematic procedure for updating the finite element modal model of a straight spur gear. The key conclusions are:
- The initial finite element model of the straight spur gear, built using standard involute geometry and APDL, yields natural frequencies that deviate from experimental data by up to 2.99%. This discrepancy necessitates model updating to improve accuracy.
- By formulating the model updating as an optimization problem, we minimized the sum of relative errors of the first six natural frequencies using the ANSYS optimization module. The design variables were the elastic modulus, density, and Poisson’s ratio.
- After 13 iterations, the optimized material parameters (E=207 GPa, ρ=7941.8 kg/m³, μ=0.2998) reduced the maximum relative error to 1.24%, and the sum of errors from 12.53% to 5.74%. The updated finite element model now matches the experimental modal characteristics with high fidelity.
- The updated model also provides a complete set of modal parameters, including damping ratios from the experiment, enabling accurate dynamic simulations of the straight spur gear in subsequent engineering analyses.
This study demonstrates that the combination of parametric finite element modeling and physical-parameter-based optimization is a practical and efficient method for obtaining a reliable dynamic model of complex structures such as the straight spur gear. The methodology can be extended to other gear types and mechanical components.