Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture 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 1,081 records · Page 60

Development of GENOA Progressive Failure Parallel Processing Software Systems

A capability consisting of software development and experimental techniques has been developed and is described. The capability is integrated into GENOA-PFA to model polymer matrix composite (PMC) structures. The capability considers the physics and mechanics of composite materials and structure by integration of a hierarchical multilevel macro-scale (lamina, laminate, and structure) and micro scale (fiber, matrix, and interface) simulation analyses. The modeling involves (1) ply layering methodology utilizing FEM elements with through-the-thickness representation, (2) simulation of effects of material defects and conditions (e.g., voids, fiber waviness, and residual stress) on global static and cyclic fatigue strengths, (3) including material nonlinearities (by updating properties periodically) and geometrical nonlinearities (by Lagrangian updating), (4) simulating crack initiation. and growth to failure under static, cyclic, creep, and impact loads. (5) progressive fracture analysis to determine durability and damage tolerance. (6) identifying the percent contribution of various possible composite failure modes involved in critical damage events. and (7) determining sensitivities of failure modes to design parameters (e.g., fiber volume fraction, ply thickness, fiber orientation. and adhesive-bond thickness). GENOA-PFA progressive failure analysis is now ready for use to investigate the effects on structural responses to PMC material degradation from damage induced by static, cyclic (fatigue). creep, and impact loading in 2D/3D PMC structures subjected to hygrothermal environments. Its use will significantly facilitate targeting design parameter changes that will be most effective in reducing the probability of a given failure mode occurring.

Abdi, Frank↗

Damage Accumulation and Failure of Plasma-Sprayed Thermal Barrier Coatings under Thermal Gradient Cyclic Conditions

Thermal barrier coatings will be more aggressively designed to protect gas turbine engine hot-section components in order to meet future engine higher fuel efficiency and lower emission goals. A fundamental understanding of the sintering and thermal cycling induced delamination of thermal barrier coating systems under engine-like heat flux conditions will potentially help to improve the coating temperature capability. In this study, a test approach is established to emphasize the real-time monitoring and assessment of the coating thermal conductivity, which can initially increase under the steady-state high temperature thermal gradient test due to coating sintering, and later decrease under the thermal gradient cyclic test due to coating cracking and delamination. Thermal conductivity prediction models have been established for a ZrO2-(7- 8wt%)Y2O3 model coating system in terms of heat flux, time, and testing temperatures. The coating delamination accumulation is then assessed based on the observed thermal conductivity response under the combined steady-state and cyclic thermal gradient tests. The coating thermal gradient cycling associated delaminations and failure mechanisms under simulated engine heat-flux conditions will be discussed in conjunction with the coating sintering and fracture testing results.

Zhu, Dongming↗

Benchmarking Mixed Mode Failure in Progressive Damage and Failure Analysis Methods

The verification and validation of progressive-damage-analysis finite element methods are difficult but critical tasks to undertake during their development. Verification exercises assess whether a predictive analysis tool produces results that are consistent with the fundamental concepts and assumptions of the tool under evaluation. Ideally, closed-form analytical solutions can be derived for which method verification results can be compared. Problems selected for computational tool verification are often simple and isolate individual features of the tool. In the case of progressive damage finite element methods, verifications should be performed to evaluate the ability of the model to predict the initiation of damage and its growth through the finite element mesh under a variety of conditions. Mabson et al. proposed a test case of a unidirectional, fiber-reinforced plate with a center crack subjected to tensile loads to evaluate matrix crack propagation predictions. The problem was modeled using the Abaqus Hashin continuum damage mechanics (CDM) model for fiber-reinforced composites. Different combinations of matrix strengths and element sizes were used in the simulations, and the results were compared to a closed-form solution based on linear elastic fracture mechanics (LEFM). It was determined that the Abaqus CDM model could predict the LEFM solution of Mode I cracks only when the finite element mesh density met specific requirements based on the material properties. This paper presents closed-form LEFM solutions for a center notch mixed mode (CNMM) verification problem. Parametric finite element analyses were developed using progressive damage analysis methods of both the Discrete Damage Mechanics (DDM) and CDM classes. The progressive damage analysis methods applied in the analyses of the CNMM problem include CompDam and the Floating Node Method. Analyses were conducted with various mode mixities and element sizes to verify that the damage models were working as intended and to identify any limits of applicability.

Leone, Frank↗

Hardness-based plasticity and fracture model for quench-hardenable boron steel (22MnB5)

A comprehensive strain hardening and fracture characterization of different grades of boron steel blanks has been performed, providing the foundation for the implementation into the modular material model (MMM) framework developed by Volkswagen Group Research for an explicit crash code. Due to the introduction of hardness-based interpolation rules for the characterized main grades, the hardening and fracture behavior is solely described by the underlying Vickers hardness. In other words, knowledge of the hardness distribution within a hot-formed component is enough to set up the newly developed computational model. The hardness distribution can be easily introduced via an experimentally measured hardness curve or via hardness mapping from a corresponding hot-forming simulation. For industrial application using rather coarse and computationally inexpensive shell element meshes, the user material model has been extended by a necking/post-necking model with reduced mesh-dependency as an additional failure mode. The present paper mainly addresses the necking/post-necking model.

36 MATERIALS SCIENCE↗

The Effect of Aging on the Delamination Fracture Energy of Glass Solar Reflectors

Thin metal and metal oxide films deposited on glass or polymer substrates have been widely applied to manufacture front surfaced solar reflectors such as those used in the concentrated solar power (CSP) systems due to the facts that they are highly transparent, resistant to ultraviolet (UV) light and abrasion, and chemically inert. However, the reliable operation over long periods of outdoor exposure of such solar reflectors depends on the adequate adhesion between each layer to prevent delamination and subsequent failure. We report results from measurements of the delamination fracture energy (DFE) of the weakest interface in both as-deposited and aged reflectors, where the reflector construction consists of a glass/copper/silver/oxide stack. In this study, DFE was measured using a width-tapered beam method and the test results showed that the weakest interface of this solar reflector is the copper-glass interface, with a delamination fracture energy of approximately 4.4 J/m 2 . In order to verify the test results, finite element (FE) simulations were run using the commercial software package ABAQUS 6.14 based on the cohesive zone model (CZM). In the simulations, the DFE derived from the delamination measurement was treated as an input and the delamination forces from the FE simulations were compared to the experimental results. The good agreement between the tested and simulated delamination force indicated that this width-tapered beam method is accurate even for very small DFE, making it appropriate to test the DFE of front surfaced glass solar reflectors. Additionally, in order to study the aging effect on the DFE, accelerated aging tests were conducted using samples with different protective coatings, one group of samples with alumina ( Al 2 O 3 ) protective coating and the other group with titania ( TiO 2 ) protective coating. The DFE after 50 and 100 hours of accelerated aging were measured. The test results showed that the copper-glass interface would still be the weakest interface for all samples after aging, and the DFE for samples with Al 2 O 3 protective coating was unchanged whereas the copper-glass interface was found to undergo statistically significant change after both 50 and 100 hours of aging for samples with TiO 2 protective coating. We found that DFE at the copper/glass interface decreased after 50 hours of aging then increased after 100 hours of aging. The origin of this non-monotonic response of TiO 2 coated reflectors under accelerated aging conditions is unknown.

36 MATERIALS SCIENCE↗

Mechanism of Decrease in Impact Toughness in a Low-Carbon MnCrMoNiCu Plate Steel with Increasing Austenitizing Temperature

In order to reveal how microscopic factors affect the toughness and the occurrence of cleavage fracture of a low-carbon MnCrMoNiCu alloyed steel, a series of thermal treatments was performed on the steel employing a thermomechanical simulator. These involved reheating samples at different temperatures (950-1250 °C), producing different prior austenite sizes, followed by a continuous cooling transformation process. The Charpy V-notch toughness was determined, and the effect of austenite grain size on the ductile-to-brittle transition temperatures of the steel was investigated. The microstructural evolution on the austenite sizes was studied, fracture features were characterized, the critical event for cleavage fracture was identified, and the local cleavage fracture stress σ{sub f} was calculated. The impact toughness decreased as the austenitizing temperature increased. A quantitative relationship between σ{sub f} and the size of the initial cleavage fracture facet (microcrack nucleus) a{sub f} in the lathy martensite + bainite microstructure has been developed.

36 MATERIALS SCIENCE↗

Mathematical Foundations for Nonlocal Interface Problems: Multiscale Simulations of Heterogeneous Materials (Final LDRD Report)

Nonlocal models provide a much-needed predictive capability for important Sandia mission applications, ranging from fracture mechanics for nuclear components to subsurface flow for nuclear waste disposal, where traditional partial differential equations (PDEs) models fail to capture effects due to long-range forces at the microscale and mesoscale. However, utilization of this capability is seriously compromised by the lack of a rigorous nonlocal interface theory, required for both application and efficient solution of nonlocal models. To unlock the full potential of nonlocal modeling we developed a mathematically rigorous and physically consistent interface theory and demonstrate its scope in mission-relevant exemplar problems.

97 MATHEMATICS AND COMPUTING↗

Hot Isostatic Pressing for Reduction of Defects in Additively Manufactured Inconel-718 Preliminary Simulation Approach and Results

Selective laser melting (SLM) of Inconel-718 (IN718) is a promising advanced manufacturing technology; however, challenges exist for fracture-critical parts where pores induced during the manufacturing process reduce fatigue life of the material. Hot isostatic pressing (HIP) is a common post processing technique used to reduce the part porosity by exposing parts to high temperature and pressures. Pore closure depends on the pore shape, trapped gas, and HIP parameters, and the process can only be optimized by understanding the effect of these variables. An analysis capability is described to better understand the conditions necessary for pore closure, and the effect of several variables on porosity are explored and compared to a baseline case. The importance of valid material response models is identified as a significant factor in model accuracy; especially the creep material response model which, in the baseline case, predicts that creep strain accounts for 87.7% of the total strain at the conclusion of the HIP cycle. Lessons learned are described, in addition to future work including model validation.

Joshua M Fody↗

Experimental and Numerical Investigations of Cyclic Plastic Deformation of Al-Mg Alloy

In this work, low-cycle fatigue (LCF) tests are performed on as-received and solution-treated Al-Mg alloy (AA 5754). Solution heat treatment of the as-received alloy is performed at 530 °C for 2 h and then quenched in water at room temperature. The low-cycle fatigue response of the as-received and solution-treated alloy is evaluated at five different strain amplitudes from 0.2 to 0.6%. The strain rate of 0.5 × 10{sup −3} s{sup −1} is kept constant in all the LCF tests. The LCF behavior of AA 5754 is modeled by finite element method (FEM) using Chaboche kinematic hardening cyclic plasticity model. Hardening behavior or fatigue life of as-received and solution-treated alloy is successfully simulated using this model. The results obtained by FEM simulations are compared with those obtained by LCF experiments. Microstructural studies are also performed for as-received and solution-treated alloy. The fracture surfaces obtained from LCF tests are analyzed by scanning electron microscopy.

36 MATERIALS SCIENCE↗

Concurrent atomistic-continuum modeling of crystalline materials

In the work reported here, we present a concurrent atomistic-continuum (CAC) method for modeling and simulation of crystalline materials. The CAC formulation extends the Irving-Kirkwood procedure for deriving transport equations and fluxes for homogenized molecular systems to that for polyatomic crystalline materials by employing a concurrent two-level description of the structure and dynamics of crystals. A multiscale representation of conservation laws is formulated, as a direct consequence of Newton's second law, in terms of instantaneous expressions of unit cell-averaged quantities using the mathematical theory of distributions. Finite element (FE) solutions to the conservation equations, as well as fluxes and temperature in the FE representation, are introduced, followed by numerical examples of the atomic-scale structure of interfaces, dynamics of fracture and dislocations, and phonon thermal transport across grain boundaries. In addition to providing a methodology for concurrent multiscale simulation of transport processes under a single theoretical framework, the CAC formulation can also be used to compute fluxes (stress and heat flux) in atomistic and coarse-grained atomistic simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multi-scale mechanical modeling of Al-steel resistance spot welds

A multi-scale finite element modeling approach was developed to study the deformation and fracture behavior of Al-steel resistance spot welds. First, a micro-scale model was applied to simulate the mechanical responses of the intermetallic compound (IMC) layer having various morphologies and thicknesses under tensile and shear loading conditions. Second, the predicted tensile and shear strength of the IMC layer, that varied along the joint interface per the IMC layer morphology and thickness variation, was then introduced into 3D macro-scale models to predict the overall mechanical performance of weld coupons under coach peel, lap shear and cross-tension testing conditions. In conclusion, the numerical predictions agreed reasonably well with the experimental data.

36 MATERIALS SCIENCE↗

Effects of Atomic-Scale Structure on the Fracture Properties of Amorphous Carbon - Carbon Nanotube Composites

The fracture of carbon materials is a complex process, the understanding of which is critical to the development of next generation high performance materials. While quantum mechanical (QM) calculations are the most accurate way to model fracture, the fracture behavior of many carbon-based composite engineering materials, such as carbon nanotube (CNT) composites, is a multi-scale process that occurs on time and length scales beyond the practical limitations of QM methods. The Reax Force Field (ReaxFF) is capable of predicting mechanical properties involving strong deformation, bond breaking and bond formation in the classical molecular dynamics framework. This has been achieved by adding to the potential energy function a bond-order term that varies continuously with distance. The use of an empirical bond order potential, such as ReaxFF, enables the simulation of failure in molecular systems that are several orders of magnitude larger than would be possible in QM techniques. In this work, the fracture behavior of an amorphous carbon (AC) matrix reinforced with CNTs was modeled using molecular dynamics with the ReaxFF reactive forcefield. Care was taken to select the appropriate simulation parameters, which can be different from those required when using traditional fixed-bond force fields. The effect of CNT arrangement was investigated with three systems: a single-wall nanotube (SWNT) array, a multi-wall nanotube (MWNT) array, and a SWNT bundle system. For each arrangement, covalent bonds are added between the CNTs and AC, with crosslink fractions ranging from 0-25% of the interfacial CNT atoms. The SWNT and MWNT array systems represent ideal cases with evenly spaced CNTs; the SWNT bundle system represents a more realistic case because, in practice, van der Waals interactions lead to the agglomeration of CNTs into bundles. The simulation results will serve as guidance in setting experimental processing conditions to optimize the mechanical properties of CNT composites.

Jensen, Benjamin D.↗

Pellet cladding mechanical interaction modeling using the extended finite element method

Fracturing of light water reactor (LWR) UO{sub 2} pellet has multiple important effects on fuel performance. One particularly important concern is that cracks in the fuel cause elevated stresses in the cladding when pellet cladding mechanical interaction (PCMI) occurs. The extended finite element method (XFEM) is a powerful method to represent arbitrary propagating discrete cracks in finite element models. The use of XFEM has been previously demonstrated for modeling propagating discrete cracks in the BISON fuel performance code. This paper demonstrates an initial application of XFEM to model stress concentrations induced by fuel fractures at the fuel/cladding interface during PCMI. This is demonstrated on a study of a pre-defined stationary crack in a 2D cross-section model of a LWR fuel rod. The results from a model with a discrete crack defined with XFEM compare favorably with the results from a model with the same crack geometry defined in the finite element mesh. This study focuses on benchmarking the use of XFEM for PCMI modeling with a stationary crack, but this technique will be readily extended in the future to consider PCMI in conjunction with arbitrary, physics-driven crack propagation. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Selection of Composites Simulation Practices at NASA Langley Research Center

One of the major areas of study at NASA Langley Research Center is the development of technologies that support the use of advanced composite materials in aerospace applications. Amongst the supporting technologies are analysis tools used to simulate the behavior of these materials. This presentation will discuss a number of examples of analysis tools and simulation practices conducted at NASA Langley. The presentation will include examples of damage tolerance analyses for both interlaminar and intralaminar failure modes. Tools for modeling interlaminar failure modes include fracture mechanics and cohesive methods, whilst tools for modeling intralaminar failure involve the development of various progressive failure analyses. Other examples of analyses developed at NASA Langley include a thermo-mechanical model of an orthotropic material and the simulation of delamination growth in z-pin reinforced laminates.

Ratcliffe, James G.↗

Flight monitor for jet engine disk cracks and the use of critical length criterion of fracture mechanics

A disk crack detector is discussed which is intended to operate under flight conditions. It monitors the disk rim for surface cracks emanating from the blade root interface. An eddy current type sensor, with a remotely located capacitance/conductance bridge and signal analyzer, can reliably detect a simulated crack 3 mm long. The sensor was tested on a spinning turbine disk at 540 C. Tests indicate that the system is useful at disk rim velocities to 460 m/sec. By using fracture mechanics, it is shown for Inconel 718 th at a crack operating under a rim stress of 34 x ten to the 7th power N/sqm has a critical length of 18 mm.

Barranger, J. P.↗

Simulation of Flow in Multi-Scale Porous Media Using the Lattice Boltzmann Method on Quadtree Grids

Abstract The unified lattice Boltzmann model is extended to the quadtree grids for simulation of fluid flow through porous media. The unified lattice Boltzmann model is capable of simulating flow in porous media at various scales or in systems where multiple length scales coexist. The quadtree grid is able to provide a high-resolution approximation to complex geometries, with great flexibility to control local grid density. The combination of the unified lattice Boltzmann model and the quadtree grids results in an efficient numerical model for calculating permeability of multi-scale porous media. The model is used for permeability calculation for three systems, including a fractured system used in a previous study, a Voronoi tessellation system, and a computationally-generated pore structure of fractured shale. The results are compared with those obtained using the conventional lattice Boltzmann model or the unified lattice Boltzmann model on rectangular or uniform square grid. It is shown that the proposed model is an accurate and efficient tool for flow simulation in multi-scale porous media. In addition, for the fractured shale, the contribution of flow in matrix and fractures to the overall permeability of the fractured shale is studied systematically.

Earth Sciences↗

Deformation and Fracture in Complex-Shaped Energetic Particles

A new technical basis on the mechanics of energetic materials at the individual particle scale has been developed. Despite these particles being in most of our Sandia non-nuclear explosive components, we have historically lacked any understanding of particle behavior. Through the novel application of nanoidentation methods to single crystal films and single particles of energetic materials with complex shapes, discovery data has been collected elucidating phenomena of particle strength, elastic and plastic deformation, and fracture. This work specifically developed the experimental techniques and analysis methodologies to distill data into relationships suitable for future integration into particle level simulations of particle reassembly. This project utilized experimental facilities at CINT and the Explosive Components Facility to perform ex-situ and in-situ nanoidentation experiments with simultaneous scanning electron microscope (SEM) imaging. Data collected by an applied axial compressive load in either force-control or displacement-control was well represented by Hertzian contact theory for linear elastic materials. Particle fracture phenomenology was effectively modeled by an empirical damage model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Nanoindentation Study of the Plastic Deformation and Fracture Mechanisms in Single-Crystalline CaFe 2 As 2

In this paper, the plastic deformation and fracture mechanisms in single-crystalline CaFe 2 As 2 has been studied using nanoindentation and density functional theory simulations. CaFe 2 As 2 single crystals were grown in a Sn-flux, resulting in homogeneous and nearly defect-free crystals. Nanoindentation along the [001] direction produces strain bursts, radial cracking, and lateral cracking. Ideal cleavage simulations along the [001] and [100] directions using density functional theory calculations revealed that cleavage along the [001] direction requires a much lower stress than cleavage along the [100] direction. This strong anisotropy of cleavage strength implies that CaFe 2 As 2 has an atomic-scale layered structure, which typically exhibits lateral cracking during nanoindentation. This special layered structure results from weak atomic bonding between the (001) Ca and Fe 2 As 2 layers.

36 MATERIALS SCIENCE↗