Raju, I. S. — Page 5 | Engineering Papers

NASA NTRS1986

Boundary force method for analyzing two-dimensional cracked bodies

The Boundary Force Method (BFM) was formulated for the two-dimensional stress analysis of complex crack configurations. In this method, only the boundaries of the region of interest are modeled. The boundaries are divided into a finite number of straight-line segments, and at the center of each segment, concentrated forces and a moment are applied. This set of unknown forces and moments is calculated to satisfy the prescribed boundary conditions of the problem. The elasticity solution for the stress distribution due to concentrated forces and a moment applied at an arbitrary point in a cracked infinite plate are used as the fundamental solution. Thus, the crack need not be modeled as part of the boundary. The formulation of the BFM is described and the accuracy of the method is established by analyzing several crack configurations for which accepted stress-intensity factor solutions are known. The crack configurations investigated include mode I and mixed mode (mode I and II) problems. The results obtained are, in general, within + or - 0.5 percent of accurate numerical solutions. The versatility of the method is demonstrated through the analysis of complex crack configurations for which limited or no solutions are known.

Tan, P. W.↗

NASA NTRS1986

Prediction of fatigue crack-growth patterns and lives in three-dimensional cracked bodies

Fatigue crack growth patterns and lives for surface cracks, surface cracks at holes, and corner cracks at holes in three dimensional bodies were predicted using linear-elastic fracture mechanics concepts that were modified to account for crack-closure behavior. The predictions were made by using stress intensity factor equations for these crack configurations and the fatigue crack-growth (delta K against rate) relationship for the material of interest. The crack configurations were subjected to constant-amplitude fatigue loading under either remote tension or bending loads. The predicted crack growth patterns and crack growth lives for aluminum alloys agreed well with test data from the literature.

Newman, J. C., Jr.↗

NASA NTRS1986

Stress-intensity factors for circumferential surface cracks in pipes and rods under tension and bending loads

The purpose of this paper is to present stress-intensity factors for a wide range of nearly semi-elliptical surface cracks in pipes and rods. The configurations were subjected to either remote tension or bending loads. For pipes, the ratio of crack depth to crack length (a/c) ranged from 0.6 to 1; the ratio of crack depth to wall thickness (a/t) ranged from 0.2 to 0.8; and the ratio of internal radius to wall thickness (R/t) ranged from 1 to 10. For rods, the ratio of crack depth to crack length also ranged from 0.6 to 1; and the ratio of crack depth to rod diameter (a/D) ranged from 0.05 to 0.35. These particular crack configurations were chosen to cover the range of crack shapes (a/c) that have been observed in experiments conducted on pipes and rods under tension and bending fatigue loads. The stress-intensity factors were calculated by a three-dimensional finite-element method. The finite-element models employed singularity elements along the crack front and linear-strain elements elsewhere. The models had about 6500 degrees of freedom. The stress-intensity factors were evaluated using a nodal-force method.

Raju, I. S.↗

NASA NTRS1986

Residual thermal and moisture influences on the strain energy release rate analysis of edge delamination

A laminated plate theory analysis is developed to calculate the strain energy release rate associated with edge delamination growth in a composite laminate. The analysis includes the contribution of residual thermal and moisture stresses to the strain energy released. The strain energy release rate, G, increased when residual thermal effects were combined with applied mechanical strains, but then decreased when increasing moisture content was included. A quasi-three-dimensional finite element analysis indicated identical trends and demonstrated these same trends for the individual strain energy release rate components, G sub I and G sub II, associated with interlaminar tension and shear. An experimental study indicated that for T300/5208 graphite-epoxy composites, the inclusion of residual thermal and moisture stresses did not significantly alter the calculation of interlaminar fracture toughness from strain energy release rate analysis of edge delamination data taken at room temperature, ambient conditions.

Obrien, T. K.↗

NASA NTRS1985

Stress-intensity factors for circumferential surface cracks in pipes and rods under tension and bending loads

The purpose of this paper is to present stress-intensity factors for a wide range of nearly semi-elliptical surface cracks in pipes and rods. The configurations were subjected to either remote tension or bending loads. For pipes, the ratio of crack depth to crack length (a/c) ranged from 0.6 to 1; the ratio of crack depth to wall thickness (a/t) ranged from 0.2 to 0.8; and the ratio of internal radius to wall thickness (R/t) ranged from 1 to 10. For rods, the ratio of crack depth to crack length also ranged from 0.6 to 1; and the ratio of crack depth to rod diameter (a/D) ranged from 0.05 to 0.35. These particular crack configurations were chosen to cover the range of crack shapes (a/c) that have been observed in experiments conducted on pipes and rods under tension and bending fatigue loads. The stress-intensity factors were calculated by a three-dimensional finite-element method. The finite-element models employed singularity elements along the crack front and linear-strain elements elsewhere. The models had about 6500 degrees of freedom. The stress-intensity factors were evaluated using a nodal-force method.

Raju, I. S.↗

NASA NTRS1985

Comparisons of various configurations of the edge delamination test for interlaminar fracture toughness

Various configurations of Edge Delamination Tension (EDT) test specimens, of both brittle (T300/5208) and toughened-matrix (T300/BP907) graphite reinforced composite laminates, were manufactured and tested. The mixed-mode interlaminar fracture toughness, G sub C, was measured using (30/30 sub 2/30/90 sub N)sub s, n=1 or 2, (35/-35/0/90) sub s and (35/0/-35/90) sub s layups designed to delaminate at low tensile strains. Laminates were made without inserts so that delaminations would form naturally between the central 90 deg plies and the adjacent angle plies. Laminates were also made with Teflon inserts implanted between the 90 deg plies and the adjacent angle (theta) plies at the straight edge to obtain a planar fracture surface. In addition, interlaminar tension fracture toughness, GIc, was measured from laminates with the same layup but with inserts in the midplane, between the central 90 deg plies, at the straight edge. All of the EDT configurations were useful for ranking the delamination resistance of composites with different matrix resins. Furthermore, the variety of layups and configurations available yield interlaminar fracture toughness measurements needed to generate delamination failure criteria. The influence of insert thickness and location, and coupon size on G sub c values were evaluated.

Obrien, T. K.↗

NASA NTRS1985

Residual thermal and moisture influences on the strain energy release rate analysis of edge delamination

A laminated plate theory analysis is developed to calculate the strain energy release rate associated with edge delamination growth in a composite laminate. The analysis includes the contribution of residual thermal and moisture stresses to the strain energy released. The strain energy release rate, G, increased when residual thermal effects were combined with applied mechanical strains, but then decreased when increasing moisture content was included. A quasi-three-dimensional finite element analysis indicated identical trends and demonstrated these same trends for the individual strain energy release rate components, G sub I and G sub II, associated with interlaminar tension and shear. An experimental study indicated that for T300/5208 graphite-epoxy composites, the inclusion of residual thermal and moisture stresses did not significantly alter the calculation of interlaminar fracture toughness from strain energy release rate analysis of edge delamination data taken at room temperature, ambient conditions.

Obrien, T. K.↗

NASA NTRS1985

Large deflection behavior of quasi-isotropic laminates under low-velocity impact type point loading

Eight-ply quasi-isotropic circular composite plates of Thornel 300 graphite in Narmco 5208 epoxy resin (T300/5208) were analyzed to obtain the large deformation behavior under low-velocity impact type point loads. A simple plate-membrane coupling model was developed. The impact type point loads were replaced by equivalent quasi-static point loads. The plate-membrane coupling model was used to obtain the large deformation shapes for the thin circular composite laminates. The analyses indicated that the large deformation shapes of the composite plates under point loads vary with the center point displacements, and hence are different for different load levels. Quasi-isotropic plates were analyzed by replacing anisotropic bending stiffness components with the equivalent flexural stiffness for the isotropic plates. The plate-membrane coupling model was verified by conducting a series of tests on clamped circular quasi-isotropic laminates. Deflected shapes for the thin composite plates were experimentally obtained. These shapes agreed well with the analytically predicted shapes.

Kelkar, A.↗

NASA NTRS1985

Large deflections of circular isotropic membranes subjected to arbitrary axisymmetric loading

Circular membranes with fixed peripheral edges, subjected to arbitrary axisymmetric loading were analyzed. A single governing differential equation in terms of radial stress was used. This nonlinear governing equation was solved using the finite difference method in conjunction with Newton-Raphson method. Three loading cases, namely: (1) uniformly loaded membrane, (2) a membrane with uniform load over an inner portion, and (3) a membrane with ring load, were analyzed. Calculated central displacement and the central and edge radial stresses for uniformly loaded membrane, agreed extremely well with the classical solution.

Kelkar, A.↗

NASA NTRS1985

Analysis of interlaminar stresses in thick composite laminates with and without edge delamination

The effect of laminate thickness on the interlaminar stresses in rectangular quasi-isotropic laminates under uniform axial strain was studied. Laminates from 8-ply to infinitely thick were analyzed. Thick laminates were synthesized by stacking (45/0/-45/90) ply groups, rather than grouping like plies. Laminates with and without delaminations were studied. In laminates without delaminations, the free-edge interlaminar normal stress distribution in the outer ply groups was insensitive to total laminate thickness. The interlaminar normal stress distribution for the interior ply groups was nearly the same as for an infinitely thick laminate. In contrast, the free-edge inter-laminar shear stress distribution was nearly the same for inner and outer ply groups and was insensitive to laminate thickness. In laminates with delaminations those delaminations near the top and bottom surfaces of a thick laminate have much larger total strain-energy-release rates (G sub t) and mode I-to-total (G sub t/G sub t) ratios than delaminations deep in the interior. Therefore, delaminations can be expected to grow more easily near the surfaces of a laminate than in the interior.

Whitcomb, J. D.↗

NASA NTRS1984

Large deflections of circular isotropic membranes subjected to arbitrary axisymmetric loading

Circular membranes with fixed peripheral edges, subjected to arbitrary axisymmetric loading are analyzed. A single governing differential equation in terms of radial stress is used. This nonlinear governing equation is solved using the finite difference method in conjunction with Newton-Raphson method. Three loading cases, namely (1) uniformly loaded membrane, (2) a membrane with uniform load over an inner portion, and (3) a membrane with ring load, are analyzed. Calculated central displacement and the central and edge radial stresses for uniformly loaded membrane, agree extremely well with the classical solution.

Kelkar, A.↗

NASA NTRS1984

Methods for analysis of cracks in three-dimensional solids

Analytical and numerical methods evaluating the stress-intensity factors for three-dimensional cracks in solids are presented, with reference to fatigue failure in aerospace structures. The exact solutions for embedded elliptical and circular cracks in infinite solids, and the approximate methods, including the finite-element, the boundary-integral equation, the line-spring models, and the mixed methods are discussed. Among the mixed methods, the superposition of analytical and finite element methods, the stress-difference, the discretization-error, the alternating, and the finite element-alternating methods are reviewed. Comparison of the stress-intensity factor solutions for some three-dimensional crack configurations showed good agreement. Thus, the choice of a particular method in evaluating the stress-intensity factor is limited only to the availability of resources and computer programs.

Raju, I. S.↗

NASA NTRS1984

Three-dimensional finite-element analysis of chevron-notched fracture specimens

Stress-intensity factors and load-line displacements were calculated for chevron-notched bar and rod fracture specimens using a three-dimensional finite-element analysis. Both specimens were subjected to simulated wedge loading (either uniform applied displacement or uniform applied load). The chevron-notch sides and crack front were assumed to be straight. Crack-length-to-specimen width ratios (a/w) ranged from 0.4 to 0.7. The width-to-thickness ratio (w/B) was 1.45 or 2. The bar specimens had a height-to-width ratio of 0.435 or 0.5. Finite-element models were composed of singularity elements around the crack front and 8-noded isoparametric elements elsewhere. The models had about 11,000 degrees of freedom. Stress-intensity factors were calculated by using a nodal-force method for distribution along the crack front and by using a compliance method for average values. The stress intensity factors and load-line displacements are presented and compared with experimental solutions from the literature. The stress intensity factors and load-line displacements were about 2.5 and 5 percent lower than the reported experimental values, respectively.

Raju, I. S.↗

NASA NTRS1984

Methods for analysis of cracks in three-dimensional solids

Various analytical and numerical methods used to evaluate the stress intensity factors for cracks in three-dimensional (3-D) solids are reviewed. Classical exact solutions and many of the approximate methods used in 3-D analyses of cracks are reviewed. The exact solutions for embedded elliptic cracks in infinite solids are discussed. The approximate methods reviewed are the finite element methods, the boundary integral equation (BIE) method, the mixed methods (superposition of analytical and finite element method, stress difference method, discretization-error method, alternating method, finite element-alternating method), and the line-spring model. The finite element method with singularity elements is the most widely used method. The BIE method only needs modeling of the surfaces of the solid and so is gaining popularity. The line-spring model appears to be the quickest way to obtain good estimates of the stress intensity factors. The finite element-alternating method appears to yield the most accurate solution at the minimum cost.

Raju, I. S.↗

NASA NTRS1984

Prediction of fatigue crack-growth patterns and lives in three-dimensional cracked bodies

Fatigue crack growth patterns and lives for surface cracks, surface cracks at holes, and corner cracks at holes in three dimensional bodies were predicted using linear-elastic fracture mechanics concepts that were modified to account for crack-closure behavior. The predictions were made by using stress intensity factor equations for these crack configurations and the fatigue crack-growth (delta K against rate) relationship for the material of interest. The crack configurations were subjected to constant-amplitude fatigue loading under either remote tension or bending loads. The predicted crack growth patterns and crack growth lives for aluminum alloys agreed well with test data from the literature.

Newman, J. C., Jr.↗

NASA NTRS1984

Stress-intensity factor equations for cracks in three-dimensional finite bodies subjected to tension and bending loads

Stress intensity factor equations are presented for an embedded elliptical crack, a semielliptical surface crack, a quarter elliptical corner crack, a semielliptical surface crack along the bore of a circular hole, and a quarter elliptical corner crack at the edge of a circular hole in finite plates. The plates were subjected to either remote tension or bending loads. The stress intensity factors used to develop these equations were obtained from previous three dimensional finite element analyses of these crack configurations. The equations give stress intensity factors as a function of parametric angle, crack depth, crack length, plate thickness, and, where applicable, hole radius. The ratio of crack depth to plate thickness ranged from 0 to 1, the ratio of crack depth to crack length ranged from 0.2 to 2, and the ratio of hole radius to plate thickness ranged from 0.5 to 2. The effects of plate width on stress intensity variation along the crack front were also included.

Newman, J. C., Jr.↗

NASA NTRS1984

Analysis of interlaminar stresses in thick composite laminates with and without edge delamination

The effect of laminate thickness on the interlaminar stresses in rectangular quasi-isotropic laminates under uniform axial strain was studied. Laminates from 8-ply to infinitely thick were analyzed. Thick laminates were synthesized by stacking (45/0/-45/90) ply groups, rather than grouping like plies. Laminates with and without delaminations were studied. In laminates without delaminations, the free-edge interlaminar normal stress distribution in the outer ply groups was insensitive to total laminate thickness. The interlaminar normal stress distribution for the interior ply groups was nearly the same as for an infinitely thick laminate. In contrast, the free-edge inter-laminar shear stress distribution was nearly the same for inner and outer ply groups and was insensitive to laminate thickness. In laminates with delaminations those delaminations near the top and bottom surfaces of a thick laminate have much larger total strain-energy-release rates (G sub t) and mode I-to-total (G sub t/G sub t) ratios than delaminations deep in the interior. Therefore, delaminations can be expected to grow more easily near the surfaces of a laminate than in the interior.

Whitcomb, J. D.↗

NASA NTRS1983

Superposition method for analysis of free-edge stresses

Superposition techniques were used to transform the edge stress problem for composite laminates into a more lucid form. By eliminating loads and stresses not contributing to interlaminar stresses, the essential aspects of the edge stress problem are easily recognized. Transformed problem statements were developed for both mechanical and thermal loads. Also, a technique for approximate analysis using a two dimensional plane strain analysis was developed. Conventional quasi-three dimensional analysis was used to evaluate the accuracy of the transformed problems and the approximate two dimensional analysis. The transformed problems were shown to be exactly equivalent to the original problems. The approximate two dimensional analysis was found to predict the interlaminar normal and shear stresses reasonably well. Previously announced in STAR as N83-32846

Whitcomb, J. D.↗