In modern industrial production, large-scale mining gear transmission systems face significant reliability challenges due to complex operational environments. To address vibration and impact issues in these critical systems, we developed a parametric dynamic model using generalized finite element principles. This model integrates gear design parameters and system layout parameters to enable comprehensive optimization of mining gear technology.
The generalized finite element model discretizes shafts into Timoshenko beam elements with 8 degrees of freedom. The stiffness matrix for each shaft segment is expressed as:
$$K_s = f(r_i, n)$$
where $r_i$ represents shaft radius and $n$ denotes the discretization parameter. For helical gear pairs, the 8-DOF mesh model generates equations of motion:
$$m_p \ddot{x}_p + (c_m \dot{\delta} + k_m \delta) \cos \beta_b \sin \phi = 0$$
$$m_p \ddot{y}_p \pm (c_m \dot{\delta} + k_m \delta) \cos \beta_b \cos \phi = 0$$
$$\vdots$$
$$I_{zg} \ddot{\theta}_{zg} \pm (c_m \dot{\delta} + k_m \delta) r_g \cos \beta_b = 0$$
The time-varying mesh stiffness $k_m$ is calculated using potential energy method:
$$k_h = \int_0^l \frac{\pi E}{4(1 – v^2)} dy$$
$$k_b = \int_0^l \int_0^{d(y)} \frac{3[(d(y) – x)\cos\alpha_1(y) – h(y)\sin\alpha_1(y)]^2}{2Eh^3x} dx dy$$
$$k = \sum_{i=1}^n \frac{1}{\frac{1}{k_{h,i}} + \frac{1}{k_{b1,i}} + \cdots + \frac{1}{k_{f2,i}}}$$
The comprehensive system stiffness matrix combines all components:
$$K = f(r_i, z_i, m_i, L_i, \alpha, \beta_b, n, n_1, n_2)$$
leading to the global dynamic equation:
$$M \ddot{X} + C \dot{X} + K_f X = P_o + F_e$$
We implemented multi-objective optimization using NSGA-II algorithm with the following objectives:
| Objective | Expression | Purpose |
|---|---|---|
| Contact Ratio | $\min(-ε) = -\left[1.88 – 3.2\left(\frac{1}{z_1} + \frac{1}{z_2}\right)\cos\beta + \frac{B\sin\beta}{\pi m_n}\right]$ | Enhance transmission stability |
| Dynamic Performance | $\min\left(\sqrt{\frac{\sum_{i=1}^n F_i^2}{n_m}}\right)$ | Reduce vibration acceleration |
| System Volume | $\min\left(\sum L_i\pi r^2 + \sum \frac{\pi (z_1 + z_2)^2 m_n^2 B}{4\cos^2\beta_b}\right)$ | Minimize material usage |
Design variables encompassed 11 key parameters:
$$X = [z_j \quad m \quad B \quad r_1 \quad r_2 \quad r_3 \quad n \quad n_1 \quad n_2 \quad \alpha \quad \beta_i]^T$$
Optimization constraints included gear interference prevention, shaft strength requirements, transmission ratio limits, and contact fatigue constraints:
$$g_{\sigma_H,i} = \max(\sigma_{H,i} – \sigma_{HP,i}, 0)$$
$$\sigma_{H,i} = Z_H Z_E Z_{\epsilon\beta} \sqrt{\frac{F_{t,i}}{B_i d_{p,i}} \frac{u_i + 1}{u_i} K_A K_V K_{H\beta} K_{H\alpha}}$$
Case study results for a mining gearbox (500 RPM, 2000 N·m load) demonstrated significant improvements:
| Performance Metric | Optimization Strategy B2 | Improvement |
|---|---|---|
| 1st Stage Mesh Force | 2,557 N | 7.6% reduction |
| 2nd Stage Mesh Force | 19,537 N | 15.7% reduction |
| System Volume | 0.53 m³ | 27.0% reduction |
| Eccentric Load Performance | 2,589.02 N | 12.0% reduction |
Frequency domain analysis revealed reduced energy at mesh frequency harmonics, validating the effectiveness of our gear technology optimization approach. The NSGA-II algorithm successfully resolved conflicts between dynamic performance and volume objectives, generating Pareto-optimal solutions for mining gear technology applications.
This parametric modeling and multi-objective optimization framework provides a theoretical foundation for designing high-performance mining gear transmission systems. Future work will incorporate tribological considerations and material fatigue models to further advance gear technology in heavy-duty applications.