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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 883 records · Page 49

All sources NASA NTRS DOE OSTI

NASA NTRS2007

Ultrahigh-Temperature Ceramics

Ultrahigh temperature ceramics (UHTCs) are a class of materials that include the diborides of metals such as hafnium and zirconium. The materials are of interest to NASA for their potential utility as sharp leading edges for hypersonic vehicles. Such an application requires that the materials be capable of operating at temperatures, often in excess of 2,000 C. UHTCs are highly refractory and have high thermal conductivity, an advantage for this application. UHTCs are potentially applicable for other high-temperature processing applications, such as crucibles for molten-metal processing and high-temperature electrodes. UHTCs were first studied in the 1960 s by the U.S. Air Force. NASA s Ames Research Center concentrated on developing materials in the HfB2/SiC family for a leading-edge application. The work focused on developing a process to make uniform monolithic (2-phase) materials, and on the testing and design of these materials. Figure 1 shows arc-jet models made from UHTC materials fabricated at Ames. Figure 2 shows a cone being tested in the arc-jet. Other variations of these materials being investigated elsewhere include zirconium based materials and fiber-reinforced composites. Current UHTC work at Ames covers four broad topics: monoliths, coatings, composites, and processing. The goals include improving the fracture toughness, thermal conductivity and oxidation resistance of monolithic UHTCs and developing oxidation-resistant UHTC coatings for thermal-protection-system substrates through novel coating methods. As part of this effort, researchers are exploring compositions and processing changes that have yielded improvements in properties. Computational materials science and nanotechnology are being explored as approaches to reduce materials development time and improve and tailor properties.

Johnson, Sylvia M.↗

NASA NTRS1978

Martian lithospheric thickness from elastic flexure theory

The thickness of the elastic lithosphere in the Tharsis region of Mars is estimated from effects due to the surface load of Olympus Mons. Deformation (vertical displacement) and stress are calculated using elastic flexure theory for a range of possible lithospheric thicknesses (T), modeling the lithosphere as a thin elastic shell and the interior as a Newtonian fluid. For T below 150 km, displacement and stress rise rapidly with decreasing thickness. For T near 100 km, deformation of the region surrounding the volcano would be clearly visible in the topography, and resulting tensional stresses exceeding 5 kbar should produce observable fracturing at the surface. In contrast, for T near 200 km deformation is minimal and the tensional stress, being less than a kilobar, would not result in extensive fracturing. Since significant deformation and fracturing are not observed, it is concluded that the Martian elastic lithosphere is at least 150 km in thickness. Seismic, tectonic, and gravity observations all suggest a thick Martian lithosphere as well.

Thurber, C. H.↗

NASA NTRS1986

Episodic volcanism of tidally heated satellites with application to Io

Io is presently considered in light of a simple model for the coupled thermal and orbital evolution of a tidally heated satellite in an orbital resonance, demonstrating quantitatively how a feedback mechanism between the orbital and thermal energy of such a satellite may yield periodic surface heatflow and orbital eccentricity variations. The model predicts that the mean motion of Io may currently be increasing, as suggested by recent estimates of Io's mean motion on the basis of eclipse data. It is further inferred that the tidal stresses in the ice shell of Europa, whose eccentricity mimics that of Io, may recently have been sufficiently great to generate the fracturing observed.

Ojakangas, G. W.↗

NASA NTRS1991

Combined thermal and bending fatigue of high-temperature metal-matrix composites: Computational simulation

The nonlinear behavior of a high-temperature metal-matrix composite (HT-MMC) was simulated by using the metal matrix composite analyzer (METCAN) computer code. The simulation started with the fabrication process, proceeded to thermomechanical cyclic loading, and ended with the application of a monotonic load. Classical laminate theory and composite micromechanics and macromechanics are used in METCAN, along with a multifactor interaction model for the constituents behavior. The simulation of the stress-strain behavior from the macromechanical and the micromechanical points of view, as well as the initiation and final failure of the constituents and the plies in the composite, were examined in detail. It was shown that, when the fibers and the matrix were perfectly bonded, the fracture started in the matrix and then propagated with increasing load to the fibers. After the fibers fractured, the composite lost its capacity to carry additional load and fractured.

Gotsis, Pascal K.↗

NASA NTRS1990

The potential for ductility enhancement from surface and interface dislocation sources in NiAl

Limited ductility and toughness of NiAl and related aluminides near room temperature pose major problems in their potential application as structural materials. An analysis of these problems is presented as part of a review of the flow and fracture behavior of binary NiAl. Following this discussion is a demonstration that conditions of elastic and plastic constraint associated with phase boundaries afforded by surface films, internal lamellae, or precipitates may introduce sufficient densities of mobile dislocations to enhance the ductility of NiAl-based materials by significant amounts. Examples of this behavior are presented for several model materials, including 001- and 123-oriented single crystals of oxide-coated NiAl, directionally solidified beta-gamma-prime (Ni70Al30) and beta-gamma (Ni50Fe30Al20) in situ composites, and several NiAl/precipitate systems. The nature of the resulting dislocation substructures and the effects of several materials variables are described.

Noebe, R. D.↗

NASA NTRS1992

Combined bending and thermal fatigue of high-temperature metal-matrix composites - Computational simulation

Gotsis, Pascal K.↗

NASA NTRS1999

Fracture Probability of MEMS Optical Devices for Space Flight Applications

A bending fracture test specimen design is presented for thin elements used in optical devices for space flight applications. The specimen design is insensitive to load position, avoids end effect complications, and can be used to measure strength of membranes less than 2 microns thick. The theoretical equations predicting stress at failure are presented, and a detailed finite element model is developed to validate the equations for this application. An experimental procedure using a focused ion beam machine is outlined, and results from preliminary tests of 1.9 microns thick single crystal silicon are presented. These tests are placed in the context of a methodology for the design and evaluation of mission critical devices comprised of large arrays of cells.

Fettig, Rainer K.↗

NASA NTRS2014

Simulating Damage Due to a Lightning Strike Event: Effects of Temperature Dependent Properties on Interlaminar Damage

A multidirectional, carbon fiber-epoxy, composite panel is subjected to a simulated lightning strike, within a finite element method framework, and the effect of material properties on the failure (delamination) response is investigated through a detailed numerical study. The numerical model of the composite panel consists of individual homogenized plies with user-defined, cohesive interface elements between them. Lightning strikes are simulated as an assumed combination of excessive heat and high pressure loadings. It is observed that the initiation and propagation of lightning-induced delamination is a significant function of the temperature dependency of interfacial fracture toughness. This dependency must be defined properly in order to achieve reliable predictions of the present lightning-induced delamination in the composite panel.

computational mechanics↗

NASA NTRS2018

Implementation of J-A Methodology Elastic-Plastic Crack Instability Analysis Capability into the WARP-3D Code

Characterization of the near crack-tip stress/strain fields is the foundation of fracture mechanics. The description of the near tip stress field and the prediction of when fracture occurs is well established for brittle materials that exhibit linear elastic behavior. However, in ductile materials or conditions that violate linear elastic assumptions (Aluminum alloys, Al 2024-T3, Al 2024- T351 etc.), the elastic-plastic crack-tip stress fields are characterized by the Hutchison-Rice-Rosengren (HRR) field. The J-integral is commonly used to characterize amplitude of the HRR field under elastic-plastic conditions. The J-integral has been demonstrated for crack-tip fields that are under high constraint conditions (i.e., small-scale plasticity where the J-dominance is maintained). However, as the external load increases, yielding changes from small- to largescale plasticity and usually a loss of constraint (i.e., reduction in the triaxial stress field along the crack front). The loss of constraint leads to the deviation of the crack-tip stress fields from that given by the HRR field. Hence, the J-dominance will be gradually lost and additional parameter(s) are required to quantify the crack-tip stress fields and predict fracture behavior. The assessment objectives were to: 1) implement a two-parameter (i.e., J-A) fracture criterion into an elastic-plastic three-dimensional (3D) finite element analysis (FEA), 2) validate the implementation by comparison with the A parameter from literature data, 3) conduct material characterization tests to quantify the material behavior and provide fracture data for validation of the J-A fracture criteria, and (4) perform evaluations to establish if the J-A criteria can be used to predict fracture in a ductile metallic material (e.g., aluminum alloys). The A parameter in these criteria is the second parameter in a three-term elastic-plastic asymptotic expansion of the neartip stress behavior. A series of extensive FEAs were performed using WARP3D software package to obtain solutions for the A parameter for different specimen configurations. The methodology needed for the estimation of the A parameter in the asymptotic expansion was developed and implemented using Matlab®. A user material (UMAT) routine was used to model the material stress-strain response using a Ramberg-Osgood power law with a hardening exponent (n) and a material coefficient (alpha). This UMAT routine was successfully implemented in WARP3D software and validated through comparison with the experimental data. Three configurations were extracted from published results: 1) center cracked plate (CCP), 2) single edge-cracked plate (SECP), and 3) double edge-cracked plate (DECP). These configurations and four other configurations (three-hole tension (THT)), three-point bend (3PTB), three-hole compact tension (3PCT), and compact tension (CT)) were analyzed to verify the methodology that was developed and implemented into WARP3D. Solutions of the A parameter were obtained for remote tension loading conditions that started with small-scale yielding and continued into the large-scale plasticity regime. The results indicate that the methodology developed can be used to calculate the elastic-plastic J-A parameters for test specimens with a range of crack geometries, material strain hardening behaviors, and loading conditions. The J-A parameters were implemented as fracture criteria and used to predict the test results. For comparison, other fracture criteria were used to predict the same test results. Major findings include: The A constraint parameter A varies with specimen type and applied load thus accurate determination is crucial in predicting the failure load, and the A parameter is asymptotic as the failure load is approached, making an accurate determination difficult (i.e., small differences in the A parameter can cause large variations in failure load) for materials exhibiting elastic-plastic behavior. The failure predictions from J-A methodology were more accurate than the traditionally used KC and J methods, and have comparable scatter to that observed when using the crack-tip opening angle (CTOA) method. However, the J-A methodology requires considerable effort (expertise level and labor) to implement and to evaluate the A parameter for different specimen types and materials, or to apply this methodology to part-through crack (e.g., 3D problems) structural applications.

Hamm, Kenneth R., Jr.↗

NASA NTRS1989

A physical model for earthquakes. III - Thermodynamical approach and its relation to nonclassical theories

An earthquake theory based on fundamental thermodynamical considerations was developed using two principles. The first principle is that physical systems, including fault systems, tend toward a state of minimal Helmholtz free energy, which corresponds to static equilibrium. The second postulate is that, if the system is not at minimal free energy, thermodynamic forces arise which drive the system back toward a minimum. According to this theory, the propagation of a fracture, or the formation of a slip on a frictional surface, are manifestations of the physics which attempt to return the system to minimal free energy. It is shown that the process of return to minimal free energy by the system corresponds to the production of entropy, thereby rendering the process irreversible.

Rundle, John B.↗

NASA NTRS1991

Effect of sinter fracture and ohmic resistance on capacity retention in the nickel oxide electrode

The lifetime of batteries which utilize the nickel oxide electrode is often limited because this electrode loses a significant portion of its capacity as it is cycled. It is asserted that this capacity loss may often be attributed to cracking or separation of the conductive nickel sinter in the electrode, which forces electronic current to pass through the poorly conducting hydrated oxide and thus imposes a significant ohmic resistance. The model indicates that the oxide develops a nearly insulating layer which prevents complete discharge in the cycled electrode at usable rates. The capacity retention can be improved by reducing the cyclic stresses or strengthening the current collecting structure, redistributing it to provide a shorter current path through the solid phase, or by increasing the conductivity of the oxide to delay the formation of an insulating layer.

Lanzi, Oscar↗

NASA NTRS2002

Intralaminar and Interlaminar Progressive Failure Analysis of Composite Panels with Circular Cutouts

A progressive failure methodology is developed and demonstrated to simulate the initiation and material degradation of a laminated panel due to intralaminar and interlaminar failures. Initiation of intralaminar failure can be by a matrix-cracking mode, a fiber-matrix shear mode, and a fiber failure mode. Subsequent material degradation is modeled using damage parameters for each mode to selectively reduce lamina material properties. The interlaminar failure mechanism such as delamination is simulated by positioning interface elements between adjacent sublaminates. A nonlinear constitutive law is postulated for the interface element that accounts for a multi-axial stress criteria to detect the initiation of delamination, a mixed-mode fracture criteria for delamination progression, and a damage parameter to prevent restoration of a previous cohesive state. The methodology is validated using experimental data available in the literature on the response and failure of quasi-isotropic panels with centrally located circular cutouts loaded into the postbuckling regime. Very good agreement between the progressive failure analyses and the experimental results is achieved if the failure analyses includes the interaction of intralaminar and interlaminar failures.

Goyal, Vinay K.↗

NASA NTRS2011

Towards Simulating a Realistic Planetary Seismic Wavefield: The Contribution of the Megaregolith and Low-Velocity Waveguides

Lunar seismograms are distinctly different from their terrestrial counterparts. The Apollo lunar seismometers recorded moonquakes without distinct P- or S-wave arrivals; instead waves arrive as a diffuse coda that decays over several hours making the identification of body waves difficult. The unusual character of the lunar seismic wavefield is generally tied to properties of the megaregolith: it consists of highly fractured and broken crustal rock, the result of extensive bombardment of the Moon. The megaregolith extends several kilometers into the lunar crust, possibly into the mantle in some regions, and is covered by a thin coating of fine-scale dust. These materials possess very low seismic velocities that strongly scatter the seismic wavefield at high frequencies. Directly modeling the effects of the megaregolith to simulate an accurate lunar seismic wavefield is a challenging computational problem, owing to the inherent 3-D nature of the problem and the high frequencies (greater than 1 Hz) required. Here we focus on modeling the long duration code, studying the effects of the low velocities found in the megaregolith. We produce synthetic seismograms using 1-D slowness integration methodologies, GEMINI and reflectivity, and a 3-D Cartesian finite difference code, Wave Propagation Program, to study the effect of thin layers of low velocity on the surface of a planet. These codes allow us generate seismograms with dominant frequencies of approximately 1 Hz. For background lunar seismic structure we explore several models, including the recent model of Weber et al., Science, 2011. We also investigate variations in megaregolithic thickness, velocity, attenuation, and seismogram frequency content. Our results are compared to the Apollo seismic dataset, using both a cross correlation technique and integrated envelope approach to investigate coda decay. We find our new high frequency results strongly support the hypothesis that the long duration of the lunar seismic codes is generated by the presence of the low velocity megaregolith, and that the diffuse arrivals are a combination of scattered energy and multiple reverberations within this layer. The 3-D modeling indicates the extreme surface topography of the Moon adds only a small contribution to scattering effects, though local geology may play a larger role. We also study the effects of the megaregolith on core reflected and converted phases and other body waves. Our analysis indicates detection of core interacting arrivals with a polarization filter technique is robust and lends the possibility of detecting other body waves from the Moon.

Schmerr, Nicholas C.↗

NASA NTRS1988

The constitutive representation of high-temperature creep damage

The elastic-viscoplastic constitutive equations of Bodner-Partom were applied to modeling creep damage in a high temperature Ni-alloy, B1900 + Hf. Both tertiary creep in bulk materials and creep crack growth in flawed materials were considered. In the latter case, the energy rate line integral was used for characterizing the crack driving force, and the rate of crack extension was computed using a local damage formulation that assumed fracture was controlled by cavitation occurring within the crack-tip process zone. The results of this investigation were used to assess the evolution equation for isotropic damage utilized in the Bodner-Partom constitutive equations.

Chan, K. S.↗

NASA NTRS1976

A theory for creep crack growth

The paper derives an analytical expression for creep crack growth rate based on a model in which the higher and concentrated stresses ahead of the crack tip enhance creep deformation, thus progressively causing stress rupture and continuous crack advance. The equation derived for creep crack growth rate is expressed in terms of a geometrical factor in the stress intensity expression and in terms of the crack tip radius whose lower limit estimate is the fracture mechanics crack-tip opening displacement. The functional features of the derived creep crack growth rate equation can be compared with those of a pertinent empirical equation. The analytical expression provides some guidelines for alloy and microstructure design.

Purushothaman, S.↗

NASA NTRS1988

Dynamic fracture toughness of ceramic composites

The dynamic fracture toughness vs crack velocity relationships of TiB2 particulate-reinforced SiC-matrix and SiC whisker-reinforced Al2O3-matrix composites were determined at both room temperature and 1200 C with impacted, single-edge notched three-point bend specimens. Rapid crack initiation and propagation were monitored by a laser interferometric-displacement gage system. A FEM model that transmitted the measured outside-impact load to the specimen within the furnace was used to characterize the entire loading system-specimen unit. Small differences were measured between the room temperature and 1200 C dynamic responses.

Yang, K. H.↗

NASA NTRS1993

Multi-disciplinary coupling effects for integrated design of propulsion systems

Effective computational simulation procedures are described for modeling the inherent multi-disciplinary interactions which govern the accurate response of propulsion systems. Results are presented for propulsion system responses including multi-disciplinary coupling effects using coupled multi-discipline thermal, structural, and acoustic tailoring; an integrated system of multi-disciplinary simulators; coupled material behavior/fabrication process tailoring; sensitivities using a probabilistic simulator; and coupled materials, structures, fracture, and probabilistic behavior simulator. The results demonstrate that superior designs can be achieved if the analysis/tailoring methods account for the multi-disciplinary coupling effects. The coupling across disciplines can be used to develop an integrated coupled multi-discipline numerical propulsion system simulator.

Chamis, C. C.↗

NASA NTRS1997

Thermomechanical Behavior of Advanced SiC Fiber Multifilament Tows

In order to relate single fiber behavior to multiple fiber behavior in composites, fast-fracture tensile strength, creep, and stress-rupture studies were conducted on advanced SiC fiber multifilament tows in the temperature range from 20 to 1400 C in air as well as in inert environments. For conditions of small fiber creep (short times and low temperatures), the tow results of this study confirm the ability of limited single fiber data to model the strength behavior of multiple fibers in a bundle. For conditions of high creep (long times and high temperatures), further studies are needed to explain tow rupture behavior being better than average single fiber behavior.

Yun, Hee Mann↗