Indiana University-Purdue University Indianapolis (IUPUI) / A procedure for calibrating a discrete element (DE) computational soil model for various moisture contents using a conventional Asperity-Spring friction modeling technique is presented in this thesis. The procedure is based on the outcomes of two physical soil experiments: (1) Compression and (2) unconfined shear strength at various levels of normal stress and normal pre-stress. The Compression test is used to calibrate the DE soil plastic strain and elastic strain as a function of Compressive stress. To calibrate the DE inter-particle friction coefficient and adhesion stress as a function of soil plastic strain, the unconfined shear test is used. This thesis describes the experimental test devices and test procedures used to perform the physical terramechanical experiments. The calibration procedure for the DE soil model is demonstrated in this thesis using two types of soil: sand-silt (2NS Sand) and silt-clay(Fine Grain Soil) over 5 different moisture contents: 0%, 4%, 8%, 12%, and 16%. The DE based models response are then validated by comparing them to experimental pressure-sinkage results for circular disks and cones for those two types of soil over 5 different moisture contents. The Mean Absolute Percentage Error (MAPE) during the compression calibration was 26.9% whereas during the unconfined shear calibration, the MAPE was calculated to be 11.38%. Hence, the overall MAPE was calculated to be 19.34% for the entire calibration phase.

Fluidized beds have widespread industrial applications ranging from chemical industries to power plants. The flow inside a fluidized bed system consists of two main phases, a particle phase and the fluid phase. The two phases are strongly coupled to each other through various forces like drag and pressure. Capturing this multiphase phenomenon requires modeling strategies that possess good fidelity over a range of scales. Discrete Element Modeling (DEM) coupled with Computational Fluid Dynamics (CFD) provides a good platform to analyze the complex coupled multiphase hydrodynamics inside fluidized bed systems. Conventional DEM-CFD framework suffers from contradictory spatial resolution requirements for the particle and fluid phases, respectively. This prevents the conventional DEM-CFD method to be applied to geometries that have features comparable to the particle diameter of the solid phase. The novelty of this work lies in the development and validation of a two-grid formulation that removes the resolution restrictions of the conventional DEM-CFD framework. The results obtained from this new framework agree reasonably well with the experiments showing the capability of the new scheme to simulate conditions not possible with conventional DEM-CFD framework. In addition, this research also focuses on performing both 2D and 3D jetting fluidized bed simulations having millions of particles; validate/compare results with experiments and to perform heat transfer studies in a jetting fluidized bed system. The results suggest convective and diffusive mixing for a single jet at higher superficial velocity to be better than the mixing obtained in a multiple jet framework. The comparison with experimental results obtained in a multiple jetting setup shows that a 2D simulation captures the essential jet characteristics near the distributor plate reasonably well while a 3D simulation is needed to capture proper bubble dynamics near the freeboard of the bed. These results give insight into the detailed dynamics of fluidized bed systems and provide a foundation for a better design of these systems. / Ph. D.

Yes / Numerical simulation of wet granulation in a continuous granulator is carried out using Discrete Element Method (DEM) to discover the possibility of formation of seeded granules in a continuous process with the aim of reducing number of experimental trials and means of process control. Simple and scooped drum granulators are utilized to attain homogenous seeded granules in which the effects of drum rotational speed, particles surface energy, and particles size ratio are investigated. To reduce the simulation time a scale-up scheme is designed in which a dimensionless number (Cohesion number) is defined based on the work of cohesion and gravitational potential energy of the particles. Also a mathematical/numerical method along with a MATLAB code is developed by which the percentage of surface coverage of each granule is predicted precisely. The results show that use of continuous granulator is promising provided that a high level of shear is considered in the granulator design, e.g. it is observed that using baffles inside the drum granulators is essential for producing seeded granules. It is observed, moreover, that the optimum surface energy for scooped granulator with rotational speed of 50 rpm is 3 J/m2 which is close to the number predicted by Cohesion number. It is also shown that increasing the seed/fine size ratio enhances the seeded granulation both quantitatively (60% increase in seeds surface coverage) and qualitatively (more homogeneous granules).

In a proposed design for a concentrated solar power tower, sand is irradiated by solar energy and transfers its energy to another fluid stream by means of a finned tube heat exchanger. To maximize heat transfer and minimize potential damage to the heat exchanger, it is desired to have a very uniform flow through the heat exchanger. However, performing full scale flow tests can be expensive, impractical, and depending upon the specific quantities of interest, unsuitable for revealing the details of what it happening inside of the flow stream. Thus, the discrete element method has been used to simulate and study particulate flows. In this project, the flow of small glass beads through a square pyramid shaped hopper and a wedge shaped hopper were studied at the lab scale. These flows were also simulated using computers running two versions of discrete element modeling software – EDEM and LIGGGHTS. The simulated results were compared against the lab scale flows and against each other. They show that, in general, the discrete element method can be used to simulate lab scale particulate flows as long as certain material properties are well known, especially the friction properties of the material. The potential for increasing the accuracy of the simulations, such as using better material property data, non-uniform particle size distributions, and non-spherical particle shapes, as well as simulating heat transfer within a granular flow are also discussed.

Compressive residual stresses (CRS) are beneficial for enhancing the fatigue life of metal components. Shot Peening (SP) is an industrial cold working process that is applied to induce a field of CRS and modify the mechanical properties of the metal component. The SP process involves impacting a surface with tiny shots with forces sufficient to create plastic deformation. The process is governed by a number of important parameters such as the shot size, angle of attack, initial velocity, mass flow rate and the distance from the shot nozzle to the surface being peened. The relationship between the optimal peening outcome, particularly the residual stress distribution of the treated surface, and the peening parameters is still unknown and needs to be investigated further. Manufacturers are interested in producing a uniform peening process for complex geometries which optimises the SP parameters. Modelling the process is complex as it involves the interaction of a metallic surface with a large number of shots of very small diameter. Conventionally, such problems are solved using finite element software to predict stresses and strains of a single shot impact then applying superposition. At the moment there are no Finite Element Method (FEM) modelling solutions involving more than tens of shots. The number of shots and elements required for such a modelling process made the approach unfeasible prior to the work described herein. The objective of this work is to develop an appropriate numerical modelling approach that can better simulate the real SP process. The model will be provided by combining Discrete Element Method (DEM) with FEM. The DEM is employed to get a distribution of impact velocities over space and time which are then implemented into a FEM analysis. A discrete element model with randomly distributed steel shots bombarding a steel component at various velocities has been developed as benchmark example. With this model the SP shot - shot interaction, the shot - target interaction, the surface coverage, angle of impingement, shot size, impact velocity and the overall shot flow can be parametrically studied in details and with little computational effort. The novel approach also proposes a new method to dynamically change the coefficient of restitution for repeated impacts during the simulation and predicts the CRS more effectively. The effects of SP on different materials of relevance to gas turbine engine components will be investigated in order to improve the understanding of the interaction between the shots and the targeted material. Initially, an uncoupled analysis was peforned, in order to assess the capabilities of the two modelling systems, DEM and FEM, to delivery an improved solutuion when combining two commercially available codes. This parametric analysis is performed using the state-of-the-art Discrete Element (DE) application EDEM. In the subsequent part of this work, a dynamic Finite Element (FE) application Abaqus will be used to investigate single shot impacts and to obtain the residual stress distribution. This gives us a prescribed residual stress distribution and peening coverage. A Combined DEM/FEM tool (DEST) is proposed that eliminates any manual pre-processing required for linking/coupling, eliminating the use of two different applications and provide an integrated solution for the simulation of the Shot Peening process. In the subsequent chapter, the implementation of essential tools for the enchanced modelling of Shot Peening process functionalities, such as the nozzle, bounding box, coverage and intensity is described. A number of computational improvements are also implemented to reduce the computation time. The existing binary search is enhanced to self-balancing search tree and further improved to allow insertion and deletion of elements. A bounding box feature which removes shots that move out of the domain during the course of the simulation is also implemented. Experiments featuring single shot impacts are performed to gain better understanding the deformation process in the target material subjected to impact conditions to those occurring in the production peening. The single shot impacts are experimentally examined using SEM and EBSD. During final chapter, case studies are performed to compare the results of the simulations with large-scale experimental work. The coverage of peening of single and multiple nozzles with different angle of impingements are assessed. Finally, possible directions for further research concerning the accurate quantification of material responses to SP are identified in the report.

Dans cette thèse, nous avons étudié la simulation numérique des phénomènes d'interaction fluide-structure entre un fluide compressible et une structure déformable. En particulier, nous nous sommes intéressés au couplage par une approche partitionnée entre une méthode de Volumes Finis pour résoudre les équations de la mécanique des fluides compressibles et une méthode d'Eléments discrets pour le solide, capable de prendre en compte la fissuration. La revue des méthodes existantes de domaines fictifs ainsi que des algorithmes partitionnés couramment utilisés pour le couplage conduit à choisir une méthode de frontières immergées conservative et un schéma de couplage explicite. Il est établi que la méthode d'Eléments Discrets utilisée permet de retrouver le comportement macroscopique du matériau et que le schéma symplectique employé assure la préservation de l'énergie du solide. Puis nous avons développé un algorithme de couplage explicite entre un fluide compressible non-visqueux et un solide indéformable. Nous avons montré des propriétés de conservation exacte de masse, de quantité de mouvement et d'énergie du système ainsi que de consistance du schéma de couplage. Cet algorithme a été étendu au couplage avec un solide déformable, sous la forme d'un schéma semi-implicite. Cette méthode a été appliquée à l'étude de problèmes d'écoulements non-visqueux autour de structures mobiles : les comparaisons avec des résultats numériques et expérimentaux existants démontrent la très bonne précision de notre méthode / This work aims at the numerical simulation of compressible fluid/deformable structure interactions. In particular, we have developed a partitioned coupling algorithm between a Finite Volume method for the compressible fluid and a Discrete Element method capable of taking into account fractures in the solid. A survey of existing fictitious domain methods and partitioned algorithms has led to choose an Embedded Boundary method and an explicit coupling scheme. We first showed that the Discrete Element method used for the solid yielded the correct macroscopic behaviour and that the symplectic time-integration scheme ensured the preservation of energy. We then developed an explicit coupling algorithm between a compressible inviscid fluid and an undeformable solid. Mass, momentum and energy conservation and consistency properties were proved for the coupling scheme. The algorithm was then extended to the coupling with a deformable solid, in the form of a semi-implicit scheme. Finally, we applied this method to unsteady inviscid flows around moving structures: comparisons with existing numerical and experimental results demonstrate the excellent accuracy of our method

La déchirure musculaire est la première cause de blessure chez les athlètes. De nombreuses études décrivent ce traumatisme musculaire sans parvenir à en identifier clairement la chronologie et ses circonstances. L’objectif de la thèse est de décrire le phénomène de déchirure musculaire avec la méthode des éléments discrets, en s’appuyant sur des essais expérimentaux pour valider les modèles numériques. Dans une première partie, une revue de littérature permet d’acquérir les propriétés mécaniques des différents éléments constituant le complexe musculo-tendineux afin de pouvoir en réaliser un modèle macroscopique. Dans une deuxième partie, la modélisation du complexe musculo-tendineux est réalisée. La validation du comportement mécanique en traction passive du modèle proposé est réalisée en comparaison des travaux de L-L. Gras sur le muscle sternocléidomastoïdien humain. L’influence des paramètres morphologiques sur le comportement mécanique global est ensuite étudiée. La rupture fait l’objet de la troisième partie. Une modélisation de l’ensemble {tendon d’Achille/triceps sural} est réalisée et soumise à un test de traction passif jusqu’à rupture. La validation des résultats est faite vis-à-vis des essais expérimentaux réalisés sur cet ensemble musculaire provenant de pièces anatomiques humaines. L’étape suivante s’attache à modéliser la contraction musculaire, implémentée au niveau des fibres musculaires. Une validation du comportement actif du complexe musculo-tendineux est réalisée. Cette dernière étape, combinée à la traction destructive permet d’étudier la faisabilité de modéliser la déchirure par la méthode des éléments discrets, mais également d’étudier les structures endommagées et les mécanismes de rupture. Cela ouvre des possibilités d’utilisation cliniques de ce modèle pour comprendre et prévenir des blessures par déchirure musculaire. / Tearing of the muscle-tendon complex is a common sport-related injury for athletes. Many studies reported description of this traumatism but mechanisms leading to such an injury are still unclear as are the site of mechanical failure and involved structures. The aim of the thesis is to describe the phenomenon of the muscle-tendon-complex’s tear using the discrete element method and validating the numerical model with experimental data. In the first part, a literature review explains the different properties of the muscle-tendon complex main components’ in order to model it at the macroscopic scale. In the second part, the muscle-tendon complex is modeled. Validation of the mechanical behavior in passive tensile test is proposed by comparison with experimental data from L.-L. Gras on human sternocleidomastodeus muscle. Then, the different influences of morphometric parameters on the mechanical behavior of the complex are investigated. The third part focuses on the rupture. A model of the complex set of {Achilles tendon/surae triceps} is built and a tensile test until rupture is applied. Model validity is assessed by comparison with in vitro experiments from human cadavers. The fourth part focuses on the muscular activation, implemented inside fibers’ behavior. Validity of its active behavior is investigated. This fifth and last presents the enrichment with destructive tensile test. This added test allows first to study the feasibility to model the tear with the discrete element method; and second to focus on damaged structures and rupture’s mechanisms. This offers possibilities for clinical applications of this model to understand and prevent injuries caused by a tear of the muscle-tendon complex

Diversos materiais de interesse tecnológico podem ser considerados como materiais heterogêneos, onde sua natureza aleatória deve ser considerada para representar corretamente o comportamento não-linear. A avaliação quantitativa do dano em materiais sujeitos a estados de tensão ou deformação tem grande importância devido ao caráter crítico desse fenômeno, que pode crescer de forma abrupta e resultar em uma falha catastrófica da estrutura. Em estudos anteriores, Carpinteri e seus colaboradores apresentaram diferentes aspectos relacionados à caracterização do dano em materiais heterogêneos. Três desses aspectos são discutidos neste trabalho: (i) a proposta do número de fragilidade como medida da fragilidade da estrutura em análise; (ii) a evolução da dimensão fractal em que o processo de dano é desenvolvido; (iii) e os índices globais obtidos pela análise de Emissão Acústica (EA). Neste trabalho, uma versão do método dos elementos discretos formado por barras é utilizada para explorar esses aspectos. Espécimes de material quase-frágil são simulados e, quando possível, os resultados numéricos são comparados com resultados experimentais. Além disso, a discussão dos resultados obtidos é realizada visando compreender o comportamento desse tipo de material, como por exemplo, o fato de que estruturas com diferentes dimensões, porém com o mesmo número de fragilidade, têm comportamentos similares. Ademais, o método numérico empregado é apresentado como uma ferramenta viável para complementar a informação obtida em ensaios experimentais na avaliação do processo de dano. Pelo sinal de EA obtido através do método numérico empregado, obtiveram-se parâmetros para a classificação do tipo de falha e para o cálculo do tensor momento, onde comparando os resultados obtidos pelo modelo numérico e pela EA se observam comportamentos coerentes. / Several materials with technological interest can be considered as heterogeneous materials, where their random nature must be accounted to correctly represent the nonlinear behavior. The quantitative evaluation of damage in materials subjected to stress or strain states have great importance due to the critical nature of this phenomenon, which abruptly can rise to catastrophic failure. In previous studies, Carpinteri and his coworkers have presented different aspects of the damage process characterization in heterogeneous materials. Three of these aspects are discussed in this work: (i) the brittleness number proposal to measure the brittleness level of the analyzed structure, (ii) the relationship of the fractal dimension in which the damage process is developed, (iii) and the global indexes obtained by the Acoustic Emission (AE) analysis. In the present work, a version of discrete element method formed by bars is used to explore these concepts. A set of quasi-brittle material specimen is simulated and, when possible, the numerical results are compared with experimental ones. Moreover, a discussion of the obtained results is carried out aiming to understand the behavior of this kind of material, for instance, the fact that structures with different dimensions, but with the same brittleness number, have similar behavior. In addition, the numerical method is presented as a viable tool to complement information from experimental test on the damage process. From the AE signal obtained by the numerical method, parameters to classify the type of crack and for calculating the moment tensor were obtained, where consistent behaviors are observed by comparing the results of the numerical model and the AE.

Lors de la première transformation du bois en scierie, les grumes sont surfacées ou équarries par des têtes de fraisages appelées slabber ou canter. Sous leur action, le copeau de bois est fragmenté en plaquettes dont la valorisation est un enjeu majeur de la filière. Débouché le plus rémunérateur de cette ressource, l’industrie de la pâte à papier impose des critères dimensionnels aux plaquettes approvisionnées notamment concernant leur épaisseur. L’objectif de l’étude est donc d’améliorer la compréhension des mécanismes mis en jeu par la coupe du bois, dans une configuration simplifiée de coupe orthogonale, afin d’optimiser la granulométrie des plaquettes produites. Une campagne expérimentale de coupe sur machine-outil à commande numérique a été réalisée sur du hêtre vert afin de déterminer les mécanismes principaux actionnés. Un modèle mécanique simplifié du comportement dynamique du bois vert est déterminé, ceci afin de développer un modèle numérique du bois vert par la Méthode des Eléments Discrets (DEM) en vue de simuler sa coupe. Une étude préliminaire pour déterminer la capacité de la méthode à modéliser à l’échelle mésoscopique des milieux fortement orthotropes a été réalisée. Sa capacité et ses limites démontrées, la démarche de calibration du modèle numérique a été élaborée et le modèle sollicité en configuration de coupe orthogonale numérique. Les premières simulations présentent des résultats encourageants. / When entering sawmills, logs are faced into cants by the mean of chipper-canters. During this machining, the ribbon produced is split into small chips whose proper valorization is a high economic stake for the industry. The paper maker industry, which is the most worthwhile chips supplier, is strongly concerned by dimensional criterions of the chips for their process optimization, especially toward their thickness. The objective of this work is to improve the comprehension of cutting and fragmentation mechanisms, in a simplified orthogonal cutting configuration, to provide cutting rules to optimize the produced chips geometry. An experimental campaign dedicated to green beech cutting on a computer numerical command machining center is done to study those mechanisms. A simple mechanical model is derived, in order to develop a Discrete Element Method (DEM) model of the material to simulate cutting operations. A preliminary study dedicated to prove the capability of DEM to model a wood-like orthotropic continuous media is presented. The numerical model is then calibrated and cutting simulations are designed to copy the experimental conditions. First results are encouraging.

In this thesis, the effects of size and stress on permeability, deformability and strength of fractured rock masses are investigated. A comparison study was carried out to examine the effects of considering, or not considering, the correlation between distributions of fracture apertures and fracture trace lengths on the hydro-mechanical behavior of fractured rocks. The basic concepts used are the fundamental principles of the general theory of elasticity, Representative Elementary Volume (REV), the tensor of equivalent permeability, and the strength criteria of the fractured rocks. Due to the size and stress dependence of the hydro-mechanical properties of rock fractures, the overall effective (or equivalent) hydro-mechanical properties of the fractured rocks are also size and stress-dependent. However, such dependence cannot be readily investigated in laboratory using small samples, and so numerical modeling becomes a necessary tool for estimating their impacts. In this study, a closed-form relation is established for representing the correlation between a truncated lognormal distribution of fracture apertures and a truncated power law distribution of trace lengths, as obtained from field mapping. Furthermore, a new nonlinear algorithm is developed for predicting the relationship between normal stress and normal displacement of fractures, based on the Bandis model and the correlation between aperture and length. A large number of stochastic Discrete Fracture Network (DFN) models of varying sizes were extracted from some generated large-sized parent realizations based on a realistic fracture system description from a site investigation programme at Sellafield, UK, for calculating the REV of hydro-mechanical properties of fractured rocks. Rotated DFN models were also generated and used for evaluation of the distributions of directional permeabilities, such that tensors of equivalent permeability could be established based on stochastically established REVs. The stress-dependence of the permeability and the stress-displacement behaviour were then investigated using models of REV sizes. The Discrete Element Method (DEM) was used for numerical simulation of the fluid flow, deformability properties and mechanical strength behavior of fractured rocks. The results show significant scale-dependency of rock permeability, deformability and strength, and its variation when the correlation between aperture and trace length of fractures are concerned, with the overall permeability and deformability more controlled by dominating fractures with larger apertures and higher transmissivity and deformability, compared with fracture network models having uniform aperture. As the second moment of aperture distribution increases, a fractured rock mass shows more discrete behavior and an REV is established in smaller value of second moment with much larger model size, compared with the models with uniform fracture aperture. When the fracture aperture pattern is more scattered, the overall permeability, Young’s modulus and mechanical strength change significantly. The effect of stress on permeability and fluid flow patterns in fractured rock is significant and can lead to the existence or non-existence of a permeability tensor. Stress changes the fluid flow patterns and can cause significant channeling and the permeability tensor, and REV may be destroyed or re-established at different applied stress conditions. With an increase in the confining stress on the DEM models, the strength is increased. Compared with the Hoek-Brown criterion, the Mohr-Coulomb strength envelope provides a better fit to the results of numerical biaxial compression tests, with significant changes of the strength characteristic parameters occurring when the second moment of the aperture distribution is increased. / QC 20100702