Tumbling Mill (Media Trajectory): Difference between revisions

From Met Dynamics
Jump to navigation Jump to search
m (1 revision imported)
imported>Scott.Munro
 
(3 intermediate revisions by 3 users not shown)
Line 14: Line 14:
== Model theory ==
== Model theory ==


{{Restricted content}}
<hide>
[[File:TumblingMillMediaTraj2.png|thumb|625px|Figure 2. Geometric diagram of an outer media element resting on a lifter bar in a tumbling mill at the point of equilibrium (after Powell, 1991).{{Powell (1991)}}]]
[[File:TumblingMillMediaTraj2.png|thumb|625px|Figure 2. Geometric diagram of an outer media element resting on a lifter bar in a tumbling mill at the point of equilibrium (after Powell, 1991).{{Powell (1991)}}]]


Line 48: Line 51:


where:
where:
* <math>\kappa</math> is the friction angle, related to the coefficient of static friction between the lifter bar and the media element, <math>\mu_s</math>, by <math>\mu_s = \tan \kappa</math>. If the effects of static friction are ignored, <math>\mu_s = 0 \implies \kappa = 0</math>.
* <math>\kappa</math> is the friction angle, related to the coefficient of static friction between the lifter bar and the media element, <math>\mu_{\rm s}</math>, by <math>\mu_{\rm s} = \tan \kappa</math>. If the effects of static friction are ignored, <math>\mu_{\rm s} = 0 \implies \kappa = 0</math>.
* <math>\Omega</math> is the angular velocity of the mill (rad/s)
* <math>\Omega</math> is the angular velocity of the mill (rad/s)
* <math>g</math> is acceleration due to gravity (m/s<sup>2</sup>)
* <math>g</math> is acceleration due to gravity (m/s<sup>2</sup>)
Line 63: Line 66:
where the distance <math>\delta</math> (m) is:
where the distance <math>\delta</math> (m) is:


:<math>\delta = (R - h_r) \cdot \sin \alpha - a</math>
:<math>\delta = (R - h_{\rm r}) \cdot \sin \alpha - a</math>


and <math>h_r</math> is the radial height of the lifter bar (m) and <math>\alpha</math> is the angle between the lifter bar face and a radial line from the centre of the mill through the tip of the lifter bar (rad).
and <math>h_{\rm r}</math> is the radial height of the lifter bar (m) and <math>\alpha</math> is the angle between the lifter bar face and a radial line from the centre of the mill through the tip of the lifter bar (rad).


The angle <math>\alpha</math> is related to the lifter bar geometry by the following equations:{{Powell (1988)}}
The angle <math>\alpha</math> is related to the lifter bar geometry by the following equations:{{Powell (1988)}}
Line 84: Line 87:
* <math>y</math> is half the width of the lifter bar (m).
* <math>y</math> is half the width of the lifter bar (m).


The radial height of the lifter bar, <math>h_r</math>, is the length of the radial line between the tip of the lifter and the base, and is defined as:
The radial height of the lifter bar, <math>h_{\rm r}</math>, is the length of the radial line between the tip of the lifter and the base, and is defined as:


:<math>h_r = R - \frac{\sqrt{R^2 - y^2} - h}{\cos \psi}</math>
:<math>h_{\rm r} = R - \frac{\sqrt{R^2 - y^2} - h}{\cos \psi}</math>


The angle <math>\lambda</math> will always be <math>\frac{\pi}{2}</math> (i.e. 90°) when the lifter bar surface is tangent to the media element, i.e. when the lifter height exceeds the media element radius, as illustrated in Figure 2. However, if the lifter height is lower than the media element radius, <math>\lambda</math> is greater than 90°, as illustrated in Figure 4, and is computed as:
The angle <math>\lambda</math> will always be <math>\frac{\pi}{2}</math> (i.e. 90°) when the lifter bar surface is tangent to the media element, i.e. when the lifter height exceeds the media element radius, as illustrated in Figure 2. However, if the lifter height is lower than the media element radius, <math>\lambda</math> is greater than 90°, as illustrated in Figure 4, and is computed as:


:<math>\lambda = \pi - \arccos \left ( \dfrac{a^2 + (R - h_r)^2 - r_0^2}{2 a (R - h_r)} \right ) - \alpha</math>
:<math>\lambda = \pi - \arccos \left ( \dfrac{a^2 + (R - h_{\rm r})^2 - r_0^2}{2 a (R - h_{\rm r})} \right ) - \alpha</math>


A different formulation of <math>\gamma_0</math> at the point of equilibrium arises from the force balance when <math>\lambda</math> is greater than 90°.
A different formulation of <math>\gamma_0</math> at the point of equilibrium arises from the force balance when <math>\lambda</math> is greater than 90°.
Line 106: Line 109:
The position of a media element at tip of the lifter bar is defined by the following parameters:
The position of a media element at tip of the lifter bar is defined by the following parameters:


:<math>s_L = (R - h_r) \cdot \cos \alpha</math>
:<math>s_{\rm L} = (R - h_{\rm r}) \cdot \cos \alpha</math>


:<math>\beta_L = \arctan \left ( \frac{\delta}{s_L} \right )</math>
:<math>\beta_{\rm L} = \arctan \left ( \frac{\delta}{s_{\rm L}} \right )</math>


where <math>s_L</math> is the position of the media element along the lifter face surface at the tip of the lifter (m), and <math>\beta_L</math> (rad) is the angle between the radial line passing through the centre of the media element at the tip of the lifter and the lifter bar face (rad).
where <math>s_{\rm L}</math> is the position of the media element along the lifter face surface at the tip of the lifter (m), and <math>\beta_{\rm L}</math> (rad) is the angle between the radial line passing through the centre of the media element at the tip of the lifter and the lifter bar face (rad).


In between the point of equilibrium and the lifter bar tip, a media element may roll, slide or be projected into free flight depending on the prevailing forces.
In between the point of equilibrium and the lifter bar tip, a media element may roll, slide or be projected into free flight depending on the prevailing forces.
Line 116: Line 119:
=== Rolling ===
=== Rolling ===


If the coefficient of static friction between the lifter and the media element, <math>\mu_s</math>, is greater than zero then the element commences movement from the point of equilibrium by pure rolling down the lifter face.
If the coefficient of static friction between the lifter and the media element, <math>\mu_{\rm s}</math>, is greater than zero then the element commences movement from the point of equilibrium by pure rolling down the lifter face.


However, if the effects of static friction are ignored, rolling does not occur and the element only experiences sliding motion (see [[Tumbling Mill (Media Trajectory)#Sliding|Sliding]] below).
However, if the effects of static friction are ignored, rolling does not occur and the element only experiences sliding motion (see [[Tumbling Mill (Media Trajectory)#Sliding|Sliding]] below).
Line 152: Line 155:
</math>
</math>


The time at which a media element undergoing pure rolling will reach the tip of the lifter, <math>t_L</math> (s), may be found by setting <math>s(t) = s_L</math> and solving for <math>t</math>. This will require a numerical solution, as <math>t</math> cannot be algebraically isolated in the equation for <math>s(t)</math> above.  
The time at which a media element undergoing pure rolling will reach the tip of the lifter, <math>t_{\rm L}</math> (s), may be found by setting <math>s(t) = s_{\rm L}</math> and solving for <math>t</math>. This will require a numerical solution, as <math>t</math> cannot be algebraically isolated in the equation for <math>s(t)</math> above.  


The velocity of the media element may be found by substituting the value of <math>t_L</math> obtained into the equation for <math>\dot s</math> under pure rolling. The angle of the lifter bar face when the media element is at the tip of the lifter is similarly <math>\gamma_L = \gamma_0 + \Omega t_L</math> under pure rolling.
The velocity of the media element may be found by substituting the value of <math>t_{\rm L}</math> obtained into the equation for <math>\dot s</math> under pure rolling. The angle of the lifter bar face when the media element is at the tip of the lifter is similarly <math>\gamma_{\rm L} = \gamma_0 + \Omega t_{\rm L}</math> under pure rolling.


==== Maximum angle for pure rolling ====
==== Maximum angle for pure rolling ====
Line 160: Line 163:
During pure rolling and prior to reaching the lifter tip, the frictional force causing angular acceleration of the element may be exceeded by the acting gravitational and centrifugal forces. At this point, the element begins to slip across the lifter face, in addition to continued rolling.
During pure rolling and prior to reaching the lifter tip, the frictional force causing angular acceleration of the element may be exceeded by the acting gravitational and centrifugal forces. At this point, the element begins to slip across the lifter face, in addition to continued rolling.


A potential transition from pure rolling to combined rolling and sliding between the point of equilibrium and the tip of the lifter bar can be detected  by finding the minimum possible value of the maximum angle for pure rolling, <math>\gamma_{m\text{(min)}}</math> (rad), i.e. where the following inequality no longer holds:
A potential transition from pure rolling to combined rolling and sliding between the point of equilibrium and the tip of the lifter bar can be detected  by finding the minimum possible value of the maximum angle for pure rolling, <math>\gamma_{\rm m(min)}</math> (rad), i.e. where the following inequality no longer holds:


:<math>\sin \gamma_{m\text{(min)}} - c_4 \mu_s \cdot \cos \gamma_{m\text{(min)}} \leq \dfrac{\Omega^2}{g} \left ( c_4 \delta \mu_s + s_L \right)</math>
:<math>\sin \gamma_{\rm m(min)} - c_4 \mu_{\rm s} \cdot \cos \gamma_{\rm m(min)} \leq \dfrac{\Omega^2}{g} \left ( c_4 \delta \mu_{\rm s} + s_{\rm L} \right)</math>


where:
where:
Line 174: Line 177:
</math>
</math>


The variable <math>\gamma_{m\text{(min)}}</math> must be solved by numerical approximation.
The variable <math>\gamma_{\rm m(min)}</math> must be solved by numerical approximation.


If the angle of the lifter bar face when the media element is at the tip of the lifter after pure rolling , <math>\gamma_L</math>, exceeds <math>\gamma_{m\text{(min)}}</math> then a transition to combined rolling and sliding is possible. Otherwise, the media element only undergoes pure rolling before being projected into free flight.
If the angle of the lifter bar face when the media element is at the tip of the lifter after pure rolling , <math>\gamma_{\rm L}</math>, exceeds <math>\gamma_{\rm m(min)}</math> then a transition to combined rolling and sliding is possible. Otherwise, the media element only undergoes pure rolling before being projected into free flight.


==== Transition to rolling and sliding ====
==== Transition to rolling and sliding ====


The time at which a transition from pure rolling to combined rolling and sliding occurs, <math>t_I</math> (s), is found by numerically solving for <math>t</math> where the following inequality no longer holds:
The time at which a transition from pure rolling to combined rolling and sliding occurs, <math>t_{\rm I}</math> (s), is found by numerically solving for <math>t</math> where the following inequality no longer holds:


:<math>c_5 . \sin (\gamma_0 + \Omega t ) - c_4 \mu_s \cdot \cos (\gamma_0 + \Omega t ) \leq  
:<math>c_5 . \sin (\gamma_0 + \Omega t ) - c_4 \mu_{\rm s} \cdot \cos (\gamma_0 + \Omega t ) \leq  


c_4.\dfrac{\Omega^2}{g} \delta \mu_s + \left ( \dfrac{s_0}{g}\Omega^2 - c_1 \sin \gamma_0 \right ) \cdot \cosh \left ( c_2 . \Omega t \right ) - c_6 \sqrt{3} \cos \gamma_0 \cdot \sinh \left ( c_2 . \Omega t \right )</math>
c_4.\dfrac{\Omega^2}{g} \delta \mu_{\rm s} + \left ( \dfrac{s_0}{g}\Omega^2 - c_1 \sin \gamma_0 \right ) \cdot \cosh \left ( c_2 . \Omega t \right ) - c_6 \sqrt{3} \cos \gamma_0 \cdot \sinh \left ( c_2 . \Omega t \right )</math>


where:  
where:  
Line 201: Line 204:
</math>
</math>


The position, <math>s_I</math> (m), and velocity, <math>\dot s_I</math> (m), of the media element at <math>t_I</math> can then be computed using the pure rolling equations for <math>s(t)</math> and <math>\dot s</math> above. The angle of the lifter bar face at <math>t_I</math>, <math>\gamma_I</math> (rad), is found from <math>\gamma_I = \gamma_0 + \Omega t_I</math>.
The position, <math>s_{\rm I}</math> (m), and velocity, <math>\dot s_{\rm I}</math> (m), of the media element at <math>t_{\rm I}</math> can then be computed using the pure rolling equations for <math>s(t)</math> and <math>\dot s</math> above. The angle of the lifter bar face at <math>t_{\rm I}</math>, <math>\gamma_{\rm I}</math> (rad), is found from <math>\gamma_{\rm I} = \gamma_0 + \Omega t_{\rm I}</math>.


=== Sliding ===
=== Sliding ===
Line 209: Line 212:
If sliding motion commences after rolling, the position and velocity of the element at time <math>\tau</math> are:
If sliding motion commences after rolling, the position and velocity of the element at time <math>\tau</math> are:


:<math>s(\tau) = \left( s_I + \mu_k \delta + \dfrac{g}{2 \Omega^2} (\mu_k \cos \gamma_I - \sin \gamma_I) \right ) \cosh (\Omega \tau) + \left ( \dfrac{\dot s_I}{\Omega} - \dfrac{g}{2 \Omega^2} (\mu_k \sin \gamma_I + \cos \gamma_I ) \right ) \sinh (\Omega \tau) - \dfrac{g}{2 \Omega^2} \big( \mu_k \cos (\gamma_I + \Omega \tau) - \sin (\gamma_I + \Omega \tau) \big) - \mu_k \delta</math>
:<math>s(\tau) = \left( s_{\rm I} + \mu_{\rm k} \delta + \dfrac{g}{2 \Omega^2} (\mu_{\rm k} \cos \gamma_{\rm I} - \sin \gamma_{\rm I}) \right ) \cosh (\Omega \tau) + \left ( \dfrac{\dot s_{\rm I}}{\Omega} - \dfrac{g}{2 \Omega^2} (\mu_{\rm k} \sin \gamma_{\rm I} + \cos \gamma_{\rm I} ) \right ) \sinh (\Omega \tau) - \dfrac{g}{2 \Omega^2} \big( \mu_{\rm k} \cos (\gamma_{\rm I} + \Omega \tau) - \sin (\gamma_{\rm I} + \Omega \tau) \big) - \mu_{\rm k} \delta</math>


:<math>\dot s = \Omega \left ( s_I + \mu_k \delta + \dfrac{g}{2 \Omega^2} (\mu_k \cos \gamma_I - \sin \gamma_I \right ) \sinh (\Omega \tau) + \left ( \dot s_I -\dfrac{g}{2 \Omega} (\mu_K \sin \gamma_I + \cos \gamma_I \right ) \cosh (\Omega \tau) - \dfrac{g}{2 \Omega} \big( -\mu_k \sin (\gamma_I + \Omega \tau) - \cos (\gamma_I + \Omega \tau) \big )</math>
:<math>\dot s = \Omega \left ( s_{\rm I} + \mu_{\rm k} \delta + \dfrac{g}{2 \Omega^2} (\mu_{\rm k} \cos \gamma_{\rm I} - \sin \gamma_{\rm I} \right ) \sinh (\Omega \tau) + \left ( \dot s_{\rm I} -\dfrac{g}{2 \Omega} (\mu_K \sin \gamma_{\rm I} + \cos \gamma_{\rm I} \right ) \cosh (\Omega \tau) - \dfrac{g}{2 \Omega} \big( -\mu_{\rm k} \sin (\gamma_{\rm I} + \Omega \tau) - \cos (\gamma_{\rm I} + \Omega \tau) \big )</math>


where <math>\mu_k</math> is the coefficient of kinematic friction between the lifter bar face and the media element.
where <math>\mu_{\rm k}</math> is the coefficient of kinematic friction between the lifter bar face and the media element.


The time at which the media element reaches the tip of the lifter, <math>\tau_L</math> (s), after the transition to sliding is found by setting <math>s(\tau)=s_L</math> and numerically solving for <math>\tau</math>. The velocity <math>\dot s_L</math> (m/s) and angular position <math>\gamma_L</math> (rad) of the media element at the tip of the lifter after sliding then follow from the value of <math>\tau_L</math> obtained (i.e., <math>\gamma_L = \gamma_I + \Omega \tau</math>).
The time at which the media element reaches the tip of the lifter, <math>\tau_{\rm L}</math> (s), after the transition to sliding is found by setting <math>s(\tau)=s_{\rm L}</math> and numerically solving for <math>\tau</math>. The velocity <math>\dot s_{\rm L}</math> (m/s) and angular position <math>\gamma_{\rm L}</math> (rad) of the media element at the tip of the lifter after sliding then follow from the value of <math>\tau_{\rm L}</math> obtained (i.e., <math>\gamma_{\rm L} = \gamma_{\rm I} + \Omega \tau</math>).


The total time for the media element to travel from the point of equilibrium to the tip of the lifter bar by rolling and sliding is <math>t_L = t_I + \tau_L</math> (s).
The total time for the media element to travel from the point of equilibrium to the tip of the lifter bar by rolling and sliding is <math>t_{\rm L} = t_{\rm I} + \tau_{\rm L}</math> (s).


==== Pure sliding ====
==== Pure sliding ====


If the coefficient of static friction (<math>\mu_s</math>) is zero, i.e. ignored, rolling does not commence and the media element undergoes sliding directly from the point of equilibrium. The above equations of motion for sliding are subsequently utilised with the following values substituted:
If the coefficient of static friction (<math>\mu_{\rm s}</math>) is zero, i.e. ignored, rolling does not commence and the media element undergoes sliding directly from the point of equilibrium. The above equations of motion for sliding are subsequently utilised with the following values substituted:


:<math>s_I = s_0, \;\;\;\;\;
:<math>s_{\rm I} = s_0, \;\;\;\;\;
\dot s_I = 0, \;\;\;\;\;
\dot s_{\rm I} = 0, \;\;\;\;\;
\gamma_I = \gamma_0</math>
\gamma_{\rm I} = \gamma_0</math>


and <math>t_L = \tau_L</math>.
and <math>t_{\rm L} = \tau_{\rm L}</math>.


=== Free-flight ===
=== Free-flight ===
Line 245: Line 248:
The polar coordinates of the media element position at the point it leaves the lifter are (m, rad):
The polar coordinates of the media element position at the point it leaves the lifter are (m, rad):


:<math>(r_L,\, \phi_L) =\left ( \dfrac{s_L}{\cos \beta_L},\; \gamma_L - \beta_L \right )</math>
:<math>(r_{\rm L},\, \phi_{\rm L}) =\left ( \dfrac{s_{\rm L}}{\cos \beta_{\rm L}},\; \gamma_{\rm L} - \beta_{\rm L} \right )</math>


Switching to a Cartesian coordinate frame of reference yields the following position <math>(x_L,y_L)</math> (m) and velocity <math>(v_{xL},v_{yL})</math> (m/s) coordinates at the tip of the lifter bar:
Switching to a Cartesian coordinate frame of reference yields the following position <math>(x_{\rm L},y_{\rm L})</math> (m) and velocity <math>(v_{x{\rm L}},v_{y{\rm L}})</math> (m/s) coordinates at the tip of the lifter bar:


:<math>(x_L,\, y_L) = \left ( r_L \cos \phi_L,\; r_L \sin \phi_L \right )</math>
:<math>(x_{\rm L},\, y_{\rm L}) = \left ( r_{\rm L} \cos \phi_{\rm L},\; r_{\rm L} \sin \phi_{\rm L} \right )</math>


:<math>(v_{xL},\, v_{yL}) = \left ( \dot s_L \cos \gamma_L - \Omega r_L \sin \phi_L,\; \dot s _L \sin \gamma_L + \Omega r_L \cos \phi_L \right )</math>
:<math>(v_{x{\rm L}},\, v_{y{\rm L}}) = \left ( \dot s_{\rm L} \cos \gamma_{\rm L} - \Omega r_{\rm L} \sin \phi_{\rm L},\; \dot s _L \sin \gamma_{\rm L} + \Omega r_{\rm L} \cos \phi_{\rm L} \right )</math>


The polar coordinates of the media element velocity as it leaves the lifter are (m/s, rad):
The polar coordinates of the media element velocity as it leaves the lifter are (m/s, rad):


:<math>(v_L,\, \phi_{vL}) =\left (\sqrt{v^2_{xL} + v^2_{yL}},\; \arctan \dfrac{v_{xL}}{v_{yL}} \right )</math>
:<math>(v_{\rm L},\, \phi_{v{\rm L}}) =\left (\sqrt{{{v_{x{\rm L}}^2}} + v^2_{y{\rm L}}},\; \arctan \dfrac{v_{x{\rm L}}}{v_{y{\rm L}}} \right )</math>


Once in free flight, the media element moves in a parabolic path under the influence of gravity, with a position defined as (m):
Once in free flight, the media element moves in a parabolic path under the influence of gravity, with a position defined as (m):


:<math>(x,\, y) = \left ( x_L + v_{xL} t,\; y_L + v_{yL} t - \frac{1}{2} g t^2 \right )</math>
:<math>(x,\, y) = \left ( x_{\rm L} + v_{x{\rm L}} t,\; y_{\rm L} + v_{y{\rm L}} t - \frac{1}{2} g t^2 \right )</math>


The media element impacts the mill shell when:
The media element impacts the mill shell when:
Line 265: Line 268:
:<math>x^2 + y^2 = r_0</math>
:<math>x^2 + y^2 = r_0</math>


the solution of which is denoted <math>(x_E, y_E)</math> (m).
the solution of which is denoted <math>(x_{\rm E}, y_{\rm E})</math> (m).


Finally, the polar coordinates of the media element velocity as it leaves the lifter are (m/s, rad):
Finally, the polar coordinates of the media element velocity as it leaves the lifter are (m/s, rad):


:<math>(v_E,\, \phi_{vE}) =\left (\sqrt{v^2_{xE} + v^2_{yE}},\; \arctan \dfrac{v_{xE}}{v_{yE}} \right )</math>
:<math>(v_{\rm E},\, \phi_{v{\rm E}}) =\left (\sqrt{v^2_{x{\rm E}} + v^2_{y{\rm E}}},\; \arctan \dfrac{v_{x{\rm E}}}{v_{y{\rm E}}} \right )</math>


where:
where:


:<math>(v_{xE},\, v_{yE}) = (v_{xL}, \; v_{yL} - gt)</math>
:<math>(v_{x{\rm E}},\, v_{y{\rm E}}) = (v_{x{\rm L}}, \; v_{y{\rm L}} - gt)</math>


== Additional notes ==
=== Additional notes ===


{| style="float:right; margin-left: 20px;"
{| style="float:right; margin-left: 20px;"
Line 282: Line 285:
|}
|}


=== Friction estimates ===
==== Friction estimates ====


Powell's (1991) experimental work separately investigated the coefficients of static  and kinetic friction between a media element and a lifter bar. The findings suggest:
Powell's (1991) experimental work separately investigated the coefficients of static  and kinetic friction between a media element and a lifter bar. The findings suggest:
* Static friction (<math>\mu_s</math>) is likely approaching zero for practical applications, and
* Static friction (<math>\mu_{\rm s}</math>) is likely approaching zero for practical applications, and
* Kinematic friction (<math>\mu_k</math>) is likely in the range of 0.17 - 0.19.
* Kinematic friction (<math>\mu_{\rm k}</math>) is likely in the range of 0.17 - 0.19.


Rolling motion does not occur when static friction is zero, greatly simplifying the trajectory calculations presented above.
Rolling motion does not occur when static friction is zero, greatly simplifying the trajectory calculations presented above.


Adequate estimates of the coefficients of friction for simulation purposes may be 0 or 0.05 for <math>\mu_s</math> and 0.19 for <math>\mu_k</math>, but is left to user discretion.
Adequate estimates of the coefficients of friction for simulation purposes may be 0 or 0.05 for <math>\mu_{\rm s}</math> and 0.19 for <math>\mu_{\rm k}</math>, but is left to user discretion.


=== Extracting lifter geometry ===
==== Extracting lifter geometry ====


Modern mill interior laser scanning devices and/or design drawings provide accurate parametric maps of lifter bar geometry, which may be used to quantify the lifter height (<math>h</math>), lifter width (<math>2y</math>) and lifter face angle (<math>\rho</math>) properties for media trajectory simulation.  
Modern mill interior laser scanning devices and/or design drawings provide accurate parametric maps of lifter bar geometry, which may be used to quantify the lifter height (<math>h</math>), lifter width (<math>2y</math>) and lifter face angle (<math>\rho</math>) properties for media trajectory simulation.  
Line 308: Line 311:
<br>
<br>
<br>
<br>
</hide>


== Excel ==
== Excel ==
Line 327: Line 331:
\text{Media type}\\
\text{Media type}\\
d\text{ (m)}\\
d\text{ (m)}\\
\mu_s\\
\mu_{\rm s}\\
\mu_k\\
\mu_{\rm k}\\
\rho\text{  (rad)}\\
\rho\text{  (rad)}\\
h\text{  (m)}\\
h\text{  (m)}\\
y\text{  (m)}\\
y\text{  (m)}\\
C_s (frac)
C_{\rm s} (frac)
\end{bmatrix},\;\;\;\;\;\;
\end{bmatrix},\;\;\;\;\;\;


Line 344: Line 348:
\end{bmatrix}\\
\end{bmatrix}\\
\begin{bmatrix}
\begin{bmatrix}
\phi_L\text{ (rad)} & -\\
\phi_{\rm L}\text{ (rad)} & -\\
\end{bmatrix}\\
\end{bmatrix}\\
\begin{bmatrix}
\begin{bmatrix}
x_L \text{ (m)} & y_L \text{ (m)}\\
x_{\rm L} \text{ (m)} & y_{\rm L} \text{ (m)}\\
v_L \text{ (m)} & \phi_{vL} \text{ (rad)}\\
v_{\rm L} \text{ (m)} & \phi_{v{\rm L}} \text{ (rad)}\\
\end{bmatrix}\\
\end{bmatrix}\\
\begin{bmatrix}
\begin{bmatrix}
\phi_E\text{ (rad)} & -\\
\phi_{\rm E}\text{ (rad)} & -\\
\end{bmatrix}\\
\end{bmatrix}\\
\begin{bmatrix}
\begin{bmatrix}
x_E \text{ (m)} & y_E \text{ (m)}\\
x_{\rm E} \text{ (m)} & y_{\rm E} \text{ (m)}\\
v_E \text{ (m)} & \phi_{vE} \text{ (rad)}\\
v_{\rm E} \text{ (m)} & \phi_{v{\rm E}} \text{ (rad)}\\
\end{bmatrix}\\
\end{bmatrix}\\
\begin{bmatrix}
\begin{bmatrix}
Line 370: Line 374:
* <math>d</math> is the media element diameter (m), i.e. <math>d = 2a</math>
* <math>d</math> is the media element diameter (m), i.e. <math>d = 2a</math>
* <math>\text{Media type}</math> is the type of media element, balls or rods (''0 = Balls, 1 = Rods'')
* <math>\text{Media type}</math> is the type of media element, balls or rods (''0 = Balls, 1 = Rods'')
* <math>C_s</math> is the fraction critical speed of the mill (frac)
* <math>C_{\rm s}</math> is the fraction critical speed of the mill (frac)


and:
and:
Line 431: Line 435:
! colspan="3" style="text-align:left;" |''Results''
! colspan="3" style="text-align:left;" |''Results''
|-
|-
| colspan="3" |''Point Of Equilibrium''
| colspan="3" |''PointOfEquilibrium''
|-
|-
| PositionAngle
| PositionAngle / Phi
|style="background: #eaecf0" | Display
|style="background: #eaecf0" | Display
| Angular position of the point of equilibrium.
| Angular position of the point of equilibrium.
|-
|-
| colspan="3" |''Leaves Lifter''
| colspan="3" |''LeavesLifter''
|-
|-
|PositionAngle
|PositionAngle / Phi
|style="background: #eaecf0" | Display
|style="background: #eaecf0" | Display
| Angular position of the point where the media element leaves the lifter.
| Angular position of the point where the media element leaves the lifter.
Line 451: Line 455:
| Cartesian ''y'' coordinate of the point where the media element leaves the lifter.
| Cartesian ''y'' coordinate of the point where the media element leaves the lifter.
|-
|-
|Velocity
|Velocity / V
|style="background: #eaecf0" | Display
|style="background: #eaecf0" | Display
| Velocity of the media element at the point where it leaves the lifter.
| Velocity of the media element at the point where it leaves the lifter.
|-
|-
|VelocityAngle
|VelocityAngle / PhiV
|style="background: #eaecf0" | Display
|style="background: #eaecf0" | Display
| Angular direction of the velocity of the media element at the point where it leaves the lifter.
| Angular direction of the velocity of the media element at the point where it leaves the lifter.
|-
|-
| colspan="3" |''Impacts Shell''
| colspan="3" |''ImpactsShell''
|-
|-
|PositionAngle
|PositionAngle / Phi
|style="background: #eaecf0" | Display
|style="background: #eaecf0" | Display
| Angular position of the point where the media element impacts the mill shell.
| Angular position of the point where the media element impacts the mill shell.
Line 473: Line 477:
| Cartesian ''y'' coordinate of the point where the media element impacts the mill shell.
| Cartesian ''y'' coordinate of the point where the media element impacts the mill shell.
|-
|-
|Velocity
|Velocity / V
|style="background: #eaecf0" | Display
|style="background: #eaecf0" | Display
| Velocity of the media element at the point where it impacts the mill shell.
| Velocity of the media element at the point where it impacts the mill shell.
|-
|-
|VelocityAngle
|VelocityAngle / PhiV
|style="background: #eaecf0" | Display
|style="background: #eaecf0" | Display
| Angular direction of the velocity of the media element at the point where it impacts the mill shell.
| Angular direction of the velocity of the media element at the point where it impacts the mill shell.
|-
|-
| style="vertical-align:top;" |State
| style="vertical-align:top;" |State
|style="background: #eaecf0" style="vertical-align:top;"| Display
|style="vertical-align:top; background:#eaecf0" | Display
| Text description of the motion of the media element. Possible values are:
| Text description of the motion of the media element. Possible values are:
* "Element is projected into free flight from point of equilibrium."
* "Element is projected into free flight from point of equilibrium."
Line 493: Line 497:
|-
|-
! colspan="3" style="text-align:left;" |''Trajectory''
! colspan="3" style="text-align:left;" |''Trajectory''
|-
| ShowTrajectory
| CheckBox
| If enabled, the (x, y) coordinate data of the element trajectory is displayed below.
|-
|-
|TrajectoryData
|TrajectoryData
Line 500: Line 508:
|Point, X, Y
|Point, X, Y
|style="background: #eaecf0" | Display
|style="background: #eaecf0" | Display
| X and Y coordinates of the media element trajectory, from the point of equilibrium to point of impact with the mill shell.
| X and Y coordinates of the media element trajectory, from the point of equilibrium to point of impact with the mill shell. Only visible if ''ShowTrajectory'' is enabled above.
|}
|}


Latest revision as of 11:19, 4 December 2024

Description

Figure 1. Example tumbling mill media trajectory, including Morrell Continuum charge position.

This article describes Powell's (1991) method for predicting the trajectory of the outer media elements in a tumbling mill.[1]

The approach predicts the point of impact of media elements projected into free flight after transit down a lifter bar surface, based on element size, lifter bar geometry, mill diameter and mill rotational speed.

The tumbling mill media trajectory model is particularly useful when combined with an estimate of charge position (e.g. Morrell's method), which allows analyses such as:

  • Identifying and mitigating the risks of shell liner damage from media impact through lifter bar (or mill) design
  • Monitoring trajectory and point of impact changes with charge level, mill speed or liner and lifter bar wear in operational mills
  • Limiting the maximum operating speed (and hence power and throughput) in a mill simulation or process control sub-systems

Model theory

Nuvola apps important blue.svg.png This content is available to registered users. Please log in to view.

Excel

The tumbling mill media trajectory model may be invoked from the Excel formula bar with the following function call: 

=mdMillPower_MillMediaTrajectory(Parameters as Range)

Invoking the function with no arguments will print Help text associated with the model, including a link to this page.

The Parameters array and model results are defined below in matrix notation, along with an example image showing the selection of the same arrays in the Excel interface:


where:

  • is the mill diameter inside liners (m), i.e.
  • is the media element diameter (m), i.e.
  • is the type of media element, balls or rods (0 = Balls, 1 = Rods)
  • is the fraction critical speed of the mill (frac)

and:

  • is a text string that describes the motion of the media element
  • is the number of trajectory points for plotting
Figure 7. Example showing the selection of the Parameters (blue frame), and Results (light blue frame) arrays in Excel.

SysCAD

The media trajectory model is an optional calculation for tumbling mill units. If selected, the input and display parameters below are shown.

Tag (Long/Short) Input / Display Description/Calculated Variables/Options
MediaTraj
HelpLink ButtonModelHelp.png Opens a link to this page using the system default web browser. Note: Internet access is required.
Requirements
MillDiameter Input/Display Diameter of the mill (inside liners).
FracCS Input/Display Fraction critical speed of the mill.
MediaType Ball/Rod Type of media element, balls or rods.
MediaSize Display Diameter of the media element.
mus Input Coefficient of static friction.
muk Input Coefficient of kinematic friction.
rho Input Lifter face angle, angle between the face surface and the base of the lifter bar.
h Input Height of the lifter bar.
2y Input Width of the lifter bar.
Results
PointOfEquilibrium
PositionAngle / Phi Display Angular position of the point of equilibrium.
LeavesLifter
PositionAngle / Phi Display Angular position of the point where the media element leaves the lifter.
X Display Cartesian x coordinate of the point where the media element leaves the lifter.
Y Display Cartesian y coordinate of the point where the media element leaves the lifter.
Velocity / V Display Velocity of the media element at the point where it leaves the lifter.
VelocityAngle / PhiV Display Angular direction of the velocity of the media element at the point where it leaves the lifter.
ImpactsShell
PositionAngle / Phi Display Angular position of the point where the media element impacts the mill shell.
X Display Cartesian x coordinate of the point where the media element impacts the mill shell.
Y Display Cartesian y coordinate of the point where the media element impacts the mill shell.
Velocity / V Display Velocity of the media element at the point where it impacts the mill shell.
VelocityAngle / PhiV Display Angular direction of the velocity of the media element at the point where it impacts the mill shell.
State Display Text description of the motion of the media element. Possible values are:
  • "Element is projected into free flight from point of equilibrium."
  • "Element slides to end of the lifter and is then projected into free flight."
  • "Element slides and is projected into free flight before the end of the lifter."
  • "Element rolls to the end of the lifter and is projected into free flight."
  • "Element rolls then slides to the end of the lifter before being projected into free flight."
  • "Element rolls then slides and is projected into free flight before the end of the lifter."
  • "Element rolls and is projected into free flight before the end of the lifter."
Trajectory
ShowTrajectory CheckBox If enabled, the (x, y) coordinate data of the element trajectory is displayed below.
TrajectoryData ButtonCopyToClipboard.png Copies the full set of trajectory coordinates to the Windows clipboard in CSV format.
Point, X, Y Display X and Y coordinates of the media element trajectory, from the point of equilibrium to point of impact with the mill shell. Only visible if ShowTrajectory is enabled above.

See also

References

  1. ↑ Powell, M.S., 1991. The effect of liner design on the motion of the outer grinding elements in a rotary mill. International Journal of Mineral Processing, 31(3-4), pp.163-193.
Retrieved from "https://wiki.metdynamics.com.au/index.php?title=Tumbling_Mill_(Media_Trajectory)&oldid=2212"