SEARCH · Search NASA
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.
An elastic potential is proposed that is capable of modeling the reversible portion of the observed nonlinear response of unidirectional graphite fiber composites. The model includes both the stiffening stress-strain behavior as well as the softening Poisson's response for loading in the fiber direction. The model is compared with experimental results for Celion 6000/PMR-15 graphite-polyimide.
In this paper, we introduce a data-driven machine learning approach for modeling one-dimensional stress–strain behavior under cyclic loading, utilizing experimental data from the nickel-based Alloy 617. The study employs uniaxial creep–fatigue test data acquired under various loading histories and compares two distinct neural network-based ODE models. The first model, known as the black-box model, comprehensively describes the strain–stress relationship using a Neural ODE equation. To interpret this black-box model, we apply the Sparse Identification of Nonlinear Dynamical Systems (SINDy) technique, transforming the black-box model into an equation-based model using symbolic regression. The second model, the Neural flow rule model, incorporates Hooke’s Law for the linear elastic component, with the nonlinear part characterized by a Neural ODE. Both models are trained with experimental data to accurately reflect the observed stress–strain behavior. We conduct a detailed comparison with the standard Chaboche model, which includes three back stresses. Our results demonstrate that the neural network-based ODE models precisely capture the experimental creep–fatigue mechanical behavior, exceeding the standard Chaboche model’s accuracy. Furthermore, an interpretable model derived from the black-box neural ODE model through symbolic regression achieves accuracy comparable to the Chaboche model, enhancing its interpretability. The results highlight the potential of neural network-based ODE models to depict complex creep–fatigue behavior, eliminating the necessity for experts to define a specific, material-focused model form.
A computerized analysis of the nonlinear behavior of fibrous composite laminates including axial loading, thermal loading, temperature dependent properties, and edge effects is presented. Ramberg-Osgood approximations are used to represent lamina stress-strain behavior and percent retention curves are employed to model the variation of properties with temperature. Balanced, symmetric laminates comprised of either boron/epoxy, graphite/epoxy, or borsic-aluminum are analyzed using a quasi-three-dimensional finite element analysis. Results are presented for the interlaminar stress distributions in cross-ply, angle-ply, and more complex laminates. Nonlinear stress-strain curves for a variety of composite laminates in tension and compression are obtained and compared to other existing theories and experimental results.
A study was made of the stress-strain response of several hybrid laminates, and the damage was correlated with nonlinear stress-strain response and ultimate strength. The fibers used in the laminates were graphite, S-glass, and Kevlar. Some laminates with graphite fibers had perforated Mylar film between plies, which lowered the interlaminar bond strength. The laminate configurations were chosen to be like those of buffer strips in large panels and fracture coupons. Longitudinal and transverse specimens were loaded in tension to failure. Some specimens were radiographed to reveal damage due to edge effects. Stress-strain response is discussed in terms of damage shown by the radiographs. Ultimate strengths are compared with simple failure criteria, one of which account for damage.
Equations are presented and described which can be used to predict bounds on the tensile and flexural strengths of nongraphitic superhybrid (NGSH) composites. These equations are derived by taking into account the measured stress-strain behavior, the lamination residual stresses and the sequence of events leading to fracture. The required input for using these equations includes constituents, properties (elastic and strength), NGSH elastic properties, cure temperature, and ply stress influence coefficients. Results predicted by these equations are in reasonably good agreement with measured data for strength and for the apparent knees in the nonlinear stress-strain curve. The lower bound values are conservative compared to measured data. These equations are relatively simple and are suitable for use in the preliminary design and initial sizing of structural components made from NGSH composites.
Using standardized uniaxial tensile test specimens is not practical due to the limited volumes in irradiation capsules, molten salt degradation facilities, or oxidation apparatus. The ASTM International Standard Test Method for Tensile Strength Estimate by Disc Compression of Manufactured Graphite (ASTM D8289) was developed to provide a convenient way for estimating tensile strength. Unlike the traditional uniaxial tensile test in ASTM C749 (International Standard Test Method for Tensile Stress–Strain of Carbon and Graphite), which uses dog-bone shaped specimens larger than 12.95 mm × 120.65 mm, the ASTM D8289 standard uses smaller discs with diameters of 6–12.7 mm. A digital image correlation (DIC) system, which uses a full-field noncontact surface displacement measurement technique, was applied along with the ASTM D8289 disc splitting test on IG-110 samples and compared with previous measurements from Mersen 2114 samples. Results confirmed that the DIC technique can measure the surface displacement/strain on these small (Ø6 mm × 3 mm) graphite specimens with good repeatability. However, Mersen 2114 and IG-110 samples exhibited strain discrepancies when DIC measurements were compared with analytical and finite element simulation values. The loading history and strain results also indicated different mechanical behaviors between Mersen 2114 and IG-110, particularly the nonlinear behavior of the IG-110 samples. Good agreement was observed by comparing the splitting tensile strengths of two superfine-grained grades with results from other work. The specimen size effect is discussed when comparing the splitting tensile strength with corresponding uniaxial tensile strength of these two graphite grades.
Equations are presented and described which can be used to predict bounds on the tensile and flexural strengths of nongraphitic superhybrid (NGSH) composites. These equations are derived by taking into account the measured stress-strain behavior, the lamination residual stresses and the sequence of events leading to fracture. The required input for using these equations includes constituent composites, elastic and strength properties, NGSH elastic properties, cure temperature, and ply stress influence coefficients. Results predicted by these equations are in reasonably good agreement with measured data for strength and for the apparent 'knees' in the nonlinear stress-strain curve. The lower bound values are conservative compared to measured data. These equations are relatively simple and should be suitable for use in the preliminary design and initial sizing of structural components made from NGSH composites.
The nonlinear hysteretic behavior of rocks under cyclic loading is a crucial area of study in geomechanics. The macroscopic response of a variety of materials has been found to be contingent upon the behavior of the micro-scale structure. This project aims to develop a functional and maintainable software package for generating a multi-phase numerical mesh and accompanying simulation files for finite element modeling used in computational mechanics solvers. Meshes generated from images often lack key preprocessing that reduces noise and prevents mesh element distortion that can increase computational cost. By incorporating user feedback throughout, grain2mesh ensures a high-fidelity mesh that can be used to model grain-scale interactions such as shearing, crack propagation, and interfacial material contrast. Scientific applications of this software include material fracturing, stress-strain analysis for natural and engineered materials, and nonlinear meso-scale analysis.
The tensile stress-strain behavior of five types of boron/aluminum laminates was investigated. Longitudinal and transverse stress-strain curves were obtained for monotonic loading to failure and for three cycles of loading to successively higher load levels. The laminate strengths predicted by assuming that the zero deg plies failed first correlated well with the experimental results. The stress-strain curves for all the boron/aluminum laminates were nonlinear except at very small strains. Within the small linear regions, elastic constants calculated from laminate theory corresponded to those obtained experimentally to within 10 to 20 percent. A limited amount of cyclic loading did not affect the ultimate strength and strain for the boron/aluminum laminates. The laminates, however, exhibited a permanent strain on unloading. The Ramberg-Osgood equation was fitted to the stress-strain curves to obtain average curves for the various laminates.
Designing fiber-reinforced polymer composites (FRPCs) with a tailored nonlinear stress-strain response is crucial for applications such as energy absorption in crash structures, flexible robotics, and impact-resistant protective gear. However, the inherent complexities of composite materials and the multitude of parameters involved, render traditional design and optimization methods inadequate for achieving effective inverse design of composites. In this paper, we present an AI-based inverse design framework that effectively and efficiently generates FRPCs with targeted nonlinear stress-strain responses. We introduce a physically constrained diffusion model (PC3D_Diffusion) capable of managing the complexities of composite materials and producing detailed, high-quality designs. We propose a loss-guided, learning-free approach to generate physically feasible microstructure designs by explicitly enforcing physical constraints during the generation process. For training purposes, 1.35 million FRPC samples were created, and their corresponding stress-strain curves were computed using established physics-based computational models. The results show that PC3D_Diffusion consistently generates high-quality designs with tailored mechanical behaviors, while guaranteeing compliance with the physical constraints. PC3D_Diffusion advances FRPC inverse design and may facilitate the inverse design of other 3D materials, offering potential applications in industries reliant on materials with custom mechanical properties.
Water released from hydrous minerals in subducting slabs reacts with the overlying plate, resulting in widespread serpentinization in the mantle wedge. Deformation of serpentinized peridotites has been invoked to explain forearc seismic anisotropy, yet studies of the mechanical properties and deformation behaviors of serpentine-bearing multiphase aggregates remain limited. Here we deformed olivine-antigorite mixtures containing 70, 50, and 20 vol.% of antigorite at 2.5 – 7.6 GPa, 673 K and strain rates of ∼10 –5 –10 –4 s –1 . Elasto-viscoplastic self-consistent simulations, constrained by synchrotron X-ray diffraction (XRD) data, were used to estimate lattice strain, stress–strain partitioning, crystallographic preferred orientations (CPO) development, and aggregate strength. Selected run products were also analyzed by electron backscatter diffraction for comparison with the CPO results obtained from XRD experiments. We found olivine transitions from A- or B-type to C-type when antigorite fraction drops to 20 vol.%, coinciding with a microstructural change from interconnected weak layers to a load-bearing framework (LBF). An additional run on a sample Atg50/Ol50 with preexisting microstructures suggested the formation of LBF was promoted by these microstructures, although the preexisting antigorite CPO has been overprinted at 20.8 % strain and could be erased completely by subsequent deformation in nature. Estimated viscosity of the two-phase mixtures suggests that low-degree serpentinization (≤20 %) in the mantle wedge may increase the strength of olivine-rich peridotite and hinder slab-mantle decoupling, whereas high-degree serpentinization (≥50–70 %) weakens the peridotite and favors decoupling if sufficient viscosity contrast (>10) develops. Seismic anisotropy shows a nonlinear dependence on antigorite fraction: antigorite CPO governs the anisotropy of the mixtures with ≥50 vol.% antigorite, whereas olivine CPO dominates at low fractions (∼20 vol.%). The presence of pre-existing microstructures reduces seismic anisotropy of the deformed mixtures, however the persistence of pre-existing CPO in actively subducting slabs remains uncertain, making their significance over geological timescales questionable.
Understanding the time-dependent effects of charged vacancies on the electromechanical responses of materials is at the forefront of research for designing materials exhibiting metal-insulator transitions and memristive behavior. A Rayleighian approach is used to develop a model for studying the nonlinear kinetics of the reaction leading to generation of vacancies and electrons via the dissociation of vacancy-electron pairs. Also, diffusion and elastic effects of charged vacancies are considered to model polarization-electric potential and strain-electric potential hysteresis loops. The model captures multiphysics phenomena by introducing couplings among polarization, the electric potential, stress, strain, and concentrations of charged (multivalent) vacancies and electrons (treated as classical negatively charged particles), where the concentrations can vary due to association-dissociation reactions. A derivation of coupled time-dependent equations based on the Rayleighian approach is presented. Three limiting cases of the governing equations are considered, highlighting the effects of (1) nonlinear reaction kinetics on the generation of charged vacancies and electrons, (2) Vegard's law (i.e., the concentration-dependent local strain) on asymmetric strain-electric potential relations, and (3) coupling between a fast component and the slow component of the net polarization on the polarization-electric-field relations. The Rayleighian approach discussed in this work should pave the way for developing a multiscale modeling framework in a thermodynamically consistent manner while capturing multiphysics phenomena in ferroelectric materials. Published by the American Physical Society 2025
Architected materials derive their properties from the geometric arrangement of their internal structural elements. Their designs rely on continuous networks of members to control the global mechanical behavior of the bulk. Here, in this study, we introduce a class of materials that consist of discrete concatenated rings or cage particles interlocked in three-dimensional networks, forming polycatenated architected materials (PAMs). We propose a general design framework that translates arbitrary crystalline networks into particle concatenations and geometries. In response to small external loads, PAMs behave like non-Newtonian fluids, showing both shear-thinning and shear-thickening responses, which can be controlled by their catenation topologies. At larger strains, PAMs behave like lattices and foams, with a nonlinear stress-strain relation. At microscale, we demonstrate that PAMs can change their shapes in response to applied electrostatic charges. The distinctive properties of PAMs pave the path for developing stimuli-responsive materials, energy-absorbing systems, and morphing architectures.
The behavior of finite element models employing different constitutive relations to describe the stress-strain behavior of soils is investigated. Three models, which assume small strain theory is applicable, include a nondilatant, a dilatant and a strain hardening constitutive relation. Two models are formulated using large strain theory and include a hyperbolic and a Tresca elastic perfectly plastic constitutive relation. These finite element models are used to analyze retaining walls and footings. Methods of improving the finite element solutions are investigated. For nonlinear problems better solutions can be obtained by using smaller load increment sizes and more iterations per load increment than by increasing the number of elements. Suitable methods of treating tension stresses and stresses which exceed the yield criteria are discussed.
The application of finite element techniques to the analytic representation of the nonlinear behavior of arbitrary two-dimensional bodies containing cracks is discussed. Specific methods are proposed using which it should be possible to obtain information concerning: the description of the maximum, minimum, and residual near-tip stress and strain fields; the effects of crack closure on the near-tip behavior of stress and strain fields during cyclic loading into the plastic range; the stress-strain and displacement field behavior associated with a nonstationary crack; and the effects of large rotation near the crack tip.
Static stiffness, strength and ultimate strain after thermal cycling were investigated for graphite/Kevlar 49/epoxy and graphite/S-glass/epoxy angle-ply laminates. Tensile stress-strain curves to failure and uniaxial tensile properties were determined, and theoretical predictions of modulus, Poisson's ratio and ultimate strain, based on linear lamination theory, constituent ply properties and measured strength, were made. No significant influence on tensile stress properties due to stacking sequence variations was observed. In general, specimens containing two 0-degree Kevlar or S-glass plies were found to behave linearly to failure, while specimens containing 4 0-degree Kevlar or S-glass plies showed some nonlinear behavior.
Three-dimensional, nonlinear, finite element structural analyses were performed for a simulated aircraft combustor liner specimen in order to assess the capability of nonlinear analyses using classical inelastic material models to represent the thermoplastic-creep response of the component. In addition, the computed stress-strain history at the critical location was input into life prediction methods in order to evaluate the ability of these procedures to predict crack initiation life. It is concluded that: (1) elastic analysis is adequate for obtaining strain range and critical location; (2) inelastic analyses did not accurately represent cyclic behavior of materials; and (3) none of the crack initiation life prediction methods were satisfactory.
Strain redistribution corrections were developed for a simplified inelastic analysis procedure to economically calculate material cyclic response at the critical location of a structure for life prediction proposes. The method was based on the assumption that the plastic region in the structure is local and the total strain history required for input can be defined from elastic finite-element analyses. Cyclic stress-strain behavior was represented by a bilinear kinematic hardening model. The simplified procedure predicts stress-strain response with reasonable accuracy for thermally cycled problems but needs improvement for mechanically load-cycled problems. Neuber-type corrections were derived and incorporated in the simplified procedure to account for local total strain redistribution under cyclic mechanical loading. The corrected simplified method was used on a mechanically load-cycled benchmark notched-plate problem. The predicted material response agrees well with the nonlinear finite-element solutions for the problem. The simplified analysis computer program was 0.3% of the central processor unit time required for a nonlinear finite-element analysis.