The effects of oxidation in the thermomechanical response of porous titanium diboride have been investigated. An in-house quasi-static material point method tool was used to perform two -dimensional plane strain simulations on unoxidized hexagonal representative volume elements (RVEs) with macroporosity volume fractions of 10%, 40% and 70% to establish a baseline for the response due to geometric effects. Compressive strains of up to 30% were applied at room temperature. The 10% and 40% RVEs showed shear banding and subsequent shear failure of the inter-pore struts, while shear banding in 70% RVE weakened the struts, which lead to buckling failure. A snapshot oxidation model was then applied to the hexagonal RVEs in place of a transient, diffusion-based oxidation solver. Compressive strain simulations were performed on RVEs with oxide layers ranging from 5 to 50 μm. In RVEs with porosity of 40% or higher, oxide percolation in the struts reduced the effective elastic modulus and compressive strength, though further oxidation beyond the percolation point did not have a significant impact. Ramped and cyclic thermal loads were applied and the damage due to thermal expansion coefficient mismatch at the oxide-substrate interface decreased as the oxide layer was increased. Finally, the snapshot oxidation modeling approach was applied to large porous RVEs derived from micro-computed tomography images of titanium diboride foam. The effective elastic modulus decreased by 47% when the 5 μm layer was applied due to many thin, flexible struts becoming fully oxidized. Subsequent oxidation did not have a significant impact on the thermomechanical response. / Master of Science / Thermal loading experienced by hypersonic flight vehicles has posed significant design challenges in the development of platforms for military and re-entry applications. The advent of hypersonic strike weapons and waveriders has led to an interest in utilizing ceramics with melting points above 3000°C, called ultra-high temperature ceramics (UHTCs), that offer improved resistance to high-temperature oxidation. Beyond load-carrying applications, UHTCs imbued with macroscale porosity have been introduced as candidates for providing thermal insulation of sensitive on-board components. This thesis presents a first pass at modeling the coupled effects of oxidation and continuum damage in the thermomechanical response of such materials. Using an in-house material point method tool, two-dimensional compressive strain simulations were performed on hexagonal representative volume elements (RVEs) of titanium diboride foam with varying levels of macroporosity, along with large porous RVEs derived from micro-computed tomography images of titanium diboride foam. A snapshot oxidation model was applied to these RVEs in place of a transient, diffusion-based oxidation solver, then simulations with applied compressive strains of up to 30% were performed on RVEs with oxide layers ranging from 5 to 50 μm. Ramped and cyclic thermal loads were applied to explore the effects of thermal expansion mismatch between the substrate and oxide phases. The oxide layers were shown to reduce the effective stiffness, compressive strength, and thermal conductivity of the RVEs, with the oxidation state of the inter-pore struts having a large impact on the overall material response.

Dans les régions de montagne, la prévision des évènements gravitaires reste un défi pour la gestion des risques. Des méthodes de calcul telles que la méthode des éléments discrets (DEM), où les particules interagissent les unes avec les autres pour restituer un comportement global d’une masse granulaire, ont été utilisées pour aborder ce type de problématique. L’application de la DEM reste normalement limitée aux évènements de petits volumes impliquant un nombre de blocs plutôt faible, puisque les temps de calcul peuvent devenir rapidement prohibitifs avec l’augmentation du nombre de particules. Les méthodes de calcul continues sont donc une alternative intéressante car elles permettent de réduire les temps de calcul. Elles nécessitent cependant la définition d’une loi de comportement macroscopique capable de représenter correctement les principaux traits de comportement mécanique du matériau au sein de la masse. L'objectif principal du travail de thèse réside dans le développement d’un outil numérique permettant de modéliser certains aléas gravitaires tels que les écoulements en masse. Notre choix s’est porté sur une méthode Lagrangienne-Eulérienne (méthode des points matériels – MPM) capable de gérer de grandes déformations tout en bénéficiant des principaux avantages de la méthode des éléments finis (FEM). La méthode utilise une grille Eulérienne fixe sur laquelle se déplacent des points matériels pendant les simulations. Un outil numérique, nommé MPMbox (2D et 3D), a été développé entièrement durant la thèse en C++. Le code a été validé à l'aide d'une série de solutions analytiques en quasi-statique (tests géotechniques standards) ainsi que par des applications de la littérature incluant des déformations importantes et rapides (tests d'affaissement). Après validation, le code a été confronté aux prédictions d’un outil de calcul DEM (DEMbox) dans le cadre de simulations numériques impliquant l'écoulement (initiation, régimes transitoires, propagation et arrêt) d'un matériau granulaire (particules sphero-polyhédriques) sur un plan incliné. Les résultats ont été comparés en termes de distance de propagation, de forme du dépôt et d'énergies dissipées à l'interface et dans la masse pendant l'écoulement. Pour les applications qui ont suivies, des éléments discrets ont été couplés à la MPM afin qu'un bloc rigide (DEM) puisse interagir avec un sol déformable (MPM). Cette application a consisté en l'analyse (2D) de la collision entre un bloc rocheux rigide (rond ou carré) et un sol bicouche élastoplastique. Les investigations ont été largement basées sur la mesure de coefficients de restitution (rapport des énergies cinétiques avant et après impact) qui reste difficile à déterminer expérimentalement. / In mountainous regions, the prediction of gravitational phenomena remains a challenge for the management of risk. Computational methods such as the Discrete Element Method (DEM) have been used for the modeling of these types of phenomena, where particles interact with each other to give an overall behavior of the mass. Its application can be somewhat restricted to small and medium number of blocks, since the computational time can easily become too large. Continuum analyses are therefore an attractive approach, which can reduce the computational times, but that rely on a constitutive law to represent the behavior within the mass. The main objective of this PhD was to develop a numerical tool that allowed the modeling of some specific gravitational hazards, such as the flowing of mass. A Lagrangian-Eulerian method such as the Material Point Method (MPM) is able to handle large deformations, while preserving most of the capabilities of the Finite Element Method (FEM). The method uses an Eulerian grid which is only used as a numerical scratch-pad, and remains fixed during simulations. A numerical tool named MPMbox (2D and 3D) was then developed from the ground up using C++. The code was validated using a series of analytical solutions for quasi-static analysis (some standard geotechnical tests), as well as simulations including large and rather rapid deformations (slump tests). After validation, the code was first used to make a numerical comparison with the DEM. In the comparison, a parametric survey was carried out during which the flow of a granular material on a sloped surface was simulated. Results were compared in terms of run-out distance, spread of the deposit and energy dissipated at the interface and within the mass during the flow. For a second study, discrete elements were coupled with MPM so that a rigid block could interact with a deformable soil. This application consisted in the (2D) analysis of the collision between a discrete block (round and squared) and a bounded elasto-plastic double-layered soil (soft over hard layers). The investigations were largely based on the measurement of the restitution coefficient (ratio of kinetic energies before to after the impact), which cannot be easily determined experimentally.

The boundary value problems (BVPs) of geomechanics are challenging due to the complexity in modeling soil behavior and difficulties in modeling large deformations. While traditional numerical schemes have struggled in realistically simulating geomechanical BVPs, new numerical methods –such as the material point method (MPM)–are increasingly being used to tackle these problems. However, algorithms in MPM have not yet been sufficiently developed, scrutinized, and validated. This thesis focuses on the development, verification, and validation of MPM for use in geomechanical BVPs. In particular, the thesis focuses on simulation of cone penetration tests in both controlled environments and in field conditions.<div><br></div><div>To efficiently simulate cone penetration, a silent boundary scheme, known as a cone boundary, was proposed in the generalized interpolation material point method (GIMP), a variant of MPM. The implementation of the cone boundary in GIMP was discussed, and the boundaries were validated by comparison against several benchmark problems. The cone boundaries were shown to be suitable in transmitting energy at the boundary. In addition, the implementation of traction boundaries in GIMP was analyzed. In GIMP, traction boundaries may be implemented either at the centroid of the material point, or at the edge of the material point domain. It was shown that the implementation of traction boundaries at the edge of the domain led to stress oscillations near the boundary, which were minimized when the traction boundaries were implemented at the edge of the domain.<br></div><div><br></div><div>During cone penetration, the soil near the cone-soil interface is pushed to large strains. At large strains, soils reach critical state, a state in which the soil shears at constant volume. Simulation of incompressible materials using low-order shape functions commonly used in GIMP leads to stiffer solutions and stress oscillations. To mitigate the constraints imposed by incompressibility, the non-linear B-bar method was implemented in GIMP. The modifications required for the implementation of the B-bar method in GIMP were discussed, and the efficacy of the method in mitigating incompressibility was demonstrated by analyzing several benchmark problems.<br></div><div><br></div><div>To simulate cone penetration in saturated soil, a coupled formulation was proposed in GIMP.A single material point was used to represent both the soil matrix and water. The governing equations were solved using an explicit scheme with the velocity of the soil matrix and the velocity of water as the primary variables. The formulation was validated through problems for which analytical or numerical solutions are available.<br></div><div><br></div><div>Finally, cone penetration analyses were performed both in dry sand and saturated clays. Two bounding surface models –one for sand and one for clay –were used for accurately capturing the soil response. Cone penetration tests were performed on Ottawa 20-30 sand under a variety of loading conditions at a large calibration chamber. The penetration resistances were measured, and the displacement fields were captured using the digital image correlation technique(DIC). The cone penetration resistances predicted by MPM were within 25% of the measured values, and the displacement fields computed using MPM were similar to those captured using DIC. For saturated clays, cone penetration test results reported in the literature for a Boston Blue Clay (BBC) test site were used. The simulated cone resistance of 650 kPa lied within the CPT resistance range of 580-730 kPa reported in the field. The results demonstrate the capability of MPM in simulating cone penetration in both sands and clays provided that sufficiently accurate algorithms and advanced constitutive models capable of reproducing realistic soil behavior are used in the analyses.<br></div>

The material point method (MPM) is efficiently applied for the simulation of structures undergoing large deformations where fracture and crack initiation are expected. The eigenfracture approach is introduced in the paper at hand for the implicit MPM to model crack development and propagation in static and dynamic fracture of brittle elastic materials. Eigenfracture is an energetic fracture formulation applied in the postprocessing step of the implicit MPM, making its implementation relatively straightforward. Furthermore, the driving energy used to check crack propagation is evaluated using the representative crack elements (RCE), by which the crack is modeled as a discrete phenomenon. The RCE approach shows more realistic results compared to other split models. Additionally, the fracture description of reinforced materials within the MPM is also presented in this article by coupling truss finite elements to the MPM, considering the bond stress-slip constitutive model. Two- and three-dimensional problems in static and dynamic applications are presented to assess the efficacy of the approach.

Ultra-high temperature ceramics (UHTCs) are a class of materials that maintain their structural integrity at high temperatures, e.g. 2000 °C. They have been limited in their aerospace applications because of their relatively high density and the difficulty involved in forming them into complex shapes, like leading edges and inlets. Recent advanced processing techniques have made significant headway in addressing these challenges, where the introduction of multiscale porosity has resulted in lightweight UHTCs dubbed multiscale porous UHTCs. The effect of multiscale porosity on material properties must be characterized to enable design, but doing so experimentally can be costly, especially when attempting to replicate hypersonic flight conditions for relevant testing of selected candidate samples. As such, this dissertation seeks to computationally characterize the thermomechanical properties of multiscale porous UHTCs, specifically titanium diboride, and validate those results against experimental results so as to build confidence in the model. An implicit quasi-static variant of the Material Point Method (MPM) is developed, whose capabilities include intrinsic treatment of large deformations and contact which are needed to capture the complex material behavior of the as-simulated porous UHTC microstructures. It is found that the MPM can successfully obtain the elastic thermomechanical properties of multiscale porous UHTCs over a wide range of temperatures. Furthermore, characterizations of post-elastic behavior are found to be qualitatively consistent with data obtained from uniaxial compression experiments and Brazilian disk experiments. / Doctor of Philosophy / This dissertation explores a class of materials called ultra-high temperature ceramics (UHTCs). These materials can sustain very high temperatures without degrading, and thus have the potential to be used on hypersonic aircraft which routinely experience high temperatures during flight. In lieu of performing experiments on physical UHTC specimens, one can perform a series of computer simulations to figure out how UHTCs behave under various conditions. This is done here, with a particular focus what happens when pores are introduced into UHTCs, thus rendering them more like a sponge than a solid block of material. Doing computer simulations instead of physical experiments is attractive because of the flexibility one has in a computational environment, as well as the significantly decreased cost associated with running a simulation vs. setting up and performing an experiment. This is especially true when considering challenging operating environments like those experienced by high-speed aircraft. The ultimate goal with this research is to develop a computational tool than can be used to design the ideal distribution of pores in UHTCs so that they can best perform their intended functions.

The current state of the art in structural health monitoring is primarily reliant on sensing deformation of structures at discrete locations using sensors and detecting damage using techniques such as X-ray, microCT, acoustic emission, impedance methods etc., primarily employed at specified intervals of service life. There is a need to develop materials and structures with self-sensing capabilities such that deformation and damage state can be identified in-situ real time. In the current work, the inherent deformation and damage sensing capabilities of carbon nanotube (CNT) based nanocomposites are explored starting from the nanoscale electron hopping mechanism to effective macroscale piezoresistive response through finite elements based computational micromechanics techniques. The evolution of nanoscale conductive electron hopping pathways which leads to nanocomposite piezoresistivity is studied in detail along with its evolution under applied deformations. The nanoscale piezoresistive response is used to evaluate macroscale nanocomposite response by using analytical micromechanics methods. The effective piezoresistive response, obtained in terms of macroscale effective gauge factors, is shown to predict the experimentally obtained gauge factors published in the literature within reasonable tolerance. In addition, the effect of imperfect interface between the CNTs and the polymer matrix on the mechanical and piezoresistive properties is studied using coupled electromechanical cohesive zone modeling. It is observed that the interfacial separation and damage at the nanoscale leads to a larger nanocomposite irreversible piezoresistive response under monotonic and cyclic loading because of interfacial damage accumulation. As a sample application, the CNT-polymer nanocomposites are used as a binding medium for polycrystalline energetic materials where the nanocomposite binder piezoresistivity is exploited to provide inherent deformation and damage sensing. The nanocomposite binder medium is modeled using electromechanical cohesive zones with properties obtained through the Mori-Tanaka method allowing for different local CNT volume fractions and orientations. Finally, the traditional implementation of Material Point Method (MPM) is extended for composite problems with large deformation (e.g. large strain nanocomposite sensors with elastomer matrix) allowing for interfacial discontinuities appropriately. Overall, the current work evaluates nanocomposite piezoresistivity using a multiscale modeling framework and emphasizes through a sample application that nanocomposite piezoresistivity can be exploited for inherent sensing in materials. / Ph. D.

Dans cette thèse, la Méthode des Points Matériels (MPM) est étendue à l’approximation de Galerkin Discontinue (DG) et appliquée aux problèmes hyperboliques en mécanique des solides. La méthode résultante (DGMPM) a pour objectif de suivre précisément les ondes dans des solides subissant de fortes déformations et dont les modèles constitutifs dépendent de l’histoire du chargement. A la croisée des méthodes de types éléments finis et volumes finis, la DGMPM s’appuie sur une grille de calcul arbitraire dans laquelle des flux sont calculés au moyen de solveurs de Riemann approximés sur les arêtes entre les éléments. L’intérêt de ce type de solveurs est qu’ils permettent l’introduction de la structure caractéristique des solutions des équations aux dérivées partielles hyperboliques directement dans le schéma numérique. Les analyses de stabilité et de convergence ainsi que l’illustration de la méthode sur des simulations de problèmes unidimensionnels et bidimensionnels montrent que le schéma numérique permet d’améliorer le suivi des ondes par rapport à la MPM. Par ailleurs, un deuxième objectif poursuivi dans cette thèse consiste à caractériser la réponse des solides élastoplastiques à des sollicitations dynamiques en deux dimensions en vue d’améliorer la résolution numérique de ces problèmes. Bien qu’un certain nombre de travaux aient déjà été menés dans cette direction, les problèmes étudiés se limitent à des cas particuliers. Un cadre unifié pour l’étude de la propagation d’ondes simples dans les solides élastoplastiques en déformations et contraintes plane est proposé dans cette thèse. Les trajets de chargement suivis à l’intérieur de ces ondes simples sont de plus analysés. / In this thesis, the material point method (MPM) is extended to the discontinuous Galerkin approximation (DG) and applied to hyperbolic problems in solid mechanics. The resulting method (DGMPM) aims at accurately following waves in finite-deforming solids whose constitutive models may depend on the loading history. Merging finite volumes and finite elements methods, the DGMPM takes advantage of an arbitrary computational grid in which fluxes are evaluated at element faces by means of approximate Riemann solvers. This class of solvers enables the introduction of the characteristic structure of the solutions of hyperbolic partial differential equations within the numerical scheme. Convergence and stability analyses, along with one and two-dimensional numerical simulations,demonstrate that this approach enhances the MPM ability to track waves. On the other hand, a second purpose has been followed: it consists in identifying the response of two-dimensional elastoplastic solids to dynamic step-loadings in order to improve numerical results on these problems. Although some studies investigated similar questions, only particular cases have been treated. Thus,a generic framework for the study of the propagation of simple waves in elastic-plastic solids under plane stress and plane strain problems is proposed in this thesis. The loading paths followed inside those simple waves are further analyzed.

Les matériaux à particules hautement déformables sont des formes complexes de matière avec de nombreuses applications en chimie, pharmacie, cosmétique et agro-alimentaire. L’effet conjugué du désordre et des grandes déformations des particules conduit à des propriétés mécaniques nouvelles par rapport aux matériaux à particules indéformables. En particulier, la compressibilité et la résistance au cisaillement sont contrôlées par une combinaison de réarrangements et de changement de forme des particules. Dans ce travail de thèse, nous avons développé une approche numérique originale pour la simulation de ces systèmes. Pour permettre aux particules de se déformer indéfiniment, nous avons modélisé chaque particule par un agrégat de particules primaires sans frottement qui interagissent par une force d’attraction de type Lennard-Jones et une contrainte de non-interpénétration. La dissipation d’énergie par collisions inélastiques entre les particules primaires confère un caractère plastique aux déformations des particules. Nous avons utilisé ce modèle pour étudier les propriétés de compaction et de cisaillement de ces systèmes. Nos résultats ont permis de mettre en évidence le caractère non-linéaire de la compressibilité lorsque la compacité progressivement augmente au-delà de celles des assemblages de particules indéformables. Sous cisaillement, un état critique est atteint avec une dilatance contrôlée par la pression de confinement. Dans cet état, nous avons exploré les distributions des formes des particules, les textures et les distributions des forces pour différentes valeurs de la pression. Nous avons également comparé la compressibilité simulée par l’approche développée avec celle obtenue par la Méthode de Points Matériels (MPM) en utilisant des particules élastiques. / Soft-particle materials are complex forms of matter that occur in numerous applications in chemical, pharmaceutical, cosmetic and food products. Joint effects of disorder and large particle deformations lead to novel mechanical properties that differ from those of rigid-particle materials. In particular, the compressibility and shear resistance depend on both particle rearrangements and their shape change. In this doctoral work, we developed an original approach for numerical simulation of these systems. To allow the particles to deform without breakage, each particle is modeled as an aggregate of frictionless primary particles interacting via a Lennard-Jones attraction force and impenetrability constraints. Energy dissipation by inelastic collisions between primary particles leads to the plastic nature of particle deformations. This model was used to investigate the compaction and shear behavior of soft-particle systems. We find that the compressibility is strongly nonlinear as the packing fraction increases beyond that of a random close packing of rigid particles. In continuous shearing, a critical state is reached with a dilatancy that depends on the confining pressure. In this state, we investigate the shear resistance, distributions of particle shapes, fabric properties and inter-particle forces as a function of the confining pressure. We also compare our results with those obtained by using the Material Point Method (MPM) with elastic particles.

Polymeric composites reinforced with bio-materials have advantages over composites with synthetic reinforcements. Bio-based composites use low-cost and renewable reinforcements, have nonabrasive properties for machining, have improved damping characteristics, and have potential for energy recycling. However, the limited use of bio-based composites is because their mechanical properties are typically much lower than those of synthetic composites. The objective of this study was to combine state-of-the-art imaging tools with emerging numerical modeling methods for an integrated, multi-level characterization of bio-based reinforcements and their composites. Digital photography (2D) will allow collection of full-field digital images of the surface of sample composites, which will be used for characterization of the morphological structure of fillers (copper wire or wood particle) and of model composites. Mechanical experiments (tension load) on isolated fillers and on model composites will allow imaging of the deformed material. By correlating relative positions of thousands of surface features between consecutive images, digital image correlation (DIC) algorithms can be used to map surface deformation fields and calculate surface strain fields. Digital imaging methods can only record deformations and strains. The interpretation of those strains in terms of material properties, such as position-dependent modulus of a heterogeneous composite material, requires simultaneous modeling. The modeling must use morphology-based methods that can handle anisotropy, heterogeneity, and the complex structure of bio-based composites such as wood plastic composites. This research used the material point method (MPM) as a modeling tool. MPM is a particle-based, meshless method for solving problems in computational mechanics. The crucial advantage of MPM over other methods is the relative ease of translating pixels from digital images into material points in the analysis. Thus digital images (2D) used in our experiments were used as direct input to the MPM software, so that the actual morphologies, rather than idealized geometries, were modeled. This procedure removes typical uncertainties connected with idealization of the internal features of modeled materials. It also removes variability of specimen to specimen due to morphology variations. Full-field imaging techniques and computer modeling methods for analysis of complex materials have developed independently. This research Coupled imaging and modeling and used inverse problem methodology for studying bio-particulate composites. The potential of coupling experiments with morphology-based modeling is a relatively new area. This work studied the morphology and mechanical properties of copper wire (for validation experiments) and wood particles used for reinforcement in polymer composites. The goal was to determine the in situ mechanical and interfacial properties of copper wire and then wood particles. By comparison of DIC results to MPM, the conclusion is MPM simulation works well by simulating 3D composite structure and using Matlab software to do qualitative and quantitative comparisons. Copper validation tests showed that copper wire is too stiff compared to polymer such that the inclusion modulus had low effect on the surface strains (DIC experimental results). Wood particle worked better because modulus of wood is much lower than copper. By qualitative comparison of the wood particle specimens, we could deduce that the in situ properties of wood particles are lower than bulk wood. Quantitative analysis concentrated on small area and got more exact results. In a 90 degree particle quantitative study, MPM simulations were shown to be capable of tracking the structure of wood particle plastic, which involved failure. The entire approach, however, is not very robust. We can get some results for mechanical properties, but it does not seem possible to extract all anisotropic properties from a few DIC tests, as some researcher have suggested. / Graduation date: 2012

[pt] Os problemas geotécnicos e geológicos envolvem a descrição do comportamento de materiais tais como solo e rocha e sua eventual interação com fluidos e estruturas. Em geral, a evolução desses problemas é caracterizada por grandes deformações e deslocamentos, descontinuidades, heterogeneidades e um comportamento constitutivo complexo. A abordagem deste tipo de problemas requer técnicas numéricas que levem em conta essas características, sem apresentar inconvenientes numéricos associados à distorção dos elementos como acontece no método dos elementos finitos (FEM). A presente tese desenvolve um algoritmo computacional baseado no método do ponto material (MPM) para aproximar a solução das equações governantes dos fenômenos mencionados. O algoritmo é baseado numa formulação dinâmica tridimensional do contínuo, considerando grandes deformações. Os amortecimentos de Rayleigh o local não viscoso são incorporados para modelar problemas dinâmicos e quase-estáticos. A geração dinâmica das poro-pressões é formulada assumindo o meio poroso saturado e um ponto material para discretizar a mistura. Diversas técnicas de suavização das pressões são avaliadas em problemas de impacto sobre meios saturados. Diferentes modelos constitutivos são implementados para modelar tanto a formação da superfície de ruptura e o processo de escoamento da massa de solo durante as instabilidades, quanto a gênese, evolução e quantificação de zona de falhas nos processos geológicos. Para a abordagem da discretização de problemas de escala geológica de grande número de partículas usando o MPM, uma metodologia é proposta e verificada com a discretização do deslizamento do Daguangbao, na China. Visando à diminuição do tempo computacional, o algoritmo é implementado segundo o paradigma de programação paralela. / [en] Geotechnical and geological problems involve the description of the behavior of materials such as soil and rock, and their eventual interaction with fluids and structures. In general, the evolution of these problems is characterized by large deformations and displacements, discontinuities, heterogeneities and complex constitutive behavior. Addressing these problems requires numerical techniques that take these characteristics into account, without numerical drawbacks associated with element distortion as occurs in the finite element method (FEM). In this thesis is developed a computational algorithm based on the material point method (MPM) to approximate the solution of the governing equations to the mentioned phenomena. The algorithm is based on a three-dimensional dynamic formulation of the continuum considering large deformations. Rayleigh damping and non-viscous local damping are incorporated to model dynamic and quasi-static problems. The dynamic generation of pore pressures is formulated assuming the saturated porous medium and a single material point to discretize the mixture. Different techniques are evaluated to mitigate spurious pressure in impact problems on saturated media. Different constitutive models are implemented to model the failure surface and the soil mass flow process during slope instabilities, as well as the genesis, evolution and failure zone quantification in geological processes. To address the discretization of large-scale geological problems using MPM, a methodology is proposed and validated with the discretization of the Daguangbao landslide, in China. In order to decrease the computational time, the algorithm is implemented according to the parallel programming paradigm.