Engineering Papers⌕ Search

Engineering topics

Shivakumar, K. N.

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

54 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Numerical fracture simulation of compact and bend specimens

An elastic-plastic finite-element analysis with a critical crack-tip-opening displacement criterion was used to simulate fracture of various size compact and bend specimens made of HY-130 steel. From the calculated load-crack-extension and load-displacement curves, J-resistance (J-R) curves were determined by several methods. The simulated 3-R curves were insensitive to specimen size up to maximum load but were sensitive to specimen configuration for crack extensions greater than 10 percent of the initial uncracked ligament length.

Shivakumar, K. N.↗
NASA NTRS1987

Strain-energy release rate analysis of a laminate with a postbuckled delamination

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 NTRS1987

Energy dissipation associated with crack extension in an elastic-plastic material

Crack extension in elastic-plastic material involves energy dissipation through the creation of new crack surfaces and additional yielding around the crack front. An analytical procedure, using a two-dimensional elastic-plastic finite element method, was developed to calculate the energy dissipation components during a quasi-static crack extension. The fracture of an isotropic compact specimen was numerically simulated using the critical crack-tip-opening-displacement (CTOD) growth criterion. Two specimen sizes were analyzed for three values of critical CTOD. Results from the analyses showed that the total energy dissipation rate consisted of three components: the crack separation energy rate, the plastic energy dissipation rate, and the residual strain energy rate. All three energy dissipation components and the total energy dissipation rate initially increased with crack extension and finally reached constant values.

Shivakumar, K. N.↗
NASA NTRS1987

Strain-energy release rate analysis of a laminate with a postbuckled delamination

Delamination growth due to local buckling of a delamination was studied. Delamination growth was assumed to be related to the total strain-energy release rate, G. In order to calculate the distribution of G along the delamination front, a new virtual crack closure technique was developed which is suitable for use with plate analysis. A check of the technique indicated that it is accurate. For square and rectangular delaminations, there is a large variation of G along the delamination front. Hence, self-similar growth is not expected. Whether a delamination grows in the load direction or perpendicular to the load direction was found to depend on the current delamination aspect ratio, the strain level, and the absolute size of the delamination.

Whitcomb, John D.↗
NASA NTRS1987

Energy dissipation associated with crack extension in an elastic-plastic material

Crack extension in elastic-plastic material involves energy dissipation through the creation of new crack surfaces and additional yielding around the crack front. An analytical procedure, using a two-dimensional elastic-plastic finite element method, was developed to calculate the energy dissipation components during a quasi-static crack extension. The fracture of an isotropic compact specimen was numerically simulated using the critical crack-tip-opening-displacement (CTOD) growth criterion. Two specimen sizes were analyzed for three values of critical CTOD. Results from the analysis showed that the total energy dissipation rate consisted of three components: the crack separation energy rate, the plastic energy dissipation rate, and the residual strain energy rate. All three energy dissipation components and the total energy dissipation rate initially increased with crack extension and finally reached constant values.

Shivakumar, K. N.↗
NASA NTRS1986

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 NTRS1985

Prediction of impact force and duration due to low-velocity impact on circular composite laminates

Two simple and improved models - energy-balance and spring-mass - were developed to calculate impact force and duration during low velocity impact of circular composite plates. Both models include the contact deformation of the plate and the impactor as well as bending, transverse shear, and membrane deformations of the plate. The plate was transversely isotropic graphite/epoxy composite laminate and the impactor was a steel sphere. Calculated impact forces from the two analyses agreed with each other. The analyses were verified by comparing the results with reported test data.

Shivakumar, K. N.↗
NASA NTRS1985

Buckling of a sublaminate in a quasi-isotropic composite laminate

The buckling of an elliptic delamination embedded near the surface of a thick quasi-isotropic laminate was predicted. The thickness of the delaminated ply group (the sublaminate) was assumed to be small compared to the total laminate thickness. Finite-element and Rayleigh-Ritz methods were used for the analyses. The Rayleigh-Ritz method was found to be simple, inexpensive, and accurate, except for highly anisotropic delaminated regions. Effects of delamination shape and orientation, material anisotropy, and layup on buckling strains were examined. Results show that: (1) the stress state around the delaminated region is biaxial, which may lead to buckling when the laminate is loaded in tension; (2) buckling strains for multi-directional fiber sublaminates generally are bounded by those for the 0 deg and 90 deg unidirectional sublaminates; and (3) the direction of elongation of the sublaminate that has the lowest buckling strain correlates with the delamination growth direction.

Shivakumar, K. N.↗
NASA NTRS1984

Delamination growth analysis in quasi-isotropic laminates under loads simulating low-velocity impact

A geometrically nonlinear finite-element analysis was developed to calculate the strain energy released by delamination plates during impact loading. Only the first mode of deformation, which is equivalent to static deflection, was treated. Both the impact loading and delamination in the plate were assumed to be axisymmetric. The strain energy release rate in peeling, G sub I, and shear sliding, G sub II, modes were calculated using the fracture mechanics crack closure technique. Energy release rates for various delamination sizes and locations and for various plate configurations and materials were compared. The analysis indicated that shear sliding (G sub II) was the primary mode of delamination growth. The analysis also indicated that the midplane (maximum transverse shear stress plane) delamination was more critical and would grow before any other delamination of the same size near the midplane region. The delamination growth rate was higher (neutrally stable) for a low toughness (brittle) matrix and slower (stable) for high toughness matrix. The energy release rate in the peeling mode, G sub I, for a near-surface delamination can be as high as 0.5G sub II and can contribute significantly to the delamination growth.

Shivakumar, K. N.↗
NASA NTRS1984

Delamination growth analysis in quasi-isotropic laminates under loads simulating low-velocity impact

A geometrically nonlinear finite-element analysis has been developed to calculate the strain energy released by delaminating plates during impact loading. Only the first mode of deformation, which is equivalent to static deflection, was treated. Both the impact loading and delamination in the plate were assumed to be axisymmetric. The strain energy release rate in peeling, GI, and shear sliding, GII, modes were calculated using the fracture mechanics crack closure technique. Energy release rates for various delamination sizes and locations and for various plate configurations and materials were compared. The analysis indicated that shear sliding was the primary mode of delamination growth. The analysis also indicated that the midplane (maximum transverse shear stress plane) delamination was more critical and would grow first before any other delamination of the same size near the midplane region. The delamination growth rate was higher (neutrally stable) for a low toughness (brittle) matrix and slower (stable) for high toughness matrix. The energy release rate in the peeling mode, GI, for a near-surface delamination can be as high as 0.5GII, and can contribute significantly to the delamination growth.

Shivakumar, K. N.↗
NASA NTRS1984

Buckling of a sublaminate in a quasi-isotropic composite laminate

The buckling of an elliptic delamination embedded near the surface of a thick quasi-isotropic laminate was predicted. The thickness of the delaminated ply group (the sublaminate) was assumed to be small compared to the total laminate thickness. Finite-element and Rayleigh-Ritz methods were used for the analyses. The Rayleigh-Ritz method was found to be simple, inexpensive, and accurate, except for highly anisotropic delaminated regions. Effects of delamination shape and orientation, material anisotropy, and layup on buckling strains were examined. Results show that: (1) the stress state around the delaminated region is biaxial, which may lead to buckling when the laminate is loaded in tension; (2) buckling strains for multi-directional fiber sublaminates generally are bounded by those for the 0 deg and 90 deg unidirectional sublaminates; and (3) the direction of elongation of the sublaminate that has the lowest buckling strain correlates with the delamination growth direction.

Shivakumar, K. N.↗
NASA NTRS1983

Prediction of impact force and duration during low velocity impact on circular composite laminates

Two simple and improved models--energy-balance and spring-mass--were developed to calculate impact force and duration during low velocity impact of circular composite plates. Both models include the contact deformation of the plate and the impactor as well as bending, transverse shear, and membrane deformations of the plate. The plate was transversely isotropic graphite/epoxy composite laminate and the impactor was a steel sphere. Calculated impact forces from the two analyses agreed with each other. The analyses were verified by comparing the results with reported test data.

Shivakumar, K. N.↗
NASA NTRS1983

Analysis of progressive damage in thin circular laminates due to static-equivalent impact loads

Clamped circular graphite/epoxy plates (25.4, 38.1, and 50.8 mm radii) with an 8-ply quasi-isotropic layup were analyzed for static-equivalent impact loads using the minimum-total-potential-energy method and the von Karman strain-displacement equations. A step-by-step incremental transverse displacement procedure was used to calculate plate load and ply stresses. The ply failure region was calculated using the Tsai-Wu criterion. The corresponding failure modes (splitting and fiber failure) were determined using the maximum stress criteria. The first-failure mode was splitting and initiated first in the bottom ply. The splitting-failure thresholds were relatively low and tended to be lower for larger plates than for small plates. The splitting-damage region in each ply was elongated in its fiber direction; the bottom ply had the largest damage region. The calculated damage region for the 25.4-mm-radius plate agreed with limited static test results from the literature.

Shivakumar, K. N.↗
NASA NTRS1983

Bolt clampup relaxation in a graphite/epoxy laminate

A simple bolted joint was analyzed to calculate bolt clampup relaxation for a graphite/epoxy (T300/5208) laminate. A viscoelastic finite element analysis of a double-lap joint with a steel bolt was conducted. Clampup forces were calculated for various steady-state temperature-moisture conditions using a 20-year exposure duration. The finite element analysis predicted that clampup forces relax even for the room-temperature-dry condition. The relaxations were 8, 13, 20, and 30 percent for exposure durations of 1 day, 1 month, 1 year, and 20 years, respectively. As expected, higher temperatures and moisture levels each increased the relaxation rate. The combined viscoelastic effects of steady-state temperature and moisture appeared to be additive. From the finite-element analysis, a simple equation was used to calculate clampup forces for the same temperature-moisture conditions as used in the finite-element analysis. The two sets of calculated results agreed well. Previously announced in STAR as N82-19320

Shivakumar, K. N.↗
NASA NTRS1982

An equation for bolt clampup relaxation in transient environments

An equation for bolt clampup relaxation for transient temperature-moisture (T-M) conditions was derived starting with a relaxation equation for steady-state conditions, and then using an incremental time approach that exploits the superposition principle for linear viscoelasticity. The resulting equation uses the initial T-M condition (at the time of clamping), the T-M history after clamping, and elastic clampup coefficients for temperature and moisture changes. For a given material and joint configuration, the clampup coefficients are constants that can be calculated by elastic analyses. The clampup equation was used to calculate the changes in clampup occurring in a T300/5208 graphite/epoxy joint exposed to a one-year history of temperature and moisture. Two cases were considered: one was a dry joint exposed to a relatively humid environment and the other was a nearly saturated joint exposed to an arid environment.

Shivakumar, K. N.↗
NASA NTRS1982

Bolt clampup relaxation in a graphite/epoxy laminate

A simple bolted joint was analyzed to calculate bolt clampup relaxation for a graphite/epoxy (T300/5208) laminate. A viscoelastic finite element analysis of a double-lap joint with a steel bolt was conducted. Clampup forces were calculated for various steady-state temperature-moisture conditions using a 20-year exposure duration. The finite element analysis predicted that clampup forces relax even for the room-temperature-dry condition. The relaxations were 8, 13, 20, and 30 percent for exposure durations of 1 day, 1 month, 1 year, and 20 years, respectively. As expected, higher temperatures and moisture levels each increased the relaxation rate. The combined viscoelastic effects of steady-state temperature and moisture appeared to be additive. From the finite-element analysis, a simple equation was developed for clampup force relaxation. This generalized equation was used to calculate clampup forces for the same temperature-moisture conditions as used in the finite-element analysis. The two sets of calculated results agreed well.

Shivakumar, K. N.↗