No description available.

The study of the shear behavior of particulate (soil) – continuum (man-made material) interfaces has received significant attention during the last three decades. The historical belief that the particulate – continuum interface represents the weak link in most geotechnical systems has been shown to be incorrect for many situations. Namely, prescribing properties of the continuum material, such as its surface roughness and hardness, can result in interface strengths that are equal to the contacting soil mass internal shear strength. This research expands the engineering implications of these findings by studying the response of interface systems in different loading conditions. Specifically, the axial and torsional shear modes are studied in detail. Throughout this thesis it is shown that taking an engineering approach to design the loading conditions induced to the interface system can result in interface strengths that exceed the previously considered limiting shear strength of the contacting soil. Fundamental experimental and numerical studies on specimens of different types of sand subjected to torsional and axial interface shear highlighted the inherent differences of these processes. Specifically, micro-scale soil deformation measurements showed that torsional shear induces larger soil deformations as compared to axial shear, as well as complex volume-change tendencies consisting of dilation and contraction in the primary and secondary shear zones. Studies on the global response of torsional and axial shear tests showed that they are affected differently by soil properties such as particle angularity and roughness. This difference in global behavior highlights the benefits of making systems that transfer load to the contacting soil in different manners available for use in geotechnical engineering. Discrete Element Modeling (DEM) simulations allowed for internal information of the specimens to be studied, such as their fabric and shear-induced loading conditions. These findings allowed for the development of links between the measured micro-scale behavior and the observed global-scale response. The understanding of the behavior of torsional and axial interfaces has allowed provides a framework for the development of enhanced geotechnical systems and applications. The global response of torsional shear found to induce larger cyclic contractive tendencies within the contacting soil mass. Therefore, this shear mode is more desirable than the conventional axial shear for the study of phenomena that depend on soil contractive behavior, such as liquefaction. A study on the influence of surface roughness form revealed that surfaces with periodic profiles of protruding elements that prevent clogging are capable of mobilizing interface friction angles that are 20 to 60% larger than the soil friction angle. These findings have direct implications in engineering design since their implementation can result in more resilient and sustainable geotechnical systems.

Doctor of Philosophy / Department of Biological & Agricultural Engineering / Mark E. Casada / Ronaldo G. Maghirang / The United States grain handling infrastructure is facing major challenges to meet worldwide customer demands for wholesome, quality, and safe grains and oilseeds for food and feed. Several challenges are maintaining grain quality during handling; reducing dust emissions for safety and health issues; growing shift from commodity-based to specialty (trait-specific) markets; proliferation of genetically modified crops for food, feed, fuel, pharmaceutical, and industrial uses; and threats from biological and chemical attacks. This study was conducted to characterize the quality of grain and feed during bucket elevator handling to meet customer demand for high quality and safe products. Specific objectives were to (1) determine the effect of repeated handling on the quality of feed pellets and corn; (2) characterize the dust generated during corn and wheat handling; (3) develop and evaluate particle models for simulating the flow of grain during elevator handling; and (4) accurately simulate grain commingling in elevator boots with discrete element method (DEM). Experiments were conducted at the research elevator of the USDA-ARS Center for Grain and Animal Health Research (CGAHR) to determine the effect of repeated handling on the quality of corn-based feed pellets and corn. Repeated handling did not significantly influence the durability indices of feed pellets and corn. The feed pellets, however, had significantly greater breakage (3.83% per transfer) than the corn (0.382% per transfer). The mass of particulate matter < 125 μm was less for feed pellets than for corn. These corn-based feed pellets can be an alternative to corn in view of their handling characteristics. Another series of experiments was conducted in the same elevator to characterize the dust generated during corn and wheat handling. Dust samples were collected from the lower and upper ducts upstream of the cyclones in the elevator. Handling corn produced more than twice as much total dust than handling wheat (185 g/t vs. 64.6 g/t). Analysis of dust samples with a laser diffraction analyzer showed that the corn samples produced smaller dust particles, and a greater proportion of small particles, than the wheat samples. Published data on material and interaction properties of selected grains and oilseeds that are relevant to DEM modeling were reviewed. Using these material and interaction properties and soybeans as the test material, the DEM fundamentals were validated by modeling the flow of soybean during handling with a commercial software package (EDEM). Soybean kernels were simulated with single- and multi-sphere particle shapes. A single-sphere particle model best simulated soybean kernels in the bulk property tests. The best particle model had a particle coefficient of restitution of 0.6; particle static friction of 0.45 for soybean-soybean contact (0.30 for soybean-steel interaction); particle rolling friction of 0.05; normal particle size distribution with standard deviation factor of 0.4; and particle shear modulus of 1.04 MPa. The single-sphere particle model for soybeans was implemented in EDEM to simulate grain commingling in a pilot-scale bucket elevator boot using 3D and quasi-2D models. Pilot-scale boot experiments of soybean commingling were performed to validate these models. Commingling was initially simulated with a full 3D model. Of the four quasi-2D boot models with reduced control volumes (4d, 5d, 6d, and 7d; i.e., control volume widths from 4 to 7 times the mean particle diameter) considered, the quasi-2D (6d) model predictions best matched those of the initial 3D model. Introduction of realistic vibration motion during the onset of clear soybeans improved the prediction capability of the quasi-2D (6d) model. The physics of the model was refined by accounting for the initial surge of particles and reducing the gap between the bucket cups and the boot wall. Inclusion of the particle surge flow and reduced gap gave the best predictions of commingling of all the tested models. This study showed that grain commingling in a bucket elevator boot system can be simulated in 3D and quasi-2D DEM models and gave results that generally agreed with experimental data. The quasi-2D (6d) models reduced simulation run time by 29% compared to the 3D model. Results of this study will be used to accurately predict impurity levels and improve grain handling, which can help farmers and grain handlers reduce costs during transport and export of grains and make the U.S. grain more competitive in the world market.

Thesis (MScEng)--Stellenbosch University, 2012. / ENGLISH ABSTRACT: In this study, tool forces from rock cutting tests were numerically simulated through a discrete element method (DEM) in association with PFC3D™. Tribological interactions such as contact, shearing, fracturing, friction and wear were presented during these cutting simulations. Particle assemblies, representing Paarl granite and Sandstone-2, were created in PFC3D™ through a material-genesis procedure. The macro-properties of these particle assemblies, namely Young’s modulus, Poisson’s ratio, uniaxial and triaxial compressive strength and Brazilian tensile strength, were calibrated by modelling the uniaxial and triaxial compressive strength test and the Brazilian tensile strength test. The calibration was done through adjustment of the micro-properties of the assembly, namely the stiffness and strength parameters of the particles and bonds. The influence of particle size on the calibration was also investigated. These assemblies were used in the rock cutting tests. Results suggested that DEM can reproduce the damage formation during calibration tests successfully. From the results obtained from the calibration tests, it was also concluded that particle size is not a free parameter but influences the macro-properties greatly. Different rock cutting tools were simulated, namely point-attack (conical) picks, chisel-shaped tools and button-shaped tools. The numerical cutting tools were treated as rigid walls to simplify the simulation and the tool forces were not influenced by wear. In each simulation the cutting tools advanced at a constant velocity. The tool forces acting on the cutting tool, in three orthogonal directions, were recorded during the numerical simulations and the peak cutting forces were predicted by theoretical equations. The damage to the Paarl granite and Sandstone-2 assemblies was revealed as broken bonds, which merge into microscopic fractures. The mean peak cutting forces of sharp cutting tools obtained from numerical, theoretical and experimental models (from the literature) were compared. Finally the influence of factors, including wear on the tool and depth of cut, on the value of tool forces was also investigated. The results from the rock cutting tests revealed that the correlation between the numerical and the experimental models as well as the theoretical and experimental models was not strong when using sharp point-attack and chisel-shaped picks. It was concluded that the influence of wear plays a substantial part in the cutting process and it has to be included during the numerical simulation for the results to be accurate and verifiable. This study also found that there is a non-linear increase in tool forces with an increase in depth of cut, since the contact area increases. At larger cutting depths, chip formation also generally increased and therefore damage to the sample as well as wear on the cutting tool will be minimized at shallow cutting depths. Overall this study concludes that DEM are capable of simulating calibration methods and rock cutting processes with different cutting tools and producing results which are verifiable with experimental data. Therefore numerical prediction of tool forces will allow the design of efficient cutting systems and the operational parameters as well as the performance prediction of excavation machines. / AFRIKAANSE OPSOMMING: In hierdie studie is die kragte wat tydens rotssny-toetse op die sny gereedskap inwerk, numeries gesimuleer met behulp van ‘n diskrete element metode (DEM) in samewerking met PFC3D™. Tribologiese interaksies soos kontak, skeer, breking, wrywing en slytasie is gedurende hiersie snytoetse voorgestel. Partikel versamelings, wat Paarl graniet en Sandsteen-2 verteenwoordig, is in PFC3D™ geskep deur middel van ‘n materiaal-skeppings prosedure. Die makro-eienskappe van die partikel versamelings, naamlik Young se modulus, Poisson se verhouding, eenassige en drie-assige druksterkte en Brasiliaanse treksterkte, is gekalibreer deur modellering van die eenassige en drie-assige druksterkte toets en die Brasiliaanse treksterkte toets. Die kalibrasie is gedoen deur aanpassing van die mikro-eienskappe, naamlik die styfheid en die sterkte parameters van die partikels en bindings. Die invloed van partikelgrootte is ook ondersoek. Daarna is hierdie versamelings in die rotssny-toetse gebruik. Resultate het daarop gedui dat DEM die kraakvorming gedurende kalibrasie toetse suksesvol kan reproduseer. Vanuit die kalibrasie is ook gevind dat die partikelgrootte nie ‘n vrye parameter is nie, maar die makro-eienskappe grotendeels beïnvloed. Verskillende rotssny gereedskap is gesimuleer, naamlik koniese, beitel-vormige en knopie-vormige instrumente. Die numeriese sny gereedskap is gesimuleer as rigiede mure om simulasies te vereenvoudig en die gereedskap-kragte is dus nie deur slytasie beïnvloed nie. Tydens elke simulasie is die sny gereedskap vorentoe beweeg teen ‘n konstante snelheid. Die gereedskap-kragte, in drie ortogonale rigtings, is aangeteken gedurende die numeriese simulasies en die piek snykragte is ook voorspel deur teoretiese vergelykings. Die skade aan die Paarl graniet en Sandsteen-2 versamelings, is voorgestel as gebreekte bindings, wat saamsmelt tot mikroskopiese frakture. Die gemiddelde piek snykragte van skerp sny gereedskap van numeriese, teoretiese en eksperimentele modelle (uit die literatuur) is vergelyk. Ten slotte is die invloed wat faktore, onder andere die slytasie van gereedskap en die snydiepte, op die grootte van die kragte het ondersoek. Die resultate van die rotssny-toetse het aan die lig gebring dat die korrelasie tussen die numeriese en eksperimentale modelle sowel as die teoretiese en eksperimentele modelle nie sterk is tydens die gebruik van skerp koniese en beitel-vormige instrumente nie. Die gevolgtrekking is gemaak dat die invloed van slytasie van sny gereedskap ‘n wesenlike rol speel in die snyproses en dat dit in die numeriese simulasie ingesluit moet word sodat die resultate akkuraat en virifieerbaar is. Hierdie studie het ook gevind dat daar ‘n nie-lineêre toename in die gereedskap-kragte is met ‘n toename in snydiepte aangesien die kontak-area toeneem met ‘n toename in die snydiepte. By groter snydieptes, het die formasie van afsplinterings verhoog en dus sal skade aan die partikel versamelings en die slytasie van die gereedskap geminimeer word by vlakker snydieptes. Algeheel het die studie tot die gevolgtrekking gekom dat DEM in staat is om kalibrasie metodes en rotssny-toetse met verskillende sny gereedskap te simuleer asook om resultate te produseer wat verifieerbaar is met eksperimentele data. Numeriese voorspellings van die gereedskap-kragte sal dus toelaat om doeltreffende sny prosesse en operasionele parameters te ontwerp sowel as om die werkverrigting van uitgrawings masjiene te voorspel.

Thesis (MScEng (Mechanical and Mechatronic Engineering))--University of Stellenbosch, 2011. / ENGLISH ABSTRACT: The removal of grape berries from the stems is an important step in the wine making process. Various problems are experienced using the destemming machines currently available, where the berries are mechanically removed and separated from the stems by a rotating beater shaft and drum. Not all berries are removed from the stems and broken stems can end up with the removed berries which can result in unwanted characters and flavours in the wine. The development of these machines is currently limited to experimental tests. In this study, the destemming process was investigated experimentally. The ability of the Discrete Element Method (DEM) to simulate this process was also investigated. A range of experiments was designed to obtain the material properties of the grapes. These experiments included the measurement of the stem stiffness and break strength, the berry stiffness, and the force needed to remove a berry from the stem. Experiments were conducted to gain further insight into the destemming process. Firstly, a simplified destemming machine with only a beater shaft and a single grape bunch was built. The influence of the bunch size and the speed of the beater shaft on the number of berries removed from the stems were investigated. Secondly, field tests on a commercial destemming machine were conducted and the performance of the machine was measured. A DEM model of both the simplified and the commercial destemming machine were built. Commercial DEM software was used with linear contact and bond models. The stems were built from spherical particles bonded together and a single spherical particle was used to represent each berry. The measured stiffnesses and break strengths were used to set the particle and bond properties. Modelling the simplified destemming machine, it was found that the DEM model could accurately predict the effect of the bunch size and the speed of the beater shaft on the number of berries removed from the stems. The model of the commercial destemming machine could accurately predict the machine’s performance in terms of the number of berries removed as well as the number of broken stems. / AFRIKAANSE OPSOMMING: Die verwydering van druiwekorrels vanaf die stingels is ʼn belangrike stap tydens die wynmaak proses. Verskeie probleme word ondervind met huidige beskikbare ontstingelaars, waar die korrels meganies verwyder en skei word vanaf die stingels deur middel van ʼn roterende klop-as en drom. Nie alle korrels word vanaf die stingels verwyder nie en gebreekte stingels kan saam met die verwyderde korrels beland, wat ongewensde karakters en geure in die wyn kan veroorsaak. Die ontwikkeling van ontstingelaars is tans beperk tot eksperimentele toetse. In hierdie studie is die ontstingel proses eksperimenteel ondersoek Die vermoë van die Diskrete Element Metode (DEM) om hierdie proses te simuleer is ook ondersoek. ʼn Reeks eksperimente is ontwikkel om die materiaal eienskappe van die druiwe te bepaal. Hierdie eksperimente sluit in die meet van die styfheid en breeksterkte van die stingel, die korrel styfheid, en die krag benodig om ʼn korrel vanaf die stingel te verwyder. Eksperimente is gedoen om verdere insig oor die ontstingel proses te bekom. Eerstens is ʼn vereenvoudigde ontstingelaar gebou, met slegs ʼn klop-as en een tros. Die invloed van die trosgrootte en die klop-as spoed op die aantal korrels wat verwyder is, is ondersoek. Tweedens is ʼn toets in die veld gedoen met ʼn kommersiële ontstingelaar om die werkverrigting van die masjien te bepaal. ʼn DEM model van beide die vereenvoudigde en kommersiële ontstingelaar is gebou. Kommersiële DEM sagteware is gebruik met lineêre kontak- en bindingsmodelle. Die stingels is gebou deur sferiese partikels aan mekaar te bind en ʼn enkele sferiese partikel is gebruik om ʼn druiwe korrel voor te stel. Die gemete styfhede en breeksterktes is gebruik om die partikel- en bindingseienskappe te spesifiseer. Die modellering van die vereenvoudigde ontstingelaar het getoon dat die DEM model akkuraat kan voorspel wat die invloed is van die trosgrootte en die klop-as spoed op die aantal korrels wat verwyder is. Die model van die kommersiële ontstingelaar kon die werkverrigting van die masjien akkuraat voorspel in terme van die aantal korrels wat verwyder is asook die aantal gebreekte stingels.

Thesis (MScEng)--Stellenbosch University, 2013. / ENGLISH ABSTRACT: Production of graphite dust inside the Pebble Bed Modular Reactor (PBMR) influences the reactor operation negatively. Graphite is used as a moderator in the reactor core and the formation and transportation of graphite dust away from the reactor core decreases the amount of moderator which in turn has a negative impact on the reactor operation. High levels of radioactive dust may also contaminate reactor components which may pose a health risk to maintenance personnel. In this study a pressure vessel was designed and used to measure the wear of a graphite pebble in helium at elevated temperatures. By means of a multi-linear regression analysis a proper mathematical function was established in order to relate graphite wear to certain tribological parameters. These parameters were identified through a literature study. Discrete Element Modelling (DEM) was used to simulate the gravitational flow of graphite pebbles through the reactor core. The experimentally determined mathematical function was incorporated into the DEM simulation to estimate the annual mass of graphite dust to be produced by the PBMR pebble bed as a result of pebble-pebble interaction and pebble-wall interaction during refuelling. / AFRIKAANSE OPSOMMING: Die vorming van grafiet stof binne die korrelbed-modulêre reaktor (PBMR) beïnvloed die werking daarvan negatief. Grafiet word gebruik as 'n moderator in die reaktor kern en die vorming en vervoer van grafietstof weg van die reaktor kern lei tot 'n afname in die hoeveelheid moderator en dit het 'n negatiewe impak op die werking van die reaktor. Hoë vlakke van radioaktiewe grafietstof kontamineer ook reaktorkomponente wat 'n gesondheidsrisiko vir onderhoudspersoneel inhou. In hierdie studie was 'n drukvat ontwerp en gebruik om die slytasie van 'n grafietkorrel in helium by verhoogde temperature te meet. 'n Multi-lineêre regressie analise is dan gebruik om 'n wiskundige funksie daar te stel wat die verband tussen grafietslytasie en die eksperimentele parameters vas stel. Hierdie parameters was met behulp van 'n literatuurstudie geïdentifiseer. Diskrete Element Modellering (DEM) was gebruik om die gravitasionele vloei van grafietkorrels in die reaktor te modelleer. Die eksperimenteel bepaalde wiskundige funksie word in die DEM simulasie ge-inkorporeer om 'n skatting te maak van die jaarlikse massa grafietstof wat gevorm sal word in die PBMR korrelbed as 'n gevolg van korrel-korrel interaksie en korrel-wand interaksie gedurende hersirkulasie.

The dissertation details work aimed toward the development and implementation of a 3-D impact fragmentation module to perform rock fall analysis by taking into account impact fragmentation. This fragmentation module is based on a database of a large set of impact simulations using a fully calibrated discrete element model (DEM), and is employed to predict impact fragmentation processes in rockfall analysis by either training a neural network model or linearly interpolating the database. A DEM was employed to model impact fragmentation in the study. A DEM code was developed from scratch. The model was first calibrated and verified with experimental results to demonstrate the capability of modeling both quasi-static and dynamic material behavior. Algorithms to calibrate the model’s micro-parameters against triaxial tests on rocks were presented. Sensitivity analyses were used to identify the deformability micro-parameters by obtaining relationships between microscopic and macroscopic deformability properties. The strength model parameters were identified by a global optimization process aimed at minimizing the difference between computed and experimental failure envelopes. When applied to the experimental results of tested granite, this calibration process produced a good agreement between simulated and experimental results for both deformability and strength properties. Dynamic compression and SHPB tests were performed to verify the dynamic model. A strain-rate-dependent dynamic strength was observed in the experimental results. This strain-rate-dependent dynamic strength was also confirmed by the numerical results. No rate-dependent constitutive model was used in the DEM to simulate dynamic behavior. This simulated rate-dependent dynamic strength can be attributed to material inertia because the inertia inhibits crack growth. Some fundamental mechanisms of impact fragmentation associated with rockfalls were then numerically investigated. The developed DEM code was coupled with a simplified impact model inspired by the theory of dynamic foundations. It has been shown that the magnitude of impact velocity, the angle of the incidence, the ground condition all play very important roles in impact fragmentation. Several case studies were performed to validate the developed impact fragmentation module in rock fall analysis. It has been demonstrated that the developed fragmentation module can reasonably predict impact fragmentation and perform some risk analysis in rock fall analysis. / text

The Distinct Element Method (DEM) and Discontinuous Deformation Analysis (DDA) are the two most commonly used discrete element methods in rock mechanics. Discrete element approaches are computationally expensive as they involve the interaction of multiple discrete bodies with continuously changing contacts. Therefore, it is very important to ensure that the method selected for the analysis is computationally efficient. In this research, a general assessment of DDA and DEM is performed from a computational efficiency perspective, and relevant enhancements to DDA are developed. The computational speed of DDA is observed to be considerably slower than DEM. In order to identify reasons affecting the computational efficiency of DDA, fundamental aspects of DDA and DEM are compared which suggests that they mainly differ in the contact mechanics, and the time integration scheme used. An in-depth evaluation of these aspects revealed that the openclose iterative procedure used in DDA which exhibits highly nonlinear behavior is one of the main reasons causing DDA to slow down. In order to improve the computational efficiency of DDA, an alternative approach based on a more realistic rock joint behavior is developed in this research. In this approach, contacts are assumed to be deformable, i.e., interpenetrations of the blocks in contact are permitted. This eliminated the computationally expensive open-close iterative procedure adopted in DDA-Shi and enhanced its speed up to four times. In order to consider deformability of the blocks in DDA, several approaches are reported. The hybrid DDA-FEM approach is one of them, although this approach captures the block deformability quite effectively, it becomes computationally expensive for large-scale problems. An alternative simplified uncoupled DDA-FEM approach is developed in this research. The main idea of this approach is to model rigid body movement and the block internal deformation separately. Efficiency and simplicity of this approach lie in keeping the DDA and the FEM algorithms separate and solving FEM equations individually for each block. Based on a number of numerical examples presented in this dissertation, it is concluded that from a computational efficiency standpoint, the implicit solution scheme may not be appropriate for discrete element modelling. Although for quasi-static problems where inertia effects are insignificant, implicit schemes have been successfully used for linear analyses, they do not prove to be advantageous for contact-type problems even in quasi-static mode due to the highly nonlinear behavior of contacts.

Computer simulation using the Discrete Element Method (DEM) has emerged as a powerful tool in studying the behaviour of particulate systems during powder flow and compaction. Contact law between particles is the most important input to the Discrete Element simulation. However, most of the present simulations employ over-simplistic contact laws which cannot capture the real behaviour of particulate systems. For example, plastic yielding, material brittleness, sophisticated particle geometry, surface roughness, and particle adhesion are all vitally important factors affecting the behaviour of particle interactions, but have been largely ignored in most of the DEM simulations. This is because it is very difficult to consider these factors in an analytical contact law which has been the characteristic approach in DEM simulations. This thesis presents a strategy for obtaining the contact laws numerically and a comprehensive study of all these factors using the numerical approach. A numerical method, named as the Material Point Method (MPM) in the literature, is selected and shown to be ideal to study the particle interactions. The method is further developed in this work in order to take into account all the factors listed above. For example, to study the brittle failure during particle impact, Weibull’s theory is incorporated into the material point method; to study the effect of particle adhesion, inter-atomic forces are borrowed from the Molecular Dynamic model and incorporated into the method. These developments themselves represent a major progress in the numerical technique, enabling the method to be applied to a much wider range of problems. The focus of the thesis is however on the contact laws between extremely fine particles. Using the numerical technique as a tool, the entire existing theoretical framework for particle contact is re-examined. It is shown that, whilst the analytical framework is difficult to capture the real particle behaviour, numerical contact laws should be used in its place.

Lunar soil, defined as the finest part of the lunar regolith which covers the entire surface of the Moon, has shown to have remarkable shear strength properties, highlighted by the clearly visible effects of soil cohesion. The main objective of this thesis is to unveil the physical explanations causing this unusual soil behaviour in a waterless, airless, lunar environment. Ultra-High Vacuum (UHV), in particular, is considered responsible for increasing the strength of surface energy forces due to lunar soil outgassing. In turn, the presence of surface energy forces, arising from van der Waals intermolecular forces, is thought to alter the mechanical properties of lunar soil. A particle-based microscopic approach by means of the Discrete Element Method (DEM) was utilised to investigate the effects of surface energy forces on the macroscopic soil be- haviour. A micro-mechanical contact model, based on the JKR theory, was selected to describe the inter-granular behaviour between lunar soil particles. Physical and geometrical parameters typical of lunar soil were employed. Several triaxial tests were run to identify a link, if any, between the microscopic surface energy parameter and the macroscopic soil cohesion, which was interpreted as a true soil cohesion. In addition, very low stress levels and high soil densities were simulated in order to take into account the low gravitational field and the high state of soil compaction caused by continuous meteorite impacts on the Moon. Results from triaxial tests were analysed at both the peak and critical state. It was found that in the ideal case of perfectly spherical grains, the presence of adhesion is a source of noticeable macroscopic soil cohesion. However, no influence was observed in terms of macroscopic friction angle. Furthermore, a brittle macroscopic soil behaviour was revealed, owing to the simulated inter-granular chemical bonds and the very low stress conditions applied. Finally, similar to the behaviour of cemented sands, very little cohesion was recorded at the critical state. Subsequently, particle shape effects were investigated by complementing the numerical model with a simple form of inter-particle rolling resistance. Simulations were also run with non-convex grains of increasing geometrical complexity in order to simulate more realistically the irregular shapes of lunar soil grains. In both cases, the interplay of surface energy forces with particle shape effects resulted in even higher shear strength, with predictions similar to the estimates of shear strength for real lunar soil. Once again, the peak strength was dominated by macroscopic cohesion which, on the other hand, was hardly observable at the critical state, confirming the tendency observed from spherical grains. Finally, the practical implications of the above findings were discussed in terms of bearing capacity, trafficability and slope stability on the lunar surface. In particular an analytical approach, based on the bearing capacity problem, was devised to study the performance of a rigid wheel rotating on a lunar terrain and operating under different dynamic conditions.