Engineering topics

Raju, I. S.

Publications and source records attributed to Raju, I. S..

At least 37 records · Page 2

NASA NTRS1999

Modeling the Influence of Stitching on Delamination Growth in Stitched Warp-Knit Composite Lap Joints

The effect of stitches on the failure of a single lap joint configuration was determined in a combined experimental and analytical study. The experimental study was conducted to determine debond growth under static monotonic loading. The stitches were shown to delay the initiation of the debond and provide load transfer beyond the load necessary to completely debond the stitched lap joint. The strain energy release rates at the debond front were calculated using a finite element-based technique. Models of the unstitched configuration showed significant values of modes I and II across the width of the joint and showed that mode III is zero at the centerline but increases near the free edge. Models of the stitched configuration showed that the stitches effectively reduced mode I to zero, but had less of an effect on modes II and III.

Glaessgen, E. H.↗

NASA NTRS1999

Debonding of Stitched Composite Joints: Testing and Analysis

The effect of stitches on the failure of a single lap joint configuration was determined in a combined experimental and analytical study. The experimental study was conducted to determine debond growth under static monotonic loading. The stitches were shown to delay the initiation ofthe debond and provide load transfer beyond the load necessary to completely debond the stitched lap joint. The strain energy release rates at the debond front were calculated using a finite element-based technique. Models of the unstitched configuration showed significant values of modes I and II across the width of the joint and showed that mode III is zero at the centerline but increases near the free edge. Models of the stitched configuration showed that the stitches effectively reduced mode I to zero, but had less of an effect on modes II and III.

Glaessgen, E. H.↗

NASA NTRS1999

Delamination and Stitched Failure in Stitched Composite Joints

The effect of stitches on the failure of a single lap joint configuration was determined in a combined experimental and finite element study. The experimental program was conducted to determine debond growth under static monotonic loading. The stitches were shown to delay the initiation of the debond and provide load transfer beyond the load necessary to completely debond the stitched lap joint. The experimentally determined debond length vs. applied load was used as an input parameter in the finite element analysis of both configurations. The strain energy release rates at the debond from were calculated using plate finite elements. Nonlinear fastener elements were used to model the stitches and multipoint constraints were used to model the contact problem. Models of the unstitched configuration showed significant values of modes I and II across the width of the joint and showed that mode III is zero at the centerline but increases near the free edge. Models of the stitched configuration showed that the stitches were effective in reducing mode I to zero, but had less of an effect on modes II and III.

Glaessgen, E. H.↗

NASA NTRS1998

Fracture Mechanics Analysis of Stitched Stiffener-Skin Debonding

An analysis based on plate finite elements and the virtual crack closure technique has been implemented to study the effect of stitching on mode I and mode II strain energy release rates for debond configurations. The stitches were modeled as discrete nonlinear fastener elements with a compliance determined by experiment. The axial and shear behavior of the stitches was considered, however, the two compliances and failure loads were assumed to be independent. Both a double cantilever beam (mode I) and a mixed mode skin-stiffener debond configuration were studied. In the double cantilever beam configurations, G(sub I) began to decrease once the debond had grown beyond the first row of stitches and was reduced to zero for long debonds. In the mixed-mode skin-stiffener configurations, G(sub I) showed a similar behavior as in the double cantilever beam configurations, however, G(sub u), continued to increase with increasing debond length.

Glaessgen, E. H.↗

NASA NTRS1998

Effect of Shear Deformation and Continuity on Delamination Modelling with Plate Elements

The effects of several critical assumptions and parameters on the computation of strain energy release rates for delamination and debond configurations modeled with plate elements have been quantified. The method of calculation is based on the virtual crack closure technique (VCCT), and models that model the upper and lower surface of the delamination or debond with two-dimensional (2D) plate elements rather than three-dimensional (3D) solid elements. The major advantages of the plate element modeling technique are a smaller model size and simpler geometric modeling. Specific issues that are discussed include: constraint of translational degrees of freedom, rotational degrees of freedom or both in the neighborhood of the crack tip; element order and assumed shear deformation; and continuity of material properties and section stiffness in the vicinity of the debond front, Where appropriate, the plate element analyses are compared with corresponding two-dimensional plane strain analyses.

Glaessgen, E. H.↗

NASA NTRS1998

Debond Analyses for Stitched Composite Structures

The effect of stitching on mode I and mode II strain energy release rates for debond configurations is studied using an analysis based on plate finite elements and the virtual crack closure technique. The stitches were modeled as discrete nonlinear fastener elements with a compliance determined by experiment. The axial and shear behavior of the stitches was considered with both the compliances and failure loads assumed to be independent. The mode I strain energy release rate, G(sub I), was shown to decrease once the debond had grown beyond the first row of stitches and was reduced to zero for long debonds, however, the mode II strain energy release rate, G(sub II), continued to be of significant magnitude over the range of debond lengths considered.

Glaessgen, E. H.↗

NASA NTRS1997

A Review of the NASA Textile Composites Research

During the past 15 years NASA has taken the lead role in exploiting the benefits of textile reinforced composite materials for application to aircraft structures. The NASA Advanced Composites Technology (ACT) program was started in 1989 to develop composite primary structures for commercial transport airplanes with costs that are competitive with metal structures. As part of this program, several contractors investigated the cost, weight, and performance attributes of textile reinforced composites. Textile composites made using resin transfer molding type processes were evaluated for numerous applications. Methods were also developed to predict resin infiltration and flow in textile preforms and to predict and measure mechanical properties of the textile composites. This paper describes the salient results of that program.

Poe, C. C., Jr.↗

NASA NTRS1993

surf3d: A 3-D finite-element program for the analysis of surface and corner cracks in solids subjected to mode-1 loadings

A computer program, surf3d, that uses the 3D finite-element method to calculate the stress-intensity factors for surface, corner, and embedded cracks in finite-thickness plates with and without circular holes, was developed. The cracks are assumed to be either elliptic or part eliptic in shape. The computer program uses eight-noded hexahedral elements to model the solid. The program uses a skyline storage and solver. The stress-intensity factors are evaluated using the force method, the crack-opening displacement method, and the 3-D virtual crack closure methods. In the manual the input to and the output of the surf3d program are described. This manual also demonstrates the use of the program and describes the calculation of the stress-intensity factors. Several examples with sample data files are included with the manual. To facilitate modeling of the user's crack configuration and loading, a companion program (a preprocessor program) that generates the data for the surf3d called gensurf was also developed. The gensurf program is a three dimensional mesh generator program that requires minimal input and that builds a complete data file for surf3d. The program surf3d is operational on Unix machines such as CRAY Y-MP, CRAY-2, and Convex C-220.

Raju, I. S.↗

NASA NTRS1993

Analysis and sizing of Mars aerobrake structure

A cone-sphere aeroshell structure for aerobraking into Martian atmosphere is studied. Using this structural configuration, a space frame load-bearing structure is proposed. To generate this structure efficiently and to perform a variety of studies of several configurations, a mesh generator that utilizes only a few configurational parameters is developed. A finite element analysis program that analyzes space frame structures was developed. A sizing algorithm that arrives at a minimum mass configuration was developed and integrated into the finite element analysis program. A typical 135-ft-diam aerobrake configuration was analyzed and sized. The minimum mass obtained in this study using high modulus graphite/epoxy composite material members is compared with the masses obtained from two other aerobrake structures using lightweight erectable tetrahedral truss and part-spherical truss configurations. Excellent agreement for the minimum mass was obtained with the three different aerobrake structures. Also, the minimum mass using the present structure was obtained when the supports were not at the base but at about 75 percent of the base diameter.

Raju, I. S.↗

NASA NTRS1993

A displacement gradient BEM for accurate stress computation near boundaries in 2-D anisotropic problems

A displacement gradient method for 2D anisotropic elasticity problems is presented, which effectively minimizes the boundary layer effect through a two-step procedure. First, the boundary integral equations are solved for the unknown boundary displacements and tractions. Second, a direct integral equation for displacement gradients is developed in terms of boundary tractions. Three methods based on different evaluation procedures and locations for determining the displacement gradients are proposed. In the first method the displacement gradients are averaged at nodes common to adjacent elements. The second method stores the gradients element-wise. In the third method, the gradients are evaluated at the nodes of discontinuous elements. The three methods are applied to near-isotropic plates with circular and elliptic cutouts. It is concluded that all three methods can yield accurate stress distributions.

Sistla, R.↗

NASA NTRS1993

Coupling finite and boundary element methods for 2-D elasticity problems

A finite element-boundary element (FE-BE) coupling method for two-dimensional elasticity problems is developed based on a weighted residual variational method in which a portion of the domain of interest is modeled by FEs and the remainder of the region by BEs. The performance of the FE-BE coupling method is demonstrated via applications to a simple 'patch test' problem and three-crack problems. The method passed the patch tests for various modeling configurations and yielded accurate strain energy release rates for the crack problems studied.

Krishnamurthy, T.↗

NASA NTRS1993

Computation of strain energy release rates for skin-stiffener debonds modeled with plate elements

An efficient method for predicting the strength of debonded composite skin-stiffener configurations is presented. This method, which is based on fracture mechanics, models the skin and the stiffener with two-dimensional (2D) plate elements instead of three-dimensional (3D) solid elements. The skin and stiffener flange nodes are tied together by two modeling techniques. In one technique, the corresponding flange and skin nodes are required to have identical translational and rotational degrees-of-freedom. In the other technique, the corresponding flange and skin nodes are only required to have identical translational degrees-of-freedom. Strain energy release rate formulas are proposed for both modeling techniques. These formulas are used for skin-stiffener debond cases with and without cylindrical bending deformations. The cylindrical bending results are compared with plane-strain finite element results. Excellent agreement between the two sets of results is obtained when the second technique is used. Thus, from these limited studies, a preferable modeling technique for skin-stiffener debond analysis using plate elements is established.

Wang, J. T.↗

NASA NTRS1993

Fracture mechanics analyses for skin-stiffener debonding

The debond configurations presently subjected to 3D FEM fracture mechanics analyses are respectively of the flange-skin strip and skin-stiffener configuration type. Two methods employing the virtual crack closure technique were used to evaluate the strain energy release rate, or 'G-value' distributions across the debond front. Both methods yielded nearly identical G-value distributions for the debond configurations studied; they were compared with plane strain and shell analyses results from the literature for the flange skin strip configuration, and found to be in good agreement. Mode II is dominant for the skin-stiffener debond configuration.

Raju, I. S.↗

NASA NTRS1992

A boundary element alternating method for two-dimensional mixed-mode fracture problems

A boundary element alternating method, denoted herein as BEAM, is presented for two dimensional fracture problems. This is an iterative method which alternates between two solutions. An analytical solution for arbitrary polynomial normal and tangential pressure distributions applied to the crack faces of an embedded crack in an infinite plate is used as the fundamental solution in the alternating method. A boundary element method for an uncracked finite plate is the second solution. For problems of edge cracks a technique of utilizing finite elements with BEAM is presented to overcome the inherent singularity in boundary element stress calculation near the boundaries. Several computational aspects that make the algorithm efficient are presented. Finally, the BEAM is applied to a variety of two dimensional crack problems with different configurations and loadings to assess the validity of the method. The method gives accurate stress intensity factors with minimal computing effort.

Raju, I. S.↗

NASA NTRS1992

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.↗

NASA NTRS1992

GENSURF: A mesh generator for 3D finite element analysis of surface and corner cracks in finite thickness plates subjected to mode-1 loadings

A computer program that generates three-dimensional (3D) finite element models for cracked 3D solids was written. This computer program, gensurf, uses minimal input data to generate 3D finite element models for isotropic solids with elliptic or part-elliptic cracks. These models can be used with a 3D finite element program called surf3d. This report documents this mesh generator. In this manual the capabilities, limitations, and organization of gensurf are described. The procedures used to develop 3D finite element models and the input for and the output of gensurf are explained. Several examples are included to illustrate the use of this program. Several input data files are included with this manual so that the users can edit these files to conform to their crack configuration and use them with gensurf.

Raju, I. S.↗

NASA NTRS1992

Coupling finite and boundary element methods for two-dimensional potential problems

A finite element-boundary element (FE-BE) coupling method based on a weighted residual variational method is presented for potential problems, governed by either the Laplace or the Poisson equations. In this method, a portion of the domain of interest is modeled by finite elements (FE) and the remainder of the region by boundary elements (BE). Because the BE fundamental solutions are valid for infinite domains, a procedure that limits the effects of the BE fundamental solution to a small region adjacent to the FE region, called the transition region (TR), is developed. This procedure involves a judicious choice of functions called the transition (T) functions that have unit values on the BE-TR interface and zero values on the FE-TR interface. The present FE-BE coupling algorithm is shown to be independent of the extent of the transition region and the choice of the transition functions. Therefore, transition regions that extend to only one layer of elements between FE and BE regions and the use of simple linear transition functions work well.

Krishnamurthy, T.↗

NASA NTRS1992

An independent refinement and integration procedure in multiregion finite element analysis

An independent refinement and integration procedure is developed to couple together independently modeled (global and local) regions in a single analysis. The models can have different levels of refinement and along the interface between them the finite element nodes need not coincide with one another. In the local model all the nodes except the nodes at the interface are statically condensed and the reduced stiffness matrix is obtained. For this static condensation a modified frontal solution technique is employed. A spline interpolation function that satisfies the linear isotropic plate bending differential equation is used to relate the local model interface nodal displacements to the global model interface displacements. The proposed independent refinement and integration procedure is evaluated by applying it to two- and three-dimensional cases involving inplane and out-of-plane deformation. The procedure yielded very accurate results for all the examples studied.

Krishnamurthy, T.↗