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 109 records · Page 6

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↗

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.↗

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.↗

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.↗

Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design

Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.

Fahrenthold, Eric P.↗

Surface-Correlated Nanophase Iron Metal in Lunar Soils: Petrography and Space Weathering Effects

Space weathering is a term used to include all of the processes that act on material exposed at the surface of a planetary or small body. In the case of the Moon, it includes a variety of processes that formed the lunar regolith, caused the maturation of lunar soils, and formed patina on rock surfaces. The processes include micrometeorite impact and reworking, implantation of solar wind and flare particles, radiation damage and chemical effects from solar particles and cosmic rays, interactions with the lunar atmosphere, and sputtering erosion and deposition. Space weathering effects collectively result in a reddened continuum slope, lowered albedo, and attenuated absorption features in reflectance spectra of lunar soils as compared to finely comminuted rocks from the same Apollo sites. Understanding these effects is critical in order to fully integrate the lunar sample collection with remotely sensed data from recent robotic missions (e.g., Lunar Prospector, Clementine, Galileo). Our objective is to determine the origin of space weathering effects in lunar soils through combined electron microscopy and microspectrophotometry techniques applied to individual soil particles from <20 pm size factions (dry-sieved) of mature lunar soils. It has been demonstrated that it is the finest size fraction (<25 pm) of lunar soils that dominates the optical properties of the bulk soils.

Keller, Lindsay P.↗

Microstructural Studies of Space Weathering Effects in Lunar Materials

Space weathering is a term used to include all of the processes which act on material exposed at the surface of a planetary or small body. In the case of the moon, it includes a variety of processes which have formed the lunar regolith, caused the maturation of lunar soils, and formed patina on rock surfaces. The processes include micrometeorite impact and reworking, implantation of solar wind and flare particles, radiation damage and chemical effects from solar particles and cosmic rays, interactions with the lunar atmosphere, and sputtering erosion and deposition. Space weathering effects collectively result in a reddened continuum slope, lowered albedo, and attenuated absorption features in reflectance spectra of lunar soils as compared to finely comminuted rocks from the same Apollo sites. However, the regolith processes that cause these effects are not well known, nor is the petrographic setting of the products of these processes fully understood. An interesting confluence of events occurred in the early 1990s. First, came the discovery of vapor deposited coatings on lunar regolith grains by Keller and McKay, who showed that amorphous coatings from 50-100 nm thick containing fine-grained Fe metal (1-10 nm in diameter) were common in the fine size fraction of several mature lunar soils. The other discovery was the recognition that the optical properties of lunar soils were dominated by fine grain sized material (less than 45 micrometer fraction) by Pieters and coworkers. These discoveries led to coordinated studies that looked at the mineralogy, chemistry, and optical properties of lunar soils as function of composition, maturity, and grain size fraction. One of the major revelations from these studies was the recognition that much of the nanophase Fe metal is surface-correlated especially in the finest size fractions, and that it was this nanophase Fe that dominated the optical properties of the soil.

Keller, L. P.↗

Computational Modeling and Experimental Characterization of Martensitic Transformations in Nicoal for Self-Sensing Materials

Fundamental changes to aero-vehicle management require the utilization of automated health monitoring of vehicle structural components. A novel method is the use of self-sensing materials, which contain embedded sensory particles (SP). SPs are micron-sized pieces of shape-memory alloy that undergo transformation when the local strain reaches a prescribed threshold. The transformation is a result of a spontaneous rearrangement of the atoms in the crystal lattice under intensified stress near damaged locations, generating acoustic waves of a specific spectrum that can be detected by a suitably placed sensor. The sensitivity of the method depends on the strength of the emitted signal and its propagation through the material. To study the transition behavior of the sensory particle inside a metal matrix under load, a simulation approach based on a coupled atomistic-continuum model is used. The simulation results indicate a strong dependence of the particle's pseudoelastic response on its crystallographic orientation with respect to the loading direction and suggest possible ways of optimizing particle sensitivity. The technology of embedded sensory particles will serve as the key element in an autonomous structural health monitoring system that will constantly monitor for damage initiation in service, which will enable quick detection of unforeseen damage initiation in real-time and during onground inspections.

Wallace, T. A.↗

Electrical characterization of plasma-grown oxides on gallium arsenide

Plasma-grown GaAs oxides and their interfaces have been characterized by measuring the electrical properties of metal-oxide-semiconductor capacitors and of Schottky junctions. The current transport mechanism in the oxide at high electrical field was found to be Frankel-Poole emission, with an electron trap center at 0.47 eV below the conduction band of the oxide. The interface-state density, evaluated from capacitance and conductance measurements, exhibits a U-shaped interface-state continuum extending over the entire band gap. Two discrete deep states with high concentration are superimposed on this continuum at 0.40 and 0.70 eV below the conduction band. The results obtained from measurements on Schottky junctions have excluded the possibility that these two deep states originate from plasma damage. Possible origins of these states are discussed in this paper.

Hshieh, F. I.↗