Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element simulation”

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

An assessment of the utility of multirate time integration for the modeling of laser powder bed fusion

Finite element simulation of the laser powder bed fusion process is made challenging by the inherently multiscale nature of the process. When using the typical global time stepping techniques, slowly-evolving regions of the domain receive the same numerical treatment as the regions with the highest temperature rates. The current work details the implementation and evaluation of an implicit multirate method which is able to advance different regions of the domain with distinct time steps, depending on their current solution rate. Previous work indicates that this representation of the temporal scales of the problem can lead to significant reductions in wall clock run time, and it is shown herein that multirate time integration, when used with a uniform process-scale mesh, can result in speedups of approximately 19, 42, and 87 times, for domains with edge lengths of approximately 1 mm, 2 mm, and 5 mm, respectively. When used in conjunction with $\ h$-refinement (limited to two levels of refinement), the resulting speedups (taken relative to the uniform mesh) are around 28, 85, and 665 times. The method is demonstrated to converge as indicated by the literature, and its use with an AM-Bench domain is demonstrated.

42 ENGINEERING↗
DOE OSTI2025

Energy absorption of architectured PrintCast interpenetrating composites in tension

Additively manufactured (AM) metal-metal composites consisting of PrintCasted 316 L austenitic stainless-steel lattice structures infiltrated with A356 casting alloy, have recently been developed for use in high energy absorption systems with potential applications ranging from static load bearing to dynamic blast containment structures. This system has a unique mechanical behavior as the volume fraction of lattice increases showing a transition from localized to de-localized failure and dramatic increase in energy absorption capability. In this work, PrintCast A356/316 L composite tensile specimens were produced with lattice volume fractions ranging from 20 % to 50 % to capture the range of this behavior. Finite element simulations support neutron diffraction measurements of stress state. Results illustrate that in tension, the reinforcement material is in tension while the matrix support material is in compression, information offering significant insight into the transition to de-localized failure. Moreover, the simulation results provide further insight into how interfacial bonding (or lack of bonding) affects the energy absorption capabilities of the PrintCast composites.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effects of edge-seal design on the mechanical and thermal performance of vacuum-insulated glazing

Although vacuum-insulated glazing (VIG) has been proposed as a promising solution towards developing energy-efficient buildings, VIGs have not become popular in the market due to several technical challenges including the complexity of the fabrication process. In particular, the edge-seal is a key component that significantly affects the thermal insulation and mechanical performance, and the development of edge-seal with adequate thermal insulation, mechanical strength, and reasonable processing cost is essential to overcome such technical issues in VIG. For this purpose, effects of edge-seal design parameters on the VIG performance should be identified. In this research, we analyzed the edge-seal for thermal transport as well as structural stresses to study the effects, and then identified and evaluated the material mixes for the edge-seal requirements. Here, the finite element simulations showed the significance of VIG corner calculation on overall thermal transmittance and the importance of seal conductivity below 1 W/m.K. The experiments with the flexible seals with different ratios of fine glass powder demonstrated that the measured shear strength values for the seal with less than 30% glass powder were more than 10 times larger than the calculated shear stress values. Based on these simulation and experimental results, a flexible sealant was developed using a proprietary mix of ceramic materials that meets the requirements of the designed VIG edge-seal, including structural as well as thermal stress resistance and a low conductivity. Moreover, the sealant is self-curing under atmospheric conditions, and thus it does not require costly inline process of laser curing or oven baking.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2016

A generalized multi-resolution expansion for uncertainty propagation with application to cardiovascular modeling

Computational models are used in a variety of fields to improve our understanding of complex physical phenomena. Recently, the realism of model predictions has been greatly enhanced by transitioning from deterministic to stochastic frameworks, where the effects of the intrinsic variability in parameters, loads, constitutive properties, model geometry and other quantities can be more naturally included. A general stochastic system may be characterized by a large number of arbitrarily distributed and correlated random inputs, and a limited support response with sharp gradients or event discontinuities. This motivates continued research into novel adaptive algorithms for uncertainty propagation, particularly those handling high dimensional, arbitrarily distributed random inputs and non-smooth stochastic responses. Here, we generalize a previously proposed multi-resolution approach to uncertainty propagation to develop a method that improves computational efficiency, can handle arbitrarily distributed random inputs and non-smooth stochastic responses, and naturally facilitates adaptivity, i.e., the expansion coefficients encode information on solution refinement. Our approach relies on partitioning the stochastic space into elements that are subdivided along a single dimension, or, in other words, progressive refinements exhibiting a binary tree representation. We also show how these binary refinements are particularly effective in avoiding the exponential increase in the multi-resolution basis cardinality and significantly reduce the regression complexity for moderate to high dimensional random inputs. The performance of the approach is demonstrated through previously proposed uncertainty propagation benchmarks and stochastic multi-scale finite element simulations in cardiovascular flow.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Multiscale topology optimization using neural network surrogate models

In this study, we are concerned with optimization of macroscale elastic structures that are designed utilizing spatially varying microscale metamaterials. The macroscale optimization is accomplished using gradient-based nonlinear topological optimization. But instead of using density as the optimization decision variable, the decision variables are the multiple parameters that define the local microscale metamaterial. This is accomplished using single layer feedforward Gaussian basis function networks as a surrogate models of the elastic response of the microscale metamaterial. The surrogate models are trained using highly resolved continuum finite element simulations of the microscale metamaterials and hence are significantly more accurate than analytical models e.g. classical beam theory. Because the derivative of the surrogate model is important for sensitivity analysis of the macroscale topology optimization, a neural network training procedure based on the Sobolev norm is described. Since the SIMP method is not appropriate for spatially varying lattices, an alternative method is developed to enable creation of void regions. Lastly, the efficacy of this approach is demonstrated via several examples in which the optimal graded metamaterial outperforms a traditional solid structure.

42 ENGINEERING↗
DOE OSTI2023

Denoising diffusion algorithm for inverse design of microstructures with fine-tuned nonlinear material properties

Here we introduce a denoising diffusion algorithm to discover microstructures with nonlinear fine-tuned properties. Denoising diffusion probabilistic models are generative models that use diffusion-based dynamics to gradually denoise images and generate realistic synthetic samples. By learning the reverse of a Markov diffusion process, we design an artificial intelligence to efficiently manipulate the topology of microstructures to generate a massive number of prototypes that exhibit constitutive responses sufficiently close to designated nonlinear constitutive behaviors. To identify the subset of microcstructures with sufficiently precise fine-tuned properties, a convolutional neural network surrogate is trained to replace high-fidelity finite element simulations to filter out prototypes outside the admissible range. Results of this study indicate that the denoising diffusion process is capable of creating microstructures of fine-tuned nonlinear material properties within the latent space of the training data. More importantly, this denoising diffusion algorithm can be easily extended to incorporate additional topological and geometric modifications by introducing high-dimensional structures embedded in the latent space. Numerical experiments are conducted on the open-source mechanical MNIST data set (Lejeune, 2020). Consequently, this algorithm is not only capable of performing inverse design of nonlinear effective media, but also learns the nonlinear structure–property map to quantitatively understand the multiscale interplay among the geometry, topology, and their effective macroscopic properties.

42 ENGINEERING↗
DOE OSTI2024

Predicting non-linear stress–strain response of mesostructured cellular materials using supervised autoencoder

Recent breakthroughs in advanced manufacturing capabilities have made it possible to design and print sophisticated topologies of cellular structures using diverse engineering materials such as metals, polymers, and ceramics. In these architectured materials, it is often desirable to tailor the mechanical properties by altering the unit cell topology. This necessitates an in-depth understanding of how the topology of the unit cell structure affects the macroscopic behavior of the material in both the linear and the non-linear regimes encountered under large compression. Here, we have developed a machine learning (ML) approach capable of accelerating the prediction of the stress–strain response of a polymer-based cellular structure under uniaxial confined compression. As part of generating the training data for ML, 60,000 mesostructures were generated using a relatively novel approach based on cellular automata, and their corresponding stress–strain responses were obtained from the finite element simulations. Principal component analysis (PCA) was used to reduce the dimensionality of the stress–strain curves. With only 20 principal components, PCA captured 99.89% of the variance in the stress–strain curves while reducing the dimensionality by 5X. ML using supervised autoencoder was able to successfully speed up the prediction of the non-linear stress–strain response of a unit cell by up to 4600X. The proposed method can serve as an efficient data generation tool and a rapid means for predicting the structure–property relationship through accelerated forward modeling of cellular materials under compaction, in cases where the macroscopic stress–strain response is governed by the unit-cell topology.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Stochastic nonlinear analysis of unidirectional fiber composites using image-based microstructural uncertainty quantification

Here, we present a data-driven nonlinear uncertainty quantification and propagation framework to study the microstructure-induced stochastic performance of unidirectional (UD) carbon fiber reinforced polymer (CFRP) composites. The proposed approach integrates (1) microscopic image characterization, (2) stochastic microstructure reconstruction, and (3) efficient multiscale finite element simulations enabled by self-consistent clustering (SCA) analysis. To model the complex microstructural variability, the proposed UQ methods take the non-Gaussian uncertainty sources into account through a distribution-free sampling approach leveraging nonparametric and asymptotic statistical tools. A hierarchical conditional sampling strategy enables the simultaneous sampling of multiple sources of uncertainties. Our approach provides insights into the impact of microstructural variabilities, which are shown to have an increasing impact on the nonlinear responses of UD CFRP parts under progressive compression loading and ultimately on the failure rate over time. We discover that before CFRP parts start to fail, a characteristic time period emerges with distinctive uncertainty distributions specific to the microstructure variability and the probability of failure. Identifying the failure time period is crucial to the reliability prediction, which is an essential component of CFRP design.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Performance assessment of 3D printed multi-material energy absorber for automotive bumper: pedestrian lower extremity protection

Designing an energy absorber for automotive bumpers involves balancing low-speed and high-speed impacts to ensure safety, reduce repair costs, and meet regulatory standards. Here, this study explores a novel design using multi-material 3D printing and structural optimization to fabricate a lightweight and cost-efficient energy absorber. The design effectively dissipates energy in low-speed collisions and minimizes force transmission in high-speed pedestrain impacts, helping to meet both safety and performance requirements. The energy absorber design combines 20% carbon fiber-reinforced acrylonitrile butadiene styrene (CF-ABS) and thermoplastic polyurethane (TPU) for optimal stiffness and flexibility. It uses 3D-printed lattice structures optimized through finite element simulations to help meet both low-speed and high-speed impact requirements. Full-scale energy absorbers were 3D-printed using optimized CF-ABS/TPU blends and tested under high-speed impact using the Flexible Pedestrian Legform Impactor (Flex-PLI). For fair comparison, a baseline bumper with a traditional triangular lattice structure, also 3D-printed from the same CF-ABS/TPU materials, was similarly tested. Interestingly, both the optimized and baseline 3D-printed energy absorbers showed nearly identical performance, successfully meeting injury limits. Their performances were also benchmarked against an injection-molded energy absorber. While both 3D-printed and injection-molded designs met injury limits, the 3D-printed absorber exhibited a higher tibia bending moment, indicating an opportunity for further optimization. A Techno-Economic Analysis compared the costs of producing energy absorbers using traditional manufacturing and 3D printing. The analysis highlighted that 3D printing offers cost benefits for low to medium production volumes, with the total cost per energy absorber at ∼ $\$$74, compared to traditional methods that become economical beyond 2000 units.

Additive manufacturing↗
DOE OSTI2015

Measurement of fracture toughness by nanoindentation methods: Recent advances and future challenges

In this study, we describe recent advances and developments for the measurement of fracture toughness at small scales by the use of nanoindentation-based methods including techniques based on micro-cantilever beam bending and micro-pillar splitting. A critical comparison of the techniques is made by testing a selected group of bulk and thin film materials. For pillar splitting, cohesive zone finite element simulations are used to validate a simple relationship between the critical load at failure, the pillar radius, and the fracture toughness for a range of material properties and coating/substrate combinations. The minimum pillar diameter required for nucleation and growth of a crack during indentation is also estimated. An analysis of pillar splitting for a film on a dissimilar substrate material shows that the critical load for splitting is relatively insensitive to the substrate compliance for a large range of material properties. Experimental results from a selected group of materials show good agreement between single cantilever and pillar splitting methods, while a discrepancy of ~25% is found between the pillar splitting technique and double-cantilever testing. It is concluded that both the micro-cantilever and pillar splitting techniques are valuable methods for micro-scale assessment of fracture toughness of brittle ceramics, provided the underlying assumptions can be validated. Although the pillar splitting method has some advantages because of the simplicity of sample preparation and testing, it is not applicable to most metals because their higher toughness prevents splitting, and in this case, micro-cantilever bend testing is preferred.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Deformation fields near a steady fatigue crack with anisotropic plasticity

In this work, from finite element simulations based on an irreversible, hysteretic cohesive interface model, a steady fatigue crack can be realized if the crack extension exceeds about twice the plastic zone size, and both the crack increment per loading cycle and the crack bridging zone size are smaller than the plastic zone size. The corresponding deformation fields develop a plastic wake behind the crack tip and a compressive residual stress field ahead of the crack tip. In addition, the Hill’s plasticity model is used to study the role of plastic anisotropy on the retardation of fatigue crack growth and the elastic strain fields. It is found that for Mode-I cyclic loading, an enhanced yield stress in directions that are inclined from the crack plane will lead to slower crack growth rate, but this retardation is insignificant for typical degrees of plastic anisotropy. Furthermore, these results provide key inputs for future comparisons to neutron and synchrotron diffraction measurements that provide full-field lattice strain mapping near fracture and fatigue crack tips, especially in textured materials such as wrought or rolled Mg alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Extremized nonlinear and linearized responses in soft metamaterials enabled by gradient-based design and grayscale digital light processing

In this study, we develop a gradient-based design approach that exploits grayscale digital light processing (DLP) 3D printing for extremizing the nonlinear and linearized response of soft metamaterials — materials that harness engineered geometric instabilities to undergo large and programmable changes in configuration. Grayscale DLP approaches modulate local mechanical properties at the pixel scale by tuning the light intensity within a single grayscale image, unlocking an exceptionally large design space. To effectively navigate this space, we develop smooth mappings between local light intensity values and global quantities of interest that characterize the behavior of soft metamaterials. Enabling these smooth mappings are robust and differentiable nonlinear finite element simulations powered by a trust region solver. A PDE-constrained optimization problem is then solved to invert these mappings and produce light intensity distributions that endow the printed part with varying stiffness and flexibility in distinctive regions. It is shown that optimizing the distribution of soft and stiff phases throughout a metamaterial structure results in markedly different buckling and self-contact configurations to drive extremized nonlinear compression and linearized vibration responses. Optimized light intensity distributions are translated to grayscale images and directly used to print soft metamaterial samples, showing remarkable agreement between the buckling and self-contact response in simulated and measured deformed configurations.

Additive manufacturing↗
DOE OSTI2023

Parametric estimation of Poisson's ratio for thin hinged-hinged plates

Cornu's method is an elegant calculation besieged by an impractical approach to obtain accurate estimates for Poisson's ratio. Conventionally, Cornu's method requires several components, and each component can adversely affect the accuracy of the measurement. Furthermore, Cornu's conventional method requires a long beam because beams with short length-to-width ratios cause the estimate of Poisson's ratio to diverge from the true value of Poisson's ratio. We believe that, with the right modifications, Cornu's method can become an attractive approach to obtaining precise estimates for Poisson's ratio from mode shapes. Here we use finite element simulations to show how to use Cornu's method to estimate Poisson's ratio from a mode shape. Our modified Cornu's method removes knife-edges and loading components for a hinged-hinged plate under steady-state excitation. Given the true value of Poisson's ratio, as a performance specification, we show that simple parametric expressions can fit estimates for Poisson's ratio for different length-to-width ratios of thin hinged-hinged plates. Additionally, we show that estimates for Poisson's ratio from higher modes align with the results from the first mode and explain our expectation for this outcome. Furthermore, our results challenge the idea that anticlastic, monoclastic, and synclastic deformation uniquely correspond to positive, zero, and negative estimates of Poisson's ratio, respectively. With the emergence of materials-by-design, we expect that this parametric technique will be able to assist in experimental qualification of thin beam and plate structures with respect to the desired value of Poisson's ratio.

42 ENGINEERING↗
DOE OSTI2025

Hyperelastic nature of the Hoek–Brown criterion

In this article, we propose a nonlinear elasto-plastic model, for which a specific class of hyperbolic elasticity arises as a straight consequence of the yield criterion invariance on the plasticity level. We superimpose this nonlinear elastic (or hyperelastic) behavior with plasticity obeying the associated flow rule. Interestingly, we find that a linear yield criterion on the thermodynamical force associated with plasticity results in a quadratic yield criterion in the stress space. This suggests a specific hyperelastic connection between Mohr–Coulomb and Hoek–Brown (or alternatively between Drucker–Prager and Pan–Hudson) yield criteria. We compare the elasto-plastic responses of standard tests for the Drucker–Prager yield criterion using either linear or the suggested hyperbolic elasticity. Notably, the nonlinear case stands out due to dilatancy saturation observed during cyclic loading in the triaxial compression test. We conclude this study with structural finite element simulations that clearly demonstrate the numerical applicability of the proposed model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

A visco-plastic framework for interface processes in sedimentary reservoir rocks at HPHT conditions

As energy operations face the challenge of reservoirs at ever-increasing depths, modelling the response of reservoir rocks at high pressure and high temperature (HPHT) conditions is a crucial step for successfully unlocking new resources. At these conditions, rocks can experience thermal and pressure sensitivity, as well as admit internal phase changes at the pore-solid interfaces that determine their macroscopic response to external loading. It is therefore expected that constitutive laws of Geomechanics should be enriched to accommodate such effects. In this work, a multi-physics constitutive theory for sedimentary rocks is proposed within the general framework of viscoplasticity. The viscosity of the material is assumed to be a function of the temperature, pore-pressure and energy required to alter the inter-granular interfaces. This energy is expressed through the chemical potentials of the phases involved during processes like chemical dissolution/precipitation of the interfaces or mechanical debonding. The resulting flow law and corresponding stress equilibrium are coupled to the energy and mass conservation laws, constituting a closed system of equations, which is solved using the Finite Element simulator REDBACK. A series of numerical tests for performance and calibration is then successfully performed against different types of reservoir rocks, confining pressures and temperatures (from room temperature to over 800K). Here, we show that the mechanical response of sedimentary porous rocks at strains usually achieved in laboratory testing can be explained by accounting for the interface processes taking place at the cementitious material bonding the grains, a process that can also determine the brittle-to-ductile transition of sedimentary rocks.

58 GEOSCIENCES↗
DOE OSTI2015

Microstructure-sensitive small fatigue crack growth assessment. Effect of strain ratio multiaxial strain state and geometric discontinuities

Fatigue crack initiation in the high cycle fatigue regime is strongly influenced by microstructural features. Research efforts have usually focused on predicting fatigue resistance against crack incubation without considering the early fatigue crack growth after encountering the first grain boundary. However, a significant fraction of the variability of the total fatigue life can be attributed to growth of small cracks as they encounter the first few grain boundaries, rather than crack formation within the first grain. Our paper builds on the framework previously developed by the authors to assess microstructure-sensitive small fatigue crack formation and early growth under complex loading conditions. Moreover, the scheme employs finite element simulations that explicitly render grains and crystallographic directions along with simulation of microstructurally small fatigue crack growth from grain to grain. The methodology employs a crystal plasticity algorithm in ABAQUS that was previously calibrated to study fatigue crack initiation in RR1000 Ni-base superalloy. Our work present simulations with non-zero applied mean strains and geometric discontinuities that were not previously considered for calibration. Results exhibit trends similar to those found in experiments for multiple metallic materials, conveying a consistent physical description of fatigue damage phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Eulerian finite element implementations of a dislocation density-based continuum model

In Eulerian finite element simulations, the mesh moves relative to the material. After every change of position between the mesh and the material, the state variables are interpolated to the new mesh position, which is referred to as advection. Large strain crystal plasticity models are based on the multiplicative decomposition of the total deformation gradient. Further, the stress is evaluated as a function of the thermoelastic strain, temperature, and other state variables. Advection of tensor quantities, such as the strain, is coupled with possibly significant advection errors. In an effort to reduce the advection errors, we develop two rate forms of an established dislocation density-based continuum model. To that end, we replace the multiplicative decomposition of the deformation gradient with the additive decomposition of the velocity gradient, and define the stress rate instead of the total stress. The Eulerian implementation is compared with Lagrangian calculations, and two numerical examples with severe deformation levels are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Effects of geometric factors and shear band patterns on notch sensitivity in bulk metallic glasses

Our recent experiments in notched bulk metallic glasses have found reduced, or insensitive, or improved strengths, while in many of these cases the ductile strain prior to final failure is enhanced. First, although the inverse notch effect is explained by a shift from shear localization to cavitation failure, it is suggested in this work that the synergistic effect between cohesive fracture at the notched area and shear bands emanating from the notch roots may extend the parametric space for the notch insensitive behavior. Second, the dependence of shear band patterns on notch geometric factors is determined by the Rudnicki-Rice theory and the free-volume-based finite element simulations. Our results suggest conditions for shear band multiplication to take place and for the shear-localization-induced failure to be delayed.

36 MATERIALS SCIENCE↗