Design, Development and Validation of An Impact Wear Testing Device For The Minerals Industry
Design, Development and Validation of An Impact Wear Testing Device For The Minerals Industry
Design, Development and Validation of An Impact Wear Testing Device For The Minerals Industry
Lokesh K. Thakur & Peter Radziszewski Mechanical Engineering Department McGill University, Montreal, Quebec H3A 2K6 e-mail : lokesh.thakur@mail.mcgill.ca
Abstract
Total steel media wear in a given mill (ball or SAG) grinding process is a product of three recognized wear mechanisms-impact, abrasion, and corrosion of which the contribution of each wear mechanism to total media wear has not been well established. A total media wear model can be dened on the assumption that the effect of each wear mechanism can be independently determined and this effect can be tied to mill charge motion as determined or estimated using a charge motion simulator, which allows for a total media wear model to be dened as the summation of the wear results of each mechanism. This necessitates the need for developing impact, abrasion and corrosion wear testers that will allow the study of media wear over a wide range of energy levels. Renement of the impact test will allow studying how impact media wear at high energies behaves as a function of those energies. This work is focused on design and development of a high-energy impact wear tester. A comprehensive illustration on designing the newest version of the impact wear tester is explained in detail in order to exhibit how theoretical, virtual and experimental analyses could be integrated while designing a mechanical system.
as a function of the operating parameters of a given mill [2, 3]. A tumbling mill, whether it is autogenous, semiautogenous, ball or rod, is a system comprises of a number of interrelated and interactive elements that work together in order to grind a given ore. This comminution process is achieved by individual balls that compose actual ball mill elements, which bring about ore breakage [4]. Jointly, these balls form the mill ball charge which, during ball mill operation, typically has a charge prole as can be found in the gures 1.
Figure 1: Typical ball charge motion prole. A ball charge motion prole shows three comminution zones (gure 1) that can be distinguished by (i) ball layers sliding over one another grinding material trapped between them, the grinding zone, (ii) balls rolling over one another breaking material in low energy impact, the tumbling zone and by (iii) ball in ight re-entering into the ball charge crushing material in high energy impact, the crushing zone. A total media wear model can be dened on the assumption that the effect of each wear mechanisms can be independently determined and this effect can be tied to mill charge motion as determined or estimated using a charge motion simulator which allows for a total media wear model to be dened as the summation of the wear result of the wear each mechanism [5].
3
1. Introduction
In mineral processing industries, grinding processes can be dened as pulverization or comminution of big particles in order to reduce them to minute particles. Depending upon the mining operation, half of the overall operating cost of comminution processes can be found in the cost of steel media wear. Metal wear is primarily associated with these processes. Grinding performance and quality can also be affected by wear, which consists of three recognized wear mechanisms-abrasion, corrosion and impact. In such a context, predictive wear models have become a necessity to determine most optimum grinding conditions that can reduce process wear while maintaining grinding performance and quality [1]. These impact and abrasion energy spectra provide an opportunity to study media wear
m total =
i=1
m i
(1)
where i=1,2,3 represents the abrasion, corrosion and impact wear mechanisms respectively. It is assumed that the contribution of each wear mechanism in equation 1 can be determined using the impact/abrasion energy spectra along with the lab test results. These lab tests were developed from some existing tests. Abrasive Wear: The abrasive wear characteristics, as a function of applied abrasion force, were determined using steel-on-steel abrasion wheel test where the abrasive was a ground ore. Corrosive Wear: The raw corrosion wear characteristics were determined from a batch test using the same ground ore and mill water in order to replicate a similar corrosion environment. Impact Wear: Tumbling and free falling grinding media explain the phenomena of impact associated with the energies dissipated in tumbling and crushing zones of the ball charge motion prole (gure 1). Three impact testers found in the literature [6, 7, 8] present promising aspect for determining grinding media impact wear. The impact chamber principle shown by the Scieszka and Dutkiewicz tester [8] presented the most exible impact wear/impact energy results for different ore-metal-environment combinations. Based on these analyses, it is possible to reformulate the equation 1 for a similar ore-metal-environment industrial context to,
nabr
int
test must be developed. A new impact tester has been developed that would generate impacts in the order of those found in the USBM media drop test [11]. The design criteria for the newer test is exibility in impact energies and impact velocities. This has led to the development of a precompressed spring cam assembly that can be schematized as in the gure 2.
mtotal =
i=1
+
j =1
As mentioned in the P9L study on exploring total media wear [10], it has been observed that the laboratory set-up although allowing impact variation, did not attain the desired impact energy levels. As the impact chamber can not generate impact energies greater than that equal to a 4 cm drop, a more aggressive 2
Figure 3: Version I of the impact tester assembly [12]. Initial Test Results: A few tests were carried on the version I of the high energy impact tester in order to vali-
date its performance against predened set of goals, which were exclusively focused on achieving maximum speed at the impact. High resolution slotted optical switches were used for measuring speed at the impact. Following assumptions taken during Observations [12]: (1) the system was assumed to be a frictionless system and (2) acceleration gradient was assumed to be equal to zero. A medium pressure die spring (color coded blue) with a spring constant () = 29.427 N/mm [13] was used. The coupling was positioned on the vertical arm in such way that the spring compression within the efcient operating range (25 %-35% of free length) could be achieved [13]. The spring was compressed for approximately 4 or 101 mm, which is 33.33% of its free length (12 or 304.8 mm). This version of the high energy impact could only produce impact energy equivalent to 9.8 cm of drop height. Following errors were found in the results: (1) optical switches were very sensitive to high amplitude of vibrations produced by the mechanical system during its operation and (2) electronic counter couldnt be calibrated properly, this is why its readings could have had moderate degree of error. Because of above mentioned errors in the observation, optical switches couldnt be used for the further measurements. This is why a high-speed camera was used for measuring speed at the impact while running tests on the version II of the high energy impact wear tester.
Figure 5: Path traced by cam-follower contact point [12]. forward direction along with the Lever arm assembly under the inuence of potential energy delivered by the compressed spring. This assembly moves until the caret, which is assembled at the forefront of mobile assembly, collides with the stationary target ball. As shown in the gure 5, the cam-follower contact point moves back and forth between locations O1 and O2 .
Derivation for Calculating Spring Displacement: In the triangle O0 , O3 and O2 (gure 5) ( O2 O0 O3 ) can be calculated as [15], (O3 O2 )2 = (O2 O0 )2 + (O0 O3 )2 2(O2 O0 )(O0 O3 ) cos (4) For a small angular displacement, in the triangle O0 , O1 and O2 , Dmax or maximum displacement of the spring for a radius rmax (gure 5) can be calculated as [15], Dmax = 2 sin(/2)rmax (5)
Since the value of is constant for a given set of lengths therefore a formula can be derived as [12], Di = kri Figure 4: Cam prole in mm [14]. As shown in the gures 4 and 5, contact point at the follower (location O1 ) rotates at a radius of rmax from the pivot point O. The contact point at the cam (location O1 ), which is situated at a distance of rmax from the pivot point Oo , rotates from the initial value of radius (59.20 mm) to the maximum value of radius (6 or 152.40 mm). After leaving the cam (location O2 ) , follower moves rapidly in the 3 i = 1, 2, 3, ... (6)
where ri is a variable radius dened as the distance between coupling axis and the axis of rotation, which passes through the pivot point O0 and k is a constant, which depends on dimension of the Lever arm and the cam [12]. Derivation for Calculating Speed at the Impact: Equation for calculating speed at the impact can be derived from law of conservation of energy. In case of reciprocating impact mechanism potential energy stored in the spring media will be delivered to the reciprocating mass and the Lever arm, out of which some amount of energy will be lost
in overcoming the friction. Potential energy stored in the spring can be expressed as, Eps = 1 2 x 2 (7)
Di = displacement or compression amplitude of the spring in mm. Now, an expression for calculating speed at the impact can be written as [12], Eps + Emhc = Emrc + Emvc + Ef va + Ef rc + Eabs (13) where Eabs is the maximum energy absorbed by the Lever arm during loading conditions, which is equal to 1/2F fmax. Where F is the load on the Lever arm in N and fmax is the maximum deection of it under the inuence of force F in mm. The equation 13 can be rewritten as, Eps + 1000 1000 2 2 (Wla sin )vi = (Wrc + Wla cos )vi 2 2 (14)
where, Eps is the potential energy stored in the spring in Nmm or Nm. This energy will be converted into momentum energy and some part of it will be lost in overcoming the friction. Momentum energy (Emrc ) for the reciprocating components can be expressed as, Emrc = 1 (W )v 2 2 rc i (8)
where, vi = Speed at the impact in m/sec. Wrc = Total weight of the reciprocating components in kg. Momentum energy (Emvc ) for the rotating components (Lever arm assembly) can be expressed as [12], Emvc = 1 2 (W cos )vi 2 la (9)
For a particular set of conditions i.e., spring constant, spring displacement etc. x and Di would be equal, therefore equation 14 can be rewritten as, vi =
(r ) cos + W D E Di 2 2[dr Wla p ds rc i abs ] + W cos W sin ] 1000[Wrc la la (15)
where, Wla = Weight of the Lever arm assembly in kg. = Angular displacement of the Lever arm in radians. As we can see from the equation 9, only vertical component (Wla cos ) of the Lever arm assembly is considered due to the fact that only this component would oppose the motion. Horizontal inertial component of the Lever arm assembly will help the spring to restore its original length, therefore momentum energy (Emhc ) associated with this component will be added to the input energy to the system, it can be expressed as [12], Emhc = 1 2 (W sin )vi 2 la (10)
Part of the input energy lost in overcoming the friction (Ef va ) for the Lever arm assembly can be expressed as [12],
Ef va = dr Wla cos (rp )
(11)
where, dr = coefcient of dynamic rolling friction for the rotating mass (Lever arm), and rp = radius of the pivot pin about which Lever arm is hanged in mm, which is equal to 5/16 or 7.94mm. Part of the input energy lost in overcoming the friction (Ef rc ) for the reciprocating components assembly can be expressed as [12],
Ef rc = ds Wrc Di
(12)
where, ds = coefcient of dynamic sliding friction for the reciprocating components assembly, and 4
The Second Simulation Run: In the second simulation run, weight of all the moving components and friction coefcient (by assuming proper lubrication and improvement in the design) were halved except weight of the Lever arm assembly, which was reduced to one forth of its original value. The maximum speed at the impact in m/sec was found to be equal to 11.26 m/sec. Verication of the Virtual Model: If maximum radius of the cam is Rcam , an expression for calculating maximum torque (Tmax ) delivered by the cam under equilibrium condition can be written as [12], Tmax = r1 (x) Rcam r Figure 6: Previous Lever arm assembly. (16)
where Rcam is the maximum radius of the cam, which is equal to 6 or 152.40mm. The maximum value of torque calculated from the ADAMS simulation was 1375.14 N m and from the theoretical calculations was 1311.01 N m. The virtual model had an error of 4.89%, which is quite reasonable after considering the fact that errors are almost inevitable during the building of virtual prototype model from a physical prototype [12]. Virtual prototyping helped to know which variables are critical to the optimization. To dene variables, we need to know which variables have an effect on the variable to be minimized. The objective was to minimize volume of the moving components which is directly related to weight of the components. Further, results from ADAMS simulations were used for determining variables for optimizing designs while using ANSYS for the nite element analysis purposes.
Figure 7: Lever arm assembly: nal design [12]. movement. (3) Very easy to change the position of coupling on the Lever arm for altering spring compression amplitude. (3) Overall weight of it, if the same material used would be about 3.10 kg lighter than the previous Lever arm assembly design (gure 6). Material Selection: Titanium ASTM Grade 5, Ti - 6Al - 4V (alpha-beta alloy) was chosen as the material to manufacture Lever arm assembly because of the following reasons [18]: (1) The most easily available and the most widely used titanium alloy. (2) It has very high strength and excellent fatigue properties. (3) The alloy is weldable. (4) It has very good mechanical properties up to 550o C. (5) Its intrinsic resistance to corrosion practically eliminates maintenance in most environments. Ti-6Al-4V (ASTM Grade 5) has a fatigue strength of 490 MPa for a smooth specimen (Kt =1) under rotating bending condition for stress to cause failure in 107 cycles [18]. Maximum Allowable Design Stress: The safety factor for structural systems, proposed Philon of Byzantium (3rd century BC) [19], is dened as follows N= capacity strength = >1 load stress (17)
The value of a central safety factor (equation 17) should not be less than two (2) for most structural applications and
should routinely be set at three (3). Uncertainty in loading, uncertainty in material properties, foreseeable abuse, and challenging service environments demand higher values of the safety factor. A long service life also requires a larger value of safety factor. High reliability applications require systems with a larger central safety factor value. Cyclic loading induces fatigue failure in structural components. The safety factor for primarily cyclic loading should be based on the endurance limit rather than on yield or tensile strength [20]. After taking consideration of uncertainties and assumptions associated with the overall reliability of the design; a value of three (3) as a FOS can reasonably be assumed. From the equation 17, maximum allowable design stress can be calculated as [12], Design Stress = 490 = 163.33 M P a 3 (18)
Figure 8: Shape factors for I-section [12]. tational nite element analyses (FEA) were used. ANSYS 7.0 (2002) was used comprehensively for this purpose.
where fatigue or endurance strength of the Titanium Grade 5, Ti - 6Al - 4V (alpha-beta alloy) is 490.0 MPa [18]. As on results, the use of Titanium ASTM Grade 5 further lightens the weight of Lever arm assembly to approximately 5.0 kg [12].
cause Titaniums manufacturing cost is higher than the material cost itself. This is why a new I beam cross section was adapted according to the available Titanium sheets dimensions for reducing the manufacturing cost after consulting a material supplier.
Figure 9: Finally adapted I beam section in mm [12]. The gure 10 shows results obtained from the nal FEA for the I section showed in the gure 9. Figure 11: An exploded view of the coupling assembly [12]. Design of the Front Part: Boundary conditions for the front part were determined by assuming that the spring force would be transmitted to the part assembled in series. Which means that two front anges of the front part would be under the inuence of a force, tensile in nature and 12, 700 N in magnitude [12]. Determination of Dimensions by Space Planning: A simple FEA was performed for determining thickness of the horizontal link. In order to reduce the manufacturing and the material cost associated with Titanium (Ti6Al-4V), high quality regular structural steel with fatigue strength=500 MPa (107 cycles) [23] was used to manufacture the horizontal link (gure 11). It was found out that a thickness of 13.00 mm would be sufcient to ensure no failure (FOS = 3.0) under the inuence of 12, 700 N of tensile force [12]. Several materials were considered for manufacturing the front part pin and the bolts, which could be used in assembling front part with the side brackets and the horizontal link (gure 11). It was found out that hex socket head shoulder screws of 5/8 or 15.88 mm diameter made from grounded and heat treated stainless steel have a minimum single shear strength of 28, 785 lbs (128, 070 N) at its shoulder part [25]. A few thorough nite element analyses were performed in order to determine remaining dimensions of the front part. Following gure shows an FEA result from ANSYS for the nally designed front part. Design of the Side Brackets: Space planning in conjunction with ANSYS FEA was used to reach the nal design of the side brackets. Assuming the symmetrical distribution of the force exerted on the left and right side brackets due to the displacement of the spring, it can be stated that the force 7
Figure 10: FEA of the nally adapted I beam section [12]. Validation of FEA Results with Theoretical Analysis: Lever arm design was used for ANSYS simulation validation purposes due to the fact that only the I beam analysis was carried out comprehensively at the theoretical level because of its regular structural geometry. And FEA was use to analyze designs of all the Lever arm assembly components comprehensively. I beam properties were calculated for the I cross section shown in the gure 9. Design of beams may be based on strength (stress) or on stiffness if deection must be limited. Maximum stress in MPa of an I beam can be calculated as, Smax = Mbmax 12700 563/4 = = 147.078 I/c 308700.613/25.4 (21)
As it can be seen from the gure 10, maximum stress in the beam calculated from ANSYS is 147.072 MPa, which is very close to the value of maximum stress calculated theoretically (equation 21).
follower assembly.
Figure 14: An exploded view of follower assembly [12]. Figure 12: Front part design: nal FEA [12]. will be equally distributed between them. A few thorough nite element analyses were performed in order to determine remaining dimensions of the side brackets. Following gure shows an FEA result from ANSYS for the nally designed side bracket. While determining dimensions of the main body, a few interrelated parameters like outer diameter of the follower roller, its position on the main body and the main bodys position on the Lever arm were determined in order to insure no contact between tip of the cam and the bottom part of the Lever arm, when the cam leaves the follower roller. A few trigonometrical geometries were sketched for determining the path traced by the cam-follower roller contact point. Later ADAMS simulation software in conjunction with ANSYS was used to simulate various position of the main body on the Lever arm in order to reduce stress concentration at the Lever arm-main body interface. It was found out that an outer diameter of 32.00 mm for the follower roller will help to insure no contact between the cam and the bottom part of the Lever arm [12]. A few thorough nite element analyses were performed in order to determine remaining dimension of the main body. The gure 15 shows an FEA result from ANSYS for the nally designed main body. Figure 13: Side bracket design: nal FEA [12].
Figure 16: A photo of the Lever arm assembly [12]. green) with a spring constant () = 124.717 N/mm [13] was used. The coupling was positioned on the vertical arm in such way that the spring compression within the efcient operating range (15%-25% of the free length) could be achieved [13]. The spring was compressed for 2.36 or 60.00 mm, which is 19.66% of its free length (12 or 304.80 mm) [12]. Test Results and Validation: A high speed camera with 1000 frames/sec was used to measure speed at the impact. The speed at the impact was found to be equal to 6.0 m/sec for a compression of 60.00 mm, which compared well with that calculated from theory for the same compression. Speed at the impact was calculated using equation 15 was found to be equal to 6.88 m/sec, which overestimates the observed impact speed by 12.80%, which is a reasonable error [12].
Table 2: Weights after upgrading impactor prototype [12]. The gure 16 shows an actual photograph taken after assembling the complete Lever arm assembly along with the horizontal link to the impact tester version II. After having all the parameters related to the impact testers dynamics e.g., dimensions and weights of the all the reciprocating components, the equation 6 for calculating spring displacement (Di ) was modied [12]. The value of angle was calculated for the Lever arm posture shown in the gure 5, the equation 4 gives = 13.21o. Now the equation 6 can be modied as, Di = 2 sin(13.21/2)ri = 0.23ri (22)
where is a correction factor, which was introduced after assembling the impactor prototype. Based on observations, a value of 0.735 was introduced as a correction factor. Tests on the Impactor Performance: For impact testers validation purposes, a high pressure die spring (color coded 9
Acknowledgments
This work is co-sponsored by AMIRA-JKMRC P9M project and a Natural Sciences and Engineering Research Council (NSERC) of Canada collaborative research and development matching grant.
[12] L. K. Thakur, Design and Development of an Impact Wear Testing Device for the Minerals Industry, Master of Engineering Thesis, Mecahanical Engineering Department, McGill University, Montreal, Quebec, 2004. [13] Producto, The Complete Catalogue of Springs for All Your Needs, Producto Diemakers Suppliers Limited, Montreal, Quebec, 1996, pp. 1-6. [14] M. Girard, Simulateur DImpact. Rapport De Conception, Universite du Que bec en Abitibi-Te miscamingue, Montr eal, Quebec, 2000. [15] H. P. Greenspan and D. J. Benney, 1973, Calculus. An Introduction to Applied Mathematics, McGraw-Hill, Inc., New York, 1973. [16] ADAMS/Solver, Using ADAMS/Solver, Automatic Dynamic Analysis of Mechanical Systems, MSC Software, Santa Ana, CA, 2001. [17] ADAMS/View 12.0 Help, Intranet HTML Help On Using ADAMS/View 12.0, Automatic Dynamic Analysis of Mechanical Systems, MSC Software, Santa Ana, CA, 2001. [18] S. D. Henry, G. M. Davidson, S. R. Lampman, F. Reidenbach, R. L. Boring and W.W. Jr. Scott, Fatigue Data Book: Light Structural Alloys, ASM International, OH, 1995, pp. 181-385 & pp. 1-110. [19] J. E. Shigley and C. R. Mischke, Mechanical Engineering Design, Sixth Edition, McGraw Hill, Inc., New York, 2001. [20] R. C. Rice, B. N. Leis and D. V. Nelson, Fatigue Design Handbook/prepared under the auspices of the Design Handbook Division of the SAE Fatigue Design and Evaluation Technical Committee, AE-10, Second Edition, Society of Automotive Engineers, Warrendale, PA, 1988. [21] T. Baumeister, E. A. Avallone and T. BaumeisterIII, Marks Standard Handbook for Mechanical Engineers, Eighth Edition, McGraw-Hill Book Company, New York, 1988, pp. 5.25.37. [22] V. Adams and A. Askenazi, Building Better Products With Finite Element Analysis, OnWord Press, Santa Fe, NM, 1999. [23] MetWeb, Material Property Data, Automation Creations, Inc., Blacksburg, VA, 2003, URL: http://www.matls.com/index.asp?ckck=1. [24] J. Fernandes, Technical Service & Sales Representative, Private Consulting, Elimetal, Inc., St. Laurent, Quebec, 2003. [25] Smith Fasteners, Socket Shoulder Screws Dimensions and Specications, Smith Fastener Company, Bell, CA, 2003, URL: http://www.smithfast.com/shoulderscrew.html. [26] Titanium Industries, Data & Reference Manual, Titanium Industries, Inc., Montreal, QC, 2003.
Greek Symbols
: Unit, angular; S : Coefcient of variation; : Angle; : Spring constant; : Correction factor; : Coefcient of friction; p : Poissons ratio; : Metal density; : Angular displacement.
References
[1] P. Radziszewski, Predictive Model For Ball Mill Wear, Canadian Metallurgical Quarterly, 1997a, Vol. 36, No. 2, pp. 87-93. [2] P. Radziszewski and S. Tarasiewicz, Modelling and Simulation of Ball Mill Wear, Wear, 1993a, Vol. 160, pp. 309-316. [3] P. Radziszewski and S. Tarasiewicz, Simulation of Ball Charge and Liner Mill, Wear, 1993b, Vol. 169, pp. 77-85. [4] P. Radziszewski, Exploring Total Media Wear, Mechanical Engineering Department Technical Report, McGill University, Montreal, Quebec, 2001. [5] P. Radziszewski, Chapter 4: Ball Mill Media Wear and Prediction, JKMRC/AMIRA P9L Technical Project Report, Queensland, Australia, 1997b, pp. 35-51. [6] Blickensderfer and Tylczak, Evaluation of Commercial US Grinding Balls by Laboratory Impact and Abrasion Tests, Miner. Metall. Process, 1989, Vol. 6, pp. 60-66. [7] Xu, Tan and Zhou, Impact Wear Mechanisms of Medium Carbon Steel Under Various Dry and Wet Conditions, The ASME J. of Wear of Materials, 1991, pp. 33-40. [8] Scieszka and Dutkiewicz, Testing Abrasive Wear in Mineral Comminution, The International J. of Mineral Processing, 1991, Vol. 32, pp. 81-110. lEcole [9] Dorlot, Ba lon and Masounave, Des Mat eriaux, Ed. Polytechnique, Montr eal, Quebec, 1986. [10] P. Radziszewski, Exploring Total Media Wear, JKMRC/AMIRA P9L Technical Project Report, Queensland, Australia, 1999. [11] P. Radziszewski, An Update on Mill Media Impact Test Improvement, JKMRC/AMIRA P9L Technical Project Report, Queensland, Australia, 2000.
10