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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Mapping Process Model Results to Fracture Model Initial Conditions for Adhesive Bonding

A systematic approach is proposed for mapping adhesive bonding process outcomes to initial conditions for progressive damage analysis (PDA). Herein, two mapping procedures are developed: 1) direct mapping for residual stresses, strains, and deformations obtained from a process model; and 2) functional mapping for quantities such as bondline thickness, degree of cure, and porosity that are assumed to have a functional relationship with fracture toughness but are not discretely modeled in the PDA. The spatial distribution of functionally mapped quantities may be determined by process modeling or inspection. The functional maps between the quantity of interest and fracture toughness may be obtained from independent lower-scale numerical simulations or empirical test campaigns. The proposed mapping procedure links the adhesive bonding process to the structural performance. Therefore, by accounting for the effects of bondline nonuniformities in structural analysis, analytical predictive accuracy can be improved, and off-nominal conditions can be evaluated. The mapping procedures are demonstrated and verified with example problems.

Andrew C Bergan↗
NASA NTRSDate not supplied

Mapping Process Model Results to Fracture Model Initial Conditions for Adhesive Bonding

A systematic approach is proposed for mapping adhesive bonding process outcomes to initial conditions for progressive damage analysis (PDA). Herein, two mapping procedures are developed: 1) direct mapping for residual stresses, strains, and deformations obtained from a process model; and 2) functional mapping for quantities such as bondline thickness, degree of cure, and porosity that are assumed to have a functional relationship with fracture toughness but are not discretely modeled in the PDA. The spatial distribution of functionally mapped quantities may be determined by process modeling or inspection. The functional maps between the quantity of interest and fracture toughness may be obtained from independent lower-scale numerical simulations or empirical test campaigns. The proposed mapping procedure links the adhesive bonding process to the structural performance. Therefore, by accounting for the effects of bondline nonuniformities in structural analysis, analytical predictive accuracy can be improved, and off-nominal conditions can be evaluated. The mapping procedures are demonstrated and verified with example problems.

Andrew Bergan↗
NASA NTRS1993

Explicit 3D continuum fracture modeling with smooth particle hydrodynamics

Impact phenomena shaped our solar system. As usual for most solar system processes, the scales are far different than we can address directly in the laboratory. Impact velocities are often much higher than we can achieve, sizes are often vastly larger, and most impacts take place in an environment where the only gravitational force is the mutual pull of the impactors. The Smooth Particle Hydrodynamics (SPH) technique has been applied in the past to the simulations of giant impacts. In these simulations, the colliding objects were so massive (at least a sizeable fraction of the Earth's mass) that material strength was negligible compared to gravity. This assumption can no longer be made when the bodies are much smaller. To this end, we have developed a 3D SPH code that includes a strength model to which we have added a von Mises yielding relation for stresses beyond the Hugoniot Elastic Limit. At the lower stresses associated with brittle failure, we use a rate-dependent strength based on the nucleation of incipient flaws whose number density is given by a Weibull distribution. Following Grady and Kipp and Melosh et al., we introduce a state variable D ('damage'), 0 less than D less than 1, which expresses the local reduction in strength due to crack growth under tensile loading. Unfortunately for the hydrodynamics, Grady and Kipp's model predicts which fragments are the most probable ones and not the ones that are really formed. This means, for example, that if a given laboratory experiment is modeled, the fragment distribution obtained from the Grady-Kipp theory would be equivalent to a ensemble average over many realizations of the experiment. On the other hand, the hydrodynamics itself is explicit and evolves not an ensemble average but very specific fragments. Hence, there is a clear incompatibility with the deterministic nature of the hydrodynamics equations and the statistical approach of the Grady-Kipp dynamical fracture model. We remedy these shortcomings by making the incipient flaw distribution explicit, i.e., particles carry activation strains which are distributed at random with a probability of occurrence given by the Weibull distribution. If the local principal axis strain exceeds this limit, damage starts to grow. By growing explicit cracks together with statistical cracks (damage) at the sub-particle scale, we ensure that material strength and fragmentation is independent of model resolution. We tested our scheme by simulating laboratory impact experiments on basalt spheres.

Benz, W.↗
NASA NTRS2013

Model Verification and Validation Concepts for a Probabilistic Fracture Assessment Model to Predict Cracking of Knife Edge Seals in the Space Shuttle Main Engine High Pressure Oxidizer

Physics-based models are routinely used to predict the performance of engineered systems to make decisions such as when to retire system components, how to extend the life of an aging system, or if a new design will be safe or available. Model verification and validation (V&V) is a process to establish credibility in model predictions. Ideally, carefully controlled validation experiments will be designed and performed to validate models or submodels. In reality, time and cost constraints limit experiments and even model development. This paper describes elements of model V&V during the development and application of a probabilistic fracture assessment model to predict cracking in space shuttle main engine high-pressure oxidizer turbopump knife-edge seals. The objective of this effort was to assess the probability of initiating and growing a crack to a specified failure length in specific flight units for different usage and inspection scenarios. The probabilistic fracture assessment model developed in this investigation combined a series of submodels describing the usage, temperature history, flutter tendencies, tooth stresses and numbers of cycles, fatigue cracking, nondestructive inspection, and finally the probability of failure. The analysis accounted for unit-to-unit variations in temperature, flutter limit state, flutter stress magnitude, and fatigue life properties. The investigation focused on the calculation of relative risk rather than absolute risk between the usage scenarios. Verification predictions were first performed for three units with known usage and cracking histories to establish credibility in the model predictions. Then, numerous predictions were performed for an assortment of operating units that had flown recently or that were projected for future flights. Calculations were performed using two NASA-developed software tools: NESSUS(Registered Trademark) for the probabilistic analysis, and NASGRO(Registered Trademark) for the fracture mechanics analysis. The goal of these predictions was to provide additional information to guide decisions on the potential of reusing existing and installed units prior to the new design certification.

Pai, Shantaram S.↗
NASA NTRS2011

Investigating Some Technical Issues on Cohesive Zone Modeling of Fracture

This study investigates some technical issues related to the use of cohesive zone models (CZMs) in modeling fracture processes. These issues include: why cohesive laws of different shapes can produce similar fracture predictions; under what conditions CZM predictions have a high degree of agreement with linear elastic fracture mechanics (LEFM) analysis results; when the shape of cohesive laws becomes important in the fracture predictions; and why the opening profile along the cohesive zone length needs to be accurately predicted. Two cohesive models were used in this study to address these technical issues. They are the linear softening cohesive model and the Dugdale perfectly plastic cohesive model. Each cohesive model constitutes five cohesive laws of different maximum tractions. All cohesive laws have the same cohesive work rate (CWR) which is defined by the area under the traction-separation curve. The effects of the maximum traction on the cohesive zone length and the critical remote applied stress are investigated for both models. For a CZM to predict a fracture load similar to that obtained by an LEFM analysis, the cohesive zone length needs to be much smaller than the crack length, which reflects the small scale yielding condition requirement for LEFM analysis to be valid. For large-scale cohesive zone cases, the predicted critical remote applied stresses depend on the shape of cohesive models used and can significantly deviate from LEFM results. Furthermore, this study also reveals the importance of accurately predicting the cohesive zone profile in determining the critical remote applied load.

Wang, John T.↗
NASA NTRS2017

Sensitivity Analysis of the Bone Fracture Risk Model

Introduction: The probability of bone fracture during and after spaceflight is quantified to aid in mission planning, to determine required astronaut fitness standards and training requirements and to inform countermeasure research and design. Probability is quantified with a probabilistic modeling approach where distributions of model parameter values, instead of single deterministic values, capture the parameter variability within the astronaut population and fracture predictions are probability distributions with a mean value and an associated uncertainty. Because of this uncertainty, the model in its current state cannot discern an effect of countermeasures on fracture probability, for example between use and non-use of bisphosphonates or between spaceflight exercise performed with the Advanced Resistive Exercise Device (ARED) or on devices prior to installation of ARED on the International Space Station. This is thought to be due to the inability to measure key contributors to bone strength, for example, geometry and volumetric distributions of bone mass, with areal bone mineral density (BMD) measurement techniques. To further the applicability of model, we performed a parameter sensitivity study aimed at identifying those parameter uncertainties that most effect the model forecasts in order to determine what areas of the model needed enhancements for reducing uncertainty. Methods: The bone fracture risk model (BFxRM), originally published in (Nelson et al) is a probabilistic model that can assess the risk of astronaut bone fracture. This is accomplished by utilizing biomechanical models to assess the applied loads; utilizing models of spaceflight BMD loss in at-risk skeletal locations; quantifying bone strength through a relationship between areal BMD and bone failure load; and relating fracture risk index (FRI), the ratio of applied load to bone strength, to fracture probability. There are many factors associated with these calculations including environmental factors, factors associated with the fall event, mass and anthropometric values of the astronaut, BMD characteristics, characteristics of the relationship between BMD and bone strength and bone fracture characteristics. The uncertainty in these factors is captured through the use of parameter distributions and the fracture predictions are probability distributions with a mean value and an associated uncertainty. To determine parameter sensitivity, a correlation coefficient is found between the sample set of each model parameter and the calculated fracture probabilities. Each parameters contribution to the variance is found by squaring the correlation coefficients, dividing by the sum of the squared correlation coefficients, and multiplying by 100. Results: Sensitivity analyses of BFxRM simulations of preflight, 0 days post-flight and 365 days post-flight falls onto the hip revealed a subset of the twelve factors within the model which cause the most variation in the fracture predictions. These factors include the spring constant used in the hip biomechanical model, the midpoint FRI parameter within the equation used to convert FRI to fracture probability and preflight BMD values. Future work: Plans are underway to update the BFxRM by incorporating bone strength information from finite element models (FEM) into the bone strength portion of the BFxRM. Also, FEM bone strength information along with fracture outcome data will be incorporated into the FRI to fracture probability.

mathematical models↗
NASA NTRSDate not supplied

Multiscale Modeling of Fracture Strength in Fibrous Thermal Protection System Materials

This work presents a multiscale modeling approach to predict the fracture strength of fibrous Thermal Protection System (TPS) materials. The model assumes that system failure is initiated at the joints between individual fibers. We investigated three distinct TPS compositions: amorphous silica, alumina and aluminosilicate fibers. Molecular dynamics (MD) simulations were employed to determine the fracture strength values of these fiber joints for both material systems. These fracture strength values were then integrated into simulations of 3D randomly populated fiber structures, where tensile load transfer occurs through the fiber joints. These microscale properties are upscaled through a renormalization approach [1] to predict macroscale tensile strength of 3D random fiber networks, accounting for joint-dominated failure and effective load-bearing area. The study concludes by demonstrating the resulting strength variation as a function of material composition, fiber density, and morphology. We also show validation of results by comparing them against explicit fiber finite element (FE) modeling [2] where fiber joint fracture is represented by cohesive elements.

Jaehyun Cho↗
NASA NTRS2006

Multiscale Modeling of Grain-Boundary Fracture: Cohesive Zone Models Parameterized From Atomistic Simulations

A multiscale modeling strategy is developed to study grain boundary fracture in polycrystalline aluminum. Atomistic simulation is used to model fundamental nanoscale deformation and fracture mechanisms and to develop a constitutive relationship for separation along a grain boundary interface. The nanoscale constitutive relationship is then parameterized within a cohesive zone model to represent variations in grain boundary properties. These variations arise from the presence of vacancies, intersticies, and other defects in addition to deviations in grain boundary angle from the baseline configuration considered in the molecular dynamics simulation. The parameterized cohesive zone models are then used to model grain boundaries within finite element analyses of aluminum polycrystals.

Glaessgen, Edward H.↗
NASA NTRSDate not supplied

Macro-scale Testing and Micromechanics Modeling of Fracture Behaviors for Boron Carbide Composites with Hierarchical Microstructures

With a unique combination of properties including high hardness, low density, wear and corrosion resistance, thermal stability, high neutron absorption, and semi-conductivity, boron carbide (B4C) is a candidate material for various engineering applications that involve extreme environment. The current applications of boron carbide, however, are limited by its intrinsic brittleness due to its strong covalent bonding. To toughen boron carbide, in this work hierarchical microstructure designs was used to provide multiple toughening mechanisms including crack deflection/bridging, micro-crack toughening, etc. Using field assisted sintering, B4C composites with hierarchical microstructure features including graphite platelets, micron and sub-micron sized TiB2 reinforcements were fabricated. The fracture toughness of fabricated B4C composites were previously measured at micro-scale using micro-indentation followed by post-testing microstructure inspection. However, questions including whether the fracture toughness enhancement measured at micro-scale can translate to macro-scale mechanical properties, and what are the fundamental mechanisms behind observed fracture toughness enhancement, remain to be answered. In this study, the fabricated B4C composites were tested using standardized four-point bending method to obtain fracture toughness at macro-scale. In addition, micromechanics modeling was conducted using MAC/GMC code and crack-band model to study the effect of residual stress and weak interphases on fracture behaviors of B4C composites reinforced with TiB2 particles. Through standardized four-point bending tests, fracture toughness enhancements up to 2.85, 3.32, and 3.65 MPa∙m(exp 1/2) (from 2.38 MPa∙m(exp 1/2)) were achieved for B4C composites with graphite platelets addition (micro/nano B4C), with TiB2 formation (micro B4C-TiB2), and with both graphite and TiB2 addition (micro/nano B4C-TiB2) respectively. Micromechanics modeling indicated that introduction of thermal residual stress and weak interphases caused enhanced micro-cracking behavior and resulted in the observed fracture toughness enhancement. These results furthered understanding about the mechanical behaviors at macro-scale and the mechanisms behind observed fracture toughness enhancement for B4C composites with hierarchical microstructures and can provide reference data for the future design of B4C composites with optimized microstructures for further fracture toughness enhancement.

Boron Carbide↗
NASA NTRS2002

Micromechanics Modeling of Fracture in Nanocrystalline Metals

Nanocrystalline metals have very high theoretical strength, but suffer from a lack of ductility and toughness. Therefore, it is critical to understand the mechanisms of deformation and fracture of these materials before their full potential can be achieved. Because classical fracture mechanics is based on the comparison of computed fracture parameters, such as stress intlmsity factors, to their empirically determined critical values, it does not adequately describe the fundamental physics of fracture required to predict the behavior of nanocrystalline metals. Thus, micromechanics-based techniques must be considered to quanti@ the physical processes of deformation and fracture within nanocrystalline metals. This paper discusses hndamental physicsbased modeling strategies that may be useful for the prediction Iof deformation, crack formation and crack growth within nanocrystalline metals.

Glaessgen, E. H.↗
NASA NTRS2002

Prediction of Long-Term Strength of Thermoplastic Composites Using Time-Temperature Superposition

Accelerated tests for composite failure were investigated. Constant ramp transverse strength tests on thermoplastic composite specimens were conducted at four temperatures from 300 F to 450 F and five duration times from 0.5 sec to 24 hrs. Up to 400 F, the time-temperature-superposition method produces a master curve allowing strength at longer times to be estimated from strength tests conducted over shorter times but at higher temperatures. The shift factors derived from compliance tests applied well to the strength data. To explain why strength behaved similar to compliance, a viscoelastic fracture model was investigated based on the hypothesis that the work of fracture for crack initiation at some critical flaw remains constant with time and temperature. The model, which used compliance as input, was found to fit the strength data only if the critical fracture energy was allowed to vary with stress rate. Fracture tests using double cantilever beam specimens were conducted from 300 F to 450 F over time scales similar to the strength study. The toughness data showed a significant change with loading rate, less variation with temperature, did not form a master curve, and could not be correlated with the fracture model. Since the fracture model did not fit the fracture data, an alternative explanation based on the dilatational strain energy density was proposed. However the usefulness of this model is severely limited because it relies on a critical parameter which varies with loading rate.

Reeder, James R.↗
NASA NTRS1992

Floor-fractured crater models for igneous crater modification on Venus

Although crater modification on the Earth, Moon, and Mars results from surface erosion and crater infilling, a significant number of craters on the Moon also exhibit distinctive patterns of crater-centered fracturing and volcanism that can be modeled as the result of igneous crater modification. Here, we consider the possible effects of Venus surface conditions on this model, describe two examples of such crater modification, and then briefly discuss the constraints these craters place on conditions at depth.

Wichman, R. W.↗
NASA NTRS1982

Factors influencing the thermally-induced strength degradation of B/Al composites

Literature data related to the thermally-induced strength degradation of B/Al composites were examined in the light of fracture theories based on reaction-controlled fiber weakening. Under the assumption of a parabolic time-dependent growth for the interfacial reaction product, a Griffith-type fracture model was found to yield simple equations whose predictions were in good agreement with data for boron fiber average strength and for B/Al axial fracture strain. The only variables in these equations were the time and temperature of the thermal exposure and an empirical factor related to fiber surface smoothness prior to composite consolidation. Such variables as fiber diameter and aluminum alloy composition were found to have little influence. The basic and practical implications of the fracture model equations are discussed.

Dicarlo, J. A.↗
NASA NTRS1983

Factors influencing the thermally-induced strength degradation of B/Al composites

Literature data related to the thermally-induced strength degradation of B/Al composites were examined in the light of fracture theories based on reaction-controlled fiber weakening. Under the assumption of a parabolic time-dependent growth for the interfacial reaction product, a Griffith-type fracture model was found to yield simple equations whose predictions were in good agreement with data for boron fiber average strength and for B/Al axial fracture strain. The only variables in these equations were the time and temperature of the thermal exposure and an empirical factor related to fiber surface smoothness prior to composite consolidation. Such variables as fiber diameter and aluminum alloy composition were found to have little influence. The basic and practical implications of the fracture model equations are discussed. Previously announced in STAR as N82-24297

Dicarlo, J. A.↗
NASA NTRS1992

Floor-fractured crater models of the Sudbury structure, Canada

The Sudbury structure in Ontario, Canada, is one of the oldest and largest impact structures recognized in the geological record. It is also one of the most extensively deformed and volcanically modified impact structures on Earth. Although few other terrestrial craters are recognized as volcanically modified, numerous impact craters on the Moon have been volcanically and tectonically modified and provide possible analogs for the observed pattern of modification at Sudbury. We correlate the pattern of early deformation at Sudbury to fracture patterns in two alternative lunar analogs and then use these analogs both to estimate the initial size of the Sudbury structure and to model the nature of early crater modification at Sudbury.

Wichman, R. W.↗
NASA NTRS2011

Statistical Models of Fracture Relevant to Nuclear-Grade Graphite: Review and Recommendations

The nuclear-grade (low-impurity) graphite needed for the fuel element and moderator material for next-generation (Gen IV) reactors displays large scatter in strength and a nonlinear stress-strain response from damage accumulation. This response can be characterized as quasi-brittle. In this expanded review, relevant statistical failure models for various brittle and quasi-brittle material systems are discussed with regard to strength distribution, size effect, multiaxial strength, and damage accumulation. This includes descriptions of the Weibull, Batdorf, and Burchell models as well as models that describe the strength response of composite materials, which involves distributed damage. Results from lattice simulations are included for a physics-based description of material breakdown. Consideration is given to the predicted transition between brittle and quasi-brittle damage behavior versus the density of damage (level of disorder) within the material system. The literature indicates that weakest-link-based failure modeling approaches appear to be reasonably robust in that they can be applied to materials that display distributed damage, provided that the level of disorder in the material is not too large. The Weibull distribution is argued to be the most appropriate statistical distribution to model the stochastic-strength response of graphite.

Nemeth, Noel N.↗
NASA NTRS1985

Notched strength of composite laminates: Predictions and experiments - A review

A self-contained review of several semiempirical fracture models for predicting notched strength of composite laminates is presented, based on notched strength data on 70 different laminate configurations of graphite/epoxy, boron/aluminum, and graphite/polyimide. Emphasis is placed on experimental results concerning such failure factors as delamination, splitting, and size of damage zone. Moreover, the fracture model parameters are correlated with the notch sensitivity of composite laminates, and the applicability of the correlations in describing the material notch sensitivity is evaluated. The predictions provided by the different models were found to be identical for all practical purposes.

Awerbuch, J.↗