The Discrete Element Method was first introduced by Cundall and Strack (1979) to model granular soils within the context of geotechnical engineering. The material is modeled as a random assembly of discrete elements. Each particle interacts with neighboring particles through contact forces that can be built up and broken at any time. The particles were modeled as discs in 2-D or as spheres in 3-D. Research studies have been conducted to improve the simulation of actual grain shapes. Ashmawy et al. (2003) developed the overlapping rigid clusters (ORC) method to accurately model irregular particle shapes. The idea relies on clumping a number of overlapping discs such that their coincides with that of the actual particle. In this dissertation, experimental verification program is presented. An experimental setup was built and model-grains were manufactured in the laboratory. A numerical simulation for the experimental test was carried out. The numerical and experimental results were compared qualitatively and quantitatively. A good agreement was observed within small displacements ranges. However, results were considerably different at large displacements. Numerical results utilizing the ORC method were closer to the experimental results than those of discs. A sequential and operator-independent procedure, which relies on the ORC concept, was developed. Identical inertial properties between the actual particle and the model were ensured. The new procedure was implemented for rounded and angular particles. The effect of particle shape and angularity on the strength and dilatancy characteristics of granular soils was investigated. A modified shape factor, which relies on the work introduced by Sukumaran and Ashmawy (2001), was developed. A series of pure shear testing simulations was performed on different shape and angularity particle groups. Angularity had a remarkable effect on strength and dilatancy properties compared to shape. The effect of interparticle friction on dilatancy was studied. An attempt was made to use an equivalent interparticle friction to model different particle shapes. It was concluded that there is no one-to-one equivalency between interparticle friction and shape or angularity. Instead, the interparticle friction must be continuously altered as a function of confining pressure and void ratio to achieve the required effect.

The general objective of the research effort reported in this thesis was to develop the knowledge required to optimize the design and operation of solid and semi-solid manure handling and land application equipment. Selected physical and rheological properties of manure products deemed to have an influence on the performances of manure handling and land application equipment were measured and general trends were identified among the measured properties. Relationships were also established between the measured properties and the type of manure as well as its total solids concentration. Field experiments were carried out to evaluate the effects of selected mechanical configurations, operating parameters and product properties on the discharge of manure spreaders. The influence of the type of conveying system (scraper conveyor and system of four augers) and the velocity at which it is operated, the geometry of the holding system and the position of a flow-control gate were all included in the analysis. The discharge rates of the machines as well as the specific energy required by the unloading operations were measured. A numerical modeling method called discrete element method (DEM) was used to create virtual manure, a numerical model of the real product. The measured physical and flow properties were used to develop and validate the virtual manure models. It was found that manure products could successfully be represented in a DE framework and that several parameters defining the contact constitutive model in the DEM had an influence on the behaviour of the virtual products. The DEM was then used to study machine-product interactions taking place in handling and land application equipment. Results from field experiments carried out using various land application equipment were used in the development and validation of the interaction models. The predicted flow rates and power requirements were in good agreement with measured data. The results obtained allowed for a better understanding of the flow of manure products in manure handling and land application equipment. It was found that the constitutive model used for the product influenced the results of the machine-product interactions models. A precision banded applicator under development at the University of Saskatchewan was also modeled. The discharge rate of this equipment is influenced by a number of parameters. The predicted mass distribution across the width of the banded applicator was well correlated to the experimental results. The models developed in this thesis have the potential to become powerful engineering tools for the design of improved machines for the handling and land application of solid and semi-solid manure.

The objective of the present study was to develop a fluid flow-coupled distinct element model capable of capturing the undrained behaviour of granular soils by considering fundamental physical mechanisms that involve fluid flow and particle interaction. The method considers granular media as assemblies of ellipsoidal particles arranged on a plane and interacting by means of contact forces. Saturation effects are incorporated by assuming that particles are immersed in fluid, the flow of which is simulated as occurring through a network of conduits. The flow through conduits is according to a Hagen-Poiseuille relation; a transient solution is obtained by solving a system of differential equations. The developed fluid-flow coupled distinct element was used to conduct various numerical simulations and the mechanisms of undrained deformations were examined from a micromechanical point of view. The dissertation begins with a literature review on the undrained behaviour of granular materials as observed in laboratory experiments. A review of previous attempts to simulate undrained tests micromechanically is also presented, and the advantages and disadvantages of various methods are examined. The capability of the developed model to simulate two-dimensional fluid-flow and pressure dissipation problems is demonstrated by means of comparisons with analytical solutions. Fluid pressure dissipation problems are qualitatively compared with Terzaghi's one-dimension theory of consolidation. It is shown that transient flow problems are accurately modelled by the fluid flow network approach. Simulated compression tests were carried out to examine the effects of different confining pressures and initial densities on the macroscopic response. The results compare favorably with those commonly observed in undrained laboratory experiments. Simulated tests are analyzed from a micromechanical point of view. It is shown that macroscopic behaviour can be traced to changes in micromechanical fabric descriptors. The effects of the interparticle friction angle on the undrained behaviour of the assemblies are investigated. The undrained strength is considerably increased by increasing interparticle friction. The main mechanism found to be responsible for the development of higher strength is the tendency of the specimens to dilate during shear distortion. The effects of the principal stress direction on the macroscopic response are examined. The behaviour of initially anisotropic samples is significantly altered by the direction of the principal stresses relative to the anisotropy direction. It is demonstrated that macroscopic permeability of the media has a considerable effect on the strength. This behaviour is attributed to the inhomogeneity of pore pressure distributions which increases with decreased permeability. The results presented are generally in agreement with observations previously reported from laboratory experiments. The possible applications of the model for future research are also discussed.

A two dimensional continuous numerical model based on Discrete Element Method is used to investigate the behaviour of accretionary wedges with different basal frictions. The models are based on elastic-plastic, brittle material and computational granular dynamics, and several characteristics of the influence of the basal friction are analysed. The model results illustrate that the wedge’s deformation and geometry, for example, fracture geometry, the compression force, area loss, displacement, height and length of the accretionary wedge etc., are strongly influenced by the basal friction. In general, the resulting wedge grows steeper, shorter  and higher, and the compression force is larger when shortened  above a larger friction basement.  Especially, when there is no basal friction, several symmetrical wedges will distribute symmetrically in the domain. The distribution of the internal stress when a new accretionary prime is forming is also studied. The results illustrate that when the stress in a certain zone is larger than a critical number, a new thrust will form there.

The objective of the present study was to develop a fluid flow-coupled distinct element model capable of capturing the undrained behaviour of granular soils by considering fundamental physical mechanisms that involve fluid flow and particle interaction. The method considers granular media as assemblies of ellipsoidal particles arranged on a plane and interacting by means of contact forces. Saturation effects are incorporated by assuming that particles are immersed in fluid, the flow of which is simulated as occurring through a network of conduits. The flow through conduits is according to a Hagen-Poiseuille relation; a transient solution is obtained by solving a system of differential equations. The developed fluid-flow coupled distinct element was used to conduct various numerical simulations and the mechanisms of undrained deformations were examined from a micromechanical point of view. The dissertation begins with a literature review on the undrained behaviour of granular materials as observed in laboratory experiments. A review of previous attempts to simulate undrained tests micromechanically is also presented, and the advantages and disadvantages of various methods are examined. The capability of the developed model to simulate two-dimensional fluid-flow and pressure dissipation problems is demonstrated by means of comparisons with analytical solutions. Fluid pressure dissipation problems are qualitatively compared with Terzaghi's one-dimension theory of consolidation. It is shown that transient flow problems are accurately modelled by the fluid flow network approach. Simulated compression tests were carried out to examine the effects of different confining pressures and initial densities on the macroscopic response. The results compare favorably with those commonly observed in undrained laboratory experiments. Simulated tests are analyzed from a micromechanical point of view. It is shown that macroscopic behaviour can be traced to changes in micromechanical fabric descriptors. The effects of the interparticle friction angle on the undrained behaviour of the assemblies are investigated. The undrained strength is considerably increased by increasing interparticle friction. The main mechanism found to be responsible for the development of higher strength is the tendency of the specimens to dilate during shear distortion. The effects of the principal stress direction on the macroscopic response are examined. The behaviour of initially anisotropic samples is significantly altered by the direction of the principal stresses relative to the anisotropy direction. It is demonstrated that macroscopic permeability of the media has a considerable effect on the strength. This behaviour is attributed to the inhomogeneity of pore pressure distributions which increases with decreased permeability. The results presented are generally in agreement with observations previously reported from laboratory experiments. The possible applications of the model for future research are also discussed.

The general objective of the research effort reported in this thesis was to develop the knowledge required to optimize the design and operation of solid and semi-solid manure handling and land application equipment. Selected physical and rheological properties of manure products deemed to have an influence on the performances of manure handling and land application equipment were measured and general trends were identified among the measured properties. Relationships were also established between the measured properties and the type of manure as well as its total solids concentration. Field experiments were carried out to evaluate the effects of selected mechanical configurations, operating parameters and product properties on the discharge of manure spreaders. The influence of the type of conveying system (scraper conveyor and system of four augers) and the velocity at which it is operated, the geometry of the holding system and the position of a flow-control gate were all included in the analysis. The discharge rates of the machines as well as the specific energy required by the unloading operations were measured. A numerical modeling method called discrete element method (DEM) was used to create virtual manure, a numerical model of the real product. The measured physical and flow properties were used to develop and validate the virtual manure models. It was found that manure products could successfully be represented in a DE framework and that several parameters defining the contact constitutive model in the DEM had an influence on the behaviour of the virtual products. The DEM was then used to study machine-product interactions taking place in handling and land application equipment. Results from field experiments carried out using various land application equipment were used in the development and validation of the interaction models. The predicted flow rates and power requirements were in good agreement with measured data. The results obtained allowed for a better understanding of the flow of manure products in manure handling and land application equipment. It was found that the constitutive model used for the product influenced the results of the machine-product interactions models. A precision banded applicator under development at the University of Saskatchewan was also modeled. The discharge rate of this equipment is influenced by a number of parameters. The predicted mass distribution across the width of the banded applicator was well correlated to the experimental results. The models developed in this thesis have the potential to become powerful engineering tools for the design of improved machines for the handling and land application of solid and semi-solid manure.

In the past, most research on soil mechanical properties was carried out for cohesionless soils in the fields of civil and geotechnical engineering. Little research has been carried out for mechanical properties of agricultural soil, which are essential for designing soil engaging tools in agriculture. In this study, unconsolidated undrained triaxial compression tests were performed to study the effects of moisture level and confining pressure on a sandy loam soil. The soil specimens tested had three moisture levels, and they were high (27-29% d.b.), medium (19-21% d.b.) and low (9-11% d.b.). The confining pressures used for the triaxial tests were 50, 100, and 150 kPa. Soil specimen was loaded at a strain rate of 1%/min. Measurements from the tests included stress-strain curve, shear strength, Young’s modulus, Poisson’s ratio, angle of internal friction, and cohesion. A model was developed using the Discrete Element Method (DEM) and computed by Particle Flow Code in three dimensions (PFC3D), a common DEM software. The model simulated the triaxial compression tests, and the model specimen was an assembly of 5-mm spherical particles which were defined by a set of micro parameters. During simulations, soil shear strength was monitored under different confining pressures. Through sensitivity analysis, it was found that most of the micro parameters affected the simulated soil shear strengths and the stress-strain behaviours. The most influential micro parameter was particle friction coefficient. This micro parameter was calibrated with the data from triaxial tests for different combinations of soil moisture levels and confining pressures. The calibrated particle friction coefficients varied from 0.2 to 1.0. The calibrations were done through matching the shear strengths between simulations and measurements, and the relative errors ranged between 0 and 6 %. / October 2015

Tobulėjant skaitiniams metodams ir kompiuterinei technikai atsivėrė galimybė naujų, sudėtingesnių mechaninių objektų modeliavimui. Turėdami naujus sudėtingesnių objektų modelius tyrėjai gali pritaikyti aprašytas ir sumodeliuotas šių objektų savybes su mikro struktūros ypatumais esamų ar busimų savybių nustatymui bei naujų medžiagų kūrimui. Tam intensyviai naudojami kaip eksperimentiniai taip ir skaitiniai metodai, kurių tobulinimui šiuo metu yra skiriamas labai didelis dėmesys. Skaitinis eksperimentas, kaip medžiagos tyrimo būdas pasitelkiamas dar ir todėl kad yra pigesnis ir leidžia interpretuoti jau žinomus eksperimentų rezultatus, o taip suteikia žinių naujiems tyrimams. Vienas iš metodų, kuris modeliuoja makroskopines medžiagų savybes remdamasis medžiagos mikro savybėmis yra diskrečiųjų elementų metodas (DEM). DEM metodas remiasi idėja, kad bet kokia fizikinė struktūra gali būti aprašyta kaip judančių dalelių sistema. Ši idėja pradėta taikyti ir vientisam deformuojamam kūnui aprašyti. Skirtingai nuo biriųjų medžiagų, vientiso kūno dalelės ir tarp jų egzistuojančios sąveikos yra kitokios prigimties, o jų modeliai yra fizikinės ir matematinės abstrakcijos rezultatas. Vientiso deformuojamo kūno modeliavimas diskrečiais elementais yra tik pradinėje stadijoje, o vientiso požiūrio į diskrečių elementų modelius dar nėra. Yra kelios hipotetinės versijos, grindžiamos skirtingais požiūriais. Taikant DEM kūnui, pirmas žingsnis būtu tampriųjų savybių modeliavimas. Tai yra... [toliau žr. visą tekstą] / Development of numerical methods and computation environments opened the possibility of new, more sophisticated mechanical objects modelling. In this context it is natural desire of the researchers to describe macroscopic mechanical characteristics of the materials by their microstructure, which can be adapted for simulation of the existing and future materials. For this purpose researchers are using intensively experimental and numerical methods for the development of which the highest priority is given. Numerical experiments are used because they are cheaper and allow the interpretation of already known results of experiments and provide information to new investigations. One of the methods used for modelling of macroscopic properties modelling is based on microscopic properties of material is discrete element method (DEM). The DEM traditionally was applied for the granular materials. The basic idea of DEM is that any physical structure could be described as a system of moving particles. This idea could be also applied to the description of solid deformable body. Particles forming solid body and existing interaction between them are of different nature than the granular materials because their models are often the result of physical and mathematical abstraction. The modelling of solid deformable body with the discrete elements is just at the initial stage and the unified approach to discrete elements models doesn’t exist. There are several versions of models, based... [to full text]

Development of numerical methods and computation environments opened the possibility of new, more sophisticated mechanical objects modelling. In this context it is natural desire of the researchers to describe macroscopic mechanical characteristics of the materials by their microstructure, which can be adapted for simulation of the existing and future materials. For this purpose researchers are using intensively experimental and numerical methods for the development of which the highest priority is given. Numerical experiments are used because they are cheaper and allow the interpretation of already known results of experiments and provide information to new investigations. One of the methods used for modelling of macroscopic properties modelling is based on microscopic properties of material is discrete element method (DEM). The DEM traditionally was applied for the granular materials. The basic idea of DEM is that any physical structure could be described as a system of moving particles. This idea could be also applied to the description of solid deformable body. Particles forming solid body and existing interaction between them are of different nature than the granular materials because their models are often the result of physical and mathematical abstraction. The modelling of solid deformable body with the discrete elements is just at the initial stage and the unified approach to discrete elements models doesn’t exist. There are several versions of models, based on... [to full text] / Tobulėjant skaitiniams metodams ir kompiuterinei technikai atsivėrė galimybė naujų, sudėtingesnių mechaninių objektų modeliavimui. Turėdami naujus sudėtingesnių objektų modelius tyrėjai gali pritaikyti aprašytas ir sumodeliuotas šių objektų savybes su mikro struktūros ypatumais esamų ar busimų savybių nustatymui bei naujų medžiagų kūrimui. Tam intensyviai naudojami kaip eksperimentiniai taip ir skaitiniai metodai, kurių tobulinimui šiuo metu yra skiriamas labai didelis dėmesys. Skaitinis eksperimentas, kaip medžiagos tyrimo būdas pasitelkiamas dar ir todėl kad yra pigesnis ir leidžia interpretuoti jau žinomus eksperimentų rezultatus, o taip suteikia žinių naujiems tyrimams. Vienas iš metodų, kuris modeliuoja makroskopines medžiagų savybes remdamasis medžiagos mikro savybėmis yra diskrečiųjų elementų metodas (DEM). DEM metodas remiasi idėja, kad bet kokia fizikinė struktūra gali būti aprašyta kaip judančių dalelių sistema. Ši idėja pradėta taikyti ir vientisam deformuojamam kūnui aprašyti. Skirtingai nuo biriųjų medžiagų, vientiso kūno dalelės ir tarp jų egzistuojančios sąveikos yra kitokios prigimties, o jų modeliai yra fizikinės ir matematinės abstrakcijos rezultatas. Vientiso deformuojamo kūno modeliavimas diskrečiais elementais yra tik pradinėje stadijoje, o vientiso požiūrio į diskrečių elementų modelius dar nėra. Yra kelios hipotetinės versijos, grindžiamos skirtingais požiūriais. Taikant DEM kūnui, pirmas žingsnis būtu tampriųjų savybių modeliavimas. Tai yra... [toliau žr. visą tekstą]

The shear rupture of massive (intact non-jointed) brittle rock in underground high stress mines occurs under a variety of different boundary conditions ranging from constant stress (no resistance to deformation) to constant stiffness (resistance to deformation). While a variety of boundary conditions exist, the shear rupture of massive rock in the brittle field is typically studied under constant stress boundary conditions. According to the theory, the fracturing processes leading to shear rupture zone creation occur at or near peak strength with a shear rupture surface created in the post-peak region of the stress-strain curve. However, there is evidence suggesting that shear rupture zone creation can occur pre-peak. Limited studies of shear rupture in brittle rock indicate pre-peak shear rupture zone creation under constant stiffness boundary conditions. This suggests that the boundary condition influences the shear rupture zone creation characteristics. In this thesis, shear rupture zone creation in brittle rock is investigated in direct shear under constant normal stress and normal stiffness boundary conditions. It is hypothesized that the boundary condition under which a shear rupture zone is created influences its characteristics (i.e., shear rupture zone geometry, load-displacement response, shear rupture zone creation relative to the load-displacement curve, and peak and ultimate strengths). In other words, it is proposed that the characteristics of a shear rupture zone are not only a function of the rock or rock mass properties but the boundary conditions under which the rupture zone is created. The hypothesis is tested and proven through a series of simulations using a two dimensional particle based Distinct Element Method (DEM) and its embedded grain based method. The understanding gained from these simulations is then used in the analysis and re-interpretation of rupture zone creation in two mine pillars. This is completed to show the value and practical application of the improved understanding gained from the simulations. The re-interpretation of these case histories suggests that one pillar ruptured predominately under a constant stress boundary condition while the other ruptured under a boundary condition changing from stiffness to stress control.