Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum damage”

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 163 records · Page 9

The Development of Computational Techniques for the Nonlinear Analysis of Composite Structures at High Temperature

The objective of this research project was to develop robust and efficient computational tools for deformation and life analysis techniques. This included: development of viscoplastic continuum models, associated numerical integration techniques, suitable micromechanics models/methods, and coupled deformation and damage algorithms. The work was involved with models/methods applicable to composite structures on both microscale and macroscale levels. The microscale analysis utilized Aboudi's generalized method of cells (GMS) micromechanical model as well as finite element models of similar composite architectures, i.e., square and hexagonal fiber pack models.

Keith, Theo G., Jr.↗

Modeling Heatshield Erosion due to Dust Particle Impacts for Martian Entries

Heatshield design for spacecraft entering the atmosphere of Mars may be affected by the presence of atmospheric dust. Particle impacts with sufficient kinetic energy can cause spallation damage to the heatshield that must be estimated. The dust environment in terms of particle size distribution and number density can be inferred from ground-based or atmospheric observations at Mars. Using a Lagrangian approach, the particle trajectories through the shock layer can be computed using a set of coupled ordinary differential equations. The dust particles are small enough that non-continuum effects must be accounted for when computing the drag coefficient and heat transfer to the particle surface. Surface damage correlations for impact crater diameter and penetration depth are presented for fused-silica, AVCOAT, Shuttle tiles, cork, and Norcoat Liège. The cork and Norcoat Liège correlations are new and were developed in this study. The modeling equations presented in this paper are applied to compute the heatshield erosion due to dust particle impacts on the ExoMars Schiaparelli entry capsule during dust storm conditions.

Palmer, Grant↗

Modeling Heatshield Erosion Due to Dust Particle Impacts for Martian Entries

Heatshield design for spacecraft entering the atmosphere of Mars may be affected by the presence of atmospheric dust. Particle impacts with sufficient kinetic energy can cause spallation damage to the heatshield that must be estimated. The dust environment in terms of particle size distribution and number density can be inferred from ground-based or atmospheric observations at Mars. Using a Lagrangian approach, the particle trajectories through the shock layer can be computed using a set of coupled ordinary differential equations. The dust particles are small enough that non-continuum effects must be accounted for when computing the drag coefficient and heat transfer to the particle surface. Surface damage correlations for impact crater diameter and penetration depth are presented for fused-silica, AVCOAT, Shuttle tiles, cork, and Norcoat Liège. The cork and Norcoat Liège correlations are new and were developed in this study. The modeling equations presented in this paper are applied to compute the heatshield erosion due to dust particle impacts on the ExoMars Schiaparelli entry capsule during dust storm conditions.

Palmer, Grant↗

Micromechanics Fatigue Damage Analysis Modeling for Fabric Reinforced Ceramic Matrix Composites

A micromechanics analysis modeling method was developed to analyze the damage progression and fatigue failure of fabric reinforced composite structures, especially for the brittle ceramic matrix material composites. A repeating unit cell concept of fabric reinforced composites was used to represent the global composite structure. The thermal and mechanical properties of the repeating unit cell were considered as the same as those of the global composite structure. The three-phase micromechanics, the shear-lag, and the continuum fracture mechanics models were integrated with a statistical model in the repeating unit cell to predict the progressive damages and fatigue life of the composite structures. The global structure failure was defined as the loss of loading capability of the repeating unit cell, which depends on the stiffness reduction due to material slice failures and nonlinear material properties in the repeating unit cell. The present methodology is demonstrated with the analysis results evaluated through the experimental test performed with carbon fiber reinforced silicon carbide matrix plain weave composite specimens.

Min, J. B.↗

Next Generation Durability and Damage Tolerance to Support Certification of Flight Hardware

Durability and Damage Tolerance (D&DT), as currently applied to flight hardware throughout the Agency, is based on continuum and similitude assumptions and does not consider local material properties, environments and responses. These limitations impact our ability to support certification of new component designs, material systems and manufacturing approaches. The primary concern in this Engineering R&A plan is related to the imminent proliferation of parts produced by additive manufacturing (AM). These parts are being advocated by NASA’s suppliers for use on a myriad of flight hardware because of their design flexibility and cost advantages. AM provides opportunities to reduce part counts through complex geometry and reduce manufacturing costs of low volume parts. However, AM materials have some notable metallurgical and microstructural differences compared to traditionally fabricated materials. One of the challenges for the acceptance of AM is the greater tendency for a deleterious defect state, most commonly in the form of porosity, to exist in AM parts. Though this defect state is typically reduced through a hot isostatic pressing (HIP) processing step, the structural performance risks associated with the remaining defects and HIP-healed features is not known. Potential fracture control issues that must be resolved stem from the real possibility that an unhealed defect or a closed defect with less than pristine strength remains at a fracture critical location after HIP. As a result, NASA and the entire aerospace community have been confronted with the need to develop a robust and relevant certification methodology to enable safe implementation of these components. The present work is an important step toward positioning NASA to credibly respond to vendors’ push to implement this new AM materials technology by improving our understanding of AM processing and performance and transitioning that research-based understanding to next-generation engineering capabilities.

Edward H Glaessgen↗

Finite Element Modeling of Plastic Deformation During Spin Forming of Aluminum 6061-O

Spin and flow forming are metal deformation techniques in which a disk or tube of material is radially thinned and axially lengthened over a rotating mandrel. A finite-element continuum model (FEM) based on commercial software (DEFORM®) was developed to determine the stress/strain distribution and damage accumulation in a flow-formed, near-net-shape part. Computational validation and experimental verification of the model was leveraged to assess the effectiveness of computer simulations. Quantitative comparisons between experimental and computational results show promise for applying a FEM to the flow forming process. Experimental failure locations correlated well with computed damage gradients, and there was good agreement between the measured and predicted roller forces during forming.

Elizabeth Urig↗

Novel Technology of Non-Contact Real-Time Radiation Damage Sensor for High Power Targets

High-power proton beams planned for forthcoming long-baseline neutrino experiments will subject solid targets to unprecedented radiation damage, threatening reliability and increasing costs. To address this challenge, we are building a real-time, non-contact sensor that tracks damage by measuring broadband laser reflectivity changes from the target surface. A low-power super-continuum Class 3B laser illuminates the sample inside a vacuum test chamber while a high-resolution fiber-coupled spectrometer captures S- and P-polarized light; spectral shifts reveal defect-driven variations in optical constants. My internship goal is to design, build, and commission this prototype by implementing a laser-safety interlock and light shield, integrating remote-operation cameras, and fabricating modular 3-D-printed mounts that enable tool-free swaps without disturbing alignment. The laser, spectrometer, and vacuum chamber have been delivered; interlock hardware, cameras, and mounts are in final assembly, and leak testing of the chamber is underway. Upcoming work will focus on initial calibration of the sensor and executing first beam-irradiation studies.

Pumarino, Rafael↗

Enhancing Lithium-Ion Battery Aging Simulations: Coupling a High-Resolution, 3D, Grain-Scale Electromechanical Model to a Single-Particle Model

One of the main goals in modeling lithium-ion batteries is to improve/predict longevity and resilience of new chemistries. Unfortunately, this requires simulation of thousands of charge/discharge cycles, which can be rather time consuming depending on the fidelity of the simulation. The purpose of this talk is to discuss a new modeling framework that couples a high-resolution, continuous damage model (CDM) to a single particle model (SPM) resulting in a good combination of speed and fidelity. In previous work, a 3D, continuum-level chemo-mechanical model was developed to investigate cracking within a single cathode particle comprised of hundreds to thousands of randomly oriented grains. The CDM predicted that particle fracture is primarily due to non-ideal grain interactions with slight dependence on high-rate charge demands. Essentially, when neighboring grains were misaligned, they expanded different rates relative to one another leading to high stresses and ultimately the formation of intraparticle cracks. The model predicted that small particles with large grains develop significantly less damage than larger particles with small grains. Finally, the model predicted most of the chemo-mechanical damage accumulates in the first charge after formation. This chemo-mechanical "damage saturation" effect indicated that initial particle fracture occurs within the first few cycles, while long-term cathode degradation is not solely chemo-mechanically induced. This led to a need for simulating fatigue-like mechanism that degrade the battery over longer time scales. In order to reach the time scales, need to resolve fatigue-like degradation, the CDM needs to be complemented by a faster model. Therefore, recent efforts have been focused on using results from the CDM to inform parameters within the SPM. These parameters are homogenization factors that are associated with diffusion, particle radius, and/or exchange current density. By coupling the CDM to the SPM, the aging simulation is broken up into two domains: short-term and long-term degradation. The short-term degradation occurs over a single cycle and is handled by the CDM due of its high fidelity, but relatively expensive computational cost. Such mechanism include break-in crack caused by mismatches in grain orientation. The long-term degradation occurs over tens of cycles and is handled by the SPM due it's computational efficiency. Most of the fatigue-like mechanism fall into this category. The eventual goal of this modeling framework is to upscale to a psudo-2D model allowing for full electrode simulation, which are informed by high-fidelity grain-scale simulations.

electrochemistry↗

Deformation mechanisms in negative Poisson's ratio materials - Structural aspects

Poisson's ratio in materials is governed by the following aspects of the microstructure: the presence of rotational degrees of freedom, non-affine deformation kinematics, or anisotropic structure. Several structural models are examined. The non-affine kinematics are seen to be essential for the production of negative Poisson's ratios for isotropic materials containing central force linkages of positive stiffness. Non-central forces combined with pre-load can also give rise to a negative Poisson's ratio in isotropic materials. A chiral microstructure with non-central force interaction or non-affine deformation can also exhibit a negative Poisson's ratio. Toughness and damage resistance in these materials may be affected by the Poisson's ratio itself, as well as by generalized continuum aspects associated with the microstructure.

Lakes, R.↗

Study of Plume Impingement Effects in the Lunar Lander Environment

Plume impingement effects from the descent and ascent engine firings of the Lunar Lander were analyzed in support of the Lunar Architecture Team under the Constellation Program. The descent stage analysis was performed to obtain shear and pressure forces on the lunar surface as well as velocity and density profiles in the flow field in an effort to understand lunar soil erosion and ejected soil impact damage which was analyzed as part of a separate study. A CFD/DSMC decoupled methodology was used with the Bird continuum breakdown parameter to distinguish the continuum flow from the rarefied flow. The ascent stage analysis was performed to ascertain the forces and moments acting on the Lunar Lander Ascent Module due to the firing of the main engine on take-off. The Reacting and Multiphase Program (RAMP) method of characteristics (MOC) code was used to model the continuum region of the nozzle plume, and the Direct Simulation Monte Carlo (DSMC) Analysis Code (DAC) was used to model the impingement results in the rarefied region. The ascent module (AM) was analyzed for various pitch and yaw rotations and for various heights in relation to the descent module (DM). For the ascent stage analysis, the plume inflow boundary was located near the nozzle exit plane in a region where the flow number density was large enough to make the DSMC solution computationally expensive. Therefore, a scaling coefficient was used to make the DSMC solution more computationally manageable. An analysis of the effectiveness of this scaling technique was performed by investigating various scaling parameters for a single height and rotation of the AM. Because the inflow boundary was near the nozzle exit plane, another analysis was performed investigating three different inflow contours to determine the effects of the flow expansion around the nozzle lip on the final plume impingement results.

Marichalar, Jeremiah↗

Analysis Methods for Progressive Damage of Composite Structures

This document provides an overview of recent accomplishments and lessons learned in the development of general progressive damage analysis methods for predicting the residual strength and life of composite structures. These developments are described within their State-of-the-Art (SoA) context and the associated technology barriers. The emphasis of the authors is on developing these analysis tools for application at the structural level. Hence, modeling of damage progression is undertaken at the mesoscale, where the plies of a laminate are represented as a homogenous orthotropic continuum. The aim of the present effort is establish the ranges of validity of available models, to identify technology barriers, and to establish the foundations of the future investigation efforts. Such are the necessary steps towards accurate and robust simulations that can replace some of the expensive and time-consuming "building block" tests that are currently required for the design and certification of aerospace structures.

Rose, Cheryl A.↗

Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures

Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.

97 MATHEMATICS AND COMPUTING↗

The dynamic tensile strength of ice and ice-silicate mixtures

The dynamic tensile strength of icy media is measured at strain rates on the order of 10,000/sec to aid in the understanding of impact and cratering phenomena. Compressed samples consisting of ice and ice-silicate mixtures with 5 and 30 wt % sand were impacted at temperatures between 230 and 250 K by projectile plexiglas plates imparting the required strain rates in less than 0.75 microsec. Taking the tensile stress corresponding to the transition from intact to spalled or fragmented samples as the dynamic tensile strength, strengths of 17, 20 and 22 MPa were obtained for the pure ice, 5 wt % sand, and 30 wt % sand specimens, respectively. The values lie considerably above those observed in static testing. A continuum fracturing model is used to obtain relations between tensile strength and stress rate as well as to derive stress and damage histories during tensile loading and the size distribution of icy fragments as a function of strain rate.

Lange, M. A.↗

Computational Implementation of a Thermodynamically Based Work Potential Model For Progressive Microdamage and Transverse Cracking in Fiber-Reinforced Laminates

A continuum-level, dual internal state variable, thermodynamically based, work potential model, Schapery Theory, is used capture the effects of two matrix damage mechanisms in a fiber-reinforced laminated composite: microdamage and transverse cracking. Matrix microdamage accrues primarily in the form of shear microcracks between the fibers of the composite. Whereas, larger transverse matrix cracks typically span the thickness of a lamina and run parallel to the fibers. Schapery Theory uses the energy potential required to advance structural changes, associated with the damage mechanisms, to govern damage growth through a set of internal state variables. These state variables are used to quantify the stiffness degradation resulting from damage growth. The transverse and shear stiffness of the lamina are related to the internal state variables through a set of measurable damage functions. Additionally, the damage variables for a given strain state can be calculated from a set of evolution equations. These evolution equations and damage functions are implemented into the finite element method and used to govern the constitutive response of the material points in the model. Additionally, an axial failure criterion is included in the model. The response of a center-notched, buffer strip-stiffened panel subjected to uniaxial tension is investigated and results are compared to experiment.

Pineda, Evan J.↗

Machine learning-enabled multiscale modeling of mechanical deformation of aluminum and Al-SiC nanocomposites

A machine learning-enabled multiscale framework is developed for modeling the mechanical response of both pure metal and nanoparticle-reinforced metal matrix nanocomposites (MMNCs). Using aluminum–silicon carbide (Al-SiC) as an example MMNC, atomistic simulations reveal three distinct deformation mechanisms (i.e., defect-free, dislocation-based, and interface separation) governed by the interfaces between the Al matrix and SiC nanoparticles. As compared with single crystal Al, the lattice undergoes a more abrupt failure once the dislocation network becomes extensive and void nucleation initiates, whereas in Al-SiC, nanoparticle interfaces enable a more gradual progression of damage. These mechanisms are captured through a combined classification-regression neural network surrogate model that bridges atomic-scale insights with continuum-scale finite element analysis. Machine learning-enabled multiscale modeling of pure Al accurately predicted strain localization and confirmed by in-situ scanning electron microscopic tensile testing on perforated Al specimens. This study underscores the promise of integrating physics-informed machine learning with hierarchical modeling to capture the interface dominated phenomena and guide the design of advanced MMNCs.

Al-SiC↗

Accelerating high-strain continuum-scale brittle fracture simulations with machine learning

Failure in brittle materials under dynamic loading conditions is a result of the propagation and coalescence of microcracks. Simulating this discrete crack evolution at the continuum level is computationally expensive or, in some cases, intractable, resulting in the need to make broad assumptions or neglect key physics. In this work, we have developed an approach using machine learning that overcomes the current inability to represent meso-scale physics at the macro-scale. Our approach leverages damage and stress data from a computationally expensive high-fidelity model that explicitly resolves microcrack behavior to build an inexpensive machine learning emulator. Once trained, the machine learning emulator is used to predict the evolution of crack length statistics, which then informs a continuum-scale constitutive model. This results in a significant speed-up of the workflow by four orders of magnitude. Both the machine learning emulator and the continuum-scale model are validated against the high-fidelity model and experimental data, respectively, showing excellent agreement. There are two key findings. The first is that we can reduce the dimensionality of the problem, establishing that the machine learning emulator only needs the length of the longest crack and one of the maximum stress components to capture the necessary physics. Another compelling finding is that the emulator can be trained in one experimental setting and transferred successfully to predict behavior in a different setting.

36 MATERIALS SCIENCE↗

Plasticity and Damage Modeling of Stress Asymmetry and Dynamic Behavior of AFS Additive Manufactured Aluminum Alloy 2219

The Solid State Additive Manufacturing (AM) process referred as MELD that fabricated the samples in this study, provides a new path for repairing, coating, joining and additive manufacturing metals and metal matrix composites. This research will be the first application of a physics-based microstructure dependent internal state variable (ISV) plasticity and damage material model to capture the mechanical response of an AM Aluminum Alloy (AA) 2219 via the MELD process. In this research, a microstructure-based internal state variable (ISV) plasticity-damage model was used to capture the mechanical behavior of AFS 2219 aluminum alloy. Aeroprobe Corporation, creator and patent holder for the MELD process, fabricated the material by pushing a solid filler rod of AA2219-T861 material through a hollow rotating tool onto an AA2219 T851 plate substrate. As feedstock, solid or powder precursor metals are pushed through a nonconsumable rotating cylindrical tool. Herein, added layers are deposited and metallurgically bonded to substrate material or previously deposited layers by the heat generated from the rotating tool through plastic deformation of the filler material. Once a layer has been added, the tool height increases, and starts the deposition of the next layer. This process results in beneficial properties such as grain refinement, homogenization and reduced porosity (fully dense). This process will experience temperatures similar to those in the weld nugget zone (WNZ) in friction stir welding (FSW), ranging from 0.6-0.9 Tm, with Tm being the melting point of the material. MELD is highly scalable with AA deposition rates reaching over 1000 cm3/hr, which allows for MELD being used for repairs, coatings, and building components. A motivating factor driving the research for physics-based history dependent material modeling of MELD components is the ability to accurately capture the stress-state and strain rate dependence in the material caused by variations in material microstructure from the MELD processing of new or repaired components. The ISV model incorporates microstructural content and is consistent with continuum level kinematics, kinetics, and thermodynamics. These features allow the ISV model to capture large deformations at the structural scale using the kinematic and isotropic hardening, while microscale damage is obtained from the microstructural features. The benefits of the ISV model arise from the inclusion of structure-property relationships identified from microstructural characterization and experimentation. The Bauschinger effect (BE) is an important concept, vital in the accurate prediction of cyclic stress-strain response of ductile materials such as metals. The ISV model has been successfully used to capture the behavior and damage, and the BE of different aluminum alloys and steels. The ISV model uses kinematic and isotropic hardening to help capture deformations of the material at the macro scale. To understand this hardening relationship, calculating the kinematic and isotropic hardening relationship in the material is warranted for a high-fidelity model. Electron Backscattered Diffraction (EBSD) was used to characterize the as-fabricated microstructure, where a fully-dense equiaxed grain morphology with average grain size of 2.5 μm was observed. Microhardness mapping of the as-built structures, monotonic tension and compression experiments at both quasi-static (0.001/s) strain rates, tension-followed-by-compression and compression-followed-by-tension experiments were performed to obtain the set of plasticity and damage constants necessary to capture strain rate and stress state behavior of this additive material. To calibrate the plasticity-damage model, a single set of constants were determined to capture the different stress states the MELD AA2219. One set of the constants was determined from experimental true stress-strain curves for the tension and compression data. Additionally, microstructural information and data from the open literature were used as the other model constants. This research is a first of its kind for AFS AA2219, includes correlating the ISV model to the monotonic experimental results that capture the isotropic and kinematic plasticity mechanical response.

Rivera, O. G.↗

A methodology to predict damage initiation, damage growth and residual strength in titanium matrix composites

In this research, a methodology to predict damage initiation, damage growth, fatigue life, and residual strength in titanium matrix composites (TMC) is outlined. Emphasis was placed on micromechanics-based engineering approaches. Damage initiation was predicted using a local effective strain approach. A finite element analysis verified the prevailing assumptions made in the formulation of this model. Damage growth, namely, fiber-bridged matrix crack growth, was evaluated using a fiber bridging (FB) model which accounts for thermal residual stresses. This model combines continuum fracture mechanics and micromechanics analyses yielding stress-intensity factor solutions for fiber-bridged matrix cracks. It is assumed in the FB model that fibers in the wake of the matrix crack are idealized as a closure pressure, and an unknown constant frictional shear stress is assumed to act along the debond length of the bridging fibers. This frictional shear stress was used as a curve fitting parameter to the available experimental data. Fatigue life and post-fatigue residual strength were predicted based on the axial stress in the first intact 0 degree fiber calculated using the FB model and a three-dimensional finite element analysis.

Bakuckas, J. G., Jr.↗