Search NASA⌕ Search

SEARCH · Search NASA

Results for “triaxial compression”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

HP-TACO: A high-pressure triaxial compression apparatus for in situ x-ray measurements in geomaterials

Triaxial compression experiments are commonly used to characterize the elastic and inelastic behavior of geomaterials. In situ measurements of grain kinematics, particle breakage, stresses, and other microscopic phenomena have seldom been made during such experiments, particularly at high pressures relevant to many geologic and man-made processes, limiting our fundamental understanding. To address this issue, we developed a new triaxial compression device called HP-TACO (High-Pressure TriAxial COmpression Apparatus). HP-TACO is a miniaturized, conventional triaxial compression apparatus permitting confining pressures up to 50 MPa and deviatoric straining of materials, while also allowing in situ x-ray measurements of grain-scale kinematics and stresses. Here, in this study, we present the design of and first results from HP-TACO during its use in laboratory and synchrotron settings to study grain-scale kinematics and stresses in triaxially compressed sands subjected to 15 and 30 MPa confining pressures. The data highlight the unique capabilities of HP-TACO for studying the high-pressure mechanics of sands, providing new insight into micromechanical processes occurring during geologic and man-made processes.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Volumetric and shear strain localization throughout triaxial compression experiments on rocks

Deformation localization is a widely observed, but rarely quantified process in the crust. Recent observations suggest that the localization of seismicity and fracture networks can help identify the approach to catastrophic failure. Here, we quantify the localization processes of the volumetric and deviatoric strain components in twelve triaxial compression experiments imaged with X-ray tomography. We capture three-dimensional images of the rock cores during triaxial compressing toward failure, and then calculate the local strain components using digital volume correlation. The divergence and curl of the incremental displacement vector field provide the volumetric and deviatoric components of the strain field. We quantify localization using the proportion of the rock occupied by high magnitudes of the volumetric and deviatoric strains, and the Gini coefficient of these high magnitude strains, which measures the deviation from a uniform process. We find that the vast majority, but not all, of the experiments experience strain localization toward failure. The rocks typically experience their maximum degree of strain localization not immediately preceding failure, but on average at 90% of the failure stress. The volumetric strain tends to localize more than the deviatoric strain. Finally, these observations support using the localization of the volumetric strain, along with the deviatoric strain, to identify the evolution of the precursory phase preceding earthquakes.

58 GEOSCIENCES↗
DOE OSTI2024

Intermediate Principal Stress Effects on the 3D Cracking Behavior of Flawed Rocks Under True Triaxial Compression

Crack initiation, growth, and coalescence in flawed rocks have been extensively studied for 2D (planar, penetrating) flaws under uniaxial/biaxial compression. However, little is known as to the mechanisms and processes of cracking from 3D flaws under true triaxial compression, where the intermediate principal stress (σ 2 ) is distinguished from the major and minor principal stresses. In this work, we systematically investigate the effects of σ 2 on the 3D cracking behavior of rock specimens with preexisting flaws, through the use of mechanistic simulations of mixed-mode fracture in rocks. We explore how two characteristics of σ 2 , namely, (i) its orientation with respect to the flaw and (ii) its magnitude, affect two aspects of the cracking behavior, namely, (i) the cracking pattern and (ii) the peak stress. Results show that the orientation of σ 2 exerts more control over the cracking pattern than the flaw inclination angle. The peak stress becomes highest when σ 1 is parallel to the flaw, whereas it becomes lowest when σ 2 is parallel to the flaw. Also, the effects of σ 2 magnitude are more significant when σ 2 becomes more oblique to the flaw plane. On the basis of our observations, we propose mechanisms underlying the cracking behavior of 3D flawed rocks under true triaxial compression.

58 GEOSCIENCES↗
DOE OSTI2021

Fracturing process of micro-concrete under uniaxial and triaxial compression: Insights from in-situ X-ray mechanical tests

This paper presents an experimental study of concrete at meso-scale (aggregates, macro-pores and mortar matrix) in order to get a better understanding of the local failure mechanisms known to drive the macroscopic mechanical behaviour of the material. The main originality comes from conducting in-situ X-ray mechanical tests on micro-concrete samples of realistic composition (including cement, sand, aggregates and water), under uniaxial compression and, for the first time, under triaxial compression at 5 MPa, 10 MPa and 15 MPa confining pressures. A timeseries analysis of the set of 3D images coming from each test allows for the measurement of the 3D kinematic fields (displacement and strain fields) throughout the experiments. The different failure patterns observed for each loading path are discussed, along with a quantification of the 3D fracturing processes at the scale of the largest heterogeneities (aggregates and macro-pores). With an increasing level of confinement, the transition from brittle to ductile response is observed, as well as an increase of the strength of the material. The pronounced impact of the meso-scale heterogeneities of concrete on their local failure mechanisms is highlighted. It is shown that strain localisation mainly originates between aggregates and mortar matrix, with the shape and location of the largest aggregates and macro-pores essentially driving the propagation of the cracking network.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Fracture Network Localization Preceding Catastrophic Failure in Triaxial Compression Experiments on Rocks

We quantify the spatial distribution of fracture networks throughout six in situ X-ray tomography triaxial compression experiments on crystalline rocks at confining stresses of 5–35 MPa in order to quantify how fracture development controls the final macroscopic failure of the rock, a process analogous to those that control geohazards such as earthquakes and landslides. Tracking the proportion of the cumulative volume of fractures with volumes >90th percentile to the total fracture volume, ∑ v 90 / v t o t indicates that the fracture networks tend to increase in localization toward these largest fractures for up to 80% of the applied differential stress. The evolution of this metric also matches the evolution of the Gini coefficient, which measures the deviation of a population from uniformity. These results are consistent with observations of localizing low magnitude seismicity before large earthquakes in southern California. In both this analysis and the present work, phases of delocalization interrupt the general increase in localization preceding catastrophic failure, indicating that delocalization does not necessarily indicate a reduction of seismic hazard. However, the proportion of the maximum fracture volume to the total fracture volume does not increase monotonically. Experiments with higher confining stress tend to experience greater localization. To further quantify localization, we compare the geometry of the largest fractures, with volumes >90th percentile, to the best fit plane through these fractures immediately preceding failure. The r 2 scores and the mean distance of the fractures to the plane indicate greater localization in monzonite than in granite. The smaller mean mineral diameter and lower confining stress in the granite experiments may contribute to this result. Tracking these various metrics of localization reveals a close association between macroscopic yielding and the acceleration of fracture network localization. Near yielding, ∑ v 90 / v t o t and the Gini coefficient increase while the mean distance to the final failure plane decreases. Macroscopic yielding thus occurs when the rate of fracture network localization increases.

58 GEOSCIENCES↗
DOE OSTI2021

Characterization of damage processes in Montney siltstone under triaxial compression using acoustic emission and diagnostic imaging

SUMMARY Crack nucleation and rock failure processes in a fine-grained siltstone (Montney Formation) under triaxial compression are investigated using combined diagnostic techniques, including ultrasonic-wave measurement, acoustic-emission (AE) monitoring, computed tomography (CT) scanning, and thin-section imaging. The sample displays a weak-to-moderate inherent seismic anisotropy and noticeable stress-induced anisotropy prior to failure. No AE event was detected until the applied axial stress reached 95 per cent of the peak value. The signal-to-noise ratio is relatively low, however, and detectable AE events are more diffuse than those observed in highly brittle rocks. The AE locations correlate with a shear fracture zone imaged by CT scanning. AE moment-tensor analysis reveals that events with larger relative magnitudes are characterized by high volumetric (tensile or compressive) components, and the initiation of the failure zone is dominated by combined shear-tensile (ST) failure. Stress inversion of the AE events with high tensile components is in good agreement with the known applied stress. Microscopic imaging of thin sections from the failed sample shows that the failure zone is an en echelon structure consisting of a major fracture with branching micro and minor cracks. This failure mechanism is consistent with a ST source mechanism and is interpreted to be associated with the fine granular structure and mineral composition of Montney siltstone.

Jia, Suzie Q. (ORCID:0000000326496190)↗
DOE OSTI2023

Assessing continuum plasticity postulates with grain stress and local strain measurements in triaxially compressed sand

Critical state and continuum plasticity theories have been used in research and engineering practice in soil and rock mechanics for decades. These theories rely on postulated relationships between material stresses and strains. Some classical postulates include coaxiality between stress and strain rates, stress–dilatancy relationships, and kinematic assumptions in shear bands. Although numerical and experimental data have quantified the strains and grain kinematics in such experiments, little data quantifying grain stresses are available. Here, we report the first-known grain stress and local strain measurements in triaxial compression tests on synthetic quartz sands using synchrotron X-ray tomography and 3D X-ray diffraction. We use these data to examine the micromechanics of shear banding, with a focus on coaxiality, stress-dilatancy, and kinematics within bands. Our results indicate the following: 1) elevated deviatoric stress, strain, and stress ratios in shear bands throughout experiments; 2) coaxial principal compressive stresses and strains throughout samples; 3) significant contraction along shear bands; 4) vanishing volumetric strain but nonvanishing stress fluctuations throughout samples at all stages of deformation. Our results provide some of the first-known in situ stress and strain measurements able to aid in critically evaluating postulates employed in continuum plasticity and strain localization theories for sands.

42 ENGINEERING↗
DOE OSTI2024

Geophysical Signatures of Crack Network Coalescence in Rocks at Multiple Length Scales

The main goal of the research project was to identify the geophysical signatures of fracture growth in natural rocks by utilizing novel geophysical techniques. The research objectives were to (a) investigate the potential for geophysical methods to determine when cracks initiate, the types and locations of propagated cracks, and the coalescence of networks of cracks in natural rocks at multiple scales, (b) determine how damage at the microscale evolved into damage at the macroscale and then link the microscopic and macroscopic observations, (c) quantify crack coalescence in rocks under realistic stress conditions using coupled mechanical-geophysical-optical visualization, and (d) identify the precursors in geophysical signals to crack coalescence. The following research thrusts were explored to achieve the research objectives: (1) uniaxial compression testing of rock specimens with and without a set of pre‐existing flaws and (2) triaxial compression testing of natural rock specimens. These thrusts allowed for exploring fracturing in rocks under realistic in situ environments and at multiple scales. This project provided educational opportunities for nine graduate and undergraduate students and resulted in 27 peer-reviewed publications. This first research thrust focused on investigating the micromechanics of fractures in rocks through uniaxial compression testing combined with advanced geophysical and imaging techniques, specifically acoustic emission (AE) monitoring, ultrasonic imaging, and 2-dimensional Digital Image Correlation (2D-DIC). By examining damage processes under time-independent and time-dependent loading conditions, insights into damage localization, crack initiation, and fracturing mechanisms were gained. It was observed that the AE signals and the strain-based measurements directly reflect the state of damage in the rock specimen and could be used to identify the cracking levels, such as the crack initiation (CI) and crack damage (CD), and the mode of deformation. A novel calibration apparatus was developed to enhance the accuracy of AE sensors, allowing for the estimation of key parameters such as magnitude, source dimension, stress drop, and radiated seismic energy associated with the fractures. The findings highlighted significant variations in the temporal evolution of AE source parameters during the primary, secondary, and tertiary stages of creep, identifying tensile cracking as the primary deformation mode. The second research thrust focused on enhancing the understanding of fracturing processes in natural rocks through triaxial compression testing, real-time AE monitoring, and ultrasonic monitoring. We investigated the impact of various factors such as fracture propagation regimes, injection parameters, rock types, and pre-existing conditions on the hydraulic fracture (HF) behavior using scaled true-triaxially loaded specimens of Barre granite and Lyons sandstone. Custom sensor housing facilitated concurrent active and passive monitoring to analyze hydro-mechanical responses and microseismicity associated with different HF scenarios. A coupled investigation of passive microseismicity and active signal attributes permitted a detailed comprehension of the various HF processes (aseismic deformation, fracture initiation and propagation, fluid permeation, and leak-off) and their dependence on the specific rock type. The findings of this research demonstrated the effectiveness of AE monitoring techniques in providing valuable insights into the impact of various factors on the behavior and dynamics of HF processes. The advancements in monitoring techniques, offering a more thorough and precise approach, represent a significant step towards optimizing HF practices and ensuring sustainable resource extraction.

58 GEOSCIENCES↗
DOE OSTI2023

Exploring pressure-dependent inelastic deformation and failure in bonded granular composites: An energetic materials perspective

In polymer-filled granular composites, damage may develop in mechanical loading prior to material failure. Damage mechanisms such as microcracking or plastic deformation in the binder phase can substantially alter the material’s mesostructure. For energetic materials, such as solid propellants and plastic bonded explosives, these mesostructural changes can have far reaching effects including degraded mechanical properties, potentially increased sensitivity to further insults, and changes in expected performance. Unfortunately, predicting damage is nontrivial due to the complex nature of these composites and the entangled interactions between inelastic mechanisms. In this work, we assess the current literature of experimental knowledge, focusing on the pressure-dependent shear response, and propose a simple simulation framework of bonded particles to study four limiting-case material formulations at both meso- and macro-scales. Further, to construct the four cases, we systematically vary the relative interfacial strength between the polymer binder and granular filler phase and also vary the polymer’s glass transition temperature relative to operating temperature which determines how much the binder can plastically deform. These simulations identify key trends in global mechanical response, such as the emergence of strain hardening or softening regimes with increasing pressure which qualitatively resemble experimental results. By quantifying the activation of different inelastic mechanisms, such as bonds breaking and plastically straining, we identify when each mechanism becomes relevant and provide insight into potential origins for changes in mechanical responses. The locations of broken bonds are also used to define larger, mesoscopic cracks to test various metrics of damage. We primarily focus on triaxial compression, but also test the opposite case of triaxial extension to highlight the impact of Lode angle on mechanical behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Predicting Fracture Network Development in Crystalline Rocks

The geometric properties of fractures influence whether they propagate, arrest, or coalesce with other fractures. Thus, quantifying the relationship between fracture network characteristics may help predict fracture network development, and perhaps precursors to catastrophic failure. To constrain the relationship and predictability of fracture characteristics, we deform eight one centimeter tall rock cores under triaxial compression while acquiring in situ X-ray tomograms. The tomograms reveal precise measurements of the fracture network characteristics above the spatial resolution of 6.5 µm. We develop machine learning models to predict the value of each characteristic using the other characteristics, and excluding the macroscopic stress or strain imposed on the rock. The models predict fracture development more accurately in the experiments performed on granite and monzonite, than the experiments on marble. Fracture network development may be more predictable in these igneous rocks because their microstructure is more mechanically homogeneous than the marble, producing more systematic fracture development that is not strongly impeded by grain contacts and cleavage planes. The varying performance of the models suggest that fracture volume, length, and aperture are the most predictable of the characteristics, while fracture orientation is the least predictable. Orientation does not correlate with length, as suggested by the idea that the orientation evolves with increasing differential stress and thus fracture length. This difference between the observed and expected relationship between orientation and length highlights the influence of mechanical heterogeneities and local stress perturbations on fracture growth as fractures propagate, link, and coalesce.

58 GEOSCIENCES↗
DOE OSTI2025

Influence of gradation on shear behavior and pore water response of sand

Here, this paper examines the influence of particle size distribution (PSD) on the shearing behavior of sand using axisymmetric triaxial compression (ATC) testing in conjunction with 3D imaging using synchrotron micro-computed tomography (SMT). Angular sand was sieved to obtain a narrow gradation (uniform) and a wide gradation (non-uniform) with uniformity coefficients (C u ) of 1.2 and 3.0, respectively. Experiments were conducted under drained conditions and a constant effective cell pressure (σ$^{'}_3$) of 50 kPa at low and high back pressures (BP). Results of paired ATC-SMT were analyzed in terms of global constitutive behavior, strain localization, and localized flow characteristics. A new particle tracking algorithm is developed to expose strain localization. Relative particle translation gradient (RPTG) results revealed that gradation affects how early the shear band develops and a tendency for narrow gradation to exhibit varying degrees of diffused shear band at high axial strain (ε 1 ) even though a clear single shear band was observed at the peak strength. Geometric tortuosity (τ) is used as a quantitative parameter to assess flow changes and the efficiency of the flow path. Representative elementary volumes (REVs) inside and outside the shear band were analyzed and the results showed a variation in flow path characteristics where wide gradation specimens had preferential flow paths.

Computed tomography↗
DOE OSTI2022

Critical review of the Mogi failure criterion based on true-triaxial laboratory data analysis and theoretical considerations

The influence of the intermediate principal stress, σ 2 , is disregarded when rock strength is analyzed using the Mohr-Coulomb or Hoek-Brown criteria. Several polyaxial strength criteria have been proposed to account for the influence of σ 2 . One of the methods that received considerable attention is the Mogi criterion. The criterion modifies the Nadai or Drucker-Prager approaches by exchanging the normal octahedral stress with what is introduced as the mean effective stress, σ m,2 . This change is based on experimental evidence that the influence of the intermediate principal stress on strength is less pronounced than the influence of the minimum principal stress. However, the Mogi criterion was derived somewhat empirically and consequently it is essential to analyze its potential theoretical deficiencies. In the current note, we analyzed the Mogi failure criterion, describing its characteristics in both τ oct –σ m,2 and σ 1 -σ 2 spaces. We investigated the relation between Mohr-Coulomb and Mogi criteria, which marks the basis for the commonly used Mogi-Coulomb approach. Here, we compared the performance of the Mogi criterion against available polyaxial strength data using three fitting approaches: Mogi-Coulomb based on triaxial compression results only, linear Mogi based on all polyaxial strength results, and power law Mogi also based on all data. The analysis reveals several serious deficiencies of the Mogi approach which all suggest that the Mogi criterion cannot properly represent the σ 2 dependence of rock strength in many cases.

58 GEOSCIENCES↗
DOE OSTI2021

Evolution of anisotropy with saturation and its implications for the elastoplastic responses of clay rocks

Many clay rocks have distinct bedding planes. Experimental studies have shown that their mechanical properties evolve with the degree of saturation (DOS), often with higher stiffness and strength after drying. For transversely isotropic rocks, the effects of saturation can differ between the bed-normal (BN) and bed-parallel (BP) directions, which gives rise to saturation-dependent stiffness and strength anisotropy. Accurate prediction of the mechanical behavior of clay rocks under partially saturated conditions requires numerical models that can capture the evolving elastic and plastic anisotropy with DOS. In this study, we present an anisotropy framework for coupled solid deformation-fluid flow in unsaturated elastoplastic media. Here we incorporate saturation-dependent strength anisotropy into an anisotropic modified Cam-Clay (MCC) model and consider the evolving anisotropy in both the elastic and plastic responses. The model was calibrated using experimental data from triaxial tests to demonstrate its capability in capturing strength anisotropy at various levels of saturation. Through numerical simulations, we demonstrate the role of evolving stiffness and strength anisotropy in the mechanical behavior of clay rocks. Plane strain simulations of triaxial compression tests were also conducted to demonstrate the impacts of material anisotropy and DOS on the mechanical and fluid flow responses.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Graphene nanoplatelets reinforced cement as a solution to leaky wellbores reinforcing weak points in hydrated Portland cement with graphene nanoparticles improves mechanical and chemical durability of wellbore cements

In this work, to improve the performance of wellbore cement under subsurface environments, graphene nanoplatelets (GNPs) were added in various percentages to Class-H cement slurry. Microstructural characterization of cement slurries cured at 90 °C and 95% RH indicates that GNP modifies the microstructure of hydrated cement by reinforcing pore spaces. As a result, the mechanical properties, such as Young's modulus and axial peak stress, obtained from high temperature triaxial compression tests are significantly improved based on different percentages of added GNPs. Furthermore, the hydrated 1 × 2 inch GNP- Portland cement cores tested under simulated deep wellbore conditions of high-temperature and high pressure (HTHP), appear to have a ductile-like behavior, when compared to a typical brittle nature of Portland cement pastes. From our observations and published data on graphene resistance to fracture, GNP addition to cement is evidently enhancing the flexibility of cement. These improvements would reduce the risks associated with wellbore cement deterioration and a consequent leakage in fossil fuel production, geothermal energy production, CO 2 storage as well as long-term sealing materials in plugging and abandonment of all wellbores at the end of their service life.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hyperelastic nature of the Hoek–Brown criterion

In this article, we propose a nonlinear elasto-plastic model, for which a specific class of hyperbolic elasticity arises as a straight consequence of the yield criterion invariance on the plasticity level. We superimpose this nonlinear elastic (or hyperelastic) behavior with plasticity obeying the associated flow rule. Interestingly, we find that a linear yield criterion on the thermodynamical force associated with plasticity results in a quadratic yield criterion in the stress space. This suggests a specific hyperelastic connection between Mohr–Coulomb and Hoek–Brown (or alternatively between Drucker–Prager and Pan–Hudson) yield criteria. We compare the elasto-plastic responses of standard tests for the Drucker–Prager yield criterion using either linear or the suggested hyperbolic elasticity. Notably, the nonlinear case stands out due to dilatancy saturation observed during cyclic loading in the triaxial compression test. We conclude this study with structural finite element simulations that clearly demonstrate the numerical applicability of the proposed model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Concurrent measurement of strain and chemical reaction rates in a calcite grain pack undergoing pressure solution: Evidence for surface-reaction controlled dissolution

Pressure solution is inferred to be a significant contributor to sediment compaction and lithification, especially in carbonate sediments. For a sediment deforming primarily by pressure solution, the compaction rate should be directly related to the rate of calcite dissolution, transport along grain contacts, and calcite reprecipitation. Previous experimental work has shown that there is evidence that deformation in wet calcite grain packs is consistent with control by pressure solution, but considerable ambiguity remains regarding the rate limiting mechanism. We present the results of laboratory compaction experiments designed to directly measure calcite dissolution and precipitation rates (recrystallization rates) concurrently with strain rate to test whether measured rates are consistent with predicted rates both in absolute magnitude and time evolution. Recrystallization rates are measured using trace element chemistry (Sr/Ca, Mg/Ca) and isotopes (87Sr/86Sr) of fluids flowing slowly through a compacting grain pack as it is being triaxially compressed. Imaging techniques are used to characterize the grain contacts and strain effects in the post-experiment grain pack. Our data show that calcite recrystallization rates calculated from all three geochemical parameters are in approximate agreement and that the rates closely track strain rate. The geochemically inferred rates are close to predicted rates in absolute magnitude. Uncertainty in grain contact dimensions makes distinguishing between surface reaction control and diffusion control difficult. Measured reaction rates decrease faster than predicted from standard pressure solution creep flow laws. This inconsistency may indicate that calcite dissolution rates at grain contacts are more complex, and more time-dependent, than suggested by geometric models designed to predict grain contact stresses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Influence of gradation on failure mode of saturated sand at particle scale

Extensive research has been reported in the literature to characterize the failure mode of dry sand using various experimental techniques such as surface optical imaging, photo-elastic materials, three-dimensional (3D) computed tomography (CT), and 3D synchrotron micro-computed tomography (SMT). However, there is a limited literature about the behavior of saturated sand. This paper presents the results of axisymmetric triaxial compression (ATC) experiments that were conducted on saturated sand specimens. The behavior of specimens composed of a uniform sand with grain size between US sieves #40 and #50 is compared to specimens conducted on the same sand that has a wider gradation. 3D SMT technique was used to acquire 3D scans while shearing the specimens to probe localized events that are completely missed or misinterpreted when analyzing ATC measurements based on global standard measurements. The results show a higher effective principal stress ratio (EPSR) for the non-uniform specimen and a thicker shear band when compared to uniform specimen.

Elnur, Mohammed [University of Tennessee]↗
DOE OSTI2025

Machine-Learning-Based Multiscale Methods for 3D Modelling of Granular Materials by Incorporating History-Dependent State Variables

Over the past decades, the prevalence of machine learning (ML) methods has made the development of ML-based constitutive models for granular materials undoubtedly a popular subject. Numerous studies have been made to feature the loading path or history-dependent stress-strain response of granular media using neural networks. In this work, a novel finite element method (FEM)–ML multiscale approach was developed by incorporating internal variables to improve the simulation accuracy of 3D history-dependent granular materials for the first time. To this end, a surrogate constitutive model based on the single-step-based multi-layer perceptron (MLP) neural network was used to replace representative volume element (RVE) simulations conducted by the discrete element method (DEM) in the multiscale FEM–DEM approach. Although the prediction principle of the MLP aligns with the FEM algorithm, artificially added internal variables are required to differentiate the loading history. To address this issue, history variables associated with the Frobenius norm are proposed to be fed into the MLP coupled with the strain tensor to extract the history-dependent behaviour of granular assemblies. The developed FEM–ML approach was demonstrated in 3D conventional triaxial compression (CTC) simulations. Compared to the multiscale FEM–DEM approach, the proposed FEM–ML method exhibits a significantly improved computational efficiency.

granular materials↗