Chapter 8: Tutorial - Cone Crusher

(Part A) Define the mantle and shaft movements required for a Cone Crusher, set up Tavares breakage parameters, and collect Boundary Collision Statistics for later post-processing.

(Part B) Track the position of tagged particles, compute the resulting PSD of resulting fragments, isolate broken fragments from whole particles, and calculate the average power draw.

8.1. Part A: Project Setup and Processing

8.1.1. OBJECTIVES

The two main purposes of this tutorial are to learn how to define the mantle and shaft movements required for a Cone Crusher simulation, and how to set up Tavares breakage model parameters.

You will learn how to:

  • Turn on the collection of intensities data for later post-processing of power draw

  • Define mantle and shaft movements for a cone crusher simulation

  • Turn on random particle orientation

  • Set up Tavares breakage model parameters

And you will use these features:

  • Boundary Collision Statistics Module

  • Cone Crusher Motion Frame

  • Particle Orientation

  • Tavares Particle Breakage

8.1.2. PREREQUISITES


Important:  This ADVANCED tutorial contains fewer details, screenshots, and procedures than other Rocky tutorials.

  • An ADVANCED tutorial is designed for users who are more familiar with the Rocky user interface (UI), and already have a good understanding of the common setup and post-processing tasks.

  • If you do not already have this level of familiarity, it is recommended that you complete at least Tutorials 01- 05 before beginning this one.


8.1.3. GEOMETRY

 

The geometries in this tutorial are composed of:

  • (1) Feed Hopper (Choke Feed)

  • (2) Housing

  • (3) Transmission Axis

  • (4) Mantle Cone

  • (5) Mantle Shaft

  • (6) Concave

The .stl files can be found in the tutorial directory.

8.1.4. PROJECT CREATION

To get started with this tutorial, do the following:

  1. Download the dem_tut08_files.zip file here .

  2. Unzip dem_tut08_files.zip to your working directory.

  3. Open Rocky 2025 R2.

  4. Create a new project.

  5. Save the empty project to a location of your choosing.

8.1.5. PROJECT SETUP

For this tutorial, information about Power on the boundary will be used.

Therefore, we must enable the collection of Boundary Collision Statistics during the Modules step.

  1. Use the information in the table that follows to start setting up your Rocky project.


    Tip:  If you run into settings or procedures in these tables that you are not yet familiar with, please refer to the Rocky User Manual and/or other Tutorials (via the Introductory Tutorials and Advanced Tutorials) to find the detailed instructions you need.


    StepEntityEditors LocationParameter or ActionSettings
    AStudy 01StudyStudy NameCone Crusher
    BModulesModulesBoundary Collision Statistics(Enabled)
    CModules

    ﹂Boundary Collision Statistics

    Boundary Collision StatisticsIntensities(Enabled)


Note:  With Intensities enabled, Rocky will collect the average dissipation and impact power values measured by each individual geometry triangle.


8.1.6. GEOMETRY IMPORT

For the Geometries step, we will import the following six separate geometry files in .stl format, disable the Feed Hopper geometry during initial particle release, and then add a surface (to associate with an inlet) to release particles into the domain.

 

  1. Use the information in the table that follows to continue setting up your Rocky project.

    StepData EntityEditors LocationParameter or ActionSettings
    AGeometriesImport WallAll six .stl geometry files with "mm" for Import Unit
    BGeometries

    ﹂Feed Hopper

    WallEnable Time3 [s]
    CGeometriesCreate a Circular Surface
    DGeometries

    ﹂Circular Surface <01>

    Circular SurfaceCenter Coordinate0, 1.65, 0 [m]
    Max Radius0.46 [m]

8.1.7. MOTION FRAME CONE CRUSHER

In this tutorial, we will be creating cone crusher movements, which are defined by the following motions:

  • The Mantle rotates eccentrically around the vertical axis, pressing the particles on the Concave and promoting breakage. At the same time, it is free to rotate around the Shaft.

  • The intersection between the Mantle Axis and the Rotation Axis characterizes the Pivot Point.

  • It is important to have your mantle geometry rotated by the desired Mantle Angle and to know in which direction this is done to properly set up the Cone Crusher movement.

These motions can be created in Rocky with a special Cone Crusher Frame that is added during the Motion Frames step.

 

When the new Cone Crusher Frame is added to your project, the following options are available in the Data Editors panel:

  • Pivot Point: Coordinate of the point around which the Mantle and Shaft will pivot.

  • Rotation Axis: The axis around which the Mantle and Shaft are rotated.

  • Rotational Velocity: The angular velocity along the Rotation Axis.

  • Initial Orientation: The vector about which the Mantle and Shaft begin their rotation.

  • Start/Stop Time: Time in which the motion starts and stops.

In addition, the Cone Crusher Frame will automatically create two separate motions that are designed to be applied to the appropriate geometries.


Important:  These motions appear only in the Motion Frame list for the imported Geometry and are not listed under Motion Frames in the Data panel.


These two motions include the following:

  • Cone Crusher <01> (Mantle) for the Mantle Cone geometry.

  • Cone Crusher <01> (Shaft) for the Mantle Shaft geometry.

 

How to calculate Initial Orientation:

  • Mantle Angle:

    (8–)


  • Initial Orientation:

    (8–)


For the Motion Frames step, we will add a type of motion frame that is specific to the cone crusher:

  1. From the Data panel, right-click Motion Frames, and then select Create Cone Crusher Frame.

     

  2. Use the information in the table below to define the new frame.

    StepData EntityEditors LocationParameter or ActionSettings
    AMotion Frames

    ﹂Cone Crusher <01>

    Cone Crusher FramePivot Point0, 0.687, 0 [m]
    Rotational Velocity15.708 [rad/s]
    Initial Orientation0.03943, 0.9992, 0 [m]
    Start Time2.5 [s]

The frame's two components can now be assigned to their respective geometries.

  1. Use the information in the table that follows to assign the new motions.

    StepData EntityEditors LocationParameter LocationSettings
    AGeometries

    ﹂Mantle Cone

    WallMotion FrameCone Crusher <01> (Mantle)
    BGeometries

    ﹂Mantle Shaft

    WallMotion FrameCone Crusher <01> (Shaft)

At this point, the movement can be previewed.

  1. From the Data panel, click Motion Frames and then click Preview.

    A new Motion Preview window will appear. To better see the results:

  2. Use the Data panel eye icons to hide all the geometries except Mantle Cone and Mantle Shaft.

  3. Also, to better see the two separate frame axes, enable the Transparency checkbox for both the Mantle Cone and Mantle Shaft geometries.


Note:  When you play the preview, no motions will be seen before 2.5 s.


 

8.1.8. MATERIAL DEFINITION

For the Materials step, two materials will be used:

  • One for all wall parts (Default Boundary), which we will use as defined by default.

  • Another for the particles (Default Particles), which we will modify.

  1. Use the information in the table below to define these materials and their interactions.

    StepData EntityEditors LocationParameter or ActionSettings
    AMaterials

    ﹂Default Particles

    MaterialUse Bulk Density(Cleared)
    Density2500 [kg/m3]
    Young's Modulus1e+07 [N/m2]
    BMaterialsMaterials Interactions |

    Default Particles ⯆

    Default Boundary ⯆

    Static Friction0.7 [s]
    Dynamic Friction0.7 [s]


Important:  The Young's Modulus value was lowered to speed up tutorial processing. For breakage simulations, it is recommended to set these values between 5e+08 and 1e+09 Pa.


8.1.9. PARTICLE GROUP DEFINITION

For the Particles step, we will create a new (rock-like) polyhedron-shaped particle group, in a range of sizes, with random orientation, and with the Tavares breakage model enabled.

  1. Use the information in the table below to start defining your particles.

    StepData EntityEditors LocationParameter or ActionSettings
    AParticlesCreate Particle
    BParticles

    ﹂Particle <01>

    ParticleShapePolyhedron
    Particle | SizeAdd size row (x3)
    (1) Size | Cumulative %0.2 [m] @ 100 [%}
    (2) ...0.15 [m] @ 70 [%]
    (3) ...0.1 [m] @ 30 [%]
    (4) ...0.08 [m] @ 10 [%]

8.1.10. RANDOM PARTICLE ORIENTATION

For this tutorial, we will change the orientation of the particles as they are injected into the inlet.

 

  • By default, particles are injected into the domain perpendicular to the inlet (top image).

  • You can change the orientation to have all particles come in at a defined angle (middle image), come in at random angles (bottom image), and many options in between.

  • For this tutorial, we will set a random orientation, as specified below.

8.1.11. PARTICLE GROUP DEFINITION

To set the orientation of the particle:

  1. From the Orientation sub-tab, enable the Random Angle checkbox, and then define the Half Range (x,y,z) values (as shown).

     


Tip:
  • For default particle shapes, or other symmetrical custom shapes, setting your values to 90, 90, 90 ensures a completely random orientation in all directions.

  • For non-symmetrical custom particle shapes (such as a bottle), you will need to set your values to 180, 180, 180 to ensure a completely random orientation.


8.1.12. TAVARES BREAKAGE MODEL

For this tutorial, we will enable the Tavares breakage model.

Like Ab-T10, the Tavares model is based upon the Voronoi fracture particle subdivision algorithm.

However, the breakage energy probability and resulting fragment size distribution are based upon the Tavares et al. (UFRJ) approach, which:

  • Models fractures by low-energy stressing, which is most relevant when particles are subject to a complex series of loading events, as within crushers.

  • Has been able to describe the progressive growth of crack-like damage that ultimately leads to the fracture of a particle under stresses significantly lower than those required for breakage in a first event.

 


Tip:  Refer to the Tavares Breakage Model Description and Calibration Guidelines white paper for more information.


8.1.13. PARTICLE GROUP DEFINITION

To set the Breakage Model:

  1. From the Breakage sub-tab, select the Enable Breakage checkbox.

  2. From the Breakage | Criteria tab, define the Model and leave all the model's parameters as default.

     

  3. From the Breakage| Fragments sub-tab, define Minimum Absolute Size and Distribution model.

     


Important:  The Minimum Absolute Size value used in this tutorial was increased to speed up processing time. Refer to the Fragment size distribution models section of the DEM Technical Manual for details on how best to set this value.


8.1.14. MASS FLOW DEFINITION

For the Inlets and Outlets step, we will create a Particle Inlet and then set our Circular Surface as the location from which we want particles to enter the simulation.

For the Solver step, we will also define our breakage start and delay times.

  1. Use the information in the table below to finish setting up your project.

    StepData EntityEditors LocationParameter or ActionSettings
    AInlets and OutletsCreate Particle Inlet
    BInputs

    ﹂Particle Inlet <01>

    Particle InletEntry PointCircular Surface <01>⯆
    Particle Inlet |

    Particles

    Add row (x1)
    (1) Particle | Mass Flow RateParticle <01>⯆ @ 350 [t/h]
    … | TimeStop3 [s]
    CSolverSolver | TimeSimulation Duration12 [s]
    Breakage | Start2.5 [s]
    … | Delay After Release1 [s]
    Solver | GeneralSimulation TargetCPU ⯆

8.1.15. SETUP CONFIRMATION

With a 3D View window opened, your Data panel and Workspace should look similar to the image below.

 

8.1.16. SIMULATE PROJECT

  1. From the Solver entity, click Start.

The Simulation Summary screen appears (as shown), then processing begins.

 

 


Tip:  You can use the Auto Refresh checkbox to view in a 3D View window the results during processing.


8.1.17. CONCLUSION

This completes Part A of this tutorial, in which Rocky was used to set up and process a cone crusher simulation.

During this tutorial, it was possible to:

  • Collect Boundary Collision Statistics data for later post-processing of power draw.

  • Use the specialized Cone Crusher Frame to easily set up and apply specific mantle and shaft motions.

  • Set particles to be released in a random orientation.

  • Enable the simulation of Breakage phenomena using the Tavares model.

What's Next?

  • If you completed this tutorial successfully, then you are ready to move on to Part B and post-process this project.

8.2. Part B: Post-Processing

8.2.1. OBJECTIVES

The main purpose of this tutorial is to learn how to analyze the results from the Cone Crusher simulation that we set up and processed in Part A.

You will learn how to:

  • Track the position of tagged particles

  • Compute the PSD of the resulting fragments

  • Isolate the broken fragments from the whole particles

  • Calculate the Average Power Draw

And you will use these features:

  • Histograms

  • Time Plots

  • User Processes, including:

    • Cylinder

    • Cube

    • Particles Time Selection

    • Property

  • Particle Calculations, including:

    • Divisions Tagging

8.2.2. PREREQUISITES


Important:  This ADVANCED tutorial contains fewer details, screenshots, and procedures than other Rocky tutorials.

  • An ADVANCED tutorial is designed for users who are more familiar with the Rocky user interface (UI), and already have a good understanding of the common setup and post-processing tasks.

  • If you do not already have this level of familiarity, it is recommended that you complete at least Tutorials 01- 05 before beginning this one.


8.2.3. OPEN PROJECT

  1. If you completed Part A of this tutorial, ensure that Rocky project is open. (Part B will continue from where Part A left off.)

  2. If you did not complete Part A, do all of the following:

    1. Download the dem_tut08_files.zip file here .

    2. Unzip dem_tut08_files.zip to your working directory.

    3. Open Rocky 2025 R2. (Look for Rocky 2025 R2 in the Program Menu or use the desktop shortcut.)


      Important:  To make use of the Rocky project file provided, you must have Rocky 2025 R2 or later. If you have an earlier version of Rocky, please upgrade Rocky to the latest version, or complete Part A from scratch.


    4. From the Rocky program, click the Open Project button, find the tutorial_08_A_pre-processing folder, and then open the tutorial_08_A_pre-processing.rocky file.

    5. Process the simulation. (From the Data panel, select Solver and then from the Data Editors panel, click the Start button.)

8.2.4. POST PROCESS PSD HISTOGRAM

Once processing is complete, you can analyze the resulting Particle Size Distribution (PSD) of the fragments.

We will do this by:

  • Using a Cube User Process to define the area beneath the cone.

  • Applying a Particles Time Selection User Process to the Cube in order to capture data for the whole simulation length.

  • Plotting the resulting Particle Size data in a Histogram.

 

  1. Use the information in the table below to set up the Cube and Particles Time Selection User Processes.

    StepData EntityEditors LocationParameter or ActionSettings
    AParticlesCreate a Cube User Process
    BUser Process

    ﹂Cube <01>

    CubeCenter0, 0.08, 0 [m]
    Magnitude1.4, 1.1, 1.4 [m]
    CUser Processes

    ﹂Cube <01>

    Create a Particles Time Selection User Process
    DUser Processes

    ﹂Particles Time Selection <01>

    Time SelectionDomain RangeAll ⯆


    Tip:  If you run into settings or procedures in these tables that you are not yet familiar with, please refer to the Rocky User Manual and/or other Tutorials to find the detailed instructions you need.


  2. Use the information in the table below to create and configure the Histogram.

    StepItemParameter or ActionSettings
    AUser Processes

    ﹂Particles Time Selection <01>

    Show in new Histogram by Particle Size
    BHistogram (window)Configure histogram (button)
    CConfigure Histogram (dialog box)WeightParticle Mass ⯆
    Number of Bins20
    Cumulative Bins(Enabled)
    Percent Values(Enabled)

 

8.2.5. POST PROCESS - DIVISIONS TAGGING

In the last Timestep of the simulation, you can observe that some Particles are accumulating in an open space of the transmission housing (in the Blue cube) (as shown).

Using the Divisions Tagging Calculation of the Particles in the hopper, we can identify from what section of the hopper these particles are coming from, which is helpful information for improving the equipment design.

 

For this case, we will create a Cylinder User Process at an Output Time far enough into the simulation that the Hopper is full of particles.

  1. From the Time toolbar, select the [60] 3 s output time.

  2. Use the information in the table below to create and configure the Cylinder.

    StepData EntityEditors LocationParameter or ActionSettings
    AParticlesCreate a Cylinder User Process
    BUser Processes

    ﹂Cylinder <01>

    CylinderSize1.01, 1.01, 1.01 [m]
    Center0, 1.2, 0 [m]

  3. With the Cylinder set up, right-click Particles, point to Particles Calculations, point to Divisions Tagging, and then click Cylinder <01>.

    This will generate divisions within the cylinder that was selected, the result of which can be seen in a 3D View window.

     

  4. From the Data panel, under Calculations, select the new Divisions Tagging (Cylinder <01>) entry.

  5. Use the information in the table to define the Tagging values.

    StepData EntityEditors LocationParameter or ActionSettings
    ACalculations

    ﹂Divisions Tagging (Cylinder <01>)...

    TaggingTime Ranger Filter | Domain RangeAfter Time
    Time Range Filter | Initial3 [s]
    Tangential Divisions12 [s]


Tip:   With a 3D View window open and selected, you can visualize the tagged particles by enabling the Transparency checkbox from the Coloring tab.


Color Particles using the newly added Property:

  1. From the Data panel, select Particles.

  2. From the Data Editors panel, select the Coloring tab.

  3. Under Nodes, change the Property to Divisions Tagging (Cylinder <01>)


    Tip:  Use the eye icons next to Cylinder and Cube (under User Processes in the Data panel) to hide the Cylinder and Cube coloring and show only the Particles coloring.


  4. Display the final Output Time.

You should now see all of the following:

  • Hopper divided in 12 circular sectors.

  • Particles and Fragments colored by the division they started in.

  • The accumulated Fragments colored by their original position in the hopper.

 


Tip:  You can duplicate a 3D view by pressing Ctrl + D, selecting View > Window, or right-clicking in the 3D view and choosing Duplicate.


Lastly for this analysis, a Cube will be added to separate the accumulated Fragments from the rest of the Particles.

  1. Use the information in the table to define this new Cube.

    StepData EntityEditors LocationParameter or ActionSettings
    AParticlesCreate a Cube User Process
    BUser Processes

    ﹂Cube <02>

    CubeCenter0.77, -0.017, 0 [m]
    Magnitude0.24, 0.24, 0.54 [m]

Then, another Histogram will be created to show the hopper region from which the accumulated Fragments are originating.

  1. Use the information in the table to define the Histogram.

    StepItemLocationParameter or ActionSettings
    AUser Processes

    ﹂Cube <02>

    Properties |

    Divisions Tagging (Cylinder <01>)...

    Show in new Histogram
    BHistogram (window)Configure histogram (button)
    CConfigure Histogram (dialog box)Number of Bins12
    PropertiesDivisions Tagging (Cylinder <01>)...
    LimitsUser Defined ⯆
    Min0.5 [-]
    Max12.5 [-]

 

 

From the Histogram, it can be seen that most of the accumulated Fragments are coming from bins 1 and 12 of the hopper.

This information may not have been obvious before the Divisions Tagging analysis. Now, an appropriate solution can be devised if this particular result is undesirable.

8.2.6. POST PROCESS - ISOLATE FRAGMENTS

It might also be useful to filter out the particles that actually broke; in Rocky, these broken particles are called fragments.

One way to isolate these fragments is by using the Laguerre-Vonoroi Size property, which has the following characteristics:

  • It is measured only for breakage fragments.

  • It is equal to the diameter of a sphere used to generate the fragments using the Laguerre-Voronoi tessellation.

  • This diameter, , can be defined as twice the minimum distance between a particle's center of gravity, , and the particle's sides, , as shown in the following equation:

    (8–1)

Using this property, the particles that did not break will have values equal to zero.

Therefore, to filter out only the fragments, we will select only the particles that have a Laguerre-Vonoroi Size value greater than zero.

8.2.7. POST PROCESS - ISOLATE FRAGMENTS

Let's start by creating the Filter User Process.

  1. Use the information in the table below.

    StepData EntityEditors LocationParameter or ActionSettings
    AParticlesCreate a Filter User Process
    BUser Processes

    ﹂Filter <01>

    PropertyPropertyLaguerre-Voronoi Size ⯆
    TypeRange ⯆
    Minimum Value1e-10 [m]
    Maximum Value0.2 [m]

Then, to visualize in a 3D View window only the filtered fragments:

  1. Select a 3D View window (or create a new one by pressing Ctrl + D).

  2. Ensure the geometries are transparent (or hidden).

  3. From the Data panel, hide Particles and all other Particles Calculations and User Processes (besides Filter <01>) using the eye icons.

  4. From the Data Editors panel, on the Filter <01> entity, select the Properties tab, and then drag and drop Laguerre-Vonoroi Size onto the 3D View window. (Results shown below.)

8.2.8. POST PROCESS - ISOLATE FRAGMENTS

 

The view now shows only the particles that have broken (fragments).

8.2.9. POST PROCESS - POWER DRAW

Earlier in Part A, we turned on the collection of some Boundary Collision Statistics; specifically, those regarding Intensities.

We can now use the resulting data to estimate how much power is required by the equipment, which is commonly known as Power Draw.

The Average Power Draw of the geometries can be estimated by the following expression:

(8–2)

where:

  • is the output time.

  • is the Average Power Draw at the current output time.

  • is the output time in which the operating conditions begin.

  • is the output time in which the analysis ends.

  • is the Output Frequency.

First, let's define .

  1. Use the information in the table below to plot the power curves for the cone and shaft geometries.

    StepData EntityEditors LocationParameter or ActionSettings
    AGeometries

    ﹂Mantle Cone

    Curves | PowerShow Curve in New Plot
    BGeometries

    ﹂Mantle Shaft

    Drag and Drop onto Time Plot <01> window

The plot shows the operation conditions start after 2.5 seconds (in which the breakage starts), so = 2.5s.

 

Now, let's apply the Average Power Draw formula to the Power curve.

  1. Use the information in the table below to continue with the analysis.

    StepItemLocationParameter or ActionSettings
    AGeometries

    ﹂Mantle Shaft

    CurvesAdd new custom curve (button)
    BAdd new

    (dialog box)

    NameAverage Power
    Output unitW
    Inputs | Power(Enabled)
    CCustom Curves (dialog box)Expression

    cumsum(A)*OUTPUT_FREQUENCY/(TIME_ELAPSED-2.5)

The new Average Power (Custom) curve appears on the Curves tab for the selected geometry, but is also made available for all other geometries automatically.

Now, let's sum the Average Power of the Mantle Cone and the Mantle Shaft geometries to estimate the Average Power Draw of the equipment:

  1. Use the information in the table below to continue with the analysis.

    StepItemLocationParameter or ActionSettings
    AGeometries

    ﹂Mantle Cone

    Curves | Average Power (Custom)Show curve in new plot
    BGeometries

    ﹂Mantle Shaft

    Drag and drop onto open Time Plot <02> window
    CTime Plot <02> (window)Table (tab)Add Formula
    DAdd Expression (dialog box)Curve CaptionAverage Power Draw
    Curve ExpressionB+C

 

Each row of this table shows the average power measured from 2.5s to the time shown in Column A.

Thus, the Average Power Draw of the simulation is in the last row of this table and it is 4840.28 W (as shown).


Note:  The values you end up with in your project may vary slightly from the ones shown in this Tutorial.


8.2.10. CONCLUSION

This completes Part B of this tutorial, in which Rocky was used to analyze the results from a cone crusher simulation.

During this tutorial, it was possible to:

  • Use Cube and Particles Time Selection User Processes to plot the PSD of fragments in a Histogram.

  • Use a Cylinder User Process and Divisions Tagging Particle Calculations to track specific particles and fragments through the equipment.

  • Use the Property User Process and the Laguerre-Voronoi Size Property to isolate the particle fragments from the whole particles.

  • Use the Intensities Boundary Collision Statistics data we collected to calculate the power draw.

What's Next?

  • If you completed this tutorial successfully, then you are ready to move on to next tutorial.