SciELO - Scientific Electronic Library Online

 
vol.22 número3Funciones agua rendimiento para 14 cultivos agrícolas en condiciones del sur de La HabanaIdentificación de aspectos ambientales en la comunidad Loma de Cruz, Santiago de Cuba índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Journal

Artigo

Indicadores

  • Não possue artigos citadosCitado por SciELO

Links relacionados

  • Não possue artigos similaresSimilares em SciELO

Compartilhar


Revista Ciencias Técnicas Agropecuarias

versão On-line ISSN 2071-0054

Rev Cie Téc Agr vol.22 no.3 San José de las Lajas set. 2013

 

ARTÍCULO ORIGINAL

 

Formulation of a frictional-cohesive soil model by the Discrete Element method

 

Formulación de un modelo friccionante-cohesivo de suelo por el método de elementos discretos

 

 

Dr. C. Elvis López Bravo I, Dr. C. Miguel Herrera Suárez I, Dr. C. Omar Gonzalez Cueto I, PhD. Engelbert Tijskens II, PhD. Herman RamonII

I Universidad Central de Las Villas. Departamento Ing. Agrícola, Villa Clara, Cuba.
II Universidad de Leuven,  Departamento Bioestadistica y Sensores, Bélgica.

 

 


ABSTRACT

Friction and cohesion in efforts that take place during soil deformation calculations. Para interactions are considered normal force, tangential cohesive friction. The maximum shear is governed by the criterion of Mohr-Coulomb failure. Model calibration is performed by obtaining the relationship between micro and macroscopic parameters used in the model. The calibration of the friction was carried out by simulating the direct shear test of a sample of soil micro different friction coefficients. To calibrate cohesion used different values ​​of cohesion at the micro and the relationship was obtained by modifying the micro friction values ​​by simulating the biaxial compression test. The friction results showed a linear relationship regarding the variation of the micro value. Instead cohesion was linearly affected by the change inthe quadratic micro-cohesion and friction over micro.

Keywords: simulation, tillage, modeling, discrete element method, draft force.


RESUMEN

En el presente trabajo se propone un modelo de elementos discretos  que incluye el efecto de la fricción y la cohesión en los esfuerzos que tienen lugar durante la deformación del suelo. Para los cálculos se consideran las interacciones entre la fuerza normal, tangencial, cohesiva y friccionante. El valor máximo de fuerza cortante es regulado mediante el criterio de falla de Mohr-Coulomb. La calibración del modelo se realizó mediante la obtención de la relación existente entre los parámetros micro y macroscópicos usados en el modelo. La calibración de la fricción se realizó por medio de la simulación del ensayo de corte directo de una muestra de suelo a diferentes valores de fricción micro. Para calibrar la cohesión se utilizaron diferentes valores de cohesión a escala micro y se obtuvo la relación existente al modificar los valores de fricción micro mediante la simulación del ensayo de compresión biaxial. Los resultados de la fricción mostraron una relación lineal respecto a la variación de su valor micro. En cambio la cohesión estuvo afectada de forma lineal por la variación en la micro-cohesión y cuadrática respecto a la fricción micro.

Palabras clave: simulación, labranza, modelos, método de elementos discretos, fuerza de tiro.


 

 

INTRODUCTION


Soil in a discontinuous model is commonly represented by assemble of spherical particles with bonds that can break up and reform under external loading. Developed by (Cundall, 1971), the method was developed for simulated the rock slope stability and soil structures. When suitable contact force models are used and microscopic model parameters are carefully selected, the DEM procedure can offer a realistic description of soil mechanical behavior.

However, as a real soil is composed of a huge number of small particles with complex shapes, it is inevitable to represent such a granular material as an assembly of larger idealized particles. Micro-macro parameter relationships resolve the difference in particles sizes and shapes, but introduce more complexity into the method (Asaf et al., 2007).

Simulations of soil mechanical behaviour generally showed good qualitative agreement with experimental tests, but accurate quantitative results were not yet obtained (Shmulevich, 2010). The differences were related to the differences between the size and shape of the soil particles with those in the model.

Considering the characteristics of the clay soil obtained by (Lopez et al., 2012), especially its considerable cohesion, the present study is focused on the virtual tests to calibrate the model.To obtain the relationship between soil micro-macro friction and cohesion the simulation of direct shear and biaxial compression test were performed in the standard configuration.

 

METHODS

Model formulation

The soil-tillage tool interaction dynamics were implemented in the DEMeter++ Software (Tijskens et al., 2003). This program provides tools for building and computing the interaction between simple and complex bodies, including contact detections, relative bodies’ velocities, contact force calculation, time integration, etc. Written in C++ this software has the possibility to utilise open source libraries, compile own applications and create executable files. The output of the calculations can be saved according to the kind of application providing easy and comfortable management of the entire system. The Mohr-Coulomb criterion of shear failure is used to regulate the force in tangential direction: equation 1

where:

Ftma = tangential force limit [N],

Fn= normal force [N],

fm = micro-friction angle [°],

cm = micro-cohesion [N].

The variation in normal force DFn at the contact point has elastic and a viscous component:equation 2

where:

kn= normal stiffness [N/m],

un = normal displacement [mm],

hn = viscous damping coefficient [kg/s],

Dt = time step [s].

The stiffness parameter is calculated by the model proposed by Liao (1997) based on the stress-strain relationship. To obtainthe values of the dimensionless parameters (ak=2.65, = 0.65, =1.0) in the Liao model, a set of compression tests was carried out by Hentz et al.(2004). The model was also verified in a mathematical study by Feng et al.(2007)with the purpose of validating the DEM prediction capacity. From this modelthe stiffness in normal direction was determined as:equation 3

where:

Eab= equivalent Young’s modulus of the materials in contact [Pa];

Aint= interaction surface [m2]:

Deq= equivalent distance between the two particles [m];

n= Poison’s  ratio;

ak,= loading path;

bk, = softening factor;

g =interaction range.

Similar to the normal force, the variation of tangential force DFs was calculated considering the elastic and plastic effects:equation 4

where:

ks= tangential stiffness [N/m];

Dus= increment of tangential overlapping [m];

hs= viscous damping coefficient in tangential direction [kg/s].

The magnitude of the tangential stiffness also depends on the value of the normal one Hentz et al., (2004) and is calculated by the following equation:equation 5

The viscous damping is obtained as: equation 6

where:

i = right subscript representing normal and tangential direction;

b  = coefficient of viscous damping;

ma, mb = mass of the objects in contact [kg].

Finally, the equations to calculate the model parameters of micro-friction fm and micro-cohesion c are obtained by modeling the soil macro-behaviour during the simulation of biaxial compression test and direct shear test.

Geometrical model of virtual shear test

In order to reproduce the direct shear test with a DEM model a virtual reproduction of the direct shear test apparatus (Figure 1) was made. A set of samples for the simulation were also generated by spherical particles to test the micro-friction fm effect over macro samples strength. The geometrical setup for the simulation (Figure 1a) was divided in two main sections: upper and lower box; each one with a cylindrical cavity at the centre.

The soil specimen is obtained by sequential steps, starting with spheres generated inside the cylindrical cavity of the shear apparatus following a hexagonal compacted spatial distribution. The centre of each particle is obtained by assigning an independent position in the (x, y, z) axes for an initial number of 4,000 particles. Dimensions of the particles were calculated by random distribution of the radii between 1.0 to 1.5 mm. In the next step the particles falling under the action of the gravity force.

The particles number is fixed at 3,680 to obtain the desire height of the sample. The generated soil sample has a cylindrical shape of 70 mm diameter by 60 mm height. From this point, the geometrical model is ready for simulation of direct shear test. The model macro-parameters (Table 1) were calculated by means of the regression equations defined by Lopez et al. (2012) at the soil intermediate physical condition according to soil moisture w = 21.3% and dry bulk densityd= 1.18 g/cm3.

During the simulation the sample is cutting in the transversal direction as shown in Figure 1b.The shear operation takes place moving the upper box over a length of 20 mm at a constant velocity of 1.3 mm/s. The press surface ensures that shear take places at different normal load.

The relationship between micro and macro friction was obtained by changing the value of the micro-friction angle in several tests computed with the same soil parameters defined on the Table 1. Accordingly, four values of fm were tested (2, 4.2, 6.4 and 8.6°), each one at a normal pressure of 30, 50, 70 and 90 kPa resulting in a total of sixteen simulation tests.

Geometrical model of biaxial compression test

Soil by nature is a frictional and cohesive material; its strength is a result of both contributions. This behaviour is introduced in the model by the micro-friction and micro-cohesion parameters acting at the particle level.

By means of the simulation of biaxial compression test the relationship between real soil cohesion and micro-cohesion for the model is obtained. In order to find the variation of soil cohesion the simulation is carried out in unconfined conditions, which means without lateral pressure.

With the compression-displacement curve resulting from the test, the cohesion calibration procedure is made following the Mohr-Coulomb criterion of failure.

The procedure to form the virtual samples was the same as the one used for direct shear test, i.e. coordinates generation, particles falling, compression, and decompression. An auxiliary cylinder of 50 mm diameter is used as a mould (Figure 2a) for a virtual sample made by 5,500 spherical particles. Before running the compression test the cylinder (mould) is removed from the installation and a total of 5,350 particles form the final sample after the height is fixed. The particles follow a random radius distribution between 1.5 to 2 mm, sized at 50 mm diameter by 100 mm height.

Finally, the geometrical model of the biaxial compression test is composed by the press surface, the base, and the soil sample (Figure 2a). In the experiment the press surface is moved at a constant velocity of 1 mm/min in a downward direction (Figure 2b). The reaction forces are computed until the sample reaches a deformation of 20% in vertical direction.

The micro-macro relationship for the cohesion parameter is obtained by testing the compression virtual model at three different values of micro-friction c= 40, 80 and 120 kPa;

To analyse the effect of particle friction over micro-cohesion  the  biaxial compression test was executed four  times  changing only  the micro-friction angle at fm = (2, 4.2, 6.4 and 8.6°) and keeping constant the inter-particle cohesion.

 

RESULTS AND DISCUSSION

Macro and micro-friction relationship

The shear strength at failure is selected as the maximum tension reached during the test. After this point the sample strength decreases. When the normal pressure increases, a new value of shear strength is attained; the results for each combination are plotted in the shear-normal plane as shown in Figure 3.

Shear strength at failure increases linearly with respect to normal pressure (Figure 4), in agreement with the laboratory experiments. In addition, an increment in the micro-friction parameter enhances the values of shear tension.

The slope of the regression line of shear strengths gives the internal friction angle of the material corresponding with the micro-friction tested in the model. The relationship between macro and micro friction is obtained from the regression line describing by the soil internal friction angle versus the micro-friction parameter (Figure 5).

 

The linear regression equation characterises the particular condition tested, influenced mainly by the particles size distribution and sample porosity. The equation to define the relationship between internal frictionf and micro-frictionfm is written as:ecuacion 7

By this equation the values of inter-particles friction
fm in the DEM model can be calculated as a function of the real soil internal friction angle f.

The model should be able to reproduce the soil shear strength according to the normal pressure applied.

Macro and micro cohesion relationship

Shear strength at failure is localized in the compression-displacement curve as the maximum value achieved during the test. Particle accommodation takes place at the start of the test without substantial increases in forces (Figure 6). After the first millimetre of deformation the sample strength increases until reaches the maximum value. After this point the tension gradually decreases. These results are in agree with those obtained in clay soils by experimental mechanical tests (Herrera et al., 2008; Gonzalez, 2011).

 


 

The relationship between micro and macro cohesion for different incremental values of internal friction is shown in Figure 7.

Increasing soil micro-cohesion results in a linear increment of soil cohesion determined by the slope in the regression equation obtained for each value of micro-friction under test. In order to find an equation to involve both parameters of soil strength, the slope of c-cμ versus the corresponding values of micro-friction are plotted.

As shown in Figure 8, non-linear behaviour distinguishes this relationship. As the slope of c-c increases, the micro-friction angle increases too. The regression equation is written as:

As shown in Figure 8, non-linear behaviour distinguishes this relationship. As the slope of   c-c increases, the micro-friction angle increases too. The regression equation is written as: equation 8

where:

fm = micro-friction [°]:

a =  0.15;

b =  0.18[ ].

As shown in the equation 8, the magnitude of macro-cohesion resulting from the model increases proportionally with the increment of inter-particles cohesion and friction angle, a quadratic relationship was also found to respect of micro-friction angle . This results agree with those obtained by Delenne (2004) in a study of the influence of the friction parameter in a cohesive bonds.

 

CONCLUSIONS

To calculate the mechanical response of a cohesive soil at particle level in the DEM model the system forcescan be formed by normal force, shear force, friction force, cohesion force and gravity force. Shear force at failure is defined by the Mohr-Coulomb criterion.

Using the geometrical elements from DEMeter++ software a virtual reproduction of the direct shear test and compression was made. The strength obtained in the virtual soil samples during the simulation are used to find the relationship between micro and macro parameters.

The shear strength at failure during the simulation of shear test was found to increases linearly with respect to the normal pressure. The same behaviour is shown as the micro-friction parameter increase. A linear regression equation characterises the relationship between micro-friction and macro-friction, influenced by factors associated to particles size distribution and sample porosity.

Alinear increment in soil macro-cohesion was the result of increases micro-cohesion parameter in the model. However a quadratic behaviour was found as the micro-friction increase.

The statistical equations obtained through the simulation allow calculating the micro values for the DEM model started from the macro-parameters of the soil valid for soil models with the same configuration of particle arrangement and force system.

 

 

REFERENCES

1. ASAF, Z.; D. RUBINSTEIN y I. SHMULEVICH: "Determination of discrete element model parameters required for soil tillage", Soil and Tillage Research, 92(1-2): 227-242, 2007.

2. CUNDALL, P. A.: "A computer model for simulating progressive, large-scale movements in blocky rock systems ", Pro. Symp. Int. Soc. Rock Mech, Nancy 2(8): 1971.

3. FENG, C.; E. C. DRUMM y G. GUIOCHON: "Prediction/Verification of Particle Motion in One Dimension with the Discrete-Element Method", International Journal of Geomechanics, 7(5): 344-352, 2007.

4. GONZÁLEZ, O.: Modelación de la Compactación Provocada por el Tráfico de Neumáticos de los Vehículos Agrícolas en Suelos en Condiciones de Laboratorio, 209pp., Tesis (en opción al grado científico de Doctor en Ciencias Técnicas Agropecuarias), Universidad Agraria de La Habana, La Habana, Cuba, 2011.

5. HENTZ, S. B.; L. DAUDEVILLE y F. D. R. V. DONZE: "Identification and Validation of a Discrete Element Model for Concrete", Journal of Engineering Mechanics, 130(6): 709-719, 2004.

6. HERRERA, M.; C. E. IGLESIAS; O. GONZÁLEZ; E. LÓPEZ y A. L. SÁNCHEZ: "Simulación mediante el Método de Elementos Finitos de la respuesta mecánica de un Oxisol", Revista Ciencias Técnicas Agropecuarias, 16(4): 55-61, 2008.

7. JEAN-YVES DELENNE, M. S. E. Y. F. C. J. C. B.: "Mechanical behaviour and failure of cohesive granular materials", International Journal for Numerical and Analytical Methods in Geomechanics, 28(15): 1577-1594, 2004.

8. LIAO, C. L.; T. P. CHANG; D. H. YOUNG y C. S. CHANG: "Stress-strain relationship for granular materials based on the hypothesis of best fit", International Journal of Solids and Structures, 34(31-32): 4087-4100, 1997.

9. LÓPEZ, E.; M. HERRERA; O. GONZALEZ; E. TIJSKENS y H. RAMON: "Determination of basics mechanical properties in a tropical clay soil as a function of dry bulk density and moisture", Technic Sciences in Agriculture, 21(3): 12-16, 2012.

10. SHMULEVICH, I.: "State of the art modeling of soil-tillage interaction using discrete element method", Soil and Tillage Research, 111(1): 41-53, 2010.

11. TIJSKENS, E.; H. RAMON y J. D. BAERDEMAEKER: "Discrete element modelling for process simulation in agriculture", Journal of Sound and Vibration, 266 493-514, 2003.

 

 

Recibido: 10 de septiembre 2012.
Aprobado: 14 de junio de 2013

 

 

Elvis López Bravo, Universidad Central de Las Villas, Dpto. Ing. Agrícola, Villa Clara, Cuba, CP 54830, Correo electrónico; elvislb@uclv.edu.cu