SEARCH · Engineering Papers

Results for “inelastic deformations”

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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Materials with periodic internal structure: Computation based on homogenization and comparison with experiment

The combination of thermal and mechanical loading expected in practice means that constitutive equations of metal matrix composites must be developed which deal with time-independent and time-dependent irreversible deformation. Also, the internal state of composites is extremely complicated which underlines the need to formulate macroscopic constitutive equations with a limited number of state variables which represent the internal state at the micro level. One available method for calculating the macro properties of composites in terms of the distribution and properties of the constituent materials is the method of homogenization whose formulation is based on the periodicity of the substructure of the composite. A homogenization procedure was developed which lends itself to the use of the finite element procedure. The efficiency of these procedures, to determine the macroscopic properties of a composite system from its constituent properties, was demonstrated utilizing an aluminum plate perforated by directionally oriented slits. The selection of this problem is based on the fact that, extensive experimental results exist, the macroscopic response is highly anisotropic, and that the slits provide very high stress gradients which severely test the effectiveness of the computational procedures. Furthermore, both elastic and plastic properties were investigated so that the application to practical systems with inelastic deformation should be able to proceed without difficulty. The effectiveness of the procedures was rigorously checked against experimental results and with the predictions of approximate calculations. Using the computational results it is illustrated how macroscopic constitutive equations can be expressed in forms of the elastic and limit load behavior.

Jansson, S.↗

NASA NTRS2005

Surface Collisions Involving Particles and Moisture (SCIP'M)

Experiments were performed on the collision of a solid sphere with a nearly horizontal flat surface covered with a thin layer of viscous liquid. High-speed collisions were obtained by dropping the ball onto the surface from various heights, using gravitational acceleration. Low-speed collisions were obtained using pendulums with long strings or by launching the balls at low velocities in the reduced-gravity environment of parabolic flight. The sphere bounces only when the impact velocity exceeds a critical value. The coefficient of restitution (ratio of rebound velocity to impact velocity) increases with increasing impact velocity above the critical value, indicating the increasing relative importance of elastic deformation to viscous dissipation. The critical impact velocity increases, and the coefficient of restitution decreases, with increasing viscosity or thickness of the liquid layer and with decreasing density or size of the sphere. The ratio of the wet and dry coefficients is expressed as a function of the Stokes number (ratio of particle inertia and viscous forces), showing good agreement between theory and experiment. Similar experiments were performed with the flat surface inclined at various angles to the approaching sphere. A modified Stokes number, which is a measure of the ratio of inertia of the sphere in the normal direction to the viscous forces exerted by the fluid layer, was used for the analysis of oblique collisions. Even for these oblique collisions, it was found that no rebound of the ball was observed below a certain critical Stokes number. The coefficient of normal restitution, defined as a ratio of normal rebound velocity to normal approach velocity, was found to increase beyond the critical Stokes number and even out as it approaches the value for dry restitution at high Stokes numbers. It was also found that, for smooth spheres like steel, the normal restitution at the same modified Stokes number is independent of the angle of impact. The tangential coefficient of restitution, defined as the ratio of tangential rebound velocity to tangential approach velocity, is found to be nearly unity, except for very low approach velocities. Thus, as a first approximation, the theories that predict the coefficient of restitution for head-on wet collisions can be extended to predict the coefficient of normal restitution for oblique wet collisions. Additional experiments were performed with soft surfaces in which a porous cloth or sponge layer was placed over the hard, flat surface. In these experiments, the coefficient of restitution was found to decrease with increasing impact velocity, due to inelastic losses in the soft material. A model combining inelastic deformation and flow through porous media was developed to describe these findings.

Davis, Robert H.↗

DOE OSTI2020

A multiscale cohesive law for carbon fiber networks

Better predictive models of mechanical failure in low-weight heat shield composites would aid material certification for missions with aggressive atmospheric entry conditions. In this study, we develop such a model for the rapid engineering analysis of the failure limits of phenolic impregnated carbon ablator (PICA) - a leading heat shield material whose structural component is a carbon fiber network. We hypothesize inelastic deformation failure mechanisms and model their behavior using molecular dynamics simulations to calculate the binding energy. We then upscale this binding energy to the macroscale using a renormalization argument. The approach delivers insightful and reasonably accurate macroscale predictions that compare favorably to experiments. In application, the model is validated for a particular variety of PICA by comparison to experiment and would then be used to study design scenarios in different entry conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Assessment of thermal and radiation induced creep in the dual cooled lead lithium blanket

The creep in recently designed blanket for the fusion nuclear science facility is described by decomposing it into thermal, irradiation and cavity swelling creep. Initial estimates of thermal creep using elastic high temperature rules showed that the recent heat transfer improvements at the first wall (FW) is feasible, and the FW can withstand heat flux of 0.35 MW/m 2 at 3000 hrs without rupture. A viscoplastic model that is based on Norton's power law and relies on Multiphysics coupling of solid mechanics and heat transfer modules is used to capture the inelastic deformation in the blanket. Results showed that the design peak heat flux of 0.25 MW/m 2 produced maximum thermal creep strain of 0.45 % and the relaxation of the thermal stress at the FW. Irradiation creep is prescribed to be proportional to the displacement damage dose and applied stress. The displacements from irradiation creep radially decrease from the FW to the back wall, but the maximum deformation is found at the back wall that is connected to the Helium manifold due to the high stresses at the region. In conclusion, cavity swelling creep is dominant at steady state and combined with radiation creep to produce displacement of about 8 mm at the FW.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Training material models using gradient descent algorithms

High temperature design requires accurate constitutive models to describe material inelastic deformation and failure behavior. Oftentimes, calibrating accurate models devolves into the problem of fitting the model parameters against experimental test data. Here, we present the pyopmat package, an open source framework for calibrating constitutive models against experiment data subjected to various loading conditions using machine learning techniques. The package calculates the exact gradient of the model response with respect to the parameters using a combination of automatic differentiation and the adjoint method. Given this exact gradient, we compare the performance of several gradient-based optimization techniques in fitting realistic constitutive models against data. Here, we demonstrate the efficiency and accuracy of our package through example problems using both synthetic data, generated using known parameter sets, under monotonic and cyclic loading conditions and also with an example applying the techniques developed here to actual high temperature creep-fatigue test data.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Coupled chemo-mechanical modeling of point-defect diffusion in a crystal plasticity fast Fourier transform framework

Below the yield strength and at moderate-to-high homologous temperatures, the inelastic deformation of metals is mostly governed/rate-controlled by vacancy diffusion-mediated processes. As a function of grain size, stress, temperature and dislocation content, vacancies (or atoms) can adopt preferential diffusion paths across grain interiors, along grain boundaries, or towards and along dislocations, resulting in climb and self-climb. In the steady state and under constant load, grain boundary and grain bulk vacancy diffusion-mediated plasticity have been described in seminal works by Coble and by Nabarro and Herring, respectively. Yet, the interplay between all aforementioned potential diffusion pathways has not been comprehensively mapped. This work presents a thermodynamically-consistent full-field model integrated within a voxel-based elasto-viscoplastic fast Fourier transform framework, which considers the coupling between the diffusion-mediated plasticity mechanisms. In the proposed approach, the kinetics and kinematics of plastic deformation due to vacancy diffusion along grain boundaries and grain bulk, as well as the exchange between grain boundaries and bulk are described explicitly. A homogenization approach at the voxel level is further introduced to simultaneously consider bulk and grain boundary diffusion in a numerically efficient fashion. The new formulation predicts the expected strain rate dependencies and the scaling of the steady-state creep rate with respect to grain size, temperature, and stress. Finally, the model predicts the transition from grain bulk to grain boundary-dominated diffusion with reduction in grain size, a significant step towards capturing transitions in deformation behavior without any phenomenological or ad-hoc adjustments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Spherical Indentation and Implementation of S 3 /P for yield stress determination of brittle materials

Here, a mathematically transparent and robust experimental method has been developed to estimate the yield stress of brittle materials through the analysis of depth-sensing spherical indentation. Employing Hertzian contact mechanics, an elastically invariant ratio based on the simple equation S 3 / P = 6RE r 2 , (where S and P are contact stiffness and indentation load, respectively) has been derived that enables more accurate and confident determination of the transition from elastic to inelastic deformation; a transition that the yield stress dictates and represents. Using two diamond spheres with radii of 3.2 and 8.6 μm, the indentation test method and analyses are applied to two vitreous silicates: Corning's HPFS 7980® fused silica and Vitro's Starphire® soda lime silicate. The estimated yield strengths are 8.15 GPa ± 2.5% for the fused silica and 6.1 GPa ± 3.3% for the soda-lime silicate, and both were independent of indenter radius. Verification of this new experimental method is demonstrated with an as-drawn titanium by showing equivalence of measured yield stress by its spherical indentation and that from uniaxial compression testing. This method will enable easier and more confident estimation of yield stress in brittle materials - a property that historically has been elusive to measure for these materials using common laboratory mechanical test methods.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Achieving 5.9% elastic strain in kilograms of metallic glasses: Nanoscopic strain engineering goes macro

The ideal elastic limit is the upper bound of the achievable strength and elastic strain of solids. However, the elastic strains that bulk materials can sustain are usually below 2%, due to the localization of inelastic deformations at the lattice scale. In this study, we achieved >5% elastic strain in bulk quantity of metallic glass, by exploiting the more uniform and smaller-magnitude atomic-scale lattice strains of martensitic transformation as a loading medium in a bulk metallic nanocomposite. The self-limiting nature of martensitic transformation helps to prevent lattice strain transfer that leads to the localization of deformation and damage. This lattice strain egalitarian strategy enables bulk metallic materials in kilogram-quantity to achieve near-ideal elastic limit. This concept is verified in a model in situ bulk amorphous (TiNiFe)-nanocrystalline (TiNi(Fe)) composite, in which the TiNiFe amorphous matrix exhibits a maximum tensile elastic strain of similar to 5.9%, which approaches its theoretical elastic limit. As a result, the model bulk composite possesses a large recoverable strain of similar to 7%, a maximum tensile strength of above 2 GPa, and a large elastic resilience of similar to 79.4 MJ/m 3 . The recoverable strain and elastic resilience are unmatched by known high strength bulk metallic materials. This design concept opens new opportunities for the development of high-performance bulk materials and elastic strain engineering of the physiochemical properties of glasses.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Evaluation of external tubular configurations for a high-temperature chloride molten salt solar receiver operating above 700°C

Next-generation concentrating solar power (CSP) tower technologies target operating temperatures exceeding 700°C to increase the thermal-to-electric conversion efficiency. Molten chloride salts are one possible alternative to current commercial molten nitrate salts to enable the higher operating temperature. This paper analyzes the predicted optical, thermal-fluids, and structural performance of traditional external tubular solar receiver configurations applied with a chloride salt heat transfer fluid (HTF) and inlet/outlet temperatures of 500°C/735°C, and considers sensitivity analysis and optimization relative to receiver sizing, tube sizing, number of panels, flow circuit configurations, and solar flux concentration under constraints on internal velocity, pressure drop, wall thickness, and required creep-fatigue lifetime. The high temperature conditions increase the significance of inelastic deformation mechanisms such as creep relative to that expected in commercial 565°C nitrate salt designs. High-temperature creep and creep-fatigue damage in the metal alloy tubes are the key factors that limit allowable solar flux concentration and achievable receiver thermal efficiency at the near-800°C wall temperature conditions. For a traditional external cylindrical receiver configuration, the design parameters and conditions capable of satisfying all constraints produced, at best, a design point receiver efficiency of 78.2%, or 80.5% when excluding receiver intercept efficiency. Variation in the optimal receiver performance relative to uncertainty in the binding maximum velocity, minimum wall thickness, and minimum lifetime constraints is presented.

14 SOLAR ENERGY↗

DOE OSTI2024

New Nonreactive Force Field for Accurate Molecular Dynamics Simulations of TATB at Extreme Conditions

Insensitive high explosives based on TATB (1,3,5-triamino-2,4,6-trinitrobenzene) are needed in applications when safety is of paramount importance, but the basic material properties that give rise to its insensitivity are not fully understood. Molecular dynamics modeling using empirical force fields (FFs) has been the main route to characterize many complicated dynamical properties of TATB single crystal, but these FFs have not been comprehensively tested at extreme conditions typical of detonation. We collect a benchmark data set of (quasi)static TATB physical properties as determined by experiments and electronic structure calculations and apply this data set to validate four existing TATB FFs along with a new TATB FF that we develop here and denote as the CEA-LLNL-Missouri (CLM) FF. Benchmark data include vibrational spectra, the TATB crystal temperature–pressure–volume equation of state and lattice parameters, properties of TATB crystal polymorphs and transitions to the gaseous and liquid states, dimer energy landscapes, the pressure-dependent elastic tensor, and the energy landscape for inelastic deformation via sliding of TATB crystal layers. As a general assessment, we find that the two existing nonreactive FFs are more accurate in describing TATB’s physical properties compared to the two variants of the ReaxFF reactive FF considered. The new CLM FF is found to consistently yield similar or better agreement with experiments and electronic structure theory than any of the existing FF models, and it presents a distinct improvement in accurately modeling TATB elasticity and equation of state. So this work is expected to help improve the accuracy of FF-based modeling of complicated dynamic responses that ultimately govern the safety and performance characteristics of this material.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Evidence of non-isentropic release from high residual temperatures in shocked metals measured with ultrafast x-ray diffraction

Shock experiments are widely used to understand the mechanical and electronic properties of matter under extreme conditions. However, after shock loading to a Hugoniot state, a clear description of the post-shock thermal state and its impacts on materials is still lacking. We used diffraction patterns from 100-fs x-ray pulses to investigate the temperature evolution of laser-shocked Al–Zr metal film composites at time delays ranging from 5 to 75 ns driven by a 120-ps short-pulse laser. We found significant heating of both Al and Zr after shock release, which can be attributed to heat generated by inelastic deformation. A conventional hydrodynamic model that employs (i) typical descriptions of Al and Zr mechanical strength and (ii) elevated strength responses (which might be attributed to an unknown strain rate dependence) did not fully account for the measured temperature increase, which suggests that other strength-related mechanisms (such as fine-scale void growth) could play an important role in thermal responses under shock wave loading/unloading cycles. Our results suggest that a significant portion of the total shock energy delivered by lasers becomes heat due to defect-facilitated plastic work, leaving less converted to kinetic energy. This heating effect may be common in laser-shocked experiments but has not been well acknowledged. High post-shock temperatures may induce phase transformation of materials during shock release. Another implication for the study is the preservability of magnetic records from planetary surfaces that have a shock history from frequent impact events.

42 ENGINEERING↗

DOE OSTI2021

The mechanics of brittle granular materials with coevolving grain size and shape

The influence of particle shape on the mechanics of sand is widely recognized, especially in mineral processing and geomechanics. However, most existing continuum theories for engineering applications do not encompass the morphology of the grains and its evolution during comminution. Similarly, the relatively few engineering models accounting for grain-scale processes tend to idealize particles as spheres, with their diameters considered as the primary and sole geometric descriptor. This paper inspires a new generation of constitutive laws for crushable granular continua with arbitrary, yet evolving, particle morphology. We explore the idea of introducing multiple grain shape descriptors into Continuum Breakage Mechanics (CBM), a theory originally designed to track changes in particle size distributions during confined comminution. We incorporate the influence of these descriptors on the elastic strain energy potential and treat them as dissipative state variables. In analogy with the original CBM, and in light of evidence from extreme fragmentation in nature, the evolution of the additional shape descriptors is postulated to converge towards an attractor. Comparisons with laboratory experiments, discrete element analyses and particle-scale fracture models illustrate the encouraging performance of the theory. The theory provides insights into the feedback among particle shape, compressive yielding and inelastic deformation in crushable granular continua. These results inspire new questions that should guide future research into crushable granular systems using particle-scale imaging and computations.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Validation of Fault Displacements from Dynamic Rupture Simulations against the Observations from the 1992 Landers Earthquake

ABSTRACT Coseismic fault displacements in large earthquakes have caused significant damage to structures and lifelines on and near fault lines. Coseismic displacements represent a real threat, especially to distributed infrastructure systems. For infrastructure systems that can not avoid active faults, engineering displacement demands are defined using probabilistic fault-displacement hazard analyses (PFDHA). However, PFDHA models are sparse and poorly constrained partly due to the scarcity of detailed fault-displacement observations. Advancements in dynamic rupture simulation methods make them an attractive approach to address this important issue. Because fault displacements can be simulated for various geologic conditions as constrained by current knowledge about earthquake processes, they can be used to supplement the observation datasets. In addition to providing on-fault displacements, when used with appropriate constitutive models for the bulk medium, they can capture off-fault distributed inelastic deformations as well. For viable extrapolation, simulations must first be validated against data. In this article, we summarize the calibration and validation of the dynamic rupture model against the observations of the well-documented 1992 Landers earthquake. We defined a preferred model that reproduces several first-order fault-displacement metrics such as the on-fault partition of the total displacement, the mean fault-zone width, and the location of the peak displacement. Simulated ground motions consistent with the observations ensure that all physics important to modeling have been properly parameterized. For the extrapolation, we generated a suite of dynamic rupture models to quantify expected fault-displacement metrics, their intercorrelations, and magnitude dependencies, which are in part supported by the Landers and other recent earthquakes. Our validation and extrapolation exercise paves the way for using dynamic rupture modeling to quantitatively address fault-displacement hazard on a broader scale. The results are promising and are expected to be useful to inform PFDHA model development.

Geochemistry & Geophysics↗

DOE OSTI2022

Analysis of Microwave Resonant Cavity Transducer Degradation Mechanisms

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Second Annual Report on Development of Microwave Resonant Cavity Transducer for Fluid Flow Sensing

We are investigating a microwave resonant cavity transducer for flow sensing in the vessel of a high temperature fluid advanced reactor (AR), such as a molten salt cooled reactor (MSCR) or a sodium fast reactor (SFR). This transducer is a hollow metallic cylindrical cavity, with the flat wall of the cylinder flexible enough to undergo microscopic deflection due to dynamic fluid pressure. Membrane deflection leads to a shift in the resonant frequency, which can be detected with a spectrum analyzer. Because the transducer is intended for immersion in a high temperature corrosive fluid, understanding of material degradation is crucial for estimation of transducer performance lifetime, and development of measurement interpretation algorithms. We conducted a preliminary computational investigation of relevant damage mechanisms of a stainless steel 316 cylindrical resonator in FLiBe salt. The two main damage mechanisms, creep and corrosion, were modeled using multiphysics COMSOL software. Degradation was modeled for a temperature range 500°C to 700°C. Coupling of the damage mechanisms was not considered. These models predict significant inelastic deformation at most temperatures due to creep, and qualitatively predict chromium depletion both along the liquid/solid interface and along the grain boundaries. An algorithmic approach for compensation of these degradation effects during fluid flow measurements will be developed in the future work. To validate sensor physics, we have performed proof-of-principle test of flow sensing in water. For this test, we have developed a cylindrical resonator for K-band, which was machined from brass. The cavity was excited through WR-42 waveguide through a subwavelength hole on the side of the wall of the cylinder. To increase the spectral signal visibility, we developed a signal processing method for baseline subtraction. A flow loop for proof-of-principle test of transducer performance in water was assembled. A commercial flow meter was installed in the loop for reference measurements. Cylindrical cavity was excited in the TE 011 mode with resonant frequency f ≈ 17.8GHz. Frequency shift of cavity spectral response was obtained by gradually increasing water flow rate from 0 to 60gpm. Corresponding monotonic increase of resonant frequency shift by several MHz was observed. Approximate figure of merit of sensitivity to flow rate is 100KHz/GPM. In addition, we have identified an existing liquid sodium experimental setup for demonstration of flow sensing in environment similar to that of an advanced reactor. The setup consists of a cylindrical vessel and center feed line, where transducer inserted through the lid will measure velocity of the impinging liquid jet. As a calibration experiment, we have assembled a water vessel with center feed with the same dimensions as those of the liquid sodium setup. We have also developed and insertion probe consisting of a 50cm brass waveguide enclosed in protective SS316 tube. Using the water loop, we have demonstrated feasibility of sensing the impinging liquid jet in the vessel.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

A Multiphysics Multiscale Simulation Platform for Damage, Environmental Degradation, and Life Prediction of CMCs in Extreme Environments

This project successfully developed a multiphysics, multiscale computational framework to enhance the design and development of CMCs, with a focus on modeling highly nonlinear, time-dependent damage mechanisms and material degradation under extreme conditions, such as those experienced in turbine service environments. The project made significant advances in improving our understanding of progressive damage, oxidative degradation, and time-dependent inelastic deformation in CMCs, with particular attention to the role of uncertainties in predictions. Key outcomes include the integration of advanced material characterization, uncertainty quantification, and multiphysics constitutive models to predict the behavior of CMCs over their service life. A novel multiscale methodology was employed, which integrated microscale constituent behaviors with structural-scale responses, enabling the manufacturing defects in the microstructure that are prone to damage nucleation. Through the development of DL algorithms, the project advanced the prediction of damage initiation and crack propagation, taking into account the defect morphology and statistical variations across multiple scales. The framework was rigorously validated using thermomechanical experiments, which tested CMCs under various mechanical loadings at elevated temperatures, further enhancing the model's predictive capability. Overall, the research outcomes have provided a more accurate, reliable method for predicting CMC component life, significantly advancing material design, and improving component reliability in extreme environments. This work has strong implications for the optimization of turbine components and other high-performance applications where CMCs are used.

03 NATURAL GAS↗

DOE OSTI2025

Physics-Informed Machine Learning Model for Ceramic Matrix Composite Creep

A physics-informed recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear, time-dependent constitutive behavior of ceramic matrix composites (CMCs) driven by matrix damage and constituent creep at the microscale. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the high-fidelity generalized method of cells (HFGMC) approach which calls appropriate creep and damage models for each of the constituents. This coupling permits simulating the nonlinear behavior of CMCs based on constituent response at the microscale along with microstructural features such as fiber and porosity volume fraction and fiber radius. The microscale repeating unit cell is loaded under creep fatigue conditions to replicate the material loading experienced in a turbine engine. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input stress sequence, temperature, and microstructural features, the resulting strain history response while satisfying physical constraints related to creep rate, isochoric inelastic deformation, and strain energy density. The trained surrogate model is shown to effectively match the strain history over quantified distributions of microstructural features and relevant loading regimes and temperatures. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore, the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex microstructures, nonlinear time-dependent material response, and under non-monotonic loading conditions.

ceramic matrix composites↗