Search NASA⌕ Search

SEARCH · Search NASA

Results for “representative volume element”

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 19 records

Observation and Modeling of Dynamic Fracture Behaviors of Battery Cell Under Impact Loading Using Enhanced Representative Volume Element Concept

The burgeoning electric automobile industry has increased interest in battery safety. Battery cells experience significant mechanical stress during operation, including the impact of accidents and vibrations from driving. The potential for thermal runaway reactions in battery cells raises safety concerns. Although numerous researchers have defined the dynamic behavior of battery cells and proposed numerical models to describe it, few studies have focused on the high-strain rate mechanical impact phase correlated with the onset of fracture. In this study, we describe the dynamic behavior of pouch battery cells and propose a modeling method to study their mechanical failure under impact situations. Impact tests are conducted at various velocities and heights. To overcome numerical issues commonly encountered under rapid deformation scenarios, a new finite element model is developed based on the representative volume element model. The proposed approach efficiently simulates continuous crack propagation and brittleness behavior during impact by permitting the individual behavior of the cell components. Therefore, engineers can reliably design safer electric vehicle battery cells by measuring the properties of the cell components.

ENERGY STORAGE↗

Geometric Requirements for Modeling Polycrystalline Representative Volume Elements Using the Generalized Method of Cells

To accurately predict the mechanical response of advanced metallic parts, the grain structure of the metal must be considered and coupled to the macroscale in a multiscale model. The size of the microscale representative volume element must be sufficiently large, and the grid used to discretize the microscale must be adequately refined, in order to ensure an accurate and converged solution. The generalized method of cells is used to predict the mechanical properties of polycrystalline representative volume elements of copper containing randomly oriented grains with the properties of a single face-centered-cubic crystal. Moreover, the results of a parametric study are used to conclude the geometric requirements for an objective microscale model, with respect to the elastic stiffness.

Multiscale Modeling↗

Modeling the Stochastic Response of Fiber Reinforced Composites with Varied Representative Volume Element Sizes

Fiber reinforced composites are desirable in applications where low weight and high strength are needed, but are susceptible to variability and flaws during manufacturing, making failure predictions difficult. These flaws may occur at the microscale where mechanical properties vary locally due to regions of fiber clusters and matrix pockets. In this study, a multiscale approach was taken to model 3-point bend, 4-point bend, and tensile experiments of a unidirectional composite from only having microstructure scans of these samples and constituent properties from literature. These scans were sampled with different sized windows, and statistically equivalent microstructures were generated, then simulated for stiffness, strength, and fracture toughness using a reduced order micromechanical model and NASA’s Multiscale Analysis Tool (NASMAT). Mesoscale models were created with equivalent element sizes to microstructures and properties sampled from microscale simulation results. Results showed how microscale size affects certain mechanical properties. Also shown is how well mesoscale models agree to experiments when using stochastic element properties and varying element size.

microstructure↗

An Algorithm for Characterization of Fiber Aggregation in Composite Microstructures

Composite structures are susceptible to localized flaws, or variability, that drive global failure. To capture this variance, a multiscale model needs to be introduced that not only accurately represents the statistical nature of the composite microstructure but is also efficient. At the microscale, fiber aggregation creates local stress concentrations, where failure is likely to occur sooner than expected. A method of rapid microstructure generation, in conjunction with cluster characterization, is examined to develop accurate reproductions of 2D composite cross-sections. Parameterization of shape and size of fiber clusters is used to characterize representative volume elements. The discrete element method is used to generate pseudo-microstructures. The clustering parameters from the pseudo- and actual microstructure arrangements (obtained from micrographs) can be compared to determine the validity of the representative volume element generated with the discrete element method. Results show a promising approach to evaluating randomness of fiber distributions and how to accurately recreate microstructures for strength analysis.

Carbon fiber↗

Micromechanics Modeling of Functionally Graded Interphase Regions in Carbon Nanotube-Polymer Composites

The effective elastic properties of a unidirectional carbon fiber/epoxy lamina in which the carbon fibers are coated with single-walled carbon nanotubes are modeled herein through the use of a multi-scale method involving the molecular dynamics/equivalent continuum and micromechanics methods. The specific lamina representative volume element studied consists of a carbon fiber surrounded by a region of epoxy containing a radially varying concentration of carbon nanotubes which is then embedded in the pure epoxy matrix. The variable concentration of carbon nanotubes surrounding the carbon fiber results in a functionally graded interphase region as the properties of the interphase region vary according to the carbon nanotube volume fraction. Molecular dynamics and equivalent continuum methods are used to assess the local effective properties of the carbon nanotube/epoxy comprising the interphase region. Micromechanics in the form of the Mori-Tanaka method are then applied to obtain the global effective properties of the graded interphase region wherein the carbon nanotubes are randomly oriented. Finally, the multi-layer composite cylinders micromechanics approach is used to obtain the effective lamina properties from the lamina representative volume element. It was found that even very small quantities of carbon nanotubes (0.36% of lamina by volume) coating the surface of the carbon fibers in the lamina can have a significant effect (8% increase) on the transverse properties of the lamina (E22, k23, G23 and G12) with almost no affect on the lamina properties in the fiber direction (E11 and v12).

Seidel, Gary D.↗

Graph neural networks for mechanical property prediction of 2D fiber composites

This work investigates the ability of graph neural networks (GNNs) to homogenize 2D fiber composite microstructures. We use different inhomogeneity and anisotropy indices to motivate and show that the Volume Elements (VEs) used in ML methods should ideally be far from their Representative Volume Element (RVE) size limit and, consequently, are notably anisotropic. Hence, training only the isotropic limit properties may not be acceptable. Another aspect is the need to normalize elastic stiffness values for ML, especially when high elastic contrast ratios are encountered between composite phases or in the material set. We introduce a normalization technique based on the mean-field method (MFM) to handle such high contrast ratios and train for the entire stiffness tensor. We show that the proposed GNN approaches exhibit high accuracy and efficiency compared to traditional methods and convolutional neural networks, utilizing unstructured graphs constructed from microstructure topology. Our model successfully predicts the stiffness tensor, peak strength under bulk damage, and brittle fracture initiation strength across diverse microstructure configurations while maintaining high accuracy even for extreme material contrasts and volume fractions. We also present a method to improve prediction accuracy for small dataset sizes using Voronoi partitioning.

Brittle strength↗

A "Tread" Ahead of the Competition

Through a consortium formed by NASA's Glenn Research Center and the Great Lakes Industrial Technology Center, Goodyear acquired a software use agreement for NASA's Micromechanics Analysis Code (MAC), in an effort to design stronger and safer tires. Based on the Generalized Method of Cells micromechanics methodology, the MAC technology has allowed Goodyear scientists and engineers to perform structural analysis of composite laminates for tires all in one step, rather than the several steps previously required. Currently, the MAC code can be used to assess and optimize composite laminates; set cord-spacing guidelines (the diameter and distance between cords may impact the overall durability of a tire); explain structural impact of geometrical configurations; and compare the behavior of the different wire constructions that make up a tire. Furthermore, the MAC software contains a built-in material database and several representative volume elements for a wider choice of composites to better represent the plies and belts in a tire. The end result of the process ultimately reduces the time spent on building, testing, and adjusting tires.

Source record↗

On the Representativity of Electrode Microstructure Parameters and Their Electrochemical Response for Lithium Ion Batteries

Lithium-ion battery electrochemical models require an accurate description of the electrodes microstructures to be predictive that can be achieved through nanoscale imaging. Such observations are however limited by their field of view (FOV), as they provide only a subset of the whole electrode volume that does not necessarily represent the whole electrode microstructure heterogeneity, and therefore can bias the microstructure analysis. A microstructure scale electrochemical model was used to investigate lithium plating onset, material non-uniform utilization, and in-plane heterogeneities for an NMC-graphite full cell. To evaluate the representativeness, and thus relevance, of these model predictions, a coupled representativity analysis has been performed on the microstructure parameters and, in a novel way, on the full cell electrochemical response. Electrode microstructure parameters representativeness has been first quantified using the representative volume element (RVE) methodology. The RVE major flaw is that ultimately it can only conclude if a FOV contains representative subvolumes of the FOV, but not if the FOV itself is representative of the electrode volume. Analysis can conclude negatively ('FOV is not representative'), but not positively ('FOV is representative'). One major contribution of this work was to quantify the convergence of the RVE size with the FOV, to actually investigate the FOV representativeness and thus partly remedy this intrinsic limitation. The analysis determined that performing a standard RVE calculation, without exploring its FOV convergence, is likely to strongly underestimate the actual RVE size. The new RVE methodology has been automated in the NREL open-source Microstructure Analysis Toolbox (MATBOX) and is available to the battery community. Representativeness of microstructure parameters is however only an intermediate step, as the end-results of an electrochemical model are performances predictions. Indeed, what is the practical consequence of a given deviation for a microstructure parameter? The microstructure parameter deviation propagations to the 3D microstructure scale electrochemical response have been then quantified for different charge rates. This defines a threshold for the microstructure parameters FOV for a desired maximum deviation of the electrochemical response. Such deviation propagation analysis is analogous to error propagation analysis and is necessary to determine the relevance of microstructure scale model predictions for macroscale predictions. Electrochemical model shows cell representative section areas are increasing with C-rate, due to higher in-plane heterogeneities, indicating larger FOVs are required specifically for fast charge modeling. Therefore, we introduced the novel concept of electrochemical RVE (eRVE) that is a function of the operating conditions (thus defined as a dynamic RVE), with an increasing dependence with the C-rate. Representativity analysis of the investigated cell determined a FOV of 144.4 x 54.4 m2 is large enough to establish a convergence on the representative section areas for low to intermediate C-rate (=2.5C), but not large enough to conclude for higher rates. This work aims to emphasize the importance of representativity analysis for LIB electrode microstructures, as it is required to estimate the error, and thus the relevance, of microstructure parameters intended to be used in macroscale models. The methodology and results can help researchers to select the relevant imaging and associated FOV required to provide accurate enough microstructure parameters.

ADVANCED PROPULSION SYSTEMS↗

The effect of weak interface on transverse properties of a ceramic matrix composite

Experimental studies conducted at NASA Lewis on SiC reaction-bonded Si3N4 composite system showed that transverse stiffness and strength were much lower than those predicted from existing analytical models based on good interfacial bonding. It was believed that weakened interfaces were responsible for the decrease in tranverse properties. To support this claim, a two-dimensional FEM analysis was performed for a transverse representative volume element. Specifically, the effect of fiber/matrix displacement compatibility at the interface was studied under both tensile and compressive transverse loadings. Interface debonding was represented using active gap elements connecting the fiber and matrix. The analyses show that the transverse tensile strength and stiffness are best predicted when a debonded interface is assumed for the composite. In fact, the measured properties can be predicted by simply replacing the fibers by voids. Thus, it is found that little or no interfacial bonding exists in the composite, and that an elastic analysis can predict the transverse stiffness and strength.

Shimansky, R. A.↗

A Convolutional Neural Network for Multiscale Modeling of Composite Materials

An artificial convolutional neural network was created to efficiently mimic a micromechanics model, the High Fidelity Generalized Method of Cells, for use in multi-scale structural finite element analysis. The network was found to quickly and accurately replicate the stiffness predicted by the micromechanics model using a 2D image of an idealized representative volume element of a fiber/matrix microstructure. The long-term goal of this work is to efficiently apply multi-scale methods for predicting the damage progression of a composite structure.

Composites↗

Simulated effect of defect volume and location on very high cycle fatigue of laser beam powder bed fused AlSi10Mg

This study quantifies the interaction between volumetric defect location and size on the very high cycle fatigue (VHCF) of laser beam powder bed fused (LB-PBF) AlSi10Mg. Crystal plasticity finite element method (CPFEM) simulations were used to investigate the effects of defect location and size on the driving force for crack initiation. The CPFEM model was calibrated against uniaxial and cyclic experimental data of LB-PBF AlSi10Mg. Defect characteristics were informed by experimental data from the specimens produced in various geometries to create realistic representative volume elements (RVEs) with equivalent volume fractions of defects. By embedding defects of varying sizes and locations within the RVEs, fatigue indicator parameters (FIPs) were calculated to analyze the impact of defects’ characteristics on fatigue performance. Different combinations of defect volume and locations were generated for various microstructure instantiations, providing insight into extreme value fatigue responses. Larger defect volumes located on free surfaces consistently generated the highest FIPs, suggesting defect size and boundary proximity intensify stress concentration effects. RVEs with multiple smaller defects produced lower FIPs than those with single large critical defects. These findings underscore the critical role of defect characteristics on fatigue life, providing a foundation for future predictive modeling in fatigue-sensitive AM applications.

AlSi10Mg↗

Machine-Learning-Based Multiscale Methods for 3D Modelling of Granular Materials by Incorporating History-Dependent State Variables

Over the past decades, the prevalence of machine learning (ML) methods has made the development of ML-based constitutive models for granular materials undoubtedly a popular subject. Numerous studies have been made to feature the loading path or history-dependent stress-strain response of granular media using neural networks. In this work, a novel finite element method (FEM)–ML multiscale approach was developed by incorporating internal variables to improve the simulation accuracy of 3D history-dependent granular materials for the first time. To this end, a surrogate constitutive model based on the single-step-based multi-layer perceptron (MLP) neural network was used to replace representative volume element (RVE) simulations conducted by the discrete element method (DEM) in the multiscale FEM–DEM approach. Although the prediction principle of the MLP aligns with the FEM algorithm, artificially added internal variables are required to differentiate the loading history. To address this issue, history variables associated with the Frobenius norm are proposed to be fed into the MLP coupled with the strain tensor to extract the history-dependent behaviour of granular assemblies. The developed FEM–ML approach was demonstrated in 3D conventional triaxial compression (CTC) simulations. Compared to the multiscale FEM–DEM approach, the proposed FEM–ML method exhibits a significantly improved computational efficiency.

granular materials↗

Finite element based micro-mechanics modeling of textile composites

Textile composites have the advantage over laminated composites of a significantly greater damage tolerance and resistance to delamination. Currently, a disadvantage of textile composites is the inability to examine the details of the internal response of these materials under load. Traditional approaches to the study fo textile based composite materials neglect many of the geometric details that affect the performance of the material. The present three dimensional analysis, based on the representative volume element (RVE) of a plain weave, allows prediction of the internal details of displacement, strain, stress, and failure quantities. Through this analysis, the effect of geometric and material parameters on the aforementioned quantities are studied.

Glaessgen, E. H.↗

Seamlessly joining length scales: From atomistic thermal graphs to anisotropic continuum conductivity

Thermal transport in complex solids is governed by local structure, defects, and anisotropy, yet most continuum models still rely on oversimplified and homogenized conductivities. Here, we bridge atomistic and continuum descriptions by building finite element (FE) models directly from the site-projected thermal conductivity (SPTC), an atomic-level decomposition of the Green–Kubo thermal conductivity. We introduce a toolkit, the “Simulator Collection for Atomic-to-Continuum Scales (SCACS)”, which uses a graph neural network to predict SPTC on large atomic structures, coarse-grains these fields into anisotropic conductivity tensors, and embeds them into the heat-flow FE equation with a customized, anisotropy-aware adaptive mesh refinement scheme. Applied to silicon nanostructures, the resulting FE models act as representative volume elements, reproduce bulk conductivities, and capture interfacial and defect-driven anisotropy while maintaining thermodynamic consistency. Additionally, SCACS predicts experimental conductance trends and fields. This work demonstrates a general route for transferring atomistic transport information into device-scale thermal simulations with physics-based approximations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

3D Finite Element Analysis of Particle-Reinforced Aluminum

Deformation in particle-reinforced aluminum has been simulated using three distinct types of finite element model: a three-dimensional repeating unit cell, a three-dimensional multi-particle model, and two-dimensional multi-particle models. The repeating unit cell model represents a fictitious periodic cubic array of particles. The 3D multi-particle (3D-MP) model represents randomly placed and oriented particles. The 2D generalized plane strain multi-particle models were obtained from planar sections through the 3D-MP model. These models were used to study the tensile macroscopic stress-strain response and the associated stress and strain distributions in an elastoplastic matrix. The results indicate that the 2D model having a particle area fraction equal to the particle representative volume fraction of the 3D models predicted the same macroscopic stress-strain response as the 3D models. However, there are fluctuations in the particle area fraction in a representative volume element. As expected, predictions from 2D models having different particle area fractions do not agree with predictions from 3D models. More importantly, it was found that the microscopic stress and strain distributions from the 2D models do not agree with those from the 3D-MP model. Specifically, the plastic strain distribution predicted by the 2D model is banded along lines inclined at 45 deg from the loading axis while the 3D model prediction is not. Additionally, the triaxial stress and maximum principal stress distributions predicted by 2D and 3D models do not agree. Thus, it appears necessary to use a multi-particle 3D model to accurately predict material responses that depend on local effects, such as strain-to-failure, fracture toughness, and fatigue life.

Shen, H.↗

An evaluation of a coupled microstructural approach for the analysis of functionally graded composites via the finite-element method

A comparison is presented between the predictions of the finite-element analysis and a recently developed higher-order theory for functionally graded materials subjected to a thorough-thickness temperature gradient. In contrast to existing micromechanical theories that utilize classical (i.e., uncoupled) homogenization schemes to calculate micro-level and macro-level stress and displacement fields in materials with uniform or nonuniform fiber spacing (i.e., functionally graded materials), the new theory explicitly couples the microstructural details with the macrostructure of the composite. Previous thermo-elastic analysis has demonstrated that such coupling is necessary when: the temperature gradient is large with respect to the dimension of the reinforcement; the characteristic dimension of the reinforcement is large relative to the global dimensions of the composite and the number of reinforcing fibers or inclusions is small. In these circumstances, the standard micromechanical analyses based on the concept of the representative volume element used to determine average composite properties produce questionable results. The comparison between the predictions of the finite-element method and the higher-order theory presented herein establish the theory's accuracy in predicting thermal and stress fields within composites with a finite number of fibers in the thickness direction subjected to a thorough-thickness thermal gradient.

Pindera, Marek-Jerzy↗

Non-conformal interface-cohesive modeling with the shifted boundary method

The accurate simulation of boundary- and interface-dominated problems on complex geometries remains challenging when boundary- or interface-fitted meshes are difficult to generate, particularly for curved boundaries, polycrystalline microstructures, and dense interface networks. The Shifted Boundary Method (SBM) alleviates this meshing burden by shifting the enforcement of boundary conditions from the true boundary to a nearby surrogate boundary and recovering the effect of the true boundary through geometric correction terms, thereby enabling standard finite element spaces on non-boundary-fitted meshes. In this report, we develop a general shiftedboundary and shifted-interface framework within the open-source MOOSE framework. We first present a general SBM implementation for complex geometries on non-boundary-fitted meshes. We then adopt the Shifted Interface Method (SIM) for internal interfaces and develop a unified shifted-interface treatment in which the interface law is enforced on a surrogate interface and the effect of the true interface is recovered through shifted jumps, fluxes, and tractions. This perspective brings scalar thermal-contact and vector-valued cohesive-zone mechanics into a single framework, the latter realized as the Shifted Cohesive Zone Method (SCZM) and coupled with history-dependent constitutive models from NEML2. We further extend the MOOSE mesh infrastructure to support cohesive-zone calculations on distributed meshes. The framework is verified and demonstrated through three progressive studies: Poisson’s equation on a smoothed starshaped domain, a manufactured thermal-contact problem on a non-interface-fitted mesh, and a two-dimensional polycrystalline representative volume element combining crystal plasticity with cohesive grain-boundary interfaces. Across these studies, the shifted formulations reproduce boundary- and interface-fitted reference solutions with high fidelity, indicating that the proposed framework provides an accurate and efficient route to boundary- and interface-dominated simulations on arbitrary geometries without requiring fitted meshes.

Yang, Cheng-Hau↗