AG/SAG Mill (Variable Rates): Difference between revisions
imported>Scott.Munro m (→See also) |
imported>Scott.Munro |
||
| (One intermediate revision by the same user not shown) | |||
| Line 56: | Line 56: | ||
* <math>D_{\rm L}</math> is the discharge rate of liquids from the mill, normally assumed to equal <math>D_n</math> (h<sup>-1</sup>). | * <math>D_{\rm L}</math> is the discharge rate of liquids from the mill, normally assumed to equal <math>D_n</math> (h<sup>-1</sup>). | ||
</hide><div class="user-show"> | |||
=== Calculation sequence === | === Calculation sequence === | ||
</div><hide> | |||
The Variable Rates AG/SAG model uses a range of sub-models to quantify the breakage rate (<math>R</math>), appearance function (<math>A</math>), and discharge function (<math>D</math>) terms of the perfect mixing model, and hence compute the mill load (<math>s</math>) and product (<math>p</math>) at steady-state. | The Variable Rates AG/SAG model uses a range of sub-models to quantify the breakage rate (<math>R</math>), appearance function (<math>A</math>), and discharge function (<math>D</math>) terms of the perfect mixing model, and hence compute the mill load (<math>s</math>) and product (<math>p</math>) at steady-state. | ||
| Line 64: | Line 66: | ||
The sub-models are described in further detail below. | The sub-models are described in further detail below. | ||
</hide><div class="user-show"> | |||
=== Breakage rates === | === Breakage rates === | ||
</div><hide> | |||
[[File:AGSAGVariableRates12.png|thumb|450px|Figure 4. Breakage rate distribution characterised by cubic spline interpolation.]] | [[File:AGSAGVariableRates12.png|thumb|450px|Figure 4. Breakage rate distribution characterised by cubic spline interpolation.]] | ||
{{Model theory (Text, AG/SAG Mill, Variable Rates, Breakage Rates)|4}} | |||
</hide><div class="user-show"> | |||
=== Discharge rates === | === Discharge rates === | ||
</div><hide> | |||
The discharge rates (<math> | The discharge rates (<math>D_i</math>) are related to the hold-up of slurry in the mill and particle classification at the discharge grates. | ||
</hide><div class="user-show"> | |||
==== Slurry hold-up ==== | ==== Slurry hold-up ==== | ||
</div><hide> | |||
[[File:AGSAGVariableRates9.png|thumb|450px|Figure 5. Principal dimensions of an AG/SAG mill.]] | [[File:AGSAGVariableRates9.png|thumb|450px|Figure 5. Principal dimensions of an AG/SAG mill.]] | ||
| Line 182: | Line 131: | ||
* <math>r_{\rm m}</math> is the radius of the mill inside the liners (m) | * <math>r_{\rm m}</math> is the radius of the mill inside the liners (m) | ||
In addition to fine slurry, particles up to | In addition to fine slurry, particles up to the grate aperture size will also discharge from the mill. To estimate total discharge flow rate, <math>Q</math> (m<sup>3</sup>/h), Morrell and Stephenson (1996) suggest the following correction:{{Morrell and Stephenson (1996)}} | ||
:<math>Q = k_{\rm g} Q_{\rm m}</math> | :<math>Q = k_{\rm g} Q_{\rm m}</math> | ||
| Line 218: | Line 167: | ||
Thus, slurry hold-up, <math>L_{\rm V}</math>, can be computed for a given feed/discharge flow rate, grate design and mill. | Thus, slurry hold-up, <math>L_{\rm V}</math>, can be computed for a given feed/discharge flow rate, grate design and mill. | ||
</hide><div class="user-show"> | |||
==== Classification and discharge ==== | ==== Classification and discharge ==== | ||
</div><hide> | |||
[[File:AGSAGVariableRates10.png|thumb|450px|Figure 6. Classification function, <math>C_i</math>, with pebble port open are fraction, <math>f_p</math>, specified.]] | [[File:AGSAGVariableRates10.png|thumb|450px|Figure 6. Classification function, <math>C_i</math>, with pebble port open are fraction, <math>f_p</math>, specified.]] | ||
| Line 224: | Line 175: | ||
[[File:AGSAGVariableRates11.png|thumb|450px|Figure 7. Classification function, <math>C_i</math>, where pebble port open are fraction, <math>f_p</math>, is zero, i.e grates only.]] | [[File:AGSAGVariableRates11.png|thumb|450px|Figure 7. Classification function, <math>C_i</math>, where pebble port open are fraction, <math>f_p</math>, is zero, i.e grates only.]] | ||
{{Model theory (Text, AG/SAG Mill, Variable Rates, Classification)}} | |||
Figure 6 shows an example classification function with pebble ports included, whilst Figure 7 shows the same function with a grate-only mill. | Figure 6 shows an example classification function with pebble ports included, whilst Figure 7 shows the same function with a grate-only mill. | ||
The value of <math>d_{\rm max}</math> is adjusted during the calculation sequence (Figure 3) to ensure the fraction of solids less than <math>x_{\rm g}</math> plus water retained in the mill load computed by the perfect mixing population balance matches the slurry hold-up determined by the [[AG/SAG Mill (Variable Rates)#Slurry_flow|slurry flow]] calculations. | The value of <math>d_{\rm max}</math> is adjusted during the calculation sequence (Figure 3) to ensure the fraction of solids less than <math>x_{\rm g}</math> plus water retained in the mill load computed by the perfect mixing population balance matches the slurry hold-up determined by the [[AG/SAG Mill (Variable Rates)#Slurry_flow|slurry flow]] calculations. | ||
</hide><div class="user-show"> | |||
=== Appearance function === | === Appearance function === | ||
{{Model theory (Text, AG/SAG Mill, Variable Rates, Appearance Function)}} | |||
=== Mill power === | === Mill power === | ||
</div><hide> | |||
The Variable Rates AG/SAG model includes an implementation of the [[Tumbling Mill (Power, Morrell Continuum)|Morrell Continuum]] tumbling mill power model. The predicted mill power draw is not utilised by the Variable Rates model formulation in any manner, and is provided for information only. | The Variable Rates AG/SAG model includes an implementation of the [[Tumbling Mill (Power, Morrell Continuum)|Morrell Continuum]] tumbling mill power model. The predicted mill power draw is not utilised by the Variable Rates model formulation in any manner, and is provided for information only. | ||
</hide><div class="user-show"> | |||
==== Charge properties ==== | ==== Charge properties ==== | ||
</div><hide> | |||
The power draw prediction requires an estimate of <math>J_{\rm t}</math> (v/v), the fraction of mill volume occupied by the charge, which includes coarse ore, balls, slurry, and void spaces. | The power draw prediction requires an estimate of <math>J_{\rm t}</math> (v/v), the fraction of mill volume occupied by the charge, which includes coarse ore, balls, slurry, and void spaces. | ||
| Line 459: | Line 229: | ||
:<math>\rho_{\rm c} = \frac{J_{\rm t} \rho_{\rm S} (1 - \varepsilon + \varepsilon U S) + J_{\rm B}( \rho_{\rm B} - \rho_{\rm S})(1 - \varepsilon) + J_{\rm t} \varepsilon U (1 - S)}{J_{\rm t}}, \quad U \leq 1</math> | :<math>\rho_{\rm c} = \frac{J_{\rm t} \rho_{\rm S} (1 - \varepsilon + \varepsilon U S) + J_{\rm B}( \rho_{\rm B} - \rho_{\rm S})(1 - \varepsilon) + J_{\rm t} \varepsilon U (1 - S)}{J_{\rm t}}, \quad U \leq 1</math> | ||
</hide><div class="user-show"> | |||
==== Power draw ==== | ==== Power draw ==== | ||
</div><hide> | |||
:''Main article'': [[Tumbling Mill (Power, Morrell Continuum)]] | :''Main article'': [[Tumbling Mill (Power, Morrell Continuum)]] | ||
| Line 477: | Line 249: | ||
The complete equations are excluded here for brevity and are available at the article link above. | The complete equations are excluded here for brevity and are available at the article link above. | ||
</hide><div class="user-show"> | |||
=== Internal mesh series === | === Internal mesh series === | ||
</div><hide> | |||
{{Model theory (Text, AG/SAG Mill, Variable Rates, Internal Mesh)}} | |||
</hide><div class="user-show"> | |||
=== Multicomponent modelling === | === Multicomponent modelling === | ||
</div><hide> | |||
{{Model theory (Text, AG/SAG Mill, Variable Rates, Multicomponent)}} | |||
</hide><div class="user-show"> | |||
=== Additional notes === | === Additional notes === | ||
</div><hide> | |||
</hide><div class="user-show"> | |||
==== Breakage rates and mill load ==== | ==== Breakage rates and mill load ==== | ||
</div><hide> | |||
An important, and potentially overlooked, limitation of the Variable Rates AG/SAG mill model is the insensitivity of the breakage rate relationships to mill load. Mill simulations should therefore use mill loads close or equal to the load observed during model fitting, or 25% for design activities.{{Bailey et al. (2009)}} | An important, and potentially overlooked, limitation of the Variable Rates AG/SAG mill model is the insensitivity of the breakage rate relationships to mill load. Mill simulations should therefore use mill loads close or equal to the load observed during model fitting, or 25% for design activities.{{Bailey et al. (2009)}} | ||
</hide><div class="user-show"> | |||
==== Slurry pool ==== | ==== Slurry pool ==== | ||
</div><hide> | |||
Various published descriptions of the Variable Rates AG/SAG mill suggest that slurry pooling phenomena are excluded from slurry hold-up and power draw estimations.{{Morrell and Morrison (1996)}}{{Kojovic et al. (2012)}}{{Bueno et al. (2013)}} | Various published descriptions of the Variable Rates AG/SAG mill suggest that slurry pooling phenomena are excluded from slurry hold-up and power draw estimations.{{Morrell and Morrison (1996)}}{{Kojovic et al. (2012)}}{{Bueno et al. (2013)}} | ||
| Line 553: | Line 307: | ||
x_{\rm m}\text{ (mm)}\\ | x_{\rm m}\text{ (mm)}\\ | ||
\gamma\text{ (m/m)}\\ | \gamma\text{ (m/m)}\\ | ||
k_{\rm m}\text{ (-)}\\ | |||
J_{\rm B}\text{ (v/v)}\\ | J_{\rm B}\text{ (v/v)}\\ | ||
\rho_{\rm B}\text{ (t/m}^3\text{)}\\ | \rho_{\rm B}\text{ (t/m}^3\text{)}\\ | ||
| Line 802: | Line 556: | ||
{{SysCAD (Page, Mill, DLL*Mill)|PowerModels=true|MediaTraj=true|MediaStrings=true}} | {{SysCAD (Page, Mill, DLL*Mill)|PowerModels=true|MediaTraj=true|MediaStrings=true}} | ||
{{SysCAD (Page, AG/SAG Mill, Variable Rates, Mill)}} | |||
{{SysCAD ( | |||
{{SysCAD (Page, AG/SAG Mill, Variable Rates, Ore)}} | |||
{{SysCAD (Page, AG/SAG Mill, Variable Rates, Results)|Steady-state=true}} | |||
{{ | {{SysCAD (Page, AG/SAGMill, Variable Rates, Ri/Di)}} | ||
{{SysCAD ( | {{SysCAD (Page, AG/SAGMill, Variable Rates, Load)}} | ||
{{SysCAD (Page, Tumbling Mill, Power)|modelpage={{SysCAD (Text, UnitType Prefix)}}Mill|HF=true|MorrellC=true|MorrellE=true|MorrellD=true|HildenPowell=true}} | {{SysCAD (Page, Tumbling Mill, Power)|modelpage={{SysCAD (Text, UnitType Prefix)}}Mill|HF=true|MorrellC=true|MorrellE=true|MorrellD=true|HildenPowell=true}} | ||
| Line 1,138: | Line 576: | ||
== See also == | == See also == | ||
* [[ | * [[AG/SAG Mill (Variable Rates, Dynamic)]] | ||
* [ | |||
== External links == | |||
* [https://help.syscad.net/Example_-_09_Met_Dynamics_Projects#Crushing_and_Grinding_Example Crushing and Grinding Example] | |||
* [https://help.syscad.net/Met_Dynamics_-_Mill Met Dynamics - Mill (help.syscad.net)] | |||
* [https://help.syscad.net/Example_-_09_Met_Dynamics_Projects#Crushing_and_Grinding_Example Crushing and Grinding Example project (help.syscad.net)] | |||
* [https://help.syscad.net/Example_-_09_Met_Dynamics_Projects#Crushing,_Grinding_and_Flotation_Example Crushing, Grinding, and Flotation Example project (help.syscad.net)] | |||
== References == | == References == | ||
Latest revision as of 09:20, 9 July 2025
Description
This article describes an implementation of the Autogenous (AG) and Semi-Autogenous (SAG) mill model originated by Leung (1987) and extended with variable breakage rates by Morrell and Morrison (1996).[1][2][3]
The formulation is referred to in the associated literature as the "Variable Rates" model (Morrell et al., 2001).[4]
Model theory
Calculation sequence
Breakage rates
Discharge rates
Slurry hold-up
Classification and discharge
Appearance function
High energy
Ball load
Low energy
Combined appearance function
Mill power
Charge properties
Power draw
Internal mesh series
Multicomponent modelling
Additional notes
Breakage rates and mill load
Slurry pool
Excel
The Variable Rates AG/SAG mill model may be invoked from the Excel formula bar with the following function call:
=mdUnit_AGSAG_VariableRates(Parameters as Range, Size as Range, MillNewFeed as Range, OreSG as Range, BallSizing as Range, RConst as Range, OreBreakageParams as Range, Optional MillRecycleFeed as Range = Nothing)
Invoking the function with no arguments will print Help text associated with the model, including a link to this page.
Inputs
The required inputs are defined below in matrix notation with elements corresponding to cells in Excel row () x column () format:
where:
- is the mass flow feed rate of liquids into the mill (t/h)
- is the Specific Gravity or density of liquids in the feed (- or t/m3)
- is an index of the Appearance function to view in the results
- is an index of the Appearance function to view in the results
- is the number of ore types
- is the number of intervals of the external mesh series
- is the number of intervals of the ball mesh series below the top size, including the submesh
- is the size of the external square mesh interval that feed mass is retained on (mm)
- , i.e. descending size order from top size () to sub mesh ()
- is the size of the square mesh interval that balls are retained on (mm)
- is the mass fraction of balls retained on ball mesh series interval (% w/w)
- indicates the array is an optional input parameter, and is set to null if omitted
Results
The results are displayed in Excel as an array corresponding to the matrix notation below:
where:
- is the number of internal computation steps required to converge the load
- is the numerical error of the converged load approximation
- is the flow rate of pulp into the mill (m3/h)
- is the rotational rate of the mill (rpm)
- is the mass of ore solids in the mill (t)
- is the mass of liquids in the mill (t)
- is the mass of balls in the mill (t)
- is the total mass of ore, liquids and balls in the mill (t)
- is product mass flow rate (t/h)
- is the geometric mean size of the internal mesh series interval that mass is retained on (mm)
Example
The images below show the selection of input arrays and output results in the Excel interface.
SysCAD
The sections and variable names used in the SysCAD interface are described in detail in the following tables.
MD_Mill page
The first tab page in the access window will have this name.
| Tag (Long/Short) | Input / Display | Description/Calculated Variables/Options |
|---|---|---|
| Tag | Display | This name tag may be modified with the change tag option. |
| Condition | Display | OK if no errors/warnings, otherwise lists errors/warnings. |
| ConditionCount | Display | The current number of errors/warnings. If condition is OK, returns 0. |
| GeneralDescription / GenDesc | Display | This is an automatically generated description for the unit. If the user has entered text in the 'EqpDesc' field on the Info tab (see below), this will be displayed here.
If this field is blank, then SysCAD will display the unit class ID. |
| Requirements | ||
| On | CheckBox | This enables the unit. If this box is not checked, then the material will pass straight through the mill with no change to the size distribution. |
| NumParallelUnits | Input | The number of parallel, identical units to simulate:
|
| Method | Fixed Discharge | The discharge particle size distribution is user defined. Different distributions can be used for different solids. |
| AG/SAG (Variable Rates) | The Variable Rates AG/SAG mill model (steady-state or dynamic) is used to determine the mill product size distribution. Different parameters can be used for different solids. | |
| Rod Mill (Lynch) | The Lynch rod mill model is used to determine the mill product size distribution. Different parameters can be used for different solids. | |
| Ball (Perfect Mixing) | The Perfect Mixing ball mill model (steady-state or dynamic) is used to determine the mill product size distribution. Different parameters can be used for different solids. | |
| Stirred (Perfect Mixing) | The Perfect Mixing stirred mill model (steady-state or dynamic) is used to determine the mill product size distribution. Different parameters can be used for different solids. | |
| Mill (Herbst-Fuerstenau) | The Herbst-Fuerstenau model is used to determine the mill product size distribution. Different parameters can be used for different solids. | |
| PowerModels | CheckBox | Show alternative mill power model calculations on the Power page. |
| MediaTrajectory | CheckBox | Show mill media rolling, sliding and free flight trajectory computations on the MediaTraj page. |
| MediaStrings | CheckBox | Show media size distributions at recharge equilibrium on the MediaStrings page. |
| Options | ||
| ShowQFeed | CheckBox | QFeed and associated tab pages (eg Sp) will become visible, showing the properties of the combined feed stream. |
| ShowQProd | CheckBox | QProd and associated tab pages (eg Sp) will become visible, showing the properties of the products. |
| SizeForPassingFracCalc | Input | Size fraction for % Passing calculation. The size fraction input here will be shown in the Stream Summary section. |
| FracForPassingSizeCalc | Input | Fraction passing for Size calculation. The fraction input here will be shown in the Stream Summary section. |
| Stream Summary | ||
| MassFlow / Qm | Display | The total mass flow in each stream. |
| SolidMassFlow / SQm | Display | The Solids mass flow in each stream. |
| LiquidMassFlow / LQm | Display | The Liquid mass flow in each stream. |
| VolFlow / Qv | Display | The total Volume flow in each stream. |
| Temperature / T | Display | The Temperature of each stream. |
| Density / Rho | Display | The Density of each stream. |
| SolidFrac / Sf | Display | The Solid Fraction in each stream. |
| LiquidFrac / Lf | Display | The Liquid Fraction in each stream. |
| Passing | Display | The mass fraction passing the user-specified size (in the field SizeForPassingFracCalc) in each stream. |
| Passes | Display | The user-specified (in the field FracForPassesSizeCalc) fraction of material in each stream will pass this size fraction. |
Mill page
The Mill page is used to specify the input parameters for the mill model.
Ore page
This page is used to define the comminution properties of SysCAD species with the size distribution quality in the project.
Results page
This page is used to display the model results.
| Tag (Long/Short) | Input / Display | Description/Calculated Variables/Options |
|---|---|---|
| Results | ||
| Solver | ||
| Iterations | Display | Number of iterations to converge internal load solver. |
| IterationError | Display | Numerical approximation error of internal load solver. |
| Mill Properties | ||
| MillVolume | Display | Internal volume of the mill. |
| MillSpeed | Display | Rotational speed of the mill. |
| MillFeedRate / Feed.SLQv | Display | Volumetric feed rate of pulp into the mill. |
| Mill Contents | ||
| OreMass | Display | Mass of ore (solids with PSD) in the mill. |
| LiquidMass | Display | Mass of liquids in the mill. |
| BallMass | Display | Mass of ball media in the mill. |
| TotalChargeMass | Display | Total mass of ore, liquids and balls in the mill. |
| VolTotalLoad | Display | Volumetric fraction of mill volume of total charge (ore, liquids, balls and void space). |
| Mill Discharge | ||
| m1 | Display | Parameter of the Austin mill holdup relationship. |
| m2 | Display | Parameter of the Austin mill holdup relationship. |
| dMax | Display | Maximum discharge rate of load volume through the grate. |
| Charge Properties | ||
| S20 | Display | Size of the top (largest) 20% of the load. |
| ChargeDensity | Display | Density of the charge. |
| U | Display | Fraction of charge void space filled with slurry. |
| ThetaShoulder | Display | Angular position of the charge shoulder. |
| ThetaToe | Display | Angular position of the charge toe. |
| ChargeSurfaceRadius | Display | Radius of the inner charge surface. |
| Power | ||
| NoLoadPower | Display | No-load power draw of the mill. |
| NetPower | Display | Net power draw of the mill. |
| GrossPower | Display | Gross power draw of the mill. |
RiDi page
This page displays the breakage and discharge rates for each size interval computed by the model.
| Tag (Long/Short) | Input / Display | Description/Calculated Variables/Options |
|---|---|---|
| Rates | ||
| Size | Display | Size of each interval in internal mesh series. |
| MeanSize | Display | Geometric mean size of each interval in internal mesh series. |
| R | Display | Value of breakage rate, , for each size interval, for each ore species. |
| D | Display | Value of discharge rate, , for each size interval. |
| Ecs | Display | Value of the specific comminution energy for each size interval. |
Load page
This page displays information about the balls, solids and liquids that currently comprise the mill load.
| Tag (Long/Short) | Input / Display | Description/Calculated Variables/Options |
|---|---|---|
| Distribution | ||
| Name | Display | Shows the name of the SysCAD Size Distribution (PSD) quality associated with the feed stream. |
| IntervalCount | Display | Shows the number of size intervals in the SysCAD Size Distribution (PSD) quality associated with the feed stream. |
| SpWithPSDCount | Display | Shows the number of species in the feed stream assigned with the SysCAD Size Distribution (PSD) quality. |
| Load | ||
| SolidMass / SMt | Display | The mass of solids with the SysCAD size distribution property currently in the mill. |
| LiquidMass / LMt | Display | The mass of liquids currently in the mill. |
| BallMass / BMt | Display | The mass of ball media in the mill. |
| Size | Display | Size of each interval in the external mesh series. |
| MeanSize | Display | Geometric mean size of each interval in the external mesh series. |
| Load | Display | The mass of solids with the SysCAD size distribution property currently in the mill, by size and species.}} |
Power page
This optional page displays the inputs and results for alternative mill power models. The page is only visible if PowerModels is selected on the MD_Mill page.
| Tag (Long/Short) | Input / Display | Description/Calculated Variables/Options |
|---|---|---|
| Power | ||
| HoggFuerstenau | CheckBox | Shows inputs and results for tumbling mill power calculations using the Hogg and Fuerstenau method. |
| MorrellC | CheckBox | Shows inputs and results for tumbling mill power calculations using the Morrell Continuum method. |
| MorrellE | CheckBox | Shows inputs and results for tumbling mill power calculations using the Morrell Empirical method. |
| MorrellD | CheckBox | Shows inputs and results for tumbling mill power calculations using the Morrell Discrete Shell method. |
| HildenPowell | CheckBox | Shows inputs and results for tumbling mill power calculations using the Hilden and Powell method. |
MediaStrings page
This page displays the inputs and results for grinding mill media string calculations. The page is only visible if MediaStrings is selected on the MD_Mill page.
MediaTraj page
This page displays the inputs and results for tumbling mill media trajectory calculations. The page is only visible if MediaTrajectory is selected on the MD_Mill page.
About page
This page is provides product and licensing information about the Met Dynamics Models SysCAD Add-On.
| Tag (Long/Short) | Input / Display | Description/Calculated Variables/Options |
|---|---|---|
| About | ||
| HelpLink | Opens a link to the Installation and Licensing page using the system default web browser. Note: Internet access is required. | |
| Information | Copies Product and License information to the Windows clipboard. | |
| Product | ||
| Name | Display | Met Dynamics software product name |
| Version | Display | Met Dynamics software product version number. |
| BuildDate | Display | Build date and time of the Met Dynamics Models SysCAD Add-On. |
| License | ||
| File | This is used to locate a Met Dynamics software license file. | |
| Location | Display | Type of Met Dynamics software license or file name and path of license file. |
| SiteCode | Display | Unique machine identifier for license authorisation. |
| ReqdAuth | Display | Authorisation level required, MD-SysCAD Full or MD-SysCAD Runtime. |
| Status | Display | License status, LICENSE_OK indicates a valid license, other messages report licensing errors. |
| IssuedTo | Display | Only visible if Met Dynamics license file is used. Name of organisation/seat the license is authorised to. |
| ExpiryDate | Display | Only visible if Met Dynamics license file is used. License expiry date. |
| DaysLeft | Display | Only visible if Met Dynamics license file is used. Days left before the license expires. |
See also
External links
- Met Dynamics - Mill (help.syscad.net)
- Crushing and Grinding Example project (help.syscad.net)
- Crushing, Grinding, and Flotation Example project (help.syscad.net)
References
- ↑ Napier-Munn, T.J., Morrell, S., Morrison, R.D. and Kojovic, T., 1996. Mineral comminution circuits: their operation and optimisation. Julius Kruttschnitt Mineral Research Centre, Indooroopilly, QLD.
- ↑ Leung, K., Morrison, R.D. and Whiten, W.J., 1987. An Energy Based Ore Specific Model for Autogenous and Semi-autogenous Grinding, Copper 87, Vina del Mar, Vol. 2, pp 71 - 86
- ↑ Morrell, S. and Morrison, R.D., 1996. AG and SAG mill circuit selection and design by simulation. In International Conference on Autogenous and Semiautogenous Grinding Technology (Vol. 2, pp. 769-790).
- ↑ Morrell, S., Valery, W., Banini, G. and Latchireddi, S., 2001. Developments in AG/SAG mill modelling. Proceedings of Autogenous and Semiautogenous Grinding Technology, Vancouver, pp.71-84.





