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 109 records · Page 6

Elastic-Plastic J-Integral Solutions or Surface Cracks in Tension Using an Interpolation Methodology

No closed form solutions exist for the elastic-plastic J-integral for surface cracks due to the nonlinear, three-dimensional nature of the problem. Traditionally, each surface crack must be analyzed with a unique and time-consuming nonlinear finite element analysis. To overcome this shortcoming, the authors have developed and analyzed an array of 600 3D nonlinear finite element models for surface cracks in flat plates under tension loading. The solution space covers a wide range of crack shapes and depths (shape: 0.2 less than or equal to a/c less than or equal to 1, depth: 0.2 less than or equal to a/B less than or equal to 0.8) and material flow properties (elastic modulus-to-yield ratio: 100 less than or equal to E/ys less than or equal to 1,000, and hardening: 3 less than or equal to n less than or equal to 20). The authors have developed a methodology for interpolating between the goemetric and material property variables that allows the user to reliably evaluate the full elastic-plastic J-integral and force versus crack mouth opening displacement solution; thus, a solution can be obtained very rapidly by users without elastic-plastic fracture mechanics modeling experience. Complete solutions for the 600 models and 25 additional benchmark models are provided in tabular format.

Allen, P. A.↗

Elastic-Plastic J-Integral Solutions or Surface Cracks in Tension Using an Interpolation Methodology. Appendix C -- Finite Element Models Solution Database File, Appendix D -- Benchmark Finite Element Models Solution Database File

No closed form solutions exist for the elastic-plastic J-integral for surface cracks due to the nonlinear, three-dimensional nature of the problem. Traditionally, each surface crack must be analyzed with a unique and time-consuming nonlinear finite element analysis. To overcome this shortcoming, the authors have developed and analyzed an array of 600 3D nonlinear finite element models for surface cracks in flat plates under tension loading. The solution space covers a wide range of crack shapes and depths (shape: 0.2 less than or equal to a/c less than or equal to 1, depth: 0.2 less than or equal to a/B less than or equal to 0.8) and material flow properties (elastic modulus-to-yield ratio: 100 less than or equal to E/ys less than or equal to 1,000, and hardening: 3 less than or equal to n less than or equal to 20). The authors have developed a methodology for interpolating between the goemetric and material property variables that allows the user to reliably evaluate the full elastic-plastic J-integral and force versus crack mouth opening displacement solution; thus, a solution can be obtained very rapidly by users without elastic-plastic fracture mechanics modeling experience. Complete solutions for the 600 models and 25 additional benchmark models are provided in tabular format.

Allen, Phillip A.↗

A nondestructive technique for predicting the strength remaining in filament wound composites subjected to low-level impact

A model for predicting the fracture strength of homogeneous materials is proposed. Impacted FWC samples were evaluated using ultrasonic testing and an X-ray dye penetration method. The ability of the model to measure fracture strength was also examined. The relation between attenuation and velocity measurements is studied. It is observed that the X-ray method is not useful for predicting fracture strength because the dye could not penetrate the matrix. It is noted that fracture strength predictions derived from the fracture mechanical model and the ultrasonic measurements correlate well with actual measured fracture strengths.

Madaras, E. I.↗

On representation of fracture profiles by fractional integrals of a Wiener process

An enhancement of the crack diffusion model of brittle fracture is proposed. In the model, certain probabilities are evaluated as averages over the set of all virtual crack paths, and the latter are modeled as realization of a Wiener process. However, commonly observed fracture profiles are smoother. In fact, one can show that their fractal dimension is smaller than the value 3/2 for a Wiener process. Accordingly, it is proposed to represent fracture profiles by graphs of fractional integrals of a Wiener process. This allows one to employ additional fractographic information obtained from fracture profiles, e.g. their fractal dimension, in a predictive formalism, evaluating probabilities of failure.

Kunin, B.↗

Apollo 14 and 16 Active Seismic Experiments, and Apollo 17 Lunar Seismic Profiling

Seismic refraction experiments were conducted on the moon by Apollo astronauts during missions 14, 16, and 17. Seismic velocities of 104, 108, 92, 114 and 100 m/sec were inferred for the lunar regolith at the Apollo 12, 14, 15, 16, and 17 landing sites, respectively. These data indicate that fragmentation and comminution caused by meteoroid impacts has produced a layer of remarkably uniform seismic properties moonwide. Brecciation and high porosity are the probable causes of the very low velocities observed in the lunar regolith. Apollo 17 seismic data revealed that the seismic velocity increases very rapidly with depth to 4.7 km/sec at a depth of 1.4 km. Such a large velocity change is suggestive of compositional and textural changes and is compatible with a model of fractured basaltic flows overlying anorthositic breccias. 'Thermal' moonquakes were also detected at the Apollo 17 site, becoming increasingly frequent after sunrise and reaching a maximum at sunset. The source of these quakes could possibly be landsliding.

Source record↗

Status of the Ford program to evaluate ceramics for stator applications in automotive gas turbine engines

The paper reviews the progress of the major technical tasks of the DOE/NASA/Ford program Evaluation of Ceramics for Stator Applications in Automotive Gas Turbine Engines: reliability prediction, stator fabrication, material characterization, and stator evaluation. A fast fracture reliability model was prepared for a one-piece ceramic stator. Periodic inspection results are presented.

Trela, W.↗

Possible role of crustal flexure in the initial detachment of extensional allochthons

The existence of low-angle normal faults indicates that the ratio of shear stress (TAU) to normal stress (sigma sub N) needed to cause slip on faults is substantially less than would be predicted based on experimental data. Because the tensional strength of rock at a large scale is exceedingly low, the upper plate of a low-angle normal fault cannot be pulled down the fault ramp, but must be driven down it by its own weight. The active or recently active Sevier Desert detachment fault in western Utah dips regionally at 12 deg. The radio of shear stress to normal stress due to the weight of the upper plate on a 12 deg dipping fault surface is 0.2. In contrast, laboratory experiments indicate that slip on fracture surfaces occurs with almost all rock types when (tau/sigma) reaches values of 0.6 to 0.85, corresponding to normal-fault dips of 30 deg to 40 deg. Seismological data indicate that low deviatoric stresses are associated with movement on faults of other geometrics and are not unique to low-angle normal faults. It thus appears that approximately planar fault zones with surface areas of hundreds to thousands of square kilometers have different mechanical properties than would be predicted based on laboratory studies of fractured rock. Modeling of stresses indicates that deviatoric stresses greater than 5kb exist and are sustainable in oceanic lithosphere, and that failure occurs when tau/sigma sub N approaches 0.6 to 0.85.

Spencer, J. E.↗

Sublaminate analysis of interlaminar fracture in composites. I - Analytical model

A simple analysis method based on a transverse shear-deformation theory and a sublaminate approach is utilized to analyze a mixed-mode edge delamination specimen. The analysis provides closed-form expressions for the interlaminar shear stresses ahead of the crack, the total strain-energy release rate, and the strain-energy release-rate components. The parameters controlling the behavior are identified. The effect of specimen stacking sequence and delamination interface on the strain energy release rate components is investigated. Results are compared with a finite-element simulation for reference. The simple nature of the method makes it suitable for preliminary design analyses which require a large number of configurations to be evaluated quickly and economically.

Armanios, Erian A.↗

Matrix cracking in composite laminates with resin-rich interlaminar layers

Fracture mechanics analysis and test data for a toughened composite material that has a resin-rich interlaminar layer (RIL) were used to investigate in situ strength. Exposure to a range of environmental conditions was considered. A parametric analysis study was performed to judge the effects of laminate and material variables. A finite thickness effect, indicating an interaction between ply group thickness and effective flaw size, was found dominant. The magnitude of the effect was directly related to RIL stiffness. In situ strength was found to decrease with decreasing RIL stiffness. This work indicates the need to use a fracture mechanics model of actual lamina microstructure and heterogeneous properties to predict in situ strength in materials with RIL.

Ilcewicz, Larry B.↗

Damage and strength of composite materials: Trends, predictions, and challenges

Research on damage mechanisms and ultimate strength of composite materials relevant to scaling issues will be addressed in this viewgraph presentation. The use of fracture mechanics and Weibull statistics to predict scaling effects for the onset of isolated damage mechanisms will be highlighted. The ability of simple fracture mechanics models to predict trends that are useful in parametric or preliminary designs studies will be reviewed. The limitations of these simple models for complex loading conditions will also be noted. The difficulty in developing generic criteria for the growth of these mechanisms needed in progressive damage models to predict strength will be addressed. A specific example for a problem where failure is a direct consequence of progressive delamination will be explored. A damage threshold/fail-safety concept for addressing composite damage tolerance will be discussed.

Obrien, T. Kevin↗

Crack Turning Mechanics of Composite Wing Skin Panels

The safety of future composite wing skin integral stiffener panels requires a full understanding of failure mechanisms of these damage tolerance critical structures under both in-plane and bending loads. Of primary interest is to derive mathematical models using fracture mechanics in anisotropic cracked plate structures, to assess the crack turning mechanisms, and thereby to enhance the residual strength in the integral stiffener composite structures. The use of fracture mechanics to assess the failure behavior in a cracked structure requires the identification of critical fracture parameters which govern the severity of stress and deformation field ahead of the flaw, and which can be evaluated using information obtained from the flaw tip. In the three-year grant, the crack-tip fields under plane deformation, crack-tip fields for anisotropic plates and anisotropic shells have been obtained. In addition, methods for determining the stress intensity factors, energy release rate, and the T-stresses have been proposed and verified. The research accomplishments can be summarized as follows: (1) Under plane deformation in anisotropic solids, the asymptotic crack-tip fields have been obtained using Stroh formalism; (2) The T-stress and the coefficient of the second term for sigma(sub y), g(sub 32), have been obtained using path-independent integral, the J-integral and Betti's reciprocal theorem together with auxiliary fields; (3) With experimental data performed by NASA, analyses indicated that the mode-I critical stress intensity factor K(sub Q) provides a satisfactory characterization of fracture initiation for a given laminate thickness, provided the failure is fiber-dominated and crack extends in a self-similar manner; (4) The high constraint specimens, especially for CT specimens, due to large T-stress and large magnitude of negative g(sub 32) term may be expected to inhibit the crack extension in the same plane and promote crack turning; (5) Crack turning out of crack plane in generally anisotropic solids under plane deformation has been studied; (6) The role of T-stress and the higher-order term of sigma(sub y) on the crack turning and stability of the kinked crack has been quantified; (7) Asymptotic crack-tip fields including the effect of transverse shear deformation (Reissner plate theory) in an anisotropic plate under bending, twisting moments, and transverse shear loads has been presented; (8) The expression of the path-independent J-integral in terms of the generalized stress and strain has been derived; (9) Asymptotic crack-tip fields including the effect of transverse shear deformation (Reissner shallow shell theory) in a general anisotropic shell has been developed; (10) The Stroh formalism was used to characterize the crack tip fields in shells up to the second term and the energy release rate was expressed in a very compact form.

Yuan, F. G.↗

Failure Criteria for FRP Laminates in Plane Stress

A new set of six failure criteria for fiber reinforced polymer laminates is described. Derived from Dvorak's fracture mechanics analyses of cracked plies and from Puck's action plane concept, the physically-based criteria, denoted LaRC03, predict matrix and fiber failure accurately without requiring curve-fitting parameters. For matrix failure under transverse compression, the fracture plane is calculated by maximizing the Mohr-Coulomb effective stresses. A criterion for fiber kinking is obtained by calculating the fiber misalignment under load, and applying the matrix failure criterion in the coordinate frame of the misalignment. Fracture mechanics models of matrix cracks are used to develop a criterion for matrix in tension and to calculate the associated in-situ strengths. The LaRC03 criteria are applied to a few examples to predict failure load envelopes and to predict the failure mode for each region of the envelope. The analysis results are compared to the predictions using other available failure criteria and with experimental results. Predictions obtained with LaRC03 correlate well with the experimental results.

Davila, Carlos G.↗

Physically Based Failure Criteria for Transverse Matrix Cracking

A criterion for matrix failure of laminated composite plies in transverse tension and in-plane shear is developed by examining the mechanics of transverse matrix crack growth. Matrix cracks are assumed to initiate from manufacturing defects and can propagate within planes parallel to the fiber direction and normal to the ply mid-plane. Fracture mechanics models of cracks in unidirectional laminates, embedded plies and outer plies are developed to determine the onset and direction of propagation for unstable crack growth. The models for each ply configuration relate ply thickness and ply toughness to the corresponding in-situ ply strength. Calculated results for several materials are shown to correlate well with experimental results.

Davila, Carlos G.↗

Thermal and Environmental Barrier Coatings for Advanced Turbine Engine Applications

Ceramic thermal and environmental barrier coatings (T/EBCs) will play a crucial role in advanced gas turbine engine systems because of their ability to significantly increase engine operating temperatures and reduce cooling requirements, thus help achieve engine low emission and high efficiency goals. Advanced T/EBCs are being developed for the low emission SiC/SiC ceramic matrix composite (CMC) combustor applications by extending the CMC liner and vane temperature capability to 1650 C (3000 F) in oxidizing and water vapor containing combustion environments. Low conductivity thermal barrier coatings (TBCs) are also being developed for metallic turbine airfoil and combustor applications, providing the component temperature capability up to 1650 C (3000 F). In this paper, ceramic coating development considerations and requirements for both the ceramic and metallic components will be described for engine high temperature and high-heat-flux applications. The underlying coating failure mechanisms and life prediction approaches will be discussed based on the simulated engine tests and fracture mechanics modeling results.

Zhu, Dong-Ming↗

Distribution of Inclusion-Initiated Fatigue Cracking in Powder Metallurgy Udimet 720 Characterized

In the absence of extrinsic surface damage, the fatigue life of metals is often dictated by the distribution of intrinsic defects. In powder metallurgy (PM) alloys, relatively large defects occur rarely enough that a typical characterization with a limited number of small volume fatigue test specimens will not adequately sample inclusion-initiated damage. Counterintuitively, inclusion-initiated failure has a greater impact on the distribution in PM alloy fatigue lives because they tend to have fewer defects than their cast and wrought counterparts. Although the relative paucity of defects in PM alloys leads to higher mean fatigue lives, the distribution in observed lives tends to be broader. In order to study this important failure initiation mechanism without expending an inordinate number of specimens, a study was undertaken at the NASA Glenn Research Center where known populations of artificial inclusions (seeds) were introduced to production powder. Fatigue specimens were machined from forgings produced from the seeded powder. Considerable effort has been expended in characterizing the crack growth rate from inclusion-initiated cracks in seeded PM alloys. A rotating and translating positioning system, with associated software, was devised to map the surface inclusions in low-cycle fatigue (LCF) test bars and to monitor the crack growth from these inclusions. The preceding graph illustrates the measured extension in fatigue cracks from inclusions on a seeded LCF test bar subjected to cyclic loading at a strain range of 0.8 percent and a strain ratio (max/min) of zero. Notice that the observed inclusions fall into three categories: some do not propagate at all (arrest), some propagate with a decreasing crack growth rate, and a few propagate at increasing rates that can be modeled by fracture mechanics. The following graph shows the measured inclusion-initiated crack growth rates from 10 interrupted LCF tests plotted against stress intensities calculated for semi-elliptical cracks with the observed surface lengths. The expected scatter in the crack growth rates for stress intensity ranges near threshold is observed. These data will be used to help determine the distribution in growth rates of cracks emanating from inclusions as well as the proportion of cracks that arrest under various loading conditions.

Bonacuse, Peter J.↗

A Continuum-Atomistic Analysis of Transgranular Crack Propagation in Aluminum

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain is used to study plastic processes at a crack tip in a single crystal of aluminum. The case of mode I loading is studied. A transition from deformation twinning to full dislocation emission from the crack tip is found when the crack plane is rotated around the [111] crystallographic axis. When the crack plane normal coincides with the [112] twinning direction, the crack propagates through a twinning mechanism. When the crack plane normal coincides with the [011] slip direction, the crack propagates through the emission of full dislocations. In intermediate orientations, a transition from full dislocation emission to twinning is found to occur with an increase in the stress intensity at the crack tip. This finding confirms the suggestion that the very high strain rates, inherently present in MD simulations, which produce higher stress intensities at the crack tip, over-predict the tendency for deformation twinning compared to experiments. The present study, therefore, aims to develop a more realistic and accurate predictive modeling of fracture processes.

Yamakov, V.↗

Floating Node Method and Virtual Crack Closure Technique for Modeling Matrix Cracking-Delamination Interaction

A novel approach is proposed for high-fidelity modeling of progressive damage and failure in composite materials that combines the Floating Node Method (FNM) and the Virtual Crack Closure Technique (VCCT) to represent multiple interacting failure mechanisms in a mesh-independent fashion. In this study, the approach is applied to the modeling of delamination migration in cross-ply tape laminates. Delamination, matrix cracking, and migration are all modeled using fracture mechanics based failure and migration criteria. The methodology proposed shows very good qualitative and quantitative agreement with experiments.

DeCarvalho, N. V.↗

Floating Node Method and Virtual Crack Closure Technique for Modeling Matrix Cracking-Delamination Migration

A novel approach is proposed for high-fidelity modeling of progressive damage and failure in composite materials that combines the Floating Node Method (FNM) and the Virtual Crack Closure Technique (VCCT) to represent multiple interacting failure mechanisms in a mesh-independent fashion. In this study, the approach is applied to the modeling of delamination migration in cross-ply tape laminates. Delamination, matrix cracking, and migration are all modeled using fracture mechanics based failure and migration criteria. The methodology proposed shows very good qualitative and quantitative agreement with experiments.

DeCarvalho, Nelson V.↗