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 415 records · Page 23

A stochastic damage model for the rupture prediction of a multi-phase solid. I - Parametric studies. II - Statistical approach

A stochastic damage model for predicting the rupture of a brittle multiphase material is developed, based on the microcrack-macrocrack interaction. The model, which incorporates uncertainties in locations, orientations, and numbers of microcracks, characterizes damage by microcracking and fracture by macrocracking. A parametric study is carried out to investigate the change of the stress intensity at the macrocrack tip by the configuration of microcracks. The inherent statistical distribution of the fracture toughness arising from the intrinsic random nature of microcracks is explored using a statistical approach. For this purpose, a computer simulation model is introduced, which incorporates a statistical characterization of geometrical parameters of a random microcrack array.

Lua, Yuan J.↗

Continuum modeling of catastrophic collisions

A two dimensional hydrocode based on 2-D SALE was modified to include strength effects and fragmentation equations for fracture resulting from tensile stress in one dimension. Output from this code includes a complete fragmentation summary for each cell of the modeled object: fragment size (mass) distribution, vector velocities of particles, peak values of pressure and tensile stress, and peak strain rates associated with fragmentation. Contour plots showing pressure and temperature at given times within the object are also produced. By invoking axial symmetry, three dimensional events can be modeled such as zero impact parameter collisions between asteroids. The code was tested against the one dimensional model and the analytical solution for a linearly increasing tensile stress under constant strain rate.

Ryan, Eileen V.↗

Strain energy release rates of composite interlaminar end-notch and mixed-mode fracture: A sublaminate/ply level analysis and a computer code

A computer code is presented for the sublaminate/ply level analysis of composite structures. This code is useful for obtaining stresses in regions affected by delaminations, transverse cracks, and discontinuities related to inherent fabrication anomalies, geometric configurations, and loading conditions. Particular attention is focussed on those layers or groups of layers (sublaminates) which are immediately affected by the inherent flaws. These layers are analyzed as homogeneous bodies in equilibrium and in isolation from the rest of the laminate. The theoretical model used to analyze the individual layers allows the relevant stresses and displacements near discontinuities to be represented in the form of pure exponential-decay-type functions which are selected to eliminate the exponential-precision-related difficulties in sublaminate/ply level analysis. Thus, sublaminate analysis can be conducted without any restriction on the maximum number of layers, delaminations, transverse cracks, or other types of discontinuities. In conjunction with the strain energy release rate (SERR) concept and composite micromechanics, this computational procedure is used to model select cases of end-notch and mixed-mode fracture specimens. The computed stresses are in good agreement with those from a three-dimensional finite element analysis. Also, SERRs compare well with limited available experimental data.

Valisetty, R. R.↗

Small Scale Polygons and the History of Ground Ice on Mars

Recent progress on polygon modeling has focused on the diameter and surface relief that we expect of thermal-contraction polygons in martian permafrost. With this in mind, we developed a finite-element model of thermal-contraction-crack behavior in permafrost in a martian climate. This model was generated from a finite element code by Jay Melosh (called TECTON) originally developed for terrestrial and planetary crustal-deformation studies. We adapted this model to martian permafrost by including time (and temperature) dependent rheologies, boundary conditions, and isotropic thermal-contraction, as well as several small adaptations to a martian environment. We tested our model extensively, including comparison to an analytic solution of pre-fracture stress. We recently published an analysis of two potential sources of water for forming the recent gullies. In this work we first evaluated the potential for near-surface ground ice (in the top meter or so of soil) to melt under conditions of solar heating on sloped surfaces at high obliquity, utilizing both thermal and diffusion-based ground-ice-stability models; our results suggested that the ground ice will sublimate, and the ice table will recede to greater depths before the melting temperature can be reached. An exception can occur only for extremely salt-rich ice, depressing the freezing point.

Mellon, Michael T.↗

What is (and Isn't) Wrong with Both the Tension and Shear Failure Models for the Formation of Lineae on Europa

An unresolved problem in the interpretation of lineae on Europa is whether they formed as tension- or shear-fractures. Voyager image analyses led to hypotheses that Europan lineaments are tension cracks induced by tidal deformation of the ice crust. This interpretation continued with Galileo image analyses, with lineae being classified as crust- penetrating tension cracks. Tension fracturing has also been an implicit assumption of nonsynchronous rotation (NSR) studies. However, recent hypotheses invoke shear failure to explain lineae development. If a shear failure mechanism is correct, it will be necessary to re-evaluate any models for the evolution of Europa's crust that are based on tensile failure models, such as NSR estimates. For this reason, it is imperative that the mechanism by which fractures are initiated on Europa be unambiguously unraveled. A logical starting point is an evaluation of the pros and cons of each failure model, highlighting the lines of evidence that are needed to fully justify either model.

Kattenhorn, S. A.↗

Numerical Modeling of Shatter Cones Development in Impact Craters

Shatter cones are the characteristic forms of rock fractures in impact structures. They have been used for decades as unequivocal fingerprints of meteoritic impacts on Earth. The abundant data about shapes, apical angles, sizes and distributions of shatter cones for many terrestrial impact structures should provide insights for the determination of impact conditions and characteristics of shock waves produced by high-velocity projectiles in geologic media. However, previously proposed models for the formation of shatter cones do not agree with observations. For example, the widely accepted Johnson-Talbot mechanism requires that the longitudinal stress drops to zero between the arrival of the elastic precursor and the main plastic wave. Unfortunately, observations do not support such a drop. A model has been also proposed to explain the striated features on the surface of shatter cones but can not invoked for their conical shape. The mechanism by which shatter cones form thus remains enigmatic to date. In this paper we present a new model for the formation of shatter cones. Our model has been tested by means of numerical simulations using the hydrocodes SALE 2D enhanced with the Grady-Kipp-Melosh fragmentation model.

Baratoux, D.↗

Probabilistic Structural Analysis Methods (PSAM) for select space propulsion system components, part 2

The technical effort and computer code enhancements performed during the sixth year of the Probabilistic Structural Analysis Methods program are summarized. Various capabilities are described to probabilistically combine structural response and structural resistance to compute component reliability. A library of structural resistance models is implemented in the Numerical Evaluations of Stochastic Structures Under Stress (NESSUS) code that included fatigue, fracture, creep, multi-factor interaction, and other important effects. In addition, a user interface was developed for user-defined resistance models. An accurate and efficient reliability method was developed and was successfully implemented in the NESSUS code to compute component reliability based on user-selected response and resistance models. A risk module was developed to compute component risk with respect to cost, performance, or user-defined criteria. The new component risk assessment capabilities were validated and demonstrated using several examples. Various supporting methodologies were also developed in support of component risk assessment.

Source record↗

The Digital Astronaut Project Computational Bone Remodeling Model (Beta Version) Bone Summit Summary Report

Under the conditions of microgravity, astronauts lose bone mass at a rate of 1% to 2% a month, particularly in the lower extremities such as the proximal femur [1-3]. The most commonly used countermeasure against bone loss in microgravity has been prescribed exercise [4]. However, data has shown that existing exercise countermeasures are not as effective as desired for preventing bone loss in long duration, 4 to 6 months, spaceflight [1,3,5,6]. This spaceflight related bone loss may cause early onset of osteoporosis to place the astronauts at greater risk of fracture later in their lives. Consequently, NASA seeks to have improved understanding of the mechanisms of bone demineralization in microgravity in order to appropriately quantify this risk, and to establish appropriate countermeasures [7]. In this light, NASA's Digital Astronaut Project (DAP) is working with the NASA Bone Discipline Lead to implement well-validated computational models to help predict and assess bone loss during spaceflight, and enhance exercise countermeasure development. More specifically, computational modeling is proposed as a way to augment bone research and exercise countermeasure development to target weight-bearing skeletal sites that are most susceptible to bone loss in microgravity, and thus at higher risk for fracture. Given that hip fractures can be debilitating, the initial model development focused on the femoral neck. Future efforts will focus on including other key load bearing bone sites such as the greater trochanter, lower lumbar, proximal femur and calcaneus. The DAP has currently established an initial model (Beta Version) of bone loss due to skeletal unloading in femoral neck region. The model calculates changes in mineralized volume fraction of bone in this segment and relates it to changes in bone mineral density (vBMD) measured by Quantitative Computed Tomography (QCT). The model is governed by equations describing changes in bone volume fraction (BVF), and rates of changes in bone cell populations that remove and replace bone in packets within the bone region. The DAP bone model is unique in several respects. In particular in takes former models of volume fraction changes one step higher in fidelity and separates BVF into separate equations for mineralized and osteoid volume fractions governed by a mineralization rate. This more closely follows the physiology of the remodeling unit cycles where bone is first resorbed and then followed by the action of osteoblasts to lay down collagen matrix which eventually becomes mineralized. In another respect, the modules allow the functional description of the time rate of change of other parameters and variables in the model during a computational simulation. More detailed description of the model, preliminary validation results, current limitation and caveats, and planned advancements are provided in sections 2 through 5. The DAP bone model is being developed primarily as a research tool, and not as a clinical tool like QCT. Even if it transitions to a clinical tool, it is not intended to replace QCT or any other clinical tool. Moreover, the DAP bone model does not predict bone fracture. Its purpose is to provide valuable additional data via "forward prediction" simulations for during and after spaceflight missions to gain insight on, (1) mechanisms of bone demineralization in microgravity, and (2) the volumetric changes at the various bone sites in response to in-flight and post-flight exercise countermeasures. This data can then be used as input to the Keyak [8] (or equivalent) FE analysis method to gain insight on how bone strength may change during and after flight. This information can also be useful to help optimize exercise countermeasure protocols to minimize changes in bone strength during flight, and improve regain of bone strength post-flight. To achieve this goal, the bone model will be integrated with DAP's exercise countermeasure models to simulate the effect of exercise prescriptions on preserving bone. More specifically, the model will accept loading history due to muscle and joint force on bone and produce quantified remodeling within the bone region under influence of the applied stress. Furthermore, because they tend to respond differently, the bone remodeling model includes both trabecular bone and cortical bone.

Pennline, James↗

Fundamental analysis of the failure of polymer-based fiber reinforced composites

A mathematical model predicting the strength of unidirectional fiber reinforced composites containing known flaws and with linear elastic-brittle material behavior was developed. The approach was to imbed a local heterogeneous region surrounding the crack tip into an anisotropic elastic continuum. This (1) permits an explicit analysis of the micromechanical processes involved in the fracture, and (2) remains simple enough to be useful in practical computations. Computations for arbitrary flaw size and orientation under arbitrary applied loads were performed. The mechanical properties were those of graphite epoxy. With the rupture properties arbitrarily varied to test the capabilities of the model to reflect real fracture modes, it was shown that fiber breakage, matrix crazing, crack bridging, matrix-fiber debonding, and axial splitting can all occur during a period of (gradually) increasing load prior to catastrophic failure. The calculations also reveal the sequential nature of the stable crack growth process proceding fracture.

Kanninen, M. F.↗

Microscale Fracture Testing of HfO2-Si Environmental Barrier Coatings

Environmental barrier coatings (EBCs) play a vital role in protecting advanced turbine components constructed from ceramic matrix composites. EBC performance is evaluated based on many factors such as mechanical and thermal properties and chemical stability. This study focuses on the mechanical properties of individual phases within a HfO2-Si EBC through the combined use of chevron-notched microcantilevers and finite element modeling. The Si- and HfO2-rich phases were found to have fracture toughness K(sub IC) values of 1.36+/-0.07 MPa·m1/2 and 1.99+/-0.36 MPa·m1/2, respectively. Despite irregular notch geometry, calculated K(sub IC) values agreed reasonably with literature values of monolithic samples for the respective phases. The relatively large scatter in the HfO2-rich samples is attributed to heterogeneity in fracture feature size and distribution of the HfO2 and HfSiO4 phases. Comparisons with bulk single-edge V-notched beams suggest bulk behavior is dictated by the fracture behavior of the Sirich phase with no contribution from the HfO2-rich phase.

Robertson, Andrew L.↗

Fracture behavior of unidirectional boron/aluminum composite laminates

An experiment was conducted to verify the results of mathematical models which predict the stresses and displacements of fibers and the amount of damage growth in a center-notched lamina as a function of the applied remote stress and the matrix and fiber material properties. A brittle lacquer coating was used to detect the yielding in the matrix while X-ray techniques were used to determine the number of broken fibers in the laminate. The notched strengths and the amounts of damage found in the specimens agree well with those predicted by the mathematical model. It is shown that the amount of damage and the crack opening displacement does not depend strongly on the number of plies in the laminate for a given notch width. By heat-treating certain laminates to increase the yield stress of the alumina matrix, the effect of different matrix properties on the fracture behavior was investigated. The stronger matrix is shown to weaken the notched laminate by decreasing the amount of matrix damage, thereby making the laminate more notch sensitive.

Goree, J. G.↗

Differential continuum damage mechanics models for creep and fatigue of unidirectional metal matrix composites

Three multiaxial isothermal continuum damage mechanics models for creep, fatigue, and creep/fatigue interaction of a unidirectional metal matrix composite volume element are presented, only one of which will be discussed in depth. Each model is phenomenological and stress based, with varying degrees of complexity to accurately predict the initiation and propagation of intergranular and transgranular defects over a wide range of loading conditions. The development of these models is founded on the definition of an initially transversely isotropic fatigue limit surface, static fracture surface, normalized stress amplitude function and isochronous creep damage failure surface, from which both fatigue and creep damage evolutionary laws can be obtained. The anisotropy of each model is defined through physically meaningful invariants reflecting the local stress and material orientation. All three transversely isotropic models have been shown, when taken to their isotropic limit, to directly simplify to previously developed and validated creep and fatigue continuum damage theories. Results of a nondimensional parametric study illustrate (1) the flexibility of the present formulation when attempting to characterize a large class of composite materials, and (2) its ability to predict anticipated qualitative trends in the fatigue behavior of unidirectional metal matrix composites. Additionally, the potential for the inclusion of various micromechanical effects (e.g., fiber/matrix bond strength, fiber volume fraction, etc.), into the phenomenological anisotropic parameters is noted, as well as a detailed discussion regarding the necessary exploratory and characterization experiments needed to utilize the featured damage theories.

Arnold, S. M.↗

Metal matrix composites microfracture: Computational simulation

Fiber/matrix fracture and fiber-matrix interface debonding in a metal matrix composite (MMC) are computationally simulated. These simulations are part of a research activity to develop computational methods for microfracture, microfracture propagation and fracture toughness of the metal matrix composites. The three-dimensional finite element model used in the simulation consists of a group of nine unidirectional fibers in three by three unit cell array of SiC/Ti15 metal matrix composite with a fiber volume ration of 0.35. This computational procedure is used to predict the fracture process and establish the hierarchy of fracture modes based on strain energy release rate. It is also used to predict stress redistribution to surrounding matrix-fibers due to initial and progressive fracture of fiber/matrix and due to debonding of fiber-matrix interface. Microfracture results for various loading cases such as longitudinal, transverse, shear and bending are presented and discussed. Step-by-step procedures are outlined to evaluate composite microfracture for a given composite system.

Mital, Subodh K.↗

Metal matrix composites microfracture - Computational simulation

Fiber/matrix fracture and fiber-matrix interface debonding in a metal matrix composite (MMC) are computationally simulated. These simulations are part of a research activity to develop computational methods for microfracture, microfracture propagation and fracture toughness of the metal matrix composites. The three-dimensional finite element model used in the simulation consists of a group of nine unidirectional fibers in three by three unit cell array of SiC/Ti15 metal matrix composite with a fiber volume ration of 0.35. This computational procedure is used to predict the fracture process and establish the hierarchy of fracture modes based on strain energy release rate. It is also used to predict stress redistribution to surrounding matrix-fibers due to initial and progressive fracture of fiber/matrix and due to debonding of fiber-matrix interface. Microfracture results for various loading cases such as longitudinal, transverse, shear and bending are presented and discussed. Step-by-step procedures are outlined to evaluate composite microfracture for a given composite system.

Mital, S. K.↗

Design and construction of models for the National Transonic Facility, part 1

The design and construction of models for the National Transonic Facility (NTF) has resulted in significant technology developments in many areas. This lecture covers the development of design criteria and major research and development work that has contributed to the successful design and fabrication models for testing at full scale Reynolds number the NTF. Emphasis is placed on the materials aspect of the design and fabrication proces, including metallic materials, mechanical properties characterization, new steel alloy development, fracture toughness enhancement, and identification of fillers and solders suitable for use in cryogenic models. Quantitative data are provided which will be of value to the potential user of NTF or for application to the design and fabrication of model systems for other cryogenic wind tunnels.

Young, C. P., Jr.↗

Impact Response of Thermoplastic Composites – Experiments and Modeling

This study investigates the impact response of thermoplastic composite materials. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK reinforced with T700G (T700/LMPAEK ). Several High Energy Dynamic Impact (HEDI) tests were performed wherein an aluminum projectile impacted T700/LMPAEK panels at different velocities. In the impact experiment studied herein, the projectile impacted the T700/LMAPEK panel at 73.2 m/s and rebounded. The ensuing damage to the T700/LMPAEK panel was characterized via ultrasonic C-scans, visualizing damage modes including matrix cracking and delaminations. Since the projectile rebounded, the fiber breakage damage mode was not evident. The impact event was modeled using a combined continuum damage mechanics and cohesive zone analysis approach with the commercially available hydrocode LS-DYNA®. The dynamic deformation and damage of the continuum, i.e., the plies in the T700/LMPAEK laminate, was modeled using a Deformation Gradient Decomposition (DGD) material model, i.e., MAT299 available within LS-DYNA. The interlaminar regions were treated as zero-thickness cohesive zones, wherein the dynamic delamination mechanics was modeled using a traction-separation formulation. Mixed-mode delamination growth was captured using the Benzegaggh-Kenane (B-K) traction-separation interpolation law. Strain-rate sensitivity of interlaminar fracture toughness was included using approximations based on published literature. The outcomes of the modeling effort are discussed herein which demonstrate good correlations with test data regarding panel deflection, panel-projectile interactions, damage modes and damage extent.

Impact modeling↗

On deformation and strength.

Mechanical strength of medium formulated as microscopic behavior in continuum obtaining time- dependent macroscopic fracture strength

TENSILE STRENGTH↗

A study of the Tyrone-Mount Union lineament by remote sensing techniques and field methods

The author has identified the following significant results. This study has shown that subtle variations in fold axes, fold form, and stratigraphic thickness can be delineated. Many of the conclusions were based on extrapolation in similitude to different scales. A conceptual model was derived for the Tyrone-Mount Union lineament. In this model, the lineament the morphological expression of a zone of fracture concentrations which penetrated basement rocks and may have acted as a curtain to regional stresses or as a domain boundary between uncoupled adjacent crustal blocks.

Gold, D. P.↗