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 271 records · Page 15

A Fragment-Cloud Model for Breakup of Asteroids with Varied Internal Structures

As an asteroid descends toward Earth, it deposits energy in the atmosphere through aerodynamic drag and ablation. Asteroid impact risk assessments rely on energy deposition estimates to predict blast overpressures and ground damage that may result from an airburst, such as the one that occurred over Chelyabinsk, Russia in 2013. The rates and altitudes at which energy is deposited along the entry trajectory depend upon how the bolide fragments, which in turn depends upon its internal structure and composition. In this work, we have developed an analytic asteroid fragmentation model to assess the atmospheric energy deposition of asteroids with a range of structures and compositions. The modeling approach combines successive fragmentation of larger independent pieces with aggregate debris clouds released with each fragmentation event. The model can vary the number and masses of fragments produced, the amount of mass released as debris clouds, the size-strength scaling used to increase the robustness of smaller fragments, and other parameters. The initial asteroid body can be seeded with a distribution of independent fragment sizes amid a remaining debris mass to represent loose rubble pile conglomerations, can be given an outer regolith later, or can be defined as a coherent or fractured monolith. This approach enables the model to represent a range of breakup behaviors and reproduce detailed energy deposition features such as multiple flares due to successive burst events, high-altitude regolith blow-off, or initial disruption of rubble piles followed by more energetic breakup of the constituent boulders. These capabilities provide a means to investigate sensitivities of ground damage to potential variations in asteroid structure.

Asteroid↗

Three-Dimensional Gear Crack Propagation Studied

Gears used in current helicopters and turboprops are designed for light weight, high margins of safety, and high reliability. However, unexpected gear failures may occur even with adequate tooth design. To design an extremely safe system, the designer must ask and address the question, "What happens when a failure occurs?" With gear-tooth bending fatigue, tooth or rim fractures may occur. A crack that propagates through a rim will be catastrophic, leading to disengagement of the rotor or propeller, loss of an aircraft, and possible fatalities. This failure mode should be avoided. A crack that propagates through a tooth may or may not be catastrophic, depending on the design and operating conditions. Also, early warning of this failure mode may be possible because of advances in modern diagnostic systems. One concept proposed to address bending fatigue fracture from a safety aspect is a splittooth gear design. The prime objective of this design would be to control crack propagation in a desired direction such that at least half of the tooth would remain operational should a bending failure occur. A study at the NASA Lewis Research Center analytically validated the crack-propagation failsafe characteristics of a split-tooth gear. It used a specially developed three-dimensional crack analysis program that was based on boundary element modeling and principles of linear elastic fracture mechanics. Crack shapes as well as the crack-propagation life were predicted on the basis of the calculated stress intensity factors, mixed-mode crack-propagation trajectory theories, and fatigue crack-growth theories. The preceding figures show the effect of the location of initial cracks on crack propagation. Initial cracks in the fillet of the teeth produced stress intensity factors of greater magnitude (and thus, greater crack growth rates) than those in the root or groove areas of the teeth. Crack growth was simulated in a case study to evaluate crack-propagation paths. Tooth fracture was predicted from the crackgrowth simulation for an initial crack in the tooth fillet region. This was the desired failure mode for an ultrasafe design. Lastly, tooth loads on the uncracked mesh of the split-tooth design were up to five times greater than those on the cracked mesh if equal deflections of the cracked and uncracked teeth were considered. This effect needs to be considered in the design of a split-tooth configuration. This work was done in-house at Lewis in support of the National Rotorcraft Technology Center project, Ultra-Safe Gear Design, with the Boeing Defense and Space Group. The crack-propagation package was developed by the Cornell Fracture Group at Cornell University. The reported results, which are the initial findings of Lewis gear-crackpropagation research for the Rotorcraft Base Program, will be further investigated to develop generalized gear design guidelines.

Lewicki, David G.↗

R-Curve Instability Calculations Of Crack Growth

Report discusses use of instability method of calculation and R-curve mathematical models to analyze growth of cracks in fracture-mechanics specimens. In case of single material and structure, such analysis sometimes simple enough to be done on pocket calculator. Where microcomputer or larger computer available, comprehensive program includes libraries of driving-force equations for various configurations and R-curve mathematical models for different materials. Author concludes instability method simple and effective and model equations studied all viable in sense at lease one of them should fit almost any applicable set of crack-growth data. Method and models constitute powerful mathematical tools for analysis of fractures.

Orange, Thomas W.↗

Some fundamental aspects of the fatigue and fracture response of composite materials

A simple analytical model, based on the concept that the constraining plies surrounding cracks through the thickness of a given ply or plies control the spacing of those cracks by controlling the surrounding stress state, appears to agree well with experimental data. Such a model can be used to evaluate critically various constraint situations including the prediction of the ability of a given constraint layer to arrest flaw development. The analytical model, which uses an equilibrium element approach in the present case, is sensitive to lamina and laminate elastic properties, flaw size, stacking sequence, and an interlaminar stress transfer distance which is estimated from experimental observations. The analysis is easily extended to study flaw growth and plastic region formation.

Reifsnider, K. L.↗

Analysis of interlaminar fracture in composites under combined loading

Delamination is a predominant failure mode in continuous fiber reinforced laminated composite structures. One type of delamination is the transverse crack tip delamination which originates at the tip of transverse matrix cracks. An analytical model based on the sublaminate approach and fracture mechanics is developed to study the growth of such delaminations. Plane strain conditions are assumed and estimates are provided for the total strain energy release rate as well as the mode 1 and mode 2 contribution. The energy release rate estimates are used in combination with a simple failure law to predict critical delamination growth strains and stresses. These predictions are compared with experimental data on T300/934 Graphite Epoxy (+ or - 25/90 nano seconds) laminates in the range n=.5 to 8. A good agreement is demonstrated for the range of n where the experimental observations indicate transverse crack tip delamination to be the predominant failure mode.

Armanios, Erian A.↗

Fracture Mechanics Analysis of an Annular Crack in a Three-concentric-cylinder Composite Model

A boundary-value problem governing a three-phase concentric-cylinder model was analytically modeled to analyze annular interfacial crack problems with Love s strain functions in order to find the stress intensity factors (SIFs) and strain energy release rates (SERRs) at the tips of an interface crack in a nonhomogeneous medium. The complex form of a singular integral equation (SIE) of the second kind was formulated using Bessel s functions in the Fourier domain, and the SIF and total SERR were calculated using Jacoby polynomials. For the validity of the SIF equations to be established, the SIE of the three-concentric-cylinder model was reduced to the SIE for a two-concentric-cylinder model, and the results were compared with the previous results of Erdogan. A preliminary set of parametric studies was carried out to show the effect of interphase properties on the SERR. The method presented here provides insight about the effect of interphase properties on the crack driving force.

Kuguoglu, Latife H.↗

Prediction of Size Effects in Notched Laminates Using Continuum Damage Mechanics

This paper examines the use of a continuum damage model to predict strength and size effects in notched carbon-epoxy laminates. The effects of size and the development of a fracture process zone before final failure are identified in an experimental program. The continuum damage model is described and the resulting predictions of size effects are compared with alternative approaches: the point stress and the inherent flaw models, the Linear-Elastic Fracture Mechanics approach, and the strength of materials approach. The results indicate that the continuum damage model is the most accurate technique to predict size effects in composites. Furthermore, the continuum damage model does not require any calibration and it is applicable to general geometries and boundary conditions.

Camanho, D. P.↗

Skeletal Structural Consequences of Reduced Gravity Environments

The overall goal of this project is to provide structurally meaningful data on bone loss after exposure to reduced gravity environments so that more precise estimates of fracture risk and the effectiveness of countermeasures in reducing fracture risk can be developed. The project has three major components: (1) measure structural changes in the limb bones of rats subjected to complete and partial nonweightbearing, with and without treatment with ibandronate and periodic full weightbearing; (2) measure structural changes in the limb bones of human bedrest subjects, with and without treatment with alendronate and resistive exercise, and Russian cosmonauts flying on the Mir Space Station; and (3) validate and extend the 2-dimensional structural analyses currently possible in the second project component (bedrest and Mir subjects) using 3-dimensional finite element modeling techniques, and determine actual fracture-producing loads on earth and in space.

Ruff, Christropher B.↗

Computer Simulation of Fracture in Aerogels

Aerogels are of interest to the aerospace community primarily for their thermal properties, notably their low thermal conductivities. While the gels are typically fragile, recent advances in the application of conformal polymer layers to these gels has made them potentially useful as lightweight structural materials as well. In this work, we investigate the strength and fracture behavior of silica aerogels using a molecular statics-based computer simulation technique. The gels' structure is simulated via a Diffusion Limited Cluster Aggregation (DLCA) algorithm, which produces fractal structures representing experimentally observed aggregates of so-called secondary particles, themselves composed of amorphous silica primary particles an order of magnitude smaller. We have performed multi-length-scale simulations of fracture in silica aerogels, in which the interaction b e e n two secondary particles is assumed to be described by a Morse pair potential parameterized such that the potential range is much smaller than the secondary particle size. These Morse parameters are obtained by atomistic simulation of models of the experimentally-observed amorphous silica "bridges," with the fracture behavior of these bridges modeled via molecular statics using a Morse/Coulomb potential for silica. We consider the energetics of the fracture, and compare qualitative features of low-and high-density gel fracture.

Good, Brian S.↗

Establishment of redox conditions during planetary collisions as an origin of chondrites

Collisions between a 'cometlike' body (mixtures of chondritic materials and ice) and a slightly differentiated body were proposed for shock origin of ordinary chondrites. In this model, chondrules were formed with shock melting, and matrices were formed both by fracturing of materials and by recondensation of evaporated materials. This model can explain different redox conditions of chondrite formations by ice evaporation. Although this model was originally proposed for ordinary chondrites, we assume here that the model can be extended to chondrite formation in general. In this paper, redox conditions during chondrite formation by collisions will be discussed in the light of phase diagrams for solid-vapor equilibria.

Tsuchiyama, A.↗

A New Sea Ice Model: Pancakes to Metrics

We are developing, testing, and validating a new sea ice dynamics model that treats the ice cover as an elastic/decohesive material in the permanent pack and includes the correct frazil/pancake behavior in the marginal zone. Two salient features of present ice dynamics models are that they do not: 1) reproduce the oriented fracture patterns of openings and closings in the pack ice, and 2) accurately model the effects of frazil/pancake ice formation in the ice margin. These poorly modeled areas account for a substantial portion of the ice growth, turbulent heat flux to the atmosphere, salt flux to the ocean, and energy dissipation due to slippage, ridging, and rafting, in the Arctic. Existing sea ice models have shown limited success in predicting the degree to which a lead will open for prescribed or observed forcing conditions. An important aspect of the new model we are developing is that the existence of cracks, along with their orientation, opening, and closing, is predicted. To put this effort in perspective a short history of ice dynamics modeling and data collection is presented. The RGPS data set is used to validate the model. As part of the testing and validation of the model, we are working on a new metric for comparing linear features (leads and ridges) in the data and model to be used in data assimilation for this model. The model framework is presented as well as some results showing the creation and development of leads in a simulation of ice dynamics in the Beaufort Sea. Other presentations by the authors will show other results from this effort.

Coon, M.↗

Investigation of the effects of manufacturing variations and materials on fatigue crack detection methods in gear teeth

The fatigue life of a gear tooth can be thought of as the sum of the number of cycles required to initiate a crack, N(sub i), plus the number of cycles required to propagate the crack to such a length that fracture occurs, N(sub p). The factors that govern crack initiation are thought to be related to localized stress or strain at a point, while propagation of a fatigue crack is a function of the crack tip parameters such as crack shape, stress state, and stress intensity factor. During a test there is no clear transition between initiation and propagation. The mechanisms of initiation and propagation are quite different and modeling them separately produces a higher degree of accuracy, but then the question that continually arises is 'what is a crack?' The total life prediction in a fracture mechanics model presently hinges on the assumption of an initial crack length, and this length can significantly affect the total life prediction. The size of the initial crack is generally taken to be in the range of 0.01 in. to 0.2 in. Several researchers have used various techniques to determine the beginning of the crack propagation stage. Barhorst showed the relationship between dynamic stiffness changes and crack propagation. Acoustic emissions, which are stress waves produced by the sudden movement of stressed materials, have also been successfully used to monitor the growth of cracks in tensile and fatigue specimens. The purpose of this research is to determine whether acoustic emissions can be used to define the beginning of crack propagation in a gear using a single-tooth bending fatigue test.

Wheitner, Jeffrey A.↗

The phase composition of Triton's polar caps

Triton's polar caps are modeled as permanent nitrogen deposits hundreds of meters thick. Complex temperature variations on Triton's surface induce reversible transitions between the cubic and hexagonal phases of solid nitrogen, often with two coexisting propagating transition fronts. Subsurface temperature distributions are calculated using a two-dimensional thermal model with phase changes. The phase changes fracture the upper nitrogen layer, increasing its reflectivity and thus offering an explanation for the surprisingly high southern polar cap albedo (approximately 0.8) seen during the Voyager 2 flyby. The model has other implications for the phase transition phenomena on Triton, such as a plausible mechanism for the origin of geyser-like plume vent areas and a mechanism of energy transport toward them.

Duxbury, N. S.↗

Analytical and experimental investigation of the notched strength of thick laminates with surface notches

Experimental and analytical investigations of the fracture of a graphite/epoxy laminate with part-through semielliptic surface notches have been conducted. Experimental values of notched strength were determined for a variety of notch aspect ratios and notch depths. Analytical predictions of laminate failure were obtained from three different models. The three models were equally successful in predicting the fracture of laminates with deep notches. The experimental results of laminates with shallow notches indicated more notch sensitivity (that is greater strength reduction) than did the analytical models.

Harris, Charles E.↗

Microscale Constitutive Model Sensitivity on Multiscale Modeling of Fiber Reinforced Composites

Fiber reinforced composites are desirable in applications where low weight and high strength are needed, but are susceptible to microscale variability during manufacturing, making failure predictions difficult. The impact microscale variability has on macroscale mechanical response is difficult to predict due to the computational efficiency needed to simulate many, large, high fidelity, microscale models. In this study, a multiscale approach was taken to model 3-point bend, 4-point bend, and tensile experiments of a unidirectional composite from only having microstructure scans of these samples and constituent properties from literature. These scans were sampled with different sized windows, and statistically equivalent microstructures were generated, then simulated for stiffness, strength, and fracture toughness using an efficient micromechanical. Mesoscale models were created where element sizes equaled microstructure size, and properties were assigned through sampling of microscale simulation results. First, this study showed the effect of using Weibull scaling on constituent matrix strength on macroscale response. Then, a comparison was made between different element sizes and experiments. Finally, model dimensions were fixed, and the effect of randomly distributed local properties alone was examined. Results showed that the scatter of strength and stiffness in the experiments could be predicted well using images of the microscale fiber morphologies and that using stochastic properties produced a 3% coefficient of variation of strength for all experiments.

statistical microstructure↗

Microwave signatures and surface properties of Ovda Regio and surroundings, Venus

Magellan observations of horizontally and vertically polarized emissivity and radar specific cross sections are jointly modeled to separate dielectric constants from textural effect for lobate plains, fractured plains, fracture belts, festoon flows, and tessera in Ovda Regio and surroundings. The model assumes surface emission and scattering dominate and that signatures are controlled by a combination of quasi-specular and diffuse-scale mechanisms. Textural parameters are found to correlate with geology whereas dielectric constants are found to depend on elevation. A gradual increase in dielectric constant from rock-like values (4 to 5) to a value of approximately 50 is evident over the interval from 6054 to 6056 km, where elevations are derived from Magellan altimetry observations. Stereo radargrammetric analyses demonstrate that a return to rock-like values occurs over approximately 0.5 km above the terrain with highest dielectric constants. The abrupt return to rock-like values occurs independently of geology and reinforces the hypothesis of elevation control on dielectric constant values. Variations in dielectric constant with elevation are modeled using an 1000-ppm concentration of a ferroelectric mineral with a Curie temperature of 707 K. The model reproduces the abrupt change at highest elevations, as the Curie temperature is reached, and the gradual change as lower elevations are encountered. To our knowledge no other model explains the trends with such simplicity.

Arvidson, R. E.↗

Higher-Order Approximations for Stabilizing Zero-Energy Modes in Peridynamics Crystal Plasticity Models with Large Horizon Interactions

The non-ordinary state-based peridynamics theory combines non-local dynamic techniques with a desirable correspondence material principle, allowing for the use of continuum mechanics constitutive models. Such an approach presents a unique capability for solving problems involving discontinuities (e.g., strain localization, fracture, and fragmentation). However, the correspondence-based peridynamics models often suffer from zero-energy mode instabilities in numerical implementation, primarily due to the approximations of the non-local deformation gradient tensor. This paper focuses on a computational scheme for eliminating the zero-energy mode oscillations using a choice of influence functions that improve the truncation error in a higher-order Taylor series expansion of the deformation gradient. The novelty here is a tensor-based derivation of the linear constraint equations, which can be used to systematically identify the particle interaction weight functions for various user-specified horizon radii. In this paper, the proposed higher-order stabilization scheme is demonstrated for multi-dimensional examples involving polycrystalline and composite microstructures, along with comparisons against conventional finite element methods. The proposed stabilization scheme is shown to be highly effective in suppressing the spurious zero-energy mode oscillations in all numerical examples while enabling efficient simulations of strain localizations across material interfaces.

Non-Ordinary State-Based Peridynamics↗

A Computational Model for Simulating Spaceflight Induced Bone Remodeling

An overview of an initial development of a model of bone loss due to skeletal unloading in weight bearing sites is presented. The skeletal site chosen for the initial application of the model is the femoral neck region because hip fractures can be debilitating to the overall performance health of astronauts. The paper begins with the motivation for developing such a model of the time course of change in bone in order to understand the mechanism of bone demineralization experienced by astronauts in microgravity, to quantify the health risk, and to establish countermeasures. Following this, a general description of a mathematical formulation of the process of bone remodeling is discussed. Equations governing the rate of change of mineralized bone volume fraction and active osteoclast and osteoblast are illustrated. Some of the physiology of bone remodeling, the theory of how imbalance in remodeling can cause bone loss, and how the model attempts to capture this is discussed. The results of a preliminary validation analysis that was carried out are presented. The analysis compares a set of simulation results against bone loss data from control subjects who participated in two different bed rest studies. Finally, the paper concludes with outlining the current limitations and caveats of the model, and planned future work to enhance the state of the model.

Pennline, James A.↗