Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constitutive modeling”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Bayesian optimal experimental design for constitutive model calibration

Computational simulation is increasingly relied upon for high/consequence engineering decisions, which necessitates a high confidence in the calibration of and predictions from complex material models. However, the calibration and validation of material models is often a discrete, multi-stage process that is decoupled from material characterization activities, which means the data collected does not always align with the data that is needed. To address this issue, an integrated workflow for delivering an enhanced characterization and calibration procedure—Interlaced Characterization and Calibration (ICC)—is introduced and demonstrated. Further, this framework leverages Bayesian optimal experimental design (BOED), which creates a line of communication between model calibration needs and data collection capabilities in order to optimize the information content gathered from the experiments for model calibration. Eventually, the ICC framework will be used in quasi real-time to actively control experiments of complex specimens for the calibration of a high-fidelity material model. This work presents the critical first piece of algorithm development and a demonstration in determining the optimal load path of a cruciform specimen with simulated data. Calibration results, obtained via Bayesian inference, from the integrated ICC approach are compared to calibrations performed by choosing the load path a priori based on human intuition, as is traditionally done. The calibration results are communicated through parameter uncertainties which are propagated to the model output space (i.e. stress–strain). In these exemplar problems, data generated within the ICC framework resulted in calibrated model parameters with reduced measures of uncertainty compared to the traditional approaches.

42 ENGINEERING↗
DOE OSTI2022

A Micromorphic Length-Scale Coupling Framework for the Determination of Higher-Order Constitutive Models and the Multi-Scale Simulation of Heterogeneous Materials [Thesis]

Heterogeneous materials and materials with complex microstructures pose a unique challenge in the development of accurate models of their response to external stimuli. These difficulties principally arise due to the difficulty in characterizing and modeling the constituents and their interactions. It is usually possible, though non-trivial, to construct an explicit representation of the microstructure (a direct numeric simulation or DNS) but the method by which the complex modes of deformation and other processes can be homogenized to a reduced-order approximation is frequently unclear. Many homogenization approaches are ad-hoc and lack a strong justification beyond ease of computation. Furthermore, the homogenization approach can, in some cases, not utilize the full breadth of information available in the computation of the macro-scale stresses and deformation measures. It is also noteworthy that homogenization, by its very nature, will tend to obfuscate details of processes occurring at the lower length-scale. It is therefore of interest to include as much information as possible in the construction of the reduced-order model so as to be predictive in a variety of loading environments. We here present a length-scale bridging technique based upon the micromorphic continuum mechanics of Eringen which incorporates volume and surface area averages as a part of its construction. This approach enforces the balance equations at the micro-scale and then studies the effect of the spatially varying nature on the macro-scale. This leads, naturally, to further balance equations which are solved at the macro-scale. This work details a homogenization framework which arises naturally from the micromorphic construction of Eringen attempting to introduce no definitions beyond which are justifiable from micro-structural considerations. One of the results of this effort is the so-called “micromorphic filter” which has been developed and applied to several DNS to demonstrate its effectively. In order to determine the macroscopic degrees of freedom we utilize the special case of an overlap coupling technique where the macro-scale is fully constrained to the micro-scale. This enables us to study the resulting material properties as expected but also allows us to further study the boundary conditions on the additional degrees of freedom at the macro-scale.

42 ENGINEERING↗
DOE OSTI2024

Peridynamic neural operators: A data-driven nonlocal constitutive model for complex material responses

Neural operators, which can act as implicit solution operators of hidden governing equations, have recently become popular tools for learning the responses of complex real-world physical systems. Nevertheless, most neural operator applications have thus far been data-driven and neglect the intrinsic preservation of fundamental physical laws in data. Here, in this work, we introduce a novel integral neural operator architecture called the Peridynamic Neural Operator (PNO) that learns a nonlocal constitutive law from data. This neural operator provides a forward model in the form of state-based peridynamics, with objectivity and momentum balance laws automatically guaranteed. As applications, we demonstrate the expressivity and efficacy of our model in learning complex material behaviors from both synthetic and experimental data sets. We also compare the performances with baseline models that use predefined constitutive laws. We show that, owing to its ability to capture complex responses, our learned neural operator achieves improved accuracy and efficiency. Moreover, by preserving the essential physical laws within the neural network architecture, the PNO is robust in treating noisy data. The method shows generalizability to different domain configurations, external loadings, and discretizations.

42 ENGINEERING↗
DOE OSTI2022

Experimental study of compaction localization in carbonate rock and constitutive modeling of mechanical anisotropy

Abstract Sedimentary rocks are inherently anisotropic and prone to strain localization. While the influence of rock anisotropy on the brittle/dilative regime has been studied extensively, its influence on the ductile/compactive regime is much less explored. This paper discusses the anisotropic behavior of a high‐porosity carbonate rock from central Europe (the Maastricht Tuffeau). A set of triaxial tests with concurrent x‐ray tomography has been performed at different confining pressures. The anisotropic characteristics of this rock have been investigated by testing samples cored at different inclinations of the bedding, thus revealing non‐negligible effects of the coring direction on yielding and compaction behavior. Specifically, samples cored perpendicular to bedding display higher strength and longer stages of post‐yielding deformation before manifesting re‐hardening. Despite such alterations of the inelastic response, Digital Image Correlation has revealed that the strain localization mode is independent of the coring direction, thus being primarily affected by the confinement level. To capture the observed interaction between material anisotropy and compaction behavior at the continuum‐scale, an elastoplastic constitutive law has been proposed. For this purpose, a set of tensorial bases has been introduced to replicate how the oriented rock fabric modulates the yielding and plastic flow characteristics of the material. The analyses show that the impact of the coring direction on yield function and plastic flow rule is fundamentally different, thus requiring the use of distinct projection strategies (a strategy here defined heterotopic mapping). The performance of the model, studied through parametric analyses and by calibrating the experimental results, illustrates the improved capability of the proposed constitutive approach when applied to strongly anisotropic porous rocks.

Shahin, Ghassan↗
DOE OSTI2023

Simulations of biomass compression-screw feeding using a compressible non-Newtonian constitutive model

There is global interest in the conversion of biomass into sustainable low-carbon-footprint fuels and chemicals as an alternative to non-renewable fossil feedstocks. Feeding biomass solids into pressurized reactors is one of the key steps in biomass conversion. Predicting mechanical failure and energy requirements for this step helps avoid upstream processing bottlenecks and enables efficient operation of a biorefinery. Here, in this work, we developed a predictive computational model for biomass screw feeders that capture the highly viscous, non-Newtonian and compressible behavior of biomass slurries. Biomass compressible behavior is formulated by an equation of state and the non-Newtonian rheology is represented by a density-dependent viscosity model. Experimental data from two compression screw-feeder systems are presented as a validation for our model. Our model successfully predicted the location of the compressed biomass “plug”, biomass flow rate, and the required torque at different operating conditions for the experimental conditions studied in this work.

09 BIOMASS FUELS↗
DOE OSTI2025

A visco-plastic constitutive model for accurate densification and shape predictions in powder metallurgy hot isostatic pressing

Powder metallurgy hot isostatic pressing (PM-HIP) is an advanced manufacturing process that produces near net shape parts with high material utilization and uniform microstructures. Despite being used frequently to produce small-scale components, the application of PM-HIP to large-scale components is limited due to inadequate understanding of its complex mechanisms that cause unpredictable post-HIP shape distortions. A computational model can provide necessary information about the intermediate and final stages of the HIP process that can help understand it better and make accurate predictions. Generally, two types of computational models are employed for PM-HIP of metal powders, namely, plastic and visco-plastic models. Between these, the plastic model is preferred due to its cheaper calibration approach requiring less experimental data. However, the plastic model sometimes produces incorrect predictions when slight variations of the HIP conditions are encountered in practical situations. Therefore, this work presents a visco-plastic model that addresses these limitations of the plastic model. A novel modified calibration approach is employed for the visco-plastic model that utilizes less experimental data than existing approaches. With the new approach, the data requirement is same for both plastic and visco-plastic models. This also enables a quantitative comparison of plastic and visco-plastic models, which have been only qualitatively compared in the past. When calibrated with the same experimental data, both the models are found to produce similar results. In conclusion, the calibrated visco-plastic model is applied to several complex geometries, and the predictions are found to be in good agreement with experimental observations.

Hot isostatic pressing↗
DOE OSTI2022

A new woven composite constitutive model validated by shock wave experiments

In this paper, we present results of plate impact simulations of shock compressed woven glass fiber-reinforced plastic (GRP) performed using the Arbitrary Lagrangian–Eulerian three-dimensional finite element code. A hyperelastic large-strain-based empirical Continuum Damage Mechanics (CDM) formulation is employed to describe damage initiation and growth in the shock-compressed GRP. The model parameters calibration scheme utilizes the Velocity Interferometer System for Any Reflector normal particle velocity measurements at the free surface of the GRP target plates. The impact velocity in the experiments ranged from 8.5 to 418 m/s. The finite element model considered planar 0°/90° bidirectional plies with an individual ply thickness of 0.68 mm, stacked to reach a total laminate thickness of 6.8 mm. The anisotropic elastic strains were estimated from the experimentally determined tetragonal symmetry stiffness matrix for the GRP. The strain-based damage model captures several salient features observed in the measured free surface particle wave profiles, including the shock rise time, onset of Elastic—Elastic Cracking, and the shape of the nonlinear portion of the experimental particle velocity profiles. Furthermore, the CDM model predicts the dominant damage mode to be matrix microcracking due to shear and the associated bulk expansion (bulking) under the global compressive loading in the plate impact configuration.

42 ENGINEERING↗
DOE OSTI2023

Constitutive Modeling of the Cross-Linked Polymers During the Trio-Aging Mechanism: Effects of Water, Oxygen, UV

Abstract The synergized effects of diverse environmental conditions such as humidity, temperature, UV, and oxygen, are considered through a generic model by combining the concepts of network alteration, hygrothermal, and Photo-thermo oxidation Models [1–4]. Hygrothermal and Photo-oxidation are considered dual aging mechanisms which are due to the synergized effects of moisture-oxygen, and UV-oxygen, respectively at elevated temperatures. Here Trio-aging is considered as the synergized effects of Hygrothermal, and Photo-oxidation which the main environmental elements involved in polymer degradation are oxygen, water, and UV radiation where each of which has a different micro-structural impact on the polymer matrix. While oxygen imposes chain scission and cross-link formation, water, and UV radiation mainly cause chain scission. Similar to our recent models, the dual network hypothesis of Topolosky is utilized and as a result, the strain energy of the polymer matrix is written in terms of two states namely unaged and infinity. The infinity state network, itself, is decomposed into three networks, namely brittle, flexible, and dead networks. While the brittle network models the cross-link formation as sub-structural phenomena, the flexible and dead networks are utilized to represent the cross-link reduction and chain scission. The contribution of these networks in the infinity state is based on the amount of relevant environmental elements. The model is developed for a domain with plenty of oxygen and holds true for a gradual aging process that happens in super-thin samples.

Bahrololoumi, Amir↗