Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “COMPRESSION LOADING”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Green co-solvent-assisted one-pot synthesis of high-performance flexible lignin polyurethane foam

Although lignin is a promising candidate as a source of polyol for producing eco-friendly polyurethane foam (PUF), its direct application has limitations, primarily stemming from the brittleness of the resulting products and poor processibility at high lignin loadings. To address these technical challenges, we introduced an effective co-solvent for the lignin-based polyurethane foams (LPUFs) with up to 80wt% of lignin substitution for petroleum polyols. As a green co-solvent that involves carbon dioxide fixation during its production, propylene carbonate (PC), was blended with lignin and improved the distribution of lignin by decreasing the viscosity, providing a room-temperature reaction. The LPUFs demonstrated a comparable thermal insulation performance (R-values greater than 5 in -1 ) with commercially available PUF with a uniform size of closed-cell structure induced by the reduced viscosity of lignin-based polyol. LPUFs also showed enhanced fire retardancy and moisture resistance compared to the control PUFs, benefiting from lignin’s intrinsic fire-resistant and hydrophobic properties. To mitigate the brittleness issue at high lignin loadings, flexible polyethylene glycols (PEGs) with a high molecular weight (i.e., PEG-1000) were dissolved in PC at room temperature to further improve the mechanical strength of LPUFs. The compressive strength of LPUFs with 80 wt% lignin substitution content was remarkably improved to 226 kPa by the employment of co-solvent and PEG-1000 compared to the control PUF (<140 kPa). Moreover, the flexibility of the LPUFs was successfully controlled by simply adjusting the PEG-400/PEG-1000 ratio. In conclusion, the results of this study provide valuable insights for accelerating and broadening the utilization of lignin in sustainable materials and manufacturing areas.

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

A single-zone zero-dimensional study of HCCI combustion of methanol dehydration products to enable ignition of direct-injected methanol

Methanol is an alternative fuel gaining traction in the maritime sector. Its direct adoption, however, is accompanied by a unique set of technical challenges, such as low cetane number and high latent heat of vaporization. An approach to overcome these challenges is being developed at the US Department of Energy’s Oak Ridge National Laboratory, where onboard generation of dimethyl ether (DME) via catalytic dehydration of methanol can be used to assist in the mixing controlled combustion of direct-injected (DI) methanol. The generated mixture from this dehydration process can be premixed with intake air to condition the cylinder via. homogeneous charge compression ignition (HCCI) for subsequent DI methanol. In this preliminary work, various catalyst or reactor conversion efficiencies were simulated (using bottles) at constant DME and water flow at low load on a single-cylinder marine-variant of a CAT® C18 18 L engine with a 145 mm bore. To substantiate the experimental findings, a zero-dimensional engine model was developed in Cantera using a DME mechanism with 79 species and 658 reactions. Results presented include experimental and simulation heat release rate comparisons, species evolution information, and constant volume ignition delay (ID) for DI methanol with and without background species from HCCI of the premixed products from different reactor efficiencies. The results suggest that thermal effects dominate the DI methanol ignition process, and this work provides a chemical kinetic foundation or guideline for developing future control schemes.

Tyrewala, Daanish [ORNL] (ORCID:0000000208599324)↗
DOE OSTI2024

Full-field quantitative visualization of shock-driven pore collapse and failure modes in PMMA

The dynamic collapse of pores under shock loading is thought to be directly related to hot spot generation and material failure, which is critical to the performance of porous energetic and structural materials. However, the shock compression response of porous materials at the local, individual pore scale is not well understood. This study examines, quantitatively, the collapse phenomenon of a single spherical void in PMMA at shock stresses ranging from 0.4 to 1.0 GPa. Using a newly developed internal digital image correlation technique in conjunction with plate impact experiments, full-field quantitative deformation measurements are conducted in the material surrounding the collapsing pore for the first time. The experimental results reveal two failure mode transitions as shock stress is increased: (i) the first in situ evidence of shear localization via adiabatic shear banding and (ii) dynamic fracture initiation at the pore surface. Numerical simulations using thermo-viscoplastic dynamic finite element analysis provide insights into the formation of adiabatic shear bands (ASBs) and stresses at which failure mode transitions occur. Further numerical and theoretical modeling indicates the dynamic fracture to occur along the weakened material inside an adiabatic shear band. Finally, analysis of the evolution of pore asymmetry and models for ASB spacing elucidate the mechanisms for the shear band initiation sites, and elastostatic theory explains the experimentally observed ASB and fracture paths based on the directions of maximum shear.

42 ENGINEERING↗
DOE OSTI2025

Electro-Thermo-Mechanical (ETM) Study on Submarine Dynamic Power Cables (SDPC) for Offshore Wind Transmission

Understanding the failure mechanisms of submarine dynamic power cables (SDPC) is critical for innovative design to meet the 2035 cost reduction target of U.S. DOE Floating Offshore Wind Shot. This is important because the current design suffers a significant failure rate in the field. This project conducted a systematic electro-thermo-mechanical (ETM) experimental study on the power cores extracted from a 15kV power cable with three cores of copper conductor, ethylene propylene rubber (EPR) insulation, and continuously corrugated welded aluminum armor (CCWA). An ETM testing system was developed by integrating a high voltage (HV) amplifier, two ceramic heaters, and a rod-plate transverse compression setup into a mechanical testing machine. Increasing temperature from room temperature (RT, 22 degree C) to 90 degree C resulted in the 67% decrease in the failure mechanical load as defined by the dielectric breakdown. Under creep mechanical loading, the dielectric breakdown time was decreased by 70% for a given mechanical load when the specimen temperature increased from RT to 90oC. The failure strain in both monotonic and creep mechanical loading modes was related to the maximum mechanical load applied. Although the dielectric breakdown occurred in the ETM test, the post-test measurement revealed an impressive recovery of electric resistance. The failure analysis on cross section of tested specimens indicated a sizable gap across insulation layer near the area between copper conductor and loading rod, which apparently resulted from online dielectric failure and offline elastic recovery of components including conductor and insulation layer.

17 WIND ENERGY↗
DOE OSTI2026

Siloxane-modified polycarbosilane flexible Prepregs for fabrication of ceramic matrix composites via compression molding and PIP densification

For this work, the development of ceramic matrix composites (CMCs) using 5 harness satin carbon fiber fabric impregnated with commercial polycarbosilane precursor plasticized by siloxane copolymer was investigated aiming to improve wetting behavior and shape conformability. The polycarbosilane precursor and polysiloxane plasticizer were mixed at varying weight ratios to obtain flexible preceramic resin prepreg cloths. 13 plies of preceramic polymer prepreg cloths were stacked in 0/90° layup and cured by compression molding, followed by densified via PIP process. The CMCs were densified with eight PIP cycles followed by final crystallization at 1600 °C. The CMCs made with 10 wt% polysiloxane loading showed lower viscosity of the prepreg resin as well as higher strength and displacement to failure compared to those fabricated from unplasticized polycarbosilane prepolymer. In contrast, at higher concentration of plasticizer, viscosity of the prepreg resin increased, and the CMC became more brittle; however, it exhibited considerably higher thermal conductivity.

Carbosilane↗
DOE OSTI2025

Further Development of the Tamped Richtmyer-Meshkov Instability Method and Application to Molybdenum Dynamic Strength Calibration and Tabulation

The high pressure and high strain rate dynamic strength of Mo is experimentally and computationally investigated in the 3–20 GPa stress, 50–600 C temperature, and 10 5 –10 6 /s strain rate regimes using a modified tamped Richtmyer-Meshkov instability (RMI) method. Modifications to the tamped RMI method include a method to determine loading states during strain, a new strength calibration function based on interface shape, and a robust uncertainty quantification method. These modifications improve fidelity of the tamped RMI method, allowing evaluation of the compensating effects of pressure hardening, strain rate hardening, strain hardening, and thermal softening. The new calibration function based on interface shape is not limited to sinusoidal corrugations and could be applied to additional interface shapes. Plate impact experiments are performed at Argonne National Laboratory’s Advanced Photon Source’s Dynamic Compression Sector operated by Washington State University (DCS), driving a planar shock front through a corrugated Mo-D 2 O or Mo-C 8 F 18 interface, forcing the corrugation to significantly deform. The extent of interfacial deformation, RMI growth, is experimentally observed using X-ray phase contrast imaging at the DCS. RMI jet lengths and jet shapes are extracted from the experimental radiographs, then used to calibrate numerical simulations performed with the Sandia National Laboratories (SNL) hydrocode CTH. Mo yield strength, Y, as a function of shock pressure, P, strain rate, $\dot{\varepsilon }$, accumulated strain, ϵ, relative volumetric compression, RD , and temperature, T , is determined for each impact experiment and presented. The calibrated Mo yield strength values range 1.2–1.8 GPa, with strength generally decreasing as the impact stress increases. This trend is likely caused by thermal softening or strain localization. The tabular yield strength versus loading condition data presented in this paper can be used to fit complex strength models.

Voorhees, T. J. [Sandia National Lab. (SNL-CA), Li↗
DOE OSTI2024

A prediction model of failure threshold for shear deformation in a Zr-based bulk metallic glass

The failure of bulk metallic glasses (BMGs) during plastic deformation at room temperature is abrupt and instantaneous, while the analysis of precursor information based on avalanche events helps predict catastrophic failure. An acoustic emission (AE) signal can provide accurate precursor information for material failure, due to its sensitive and high fast calculation ability. In the current study, AE monitoring tests are carried out during uniaxial compression tests of BMGs at different strain rates. The AE experimental failure threshold, E max , is proposed on the basis of AE cumulative energy, which reflects the intensity of damage evolution at different loading conditions. Compared with the critical shear band velocity (CSBV) associated with stick-slip dynamics of serrated flow, E max is a more sensitive failure parameter since it is connected with the local microscopic changes that occur during the material response process. Here, the E max is obtained prior to reaching the CSBV since the calculation of these two avalanches analysis focuses on the different stages of shear band growth. In particular, AE events are related to the “dry” friction process in the first stage, however, the CSBV is responsible for the “viscous” glide in the second stage. Therefore, Emax is not affected by the complex interactions between the shear bands during the stick-slip process. The maximum avalanche of serrated flow, S max , is proposed as the experimental failure threshold, which depends on the applied strain rate as S max ~ $\dot{ε}$ –λ . According to the relationship of E max and S max , the theoretical failure threshold, E max , follows a criterion E max = 2545$\dot{ε}$ –λ - 4468, where λ is equivalent to 0.15 for this work. Finally, combining the different calculations and AE measurements, this model gives new insights to predict the deformation failure behavior of Zr-based BMGs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Review of the Research and Development of Brayton Cycle Technology in Nuclear Power Applications with a Focus on Compressor Technology

This study reviews the integration of Brayton Cycle (BC) systems in nuclear power generation, emphasizing their potential to enhance thermal efficiency and operational flexibility over traditional Rankine Cycle (RC) systems. Key working fluids, such as helium (He), supercritical carbon dioxide (sCO 2 ), nitrogen (N 2 ), and air, are evaluated for their performance, efficiency, and compatibility with nuclear systems. He is recognized for its high thermal conductivity and inertness at elevated temperatures, while sCO 2 demonstrates advantages in compactness and efficiency in midrange temperatures. This article also highlights the importance of compressor designs in optimizing BC performance and reviews, available compressor technologies. Axial and centrifugal compressor designs enable efficient gas compression while managing the thermal and mechanical stresses associated with high-pressure operations in nuclear systems. Combined with variable geometry components and advanced materials, these technologies address the challenges posed by varying load conditions. Despite the promising features of BC systems, several challenges persist, including high leakage rates and material degradation under extreme conditions, which necessitate robust sealing technologies and thorough testing. The insights gained from operational experiences at facilities, such as the Oberhausen II plant and the High-Temperature He Test Facility (HHV), underscore the complexities involved in designing high-temperature gas turbines for nuclear applications. This review concludes that as the nuclear industry evolves, BC systems hold significant promise for contributing to a sustainable energy future, particularly in the context of small modular reactors (SMRs) and microreactors. Further exploration of combined cycle configurations that combine BCs with RCs may enhance overall efficiency and flexibility in power generation.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Computational simulation of asymmetric phase transformation in cracked Li₇La₃Zr₂O₁₂: Variant selection and chemo-mechanical implications

Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dynamic crushing of metal lattice metamaterials: Shock mode diagrams and transition to topology-independent compaction regime

Additively manufactured lattice metamaterials offer design versatility in strength and energy absorption and provide an additional degree of freedom through the selection of the lattice topology. Under quasistatic loading, the unit cell structure can strongly affect the stiffness, yield, and post-yield behavior, but whether and to what degree the effect of lattice topology persists into dynamic loading scenarios, up to the compaction shock regime, has not been established. LLNL’ s ALE3D hydrocode was used to perform a computational investigation of dynamic loading in multiple lattice types, including the gyroid, octet, Schwarz D, and rhombic dodecahedron, under impact velocities from 0.25 to 2.25 km/s. Shock Hugoniots for each lattice topology are generated and compared, suggesting that above a critical velocity, distinctions between architectures may not persevere and compacted lattices behave similarly. Here, to investigate the transition between topology-dependent quasistatic compression and the topology-independent regime above the critical velocity, a one-dimensional elastic-linear hardening plasticity-densified solid (E-LHP-DS) shock model for lattice materials was developed that relies upon confined compression to link the quasistatic and shock mechanics. Unlike similar works, the model does not assume rigid behavior prior to yield or locking behavior at densification, allowing a richer exploration of lattice mechanics. With only six parameters, the analytical model simultaneously fit quasistatic confined compression simulations for relative densities 0.1 $≤ \bar{ρ} ≤$ 0.9 and predicted dynamic compaction behavior to traverse several distinct shock modes, each defined by a critical impact speed (equivalently, critical stresses). Comparing the numerical results to the one-dimensional E-LHP-DS shock model predictions suggests that the topology-independence under strong shocks is linked to the onset of densification, which can be predicted based on quasistatic confined compression results.

Cellular material↗
DOE OSTI2025

Generalizing the compressible pairwise interaction extended point-particle model

Ejecta physics plays an important role in material interfaces that are impacted by a strong shock wave. When a shock impacts a rough surface of solid material and melts it, the Richtmyer–Meshkov instability grows perturbations on the surface, which can eject particles. After release, the ejecta travel through the post-shock compressible flow. To accurately simulate a large number of ejecta particles, an Euler–Lagrange approach is preferred, which requires modeling the subgrid-scale physics involved with fluid–particle interactions. We generalize the previous work from Hsiao et al. (2023) to consider systems of moving particles subject to any loading shock. The following improvements were made: (1) Particles are allowed to move relative to each other (2) Non-planar shocks are accounted for along with allowing for variable shock speeds. As a result, the generalized algorithm was tested with particle-resolved simulations for canonical test cases. The results of these tests are discussed and analyzed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A hereditary integral, transient network approach to modeling permanent set and viscoelastic response in polymers

An efficient numerical framework is presented for modeling viscoelasticity and permanent set of polymers. It is based on the hereditary integral form of transient network theory, in which polymer chains belong to distinct networks each with different natural equilibrium states. Chains continually detach from previously formed networks and reattach to new networks in a state of zero stress. The free energy of these networks is given in terms of the deformation gradient relative to the configuration at which the network was born . A decomposition of the kernel for various free energies allows for a recurrence relationship to be established, bypassing the need to integrate over all time history. The technique is established for both highly compressible and nearly incompressible materials through the use of neo-Hookean, Blatz–Ko, Yeoh, and Ogden-Hill material models. Multiple examples are presented showing the ability to handle rate-dependent response and residual strains under complex loading histories.

Finite element method↗
DOE OSTI2024

Influence of high-strain-rate compression and subsequent heat treatment on (TiNbZr) 89 (AlTa) 11 refractory high-entropy alloys: Dynamic-mechanical behavior and microstructural changes

This study explored the dynamic-mechanical behavior of a novel low-density (TiNbZr) 89 (AlTa) 11 refractory high-entropy alloy (RHEA) across strain rates ranging from 1.0×10 3 to 3.5×10 3 s –1 . A significant increase in the yield and ultimate compressive strengths with rising strain rates up to 3.0×10 3 s –1 was observed and attributed to enhanced dislocation activities and stress-induced microstructural transformations. The formation of the B2 phase and Zr 5 Al 3 precipitates was found to be crucial in bolstering the alloy strength at high strain rates. Beyond strain rates of 3.0×10 3 s –1 , a decrease in strength occurred due to thermal softening and strain localization. Microstructural analyses at 3.5×10 3 s –1 revealed grain refinement, the development of micro shear bands, and dislocation tangles, which were indicative of dynamic recrystallization. Besides, the findings also revealed that the post-dynamic compression heat treatment further enhanced the hardness and microstructural stability of the alloy. These results highlight the potential of the (TiNbZr) 89 (AlTa) 11 RHEA for applications requiring materials with high strength-to-weight ratios, particularly in dynamically loaded environments. It is expected that the results of this study will further advance our fundamental understanding of the behavior of RHEAs under extreme conditions, thereby opening new avenues for material innovation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

NLR HPC Facility Power Usage Effectiveness (PUE) Data

Timeseries of Energy Systems Integration Facility (ESIF) Data Center Power Usage Effectiveness (PUE) Data provided in Parquet and compressed CSV formats Power Metrics Timeseries Fields: ts: Timestamp cooling_kw: Cooling (kilowatts) - Captures the power used by fans and pipe trace heaters associated with outdoor cooling equipment. The dedicated tower filter pump power is also captured as cooling load. energy_reuse: Energy Reuse Effectiveness hvac_kw: Heating, ventilation, and air conditioning (kilowatts) - Captures fan walls, fan coils that support the data center electrical rooms, and the make-up air unit. it_power_kw: IT equipment (kilowatts) - Captures power used by the IT equipment on the data center floor. plug_and_light_kw: Lights and utility plugs (kilowatts) - Captures power associated with the data center and dedicated mechanical room. The crank-case heater for the emergency standby generator is also captured as light and plug load. pue: Power Usage Effectiveness pump_kw: Pumps (kilowatts) - Captures power from pumps that move water in the data center Energy Recover Water loop and the Tower Water loops, and also captures power used by the boost pumps that circulate water through the fan walls. Note: The tower filter pump runs constantly to filter water from the data center cooling tower system, so 2.67 kilowatts are attributed to this pump and that is not reflected in this data field. day: Day of month Outside Weather Station Timeseries Fields: ts: Timestamp outside_air_humidity: Outside air humidity - Relative humidity percent outside_air_temp: Outside air temperature - Degrees Fahrenheit day: Day of month More detail: High-Performance Computing Data Center Power Usage Effectiveness

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A solid-state extrusion pathway extensible to upcycling aluminum scrap and end-of-life permanent magnets into ductile Al-6061/SmCo 5 structural magnetic composites

Permanent magnets offer strong potential for multifunctional structural components, but their inherent brittleness limits direct integration. Here, SolidStir® Extrusion (SSE) is used to embed SmCo5 powder into AA6061 in the solid state, and the resulting microstructure, interfacial chemistry, mechanical properties, and magnetic response are systematically characterized. SSE fragments the initial SmCo 5 powder into predominantly fine particles dispersed within an ultrafine Al matrix, with a volume fraction of ∼ 4–6 % in this proof-of-concept study. Interfacial reactions form a continuous Al–Co reaction layer, with local atomic-resolution indexing consistent with Al 13 Co 4 , together with Sm–Si-rich reaction products, while the matrix contains Sm–Co–Si-rich precipitates and Mg 2 Si, producing features from nanometers to micrometers. Compared with an Al extrudate, this architecture increases yield strength by ∼ 42 % and tensile strength by ∼ 31 %, while preserving uniform elongation, and it sustains work hardening in compression from 150 to 400 °C. The composite retained hard-magnetic behavior, with axial magnetization of 12.7 A·m 2 kg −1 at 6 T and intrinsic coercivity of 9.6 kOe. These findings demonstrate that SSE can produce Al-based load-bearing magnetic composites while preserving useful SmCo5 functionality.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A conveyor-belt magneto-optical trap of CaF

Achieving high-density samples of laser-cooled molecules is a critical step toward advancing applications in precision measurements, ultracold chemistry and quantum science. We report the experimental realization of a high-density conveyor-belt magneto-optical trap for calcium monofluoride (CaF) molecules. The obtained highly-compressed cloud has a mean radius of 64(5) μm and a peak number density of 3.6(5) × 10 10 cm −3 , a 600-fold increase over the conventional red-detuned MOTs of CaF, and the densest molecular MOT observed to date. Subsequent loading of these molecules into an optical dipole trap yields up to 2.6 × 10 4 trapped molecules at a temperature of 14(2) μK with a peak phase-space density of ~ 2.4 × 10 −6 . This opens new possibilities for a range of applications utilizing high-density, optically trapped ultracold molecules.

Yu, Scarlett S. [Harvard Univ., Cambridge, MA (Uni↗
DOE OSTI2025

Dynamic Strength and Equation of State of Epon 828 and Diethanolamine (DEA) Polymer Epoxy Under Shock Loading

Polymers are increasingly utilized in engineering applications that can experience high loading rates, necessitating increased understanding of their response under such conditions. The tamped Richtmyer-Meshkov Instability (RMI) method was used to characterize the equation of state and dynamic strength of the polymer Epon 828 cured with Diethanolamine (DEA). Plate impact experiments that drove a uniaxial shock compression wave across a sinusoidally corrugated metal-polymer interface were performed to generate shock stresses from 4-12 GPa and strain rates of approximately 1/s in the polymer. X-ray phase contrast imaging recorded the shock motion in the polymer and subsequent interface evolution. Analysis of this data yielded the polymer equation of state and, in conjunction with numerical modeling, the dynamic strength. The equation of state was validated against one-dimensional plate impact experiments from existing literature. The dynamic strength was compared to prior data for Epon 828 and related polymers at lower strain rates and found to exhibit significant strain rate and pressure-hardening effects. The strength of the Epon 828 polymer at 10 6 1/s was found to be approximately 1.5 GPa, suggesting that it is comparable to the strength of high strength metals at these dynamic conditions.

Equation of state↗
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↗