Dynamic Modeling of a Vibrating Screen Considering the Ore Inertia and Force of the Ore over the Screen Calculated with Discrete Element Method

Vibrating screens are critical machines used for size classification in mineral processing. )eir proper operation, including accurate vibration movement and slope angle, can provide the benefits of energy savings and cost reductions in the screening process and the whole mining process. Dynamic models of the vibrating screen movement available in the literature do not simulate ore motion or its interaction with screen decks. )e discrete element method (DEM) allows for the calculation of the dynamic of the ore. In this paper, two 2D three-degrees-of-freedom dynamic models for a vibrating screen are tested, using linear and nonlinear approaches for angular displacement. )ese models consider the inertia of the ore and the ore force calculated with DEM. A double-deck linear motion vibrating screen is simulated using the DEM software LIGGGTHS. DEM is used to obtain the ore parameters in the steady state and the force on the screen decks. Two cases are compared: Case 1 considers the ore as moving together with the vibrating screen, and Case 2 considers the ore force on the screen deck as calculated by DEM. Simulations are carried out with data for an industrial vibrating screen used in copper mining. )e force over the screen is significantly different between the cases. Case 1 produces a force that is unrealistic because the ore cannot produce a high-amplitude adhesion force over the screen decks. In Case 2, no adhesion force acts between the ore and deck. It is concluded that the linear dynamic model used in Case 2 is adequate to evaluate the influence of the ore on the movement of the vibrating screen. )e linear dynamic model considering the force as in Case 1 can be used to simulate a vibrating screen, as long as a correct calibration parameter is included to obtain an accurate motion amplitude.


Introduction
Size classification of particulate materials is an important process in mineral processing [1], particularly in the copper industry.Vibrating screens are frequently used to separate granulated ore materials based on particle size for particles with diameters greater than 0.5 mm [2].In the copper industry, the most commonly used vibrating screens have linear motion and horizontal, sloped, or multisloped (banana) screens [3].A vibrating screen has one or more screening surfaces (decks) with square or rectangular openings, and its vibratory movement depends on a system of unbalanced masses.e appropriate type of vibration allows for better movement of ore and stratification of minerals [4].
Vibrating screens are critical machines prone to successive failures that can result in huge economic losses [5] and must be constantly improved in order to meet the requirements of the mining industry [6].Proper operation of this machine has important benefits for the whole mineral process [7], and for this reason, many studies have been conducted to investigate the behavior and operational parameters for successful operation.For example, a change in the slope of 1 °can cause a decrease of 1-2% in the screening efficiency [1,8,9], whereas a change of 1 mm in the vibration amplitude can cause a loss of between 5% and 10% of the screening efficiency [1,9,10].ese efficiency losses depend on the vibrating screen design, ore characteristics, and operational conditions.Inadequate classification produces undersize particles in an oversize stream (overflow) entering a comminution process, resulting in extra energy expenditure and obstruction of the grinding by packing of fines [11].Inadequate classification in a closed circuit produces a short circuit of the undersize stream and results in recirculating a constant percentage of fine ore in the circuit.To ensure proper operation and high performance, it is necessary to have a dynamic model that is able to predict movement of the deck to maintain or increase throughput.
e literature presents several models and studies covering different aspects of vibrating screens.Models of screens using a probabilistic approach [12,13], stratification and passage in the screening process [14], particle movement [15][16][17][18][19], screen blocking [20], finite elements [5,[21][22][23][24], crushing plants [7,25], and phenomenological models [26] can be found in the literature.Dynamic models of vibrating screens can simulate motion of the vibrating screen structure and show good agreement with experimental measurements [27] and finite element method (FEM) results [22].e linear model proposed by He and Liu [28] considers three degrees of freedom.e excitation force is circular, and the vibrating screen structure is supported by symmetrical damping springs with equal stiffness.Liu et al. [29] developed a linear model with three degrees of freedom in 3D that has an excitation force in the vertical direction and different damping spring stiffnesses in each support position.eir study focused on vibrating screen fault diagnosis by performing a dynamical simulation on when the supports lose stiffness.Liu et al. [30] proposed a linear model similar to that developed by He and Liu [28] by incorporating a quadruple excitation mechanism into the model.e model proposed by Baragetti and Villa [22] considers motion in a plane with three degrees of freedom.
eir model was used to calculate the dynamical response and natural frequencies, and the optimization used FEM and experimental measurements.Using that model, Baragetti [31] patented a new design for a vibrating screen.For optimum screening performance, the angular displacement was set equal to zero and held constant over time, and the load eccentricity was also made null.us, the resulting equations have two degrees of freedom.Slepyan and Slepyan [32] and Zahedi and Babitsky [33] analyzed a vibrating screen operating with parametric resonance, using a linear model with damping and a tensile force [32].ey employed a nonlinear dynamic model with a system of autoresonant control [33].Peng et al. [34] developed a model with a single degree of freedom for a large vibrating screen that considered the bending and random vibration in the design of these machines, and in which the ore was simulated as a random force.Jiang et al. [27] proposed a linear dynamic model of a single-deck equal-thickness vibrating screen.In their results, the simulated amplitudes of the vertical and horizontal motion deviated by less than 5% from the corresponding experimental results.Wang et al. [19] developed a nonlinear dynamic model of a planar reciprocating vibrating screen and employed a matrix method to derive the equation of motion in order to analyze the motion of a particle on the screen.With dynamic simulation software, Jiang et al.
proposed a new design of a vibrating screen, where its screen decks are composed by rigid-flexible rods [35].
ese linear models [27][28][29][30] assume that the angular motion of the vibrating screen is low and implement linearization as sin θ ≈ θ and cos θ ≈ 1, which is useful under nominal operating conditions where angular displacement is not significant.However, in practice, significant angular motion occurs in vibrating screens during startup [31], shutdown, and under off-design operating conditions [3].
e nonlinear models [19,33] have been developed for particular types of vibrating screens and thus are not always applicable to the vibrating screens normally used in the copper industry, due to its deck material or type of movement.
Rodríguez et al. [3] developed a 2D nonlinear dynamic model of a vibrating screen with three degrees of freedom that allows for significant angular motion and damping in which the nonlinearity is geometric due to angular displacement.
ey proposed a range of admissible loss of stiffness as a percentage of nominal stiffness in order to ensure proper operation using orbital analysis.e calculated stiffness range was found to be 38% on the feed side and 46% on the discharge side.Moncada and Rodríguez [36] used a nonlinear model [3] to calculate the effects of loss of stiffness in the supporting positions on the steady and transient responses of a vibrating screen.For the steady response, the change in orbital direction was analyzed, while the transient response analyzed the change in the natural frequencies due to loss of stiffness in the supporting positions.
Simulations carried out with full-ore loads [3,19,36] assume that the ore has the same movement as the vibrating screen, i.e., the ore is added to the decks, approaching this force to the ore inertial force −M ore a cm , where M ore is the mass of the ore and a cm is the acceleration of the center of mass of the ore. is produces an unrealistic attractive force when the ore is in its highest position (ore cannot pull the vibrating screen deck) or in the free-fall phase or when throwing index is equal to zero [15,17].
Because the force generated by the ore material over the decks is inaccurate owing to the impossibility of traction or a pulling force over the deck, this study aims to calculate the ore force over a vibrating screen deck by means of simulations using the discrete element method (DEM). is method considers the interaction of each particle with each other and with decks and allows forces for different simulation conditions to be obtained.Furthermore, this method is widely used in the literature.In the mining field, several machines and processes have been simulated using DEM, including cone crushers [37], mills [38], hopper discharge [39], jaw crushers [40], feed boxes [41], and 2 Shock and Vibration vibrating screens [1,[8][9][10][42][43][44][45][46][47][48][49][50][51].Cleary et al. [42,43] performed a DEM analysis of an industrial double-deck banana screen for a range of peak accelerations and two feed size distributions.Dong et al. [44] conducted a numerical analysis of the particle flow on a banana screen and demonstrated the importance of operational parameters like the slope angle of each deck, vibration amplitude, and frequency.Zhao et al. [9] carried out a numerical study of the motion particulates follow along a circularly vibrating screen deck using 3D-DEM.ey studied the effects of vibration amplitude, throwing index, and screen deck inclination angle on the screening process.Fernandez et al. [45] used a one-way coupled model of smoothed particle hydrodynamics (SPH) and DEM to simulate large banana screens.DEM was used for the coarse particulate flow, while SPH was employed to model the transport of the fine particle slurry.Delaney et al. [46] used DEM to investigate the flow of a granular material over a horizontal vibrating screen, and performed a quantitative comparison between laboratory scale experiments and the simulation results.Liu et al. [47] simulated the particle flow on a banana screen deck using DEM and investigated the effects of slope angle and deck length on the screening process.Dong et al. [48] simulated the screening of particles for different vibration modes (linear, circular, and elliptical) and studied the resulting motion and penetration of the particles on the screen deck.Li et al. [49] used DEM to optimize the design and operational parameters of a linear vibrating screen.Jahani et al. [10] investigated the screening performance of banana screens using the DEM solver LIGGGHTS.An industrial double-deck banana screen with five panels and two laboratory single-deck banana screens with three and five panels were simulated, and the effects of design parameters such as the slope angle of decks, vibration amplitude, and frequency were analyzed.eir results were validated with partition numbers obtained from the literature [44].Zhao et al. [50] quantitatively compared DEM results and experimental data for a specially designed circularly vibrating screening model under a range of operating conditions.Jafari and Nezhad [8] studied the effects of different parameters on process efficiency and mesh wear using LIGGGHTS.Dong et al. [51] conducted a numerical analysis of the effects of aperture shape, length, and orientation on particle flow and separation in a vibrating screen process.Zhao et al. [1] analyzed the combined effects of vibration parameters on a circularly screening processes.With a DEM simulation, the effects of various design and operating variables on the efficiency of screen were investigated using open-source LIGGGHTS solver with spheral and irregularly shaped particles [52].Particle velocity, mass of oversized material, screening efficiency, and impact force were studied to reflect the performance of vibrating screen [53].
From this review of DEM simulations, we can summarize the following: (i) there is significant evidence that vibratory parameters are important for proper screening performance [1, 8-10, 42-44, 47, 52, 53]; (ii) DEM simulations do not consider the effects of ore mass flow on amplitude movement although mass flow affects amplitude [54,55] and efficiency depends on amplitude [1,9,10,42,43,52,53]; and (iii) DEM simulations do not focus on vibration.Most models use circular motion even when there are experimental data showing that the motion is not circular.In addition, rotation is neglected [19], which significantly affects the movement described for the vibrating screen supports.
e dynamic simulations in this study make two different assumptions or considerations.In Case 1 [3,36], the ore moves together with the screen deck, which means that physically the ore adheres to the screen and there is no relative displacement between them.For modeling, parameters such as the mass of the vibrating screen with the load and the position of center of mass must be calculated.In Case 2, interaction between the ore and the screen deck is represented by a time vector force for each degree of freedom (x, y, and, θ).Movement of the ore is simulated with DEM to obtain the force over the screen deck, and this force is then used in the calculation of the dynamic model.
In this study, a 2D linear model with three degrees of freedom that considers ore inertia and the ore force over the screen calculated using discrete element method is developed in order to determine the influence of the ore on the movement of a vibrating screen.For cases 1 and 2, DEM is used to simulate a double-deck vibrating screen.
e DEM simulation is carried out in LIGGGHTS and is used to obtain ore parameters and the ore force on the screen over time.Finally, the two dynamical responses are compared.

Dynamic Model
e dynamic model for a linear motion vibrating screen has three degrees of freedom in the plane xy, as shown in Figure 1(a) in its equilibrium position.e excitation force of the system is time-dependent on amplitude, F sin(Ωt), has a constant direction, α, with respect to the line linking points A and B, is applied at a distance, p, from the center of mass, and has Ω excitation frequency.e structure of the vibrating screen and all its components is defined as a rigid body inclined at an angle, β, with a center of mass located at a distance, h, perpendicular to segment AB. Supports are located at points A and B and consist of a horizontal and vertical spring and damper, each at a distance, a and b, from a projection of the center of mass, and with a stiffness, K, and damping, C. For the simulation cases in this study, two approaches for modeling the ore are used.In Case 1, the ore is modeled as a rigid body and vibrate together with the screen deck.
erefore, the mass, M, inertia, I, and position of the center of mass of the vibrating screen with a load are the sum of the empty vibrating screen plus the ore.In the second case, the ore movement is calculated by DEM simulation.e load is represented in the dynamic model by a time vector force with three components for each degree of freedom: F x,dem (t), F y,dem (t), and F θ,dem (t), applied at the center of mass of the vibrating screen.

Parameters for Case 1.
Parameters for the ore particles in the DEM simulation must be known in order to perform the dynamical simulations in Case 1.For that purpose, the mass, m i , radius, r i , inertia, I i , and position (x i , y i ) of every particle in the N particles system are used in equations (1a)-(1f ) to obtain the load parameters.
ese load parameters include the total mass, M ore , inertia, I ore , position of the ore center of mass (x cm,ore , y cm,ore ) relative to the axes (x dem , y dem ), and position relative to the vibrating screen (a ore and h ore ). Figure 2 shows a vibrating screen modeled as a rigid body, along with the upper and lower deck, feeder chute, particles, and variables used in these equations.It also illustrates the geometry to allow comparison of the dynamic modeling and DEM, and it should be noted that the variable x in the dynamic modeling and x dem in the DEM modeling have different directions.e subscript cm represents the center of mass, ore represents the ore load over the screen decks, em represents the empty vibrating screen, msh represents the screen deck, and low and up represent the lower and upper decks, respectively.
With these load data, the parameters for the vibrating screen with a load can be obtained using equations (2a)-(2g).Figure 2 also shows this situation, illustrating the position of vibrating screen center of mass without a load or empty, em, with a ore load, ore, and the new position of the center of mass of the vibrating screen with a load, represented without subscripts.
e average values of M, I, h, and a in steady state are used in the model.

Forces in Case 2.
e force calculated with DEM, F dem , is the force of all the particles over the screen decks.Applying Newton's second law to the N particles system results in the following: where F dem ′ is the force of the screen decks over the particles, m i is the mass of a particle, g � −9.81 m/s 2 j is the gravitational acceleration, and a i is the acceleration of each particle.Shock and Vibration erefore, using Newton's third law, the force of the particles over the screen decks is As F dem includes the inertial force of the particles, Case 2 also represents the dynamic characteristics of the particles.
is force is applied to the vibration screen center of mass and parameters for the vibrating screen in empty conditions are used.

Equations of Motion.
e equations of motion are developed using Lagrange's equations.Applying the linearization sin θ θ and cos θ 1, the resulting linear equations of motion are given in equations ( 5)- (7).
ese three equations correspond to a second-order linear system of equations, and can be solved using the Newmark method with a constant time step.e nonlinear model used for simulations where the angular displacement is signi cant is presented in Appendix [3].In these equations, nonlinearity is present in the terms _ θ sin θ, _ x cos θ, and cos 2θ, among others.For the force direction, nonlinearity is present in terms that consider a change in θ, as can be seen in the term cos(θ e equations of motion are developed in an equilibrium position with an inclination angle, β.If the position of the center of mass or amount of mass changes, the vibrating screen will experience a variation in its equilibrium position, the slope β.A free-body diagram calculation allows the new slope of β to be determined before these changes occur.is is useful when a loss of sti ness exists in the support system [36].

DEM
e discrete element method was originally developed by Cundall and Strack in 1979 [56] and it has proven to be an e ective numerical technique for calculating particle movement in granular material ows.It is also useful for obtaining the force acting on each particle, which is di cult to obtain on a similar scale with physical experimentation [57].Since it is a Lagrangian method, Newton's second law is applied for each particle in the simulation domain.is explicitly determines the trajectories and kinematics of each particle at every time step by accounting for the interaction between particles and their environment with contact or eld forces.Using Newton's second law, the equations of motion for the translational and rotational motion of each particle, i, in contact with a particle or wall, j, are as follows: (3) where m i is the particle mass, v i is its velocity, g is the gravitational acceleration, F ij is the contact force exerted by particle j on i, I i is the particle inertia, ω i is its angular velocity, and T ij is the contact moment.By assuming that the particles are rigid, spherical, and of uniform material, the contact force, F ij , for particleparticle and particle-wall interactions can be calculated with the Hertz-Mindlin contact model [56]. is force is composed of an elastic and viscous component and is expressed by equations (9a)-(9c) in the normal and tangential directions, being  e n and  e t the respective unit vectors.
e tangential overlap is truncated to fulfil Coulomb friction criterion in tangential force.Overlap distance δ ij is composed of a normal and tangential component, as shown in equation ( 9d).
e choice of a contact model depends primarily on the experimental comparison and type of results expected.A simple model, such as the Hooke linear model, can also produce good agreement in some situation [58].
A constant directional torque (CDT) model is used to calculate the rolling resistance [59,60].A constant moment, T r,i , is applied to a particle to represent the rolling friction, where ω i and ω j are the angular velocities of the particles in contact, μ r is the rolling friction coefficient, r ij is the effective radius, and F n,ij is the contact normal force.e torque is always in the opposite direction to the relative rotation between the two contact particles.
Figure 3 shows the schematic of the contact model and rolling resistance model, illustrating two particles, i and j, with an overlap, δ ij , and with contact stiffness, K, and damping, C. Each particle is characterized by its Young's modulus, E i , mass, m i , and radius, r i .
Model coefficients for equation ( 9) are calculated using equations given in the LIGGGHTS manual [61].ese coefficients depend on the material parameters, such as the restitution coefficient, e, Young's modulus, E, and Poisson coefficient, v.
To numerically solve the motion equations for each particle, the simulation is divided into time steps, Δt, and the equations are solved to obtain a solution at the end of each time step.e velocity Verlet algorithm is used for these calculations, which can be derived by an approximation of the Taylor series [62].
DEM simulations are implemented in LIGGGHTS [63], an open source software based on LAMMPS used to perform massive granular simulation in parallel using MPI [64]. is software makes it possible to import CAD files for the wall geometries; assign different movements to these geometries, such as a 3D vibration with a different amplitude, frequency, and phase; use nonspherical particles [65]; perform stress and wear analysis [8]; and develop smoothed particle hydrodynamics (SHP) models [66] and CFD-DEM coupling simulations [67].For these reasons, it has been used in the literature for vibrating screen DEM simulations [8,10].LIGGGHTS 3.6.0was used to perform the simulations on a desktop computer.

Vibrating Screen Simulation
e simulations in this study make use of data obtained from a vibrating screen in use at a copper mine [3]. is screen is a double-deck linear motion vibrating screen with a size of 3.66 × 7.32 m and a 10.5 mm nominal stroke.
e CAD model of this vibrating screen is shown in Figure 4(a), where the upper and lower screen decks, feed chute, and lateral and rear walls that serve as boundaries of the particles are detailed.e lateral wall on the left side is not presented in this figure, but it is included in the simulation.Only the features necessary for DEM are geometrically modeled, as the details of the structure of the machine are not relevant; the screen decks are correctly modeled.e screen deck is composed of 27 square modules in the longitudinal direction, with a side length of each square module of 305 mm.Each deck has rectangular openings, slotted with the flow direction, of 60 × 20 mm (Figure 4(b)) and 47 × 11 mm (Figure 4(c)) for the upper and lower decks, respectively.e opening area represents 30% and 36 of the upper and lower decks, respectively.
To reduce the computational cost of the DEM simulation, the simulated vibrating screen is a 1/12th scale model in the cross direction [44] of the real machine, using the symmetry in the xy-plane corresponding to a row of modules, as shown in Figure 4. is simplification implies that the intensive variables, Shock and Vibration such as force, mass, and inertia, must be multiplied by 12.
In spite of the influence of the particle shape on stratification and passing [46], spherical particles are adopted because of their lower computational cost.e particle size distribution corresponds to nominal data provided by the manufacturer, and a diameter threshold of 6 mm is used in order to reduce the simulation time [68].Young's modulus is E p � 2.3 • 10 7 Pa for the particles, and E wall � 5 • 10 5 Pa for the boundary walls.e common simplification to reduce E [10] is not applied because it influences the contact force [69], and thus, the force between the particles and decks.With the period of vibration, T, and N � 2 15  steps per cycle, a time step is Δt � T/N, or approximately 2 • 10 −6 s. at time step is acceptable for this simulation [63] because it is smaller than the Rayleigh time step.Details of the geometrical conditions and material parameters are listed in Table 1 and are based on nominal data and previous studies [17,37,46].ese parameters correspond to the simplified vibrating screen, and thus, the mass flow and width are 1/12th the real value.e subscript p refers to the particle.
To promote stratification and provide movement to the ore, a vibrating screen vibrates with a particular frequency, direction, and amplitude.In the DEM simulations, motion is imposed on both decks, with vibration in the x-axis and y-axis, as well as angularly an amount θ with respect to its center of mass.e vibration direction is ensured by the phase difference between the vibrations in x and y. e vibrational parameters are listed in Table 2.
ese parameters are obtained from dynamical simulations, and thus the frequency, f, of the deck vibration movement is the same as that in the dynamic simulations.To the comparison of the cases be valid, both cases consider the same value of M ore .
Table 3 lists the parameters used in the dynamical simulations for Cases 1 and 2 that correspond to the nominal data without load and geometrical conditions.e excitation frequency is calculated by Ω � 2πf.
For Case 1, the DEM data, parameters for the empty vibrating screen, and equations ( 1) and ( 2) are used to obtain the parameters listed in Table 4. [70].

DEM. DEM simulation results are available in
Figure 5 presents snapshots of the DEM simulation of the vibrating screen in the steady state.e ore entering through the chute, distribution on the decks, stratification, and passing can be observed.In this figure, the red spheres have a larger diameter and the blue spheres have smaller diameters.e steady state is defined as when the ore mass on the vibrating screen remains approximately constant with time.Shock and Vibration e total force and moment over the screen exerted by the ore are presented in Figure 6 for the x, y, and, θ directions.e vertical component, F y,dem , has a greater amplitude than the horizontal component F x,dem , and its peak amplitude is 7 times greater.Both amplitudes match when they have null values, as that corresponds to when the ore is detached from the screen.
e value of F y,dem exhibits an amplitude change every cycle, i.e., between peaks (1) and ( 2), that is related to the ore movement.On average, the first is equal to 122.4 kN and the second is equal to 111.5 kN. e moment, F θ,dem , has the same tendency as components x and y, and its maximum amplitude occurs when ore comes into contact with the screen deck.
e overall screening performance of both decks can be investigated with the partition curve of the overflow, which is defined as the ratio of the number of residue particles in the overflow to that of fed particles [44], and is shown in Figure 7.In terms of mass flow, the partition number is equal to for each i particle size, where o refers to overflow stream and f to feed.ese results agree with the nominal screening efficiency of the vibrating screen, which is approximately 90%.

Comparison of Case 1 and Case 2.
e results for Case 1 are compared with those of Case 2, and the force exerted by the ore on the vibrating screens and movement of the machine in its supports are analyzed.

Inertia
Force and DEM Force. Figure 8 compares the forces in Case 1 and Case 2 with a frequency spectrum and orbit.Figure 8(a) shows the frequency spectrum of the signal in the vertical direction, y.For Case 1, only the 1X component is observed.For Case 2, a constant component at 0 Hz is observed, which represents the mean time value, along with harmonics from 1X-4X.e logarithmic spectrum shown in Figure 8(c) exhibits, in addition to these harmonics, nonsynchronous components at 0.5X and 1.5X, which are highlighted in the grass spectrum.
Figure 8(b) shows both forces in the xy-plane, allowing for a graphical understanding of the direction of each force.In Case 1, the orbit is elliptical with a nearly negligible semiminor axis, and thus, it approximately corresponds to a line.Case 2 has an ore force that is always negative in the direction of y because the ore cannot produce an adhesion force that lifts up the screen decks.In analyzing these results, it can be concluded that Case 1 does not accurately physically represent the effect of the ore on the vibrating screen, because this force must be vertical and repulsive, as in the DEM simulation.

Supports Movement.
Movement in the two cases is compared by means of orbit analysis of the movement of support A. Figure 9 presents the frequency spectrum and displacement orbit for support A. Figure 9(a) shows the frequency spectrum of vertical movement and indicates that Case 1 only exhibits a 1X component, which differs by 12% from that of Case 2. Case 2 exhibits a low amplitude harmonic at 2X.For Case 2, owing to the transient condition of F dem , a resonant zone appears near 1.875 Hz.
is corresponds closely to the values of the natural frequencies, ω nt,θ and ω nt,x , which are 1.982 Hz and 2.067 Hz, respectively.
Figure 9(b) shows the orbits.In Case 1, an elliptical orbit is clearly observed, while in Case 2, the transient condition of the excitation force, F dem , results in an elliptical trajectory that changes position in the vertical direction by an amount Δy � 1.2 mm.A difference in stroke 2λ i is also observed, equal to 2λ 1 � 9.230 mm in Case 1 and 2λ 2 � 10.592 mm in case 2 on average.ese values agree with commonly found in experimental measurements and in the literature [3,10].Furthermore, Figure 9(b) presents the nominal stroke length equal to 2λ nom � 10.5 mm.Comparing nominal stroke with simulated stroke, case 2 provides a better agreement.

Angular Displacement.
Angular displacement of the vibrating screen is obtained, and it is presented in Figure 10.Both cases present a low peak amplitude, less than 0.025 °.
Figure 10(a) presents frequency spectrum showing that both have a 1X component that differs on 0.003 °.Case 2 presents a low amplitude 2X component and, as well as vertical movement, a resonant zone.Figure 10(b) shows waveform of angular displacement.It should be noted that Case 2 has nonzero mean value equal to θ 2 .
is value depends of F θ,dem , which is calculated based on F x,dem and F y,dem and the position of the center of mass of the ore with respect to the center of mass of vibrating screen.

Discussion and Conclusion
Conclusions obtained in this study can be summarized as follows: (i) e proposed dynamic model allows for prediction of the behavior of the vibrating screen operating at both the nominal condition and high angular displacement with a nonlinear model [3]. is model simulated ore movement along with the screen deck

Shock and Vibration
(Case 1), as well as the ore force over the screen decks calculated with DEM (Case 2).In both cases, the movement of the vibrating screen supports was also obtained.(ii) DEM simulation of a double-deck vibrating screen was carried out using the open-source LIGGGHTS software.Movement of the ore center of mass and the force exerted by the ore over the screen deck were obtained.e F dem force has the same frequency as the excitation force; high oscillations in the ore do not produce significant changes in the force exerted over the screen.(iii) e partition curve and stroke of the vibrating screen motion have very good agreement with the nominal data, validating the model results.(iv) In a comparison of the results for the proposed cases, where the ore is represented as moving together with the vibrating screen (Case 1), or the ore force is obtained from DEM (Case 2), it is observed that: e force over the screen deck is completely different in both cases, both in terms of the magnitude (the peak-to-peak amplitude in Case 1 is more than twice that of Case 2) and the shape on the xy-plane.Case 1 produces an unrealistic force, because it includes a contact force of adhesion between the ore and screen deck.However, this is the hypothesis used in most of the existing dynamic models available in the literature.
e DEM model allows calculation of a force closer to reality, because it calculates the interaction of each individual particle with the screen deck.(b) Case 1 available in the literature was successfully implemented.e computation of Case 1 is faster than Case 2, because only a linear ordinary differential of three-degrees-of-freedom equation is solved, in contrast to DEM simulation that simulates the movement of 42000 particles.is is the reason why Case 1 is commonly used in the mining industry.(c) Case 2 is a new simulation approach that allows the coupling of simulation results in DEM in dynamic models.is model is able to evaluate how the ore affects the movement of the vibrating screen.(d) Notwithstanding the clear inequality in the force calculated for cases 1 and 2, the approach used in Case 1 can be used to predict the movement of vibrating screen if a correct calibration parameter is included.(e) Comparing with the nominal data and the results Case 2, frequency, direction, and inclination calculated with Case 1 are accurate.e amplitude obtained using Case 1 is not accurate and must be corrected.In this case, the parameter of mass of ore M ore must be adjusted with experimental data, decreasing its value so that the amplitude increases.To decrease the value of M ore is physically correct, since the model of Case 1 considers that all the ore is in contact with the screen decks, while there is ore in free fall that is not in contact.
(v) Under nominal operating conditions, the angular response, θ(t), has low amplitude (0.014 °), whereas the steady responses obtained for the linear model, equations ( 5)- (7), and the nonlinear model, equations ( 12)-( 14), have negligible differences.Consequently, for simplification and lower computational cost, the linear model can be used for the steady case or without deterioration in the supports.For transient signals or when the amplitude of the angular response is high, the linear model is not recommended.(vi) e proposed dynamic model allows for greater accuracy and validity in different operating conditions, which is useful for predicting the angle of operation and the vibratory amplitude, parameters that affect the screening efficiency.

Data Availability
e DEM raw data used to support the findings of this study have been deposited in the "Raw data of DEM simulation of linear motion double deck vibrating screen" repository (DOI: 10.17632/cc798dvdbn.2).

Figure 1 (
b) presents the degrees of freedom of the model-horizontal displacement, x, vertical displacement, y, and angular displacement, θ.

MFigure 2 :
Figure 2: Geometric comparison of dynamic modeling and DEM, showing the change in the center of mass due to the load, with (1) upper screen deck, (2) rigid-body model, (3) lower screen deck, and (4) feed chute.

Figure 7 :FFigure 8 :FrequencyFigure 9 :
Figure 7: Partition curves of the overflow for the upper and lower decks.

Figure 10 :
Figure 10: Comparison of the angular displacement of vibrating screen: (a) frequency spectrum and (b) waveform.

Table 1 :
Parameters used for the DEM simulations.

Table 2 :
Vibrational parameters used for DEM simulation.

Table 3 :
Parameters of empty vibrating screen used for the dynamical simulations.