Earthquake-induced soil liquefaction is a phenomenon in which saturated, sandy soil loses its strength and stiffness during earthquake shaking. Liquefaction can be extremely costly and damaging to infrastructure. The commonly used "simplified" stress-based liquefaction triggering framework is correlated with metrics computed from in-situ tests like the Cone Penetration Test (CPT). While CPT-based procedures have been shown to accurately predict liquefaction occurrence in homogenous, sandy soil profiles, they tend to over-predict the occurrence of liquefaction in challenging, highly interlayered soil profiles. One contributing factor to the over-prediction is multiple thin-layer effects in CPT data, a phenomenon in which data in interlayered zones is blurred or averaged, making it difficult to identify specific layer boundaries and associated CPT parameters like tip resistance. Multiple thin-layer correction procedures have been proposed to convert the measured tip resistance in an interlayered profile (qm) to the "true" or characteristic tip resistance (qt) that would be measured without the influence of multiple thin-layer effects. In this dissertation, the efficacy of existing multiple thin-layer correction procedures is assessed. It is shown that existing procedures are not effective for layer thicknesses equal to or less than about 1.6 times the diameter of the cone. Two new multiple thin-layer correction procedures are proposed. Furthermore, a framework for numerically simulating CPTs in interlayered soil profiles using the Material Point Method (MPM) is developed. A framework for linking uncertainties associated with the numerical analyses and the laboratory CPT calibration chamber tests used to calibrate the numerical analyses is also proposed. Finally, a database of laboratory and numerically-generated CPT data is presented. It is shown how this database can be used to improve existing, and develop new, multiple thin-layer correction procedures. Ultimately, the work detailed in this dissertation will improve the characterization of highly interlayered soil profiles using CPTs to support more accurate liquefaction hazard assessment at challenging soil sites. / Doctor of Philosophy / Earthquake-induced soil liquefaction is a phenomenon in which saturated, sandy soil loses its strength and stiffness during earthquake shaking. Liquefaction can be extremely costly and damaging to infrastructure. Existing procedures used to assess liquefaction hazard were developed specifically for homogenous, sandy soil profiles. These procedures do not perform well in challenging, highly interlayered soil profiles. One reason for this is the inadequate characterization of the soil profile by the chosen in-situ test method. For example, the cone penetration test (CPT) consists of hydraulically advancing a steel probe with a conical shaped tip ("cone") into the ground. Typically, the penetrometer is about 3.6 to 4.4 cm in diameter, and data are recorded at 1 to 5 cm depth intervals. However, data recorded at a specific depth are representative of soil that falls within a zone several times the diameter of the penetrometer ahead of and behind the tip of the cone. In a highly interlayered soil profile, this means the CPT records blurred or averaged data within interlayered zones. Typical liquefaction analyses compute a factor of safety against liquefaction at every depth in the soil profile where CPT data are recorded. Hence, having data that are blurred can result in an inaccurate factor of safety against liquefaction. To account for this blurring (called multiple thin-layer effects), correction procedures have been proposed. This dissertation evaluates the effectiveness of those procedures and develops new procedures. Additionally, a numerical simulation tool is shown to be capable of simulating CPTs in layered soil profiles. This reduces the need for costly laboratory testing to further evaluate multiple thin-layer effects. Finally, a combined laboratory and numerically-generated CPT database is developed to support the improvement of, and development of new, multiple thin-layer correction procedures. The broader impacts of this work support more accurate liquefaction evaluations in challenging soil profiles worldwide, like those in Christchurch, New Zealand, and the Groningen region of the Netherlands.

Free-Fall Penetrometer (FFP) testing is an easy and rapid test procedure for seabed sediment characterization favorable to conventional geotechnical testing mainly due to its cost-effectiveness. Yet, FFP testing results are interpreted using empirical correlations, but difficulties arise to understand soil behavior under the high-strain rate (HSR) loading effects during rapid FFP penetration. The numerical simulation of FFP-soil interaction is also challenging. This study aims to numerically analyze FFP testing of saturated loose sands using the particle-based Material Point Method (MPM). The numerical analysis was conducted by simulating calibration chamber FFP tests on saturated loose quartz sand. The numerical results using quasi-static properties resulted in a reaction of the sand softer than the actual calibration chamber test. This implied the necessity of considering HSR effects. After performing parametric analyses, it was concluded that dilation plays an important role in the response of sand-water mixtures. Comparison of dry and saturated simulations showed that FFP penetration increases when the soil is dry and tends to develop a general bearing capacity failure mechanism. This is because the pore water increases the stiffness of the system and due to the increased strength that develops in saturated dilative sands when negative pore pressures develop. Local bearing failure mechanism is observed in all saturated simulations. Finally, numerical CPT (quasi-static) and FFP tests were used to examine the strain rate coefficient used in practice (K); and a consistent range between 1 to 1.5 was obtained. / Master of Science / Accurate characterization of seabed sediments is crucial to understand sediment mobilization processes and to solve nearshore engineering problems such as scouring around offshore structures. Its portability, low testing effort, and repeatability make FreeFall Penetrometer (FFP) testing a highly cost-effective sediment characterization test. Nevertheless, due to the complex penetration mechanism of FFPs in soils (e.g., high-strain rate effects due to rapid FFP loading), converting FFP output into practical information is complicated, and it heavily relies on empirical correlations. This thesis presents a numerical analysis of FFP testing on saturated sand using the Material Point Method. First, the simulation results were compared with laboratory tests. Later, a parametric study was performed to understand the effect of different material parameters on the FFP response and to highlight in a simplified manner the effects of rapid loading on the sand behavior. Additional simulations in dry sand (without water) revealed that dry conditions provide larger FFP penetrations than saturated ones for the same material parameters. Lastly, the strain rate coefficient, which is a parameter required in one of the most common empirical methods for converting FFP output into geotechnical parameters, was back-calculated. The results were consistent with values used in practice for similar conditions.

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.