Engineering topics
Publications and source records attributed to Raju, I. S..
Q3DG computer program developed to perform quasi-three-dimensional stress analysis of composite laminates containing delaminations. Calculates strain-energy-release rates for long, rectangular composite laminates containing delaminations and subjected to any combination of mechanical, thermal, and hygroscopic loading. Written in FORTRAN V.
A 3-D finite element model was developed to analyze the fiber-resin behavior near the delamination front in a graphite-epoxy double cantilever beam (DCB) specimen. The specimen interior was analyzed using a typical one-fiber slice, represented by a local 3-D fiber-resin model. The resin stresses were computed for the resin-rich layer at the ply interface as well as for the regions between the fibers close to the delamination front. However, the computed strain energy release rate G sub I along the delamination front varied by less than two percent, and was within about four percent of the plane-strain value. The von Mises yield criterion was used to estimate the extent of yielding near the delamination front. The yielding extended ahead of the delamination and also developed between the fibers. Although the fibers had only a negligible effect on G sub I, they caused yielding within the ply and therefore could influence delamination fracture toughness. The normal and shear stresses at the fiber-resin interface were computed near the delamination front. These results suggest that multi-axial stress criteria may be required to analyze fiber-resin interfaces.
Strain energy release rates for edge delaminated composite laminates were obtained using quasi 3 dimensional finite element analysis. The problem of edge delamination at the -35/90 interfaces of an 8-ply composite laminate subjected to uniform axial strain was studied. The individual components of the strain energy release rates did not show convergence as the delamination tip elements were made smaller. In contrast, the total strain energy release rate converged and remained unchanged as the delamination tip elements were made smaller and agreed with that calculated using a classical laminated plate theory. The studies of the near field solutions for a delamination at an interface between two dissimilar isotropic or orthotropic plates showed that the imaginary part of the singularity is the cause of the nonconvergent behavior of the individual components. To evaluate the accuracy of the results, an 8-ply laminate with the delamination modeled in a thin resin layer, that exists between the -35 and 90 plies, was analyzed. Because the delamination exists in a homogeneous isotropic material, the oscillatory component of the singularity vanishes.
An elastic stress analysis was conducted for a double cantilever beam (DCB) specimen using finite-element methods. The purpose of this study was to identify the important parameters that influence stresses ahead of the delamination front. The study focused on an aluminum DCB specimen, typical of adhesively-bonded joints, and on a graphite/epoxy specimen representing a cocured composite. Opening mode sigma sub y stresses ahead of the crack tip were calculated and compared with those for a monolithic reference specimen. Beyond the singularity-dominated region very near the crack tip, the sigma sub y distribution was elevated compared to the monolithic case. Both the adhesive thickness and the adherend transverse (thickness-direction) stiffness were found to influence the elevation of sigma sub y. In contrast, adherend thickness and longitudinal stiffness has very little effect on this stress distribution. Estimates for adhesive yielding beyond the aluminum DCB crack tip showed that both the area and height of the plastic zone increased to a peak value for increasing adhesive thicknesses. Results from this study would provide insight for comparing data from different DCB specimens and for designing new DCB specimens.
The boundary-force method developed by Tan (1985) and Tan et al. (1986) is extended to two-dimensional mixed boundary-value problems in fracture mechanics. The formulation is explained in detail, and results are presented in tables and graphs for sample problems involving an edge crack emanating from a semielliptical and a semicircular notch and an edge crack emanating from a V notch. The accuracy of the approach is demonstrated by comparisons with collocation solutions for an edge-clamped center-cracked tension specimen.
Stress intensity factors are presented for a wide range of quarter-elliptical corner cracks in a rectangular bar. The bars were subjected to remote tension, out-of-plane bending, or in-plane bending loading, and the stress intensity factors were calculated from a three-dimensional finite element analysis using the nodal-force method. For all loadings considered, smaller width bars produce higher normalized stress intensity factors than the larger width bars. For tensile loading, quarter-circular corner cracks give higher stress intensity factors near the free surfaces and lower values in the interior. The results are compared with data in the literature for remote tensile loading.
Two-dimensional finite-element analysis of the end-notched flexure specimen was performed using 8-node isoparametric, parabolic elements to evaluate compliance and mode II strain energy release rates, G sub II. The G sub II values were computed using two different techniques: the virtual crack-closure technique (VCCT) and the rate of change of compliance with crack length (compliance derivative method). The analysis was performed for various crack-length-to-semi-span (a/L) ratios ranging from 0.2 to 0.9. Three material systems representing a wide range of material properties were analyzed. The compliance and strain energy release rates of the specimen calculated with the present finite-element analysis agree very well with beam theory equations including transverse shear. The G sub II values calculated using the compliance derivative method compared extremely well with those calculated using the VCCT. The G sub II values obtained by the compliance derivative method using the top or bottom beam deflections agreed closely with each other. The strain energy release rates from a plane-stress analysis were higher than the plane-strain values by only a small percentage, indicating that either assumption may be used in the analysis. The G sub II values for one material system calculated from the finte-element analysis agreed with one solution in the literature and disagreed with the other solution in the literature.
Finite element method was used to analyze a tapered glass/epoxy composite laminate subjected to a tension load, in order to determine the interlaminar stress distributions the strain-energy release rate for the delamination growth that may occur due to ply drops. In a laminate having a typical configuration of a helicopter hub, the (+/- 45)3 plies were dropped in three distinct steps, each 20-ply thickness apart, with the resulting taper angle of 5.71 deg. Delaminations were assumed to initiate at the bottom of the taper on the -45/+45 interface, and the delamination growth was simulated along the taper and into the thin region. The results of the analysis indicated that a delamination initiating at the end of the taper will grow unstably along the taper and the thin laminate simultaneously.
Various methods used for determining stress-intensity factors for cracked three-dimensional bodies are reviewed. The review is limited to the determination of mode-I stress-intensity factors. Some exact solutions for cracks in infinite solids are presented, and various approximate methods that have been used for the solution of the boundary-value problem of finite solids with cracks are addressed. The techniques used to extract the stress-intensity factors from these solutions are considered, and various methods used to calculate the stress-intensity factors for through-the-thickness cracks, semielliptical surface cracks, and quarter-elliptical corner cracks at holes are compared for the case of remote tensile loading.
Two-dimensional finite-element analysis of the end-notched flexure specimen was performed using 8-node isoparametric, parabolic elements to evaluate compliance and mode II strain energy release rates, G sub II. The G sub II values were computed using two different techniques: the virtural crack-closure technique (VCCT) and the rate of change of compliance with crack length (compliance derivative method). The analysis was performed for various crack-length-to-semi-span (a/L) ratios ranging from 0.2 to 0.9. Three material systems representing a wide range of material properties were analyzed. The compliance and strain energy release rates of the specimen calculated with the present finite-element analysis agree very well with beam theory equations including transverse shear. The G sub II values calculated using the compliance derivative method compared extremely well with those calculated using the VCCT. The G sub II values obtained by the compliance derivative method using the top or bottom beam deflections agreed closely with each other. The strain energy release rates from a plane-stress analysis were higher than the plane-strain values by only a small percentage, indicating that either assumption may be used in the analysis. The G sub II values for one material system calculated from the finite-element analysis agreed with one solution in the literature and disagreed with the other solution in the literature.
The Boundary Force Method (BFM) was formulated for the three fundamental problems of elasticity: the stress boundary value problem, the displacement boundary value problem, and the mixed boundary value problem. Because the BFM is a form of an indirect boundary element method, only the boundaries of the region of interest are modeled. The elasticity solution for the stress distribution due to concentrated forces and a moment applied at an arbitrary point in a cracked infinite plate is used as the fundamental solution. Thus, unlike other boundary element methods, here the crack face 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 a center-cracked specimen subjected to mixed boundary conditions and a three-hole cracked configuration subjected to traction boundary conditions. The results obtained are in good agreement with accepted numerical solutions. The method is then used to generate stress-intensity solutions for two common cracked configurations: an edge crack emanating from a semi-elliptical notch, and an edge crack emanating from a V-notch. The BFM is a versatile technique that can be used to obtain very accurate stress intensity factors for complex crack configurations subjected to stress, displacement, or mixed boundary conditions. The method requires a minimal amount of modeling effort.
Strain energy release rates for edge delaminated composite laminates were obtained using quasi 3 dimensional finite element analysis. The problem of edge delamination at the -35/90 interfaces of an 8-ply composite laminate subjected to uniform axial strain was studied. The individual components of the strain energy release rates did not show convergence as the delamination tip elements were made smaller. In contrast, the total strain energy release rate converged and remained unchanged as the delamination tip elements were made smaller and agreed with that calculated using a classical laminated plate theory. The studies of the near field solutions for a delamination at an interface between two dissimilar isotropic or orthotropic plates showed that the imaginary part of the singularity is the cause of the nonconvergent behavior of the individual components. To evaluate the accuracy of the results, an 8-ply laminate with the delamination modeled in a thin resin layer, that exists between the -35 and 90 plies, was analyzed. Because the delamination exists in a homogeneous isotropic material, the oscillatory component of the singularity vanishes.
Explore the source record for details and available documents.
A general finite element procedure for obtaining strain-energy release rates for crack growth in isotropic materials is presented. The procedure is applicable to two-dimensional finite element analyses and uses the virtual crack-closure method. The procedure was applied to nonsingular 4-noded (linear), 8-noded (parabolic), and 12-noded (cubic) elements and to quarter-point and cubic singularity elements. Simple formulas for strain-energy release rates were obtained with this procedure for both nonsingular and singularity elements. The formulas were evaluated by applying them to two mode I and two mixed mode problems. Comparisons with results from the literature for these problems showed that the formulas give accurate strain-energy release rates.
A general finite element procedure for obtaining strain-energy release rates for crack growth in isotropic materials is presented. The procedure is applicable to two-dimensional finite element analyses and uses the virtual crack-closure method. The procedure was applied to non-singular 4-noded (linear), 8-noded (parabolic), and 12-noded (cubic) elements and to quarter-point and cubic singularity elements. Simple formulas for strain-energy release rates were obtained with this procedure for both non-singular and singularity elements. The formulas were evaluated by applying them to two mode I and two mixed mode problems. Comparisons with results from the literature for these problems showed that the formulas give accurate strain-energy release rates.
An elastic stress analysis was conducted for a double cantilever beam (DCB) specimen using finite-element methods. The purpose of this study was to identify the important parameters that influence stresses ahead of the delamination front. The study focused on an aluminum DCB specimen, typical of adhesively-bonded joints, and on a graphite/epoxy specimen representing a cocured composite. Opening mode sigma sub y stresses ahead of the crack tip were calculated and compared with those for a monolithic reference specimen. Beyond the singularity-dominated region very near the crack tip, the sigma sub y distribution was elevated compared to the monolithic case. Both the adhesive thickness and the adherend transverse (thickness-direction) stiffness were found to influence the elevation of sigma sub y. In contrast, adherend thickness and longitudinal stiffness has very little effect on this stress distribution. Estimates for adhesive yielding beyond the aluminum DCB crack tip showed that both the area and height of the plastic zone increased to a peak value for increasing adhesive thicknesses. Results from this study would provide insight for comparing data from different DCB specimens and for designing new DCB specimens.
The Q3DG is a computer program developed to perform a quasi-three-dimensional stress analysis for composite laminates which may contain delaminations. The laminates may be subjected to mechanical, thermal, and hygroscopic loads. The program uses the finite element method and models the laminates with eight-noded parabolic isoparametric elements. The program computes the strain-energy-release components and the total strain-energy release in all three modes for delamination growth. A rectangular mesh and data file generator, DATGEN, is included. The DATGEN program can be executed interactively and is user friendly. The documentation includes sections dealing with the Q3D analysis theory, derivation of element stiffness matrices and consistent load vectors for the parabolic element. Several sample problems with the input for Q3DG and output from the program are included. The capabilities of the DATGEN program are illustrated with examples of interactive sessions. A microfiche of all the examples is included. The Q3DG and DATGEN programs have been implemented on CYBER 170 class computers. Q3DG and DATGEN were developed at the Langley Research Center during the early eighties and documented in 1984 to 1985.
Although the application of multigrid methods to the equations of elasticity has been suggested, few such applications have been reported in the literature. In the present work, multigrid techniques are applied to the finite element analysis of a simply supported Bernoulli-Euler beam, and various aspects of the multigrid algorithm are studied and explained in detail. In this study, six grid levels were used to model half the beam. With linear prolongation and sequential ordering, the multigrid algorithm yielded results which were of machine accuracy with work equivalent to 200 standard Gauss-Seidel iterations on the fine grid. Also with linear prolongation and sequential ordering, the V(1,n) cycle with n greater than 2 yielded better convergence rates than the V(n,1) cycle. The restriction and prolongation operators were derived based on energy principles. Conserving energy during the inter-grid transfers required that the prolongation operator be the transpose of the restriction operator, and led to improved convergence rates. With energy-conserving prolongation and sequential ordering, the multigrid algorithm yielded results of machine accuracy with a work equivalent to 45 Gauss-Seidel iterations on the fine grid. The red-black ordering of relaxations yielded solutions of machine accuracy in a single V(1,1) cycle, which required work equivalent to about 4 iterations on the finest grid level.