Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational fracture modeling”

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 73 records · Page 4

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.↗

Separating the Influence of Environment from Stress Relaxation Effects on Dwell Fatigue Crack Growth in a Nickel-Base Disk Alloy

Both environmental embrittlement and crack tip visco-plastic stress relaxation play a significant role in determining the dwell fatigue crack growth (DFCG) resistance of nickel-based disk superalloys. In the current study performed on the Low Solvus High Refractory (LSHR) disk alloy, the influence of these two mechanisms were separated so that the effects of each could be quantified and modeled. Seven different microstructural variations of LSHR were produced by controlling the cooling rate and the subsequent aging and thermal exposure heat treatments. Through cyclic fatigue crack growth testing performed both in air and vacuum, it was established that four out of the seven LSHR heat treatments evaluated, possessed similar intrinsic environmental resistance to cyclic crack growth. For these four heat treatments, it was further shown that the large differences in dwell crack growth behavior which still persisted, were related to their measured stress relaxation behavior. The apparent differences in their dwell crack growth resistance were attributed to the inability of the standard linear elastic fracture mechanics (LEFM) stress intensity parameter to account for visco-plastic behavior. Crack tip stress relaxation controls the magnitude of the remaining local tensile stresses which are directly related to the measured dwell crack growth rates. It was hypothesized that the environmentally weakened grain boundary crack tip regions fail during the dwells when their strength is exceeded by the remaining local crack tip tensile stresses. It was shown that the classical creep crack growth mechanisms such as grain boundary sliding did not contribute to crack growth, but the local visco-plastic behavior still plays a very significant role by determining the crack tip tensile stress field which controls the dwell crack growth behavior. To account for the influence of the visco-plastic behavior on the crack tip stress field, an empirical modification to the LEFM stress intensity parameter, Kmax, was developed by incorporating into the formulation the remaining stress level concept as measured by simple stress relaxation tests. The newly proposed parameter, Ksrf, did an excellent job in correlating the dwell crack growth rates for the four heat treatments which were shown to have similar intrinsic environmental cyclic fatigue crack growth resistance.

Dwell crack growth↗

Porous Media Approach for Modeling Closed Cell Foam

In order to minimize boil off of the liquid oxygen and liquid hydrogen and to prevent the formation of ice on its exterior surface, the Space Shuttle External Tank (ET) is insulated using various low-density, closed-cell polymeric foams. Improved analysis methods for these foam materials are needed to predict the foam structural response and to help identify the foam fracture behavior in order to help minimize foam shedding occurrences. This presentation describes a continuum based approach to modeling the foam thermo-mechanical behavior that accounts for the cellular nature of the material and explicitly addresses the effect of the internal cell gas pressure. A porous media approach is implemented in a finite element frame work to model the mechanical behavior of the closed cell foam. The ABAQUS general purpose finite element program is used to simulate the continuum behavior of the foam. The soil mechanics element is implemented to account for the cell internal pressure and its effect on the stress and strain fields. The pressure variation inside the closed cells is calculated using the ideal gas laws. The soil mechanics element is compatible with an orthotropic materials model to capture the different behavior between the rise and in-plane directions of the foam. The porous media approach is applied to model the foam thermal strain and calculate the foam effective coefficient of thermal expansion. The calculated foam coefficients of thermal expansion were able to simulate the measured thermal strain during heat up from cryogenic temperature to room temperature in vacuum. The porous media approach was applied to an insulated substrate with one inch foam and compared to a simple elastic solution without pore pressure. The porous media approach is also applied to model the foam mechanical behavior during subscale laboratory experiments. In this test, a foam layer sprayed on a metal substrate is subjected to a temperature variation while the metal substrate is stretched to simulate the structural response of the tank during operation. The thermal expansion mismatch between the foam and the metal substrate and the thermal gradient in the foam layer causes high tensile stresses near the metal/foam interface that can lead to delamination.

Ghosn, Louis J.↗

Fracture Mechanics Analyses of the Slip-Side Joggle Regions of Wing-Leading Edge Panels

The Space Shuttle Orbiter wing comprises of 22 leading edge panels on each side of the wing. These panels are part of the thermal protection system that protects the Orbiter wings from extreme heating that take place on the reentry in to the earth atmosphere. On some panels that experience extreme heating, liberation of silicon carbon (SiC) coating was observed on the slip side regions of the panels. Global structural and local fracture mechanics analyses were performed on these panels as a part of the root cause investigation of this coating liberation anomaly. The wing-leading-edge reinforced carbon-carbon (RCC) panels, Panel 9, T-seal 10, and Panel 10, are shown in Figure 1 and the progression of the stress analysis models is presented in Figure 2. The global structural analyses showed minimal interaction between adjacent panels and the T-seal that bridges the gap between the panels. A bounding uniform temperature is applied to a representative panel and the resulting stress distribution is examined. For this loading condition, the interlaminar normal stresses showed negligible variation in the chord direction and increased values in the vicinity of the slip-side joggle shoulder. As such, a representative span wise slice on the panel can be taken and the cross section can be analyzed using plane strain analysis.

Raju, Ivatury S.↗

Isotopic Constraints on the Genesis of Carbonates in Martian Meteorite ALH 84001

Oxygen isotopic analyses in approximately 20 micrometer spots in a chemically diverse suite of carbonates from ALH 84001 show highly variable delta(exp 18)O values from +5.4 to +25.3%. The isotopic data are correlated with the major element composition of the carbonate. The earliest forming (Ca-rich) carbonates have the lowest delta(exp 18)O values and the late-forming Mg-rich carbonates have the highest delta(exp 18)O values. Two models that can explain the isotopic variation were investigated. The carbonates could have formed in a water-rich environment at relatively low, but highly variable temperatures. In this open-system case the lower limit to the temperature variation is approximately 125 C, with fluctuations of over 250 C possible within the constraints of the model, depending on fluid composition. Alternatively the data can be explained by a closed-system model in which carbonates precipitated from a limited amount of a CO2-rich fluid. This scenario can reproduce the range of isotopic values observed, even at relatively high temperatures (greater than 500 C). Thus, the oxygen isotopic compositions do not provide unequivocal evidence for formation of the carbonates at low temperature. Neither of these scenarios is consistent with a biological origin of the carbonates and their associated features. Olivine from ALH 84001 occurs as clusters within orthopyroxene adjacent to fractures containing disrupted carbonate globules and feldspathic shock glass. The inclusions are irregular in shape and range in size from approximately 40 micrometers to submicrometer. The olivine exhibits a limited range of chemical composition from Fo(sub 65) to Fo(sub 66). We measured delta(exp 18)O values of the olivine to be +5.1 +/- 1.4%, indistinguishable within uncertainty from the host orthopyroxene. The data suggest that the olivine formed at high temperature (greater than 800 C), and is probably unrelated to carbonate formation. Instead the olivine probably formed by metamorphic reactions involving either dehydration of hydrous silicates or reduction of opx + spinel. If the reactions took place after carbonate formation, they could have caused devolatilization of siderite, producing the magnetite in the ALH carbonate globules. This scenario is also inconsistent with a biological origin for the features in ALH 84001.

Leshin, Laurie A.↗

Mechanical behavior of 18 Ni 200 grade maraging steel at cryogenic temperatures

A comprehensive study was conducted to characterize the mechanical behavior of 18 Ni 200 grade maraging steel for cryogenic applications in the National Transonic Facility. Tensile, fatigue, impact, and fracture toughness properties were determined at room and cryogenic temperatures. Experimental variables included product form, specimen orientation, and variation in metallurgical condition accomplished with a grain refining heat treatment. Room temperature tensile yield strength was approximately 200 ksi, while the yield strength at -275 F was approximately 260 ksi. Charpy V-notch energy absorption values and fracture toughness values of 18 Ni steel were found to be dependent on temperature, product form, orientation, and metallurgical condition. Fracture toughness values ranged from 80 to 113 ksi-sq rt in. at room temperature and from 67 to 90 ksi-sq rt in. at -275 F. The results of this study indicate that 18 Ni 200 grade maraging steel exceeds the strength criterion of 150 ksi minimum yield strength for use in cryogenic wind tunnel models but is marginal with respect to the 85 ksi sq rt in fracture toughness criterion at -275 F.

Wagner, J. A.↗

Load redistribution considerations in the fracture of ceramic matrix composites

Using a macroscopic viewpoint, composite laminae are homogeneous orthotropic solids whose directional strengths are random variables. Incorporation of these random variable strengths into failure models, either interactive or noninteractive, allows for the evaluation of the lamina reliability under a given stress state. Using a noninteractive criterion for demonstration purposes, laminate reliabilities are calculated assuming previously established load sharing rules for the redistribution of load as the failure of laminae occur. The matrix cracking predicted by ACK theory is modeled to allow a loss of stiffness in the fiber direction. The subsequent failure in the fiber direction is controlled by a modified bundle theory. Results are compared with previous models which did not permit separate consideration of matrix cracking, as well as to results obtained from experimental data. The effects of variations from the ideal physical geometry which is normally used to depict the matrix cracking are also studied.

Thomas, David J.↗

Three-dimensional hybrid-stress finite element analysis of composite laminates with cracks and cutouts

A three-dimensional hybrid-stress finite element analysis of composite laminates containing cutouts and cracks is presented. Fully three-dimensional, hexahedral isoparametric elements of the hybrid-stress model are formulated on the basis of the Hellinger-Reissner variational principle. Traction-free edges, cutouts, and crack surfaces are modeled by imposition of exact traction boundary conditions along element surfaces. Special boundary and surface elements are constructed by introducing proper constraints on assumed stress functions. The Lagrangian multiplier technique is used to enforce ply-interface continuity conditions in hybrid bimaterial composite elements for modeling the interface region in a composite laminate. Two examples are given to illustrate the capability of the present method of approach: (1) the well-known delamination problem in an angle-ply laminate, and (2) the important problem of a composite laminate containing a circular hole. Results are presented in detail for each case. Implications of interlaminar and intralaminar crack initiation, growth and fracture in composites containing cracks and cutouts are discussed.

Wang, S. S.↗

The delineation and interpretation of the Earth's gravity field

The geoid and topographic fields of the central Pacific were delineated and shown to correlate closely at intermediate wavelengths (500 to 2500 km). The associated admittance shows that anomalies having wavelengths less than about 1000 km are probably supported by the elastic strength of the lithosphere. Larger wavelength anomalies are due to dynamic effects in the sublithosphere. Direct modeling of small scale convection in the asthenosphere shows that the amplitudes of observed geoid and topographic anomalies can be independently matched, but that the observed admittance cannot. Only by imposing an initial regional variation in the thermal regime is it possible to match the admittance. It is proposed that this variation may be due to differences in the onset time of convection beneath the lithosphere of different ages. That is, convection beneath thickening lithosphere is strongly dependent on the rate of thickening (V) relative to the rise time for convection. The critical Rayleigh number contains the length scale K/V, where K is thermal diffusivity. Young, fast growing lithosphere stabilizes the underlying asthenosphere unless it has an unusually low viscosity. Lithosphere of different age, separated by fracture zones, will go unstable at different times, producing regional horizontal temperature gradient that may strongly influence convection. Laboratory and numerical experiments are proposed to study this form of convection and its influence on the geoid.

Marsh, Bruce D.↗

In-situ x-ray monitoring of damage accumulation in SiC/RBSN tensile specimens

The room-temperature tensile testing of silicon carbide fiber reinforced reaction-bonded silicon nitride (SiC/RBSN) composite specimens was monitored by using in-situ x ray film radiography. Radiographic evaluation before, during, and after loading provided data on the effect of preexisting volume flaws (high density impurities, and local density variations) on the fracture behavior of composites. Results from (0)1, (0)3, (0)5, and (0)8 composite specimens, showed that x ray film radiography can monitor damage accumulations during tensile loading. Matrix cracking, fiber-matrix debonding, and fiber pullout were imaged throughout the tensile loading history of the specimens. Further, in-situ film radiography was found to be a helpful and practical technique for estimating interfacial shear strength between the SiC fiber and the RBSN matrix by the matrix crack spacing method. It is concluded that pretest, in-situ, and post-test radiography can provide for a greater understanding of ceramic matrix composite mechanical behavior, a verification of related experimental procedures, and a validation and development of related analytical models.

Baaklini, George Y.↗

Application of fiber bridging models to fatigue crack growth in unidirectional titanium matrix composites

Several fiber bridging models were reviewed and applied to study the matrix fatigue crack growth behavior in center notched (0)(sub 8) SCS-6/Ti-15-3 and (0)(sub 4) SCS-6/Ti-6Al-4V laminates. Observations revealed that fatigue damage consisted primarily of matrix cracks and fiber matrix interfacial failure in the (0)(sub 8) SCS-6/Ti-15-3 laminates. Fiber-matrix interface failure included fracture of the brittle reaction zone and cracking between the two carbon rich fiber coatings. Intact fibers in the wake of the matrix cracks reduce the stress intensity factor range. Thus, an applied stress intensity factor range is inappropriate to characterize matrix crack growth behavior. Fiber bridging models were used to determine the matrix stress intensity factor range in titanium metal matrix composites. In these models, the fibers in the wake of the crack are idealized as a closure pressure. An unknown constant frictional shear stress is assumed to act along the debond or slip length of the bridging fibers. The frictional shear stress was used as a curve fitting parameter to available data (crack growth data, crack opening displacement data, and debond length data). Large variations in the frictional shear stress required to fit the experimental data indicate that the fiber bridging models in their present form lack predictive capabilities. However, these models provide an efficient and relatively simple engineering method for conducting parametric studies of the matrix growth behavior based on constituent properties.

Bakuckas, J. G., Jr.↗

In situ X-ray monitoring of damage accumulation in SiC/RBSN tensile specimens

The room-temperature tensile testing of silicon carbide fiber reinforced reaction-bonded silicon nitride (SiC/RBSN) composite specimens was monitored by using in-situ X-ray film radiography. Radiographic evaluation before, during, and after loading provided data on the effect of preexisting volume flaws (high density impurities, and local density variations) on the fracture behavior of composites. Results from (O)1, (O)3, (O)5, and (O)8 composite specimens showed that X-ray film radiography can monitor damage accumulations during tensile loading. Matrix cracking, fiber-matrix debonding, and fiber pullout were imaged throughout the tensile loading history of the specimens. Further, in-situ film radiography was found to be a helpful and practical technique for estimating interfacial shear strength between the SiC fiber and the RBSN matrix by the matrix crack spacing method. It is concluded that pretest, in-situ, and post-test radiography can provide for a greater understanding of ceramic matrix composite mechanical behavior, a verification of related experimental procedures, and a validation and development of related analytical models.

Baaklini, George Y.↗

Variational mechanics analysis of the stresses in microdrop debond specimens

A recently derived variational mechanics analysis of stresses in single-fiber model composites has been applied to the analysis of the stresses in the microdrop debond specimen. The new analysis is more accurate than the commonly applied shear-lag or elastic-plastic analyses. The results from a sample stress state calculation suggest that interfacial failure between the fiber and the microdrop is by mode I or opening mode failure at the beginning of the microdrop. The opening mode failure is caused by a large tensile radial stress at the fiber/matrix interface. Previous analyses of microdrop debond data have been in terms of a shear strength. It is suggested that these analyses misrepresent microdrop debond results and recommend instead a failure analysis based on energy release rate and interfacial fracture toughness. A procedure for calculating the energy release rate for the growth of an interfacial crack is described.

Scheer, Robert J.↗

Determination of Stress Intensity Factor Distributions for "Interface" Cracks in Incompressible, Dissimilar Materials

The present study was undertaken in order to develop test methods and procedures for measuring the variation of the stress intensity factor through the thickness in bimaterial specimens containing cracks within and parallel to the bond line using the frozen stress photoelastic method. Since stress freezing materials are incompressible above critical temperature, and since thick plates are to be employed which tend to produce a state of plane strain near the crack tip, the interface near tip fracture equations reduce to the classic form for homogeneous materials. Moreover, zero thickness interfaces do not exist when materials are bonded together. It was decided early on that it would be important to insure a uniform straight and accurate crack tip region through the thickness of the body to reduce scatter in the SIF distribution through the thickness. It was also observed that rubberlike materials which were desired to be modeled exhibited significant tip blunting prior to crack extension and that some blunting of the tip would provide a more realistic model. It should be noted that, in normal stress freezing photoelastic work, it is considered good practice to avoid utilizing data near bond lines in photoelastic models due to the bond line stresses which inevitably develop when two parts are bonded together. Thus, the present study involves certain exploratory aspects in deviating from standard practice in stress freezing work. With the above ideas in mind, several different test methods were investigated and are described in the following sections and appendices. The geometry selected for the program was a thick, edge cracked specimen containing a bond line.

Smith, C. W.↗

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.↗

Observation of pressure ridges in SAR images of sea ice: Scattering theory and comparison with observations

Ridges and keels (hummocks and bummocks) in sea ice flows are important in sea ice research for both scientific and practical reasons. Sea ice movement and deformation is driven by internal and external stresses on the ice. Ridges and keels play important roles in both cases because they determine the external wind and current stresses via drag coefficients. For example, the drag coefficient over sea ice can vary by a factor of several depending on the fluid mechanical roughness length of the surface. This roughness length is thought to be strongly dependent on the ridge structures present. Thus, variations in ridge and keel structure can cause gradients in external stresses which must be balanced by internal stresses and possibly fracture of the ice. Ridging in sea ice is also a sign of fracture. In a practical sense, large ridges form the biggest impediment to surface travel over the ice or penetration through sea ice by ice-strengthened ships. Ridges also play an important role in the damage caused by sea ice to off-shore structures. Hence, observation and measurement of sea ice ridges is an important component of sea ice remote sensing. The research reported here builds on previous work, estimating the characteristics of ridges and leads in sea ice from SAR images. Our objective is to develop methods for quantitative measurement of sea ice ridges from SAR images. To make further progress, in particular, to estimate ridge height, a scattering model for ridges is needed. Our research approach for a ridge scattering model begins with a survey of the geometrical properties of ridges and a comparison with the characteristics of the surrounding ice. For this purpose we have used airborne optical laser (AOL) data collected during the 1987 Greenland Sea Experiment. These data were used to generate a spatial wavenumber spectrum for height variance for a typical ridge - the typical ridge is the average over 10 large ridges. Our first-order model radar scattering includes both the quasi-specular and Bragg resonant scatter mechanisms. This model is extended to include contributions from volume scatter and scatter from discrete objects. Geometrical characteristics from the AOL survey and model calculations imply that for radar wavelengths and observation geometries that are dominated by the quasi-specular scattering mechanism radar backscatter from a ridge is a measure of peak ridge height. We present scattering model results and compare them with ridges observed during the LEADEX experiment of March-April 1992 when both X, C, and L-band aircraft SAR and the ERS-1 satellite SAR observed a region in the Beaufort Sea near 86 deg N, 10 deg W. Data were also collected documenting ridge characteristics on the surface. The surface data are used to generate a SAR signature via the scattering model described above. The predicted SAR signatures compare well with the SAR observations.

Vesecky, J. F.↗

Characterization of Low-Velocity Impact Damage in Toughened Non-Crimp Fabric Laminates

The damage response of resin-infused composites was investigated using quasi-static indentation (QSI) and instrumented drop-weight impact testing. The test specimens were composed of six layers of non-crimp fabrics where each layer contained +45°, 90°, −45°, and 0° plies. The goals of the testing were to identify damage responses that may differ from traditional thermoset materials and provide damage characterization data to evaluate current damage prediction analyses and aid improvements to current damage models. The damage was characterized using X-ray computed tomography (CT) and ultrasonic testing (UT). Incremental QSI loading was used on a single specimen to investigate damage growth in nine steps at progressively higher loads. Realtime acoustic emission data were used to determine the unloading points. Unlike typical thermoset composites, the force-displacement was smooth and showed no obvious indications of damage. The primary damage consisted of large lines of fiber fracture on the plies near the impact surface and delaminations through the thickness. Nine specimens were impacted at nominal impact energies of 25 ft-lbs, 30 ft-lbs, and 35 ft-lbs. The force histories at each energy followed nearly identical paths. Again, no obvious signs of damage were observed in any of the force histories. The general damage response in the dynamic tests was similar to the QSI-loaded specimens. However, repeat specimens impacted with a given impact energy showed variations in the damage. Some characteristics of the damage response were different than those typically observed in thermoset composites and will need to be evaluated using current analysis methods.

Composite Materials↗

Characterization of Low-Velocity Impact Damage in Toughened Non-Crimp Fabric Laminates

The damage response of resin-infused composites was investigated using quasi-static indentation (QSI) and instrumented drop-weight impact testing. The test specimens were composed of six layers of non-crimp fabrics where each layer contained +45°, 90°, −45°, and 0° plies. The goals of the testing were to identify damage responses that may differ from traditional thermoset materials and provide damage characterization data to evaluate current damage prediction analyses and aid improvements to current damage models. The damage was characterized using X-ray computed tomography (CT) and ultrasonic testing (UT). Incremental QSI loading was used on a single specimen to investigate damage growth in nine steps at progressively higher loads. Realtime acoustic emission data were used to determine the unloading points. Unlike typical thermoset composites, the force-displacement was smooth and showed no obvious indications of damage. The primary damage consisted of large lines of fiber fracture on the plies near the impact surface and delaminations through the thickness. Nine specimens were impacted at nominal impact energies of 25 ft-lbs, 30 ft-lbs, and 35 ft-lbs. The force histories at each energy followed nearly identical paths. Again, no obvious signs of damage were observed in any of the force histories. The general damage response in the dynamic tests was similar to the QSI-loaded specimens. However, repeat specimens impacted with a given impact energy showed variations in the damage. Some characteristics of the damage response were different than those typically observed in thermoset composites and will need to be evaluated using current analysis methods.

Composite Materials↗