Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Inspires”

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 883 records · Page 49

Coupled geomechanical investigation of depletion-induced fault reactivation

Fault reactivation during subsurface fluid production pose significant challenges to safe and sustainable resource extraction. Here, this study presents a three-dimensional coupled geomechanical framework to investigate the processes driving fault reactivation, capturing the interactions between reservoir dynamics and geomechanical responses. Verification against theoretical estimations based on linear poroelasticity confirms the model's capacity in representing reservoir background stress responses. However, the study reveals that relying solely on background stress states can underestimate or overestimate fault reactivation potential, emphasizing the importance of including localized stress perturbations such as differential compaction and stress redistribution. Applied to a fault (M1) inspired by the geological characteristics of the Groningen field, the model shows slip initiation at 2965 m depth with 16.0 MPa depletion, aligning with field observations where seismicity occurred at approximately 3 km depth after 15.8 MPa depletion. Parametric studies reveal: (1) inelastic reservoir compaction delays fault reactivation and mitigates fault slip by reducing stress concentration, (2) higher intermediate in-situ stress magnitudes decrease the Coulomb Failure Stress (CFS) increase rate and reduce fault slip, (3) larger fault offsets amplify shear stress near the offset zone, promoting earlier reactivation and longer rupture propagation, and (4) fault permeability significantly influences pressure diffusion, with low-permeability faults leading to sharper stress changes and earlier fault destabilization. These insights highlight the critical role of geological and mechanical parameters in fault reactivation and provide a predictive framework for mitigating induced seismicity risks.

Coupled geomechanical modeling↗

The role of biochar in the degradation of sugars during fast pyrolysis of biomass

With appropriate pretreatment, sugars can be a major product from fast pyrolysis of lignocellulosic biomass. Analytical pyrolysis of pure cellulose can produce up to 60 wt% yield of levoglucosan although yields are significantly lower in continuous pyrolysis at larger scales. Secondary reactions of vaporized levoglucosan are thought to be responsible for this loss of sugar yield, suggesting changes in the design and operation of pyrolysis reactors to minimize these reactions. Micropyrolysis experiments were performed to better understand the mechanism of sugar degradation in the presence of biochar. A 57% loss in levoglucosan yield was observed for cellulose overlain with untreated biochar powder compared to the pure cellulose control sample. The addition of biochar derived from pyrolysis of untreated corn stover to a fluidized bed pyrolyzer reduced sugar yields from cellulose from 61.3 wt% to 21.3 wt% and 41.5–11.6 wt% for conventional and autothermal operation, respectively. The significant drop in sugar yield due to biochar interaction inspired change in feeder configuration for the fluidized bed pyrolyzer to reduce vapor-char interactions. Biomass feeding was changed from in-bed to above-bed injection, which allowed significant devolatilization to occur above the layer of biochar that exists at the surface of the bed. Finally, by reducing secondary reactions, bio-oil and sugar yields increased by 7.9% and 14%, respectively, for autothermal pyrolysis.

09 BIOMASS FUELS↗

On the existence of orthorhombic martensite in a near-α titanium base alloy used for additive manufacturing

Additive manufacturing is a state-of-the-art production technology to produce tailor-made and highly complex parts. Among various other alloys, Ti base alloys are frequently used in this manufacturing technique due to their well-balanced properties and their wide range of applications. Allotropic phases and the occurrence of athermal phase transformations are the main reasons why these alloys hold a great development potential and are the basis of extensive use. High cooling rates during manufacturing lead to martensitic phases and the formation of nanometer-sized microstructures resulting in extraordinary strength. Simultaneously, such high cooling rates cause a high amount of lattice defects and the occurrence of residual stresses, which finally may result in delamination effects and cracks. Usually, a common approach to reduce residual stresses during additive manufacturing is to decrease thermal gradients by increasing the heat input or preheating the building platform. Instead of applying the typical approaches to lower thermal gradients, this study deals with the origin of the ‘softer’ orthorhombic martensite by accelerating the solidification process. The implementation of the orthorhombic phase in bulk components was inspired by a new phase transformation herein reported for the first time in the powder material, which also validates the possible occurrence of two martensitic phases in the same alloy. Various sophisticated characterization techniques like high energy and high-temperature X-ray diffraction, high-resolution transmission electron microscopy as well as atom probe tomography were applied to characterize this softer orthorhombic martensitic phase in detail aiming to highlight the opportunities accompanied by this new approach for additive manufacturing of titanium alloys.

36 MATERIALS SCIENCE↗

Magnetic and transport properties of Fe-doped Weyl semimetal Co 3 Sn 2 S 2

Ferromagnetic Weyl semimetal Co 3 Sn 2 S 2 has recently attracted a lot of interest on account of its large non-saturating magnetoresistance and the high carrier mobility common to topological semimetals. Moreover, the material also provides a platform for studying the interplay between the intrinsic magnetism and topological band structures, which generates a significant anomalous Hall effect. Chemical doping is the usual approach to adjust materials’ properties, and it can lead to fascinating results. Here, we present the influence of Fe-doping on the magnetic and transport properties of polycrystalline Co 3 Sn 2 S 2 . The Fe-doping suppresses the sample’s magnetic order, increases the electrical resistivity, results in a skew-scattering contribution to the anomalous Hall effect and Fe impurities cause the Kondo effect. In addition, we report on the impact of Fe-doping on the magnetoresistance of Co 3 Sn 2 S 2 , where Fe-doping leads to a dramatic decline in magnetoresistance. Furthermore, the hysteresis behavior in the low magnetic field region, detected previously in nanoflake samples only, is observed in both pure and Fe-doped Co 3 Sn 2 S 2 bulk samples. Finally, our findings serve as a reference and source of inspiration for future chemical doping research on topological semimetals.

36 MATERIALS SCIENCE↗

Metal hydrides: a historical perspective

Metal hydrides are known for their outstanding performance as materials for hydrogen storage and processing. These materials find applications for short- and long-term energy storage, compression and supply of hydrogen gas, thermal energy storage, as electrodes and electrolytes in rechargeable batteries, for the microstructural optimisation of functional materials, in thin film technologies, as catalysts, getters and in many other uses. After the discovery of the first binary metal hydrides back in the 19th century, their studies covered all possible binary M-H systems and expanded rapidly into the field of ternary hydrides following the recognition of the excellent hydrogen storage performance of LaNi 5 - and TiFe-based materials, which operate efficiently at room temperature and at near-ambient H 2 pressures. This review aims to provide an overview of the early works, as well as selected recent results on various classes of metal hydrides. It also covers the recent activities from the major contributing countries and continents, including USA, Europe, Japan, China and Australia. These studies relate to achieving the hydrogen storage systems goals set by the Department of Energy in the United States which inspired the research activities at the national and international level, through execution of the tasks on hydrogen-based energy storage managed by the International Energy Agency. The review is prepared by international experts in the field and covers the most important past developments and also presents the recent achievements in the field.

08 HYDROGEN↗

KEBLM: Knowledge-Enhanced Biomedical Language Models

Pretrained language models (PLMs) have demonstrated strong performance on many natural language processing (NLP) tasks. Despite their great success, these PLMs are typically pretrained only on unstructured free texts without leveraging existing structured knowledge bases that are readily available for many domains, especially scientific domains. As a result, these PLMs may not achieve satisfactory performance on knowledge-intensive tasks such as biomedical NLP. Comprehending a complex biomedical document without domain-specific knowledge is challenging, even for humans. Inspired by this observation, we propose a general framework for incorporating various types of domain knowledge from multiple sources into biomedical PLMs. We encode domain knowledge using lightweight adapter modules, bottleneck feed-forward networks that are inserted into different locations of a backbone PLM. For each knowledge source of interest, we pretrain an adapter module to capture the knowledge in a self-supervised way. We design a wide range of self-supervised objectives to accommodate diverse types of knowledge, ranging from entity relations to description sentences. Once a set of pretrained adapters is available, we employ fusion layers to combine the knowledge encoded within these adapters for downstream tasks. Each fusion layer is a parameterized mixer of the available trained adapters that can identify and activate the most useful adapters for a given input. Our method diverges from prior work by including a knowledge consolidation phase, during which we teach the fusion layers to effectively combine knowledge from both the original PLM and newly-acquired external knowledge using a large collection of unannotated texts. After the consolidation phase, the complete knowledge-enhanced model can be fine-tuned for any downstream task of interest to achieve optimal performance. Extensive experiments on many biomedical NLP datasets show that our proposed framework consistently improves the performance of the underlying PLMs on various downstream tasks such as natural language inference, question answering, and entity linking. These results demonstrate the benefits of using multiple sources of external knowledge to enhance PLMs and the effectiveness of the framework for incorporating knowledge into PLMs. Finally, while primarily focused on the biomedical domain in this work, our framework is highly adaptable and can be easily applied to other domains, such as the bioenergy sector.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

H 2 O-assisted O 2 reduction by H 2 on Pt and PtAu bimetallic nanoparticles: Influences of composition and reactant coverages on kinetic regimes, rates, and selectivities

Hydrogen peroxide (H 2 O 2 ) can replace hazardous oxidants in industrial processes but is currently too expensive for many such applications. While direct synthesis of H 2 O 2 (H 2 + O 2 → H 2 O 2 ) may reduce costs in comparison to incumbent technology, current catalysts lack the requisite stability and selectivity. Here, we examine the direct synthesis of H 2 O 2 on bimetallic Pt 1 Au x (0 ≤ x ≤ 230) and Pt catalysts at steady-state in pure water and relate kinetic parameters for H 2 O 2 and H 2 O formation to possible active site structures informed by complementary characterization methods. X-ray photoelectron spectra show significant Pt surface enrichment compared to the bulk composition. Analysis of infrared spectra of mixed monolayers of 12 CO* and 13 CO* indicate that Pt and Au form substitutional surface alloys. The Pt 1 Au x nanoparticles with the greatest mole fractions of Au predominantly expose Pt monomers (i.e., isolated Pt atoms), yet Pt atoms exposed upon all these nanoparticles possess electronic structures distinct from bulk Pt. Despite these differences, rate measurements are consistent with product formation through proton-electron transfer pathways for all Pt 1 Au x catalysts. In situ XAS indicate that Pt remains metallic during H 2 O 2 synthesis. Under the most oxidizing conditions, selectivities toward H 2 O 2 increase strongly with the Au to Pt ratio from 2% for monometallic Pt to 85% for Pt 1 Au 170 . However, selectivities are similar among all catalysts within reducing conditions. Comparisons of apparent activation enthalpies for the formation of H 2 O 2 and H 2 O across these catalysts and the range of conditions suggest that Pt monomers within Au provide the greatest selectivities for H 2 O 2 formation, because these active sites present high barriers for O-O bond rupture. Further, selectivities decrease with increasing ratios of H 2 to O 2 pressures, because Pt atoms aggregate and form oligomers that readily dissociate dioxygen intermediates. The combined use of spectroscopy, kinetics, and concepts employed in reaching these conclusions take inspiration from the legacy of Prof. Michel Boudart, and specifically his elegant methods for interrogating bimetallic catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

GenMod: A generative modeling approach for spectral representation of PDEs with random inputs

Here, we propose a method for quantifying uncertainty in high-dimensional PDE systems with random parameters, where the number of solution evaluations is small. Parametric PDE solutions are often approximated using a spectral decomposition based on polynomial chaos expansions. For the class of systems we consider (i.e., high dimensional with limited solution evaluations) the coefficients are given by an underdetermined linear system in a regression formulation. This implies additional assumptions, such as sparsity of the coefficient vector, are needed to approximate the solution. Here, we present an approach where we assume the coefficients are close to the range of a generative model that maps from a low to a high dimensional space of coefficients. Our approach is inspired be recent work examining how generative models can be used for compressed sensing in systems with random Gaussian measurement matrices. Using results from PDE theory on coefficient decay rates, we construct an explicit generative model that predicts the polynomial chaos coefficient magnitudes. The algorithm we developed to find the coefficients, which we call GenMod, is composed of two main steps. First, we predict the coefficient signs using Orthogonal Matching Pursuit. Then, we assume the coefficients are within a sparse deviation from the range of a sign-adjusted generative model. This allows us to find the coefficients by solving a nonconvex optimization problem, over the input space of the generative model and the space of sparse vectors. We obtain theoretical recovery results for a Lipschitz continuous generative model and for a more specific generative model, based on coefficient decay rate bounds. We examine three high-dimensional problems and show that, for all three examples, the generative model approach outperforms sparsity promoting methods at small sample sizes.

97 MATHEMATICS AND COMPUTING↗

Learning stochastic dynamics with statistics-informed neural network

We introduce a machine-learning framework named statistics-informed neural network (SINN) for learning stochastic dynamics from data. This new architecture was theoretically inspired by a universal approximation theorem for stochastic systems, which we introduce in this paper, and the projection-operator formalism for stochastic modeling. Here, we devise mechanisms for training the neural network model to reproduce the correct statistical behavior of a target stochastic process. Numerical simulation results demonstrate that a well-trained SINN can reliably approximate both Markovian and non-Markovian stochastic dynamics. We demonstrate the applicability of SINN to coarse-graining problems and the modeling of transition dynamics. Furthermore, we show that the obtained reduced-order model can be trained on temporally coarse-grained data and hence is well suited for rare-event simulations.

97 MATHEMATICS AND COMPUTING↗

Multi-frequency progressive refinement for learned inverse scattering

Interpreting scattered acoustic and electromagnetic wave patterns is a computational task that enables remote imaging in a number of important applications, including medical imaging, geophysical exploration, sonar and radar detection, and nondestructive testing of materials. However, accurately and stably recovering an inhomogeneous medium from far-field scattered wave measurements is a computationally difficult problem, due to the nonlinear and non-local nature of the forward scattering process. We design a neural network, called Multi-Frequency Inverse Scattering Network (MFISNet), and a training method to approximate the inverse map from far-field scattered wave measurements at multiple frequencies. We consider three variants of MFISNet, with the strongest performing variant inspired by the recursive linearization method — a commonly used technique for stably inverting scattered wavefield data — that progressively refines the estimate with higher frequency content. MFISNet outperforms past methods in regimes with high-contrast, heterogeneous large objects, and inhomogeneous unknown backgrounds.

97 MATHEMATICS AND COMPUTING↗

An explicit, energy-conserving particle-in-cell scheme

We present an explicit temporal discretization of particle-in-cell schemes for the non-relativistic Vlasov equation that results in exact energy conservation when combined with an appropriate spatial discretization. The scheme is inspired by a simple, second-order explicit scheme that conserves energy exactly in the Eulerian context. We show that direct translation to particle-in-cell does not result in strict conservation, but derive a simple correction based on an analytically solvable optimization problem that recovers conservation. While this optimization problem is not guaranteed to have a real solution for every particle, we provide a correction that makes imaginary values extremely rare and still admits $\mathcal{O}$(10 –12 ) fractional errors in energy for practical simulation parameters. We present the scheme in both electrostatic – where we use the Ampère formulation – and electromagnetic contexts. With an electromagnetic field solve, the field update is most naturally linearly implicit, but the more computationally intensive particle update remains fully explicit. Here, we also show how the scheme can be extended to use the fully explicit leapfrog and pseudospectral analytic time-domain (PSATD) field solvers. The scheme is tested on standard kinetic plasma problems, confirming its conservation properties.

Energy conservation↗

Bimetallic NiCo boride nanoparticles confined in a MXene network enable efficient ambient ammonia electrosynthesis

Ambient electrocatalytic nitrogen fixation is an emerging technology for green ammonia synthesis, but the absence of optimized, stable and performant catalysts can render its practical application challenging. Herein, bimetallic NiCo boride nanoparticles confined in MXene are shown to accomplish high-performance nitrogen reduction electrolysis. Taking advantage of the synergistic effect in specific compositions with unique electronic d and p orbits and typical architecture of rich nanosized particles embedded in the interconnected conductive network, the synthesized MXene@NiCoB composite demonstrates extensive improvements in nitrogen molecule chemisorption, active area exposure and charge transport. As a result, optimal NH 3 yield rate of 38.7 μg h -1 mgcat. -1 and Faradaic efficiency of 6.92% are acquired in 0.1 M Na 2 SO 4 electrolyte. Moreover, the great catalytic performance can be almost entirely maintained in the cases of repeatedly-cycled and long-term electrolysis. Theoretical investigations reveal that the nitrogen reduction reaction on MXene@NiCoB catalyst proceeds according to the distal pathway, with a distinctly-reduced energy barrier relative to the Co 2 B counterpart. In conclusion, this work may inspire a new route towards the rational catalyst design for the nitrogen reduction reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spectroscopic and computational investigations of organometallic complexation of group 12 transition metals by methanobactins from Methylocystis sp. SB2

Methanotrophic bacteria catalyze the aerobic oxidation of methane to methanol using Cu-containing enzymes, thereby exerting a modulating influence on the global methane cycle. To facilitate the acquisition of Cu ions, some methanotrophic bacteria secrete small modified peptides known as “methanobactins,” which strongly bind Cu and function as an extracellular Cu recruitment relay, analogous to siderophores and Fe. In addition to Cu, methanobactins form complexes with other late transition metals, including the Group 12 transition metals Zn, Cd, and Hg, although the interplay among solution-phase configurations, metal interactions, and the spectroscopic signatures of methanobactin-metal complexes remains ambiguous. In this study, the complexation of Zn, Cd, and Hg by methanobactin from Methylocystis sp. strain SB2 was studied using a combination of absorbance, fluorescence, extended x-ray absorption fine structure (EXAFS) spectroscopy, and time-dependent density functional theory (TD-DFT) calculations. We report changes in sample absorbance and fluorescence spectral dynamics, which occur on a wide range of experimental timescales and characterize a clear stoichiometric complexation dependence. Mercury L3-edge EXAFS and TD-DFT calculations suggest a linear model for Hg--S coordination, and TD-DFT suggests a tetrahedral model for Zn 2+ and Cd 2+ . We observed an enhancement in the fluorescence of methanobactin upon interaction with transition metals and propose a mechanism of complexation-hindered isomerization drawing inspiration from the wild-type Green Fluorescent Protein active site. Collectively, our results represent the first combined computational and experimental spectroscopy study of methanobactins and shed new light on molecular interactions and dynamics that characterize complexes of methanobactins with Group 12 transition metals.

59 BASIC BIOLOGICAL SCIENCES↗

Visualizing the gas channel of a monofunctional carbon monoxide dehydrogenase

Carbon monoxide dehydrogenase (CODH) plays an important role in the processing of the one–carbon gases carbon monoxide and carbon dioxide. In CODH enzymes, these gases are channeled to and from the Ni-Fe-S active sites using hydrophobic cavities. In this work, we investigate these gas channels in a monofunctional CODH from Desulfovibrio vulgaris, which is unusual among CODHs for its oxygen-tolerance. By pressurizing D. vulgaris CODH protein crystals with xenon and solving the structure to 2.10 Å resolution, we identify 12 xenon sites per CODH monomer, thereby elucidating hydrophobic gas channels. We find that D. vulgaris CODH has one gas channel that has not been experimentally validated previously in a CODH, and a second channel that is shared with Moorella thermoacetica carbon monoxide dehydrogenase/acetyl-CoA synthase (CODH/ACS). This experimental visualization of D. vulgaris CODH gas channels lays groundwork for further exploration of factors contributing to oxygen-tolerance in this CODH, as well as study of channels in other CODHs. We dedicate this publication to the memory of Dick Holm, whose early studies of the Ni-Fe-S clusters of CODH inspired us all.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bioinspired oxidation of benzyl alcohol: The role of environment and nuclearity of the catalyst evaluated by multivariate analysis

Inspired by copper-containing enzymes such as galactose oxidase and catechol oxidase, in which distinct coordination environments and nuclearities lead to specific catalytic activities, we summarize here the catalytic properties of dinuclear and mononuclear copper species towards benzyl alcohol oxidation using a multivariate statistical approach. Here, the new dinuclear [Cu 2 (μ-L 1 )(μ-pz)] 2+ (1) is compared against the mononuclear [CuL 2 Cl] (2), where (L 1 ) - and (L 2 ) - are the respective deprotonated forms of 2,6-bis((bis(pyridin-2-ylmethyl)amino)methyl)-4-methylphenol, and 3-((bis(pyridin-2-ylmethyl)amino)methyl)-2-hydroxy-5-methylbenzaldehyde and (pz) - is a pyrazolato bridge. Copper(II) perchlorate (CP) is used as control. The catalytic oxidation of benzyl alcohol is pursued, aiming to assess the role of the ligand environment and nuclearity. The multivariate statistical approach allows for the search of optimal catalytic conditions, considering variables such as catalyst load, hydrogen peroxide load, and time. Species 1, 2 and CP promoted selective production of benzaldehyde at different yields, with only negligible amounts of benzoic acid. Under normalized conditions, 2 showed superior catalytic activity. This species is 3.5-fold more active than the monometallic control CP, and points out to the need for an efficient ligand framework. Species 2 is 6-fold more active than the dinuclear 1, and indicates the favored nuclearity for the conversion of alcohols into aldehydes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A constitutive framework for rocks undergoing solid dissolution

Here, we formulate a time-dependent damage theory for rocks subjected to mechanical deformation and solid dissolution. The constitutive description is inspired by the transition state theory, which states that the rate of dissolution is a function of the reactive surface area measured through the crack density in the volume. We use a gradient-enhanced damage framework in which damage depends on the deformation of the material as well as on the amount of solid mass dissolved over time. The gradient-enhanced formulation is characterized by a three-field variational formulation with the solid displacement, nonlocal equivalent strain, and nonlocal rate of solid dissolution as the basic state variables. Traditionally, time-independent damage theories have only allowed damage to increase with increasing external load. In the proposed framework, the degree of damage may increase due to solid dissolution even when the external load is held fixed. In this way, solid dissolution is viewed as a process that is responsible for bringing about rate-dependent effects such as creep and stress-relaxation, which are two common features of geomaterial behavior.

42 ENGINEERING↗

A statistical mechanics framework for polymer chain scission, based on the concepts of distorted bond potential and asymptotic matching

To design increasingly tough, resilient, and fatigue-resistant elastomers and hydrogels, the relationship between controllable network parameters at the molecular level (bond type, non-uniform chain length, entanglement density, etc.) to macroscopic quantities that govern damage and failure must be established. Many of the most successful constitutive models for elastomers have been rooted in statistical mechanical treatments of polymer chains. Typically, such constitutive models have used variants of the freely jointed chain model with rigid links. However, since the free energy state of a polymer chain is dominated by enthalpic bond distortion effects as the chain approaches its rupture point, bond extensibility ought to be accounted for if the model is intended to capture chain rupture. To that end, a new bond potential is supplemented to the freely jointed chain model (as derived in the u FJC framework of Buche and Silberstein (2021) and Buche et al. (2022)), which we have extended to yield a tractable, closed-form model of single chain behavior that should be amenable to continuum-level constitutive model development. Inspired by the asymptotically matched u FJC model response in both the low/intermediate chain force and high chain force regimes, a simple, quasi-polynomial bond potential energy function is derived. This bond potential exhibits harmonic behavior near the equilibrium state and anharmonic behavior for large bond stretches tending to a characteristic energy plateau (akin to the Lennard-Jones and Morse bond potentials). Using this bond potential, approximate yet highly-accurate analytical functions for bond stretch and chain force dependent upon chain stretch are established. Then, using this polymer chain model, a stochastic thermal fluctuation-driven chain rupture framework is developed. This framework is based upon a force-modified tilted bond potential that accounts for distortional bond potential energy, allowing for the derivation and subsequent calculation of the dissipated chain scission energy. Here, the cases of rate-dependent and rate-independent scission are accounted for throughout the rupture framework. The impact of Kuhn segment number on chain rupture behavior is also investigated. The model is fit to single chain mechanical response data collected from atomic force microscopy tensile tests for validation and to glean deeper insight into the molecular physics taking place. Due to their analytical nature, this polymer chain model and the associated rupture framework can in the future be implemented in finite element models accounting for fracture and fatigue in polydisperse elastomer networks.

36 MATERIALS SCIENCE↗

Strongly nonlinear wave propagation in elasto-plastic metamaterials: Low-order dynamic modeling

Nonlinear elastic metamaterials are known to support a variety of dynamic phenomena that enhance our capacity to manipulate elastic waves. Since these properties stem from complex, subwavelength geometry, full-scale dynamic simulations are often prohibitively expensive at scales of interest. Prior studies have therefore utilized low-order effective medium models, such as discrete mass-spring lattices, to capture essential properties in the long-wavelength limit. While models of this type have been successfully implemented for a wide variety of nonlinear elastic systems, they have predominantly considered dynamics depending only on the instantaneous kinematics of the lattice, neglecting history-dependent effects, such as wear and plasticity. Here, to address this limitation, the present study develops a lattice-based modeling framework for nonlinear elastic metamaterials undergoing plastic deformation. Due to the history- and rate-dependent nature of plasticity, the framework generally yields a system of differential-algebraic equations whose computational cost is significantly greater than an elastic system of comparable size. We demonstrate the method using several models inspired by classical lattice dynamics and continuum plasticity theory and explore means to obtain empirical plasticity models for general geometries, thereby gaining insight into the influence of microstructural plasticity on effective material performance, which can be used to improve the design of nonlinear mechanical metamaterials.

Dynamic simulation↗