Next Article in Journal
Editorial for Special Issue “Leaching of Rare Earth Elements from Various Sources”
Next Article in Special Issue
The Effect of Hydrodynamic Conditions on the Selective Flotation of Fully Liberated Low Grade Copper-Nickel Ore
Previous Article in Journal
On Hydromechanical Interaction during Propagation of Localized Damage in Rocks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Two Mass Variables on Inertia Cone Crusher Performance and Optimization of Dynamic Balance

National Engineering Research Center for Equipment and Technology of Cold Rolling Strip, Yanshan University, Qinhuangdao 066004, Hebei, China
*
Author to whom correspondence should be addressed.
Minerals 2021, 11(2), 163; https://doi.org/10.3390/min11020163
Submission received: 17 January 2021 / Revised: 30 January 2021 / Accepted: 2 February 2021 / Published: 4 February 2021
(This article belongs to the Special Issue Process Optimization in Mineral Processing)

Abstract

:
Inertia cone crushers are widely used in complex ore mineral processing. The two mass variables (fixed cone mass and moving cone mass) affect the dynamic performance of the inertia cone crusher. Particularly the operative crushing force of the moving cone and the amplitude of the fixed cone are affected, and thus the energy consumption of the crusher. In this paper, the process of crushing steel slag is taken as a specific research object, to analyze the influence of two mass variables on the inertia cone crusher performance. A real-time dynamic model based on the multi-body dynamic (MBD) and the discrete element method (DEM) is established. Furthermore, the influence of the fixed cone mass and moving cone mass on the operative crushing force, amplitude and average power draw are explored by the design of simulation experiments. The predictive regression models of inertia cone crusher performance are obtained using response surface methodology (RSM). After increasing the fixed cone mass, the optimized amplitude, average power and moving cone mass are decreased by 37.1%, 33.1% and 10%, respectively, compared to without the adjustment. Finally, a more effective dynamic balancing mechanism of inertia cone crusher is achieved, which can utilize the kinetic energy of a balancer, and minimize the mass of the fixed and moving cone. The fixed cone mass and moving cone mass of a balancing crusher are decreased by 78.9% and 22.8%, respectively, compared to without the balancing mechanism.

1. Introduction

Inertia cone crushers are widely used in the secondary and tertiary crushing stages of complex ore processing, such as the comprehensive recovery of steel slag [1,2]. A mantle rotates and swings in the crushing chamber, which is due to an eccentric vibrator transferring the rotational motion to the main shaft. As it flows downward between the mantle and concave, the ore particle is crushed several times. The total crushing force for the inertia cone crusher is provided by the eccentric vibrator and mantle. As the concave is located above several rubber absorbers, the concave can move and roll in three-dimensional space. Therefore, the operative crushing force is less than the theoretical force, and the energy consumption increases [3]. The subgroup including the concave and subsidiary components is defined as the fixed cone, and the subgroup including the mantle and eccentric vibrator is defined as the moving cone. At the condition of keeping other parameters invariable, the fixed cone mass and moving cone mass have a great impact on the operative crushing force, amplitude of fixed cone and energy consumption, whereas the increase in moving cone mass can increase the theoretical crushing force and amplitude directly. The decrease of fixed cone mass can decrease the operating crushing force and increase energy consumption indirectly. Here, from a manufacturer’s perspective, how to determine the two mass parameters is a key problem. Furthermore, at the guarantee of reasonable crushing force achievement rate and energy consumption, minimizing the mass of the fixed and moving cone is one of the main ways to reduce manufacturing cost.
Savov et al. [4] and Xia et al. [5] contributed to an initial mathematical modeling of the crushing force achievement rate. However, the models do not have the ability to take into account the effect of ore particles on the crusher. Additionally, no research regarding the mass of inertia cone crusher optimization has yet been published to our knowledge. Cleary et al. [6] and Andre et al. [7] studied the effect of feed properties (material strength, particle friction) and machine controls (CSS, speed) on cone crusher performances (particle distribution, throughput, power) based on the particle replacement model (PRM) in the software EDEM. Chen et al. [8] took the throughput and crushing force as the multi-objective optimization, and studied the effect of the parameters of the crushing chamber and speed on gyratory crusher performances based on the bonded particle model (BPM). These above studies can provide the main ways to optimize the variables (operation, chamber shape and feed properties). However, these simulation methods do not have the ability to take into account the effect of inertia parameters (fixed cone and moving cone mass) on operation performance (crushing force, amplitude and average power) for an inertia cone crusher. Cheng et al. [9] provided a powerful method whereby the coupling multi-body dynamics (MBD) [10,11] and discrete element method (DEM) [12,13] simulate the crushing behavior response for an inertia cone crusher. Currently no research using coupled MBD–DEM dynamic models for an inertia cone crusher has been published, except for our publication [9]. Barrios et al. [14] and Chung et al. [15], regarding coupled MBD–DEM models, provided useful attempts for high-pressure grinding rolls (HPGR). Furthermore, at the same industrial scale, the mass of the inertia cone crusher is much heavier than other cone crushers, such as hydraulic crushers and spring crushers. The reason is that the eccentric vibrator leads to the violent vibration of the crusher and the increase of energy consumption. Znamenskll et al. [16], regarding the dynamic balance of an inertia cone crusher, put forward a preliminary design. However, the dynamic balance design neither completely counteracts the excitation force nor utilizes the kinetic energy of the balancer. Ren et al. [17] utilized the kinetic energy of the balancer in the design. Nevertheless, the dynamic balancing mechanism is unstable and cannot completely counteract the excitation force, so it is difficult to widely use it in industry.
As such, this paper takes the process of crushing steel slag as the analysis object, and the crushing force achievement rate, amplitude of the fixed cone and average power of the drive shaft are explored by the MBD–DEM coupling method. Moreover, we propose an approach in which the use of response surface methodology (RSM) and analysis of variance (ANOVA) optimizes the fixed and moving cone mass to achieve the optimum operation performance. The results show that the operation performance is greatly improved by increasing the fixed cone mass, which increases the manufacturing cost. Finally, in order to reduce the manufacturing cost for manufacturers and the running cost for users, a more effective dynamic balancing mechanism of inertia cone crushers is achieved. Such a mechanism not only utilizes the kinetic energy of the balancer, but also minimizes the mass of the fixed cone and the moving cone.

2. The Coupled Model for Inertia Cone Crusher

2.1. Inertia Cone Crusher Theory

The inertia cone crusher consists of a main frame, a concave, a mantle, rubber absorbers, a main shaft and an eccentric vibrator. The ore particles fall from the feed chute to the crushing chamber; then, they are squeezed by the mantle and other particles. Finally, the particles are discharged from the discharge zone. Figure 1 shows a vertical cross-section and a simplified MBD model for an inertia cone crusher, where α and θ are the mantle angle and the nutation angle, respectively. l0, l1 and l2 are the axis of the crusher, concave and mantle, respectively. B1 is the fixed cone, which is fixed to the main frame. B2 is the moving cone, which is fixed to the main shaft. O1 is a spherical joint; O2 is a spherical joint between B1 and the globe bearing (B4); O3 is a cylindrical joint between B2 and the eccentric vibrator (B3); O4 is a planar joint between B3 and B4; O5 is a ball-pin joint between B3 and the connecting shaft (B5); O6 is a universal joint between B5 and the drive shaft (B6); O7 is a revolute joint between B6 and the ground (B0).

2.2. Crusher Dynamic Model Using MBD

The generalized coordinate qi of the rigid body Bi is shown in Equation (1).
q i = r i T , Λ i T T = x i ,   y i ,   z i ,   ψ i ,   θ i ,   φ i T ,   i = 1 , ..... , 6
where ri is the matrix of independent position variables (xi, yi, zi), and Λi is the matrix of independent angle variables (ψi, θi, φi). According to [18], the kinetic formula for Bi without any joint equations is derived in Equation (2).
  M i q ¨ i = Q i
where q ¨ i is the generalized acceleration of qi, and Mi and Qi can be expressed, respectively:
M i = m i E 3 × 3 0 3 × 3 0 3 × 3 J i ( i ) D i ,   Q i = F i a + F i f + F i g + F i c + F i p J i ( i ) D ˙ i + D i Λ ˙ J i ( i ) D i Λ ˙ i + T i a ( i ) + T i f ( i ) + T i c ( i ) + T i p ( i )
where mi and Ji are the mass and the inertia matrix for Bi, respectively. E and 0 denote the identity and null matrix, respectively. Fia and Tia denote the equivalent absorber forces and torques, Fif and Tif denote the joint friction forces and torques, Fig is the gravity of Bi, Fip and Tip denote the equivalent particle forces and torques, and Di is the coordinate matrix.
The formula of MBD for the inertia cone crusher without any joints can be shown as:
M q ¨ = Q M = diag ( M 1 ,   M 2 ,   M 3 ,   M 4 ,   M 5 ,   M 6 )   Q = ( Q 1 T ,   Q 2 T ,   Q 3 T ,   Q 4 T ,   Q 5 T ,   Q 6 T ) T
where q ¨ is the generalized accelerations of the multi-body system. The joint Oj (j = 1,2,…,7) equations and driving motion constraints for inertia cone crusher are expressed as:
Φ ( q , t ) = 0 ,   Φ q = Φ q = Φ i q j   ;   i = 1 , ..... , 26 ,   j = 1 , ..... , 7
According to Equation (5), the velocity and acceleration equation can be expressed as:
Φ q q ˙ = Φ t Φ q q ¨ = ζ ,   ζ = Φ q q ˙ q q ˙ + 2 Φ q t q ˙ + Φ t t
whereΦq and Φt is the Jacobian matrix for Φ(q,t). According to the Lagrange multiplier λi (i = 1, 2, …, 26), the formula of MBD for inertia cone crusher is derived in terms of the Lagrange multiplier matrix λ and generalized coordinate matrix q, and can be shown as:
M Φ q T Φ q 0 q ¨ λ = Q ζ ,     λ = λ 1 ,   λ 2 ,     ,   λ 26 T

2.3. DEM Modeling of Slag Particles Using BPM

2.3.1. BPM Theory

The BPM consists of bonding a packed distribution of particles, forming a breakage cluster [19]. As shown in Figure 2, a parallel bonding beam is created between each particle in contact, so the forces (torques) on the bonding beam are calculated from Equation (8) and Equation (9). BPM has been used in simulating the crushing behavior of particles [13,20,21].
δ F b n = k b n A Δ U n ,   δ F b t = k b t A Δ U t ,   δ T b n = k b n J Δ Θ n ,   δ T b t = k b t J 2 Δ Θ t
where Fbn, Tbn, Fbt and Tbt denote the bond normal force, normal torque, bond tangential-directed force and torque, respectively. kbn and kbt are the normal and tangential stiffness per unit area. A and J are the area of the parallel bond cross-section and polar moment of inertia, respectively.
Δ U n = v n δ t   ,   Δ U t = v t δ t ,   A = π R b 2 ,   Δ Θ n = ω n δ t ,     Δ Θ t = ω t δ t ,     J = 1 2 π R b 4
where Rb, Vn, Vt, ωn, ωt, and δt are parallel bond radius, normal velocity, tangential velocity, normal angular velocity, tangential angular velocity and time step, respectively.
The maximum normal and tangential stress are calculated according to Equation (10).
σ b max = F b n A + 2 T b t J R b max < σ b c , τ b max = F b t A + T b n J R b max < τ b c
where σbc and τbc are critical normal and critical shear strength, respectively. If the maximum stress exceeds the critical strength, the bond beam will disappear. The particle interaction depends on the Hertz–Mindlin contact model [22].

2.3.2. BPM Calibration

The feed material is a steel slag which has a complex shape and size distribution in our industrial experiments of the inertia cone crusher. The feed particle size range is 50 mm to 70 mm, and a 3D model of slag is shown in Figure 3. When using BPM for simulating the crushing behavior of steel slag, the key is to make sure that the relevant DEM parameters of the particle are calibrated. Therefore, the Hertz–Mindlin contact model parameters were determined by the uniaxial compression deformation and repose angle test, as shown in Figure 4. Table 1 shows the contact parameters.
The BPM-relevant parameters were determined by the Brazilian test and simulation experiment, as shown in Figure 5. The calibration method had been described in detail in our publication [13]. By comparing the tensile strength simulation with the experiment values (10.6 MPa), we directly provide the BPM parameters that are shown in Table 1.

2.4. The Solution of the Coupled MBD-DEM Model in Software

Combining with Section 2.1, Section 2.2 and Section 2.3, Figure 6 shows the simulation flowchart of the inertia cone crusher using the MBD–DEM coupling method. The MBD of the geometries is calculated by RecurDyn software, and the DEM of the particles is calculated using EDEM. As slag clusters flow downward between the mantle and the concave, the size becomes successively smaller in the software EDEM, as shown in Figure 7. The operative performances (the operative crushing force, amplitude and average power) for the inertia cone crusher are obtained by the software RecurDyn.

3. Analysis of the Inertia Cone Crusher Performance

3.1. Influencing Factors and Performance Goals

Basing on Section 2 and the previous publications [3,4,5], the fixed cone mass (FM) m1 and the moving cone mass (MM) m2 are taken as the influencing factors. The operative crushing force Fo is less than the theoretical crushing force Ft under normal operating conditions. The theoretical crushing force of the inertia cone crusher is provided by the moving cone when the fixed cone is not moving. As such, the crushing force achievement rate ηf is taken as one of the performance goals, according to Equation (11).
η f = 100 % × F o / F t
when the inertia cone crusher works, the fixed cone will move horizontally and deflect around a rotation point. As such, the amplitude of the rotation point displacement and the deflection angle can give a good indication for the vibration characteristics of the crusher. The theoretical and experimental results indicate that the amplitude As and deflection angle γd have a significant positive correlation [22,23]. Ignoring the deflection angle γd, the amplitude As is taken as a performance goal. Besides this, the two mass variables have a great impact on the energy consumption of a crusher, and the average power draw of drive shaft Pa is taken as a performance goal. Based on the response surface methodology (RSM) [24], the influence of the two mass variables on the performance goals is modeled with the SPSS software. The corresponding predictive regression models can be expressed as Equation (12), and the simulated experiment scheme is listed in Appendix A, Table A1.
η f = f η ( m 1 , m 2 ) A s = f A ( m 1 , m 2 ) P a = f P ( m 1 , m 2 ) ; 10 m 1 45 3 m 2 6.5
where fη, fA, and fP are the predictive models of crushing force achievement rate, amplitude and average power, respectively. Because the driving speed is determined by the productivity, the driving speed of the crusher (model: GYP1200) should not exceed 550 rpm. As such, the driving speed is set to 550 rpm in this paper.

3.2. Crushing Force Achievement Rate Analysis

The influence of FM and MM on the crushing force achievement rate, and the prediction regression curves are shown in Figure 8. The theoretical force Ft is only determined by MM, and Figure 8a shows the influence of MM on the theoretical force Ft. Figure 8b shows the influence of FM on the crushing force achievement rate under three kinds of MM (3.5 t, 4.5 t, and 5.9 t), which indicates that the crushing force achievement rate ηf increases significantly with increasing FM. Figure 8c shows the influence of MM on ηf under three kinds of FM (15 t, 25 t, and 39 t), which indicates that ηf gradually decreases with increasing MM. At the 0.05 significance level, we can find that the influence of FM and MM on ηf is significant using the quadratic regression function.
Figure 8d shows the influence of the interaction between FM and MM on ηf using the response surface methodology (RSM), and the corresponding results of ANOVA for predictive regression models are shown in Table 2. The prediction regression model of the crushing force achievement rate is expressed as:
η f = f η ( m 1 , m 2 ) = 82 . 196 + 0 . 628 m 1 - 0 . 011 m 1 2 - 0 . 557 m 2 2 + 0 . 083 m 1 m 2 ;   10 m 1 45 ,   3 m 2 6.5
The corresponding ANOVA for predictive regression coefficients is shown in Table A2. At the 0.05 significance level, we can find that the linear term of MM has the weakest impact on ηf, and the p-value is more than 0.05. As such, the linear term of MM is ignored in the prediction regression model of crushing force achievement rate. As FM increases or MM decreases, ηf can increase. This is because the increase in FM or the decreases in MM will decrease the amplitude of the fixed cone, and increase the eccentric distance of the moving cone. However, when ηf is over 90%, the increase in ηf is very small with increasing FM, and for the different moving cone mass, the fixed cone mass required to reach 90% (ηf) is different.

3.3. Amplitude Analysis

The influence of FM and MM on the amplitude of the fixed cone, and the prediction regression curves, are shown in Figure 9. Figure 9a shows the influence of FM on the amplitude under three kinds of MM (3.5 t, 4.5 t, and 5.9 t), which indicates that the amplitude As decreases with increasing FM. Figure 9b shows the influence of MM on the amplitude under three kinds of FW (15 t, 25 t, and 39 t), which indicates that the amplitude As increases with increasing MM. Figure 9c shows the influence of the interaction between FM and MM on the crushing force achievement rate ηf using RSM, and Table 3 indicates that the prediction model of As has a good fitness using the quadratic regression function with the 0.01 significance level. The prediction regression model is shown in Equation (14).
A s = f A ( m 1 , m 2 ) = 3 . 708 + 2 . 162 m 2 - 0 . 197 m 1 + 0 . 003 m 1 2 - 0 . 039 m 1 m 2 ;   10 m 1 45 ,   3 m 2 6.5
It can be found that the quadratic term of MM has the weakest impact on amplitude, and the p-value is more than 0.01 from Table A2. So, the quadratic term of MM is ignored in the predictive regression model of amplitude. As FM increases, the amplitude As can decrease, and the fixed cone will be more difficult to move for a constant theoretical crushing force. Comparing Figure 8d and Figure 9c, it can be found that when the ηf is over 90%, the decrease in amplitude As is very small with the increasing FM. However, with the change of MM, the amplitude As will change significantly for a constant FM.

3.4. Average Power Analysis

The influence of FM and MM on the average power draw, and the prediction regression curves, are shown in Figure 10. Figure 10a shows the influence of FM on the average power under three kinds of MM (3.5 t, 4.5 t, and 5.9 t), which indicates that the average power Pa decreases with increasing FM. Figure 10b shows the influence of MM on Pa under three kinds of FM (15 t, 25 t, and 39 t), which indicates that the average power Pa increases significantly with increasing MM. Figure 10c shows the influence of the interaction between FM and MM on Pa using RSM, and Table 3 indicates that the prediction model of average power Pa has a good fitness using the quadratic regression function with the 0.05 significance level. The prediction regression model is shown in Equation (15). Table A2 shows that the quadratic term of MM has the weakest impact on the average power Pa, and the p-value is more than 0.05.
A s = f A ( m 1 , m 2 ) = - 56 . 084 + 64 . 107 m 2 - 1 . 679 m 1 + 0 . 068 m 1 2 - 0 . 979 m 1 m 2 ;   10 m 1 45 ,   3 m 2 6.5
Comparing Figure 9 and Figure 10, we can see that as FM and MM increase, the variation in average power Pa is similar to that of the amplitude As. The reason for this is that as As increases, the energy consumption of the rubber absorbers can increase, and the kinetic energy of the steel slag particles also increases, resulting in the increase in friction heat energy.

3.5. Optimization Results

Combining with the above sections, it can be seen that when the crushing force achievement rate ηf is over 90%, the decrease in amplitude As and average power Pa are very small with increasing FM. Table 1 shows the FM is 20 t and the MM is 5.5 t for the industrial inertia cone crusher (model: GYP1200), and the operative crushing force Fo is 697.92 kN from the simulated experiment (Table A1). If the prediction values of Fo and ηf are 697.92 kN and 90%, the MM is set to 4.95 t (Figure 8a). According to Equation (13), the optimized fixed cone mass is 30.54 t.
According to Equations (14) and (15), the optimized amplitude As and average power Pa are 5.29 mm and 125.38 kW, respectively. Compared with the simulated experiment (Table A1), it can be seen that the optimized amplitude, average power and MM are decreased by 37.1%, 33.1% and 10.2%, respectively. Finally, we can see that the decrease in amplitude can effectively decrease the average power, and increase the crushing force achievement rate. However, the optimized FM is increased by 52.7%, so the optimized mass of the inertia cone crusher is about three times as much as the hydraulic cone crusher or spring crusher for the same industrial scale, which increases the manufacturing cost. In this paper, we design a more effective dynamic balancing mechanism for the inertia cone crusher, which decreases the amplitude and minimizes the mass of the inertia cone crusher.

4. Design of Dynamic Balancing Mechanism

4.1. Mechanics Principle of Dynamic Balancing

The total crushing force of the single exciter GYP-type inertia cone crusher is provided by the eccentric vibrator and the mantle, as shown in Equation (16):
F cr = F ice + F ucve
where Fcr is the total crushing force, Fice is the equivalent centrifugal force generated by the mantle, and Fucve is the equivalent exciting force generated by the eccentric vibrator.
The inertia cone crusher with the dynamic balancing mechanism is shown in Figure 11a. The dynamic balancing mechanism is mainly composed of a balancer and feedback mechanism. A balancer is added to the GYP-type inertia cone crusher, which can similarly rotate on the opposite side of the vibration exciter. In this way, the vibration of the crusher can be minimized to decrease the mass of the fixed and moving cone. Through a feedback mechanism to increase the crushing force, the kinetic energy of the balancer can be utilized efficiently.
The planar layout of forces acting on the mantle cone is shown in Figure 11b. Fic, Fucv and Fbec are the centrifugal forces generated by the mantle, eccentric vibrator and balancer, respectively.
F ic = m ic ω 2 e i ;   F ucv = m ucv ω 2 e u ;   F bec = m bec ω 2 e b
where mic, mucv and mbec are the mass of the mantle, eccentric vibrator and balancer, respectively; eic, eucv and ebec are the eccentric distance of the mantle, eccentric vibrator and balancer, respectively; ω is the drive shaft speed. The value of Fbec should conform to Equation (18).
F uic = m ucv + m ic ω 2 e ui = F bec
where Fuic is the equivalent force of the centrifugal force generated by the moving cone, and eui is the equivalent eccentric distance of the moving cone. Based on the lever principle, the vector force Fbeci can be expressed as
F cr + F beci + F uic + R 1 + G c = 0
where Fbeci is the force acting on the moving cone by the feedback mechanism, R1 is the constraint reaction of spherical joint O1, and Gc is the gravity of the moving cone.
Compared with Equations (16) and (19), it can be seen that the inertia cone crusher with the dynamic balancing mechanism adds the feedback force compared with the single exciter crusher. Therefore, the dynamic balancing mechanism can obviously decrease the amplitude of the fixed cone and increase the crushing force. In this way, the mass of the fixed and moving cone is minimized, and the manufacturing cost decreases significantly.

4.2. Elementary Prototype of Laboratory Experiments

4.2.1. Experimental Devices

Corresponding experiments are validated to verify the dynamic balancing mechanism. However, it is an impossible task for us to manufacture an industrial-scale inertia cone crusher with the dynamic balancing mechanism. In this paper, a laboratory prototype with the same dynamic balancing mechanism is developed. The amplitude, power draw and product size distribution of the laboratory crusher can be collected by some experimental devices, and the results are compared with the crusher without the dynamic balancing mechanism. The crusher without the dynamic balance can provide the same theoretical crushing force as the dynamic balance prototype. The experimental devices and dynamic signal acquisition systems are shown in Figure 12. The feed material is a 7.5–10 mm white marble, and the amplitude, power draw and product size distribution for the two driving speed levels (450 and 650 rpm) are compared in the following sections.

4.2.2. Amplitudes of Test Points

The displacements of two test point are sampled by displacement sensors in Figure 12a. The experimental data of two test points in the x direction are displayed for the drive speeds of 450 rpm and 650 rpm, as shown in Figure 13. In Table 3, the test data from the balancing crusher and without balancing mechanism are summarized. The results show that the amplitude and deflection angle of the balancing crusher are decreased by 80.6% and 64.2%, compared with the crusher without a balancing mechanism. Therefore, the good vibration reduction performance of the dynamic balancing mechanism is verified by experimental comparison.

4.2.3. Power draw and Product Size Distribution

The input power of the motor is sampled by an electrical parameter test instrument in Figure 12b. The experimental data are displayed for the drive speeds of 450 rpm and 650 rpm, as shown in Figure 14. Figure 14 shows that the average power of the balancing crusher is decreased by 20.9%, compared with the crusher without a balancing mechanism. Therefore, the dynamic balancing mechanism can effectively reduce the energy consumption and running cost of the inertia cone crusher.
The product size distribution is collected by square sieves in Figure 12c. The experimental data are displayed for the drive speeds of 450 rpm and 650 rpm, as shown in Figure 15. The product size distribution displays a relatively good correspondence between the balancing mechanism and without the balance. It can be seen that the mechanism realizes the purpose of utilizing the inertia force and kinetic energy of the balancer.

4.3. Optimization Verification of Industrial-Scale Inertia Cone Crusher

To verify the optimization performances and two mass variables (FM and MM), the simulated experiments using MBD–DEM coupling are performed on the industrial-scale inertia cone crusher with the dynamic balancing mechanism. The fixed cone mass and the moving cone mass (including the dynamic balancing mechanism mass) are 6.44 t and 3.82 t, respectively. Compared with Section 3.4, the optimized crushing force achievement rate is over 95% using with the dynamic balancing mechanism. Furthermore, the optimized amplitude and average power are decreased by 33.1% and 10.2%, respectively, as shown in Figure 16.
In Figure 16c shows the product size distributions for the case balancing mechanism and the case without are compared. It can be seen that the balance case is slightly finer than the without-balance case. Due to the amplitude, the average power and product size distribution have been improved, and the good crushing performance of the dynamic balancing mechanism is verified in the industrial-scale inertia cone crusher. Furthermore, the fixed cone mass and the moving cone mass are decreased by 78.9% and 22.8%, respectively. The dynamic balancing mechanism significantly reduces the manufacturing cost.

5. Conclusions

In the inertia cone crusher design, an inevitable problem concerns how to determine the two mass variables (fixed cone mass and moving cone mass), which can affect the crusher dynamic performances. Firstly, the crushing process of the inertia cone crusher is simulated using the MBD–DEM coupling. Predictive regression models, in which the two mass variables are taken as influencing factors, and the crushing force achievement rate, amplitude and average power are taken as the performance goals, are explored by the design of simulation experiments. In addition, it is found that when the achievement rate ηf is over 90%, the decrease in amplitude As and average power Pa, and the increase in ηf, are very small with increasing FM. Due to the optimized values of ηf and FM being 90% and 30.54 t, the optimized amplitude As, average power Pa, and MM are decreased by 37.1%, 33.1% and 10.2%, respectively, compared with the original crusher.
The optimized FM is increased by 52.7%, which increases the manufacturing cost. In this paper, a new and more cost-effective dynamic balancing mechanism of inertia cone crusher is achieved in order to reduce FM. The vibration reduction and inertia force utilization of the dynamic balancing mechanism are verified by the elementary prototype of the laboratory experiment. Moreover, the effect of FM and MM reduction is verified by the MBD–DEM simulation of the industrial inertia cone crusher. Compared with the without-balance case (Section 3.4), the amplitude, average power and product size distribution have been improved, and the FM and MM are decreased by 78.9% and 22.8%. As such, from a manufacturer’s perspective, the manufacturing cost decreases significantly. In order to reduce the running costs for users, the future work will prioritize manually using the design of simulation experiments around the optimum drive shaft speed, the eccentric distance of eccentric vibrator and the discharge gap.

Author Contributions

Conceptualization, J.C. and T.R.; Methodology, J.C.; Software, J.C.; Validation, Z.Z.; Formal analysis, D.L.; Investigation, J.C. and X.J.; Data curation, J.C.; Writing—original draft, J.C.; Project administration, T.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Key Technology R&D Program of China (grant number.2011BAF15B01).

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Simulated experiment scheme of two mass variables.
Table A1. Simulated experiment scheme of two mass variables.
Simulation RunFM—m1 (t)MM—m2 (t)Achievement Rate—ηf (%)Operative Force—Fo (kN)Amplitude—As (mm)Deflection Angle—γd (deg)Average Power—Pa (kW)
1113.185.79418.437.090.213105.39
2113.582.96455.657.850.239118.06
3114.081.97515.248.750.269154.19
4114.580.04562.989.630.301171.62
5115.079.38623.8110.590.329199.30
6115.575.76653.8511.360.358231.51
7115.973.83697.7911.950.383248.23
8153.188.75432.206.310.17786.16
9153.586.36472.117.030.198100.78
10154.085.34536.387.780.225133.59
11154.584.35593.228.640.254148.01
12155.081.94643.829.390.280168.72
13155.578.48666.1510.160.304207.90
14155.977.67723.9510.770.326229.78
15203.190.87443.195.230.14579.47
16203.589.18491.925.860.16488.69
17204.088.13553.816.550.186121.53
18204.587.67616.627.220.213135.34
19205.085.67673.257.910.231149.18
20205.581.87697.928.410.254187.17
21205.980.71743.239.080.272202.72
22253.192.94452.844.210.12173.14
23253.591.15502.394.760.13576.02
24254.090.29570.545.320.155105.89
25254.589.87643.465.990.176122.67
26255.088.31694.316.450.193138.81
27255.584.91726.747.050.231165.17
28255.987.01779.007.480.229184.56
29303.194.36459.723.470.10268.53
30303.592.97510.623.910.11572.56
31304.092.12578.334.400.13193.88
32304.591.48645.524.810.149113.45
33305.089.96706.505.350.163129.62
34305.587.39754.205.830.180149.31
35305.986.84799.676.210.193155.07
36353.195.48465.182.830.08764.54
37353.593.80515.133.080.09469.68
38354.092.91585.443.590.11382.87
39354.592.96653.763.980.128109.54
40355.091.02716.724.460.141121.50
41355.591.04785.694.660.155134.31
42355.989.98827.124.910.168136.34
43393.196.68470.812.470.07962.20
44393.594.36520.172.770.88067.96
45394.094.13591.513.090.10078.30
46394.593.44657.933.390.114105.51
47395.092.38725.883.770.125114.60
48395.592.08802.483.910.138128.06
49395.990.88836.754.430.149129.78
Table A2. The results of ANOVA for prediction regression coefficients.
Table A2. The results of ANOVA for prediction regression coefficients.
Model Regression Coefficient—BStandardization Coefficient—Bet Valuep ValueConfidence Interval for B
LowerUpper
fηFM0.6281.0977.718<0.010.4640.792
MM −0.400−1.9430.59
FM·FM−0.011−0.937−8.080<0.01−0.013−0.08
MM·MM−0.557−0.878−17.980<0.01−0.619−0.494
FM·MM0.0830.7597.940<0.010.0620.104
Constant82.196 82.400<0.0180.18584.206
fAFM−0.197−0.762−18.703<0.01−0.225−0.169
MM2.1620.83658.740<0.012.0632.261
FM·FM0.0030.58317.664<0.010.00250.0034
MM·MM −0.131−2.3310.025
FM·MM−0.039−0.799−28.692<0.01−0.043−0.036
Constant3.708 19.341<0.013.1924.224
fPFM−1.679−0.340−2.6410.011−2.960−0.398
MM64.1071.29928.872<0.0159.63268.582
FM·FM0.0680.6976.689<0.010.0480.088
MM·MM 0.1740.9350.355
FM·MM−0.979−1.039−11.804<0.01−1.146−0.812
Constant−56.084 −4.850<0.01−79.391−32.777

References

  1. Safronov, A.N.; Kazakov, S.V.; Shishkin, Y.V. New trends in inertia cone crusher studies. Miner. Process. J. 2012, 5, 40–42. [Google Scholar]
  2. Babaev, R.M.; Kishchenko, V.L. The effect of KID-1500 inertia cone crusher parameters upon preset crushed stone size fraction yield. Miner. Process. J. 2012, 1, 22–24. [Google Scholar]
  3. Savov, S. Constructive-mechanical review of conical inertial crushers (KID). Bulg. J. Eng. Des. 2012, 15, 23–28. [Google Scholar]
  4. Savov, S.; Nedyalkov, P.; Minin, I. Crushing force theoretical examination in one cone inertial crusher. J. Multidiscip. Eng. Sci. Technol. 2015, 2, 430–437. [Google Scholar]
  5. Xia, X.; Chen, B. Study on dynamic characteristics of inertia cone crusher and its application. In Proceedings of the XXV International Mineral Processing Congress (IMPC 2010), Brisbane, Australia, 6–10 September 2010; pp. 1399–1403. [Google Scholar]
  6. Cleary, P.W.; Sinnott, M.D.; Morrison, R.D.; Cummins, S.; Delaney, G.W. Analysis of cone crusher performance with changes in material properties and operating conditions using DEM. Miner. Eng. 2017, 100, 49–70. [Google Scholar] [CrossRef]
  7. Andre, F.P.; Tavares, L.M. Simulating a laboratory-scale cone crusher in DEM using polyhedral particles. Powder Technol. 2020, 372, 362–371. [Google Scholar] [CrossRef]
  8. Chen, Z.; Wang, G.; Xue, D.; Bi, Q. Simulation and optimization of gyratory crusher performance based on the discrete element method. Powder Technol. 2020, 376, 93–103. [Google Scholar] [CrossRef]
  9. Cheng, J.; Ren, T.; Zhang, Z.; Liu, D.; Jin, X. A dynamic model of inertia cone crusher using the discrete element method and multi-body dynamics coupling. Minerals 2020, 10, 862. [Google Scholar] [CrossRef]
  10. Xu, L. An approach for calculating the dynamic load of deep groove ball bearing joints in planar multibody systems. Nonlinear Dyn. 2012, 70, 2145–2161. [Google Scholar] [CrossRef]
  11. Langerholc, M.; Cesnik, M.; Slavic, J.; Boltear, M. Experimental validation of a complex, large-scale, rigid-body mechanism. Eng. Struct. 2012, 36, 220–227. [Google Scholar] [CrossRef]
  12. Li, C.; Honeyands, T.; O’Dea, D.; Moreno-Atanasio, R. The angle of repose and size segregation of iron ore granules: DEM analysis and experimental investigation. Powder Technol. 2017, 320, 257–272. [Google Scholar] [CrossRef]
  13. Cundall, P.A.; Strack, O.D. A discrete numerical model for granular assemblies. Geotechnique 1979, 30, 331–336. [Google Scholar] [CrossRef] [Green Version]
  14. Barrios, G.K.P.; Tavares, L.M. A preliminary model of high pressure roll grinding using the discrete element method and multi-body dynamics coupling. Int. J. Miner. Process. 2016, 156, 32–42. [Google Scholar] [CrossRef]
  15. Chung, Y.C.; Wu, C.W.; Kuo, C.Y.; Hsiau, S.S. A rapid granular chute avalanche impinging on a small fixed obstacle: DEM modeling, experimental validation and exploration of granular stress. Appl. Math. Model. 2019, 74, 540–568. [Google Scholar] [CrossRef]
  16. Znamenskll, S.B.; Krivelev, D.M. Inertial Cone Crusher for Comminution of Ores of Varying Hardness Has Base with Cylindrical Vertical Support Element, on Which Driving Counter-Balance Is Mounted on Bearings, Counter-Balance Including Concentric Portion with V Belt Groove. Russia Patent No. 2097132-C1, 20 February 1996. [Google Scholar]
  17. Ren, T.; Huang, K.; Cheng, J.; Zhang, Z.; Jin, X. Vibration performance on novel inertia cone crushers. Chin. Mech. Eng. 2019, 30, 2704–2708. [Google Scholar]
  18. Orin, D.E.; McGhee, R.B.; Vukobratovic, M.; Hartoch, G. Kinematic and kinetic analysis of open-chain linkages utilizing Newton-Euler methods. Math. Biosci. 1972, 43, 107–130. [Google Scholar] [CrossRef]
  19. Potyondy, D.O.; Cundall, P.A. A bonded-particle model for rock. Int. J. Rock Mech. Min. Sci. 2004, 41, 1329–1364. [Google Scholar] [CrossRef]
  20. Quist, J.; Evertsson, C.M. Cone crusher modelling and simulation using DEM. Miner. Eng. 2016, 85, 92–105. [Google Scholar] [CrossRef]
  21. Wang, D.; Servin, M.; Berglund, T.; Mickelsson, K.O.; Ronnback, S. Parametrization and validation of a nonsmooth discrete element method for simulating flows of iron ore green pellets. Powder Technol. 2015, 283, 475–487. [Google Scholar] [CrossRef]
  22. Mindlin, R.D. Compliance of elastic bodies in contact. J. Appl. Mech. 1949, 16, 259–268. [Google Scholar]
  23. Ren, T.; Huang, K.; Cheng, J.; Zhang, Z.; Liu, D.; Jin, X. Multi-body dynamics Analysis and test verification for inertia cone crushers. Chin. Mech. Eng. 2020, 31, 2322–2331. [Google Scholar]
  24. Myers, R.H.; Montgomery, D.C. Response Surface Methodology: Process and Product Optimization Using Designed Experiments, 2nd ed.; Wiley: New York, NY, USA, 2008. [Google Scholar]
Figure 1. Schematics of the inertia cone crusher: (a) vertical cross-section and (b) simplified multi-body dynamic (MBD) with bonded particles.
Figure 1. Schematics of the inertia cone crusher: (a) vertical cross-section and (b) simplified multi-body dynamic (MBD) with bonded particles.
Minerals 11 00163 g001
Figure 2. Schematic illustration of particle bond: (a) before a parallel bonding beam formation and after breakage two particles interact according to the Hertz–Mindlin contact model, and (b) two particles bonded together with a parallel beam according to the bonded particle model (BPM).
Figure 2. Schematic illustration of particle bond: (a) before a parallel bonding beam formation and after breakage two particles interact according to the Hertz–Mindlin contact model, and (b) two particles bonded together with a parallel beam according to the bonded particle model (BPM).
Minerals 11 00163 g002
Figure 3. The schematic illustration of the slag particle using DEM: (a) BPM of the slag particle is formed by the EDEM software, and (b) a realistic slag shape is used creating the packing structure of a normal distribution.
Figure 3. The schematic illustration of the slag particle using DEM: (a) BPM of the slag particle is formed by the EDEM software, and (b) a realistic slag shape is used creating the packing structure of a normal distribution.
Minerals 11 00163 g003
Figure 4. The photo of a steel slag subjected to (a) the uniaxial deformation and (b) the repose angle test.
Figure 4. The photo of a steel slag subjected to (a) the uniaxial deformation and (b) the repose angle test.
Minerals 11 00163 g004
Figure 5. The photo of a steel slag subjected to (a) Brazilian test and (b) BPM calibration simulation.
Figure 5. The photo of a steel slag subjected to (a) Brazilian test and (b) BPM calibration simulation.
Minerals 11 00163 g005
Figure 6. Flowchart of a coupled MBD–DEM model.
Figure 6. Flowchart of a coupled MBD–DEM model.
Minerals 11 00163 g006
Figure 7. Series of images from MBD–DEM simulations showing the bond cluster representations.
Figure 7. Series of images from MBD–DEM simulations showing the bond cluster representations.
Minerals 11 00163 g007
Figure 8. Regression curves of the GYP inertia cone crusher for (a) the relationship between the theoretical force and MM, (b) the relationship between the force achievement rate and FM, (c) the relationship between the force achievement rate and MM, and (d) the influence of the interaction between FM and MM on the force achievement rate using response surface methodology (RSM).
Figure 8. Regression curves of the GYP inertia cone crusher for (a) the relationship between the theoretical force and MM, (b) the relationship between the force achievement rate and FM, (c) the relationship between the force achievement rate and MM, and (d) the influence of the interaction between FM and MM on the force achievement rate using response surface methodology (RSM).
Minerals 11 00163 g008
Figure 9. Amplitude regression curves of the GYP inertia cone crusher for (a) the relationship between amplitude and FM, (b) the relationship between amplitude and MM, and (c) the influence of the interaction between FM and MM on amplitude using RSM.
Figure 9. Amplitude regression curves of the GYP inertia cone crusher for (a) the relationship between amplitude and FM, (b) the relationship between amplitude and MM, and (c) the influence of the interaction between FM and MM on amplitude using RSM.
Minerals 11 00163 g009
Figure 10. Average power regression curve of the GYP inertia cone crusher for (a) the relationship between the average power and FM, (b) the relationship between the average power and MM, and (c) the influence of the interaction between FM and MM on the average power using RSM.
Figure 10. Average power regression curve of the GYP inertia cone crusher for (a) the relationship between the average power and FM, (b) the relationship between the average power and MM, and (c) the influence of the interaction between FM and MM on the average power using RSM.
Minerals 11 00163 g010
Figure 11. The dynamic balancing inertia cone crusher: (a) vertical cross-section and (b) mechanical analysis of main mechanism.
Figure 11. The dynamic balancing inertia cone crusher: (a) vertical cross-section and (b) mechanical analysis of main mechanism.
Minerals 11 00163 g011
Figure 12. The photo of laboratory prototype experiments for (a) experimental devices, (b) signal acquisition systems and (c) mechanical analysis of main mechanism.
Figure 12. The photo of laboratory prototype experiments for (a) experimental devices, (b) signal acquisition systems and (c) mechanical analysis of main mechanism.
Minerals 11 00163 g012
Figure 13. Comparison of the displacements between balancing crusher and without balancing mechanism for (a) 450 rpm case, and (b) 650 rpm case.
Figure 13. Comparison of the displacements between balancing crusher and without balancing mechanism for (a) 450 rpm case, and (b) 650 rpm case.
Minerals 11 00163 g013
Figure 14. Comparison of the input powers between balancing crusher and without balancing mechanism (a) 450 rpm case, and (b) 650 rpm case.
Figure 14. Comparison of the input powers between balancing crusher and without balancing mechanism (a) 450 rpm case, and (b) 650 rpm case.
Minerals 11 00163 g014
Figure 15. Comparison of the product size distribution between balancing mechanism and without balance (a) 450 rpm case, and (b) 650 rpm case.
Figure 15. Comparison of the product size distribution between balancing mechanism and without balance (a) 450 rpm case, and (b) 650 rpm case.
Minerals 11 00163 g015
Figure 16. Comparison of the performance of industrial-scale inertia cone crusher between balance and without dynamic balance for (a) amplitude case, (b) power draw case, and (c) product size distribution case.
Figure 16. Comparison of the performance of industrial-scale inertia cone crusher between balance and without dynamic balance for (a) amplitude case, (b) power draw case, and (c) product size distribution case.
Minerals 11 00163 g016
Table 1. Simulation parameters for MBD–DEM.
Table 1. Simulation parameters for MBD–DEM.
ParameterValueUnit
DEM material properties
RockSteel
Solid density47007800(kg/m3)
Shear stiffness1.48∙1097.0∙1010(Pa)
Poisson’s ratio0.380.3
Rock_RockRock_Steel
Coefficient of static friction0.560.7
Coefficient of restitution0.150.25
Coefficient of rolling friction0.010.01
BPM parameters
Normal stiffness556(GPa/m)
Shear stiffness250(GPa/m)
Normal critical stress32(MPa)
Shear critical stress8.5(MPa)
Bond disc radius3.2(mm)
Machine
Mantle cone angle50(deg)
Driving speed550(rpm)
Fixed cone mass20,000(kg)
Moving cone mass5500(kg)
Rubber absorber properties
Stiffness coefficient kx,ky,kz350,350,970(N/mm)
Damping coefficient cx,cy,cz20,20,40(N·s/mm)
Table 2. The results of ANOVA for predictive regression model.
Table 2. The results of ANOVA for predictive regression model.
ModelDegree FreedomMean SquareF Valuep ValueDetermination Coefficient
fηRegression4363.167767.995<0.010.993
Error440.473
Total48
fARegression475.0109625.911<0.010.999
Error440.008
Total48
fPRegression427,049.462954.059<0.010.989
Error4428.352
Total48
Table 3. The summary of experimental averaged results.
Table 3. The summary of experimental averaged results.
Experiment CaseAmplitude of First Point (mm)Amplitude of Second Point (mm)Amplitude of Rotation Point (mm)Deflection Angle (mm)
450 rpmBalance0.180.550.040.18
Without0.492.110.220.63
650 rpmBalance0.210.660.050.24
Without0.552.470.250.67
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cheng, J.; Ren, T.; Zhang, Z.; Jin, X.; Liu, D. Influence of Two Mass Variables on Inertia Cone Crusher Performance and Optimization of Dynamic Balance. Minerals 2021, 11, 163. https://doi.org/10.3390/min11020163

AMA Style

Cheng J, Ren T, Zhang Z, Jin X, Liu D. Influence of Two Mass Variables on Inertia Cone Crusher Performance and Optimization of Dynamic Balance. Minerals. 2021; 11(2):163. https://doi.org/10.3390/min11020163

Chicago/Turabian Style

Cheng, Jiayuan, Tingzhi Ren, Zilong Zhang, Xin Jin, and Dawei Liu. 2021. "Influence of Two Mass Variables on Inertia Cone Crusher Performance and Optimization of Dynamic Balance" Minerals 11, no. 2: 163. https://doi.org/10.3390/min11020163

APA Style

Cheng, J., Ren, T., Zhang, Z., Jin, X., & Liu, D. (2021). Influence of Two Mass Variables on Inertia Cone Crusher Performance and Optimization of Dynamic Balance. Minerals, 11(2), 163. https://doi.org/10.3390/min11020163

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop