Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shell buckling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

536 records · Page 30

Development of Curved-Plate Elements for the Exact Buckling Analysis of Composite Plate Assemblies Including Transverse Shear Effects

The analytical formulation of curved-plate non-linear equilibrium equations that include transverse-shear-deformation effects is presented. A unified set of non-linear strains that contains terms from both physical and tensorial strain measures is used. Using several simplifying assumptions, linearized, stability equations are derived that describe the response of the plate just after bifurcation buckling occurs. These equations are then modified to allow the plate reference surface to be located a distance z(c), from the centroid surface which is convenient for modeling stiffened-plate assemblies. The implementation of the new theory into the VICONOPT buckling and vibration analysis and optimum design program code is described. Either classical plate theory (CPT) or first-order shear-deformation plate theory (SDPT) may be selected in VICONOPT. Comparisons of numerical results for several example problems with different loading states are made. Results from the new curved-plate analysis compare well with closed-form solution results and with results from known example problems in the literature. Finally, a design-optimization study of two different cylindrical shells subject to uniform axial compression is presented.

McGowan, David M.↗

Finite Elements Analysis of a Composite Semi-Span Test Article With and Without Discrete Damage

AS&M Inc. performed finite element analysis, with and without discrete damage, of a composite semi-span test article that represents the Boeing 220-passenger transport aircraft composite semi-span test article. A NASTRAN bulk data file and drawings of the test mount fixtures and semi-span components were utilized to generate the baseline finite element model. In this model, the stringer blades are represented by shell elements, and the stringer flanges are combined with the skin. Numerous modeling modifications and discrete source damage scenarios were applied to the test article model throughout the course of the study. This report details the analysis method and results obtained from the composite semi-span study. Analyses were carried out for three load cases: Braked Roll, LOG Down-Bending and 2.5G Up-Bending. These analyses included linear and nonlinear static response, as well as linear and nonlinear buckling response. Results are presented in the form of stress and strain plots. factors of safety for failed elements, buckling loads and modes, deflection prediction tables and plots, and strainage prediction tables and plots. The collected results are presented within this report for comparison to test results.

Lovejoy, Andrew E.↗

Finite element structural redesign by large admissible perturbations

In structural redesign, two structural states are involved; the baseline (known) State S1 with unacceptable performance, and the objective (unknown) State S2 with given performance specifications. The difference between the two states in performance and design variables may be as high as 100 percent or more depending on the scale of the structure. A Perturbation Approach to Redesign (PAR) is presented to relate any two structural states S1 and S2 that are modeled by the same finite element model and represented by different values of the design variables. General perturbation equations are derived expressing implicitly the natural frequencies, dynamic modes, static deflections, static stresses, Euler buckling loads, and buckling modes of the objective S2 in terms of its performance specifications, and S1 data and Finite Element Analysis (FEA) results. Large Admissible Perturbation (LEAP) algorithms are implemented in code RESTRUCT to define the objective S2 incrementally without trial and error by postprocessing FEA results of S1 with no additional FEAs. Systematic numerical applications in redesign of a 10 element 48 degree of freedom (dof) beam, a 104 element 192 dof offshore tower, a 64 element 216 dof plate, and a 144 element 896 dof cylindrical shell show the accuracy, efficiency, and potential of PAR to find an objective state that may differ 100 percent from the baseline design.

Bernitsas, Michael M.↗

A Stability Analysis of Cylindrical Panels Using a Finite Element Formulation

A cylindrical finite element suitable for the linear stability analysis of cylindrical shells is developed. Energy principles and variational methods lead to a problem formulation which lends itself to physical interpretations of the governing matrices of the finite element. By properly grouping the terms which result from taking the second variation of the potential energy of the element, it is possible to identify three distinct types of matrices. The first matrix is the conventional stiffness matrix; the second is an initial stress stiffness matrix; and the third is an initial displacement stiffness matrix. With the assumption of linearity, the buckling problem is stated in terms of the classical linear real eigenvalue equation. This problem formulation was programmed on the CDC 6600 series computer. The computer program is used to analyze the buckling of a variety of structures. Columns, arches, flat plates and curved panels with and without cutouts are considered.

Snyder, R. E.↗

Fracture Analysis of the FAA/NASA Wide Stiffened Panels

This paper presents the fracture analyses conducted on the FAA/NASA stiffened and unstiffened panels using the STAGS (STructural Analysis of General Shells) code with the critical crack-tip-opening angle (CTOA) fracture criterion. The STAGS code with the "plane-strain" core option was used in all analyses. Previous analyses of wide, flat panels have shown that the high-constraint conditions around a crack front, like plane strain, has to be modeled in order for the critical CTOA fracture criterion to predict wide panel failures from small laboratory tests. In the present study, the critical CTOA value was determined from a wide (unstiffened) panel with anti-buckling guides. The plane-strain core size was estimated from previous fracture analyses and was equal to about the sheet thickness. Rivet flexibility and stiffener failure was based on methods and criteria, like that currently used in industry. STAGS and the CTOA criterion were used to predict load-against-crack extension for the wide panels with a single crack and multiple-site damage cracking at many adjacent rivet holes. Analyses were able to predict stable crack growth and residual strength with a few percent (5%) of stiffened panel tests results but over predicted the buckling failure load on a unstiffened panel with a single crack by 10%.

Seshadri, B. R.↗

Fracture Analysis of the FAA/NASA Wide Stiffened Panels

This paper presents the fracture analyses conducted on the FAA/NASA stiffened and unstiffened panels using the STAGS (STructural Analysis of General Shells) code with the critical crack-tip-opening angle (CTOA) fracture criterion. The STAGS code with the "plane-strain" core option was used in all analyses. Previous analyses of wide, flat panels have shown that the high-constraint conditions around a crack front, like plane strain, has to be modeled in order for the critical CTOA fracture criterion to predict wide panel failures from small laboratory tests. In the present study, the critical CTOA value was determined from a wide (unstiffened) panel with anti-buckling guides. The plane-strain core size was estimated from previous fracture analyses and was equal to about the sheet thickness. Rivet flexibility and stiffener failure was based on methods and criteria, like that currently used in industry. STAGS and the CTOA criterion were used to predict load-against-crack extension for the wide panels with a single crack and multiple-site damage cracking at many adjacent rivet holes. Analyses were able to predict stable crack growth and residual strength within a few percent (5%) of stiffened panel tests results but over predicted the buckling failure load on an unstiffened panel with a single crack by 10%.

Seshadri, B. R.↗

Derivatives of buckling loads and vibration frequencies with respect to stiffness and initial strain parameters

A uniform variational approach to sensitivity analysis of vibration frequencies and bifurcation loads of nonlinear structures is developed. Two methods of calculating the sensitivities of bifurcation buckling loads and vibration frequencies of nonlinear structures, with respect to stiffness and initial strain parameters, are presented. A direct method requires calculation of derivatives of the prebuckling state with respect to these parameters. An adjoint method bypasses the need for these derivatives by using instead the strain field associated with the second-order postbuckling state. An operator notation is used and the derivation is based on the principle of virtual work. The derivative computations are easily implemented in structural analysis programs. This is demonstrated by examples using a general purpose, finite element program and a shell-of-revolution program.

Haftka, Raphael T.↗

Analysis of Composite Panel-Stiffener Debonding Using a Shell/3D Modeling Technique

Interlaminar fracture mechanics has proven useful for characterizing the onset of delaminations in composites and has been used successfully primarily to investigate onset in fracture toughness specimens and laboratory size coupon type specimens. Future acceptance of the methodology by industry and certification authorities, however, requires the successful demonstration of the methodology on the structural level. For this purpose, a panel was selected that is reinforced with stiffeners. Shear loading causes the panel to buckle, and the resulting out-of-plane deformations initiate skin/stiffener separation at the location of an embedded defect. A small section of the stiffener foot, web and noodle as well as the panel skin in the vicinity of the delamination front were modeled with a local 3D solid model. Across the width of the stiffener foot, the mixedmode strain energy release rates were calculated using the virtual crack closure technique. A failure index was calculated by correlating the results with a mixed-mode failure criterion of the graphite/epoxy material. Computed failure indices were compared to corresponding results where the entire web was modeled with shell elements and only a small section of the stiffener foot and panel were modeled locally with solid elements. Including the stiffener web in the local 3D solid model increased the computed failure index. Further including the noodle and transition radius in the local 3D solid model changed the local distribution across the width. The magnitude of the failure index decreased with increasing transition radius and noodle area. For the transition radii modeled, the material properties used for the noodle area had a negligible effect on the results. The results of this study are intended to be used as a guide for conducting finite element and fracture mechanics analyses of delamination and debonding in complex structures such as integrally stiffened panels.

Krueger, Ronald↗

Residual Strength Analyses of Riveted Lap-Splice Joints

The objective of this paper was to analyze the crack-linkup behavior in riveted-stiffened lap-splice joint panels with small multiple-site damage (MSD) cracks at several adjacent rivet holes. Analyses are based on the STAGS (STructural Analysis of General Shells) code with the critical crack-tip-opening angle (CTOA) fracture criterion. To account for high constraint around a crack front, the "plane strain core" option in STAGS was used. The importance of modeling rivet flexibility with fastener elements that accurately model load transfer across the joint is discussed. Fastener holes are not modeled but rivet connectivity is accounted for by attaching rivets to the sheet on one side of the cracks that simulated both the rivet diameter and MSD cracks. Residual strength analyses made on 2024-T3 alloy (1.6-mm thick) riveted-lap-splice joints with a lead crack and various size MSD cracks were compared with test data from Boeing Airplane Company. Analyses were conducted for both restrained and unrestrained buckling conditions. Comparison of results from these analyses and results from lap-splice-joint test panels, which were partially restrained against buckling indicate that the test results were bounded by the failure loads predicted by the analyses with restrained and unrestrained conditions.

Seshadri, B. R.↗

Residual Strength Analyses of Monolithic Structures

Finite-element fracture simulation methodology predicts the residual strength of damaged aircraft structures. The methodology uses the critical crack-tip-opening-angle (CTOA) fracture criterion to characterize the fracture behavior of the material. The CTOA fracture criterion assumes that stable crack growth occurs when the crack-tip angle reaches a constant critical value. The use of the CTOA criterion requires an elastic- plastic, finite-element analysis. The critical CTOA value is determined by simulating fracture behavior in laboratory specimens, such as a compact specimen, to obtain the angle that best fits the observed test behavior. The critical CTOA value appears to be independent of loading, crack length, and in-plane dimensions. However, it is a function of material thickness and local crack-front constraint. Modeling the local constraint requires either a three-dimensional analysis or a two-dimensional analysis with an approximation to account for the constraint effects. In recent times as the aircraft industry is leaning towards monolithic structures with the intention of reducing part count and manufacturing cost, there has been a consistent effort at NASA Langley to extend critical CTOA based numerical methodology in the analysis of integrally-stiffened panels.In this regard, a series of fracture tests were conducted on both flat and curved aluminum alloy integrally-stiffened panels. These flat panels were subjected to uniaxial tension and during the test, applied load-crack extension, out-of-plane displacements and local deformations around the crack tip region were measured. Compact and middle-crack tension specimens were tested to determine the critical angle (wc) using three-dimensional code (ZIP3D) and the plane-strain core height (hJ using two-dimensional code (STAGS). These values were then used in the STAGS analysis to predict the fracture behavior of the integrally-stiffened panels. The analyses modeled stable tearing, buckling, and crack branching at the integral stiffener using different values of critical CTOA for different material thicknesses and orientation. Comparisons were made between measured and predicted load-crack extension, out-of-plane displacements and local deformations around the crack tip region. Simultaneously, three-dimensional capabilities to model crack branching and to monitor stable crack growth of multiple cracks in a large thick integrally-stiffened flat panels were implemented in three-dimensional finite element code (ZIP3D) and tested by analyzing the integrally-stiffened panels tested at Alcoa. The residual strength of the panels predicted from STAGS and ZP3D code compared very well with experimental data. In recent times, STAGS software has been updated with new features and now one can have combinations of solid and shell elements in the residual strength analysis of integrally-stiffened panels.

Forth, Scott↗

SCBUCKLE user's manual: Buckling analysis program for simple supported and clamped panels

The program SCBUCKLE calculates the buckling loads and mode shapes of cylindrically curved, rectangular panels. The panel is assumed to have no imperfections. SCBUCKLE is capable of analyzing specially orthotropic symmetric panels (i.e., A(sub 16) = A(sub 26) = 0.0, D(sub 16) = D(sub 26) = 0.0, B(sub ij) = 0.0). The analysis includes first-order transverse shear theory and is capable of modeling sandwich panels. The analysis supports two types of boundary conditions: either simply supported or clamped on all four edges. The panel can be subjected to linearly varying normal loads N(sub x) and N(sub y) in addition to a constant shear load N(sub xy). The applied loads can be divided into two parts: a preload component; and a variable (eigenvalue-dependent) component. The analysis is based on the modified Donnell's equations for shallow shells. The governing equations are solved by Galerkin's method.

Cruz, Juan R.↗

Development of Curved-Plate Elements for the Exact Buckling Analysis of Composite Plate Assemblies Including Transverse-Shear Effects

The analytical formulation of curved-plate non-linear equilibrium equations including transverse-shear-deformation effects is presented. A unified set of non-linear strains that contains terms from both physical and tensorial strain measures is used. Linearized, perturbed equilibrium equations (stability equations) that describe the response of the plate just after buckling occurs are derived. These equations are then modified to allow the plate reference surface to be located a distance z(sub c) from the centroidal surface. The implementation of the new theory into the VICONOPT exact buckling and vibration analysis and optimum design computer program is described. The terms of the plate stiffness matrix using both classical plate theory (CPT) and first-order shear-deformation plate theory (SDPT) are presented. The effects of in-plane transverse and in-plane shear loads are included in the in-plane stability equations. Numerical results for several example problems with different loading states are presented. Comparisons of analyses using both physical and tensorial strain measures as well as CPT and SDPT are made. The computational effort required by the new analysis is compared to that of the analysis currently in the VICONOPT program. The effects of including terms related to in-plane transverse and in-plane shear loadings in the in-plane stability equations are also examined. Finally, results of a design-optimization study of two different cylindrical shells subject to uniform axial compression are presented.

McGowan, David M.↗

Development of Curved-Plate Elements for the Exact Buckling Analysis of Composite Plate Assemblies Including Transverse-Shear Effects

The analytical formulation of curved-plate non-linear equilibrium equations including transverse-shear-deformation effects is presented. The formulation uses the principle of virtual work. A unified set of non-linear strains that contains terms from both physical and tensorial strain measures is used. Linearized, perturbed equilibrium equations (stability equations) that describe the response of the plate just after buckling occurs are then derived after the application of several simplifying assumptions. These equations are then modified to allow the reference surface of the plate to be located at a distance z(sub c) from the centroidal surface. The implementation of the new theory into the VICONOPT exact buckling and vibration analysis and optimum design computer program is described as well. The terms of the plate stiffness matrix using both Classical Plate Theory (CPT) and first-order Shear-Deformation Plate Theory (SDPT) are presented. The necessary steps to include the effects of in-plane transverse and in-plane shear loads in the in-plane stability equations are also outlined. Numerical results are presented using the newly implemented capability. Comparisons of results for several example problems with different loading states are made. Comparisons of analyses using both physical and tensorial strain measures as well as CPT and SDPF are also made. Results comparing the computational effort required by the new analysis to that of the analysis currently in the VICONOPT program are presented. The effects of including terms related to in-plane transverse and in-plane shear loadings in the in-plane stability equations are also examined. Finally, results of a design-optimization study of two different cylindrical shells subject to uniform axial compression are presented.

McGowan, David Michael↗

DESIGN AND TEST EXPERIENCES WITH INSTABILITY OF MAJOR AIRFRAME COMPONENTS

Two test incidents involving instability of large scale commercially built structures are described. Two classes of structure are discussed; the first, a fuselage with skin designed to buckle at low stress, and the second, a wing whose surface remains unbuckled to failure. The structure in the region of failure is defined and the failures described and illustrated. Insofar as possible, the stresses in the critical area at the time of failure are reported and compared to strength determined by analysis. Both fuselage and wing surfaces were observed to fail in the mode of a medium range column when adequate support was provided by ribs and frames. Initial failure in the wing example was premature due to a design deficiency in rib strength. A clear illustration is given of the effect of rib stiffness on wing surface stability. Adding stiffness to wing ribs increases the limit of surface stability to the theoretical flat panel value.

SHELL STABILITY↗