Engineering topics

Shivakumar, K. N.

Publications and source records attributed to Shivakumar, K. N..

At least 37 records · Page 2

NASA NTRS1990

Treatment of singularities in cracked bodies

Three-dimensional finite-element analyses of middle-crack tension (M-T) and bend specimens subjected to mode I loadings were performed to study the stress singularity along the crack front. The specimen was modeled using 20-node isoparametric elements. The displacements and stresses from the analysis were used to estimate the power of singularities using a log-log regression analysis along the crack front. The analyses showed that finite-sized cracked bodies have two singular stress fields of the form rho = C sub o (theta, z) r to the -1/2 power + D sub o (theta, phi) R to the lambda rho power. The first term is the cylindrical singularity with the power -1/2 and is dominant over the middle 96 pct (for Poisson's ratio = 0.3) of the crack front and becomes nearly zero at the free surface. The second singularity is a vertex singularity with the vertex point located at the intersection of the crack front and the free surface. The second term is dominant at the free surface and becomes nearly zero away from the boundary layer. The thickness of the boundary layer depends on Poisson's ratio of the material and is independent of the specimen type. The thickness of the boundary layer varied from 0 pct to about 5 pct of the total specimen thickness as Poisson's ratio varied from 0.0 to 0.45. Because there are two singular stress fields near the free surface, the strain energy release rate (G) is an appropriate parameter to measure the severity of the crack.

Shivakumar, K. N.↗

NASA NTRS1990

Evaluation of finite-element models and stress-intensity factors for surface cracks emanating from stress concentrations

Tan, P. W.↗

NASA NTRS1990

Treatment of singularities in a middle-crack tension specimen

A three-dimensional finite-element analysis of a middle-crack tension specimen subjected to mode I loading was performed to study the stress singularity along the crack front. The specimen was modeled using 20-node isoparametric elements with collapsed nonsingular elements at the crack front. The displacements and stresses from the analysis were used to estimate the power of singularities, by a log-log regression analysis, along the crack front. Analyses showed that finite-sized cracked bodies have two singular stress fields. Because of two singular stress fields near the free surface and the classical square root singularity elsewhere, the strain energy release rate appears to be an appropriate parameter all along the crack front.

Shivakumar, K. N.↗

NASA NTRS1990

An equivalent domain integral method in the two-dimensional analysis of mixed mode crack problems

An equivalent domain integral (EDI) method for calculating J-integrals for two-dimensional cracked elastic bodies is presented. The details of the method and its implementation are presented for isoparametric elements. The EDI method gave accurate values of the J-integrals for two mode I and two mixed mode problems. Numerical studies showed that domains consisting of one layer of elements are sufficient to obtain accurate J-integral values. Two procedures for separating the individual modes from the domain integrals are presented.

Raju, I. S.↗

NASA NTRS1989

Effects of T-tabs and large deflections in DCB specimen tests

A simple strength of materials analysis was developed for a double-cantilever beam (DCB) specimen to account for geometric nonlinearity effects due to large deflections and T-tabs. A new DCB data analysis procedure was developed to include the effects of these nonlinearities. The results of the analysis were evaluated by DCB tests performed for materials having a wide range of toughnesses. The materials used in the present study were T300/5208, IM7/8551-7, and AS4/PEEK. Based on the present analysis, for a typical deflection/crack length ratio of 0.3 (for AS4/PEEK), T-tabs and large deflections cause a 15 percent and 3 percent error, respectively, in the computer Mode 1 strain energy release rate. Design guidelines for DCB specimen thickness and T-tab height were also developed in order to keep errors due to these nonlinearities within 2 percent. Based on the test results, for both hinged and tabbed specimens, the effects of large deflection on the Mode 1 fracture toughness (G sub Ic) were almost negligible (less than 1 percent) in the case of T300/5208 and IM7/8551-7; however, AS4/PEEK showed a 2 to 3 percent effect. The effects of T-tabs G sub Ic were more significant for all the materials with T300/5208 showing a 5 percent error, IM7/8551-7 a 15 percent error, and, AS4/PEEK a 20 percent error.

Naik, R. A.↗

NASA NTRS1989

Strain-energy release rate analysis of plates with postbuckled delaminations

The objectives are to present the derivation of the new virtual crack closure technique, evaluate the accuracy of the technique, and finally to present the results of a limited parametric study of laminates with a postbuckled delamination. Although the new virtual crack closure technique is general, only homogeneous, isotropic laminates were analyzed. This was to eliminate the variation of flexural stiffness with orientation, which occurs even for quasi-isotropic laminates. This made it easier to identify the effect of geometrical parameters on G. The new virtual crack closure technique is derived. Then the specimen configurations are described. Next, the stress analyses is discussed. Finally, the virtual crack closure technique is evaluated and then used to calculate the distribution of G along the delamination front of several laminates with a postbuckled delamination.

Whitcomb, John D.↗

NASA NTRS1989

Treatment of singularities in cracked bodies

Three-dimensional finite-element analyses of middle-crack tension (M-T) and bend specimens subjected to mode I loadings were performed to study the stress singularity along the crack front. The specimen was modeled using 20-node isoparametric elements. The displacements and stresses from the analysis were used to estimate the power of singularities using a log-log regression analysis along the crack front. The analyses showed that finite-sized cracked bodies have two singular stress fields of the form rho = C sub o (theta, z) r to the -1/2 power + D sub o (theta, phi) R to the lambda rho power. The first term is the cylindrical singularity with the power -1/2 and is dominant over the middle 96 pct (for Poisson's ratio = 0.3) of the crack front and becomes nearly zero at the free surface. The second singularity is a vertex singularity with the vertex point located at the intersection of the crack front and the free surface. The second term is dominant at the free surface and becomes nearly zero away from the the boundary layer. The thickness of the boundary layer depends on Poisson's ratio of the material and is independent of the specimen type. The thickness of the boundary layer varied from 0 pct to about 5 pct of the total specimen thickness as Poisson's ratio varied from 0.0 to 0.45. Because there are two singular stress fields near the free surface, the strain energy release rate (G) is an appropriate parameter to measure the severity of the crack.

Shivakumar, K. N.↗

NASA NTRS1989

Implementation of equivalent domain integral method in the two-dimensional analysis of mixed mode problems

An equivalent domain integral (EDI) method for calculating J-intergrals for two-dimensional cracked elastic bodies is presented. The details of the method and its implementation are presented for isoparametric elements. The total and product integrals consist of the sum of an area of domain integral and line integrals on the crack faces. The line integrals vanish only when the crack faces are traction free and the loading is either pure mode 1 or pure mode 2 or a combination of both with only the square-root singular term in the stress field. The EDI method gave accurate values of the J-integrals for two mode I and two mixed mode problems. Numerical studies showed that domains consisting of one layer of elements are sufficient to obtain accurate J-integral values. Two procedures for separating the individual modes from the domain integrals are presented. The procedure that uses the symmetric and antisymmetric components of the stress and displacement fields to calculate the individual modes gave accurate values of the integrals for all problems analyzed. The EDI method when applied to a problem of an interface crack in two different materials showed that the mode 1 and mode 2 components are domain dependent while the total integral is not. This behavior is caused by the presence of the oscillatory part of the singularity in bimaterial crack problems. The EDI method, thus, shows behavior similar to the virtual crack closure method for bimaterial problems.

Raju, I. S.↗

NASA NTRS1989

Numerical fracture simulation of bend specimens using a CTOD criterion

Stable crack growth in 50-mm-thick side-grooved three-point bend specimens of HY-130 steel is modeled numerically. A two-dimensional elastoplastic FEM formulation incorporating a CTOD fracture criterion is developed analytically, and the theoretical predictions are compared with experimental data in graphs. Good general agreement is obtained, except for the crack extension, which was typically overpredicted. The value of the J integral at crack initiation is found to be 83 kJ/sq m.

Shivakumar, K. N.↗

NASA NTRS1989

An equivalent domain integral for analysis of two-dimensional mixed mode problems

An equivalent domain integral (EDI) method for calculating J-integrals for two-dimensional cracked elastic bodies subjected to mixed mode loading is presented. The total and product integrals consist of the sum of an area or domain integral and line integrals on the crack faces. The EDI method gave accurate values of the J-integrals for two mode I and two mixed mode problems. Numerical studies showed that domains consisting of one layer of elements are sufficient to obtain accurate J-integral values. Two procedures for separating the individual modes from the domain integrals are presented. The procedure that uses the symmetric and antisymmetric components of the stress and displacement fields to calculate the individual modes gave accurate values of the integrals for all the problems analyzed.

Raju, I. S.↗

NASA NTRS1989

Three-dimensional aspects of plasticity-induced fatigue crack closure

A three-dimensional elastoplastic finite element analysis of fatigue crack growth and closure in a center-crack tension specimen was performed in order to investigate the crack opening stress as a function of specimen thickness under constant amplitude loading conditions. Elastic-perfectly plastic conditions were assumed and cycle loading was kept at R = 0.1. Firstly, by imposing proper boundary conditions on the three-dimensional model, a plane strain analysis was carried out. The ratio of stabilized opening stress over the maximum applied stress under plane strain conditions was found to be 0.28. Secondly, the model thickness was varied under fully three-dimensional conditions, and crack-opening stresses were determined for each specimen thickness. The crack-opening stress is found to vary through the thickness for the center crack tension specimen. On the specimen surface and in the midplane the crack-opening stress levels tend to two-dimensional solutions for plane stress and plane strain conditions, respectively.

Chermahini, R. G.↗

NASA NTRS1988

Three-dimensional elastic analysis of a composite double cantilever beam specimen

Attention is given to the stresses and the strain energy release rate along the delamination front in the present three-dimensional elastic analysis of a 24-ply, cocured double-cantilever beam specimen by means of 20-noded parabolic-isoparametric finite elements. At the free surface, the strain energy release rate was found to be substantially smaller than the plane strain value; this is suggested to be due to the free-surface effect that exists where the delamination meets the surface edge.

Raju, I. S.↗

NASA NTRS1988

A re-evaluation of finite-element models and stress-intensity factors for surface cracks emanating from stress concentrations

A re-evaluation of the 3-D finite-element models and methods used to analyze surface crack at stress concentrations is presented. Previous finite-element models used by Raju and Newman for surface and corner cracks at holes were shown to have ill-shaped elements at the intersection of the hole and crack boundaries. These ill-shaped elements tended to make the model too stiff and, hence, gave lower stress-intensity factors near the hole-crack intersection than models without these elements. Improved models, without these ill-shaped elements, were developed for a surface crack at a circular hole and at a semi-circular edge notch. Stress-intensity factors were calculated by both the nodal-force and virtual-crack-closure methods. Both methods and different models gave essentially the same results. Comparisons made between the previously developed stress-intensity factor equations and the results from the improved models agreed well except for configurations with large notch-radii-to-plate-thickness ratios. Stress-intensity factors for a semi-elliptical surface crack located at the center of a semi-circular edge notch in a plate subjected to remote tensile loadings were calculated using the improved models. The ratio of crack depth to crack length ranged form 0.4 to 2; the ratio of crack depth to plate thickness ranged from 0.2 to 0.8; and the ratio of notch radius to the plate thickness ranged from 1 to 3. The models had about 15,000 degrees-of-freedom. Stress-intensity factors were calculated by using the nodal-force method.

Tan, P. W.↗

NASA NTRS1988

A virtual crack-closure technique for calculating stress intensity factors for cracked three dimensional bodies

A three-dimensional virtual crack-closure technique is presented which calculates the strain energy release rates and the stress intensity factors using only nodal forces and displacements from a standard finite element analysis. The technique is an extension of the Rybicki-Kanninen (1977) method, and it assumes that any continuous function can be approximated by a finite number of straight line segments. Results obtained by the method for surface cracked plates with and without notches agree favorably with previous results.

Shivakumar, K. N.↗

NASA NTRS1988

A fiber-resin micromechanics analysis of the delamination front in a DCB specimen

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.

Crews, J. H.↗

NASA NTRS1988

An elastic-plastic finite-element analysis of the J-resistance curve using a CTOD criterion

The fracture process in compact and bend specimens was simulated using a two-dimensional finite-element analysis of the J-resistance (J-R) curve and a CTOD criterion. The J-R curves were calculated from the numerical results for each specimen type using several different methods. In general, the J-R curves obtained for the bend specimens were found to be higher than those for the compact specimens, especially beyond maximum load. However, below the maximum load, the modified deformation theory of plasticity and the contour-integral J(r) method resulted in very similar J-R curves for both specimen types.

Newman, J. C., Jr.↗

NASA NTRS1988

Factors influencing elastic stresses in double cantilever beam specimens

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.

Crews, J. H., Jr.↗

NASA NTRS1988

Three-dimensional finite-element simulation of fatigue crack growth and closure

A three-dimensional elastic-plastic finite-element analysis of crack growth and closure under cyclic loading has been performed in order to investigate the behavior of a crack in a finite-thickness middle-crack tension specimen. The cases of both constant-amplitude loading and a single-spike overload are considered. The calculated crack-opening stresses for the exterior of the specimen agree with previous plane-stress results, and those obtained for the interior of the specimen agree with previous plane-strain results.

Chermahini, R. G.↗