EDEM SOP
EDEM SOP
of mill liner
Title: Description
Table of Contents:
2. EDEM
2.1. EDEM Creator
2.2. EDEM Simulator
2.3. EDEM Analyst
Page | 1
Title: Description
The Discrete Element Method (DEM) is a numerical method for modelling the dynamics of solid
particles which interact with each other at discrete contact points. In DEM simulations, the
properties of a stressed assembly of rigid spherical particles (e.g. position, velocity and contact
forces) are updated at every numerical iteration time step. The translational and rotational
displacements of each particle are obtained by explicitly integrating the governing differential
equations based on the Newton’s second law of motion, while the contact forces between particles
are calculated using well-defined force–displacement contact models.
In this process, the interactions between particles are monitored at all contacts and the
movement of each individual particle is traced. It is assumed that the velocity of each particle is
constant within each iteration step.
The simulation starts by detecting the contact points between particles. The contact point
and the particle overlapping distance are then used to calculate the interaction forces via the
force–displacement relationships. The resulting contact forces are applied at the centre of each
particle, causing particles to move. According to the Newton’s second law of motion, the particle
acceleration, velocity and displacement are updated at the end of each numerical iteration step.
This mechanical loop continues until the prescribed total iteration steps are reached.
The accuracy of DEM modelling depends highly on the contact models used to calculate
the particle–particle interactions. The fundamental problem is how to reproduce the correct
interactions between two solid grains in the DEM. The classical DEM model treats all solid grains
as either disks (i.e. two-dimensional (2D) model) or spheres (i.e. three-dimensional (3D) model),
as this would greatly simplify the calculation of the overlap distance between particles. Though this
model has shed some lights on the mechanical behaviour of granular materials, it is inaccurate
due to the oversimplified assumptions that all solid particles are spherical and free to rotate.
To account for the particle shape effects, three ways of DEM modelling can be used. The
first of the three is to use a clump of spherical particles based on the bonded-particle method. In
this method, the real shape of soil/rock mass is approximated by a clump of bonded spherical
grains, which behaves as a rigid body. The second approach is to model the soil/rock blocks as
polygons.
Three distinct particle contact models can be identified, namely, (a) the normal contact
model; (b) the tangential contact model; and (c) the rolling contact model. The mechanical
responses of these three different contact models are closely related to the relative displacement
between the two particles, as shown in Fig.1. The interaction forces between two spheres include
the normal (Fn) and tangential (Ft) forces, and rolling moment (Mr), respectively.
Page | 2
Title: Description
Fig. 1
In Fig.1, the normal contact force (Fn) acting on a single particle is linearly proportional to
the overlap distance between two particles, which is expressed as
Fn = Kn* Un
in which Kn is the normal contact stiffness, Un is the normal overlap distance.
Figure.1b shows that a maximum tangential force exists, before which, the tangential force (F t) can
be calculated incrementally, as
Fn=Fn-1 +Ks*dUs
where Fn and Fn-1 are tangential forces calculated at current and previous iteration steps. K s is the
shear stiffness, and dUs is the incremental tangential sliding displacement.
Fig. 2
Page | 3
Title: Description
2. EDEM
We use EDEM (which is the market-leading software for bulk material simulation) as a tool ,
powered by state-of-the-art Discrete Element Modelling (DEM) technology, EDEM quickly
and accurately simulates and analyzes the behaviour of bulk materials such as coal, mined
ores, soil, tablet, and powders.
EDEM software is made of 3 core components: EDEM Creator, Simulator and Analyst
Page | 4
Title: Description
Fig. 01
Fig. 02
Enter the title and write the description of the simulation for reference as shown in fig. 2
Page | 5
Title: Description
3. Rename the material for proper identification. Select the material and right click; you will
get the following display on your screen as shown in fig. 3. Select <Rename Material>
and name the material as per your convenience and click <OK>.
Fig. 03
4. Once you have renamed the material (e.g. Ore), select the same material and you will
see the following display on your screen.
5. Enter the values of the Poisson’s ratio, density and any one out of shear and Young’s
modulus according to the material’s properties as per the table given below
Page | 6
Title: Description
(NOTE: In all the cases we enter the value of shear modulus instead of Young’s modulus)
6. If you want to define more than one type of bulk material then simply follow the same
procedure (from step 1 to 5) as given above for adding the new material and renaming
it.
(NOTE: In the case of mill shell region simulation we create ball and ore as two different
materials)
7. Adding particle geometry (shape and size).
Select the bulk material that you have defined and right click, you will get the following
display on your screen.
Page | 7
Title: Description
11. For the single sphere prototype, enter its radius in mm by selecting the corresponding
cell as shown in the fig. 4
Fig. 04
12. Select <Size Distribution>, you will get the following display on your screen
Fig.05
13. The default size distribution scheme is of type “fixed “where all the particles are of same
size. To change the size distribution scheme , click on the dropdown as shown in fig. 6
Fig. 6
14. You will see various distribution schemes available but most often we select the <user
defined>.
(NOTE: In “USER DEFINED”, we can define the variation in size of the particles of the bulk
material with their relative composition.)
Page | 8
Title: Description
15. You will get the following display on the screen with the given default values.
Fig. 7
16. To add particles of different size, select the scale cell of default particle and right click to
<Add>the number of size distributions as per your requirement. In this case we have
considered 3 size distributions as shown in the fig. 8
Fig. 08
Here Scale 1 means the original particle size as defined in step 11. The other scales will be
defined by a reducing factor of √2 of the base scale. The %age of mass against each size
distribution should be defined accordingly so that the total sum always adds to 100.
(NOTE : Here we have shown scaling by radius but scaling by volume can be used as well.)
17. Follow the same procedure for introducing further different sized particles and defining
their percentage mass composition in the bulk material.
18. Select <Properties> as shown in the fig. 9. Then check the box against “Auto
Calculation”. This process will position the particle about its centre of mass, and it will
also populate the Mass, Volume, and Moment of Inertia boxes of the particle details.
Page | 9
Title: Description
Fig. 09
19. Select <Equipment Material> and right click, you will get the following display.
20. Select <Add Equipment Material>to define the material of the equipment and follow
the same procedure for renaming and assigning its physical properties as done in the
case of bulk material.
Page | 10
Title: Description
(NOTE: This list also includes the first material that you selected because you need to
define how particles will behave when they encounter other particles or equipment that is
made out of the same material.)
Fig. 10
5. You can continue to add interactions to your simulation by pressing the green plus
button.
(NOTE: Once you have created an interaction pair for a combination of materials, EDEM
will remove the option to choose those materials from the Create Interactions list).
(NOTE: EDEM will recognize that the interaction between material 1 and material 2 is the
same as the interaction between material 2 and material 1 and so you need to create this
interaction once)
6. If you click the green plus
button even after all
interactions have been
defined , the following
window will get appear on
your screen saying no
further interactions can be
defined.
Page | 11
Title: Description
7. Similarly, define the interactions of the other remaining materials that you have defined
under the “Bulk Material” and “Equipment Material” section.
In our case, the liner has the effective diameter of 5000mm and length 300mm, shell plate
has a thickness of 60mm and lifter bar has 210-180 SDS profile as shown below. The axis
of the liner is along the Z axis.
Fig. 12
Page | 12
Fig. 13
Title: Description
Fig. 14
This ensures that the part is created at the right scale in EDEM.
5. The next option <Merge Sections> is checked to combine together all the geometry in
a CAD file as a single EDEM geometry object, rather than individual parts. This is done
to ensure that each part of the assembly has the same kinematics.
In our case, we have merged the
geometry.
Fig. 15
7. The EDEM will create the geometry with a default name as shown in the fig. 16
To rename the geometry, select that geometry and right click and select <Rename
Geometry>
Fig. 16
Page | 13
Title: Description
8. A new window will get appear (as shown in the fig.17) where you can enter the name of
the geometry as per your convenience and click <OK>.( In our case, it is “liner” )
Fig. 17
9. Select your geometry and you will see a dropdown with two options namely
(NOTE: When a geometry part is set to Physical it means that it has a material assigned to
it and when particles come into contact with it, they will react according to the contact model
you have assigned.
When geometry is defined as Virtual, particles will be free to pass
through the surfaces as though they were not present. Note since the particles will not
interact with a virtual geometry, there is no option to assign material to a part when Virtual
has been selected. )
10. For liner, select <Physical> as it will interact with the bulk material and then select its
material accordingly.
11. To add kinematics to the liner, select <liner> from the Geometries section, right click
and then select <Add Kinematic> as shown in the fig. 18
Fig. 18
Page | 14
Title: Description
Fig. 19
13. Select <New Kinematic 1> and define the specifics of the motion on the details view.
Static factory: It will create all the requested particles in a single time-step of EDEM
simulation. Once the requested volume of material has been produced, the Factory will be
classed as complete and EDEM will move on to processing the next time-step. Static
Factories are useful for quickly filling volumes.
Dynamic factory: It acts over multiple time-steps and allows you to create material
continuously. Dynamic factories are useful for representing a feed of material into a system.
1. Select <Geometries> and right click and then select <Add Geometry>. You will see a
list of options as shown below
Page | 15
Title: Description
(You can select any of the three given geometry for the factory. In our case, we have
selected ‘Box‘.)
2. Select <Box>, you will have the following display on your screen
3. The factory will have the default name ‘ New Section 1’. You can rename it by right
click on that factory and then select <Rename Geometry>.
4. A new window will appear, then enter the name of the factory and click <OK>.
5. Select your specified factory and click on the ‘Type’ dropdown box and then select
<Virtual>.
(NOTE: Virtual because we don’t want it to interact with the particle or equipment)
6. Select <Box> and you will see the size panel of the box where you can define its centre
coordinates, its length along the three principle axes and its rotation about those axes.
Page | 16
Title: Description
In our case, the effective diameter of the liner is 5000mm. So, we have
increased the length along x axis from default 1000mm to 4000mm as shown in the fig.20
Fig. 20
We then changed the viewing direction of the geometry from positive X axis to positive Y
axis by clicking on the ‘View ‘dropdown present on the toolbar and selecting <+ Y> as
shown below
As the dimension of the factory along Z axis was exceeding the dimension of our geometry,
as shown below,
So, we decreased its length along Z axis from default 1000mm to 250mm and we got the
following display on the screen
7. Select your specified factory and right click and then select <Add Factory>, you will
have the following display on your screen
Page | 17
Title: Description
Fig. 21
( NOTE: Please ensure in the’ Material’ dropdown box, the desired bulk material is
selected. In the above figure, it is ore.)
‘Unlimited Number’s used when unlimited amount of bulk material is required to be created
by the factory.
For limited case, there are two options,
Total Number: It is used when you want to create a certain fixed number of particles of a
specified material.
Total Mass: It is used when you want to create a certain fixed mass of a specified material.
In our case, we have checked ‘Total Mass’ checkbox (Highlighted in the fig. 21) and
entered the mass and the discharging rate of the ore based on the following calculations:
Page | 18
Title: Description
Simulation domain: It is the region in which particles can exist within an EDEM simulation.
It is represented in the EDEM viewer window as the red box outline
All calculations are performed on objects within the Domain and should a particle pass out
of the Domain wall it will be removed unless periodic boundaries are turned on.
Gravity Force: By default gravity is set to a value of 9.81 meters per second squared
acting in the negative Z direction. You can change the direction gravity is acting by
adjusting the values in the boxes.
Periodic Boundaries: For some applications, it is desirable to have particles that exit one
side of the domain to reappear on the opposite side in the same position and with the same
velocity. This behaviour can be achieved using the Linear Periodic Boundary settings.
2. If you want to change the size of the Domain, deselect <Auto Update from
Geometry>and manually enter the domain size.
In our case, we had the following default domain size as seen from the positive Y direction.
Page | 19
Title: Description
After deselecting <Auto Update from Geometry>and entering the desired values in the
Domain panel (ensuring that the Domain completely encloses the factory, though it may or
may not enclose the geometry), we got the following display on our viewer window.
(NOTE: The Domain should always enclose the factory completely otherwise some of the
particles generated by the factory would get outside the Domain and would not be
considered in the simulation.)
3. In the Gravity panel, you
can activate (or deactivate)
gravity by selecting (or
deselecting) ‘Gravity’ and
then select the direction in
which it will act and enter its
magnitude.
Page | 20
Title: Description
Fig. 22
5. Click the check box that corresponds to the sides of the domain where periodic
boundaries are required(Highlighted in the fig. 22)
In our case, it is Z. Now when particles leave the domain from the selected sides, they will
be recreated on the opposite side rather than being deleted.
Once you have completed the set-up of your simulation in the EDEM Creator, you can save
your file and switch to the EDEM Simulator to begin processing.
Time Integration: At the top of the Simulator panel are the time-step options. By default,
EDEM uses Euler time integration scheme.DEM uses a Rayleigh Time-step approach to
ensure stability is achieved during processing. The Rayleigh time-step is the maximum
theoretical time-step that can be used for a simulation. It is based on the material properties
used and the minimum size that is present in your simulation set-up.
Page | 21
Title: Description
To ensure stability we recommend that you set the time-step to be between 18 to 22% of
the Rayleigh time-step. But it should never be more than 40% as mathematical instability
would incur.
Simulation Time: In this section, the
total run time of the simulation is
defined.
Data Save: It is where you can define
the rate at which results will be saved.
For ex., if you are running for 10
seconds and are saving at a rate for
0.1 seconds then when you analyze
your results you will have 100 write-out
points available.
Cell size: In this section, define the
cell size in terms of the radius of the
minimum sized particle. In most cases,
3R gives the optimum result.
Number of CPU Cores: You can
decide how many CPU cores to be
engaged in the simulation. Greater the
number of CPU cores, faster will be the
calculation.
In our case (Highlighted in the figure
above),
Deselect <Auto Time Step> and
enter against ‘Fixed Time Step’: 22%.
Enter against ‘Total Time’: 10s
Enter against ‘Target Save Interval’: 0.01.
Enter against ‘Cell Size’: 3R.
Click on the Play button on the bar below the viewer window. During processing you will
see the progress bar increasing which gives you the indication of how far the simulation
has been done.
Page | 22
Title: Description
Fig.23
You can also vary the opacity (Highlighted in fig.23) of a piece of equipment, and if you
want to make geometry object more transparent then you can reduce the opacity value.
Colouring: You can adjust the colouring of any item in your EDEM simulation using the
options that are available in the display section of the Analyst tree. You can colour your
material and equipment according to specific properties that relate to their current
behaviour. For ex., you can colour geometry on the basis of force that is acting on them, or
particles according to their velocity.
To view a particle (in our case, it is “ore”) according to their velocity, follow the
following steps:
1. Select the particle whose display you want to modify.
In the present case, we have shown for “ore’.
2. Click on the <Color by> dropdown (Highlighted in fig.
24) and you will see a list of properties based on
which you can modify the display.
3. Select <Velocity> and click on the <Component >
dropdown, you will see a list of velocity features
based on which the colour of your display will be
modified but in the present case, we have selected
<Magnitude>.
4. Check the <Auto Update> checkbox for <Min
Value> and <Max Value>. (Highlighted in the fig.24).
5. If you want to see the legend of the specified
property, then check the <Show legend> checkbox
and then click <Apply >. (Highlighted in the fig.24).
Page | 23
Fig. 24
Title: Description
2. Click on the ‘Record Animation’ button present at the bottom right corner of the Analyst
window. (Highlighted in the figure below)
3. This will open an ‘Export Video and Images’ window where you can configure the
animation options. Select the export format < Windows Media Video> format as shown
highlighted in the fig.25.
Force Analysis: One can find the forces acting on a particular lifter bar at different time-
steps on EDEM Analyst.
1. For that, one should not merge the complete geometry
while importing it (as done in the 5th step of section 3.1.4)
but should leave the <Merge Sections> checkbox
unchecked and click <OK>.You will then see a list of entries
under the <Geometries> section where each part of the geometry is given a default
name.
2. Select the part from the list which corresponds to your desired lifter bar and name it as
per your convenience.
Page | 24
Title: Description
In our case, our desired lifter bar (Highlighted as red in the figure below) has the default
name “Solid2 (6) “. So, we renamed it to “lifter1”.
3. Select any geometrical part apart from <lifter1> and right click and select <Merge
Sections >. A new window will appear showing the list of different geometrical parts
including “lifter1”
4. Select all the parts except <lifter1> and click <OK>. This will merge all the parts except
“lifter1“together.
5. The merged geometry will get a default name. You can name it as per your
convenience. (We have named it to “total liner “.)
You can export the data containing the value of force acting on your specified lifter bar at
different time-steps by following the below given steps:
1. Select <File>, then <Export> and then <Results Data …> (as shown in fig.26). A new
window will get appear as shown in fig. 27.
Page | 25
Title: Description
Fig. 26
9. Click on the <Section>dropdown (Highlighted in fig.27) and select <lifter1 > as in the
present case we want forces acting on “lifter1”.
10. Click on the <Query type> dropdown (Highlighted in the fig.27) and select <Maximum>
because we want data of maximum force acting on “lifter1”
11. Click <Export> as highlighted in the fig.2.
You can also obtain the graph showing the variation of force acting on “lifter1” against time
on EDEM Analyst by following the below given steps :
Page | 26
Title: Description
Fig.28
Page | 27
Title: Description
Fig. 29
1. In the “Bulk Material “section, define the material “Ore “(as highlighted in the fig. 29).
Since steel balls do not enter the discharger system, so no need to define it under this
section.
2. In the “Equipment Material” section, define the material “Discharger system” having the
same physical properties as that of steel.
3. Define the interaction of “Ore” as done before in the section 3.1.3. There will be only two
interactions namely Ore-Ore and Ore-(Discharger System).
4. Right click on <Geometries> and select <Import Geometry> and open the desired file
of the discharger system geometry. A new window will get displayed as shown in the fig.
30
Fig.30
5. Press the <Choose Units> dropdown and select <Millimeters >. Leave the <Merge
Sections> checkbox unchecked and click <OK>.
In the case of discharger system, the factory is not created separately like we did in the
case of shell region. But it is made as a part of the geometry of the discharger system.
6. Select the part which corresponds to your factory and rename it to “factory” as shown in
the fig. 31.
Page | 28
Title: Description
Fig. 31
9. Select the type and the material of the discharger system as highlighted in the fig. 31.
10. Select <factory> and press on <Type> dropdown and select <Virtual>.
11. Right click on <factory> and select <Add Factory>.A new factory will get created with
the default name “New Factory 1” but the default type of any factory on EDEM is
dynamic.
12. To change the factory type, right click on <New Factory 1>and select <Change
Factory Type> as shown in the fig. 32.With this, the factory type will change from
dynamic to static.
Page | 29
Title: Description
Fig. 32
In the case of discharger system, the static factory is required to create a specific amount of
bulk material in a single time-step. The amount and the time instances depend on the
angular speed of the discharger system, number of dischargers and the capacity of the
discharger system and can be determined by the below given relationships :
In our case, the angular velocity is 11.58 rpm, capacity is 34TPH and the number of
dischargers is 24. Mass to be discharged per minute = (34 * 1000) / 60 = 566.6667 kg;
Mass to be discharged per revolution of the mill = 566.6667 / 11.58 = 48.93 kg;
Mass to be discharged per discharger per revolution of mill = 48.93 / 24 = 2.038 kg
Time interval between material creation = 60 / (11.58 * 24) = 0.2158 sec.
So, the number of factories you will need for 1 revolution will be 24.In our case, we did
simulation for 2 revolutions, hence 48 factories as shown in the fig.33.
Page | 30
Fig. 33
Title: Description
Fig. 34
15. Select <Total Mass> (Highlighted in the fig.34) and enter its value against it as shown in
the fig. 34.
16. In the <Start Time> cell, enter the time instant at which you want to generate the above
specified mass.
Similarly, for each factory enter the corresponding mass and the time instances.
17. Select the <Environment> section, you will have the following display
Page | 31
Fig. 35
Title: Description
18. Deselect <Auto Update from Geometry> (highlighted in the figure 35)
and enter the domain coordinates manually such that the domain
extends to a significant distance from the exit point of the discharger
system as shown in right.
19. Check the <Gravity> checkbox highlighted in the fig. 35 and enter its
value.
20. Save the changes and click on the EDEM Simulator icon and follow
the under given steps:
Deselect <Auto Time Step> and enter against ‘Fixed Time Step’
22%.
Enter against ‘Total Time’ 10s
Enter against ‘Target Save Interval’ 0.01.
Enter against ‘Cell Size’ 7R.
Click on the Play button on the bar below the viewer window.
Page | 32
Title: Description
Page | 33