Engineering PapersSearch

Engineering topics

Shivakumar, K. N.

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

At least 19 records

Modified Mode-I Cracked Sandwich Beam (CSB) Fracture Test

Five composite sandwich panels were fabricated using vacuum assisted resin transfer molding (VARTM). Four of these panels had E-glass/vinylester facesheets and one had carbon/epoxy facesheets. The sandwich panels had different density PVC foam cores. The four E-glass panels had core densities of 80, 100, 130, 200 kg/cu m. The sandwich with carbon/epoxy 3 facesheets had a core with density of 100 kg/cu m. Fracture tests were conducted using a modified Cracked Sandwich Beam (CSB) test configuration. Load displacement curves were obtained for loading and unloading of the specimens during crack growth. Various increments of crack growth were monitored. Critical Strain Energy Release Rates (SERR) were determined from the tests using the area method. The critical values of SERR can be considered the fracture toughness of the sandwich material. The fracture toughness ranged 367 J/sq m to 1350 J/sq m over the range of core densities. These results are compared to the Mode-I fracture toughness of the PVC foam core materials and values obtained for foam-cored sandwiches using the TSD specimen. Finite-element analyses (FEA) were performed for the test configuration and Strain Energy Release Rates were calculated using the Virtual Crack Closure Technique (VCCT). The SERR values determined from the FEA were scaled to the fracture loads, or critical loads, obtained from the modified CSB tests. These critical loads were in close agreement with the test values.

Smith, S. A.

Aero-Thermo-Structural Analysis of Inlet for Rocket Based Combined Cycle Engines

NASA has been developing advanced space transportation concepts and technologies to make access to space less costly. One such concept is the reusable vehicles with short turn-around times. The NASA Glenn Research Center's concept vehicle is the Trailblazer powered by a rocket-based combined cycle (RBCC) engine. Inlet is one of the most important components of the RBCC engine. This paper presents fluid flow, thermal, and structural analysis of the inlet for Mach 6 free stream velocity for fully supersonic and supercritical with backpressure conditions. The results concluded that the fully supersonic condition was the most severe case and the largest stresses occur in the ceramic matrix composite layer of the inlet cowl. The maximum tensile and the compressive stresses were at least 3.8 and 3.4, respectively, times less than the associated material strength.

Shivakumar, K. N.

Analysis of surface cracks at hole by a 3-D weight function method with stresses from finite element method

Parallel with the work in Part-1, stress intensity factors for semi-elliptical surface cracks emanating from a circular hole are determined. The 3-D weight function method with the 3D finite element solutions for the uncracked stress distribution as in Part-1 is used for the analysis. Two different loading conditions, i.e. remote tension and wedge loading, are considered for a wide range in geometrical parameters. Both single and double surface cracks are studied and compared with other solutions available in the literature. Typical crack opening displacements are also provided.

Zhao, W.

Analysis of corner cracks at hole by a 3-D weight function method with stresses from finite element method

Stress intensity factors for quarter-elliptical corner cracks emanating from a circular hole are determined using a 3-D weight function method combined with a 3-D finite element method. The 3-D finite element method is used to analyze uncracked configuration and provide stress distribution in the region where crack is to occur. Using this stress distribution as input, the 3-D weight function method is used to determine stress intensity factors. Three different loading conditions, i.e. remote tension, remote bending and wedge loading, are considered for a wide range in geometrical parameters. The significance in using 3-D uncracked stress distribution and the difference between single and double corner cracks are studied. Typical crack opening displacements are also provided. Comparisons are made with solutions available in the literature.

Zhao, W.

Calculating 3-D Crack-Opening Behavior Of A Fatigue Crack

Method for calculating through-the-thickness, three-dimensional (3-D) crack-opening behavior in standard fatigue-crack-growth test specimens devised. Involves determination of crack-opening loads from crack-propagation data and finite-element analyses. In method, one considers variation in rate of crack growth along crack front, relationship between rate of crack growth and effective stress intensity factor, and variation in three-dimensional stress-intensity factor along crack front. Developed to provide simpler means for obtaining fatigue-crack-opening behavior loads.

Newman, J. C., Jr.

Stress intensity factors for surface and corner cracks emanating from a wedge-loaded hole

To assist analysis of riveted lap joints, stress intensity factors are determined for surface and corner cracks emanating from a wedge-loaded hole by using a 3-D weight function method in conjunction with a 3-D finite element method. A stress intensity factor equation for surface cracks is also developed to provide a closed-form solution. The equation covers commonly-encountered geometrical ranges and retains high accuracy over the entire range.

Zhao, W.

Analysis of local delaminations caused by angle ply matrix cracks

Two different families of graphite/epoxy laminates with similar layups but different stacking sequences, (0,theta,-theta) sub s and (-theta/theta/0) sub s were analyzed using three-dimensional finite element analysis for theta = 15 and 30 degrees. Delaminations were modeled in the -theta/theta interface, bounded by a matrix crack and the stress free edge. The total strain energy release rate, G, along the delamination front was computed using three different techniques: the virtual crack closure technique (VCCT), the equivalent domain Integral (EDI) technique, and a global energy balance technique. The opening fracture mode component of the strain energy release rate, Gl, along the delamination front was also computed for various delamination lengths using VCCT. The effect of residual thermal and moisture stresses on G was evaluated.

Salpekar, Satish A.

Three-dimensional elastic-plastic finite-element analyses of constraint variations in cracked bodies

Three-dimensional elastic-plastic (small-strain) finite-element analyses were used to study the stresses, deformations, and constraint variations around a straight-through crack in finite-thickness plates for an elastic-perfectly plastic material under monotonic and cyclic loading. Middle-crack tension specimens were analyzed for thicknesses ranging from 1.25 to 20 mm with various crack lengths. Three local constraint parameters, related to the normal, tangential, and hydrostatic stresses, showed similar variations along the crack front for a given thickness and applied stress level. Numerical analyses indicated that cyclic stress history and crack growth reduced the local constraint parameters in the interior of a plate, especially at high applied stress levels. A global constraint factor alpha(sub g) was defined to simulate three-dimensional effects in two-dimensional crack analyses. The global constraint factor was calculated as an average through-the-thickness value over the crack-front plastic region. Values of alpha(sub g) were found to be nearly independent of crack length and were related to the stress-intensity factor for a given thickness.

Newman, J. C., Jr.

An equivalent domain integral method for three-dimensional mixed-mode fracture problems

A general formulation of the equivalent domain integral (EDI) method for mixed mode fracture problems in cracked solids is presented. The method is discussed in the context of a 3-D finite element analysis. The J integral consists of two parts: the volume integral of the crack front potential over a torus enclosing the crack front and the crack surface integral due to the crack front potential plus the crack face loading. In mixed mode crack problems the total J integral is split into J sub I, J sub II, and J sub III representing the severity of the crack front in three modes of deformations. The direct and decomposition methods are used to separate the modes. These two methods were applied to several mixed mode fracture problems, were analyzed, and results were found to agree well with those available in the literature. The method lends itself to be used as a post-processing subroutine in a general purpose finite element program.

Shivakumar, K. N.

Stress concentrations for straight-shank and countersunk holes in plates subjected to tension, bending, and pin loading

A three dimensional stress concentration analysis was conducted on straight shank and countersunk (rivet) holes in a large plate subjected to various loading conditions. Three dimensional finite element analysis were performed with 20 node isoparametric elements. The plate material was assumed to be linear elastic and isotropic, with a Poisson ratio of 0.3. Stress concentration along the bore of the hole were computed for several ratios of hole radius to plate thickness (0.1 to 2.5) and ratios of countersink depth to plate thickness (0.25 to 1). The countersink angles were varied from 80 to 100 degrees in some typical cases, but the angle was held constant at 100 degrees for most cases. For straight shank holes, three types of loading were considered: remote tension, remote bending, and wedge loading in the hole. Results for remote tension and wedge loading were used to estimate stress concentration for simulated rivet in pin loading. For countersunk holes only remote tension and bending were considered. Based on the finite element results, stress concentration equations were developed. Whenever possible, the present results were compared with other numerical solutions and experimental results from the literature.

Shivakumar, K. N.

Parallel computation in a three-dimensional elastic-plastic finite-element analysis

A CRAY parallel processing technique called autotasking was implemented in a three-dimensional elasto-plastic finite-element code. The technique was evaluated on two CRAY supercomputers, a CRAY 2 and a CRAY Y-MP. Autotasking was implemented in all major portions of the code, except the matrix equations solver. Compiler directives alone were not able to properly multitask the code; user-inserted directives were required to achieve better performance. It was noted that the connect time, rather than wall-clock time, was more appropriate to determine speedup in multiuser environments. For a typical example problem, a speedup of 2.1 (1.8 when the solution time was included) was achieved in a dedicated environment and 1.7 (1.6 with solution time) in a multiuser environment on a four-processor CRAY 2 supercomputer. The speedup on a three-processor CRAY Y-MP was about 2.4 (2.0 with solution time) in a multiuser environment.

Shivakumar, K. N.

An inverse method for the calculation of through-thickness fatigue crack closure behavior

An inverse technique was used to calculate through-thickness fatigue crack closure behavior. The through-thickness variation in crack opening stress-intensity factor was calculated by considering the variation in the three-dimensional stress-intensity factor, the variation in crack growth rate along the crack front, and a relationship between the crack growth rate and effective stress-intensity factor range (da/dN-Delta-K(eff)). The three-dimensional stress-intensity factor variation was obtained from an elastic finite element analysis of specific crack-front profiles observed experimentally. The variation in crack growth rate along the crack front was obtained experimentally from comparison of observed crack front changes. The da/dN-Delta-K(eff) relationship was estimated from high stress ratio, constant load amplitude, and fatigue crack growth tests. The through-thickness crack opening stress-intensity factor results agreed with crack opening measurements obtained from fatigue striations, near-tip strain gages, and remote strain and displacement gages.

Dawicke, D. S.

Strain energy release rate distributions for double cantilever beam specimens

A 24-ply composite double cantilever-beam specimen under mode I (opening) loading has been analyzed by a 3D FEM code that calculated along a straight delamination starter for several different specimen materials. An isotropic specimen was found to have a strain-energy release rate distribution which varied along its delamination front due to the boundary-layer effect and another effect associated with the anticlastic curvature of the bent specimen arms. A 0-deg graphite-reinforced epoxy specimen had a nearly-uniform strain-energy release rate distribution which dropped only near the edge, due to the boundary-layer effect, and a +/- 45-deg graphite/epoxy specimen exhibited a pronounced strain-energy release rate variation across the specimen width.

Crews, J. H., Jr.

An equivalent domain integral method for three-dimensional mixed-mode fracture problems

A general formulation of the equivalent domain integral (EDI) method for mixed mode fracture problems in cracked solids is presented. The method is discussed in the context of a 3-D finite element analysis. The J integral consists of two parts: the volume integral of the crack front potential over a torus enclosing the crack front and the crack surface integral due to the crack front potential plus the crack face loading. In mixed mode crack problems the total J integral is split into J sub I, J sub II, and J sub III representing the severity of the crack front in three modes of deformations. The direct and decomposition methods are used to separate the modes. These two methods were applied to several mixed mode fracture problems, were analyzed, and results were found to agree well with those available in the literature. The method lends itself to be used as a post-processing subroutine in a general purpose finite element program.

Shivakumar, K. N.

Stress intensity factors for large aspect ratio surface and corner cracks at a semi-circular notch in a tension specimen

Stress intensity factor solutions for semielliptic surface and quarter-elliptic corner cracks emanating from a semicircular notch in a tension specimen are presented. A three-dimensional finite-element analysis in conjunction with the equivalent domain integral was used to calculate stress intensity factors (SIF). SIF solutions for surface or corner crack (crack length to depth ratio of 2) at a notch are presented for a wide range of crack sizes and notch radii. Results showed that the SIF are larger for larger crack lengths and for larger notch radii. The SIF are nearly constant all along the crack front for deep surface cracks and for all corner cracks analysed.

Shivakumar, K. N.

Parallel computation in a three-dimensional elastic-plastic finite-element analysis

CRAY 2 and CRAY Y-MP evaluation runs are undertaken for a three-dimensional, elastoplastic FEM code implementation of the 'autotasking' parallel-processing technique, which is used in all code components except the matrix-equation solver. For a typical example problem in the case of a four-processor CRAY 2, a speedup factor of 2.1 was found achievable in a dedicated environment, and 1.7 in a multiuser environment; for a three-processor CRAY Y-MP, speedup in a multiuser environment was about 2.4. Only minimal effort was required to implement autotasking in this implementation code.

Shivakumar, K. N.

A three-dimensional equivalent domain integral method for cracked solids

Three-dimensional FEM analysis is presently undertaken via the general formulation of the equivalent-domain integral method for mixed-mode fracture problems in general anisotropic and nonlinear materials. In mixed-mode integrals, the total J-integral is split into three separate modes of fracture; direct and decomposition methods are effective in this separation of the modes, as presently demonstrated by their application to several pure and mixed-mode fracture problems. The results are found to agree well with those reported in the literature, and suggest the usefulness of the method as a FEM-program postprocessing subroutines.

Shivakumar, K. N.

ZIP3D: An elastic and elastic-plastic finite-element analysis program for cracked bodies

ZIP3D is an elastic and an elastic-plastic finite element program to analyze cracks in three dimensional solids. The program may also be used to analyze uncracked bodies or multi-body problems involving contacting surfaces. For crack problems, the program has several unique features including the calculation of mixed-mode strain energy release rates using the three dimensional virtual crack closure technique, the calculation of the J integral using the equivalent domain integral method, the capability to extend the crack front under monotonic or cyclic loading, and the capability to close or open the crack surfaces during cyclic loading. The theories behind the various aspects of the program are explained briefly. Line-by-line data preparation is presented. Input data and results for an elastic analysis of a surface crack in a plate and for an elastic-plastic analysis of a single-edge-crack-tension specimen are also presented.

Shivakumar, K. N.