Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture modelling”

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 613 records · Page 34

Multiscale Modeling of Damage Processes in Aluminum Alloys: Grain-Scale Mechanisms

This paper has two goals related to the development of a physically-grounded methodology for modeling the initial stages of fatigue crack growth in an aluminum alloy. The aluminum alloy, AA 7075-T651, is susceptible to fatigue cracking that nucleates from cracked second phase iron-bearing particles. Thus, the first goal of the paper is to validate an existing framework for the prediction of the conditions under which the particles crack. The observed statistics of particle cracking (defined as incubation for this alloy) must be accurately predicted to simulate the stochastic nature of microstructurally small fatigue crack (MSFC) formation. Also, only by simulating incubation of damage in a statistically accurate manner can subsequent stages of crack growth be accurately predicted. To maintain fidelity and computational efficiency, a filtering procedure was developed to eliminate particles that were unlikely to crack. The particle filter considers the distributions of particle sizes and shapes, grain texture, and the configuration of the surrounding grains. This filter helps substantially reduce the number of particles that need to be included in the microstructural models and forms the basis of the future work on the subsequent stages of MSFC, crack nucleation and microstructurally small crack propagation. A physics-based approach to simulating fracture should ultimately begin at nanometer length scale, in which atomistic simulation is used to predict the fundamental damage mechanisms of MSFC. These mechanisms include dislocation formation and interaction, interstitial void formation, and atomic diffusion. However, atomistic simulations quickly become computationally intractable as the system size increases, especially when directly linking to the already large microstructural models. Therefore, the second goal of this paper is to propose a method that will incorporate atomistic simulation and small-scale experimental characterization into the existing multiscale framework. At the microscale, the nanoscale mechanics are represented within cohesive zones where appropriate, i.e. where the mechanics observed at the nanoscale can be represented as occurring on a plane such as at grain boundaries or slip planes at a crack front. Important advancements that are yet to be made include: 1. an increased fidelity in cohesive zone modeling; 2. a means to understand how atomistic simulation scales with time; 3. a new experimental methodology for generating empirical models for CZMs and emerging materials; and 4. a validation of simulations of the damage processes at the nano-micro scale. With ever-increasing computer power, the long-term ability to employ atomistic simulation for the prognosis of structural components will not be limited by computation power, but by our lack of knowledge in incorporating atomistic models into simulations of MSFC into a multiscale framework.

Hochhalter, J. D.↗

Thermal evolution of a differentiated Ganymede and implications for surface features

Thermodynamic models are developed for the processes which controlled the evolution of the surface Ganymede, an icy Jovian satellite assumed to have a rock-rich core surrounded by a water-ice mantle. Account is taken of a heat pulse which would have arisen from a Rayleigh-Taylor instability at a deep-seated liquid-solid water interface, rapid fracturing from global stresses imposed by warm ice diapiric upwelling, impacts by large meteorites, and resurfacing by ice flows (rather than core formation). Comparisons are made with existing models for the evolution of Callisto, and the difficulties in defining a mechanism which produced the groove terrain of Ganymede are discussed.

Kirk, R. L.↗

Face Sheet/Core Disbond Growth in Honeycomb Sandwich Panels Subjected to Ground-Air-Ground Pressurization and In-Plane Loading

Typical damage modes in light honeycomb sandwich structures include face sheet/core disbonding and core fracture, both of which can pose a threat to the structural integrity of a component. These damage modes are of particular interest to aviation certification authorities since several in-service occurrences, such as rudder structural failure and other control surface malfunctions, have been attributed to face sheet/core disbonding. Extensive studies have shown that face sheet/core disbonding and core fracture can lead to damage propagation caused by internal pressure changes in the core. The increasing use of composite sandwich construction in aircraft applications makes it vitally important to understand the effect of ground-air-ground (GAG) cycles and conditions such as maneuver and gust loads on face sheet/core disbonding. The objective of the present study was to use a fracture mechanics based approach developed earlier to evaluate the loading at the disbond front caused by ground-air-ground pressurization and in-plane loading. A honeycomb sandwich panel containing a circular disbond at one face sheet/core interface was modeled with three-dimensional (3D) solid finite elements. The disbond was modeled as a discrete discontinuity and the strain energy release rate along the disbond front was computed using the Virtual Crack Closure Technique (VCCT). Special attention was paid to the pressure-deformation coupling which can decrease the pressure load within the disbonded sandwich section significantly when the structure is highly deformed. The commercial finite element analysis software, Abaqus/Standard, was used for the analyses. The recursive pressure-deformation coupling problem was solved by representing the entrapped air in the honeycomb cells as filled cavities in Abaqus/Standard. The results show that disbond size, face sheet thickness and core thickness are important parameters that determine crack tip loading at the disbond front. Further, the pressure-deformation coupling was found to have an important load decreasing effect [6]. In this paper, a detailed problem description is provided first. Second, the analysis methodology is presented. The fracture mechanics approach used is described and the specifics of the finite element model, including the fluid-filled cavities, are introduced. Third, the initial model verification and validation are discussed. Fourth, the findings from a closely related earlier study [6] are summarized. These findings provided the basis for the current investigation. Fifth, an aircraft ascent scenario from 0 to 12192 m (0 to 40000 ft) is considered and the resulting crack tip loading at the disbond front is determined. In-plane loading to simulate maneuvers and gust conditions are also considered. Sixth, the results are shown for a curved panel, which was used to simulate potential fuselage applications. Finally, a brief summary of observations is presented and recommendations for improvement are provided.

Chen, Zhi M.↗

Microstructurally based model of fatigue initiation and growth

A model to calculate fatigue life is developed based on the assumption that fatigue life is entirely composed of crack growth from an initial microstructural inhomogeneity. Specifically, growth is considered to start from either an ellipsoidal void, a cracked particle, or a debonded particle. The capability of predicting fatigue life from material microstructure is based on linear elastic fracture mechanics principles, the sizes of the crack-initiating microstructural inhomogeneities, and an initiation parameter that is proportional to the cyclic plastic zone size. A key aspect of this modeling approach is that it is linked with a general purpose probability program to analyze the effect of the distribution of controlling microstructural features within the material. This enables prediction of fatigue stress versus life curves for various specimen geometries using distributional statistics obtained from characterizations of the microstructure. Results are compared to experimental fatigue data from an aluminum alloy.

Brockenbrough, J. R.↗

Ceramics Analysis and Reliability Evaluation of Structures (CARES). Users and programmers manual

This manual describes how to use the Ceramics Analysis and Reliability Evaluation of Structures (CARES) computer program. The primary function of the code is to calculate the fast fracture reliability or failure probability of macroscopically isotropic ceramic components. These components may be subjected to complex thermomechanical loadings, such as those found in heat engine applications. The program uses results from MSC/NASTRAN or ANSYS finite element analysis programs to evaluate component reliability due to inherent surface and/or volume type flaws. CARES utilizes the Batdorf model and the two-parameter Weibull cumulative distribution function to describe the effect of multiaxial stress states on material strength. The principle of independent action (PIA) and the Weibull normal stress averaging models are also included. Weibull material strength parameters, the Batdorf crack density coefficient, and other related statistical quantities are estimated from four-point bend bar or unifrom uniaxial tensile specimen fracture strength data. Parameter estimation can be performed for single or multiple failure modes by using the least-square analysis or the maximum likelihood method. Kolmogorov-Smirnov and Anderson-Darling goodness-of-fit tests, ninety percent confidence intervals on the Weibull parameters, and Kanofsky-Srinivasan ninety percent confidence band values are also provided. The probabilistic fast-fracture theories used in CARES, along with the input and output for CARES, are described. Example problems to demonstrate various feature of the program are also included. This manual describes the MSC/NASTRAN version of the CARES program.

Nemeth, Noel N.↗

Design of ceramic components with the NASA/CARES computer program

The ceramics analysis and reliability evaluation of structures (CARES) computer program is described. The primary function of the code is to calculate the fast-fracture reliability or failure probability of macro-scopically isotropic ceramic components. These components may be subjected to complex thermomechanical loadings, such as those found in heat engine applications. CARES uses results from MSC/NASTRAN or ANSYS finite-element analysis programs to evaluate how inherent surface and/or volume type flaws component reliability. CARES utilizes the Batdorf model and the two-parameter Weibull cumulative distribution function to describe the effects of multiaxial stress states on material strength. The principle of independent action (PIA) and the Weibull normal stress averaging models are also included. Weibull material strength parameters, the Batdorf crack density coefficient, and other related statistical quantities are estimated from four-point bend bar or uniform uniaxial tensile specimen fracture strength data. Parameter estimation can be performed for a single or multiple failure modes by using a least-squares analysis or a maximum likelihood method. Kolmogorov-Smirnov and Anderson-Darling goodness-to-fit-tests, 90 percent confidence intervals on the Weibull parameters, and Kanofsky-Srinivasan 90 percent confidence band values are also provided. Examples are provided to illustrate the various features of CARES.

Nemeth, Noel N.↗

Ceramic component reliability with the restructured NASA/CARES computer program

The Ceramics Analysis and Reliability Evaluation of Structures (CARES) integrated design program on statistical fast fracture reliability and monolithic ceramic components is enhanced to include the use of a neutral data base, two-dimensional modeling, and variable problem size. The data base allows for the efficient transfer of element stresses, temperatures, and volumes/areas from the finite element output to the reliability analysis program. Elements are divided to insure a direct correspondence between the subelements and the Gaussian integration points. Two-dimensional modeling is accomplished by assessing the volume flaw reliability with shell elements. To demonstrate the improvements in the algorithm, example problems are selected from a round-robin conducted by WELFEP (WEakest Link failure probability prediction by Finite Element Postprocessors).

Powers, Lynn M.↗

Ceramic component reliability with the restructured NASA/CARES computer program

The Ceramics Analysis and Reliability Evaluation of Structures (CARES) integrated design program on statistical fast fracture reliability and monolithic ceramic components is enhanced to include the use of a neutral data base, two-dimensional modeling, and variable problem size. The data base allows for the efficient transfer of element stresses, temperatures, and volumes/areas from the finite element output to the reliability analysis program. Elements are divided to insure a direct correspondence between the subelements and the Gaussian integration points. Two-dimensional modeling is accomplished by assessing the volume flaw reliability with shell elements. To demonstrate the improvements in the algorithm, example problems are selected from a round-robin conducted by WELFEP (WEakest Link failure probability prediction by Finite Element Postprocessors).

Powers, Lynn M.↗

Application of Video Image Correlation Techniques to the Space Shuttle External Tank Foam Materials

Results that illustrate the use of a video-image-correlation-based displacement and strain measurement system to assess the effects of material nonuniformities on the behavior of the sprayed-on foam insulation (SOFI) used for the thermal protection system on the Space Shuttle External Tank are presented. Standard structural verification specimens for the SOFI material with and without cracks and subjected to mechanical or thermal loading conditions were tested. Measured full-field displacements and strains are presented for selected loading conditions to illustrate the behavior of the foam and the viability of the measurement technology. The results indicate that significant strain localization can occur in the foam because of material nonuniformities. In particular, elongated cells in the foam can interact with other geometric or material discontinuities in the foam and develop large-magnitude localized strain concentrations that likely initiate failures. Furthermore, some of the results suggest that continuum mechanics and linear elastic fracture mechanics might not adequately represent the physical behavior of the foam, and failure predictions based on homogeneous linear material models are likely to be inadequate.

Hilburger, Mark W.↗

Application of Video Image Correlation Techniques to the Space Shuttle External Tank Foam Materials

Results that illustrate the use of a video-image-correlation-based displacement and strain measurement system to assess the effects of material nonuniformities on the behavior of the sprayed-on foam insulation (SOFI) used for the thermal protection system on the Space Shuttle External Tank are presented. Standard structural verification specimens for the SOFI material with and without cracks and subjected to mechanical or thermal loading conditions were tested. Measured full-field displacements and strains are presented for selected loading conditions to illustrate the behavior of the foam and the viability of the measurement technology. The results indicate that significant strain localization can occur in the foam because of material nonuniformities. In particular, elongated cells in the foam can interact with other geometric or material discontinuities in the foam and develop large-magnitude localized strain concentrations that likely initiate failures. Furthermore, some of the results suggest that continuum mechanics and linear elastic fracture mechanics might not adequately represent the physical behavior of the foam, and failure predictions based on homogeneous linear material models are likely to be inadequate.

Hilburger, Mark W.↗

Effect of Data Reduction and Fiber-Bridging on Mode I Delamination Characterization of Unidirectional Composites

Reliable delamination characterization data for laminated composites are needed for input in analytical models of structures to predict delamination onset and growth. The double-cantilevered beam (DCB) specimen is used to measure fracture toughness, GIc, and strain energy release rate, GImax, for delamination onset and growth in laminated composites under mode I loading. The current study was conducted as part of an ASTM Round Robin activity to evaluate a proposed testing standard for Mode I fatigue delamination propagation. Static and fatigue tests were conducted on specimens of IM7/977-3 and G40-800/5276-1 graphite/epoxies, and S2/5216 glass/epoxy DCB specimens to evaluate the draft standard "Standard Test Method for Mode I Fatigue Delamination Propagation of Unidirectional Fiber-Reinforced Polymer Matrix Composites." Static results were used to generate a delamination resistance curve, GIR, for each material, which was used to determine the effects of fiber-bridging on the delamination growth data. All three materials were tested in fatigue at a cyclic GImax level equal to 90% of the fracture toughness, GIc, to determine the delamination growth rate. Two different data reduction methods, a 2-point and a 7-point fit, were used and the resulting Paris Law equations were compared. Growth rate results were normalized by the delamination resistance curve for each material and compared to the nonnormalized results. Paris Law exponents were found to decrease by 5.4% to 46.2% due to normalizing the growth data. Additional specimens of the IM7/977-3 material were tested at 3 lower cyclic GImax levels to compare the effect of loading level on delamination growth rates. The IM7/977-3 tests were also used to determine the delamination threshold curve for that material. The results show that tests at a range of loading levels are necessary to describe the complete delamination behavior of this material.

Murri, Gretchen B.↗

Evaluation of Delamination Growth Characterization Methods Under Mode I Fatigue Loading

Reliable delamination characterization data for laminated composites are needed for input to analytical models of structures to predict delamination. The double-cantilevered beam (DCB) specimen is used with laminated composites to measure fracture toughness, G(sub Ic), delamination onset strain energy release rate, and growth rate data under cyclic loading. In the current study, DCB specimens of IM7/8552 graphite/epoxy supplied by two different manufacturers were tested in static and fatigue to compare the measured characterization data from the two sources, and to evaluate a proposed ASTM standard for generating Paris Law equations. Static results were used to generate compliance calibration constants for the fatigue data, and a delamination resistance curve, G(sub IR), which was used to determine the effects of fiber-bridging on delamination growth. Specimens were tested in fatigue at a cyclic G(sub Imax) level equal to 50, 40 or 30% of G(sub Ic), to determine a delamination onset curve and delamination growth rate. The delamination onset curve equations had similar exponents and the same trends. Delamination growth rate was calculated by fitting a Paris Law to the da/dN versus G(sub Imax) data. Both a 2-point and a 7-point data reduction method were used and the Paris Law equations were compared. To determine the effects of fiber-bridging, growth rate results were normalized by the delamination resistance curve for each material and compared to the non-normalized results. Paris Law exponents were found to decrease by 31% to 37% due to normalizing the growth data. Normalizing the data also greatly reduced the amount of scatter between the different specimens. Visual data records from the fatigue testing were used to calculate individual compliance calibration constants from the fatigue data for some of the specimens. The resulting da/dN versus G(sub Imax) plots showed much improved repeatability between specimens. Gretchen

Murri, Gretchen B.↗

Computational Analysis of Artificial Gravity as a Possible Countermeasure to Spaceflight Induced Bone Loss

During exploration class missions, such as to asteroids and Mars, astronauts will be exposed to reduced gravity for extended periods. Data has shown that astronauts lose bone mass at a rate of 1% to 2% a month in microgravity, particularly in lower extremities such as the proximal femur. Exercise countermeasures have not completely eliminated bone loss from long duration spaceflight missions, which leaves astronauts susceptible to early onset osteoporosis and greater risk of fracture. Introduction of the Advanced Resistive Exercise Device and other large exercise devices on the International Space Station (ISS), coupled with improved nutrition, has further minimized bone loss. However, unlike the ISS, exploration vehicles will have very limited volume and power available to accommodate such capabilities. Therefore, novel concepts like artificial gravity systems are being explored as a means to provide sufficient load stimulus to the musculoskeletal system to mitigate bone changes that may lead to early onset osteoporosis and increased risk of fracture. Currently, there is minimal data available to drive further research and development efforts to appropriately explore such options. Computational modeling can be leveraged to gain insight on the level of osteoprotection that may be achieved using artificial gravity produced by a spinning spacecraft or centrifuge. With this in mind, NASA's Digital Astronaut Project (DAP) has developed a bone remodeling model that has been validated for predicting volumetric bone mineral density (vBMD) changes of trabecular and cortical bone both for gravitational unloading condition and the equivalent of 1g daily load stimulus. Using this model, it is possible to simulate vBMD changes in trabecular and cortical bone under different gravity conditions. In this presentation, we will discuss our preliminary findings regarding if and how artificial gravity may be used to mitigate spaceflight induced bone loss.

Mulugeta, L.↗

Disruption of Kilometer-Sized Asteroids by Energetic Collisions

Past collisions determine the current mechanical properties of asteroids, and hence their response to future collisions. To understand the evolution of their diverse shapes, structures and spin states, we model hypervelocity impacts into small asteroids using a modern 3D smooth-particle hydrocode (SPH3D) that takes material strength and fracture into account.

Asteroids↗

Composite Bond Line Measurements Based on A Bayesian Analysis of Flash Thermography Data

For bonded composite materials, an accurate characterization of the adhesive bond line is needed to predict failure modes and fracture toughness. In this paper, bond line thickness was estimated from data obtained using through transmission flash thermography. The forward model that predicts back surface temperature is based on a three layer heat diffusion equation with varying diffusivity and flux boundary conditions. The corresponding inverse problem of estimating bond line thickness from measurement data was solved using a Bayesian approach that assumed Gaussian priors for the bond line thickness and thermal diffusivity of the adherends. Finally, the outputs of the thermography based method were compared to measurements that were collected using a micrometer and ultrasound testing.

Flash thermography↗

Composite Bond Line Measurements Based on a Bayesian Analysis of Flash Thermography Data

For bonded composite materials, an accurate characterization of the adhesive bond line is needed to predict failure modes and fracture toughness. In this paper, bond line thickness was estimated from data obtained using through transmission flash thermography. The forward model that predicts back surface temperature is based on a three layer heat diffusion equation with varying diffusivity and flux boundary conditions. The corresponding inverse problem of estimating bond line thickness from measurement data was solved using a Bayesian approach that assumed Gaussian priors for the bond line thickness and thermal diffusivity of the adherends. Finally, the outputs of the thermography based method were compared to measurements that were collected using a micrometer and ultrasound testing.

Flash thermography↗

A pressurized cylindrical shell with a fixed end which contains an axial part-through or through crack

A cylindrical shell having a very stiff and plate or a flange is considered. It is assumed that near the end the cylinder contains an axial flaw which may be modeled as a part through surface crack or a through crack. The effect of the end constraining on the stress intensity factor which is the main fracture mechanics parameter is studied. The applied loads acting on the cylinder are assumed to be axisymmetric. Thus the crack problem under consideration is symmetric with respect to the plane of the crack and consequently only the Mode 1 stress intensity factors are nonzero. With this limitation, the general perturbation problem for a cylinder with a built in end containing an axial crack is considered. Reissner's shell theory is used to formulate the problem. The part through crack problem is treated by using a line spring model. In the case of a crack tip terminating at the fixed end it is shown that the integral equations of the shell problem has the same generalized Cauchy kernel as the corresponding plane stress elasticity problem.

Yahsi, O. S.↗

Micromechanics, Fracture Mechanics and Gas Permeability of Composite Laminates for Cryogenic Storage Systems

A micromechanics method is developed to investigate microcrack propagation in a liquid hydrogen composite tank at cryogenic temperature. The unit cell is modeled using square and hexagonal shapes depends on fiber and matrix layout from microscopic images of composite laminates. Periodic boundary conditions are applied to the unit cell. The temperature dependent properties are taken into account in the analysis. The laminate properties estimated by the micromechanics method are compared with empirical solutions using constituent properties. The micro stresses in the fiber and matrix phases based on boundary conditions in laminate level are calculated to predict the formation of microcracks in the matrix. The method is applied to an actual liquid hydrogen storage system. The analysis predicts micro stresses in the matrix phase are large enough to cause microcracks in the composite. Stress singularity of a transverse crack normal to a ply-interface is investigated to predict the fracture behavior at cryogenic conditions using analytical and finite element analysis. When a transverse crack touches a ply-interface of a composite layer with same fiber orientation, the stress singularity is equal to 1/2. When the transverse crack propagates to a stiffer layer normal to the ply-direction, the singularity becomes less than 1/2 and vice versa. Finite element analysis is performed to predict the fracture toughness of a laminated beam subjected to fracture loads measured by four-point bending tests at room and cryogenic temperatures. As results, the fracture load at cryogenic temperature is significantly lower than that at room temperature. However, when thermal stresses are taken into consideration, for both cases of room and cryogenic temperatures, the difference of the fracture toughness becomes insignificant. The result indicates fracture toughness is a characteristic property, which is independent to temperature changes. The experimental analysis is performed to investigate the effect of cryogenic cycling on permeability for various composite material systems. Textile composites have lower permeability than laminated composites even with increasing number of cryogenic cycle. Nano-particles dispersed in laminated composites do not show improvement on permeability. The optical inspection is performed to investigate the microcrack propagation and void content in laminated composites and compared the microscopic results before and after cryogenic cycling.

Choi, Sukjoo↗