* developed by Dr. Charles C. Rankin and others. STAGS continued to be developed by Dr. Rankin until his death in August, 2012. Many of the unique and sophisticated capabilities of STAGS described below have been or are now being added to the capabilities of the widely used commercial software, ABAQUS.
The following text is included in the paper, “Optimization of an axially compressed ring and stringer stiffened cylindrical shell with a general buckling modal imperfection”, AIAA 48th Structures, Structural Dynamics, and Materials Conference, Paper no. AIAA20072216, 2007. Tables and figures referred to here are in that paper.
STAGS (STructural Analysis of General Shells [17]) is a finite element code for generalpurpose nonlinear analysis of stiffened shell structures of arbitrary shape and complexity. Its capabilities include stress, stability, vibration, and transient analyses with both material and geometric nonlinearities permitted in all analysis types. STAGS includes enhancements, such as a higher order thick shell element, more advanced nonlinear solution strategies, and more comprehensive postprocessing features such as a link with STAPL, a postprocessor used to generate many of the figures in this paper: figures that display the STAGS model, such as Figs. 1ac and 2, for example. (Figure 2 from the paper cited above is shown below.)
Research and development of STAGS by Rankin, Brogan, Almroth, Stanley, Cabiness, Stehlin and others, formerly of the Computational Mechanics Department of the Lockheed Martin Advanced Technology Center, has been under continuous sponsorship from U.S. government agencies for the past 40 years. During this time particular emphasis has been placed on improvement of the capability to solve difficult nonlinear problems such as the prediction of the behavior of axially compressed stiffened panels loaded far into their locally postbuckled states. STAGS has been extensively used worldwide for the evaluation of stiffened panels and shells loaded well into their locally postbuckled states. See [7], for example.
A large rotation algorithm that is independent of the finite element library has been incorporated into STAGS [4]. With this algorithm there is no artificial stiffening due to large rotations. The finite elements in the STAGS library do not store energy under arbitrary rigidbody motion, and the first and second variations of the strain energy are consistent. These properties lead to quadratic convergence during Newton iterations.
Solution control in nonlinear problems includes specification of load levels or use of the advanced RiksCrisfield path parameter [7] that enables traversal of limit points into the postbuckling regime. Two load systems with different histories (Load Sets A and B) can be defined and controlled separately during the solution process. Flexible restart procedures permit switching from one strategy to another during an analysis, including shifts from bifurcation buckling to nonlinear collapse analyses and back and shifts from static to transient and transient to static analyses with modified boundary conditions and loading. STAGS provides solutions to the generalized eigenvalue problem for buckling and vibration from a linear (Fig. 24) or nonlinear (Figs. 26, 27) stress state.
Quadric surfaces can be modeled with minimal user input as individual substructures called "shell units" in which the analytic geometry is represented exactly. "Shell units" can be connected along edges or internal grid lines with partial or complete compatibility. In this way complex structures can be assembled from relatively simple units. Alternatively, a structure of arbitrary shape can be modeled with use of an "element unit".
Geometric imperfections can be generated automatically in a variety of ways, thereby permitting imperfectionsensitivity studies to be performed. For example, imperfections can be generated by superposition of several buckling modes determined from previous linear and nonlinear STAGS analyses of a given case. (See Parts 47 of Table 9 and Figs. 24, 26, and 27, for example).
A variety of material models is available, including both plasticity and creep. STAGS handles isotropic and anisotropic materials, including composites consisting of up to 60 layers of arbitrary orientation. Four plasticity models are available, including isotropic strain hardening, the White Besseling (mechanical sublayer model), kinematic strain hardening, and deformation theory.
Two independent load sets, each composed from simple parts that may be specified with minimal input, define a spatial variation of loading. Any number of point loads, prescribed displacements, line loads, surface tractions, thermal loads, and "live" pressure (hydrostatic pressure which remains normal to the shell surface throughout large deformations) can be combined to make a load set. For transient analysis the user may select from a menu of loading histories, or a general temporal variation may be specified in a userwritten subroutine.
Boundary conditions (B.C.) may be imposed either by reference to certain standard conditions or by the use of single and multipoint constraints. Simple support, symmetry, antisymmetry, clamped, or userspecified B.C. can be defined on a "shell unit" edge. Singlepoint constraints that allow individual freedoms to be free, fixed, or a prescribed nonzero value may be applied to grid lines and surfaces in "shell units" or "element units". A useful feature for buckling analysis allows these constraints to differ for the prebuckling stress and eigenvalue analyses. Langrangian constraint equations containing up to 100 terms may be defined to impose multipoint constraints.
STAGS has a variety of finite elements suitable for the analysis of stiffened plates and shells. Simple four node quadrilateral plate elements with a cubic lateral displacement field (called "410" and "411" elements) are effective and efficient for the prediction of postbuckling thin shell response. A linear (410) or quadratic (411) membrane interpolation can be selected. For thicker shells in which transverse shear deformation is important (and for the thinshell cases described in this paper), STAGS provides the Assumed Natural Strain (ANS) ninenode element (called "480" element). A twonode beam element compatible with the four node quadrilateral plate element is provided to simulate stiffeners and beam assemblies. Other finite elements included in STAGS are described in the STAGS literature [17].
Graphic from blah blah, more blah blah.
Abstracts and Documentation

A more detailed abstract on STAGS.

STAGS (STructural Analysis of General Shells is a finite element code for generalpurpose nonlinear analysis of stiffened shell structures of arbitrary shape and complexity. Its capabilities include stress, stability, vibration, and transient analyses with both material and geometric nonlinearities permitted in all analysis types. STAGS includes enhancements, such as a higher order thick shell element, more advanced nonlinear solution strategies, and more comprehensive postprocessing features such as a link with STAPL, a postprocessor used to generate figures that display undeformed and deformed STAGS models. Research and development of STAGS by Rankin, Brogan, Almroth, Stanley, Cabiness, Stehlin and others, formerly of the Computational Mechanics Department of the Lockheed Martin Advanced Technology Center, has been under continuous sponsorship from U.S. government agencies for the past 40 years. During this time particular emphasis has been placed on improvement of the capability to solve difficult nonlinear problems such as the prediction of the behavior of axially compressed stiffened panels loaded far into their locally postbuckled states.

STAGS is fundamentally a finite element code for generalpurpose analysis of shell structures of arbitrary shape and complexity, with additional capabilities for analysis of solids and other types of structural configurations. Shells to be analyzed by STAGS may be thin or thick, unstiffened or stiffened, with stiffeners modeled either as beams or as shells. The availability of numerous wallfabrications and stiffenercrosssection options combined with a variety of material models permits tremendous flexibility in modeling a wide spectrum of construction types. Analysis capabilities include stress, stability, vibration, and transient analyses, with both material and geometric nonlinearities permissible. STAGS has a long and distinguished history starting in the late 1960’s in the Lockheed Missiles and Space Company’s Research and Development Division (R&DD) as a companysponsored research code to study the stability of shells and continuing to the present at LockheedMartin’s Advanced Technology Center. The primary goal of the STAGS group at Lockheed and at LockheedMartin has always been to achieve the best shell stability code available anywhere. Research with and development of STAGS has been under continuous sponsorship from U.S. government agencies and Lockheed’s Independent Research program from the beginning.
References
 The STAGS computer code, NASA CR2950, NASA Langley Research Center, Hampton, VA., 1978. ,

Enhancements to the STAGS computer code, NASA Contractor Report 4000, Contract NAS116723, November 1986.
,The research conducted during the last three years has greatly enhanced the power of the STAGS family of programs. Members of this family include STAGSCl and RRSYS. As a result of improvements implemented during these last three years, it is now possible to address the full collapse of a structural system, up to and beyond critical points where its resistance to the applied loads vanishes or suddenly changes. This also includes the important class of problems where a multiplicity of solutions exists at a given point (bifurcation), and where until now no solution could be obtained along any alternate (secondary) load path with any standard production finiteelement code. Formerly, when rotations for any part of a collapsing structure exceeded the small rotation regime (less than 10 degrees), finite elements in that region of the structure displayed unpredictable behavior that arose out of a violation of the basic assumptions concerning deformations within the individual elements. All elements in STAGS suffered to some degree, with the more economical elements (such as the SH410 shell element) locking almost completely as rotations grew. Now, the latest version of STAGSC1 includes an elementindependent corotational formulation for all static nonlinear collapse analyses that virtually eliminates all problems related to largerotations for small strain problems. Some of the remaining improvements to STAGSC1 include a postprocessor that generates secondary solution data (strains, stresses, and resultants) from previouslysaved primary data (displacements, velocities, eigenvectors, plastic strains). This postprocessor eliminates the need for saving bulky secondary data that may or may not later be needed for printout or display (e.g. STAPL). A new data translator formats primary solution data created on one type of machine (such as a CDC or CRAY) for simple transfer to another machine (like a VAX/VMS) for local, interactive, and inexpensive postprocessing or solution restart.

The Numerical Simulation of the collapse process of axially compressed cylindrical shells with measured imperfections, Report LR 705, Technical University of Delft, March 1993
,Thin cylindrical shells in compression buckle explosively as is known from numerous experiments. What happens visually during the buckling phase has been revealed only in some isolated instances where high speed video recordings were made of these phenomena. But as far as it is known the actual buckling process of a cylindrical shell has never been numerically simulated. We believe that the most likely reason for this is that the whole arsenal of computational tools that are necessary for this task were not yet available at the time such calculations were contemplated. But now the situation has changed. The potential of present day finite element software and the power of computer hardware have increased so much that the restrictions mentioned above no longer exist. To illustrate this point we present here the numerical computation of the collapse process of a cylindrical shell under prescribed edge displacements. The shell that we chose for our computational experiment was manufactured, measured and tested some twenty years ago by Singer, Arbocz and Babcock under carefully controlled conditions in a laboratory environment.

On the solution of mode jumping phenomena in thin walled shell structures, Computer Methods in Applied Mechanics and Engineering, First ASCE/ASM/SES Mechanics Conference, Charlottesville, VA, 1996: Vol. 136, pp. 5992
,This paper is an investigation into the merits of a hybrid procedure for the numerical simulation of transient buckling problems. The procedure consists of the combination of a classical pathfollowing method with a transient integration method where the first method is used for the quasistatic (stable) parts of the simulation and the second method for the parts of the simulation that belong to the transient domain. It is shown that the success of the procedure is guaranteed by a proper formulation of the socalled matching conditions that define the transition from one mode of operation to the other. The chosen strategy turns out to be very robust and it has as an added advantage that it can be applied with relative ease. The power of the approach is demonstrated with the presentation of two simulations: The mode jumping problem of a plate strip and the collapse of a thinwalled composite cylinder in compression.
 STAGS User Manual – Version 5.0 , Rhombus Consultants Group, Inc., Palo Alto, CA. Revised January 2005. Previously Report No. LMSC P032594, Lockheed Martin Missles and Space Company, Palo Alto, CA., June 1998. ,

Application of Linear Finite Elements to Finite Strain Using Corotation, AIAA Paper 20061751, 47th AIAA Structures, Structural Dynamics, and Materials Conference, 2006.
,Corotation methods have had great success in extending offtheshelf element kernels to handle seamlessly largedeflection, largerotation response with the limitation that strains must remain moderate. Elementindependent methods promote element reuse by separating necessary operations that filter out rigid rotation and enforce element selfequilibrium from other, elementspecific tasks. The current paper presents a method that preserves element independence while at the same time removes the limitation of moderate strains, based on polar decomposition of the deformation state at the element centroid. The paper will begin with a short review of the foundations of corotation, followed by a detailed derivation of the new method, with close attention to preserving element independence. The efficacy of the method will be demonstrated with several cases of hyperelastic response with huge strains. A comparison will be made between elements with linear straindisplacement kernels and Green’s strain, with almost identical results.

The use of shell elements for the analysis of large strain response, AIAA Paper 20072384, 48th AIAA Structures, Structural Dynamics, and Materials Conference, Hawaii, April 2007.
,Last year we demonstrated that ElementIndependent Corotation based on polar decomposition at the element centroid permits analysis of systems undergoing largestrain response with standard offtheshelf elements originally designed for moderate strains. This capability has implications for a wide variety of problems of special importance to aerospace structures, many of which consist of thin members best modeled by shell elements, where a plane stress and linear throughthethickness strain variation approximate quite closely the actual response. Given the motivation to minimize problem size and computer resource requirements, it is very desirable to have shell elements that can handle large membrane strains and at the same time faithfully reproduce plane stress response in the presence of severe thinning of the shell cross section. In this paper, I shall present our strategy for computing results presented last year for shells. I shall demonstrate superior accuracy and convergence properties for large strain response for a variety of standard shell elements with and without transverse shear deformation.

User’s manual for STAGS, Vol. 1: Theory, NASA Contractor Report 165670, Contract NAS110843, March 1978.
,This is a document for an earlier version of STAGS. Users, especially those at NASA Langley Research Center, have found it useful in addition to the more recent documents listed above.

Structural analysis of general shells: Example cases for STAGSC, Vol. III (preliminary and incomplete), , , January 1977.
,This is a document for an earlier version of STAGS. Users, especially those at NASA Langley Research Center, have found it useful in addition to the more recent documents listed above.
In Memoriam: Dr. Charles C. Rankin (August 1942 – August 2012)
Charles Rankin is known for his salient fundamental contributions over more than 30 years to the field of solid and structural mechanics. All of these contributions have been formulated by Charles and implemented into his generalpurpose nonlinear static and dynamic finite element code called STAGS (STructural Analysis of General Shells), widely used especially at NASA Langley Research Center.
Charles was no mere programmer implementing into a code the mechanics theories of others; he developed these theories mainly by himself. Many of Charles’ fundamental contributions are now finding their way into the most widely used commercial structural computer programs such as MSC_NASTRAN, ANSYS and ABAQUS. In this way Charles’ important original contributions will for the foreseeable future have a major impact on research and engineering in academia, government and industry.
Specifically Charles Rankin’s technical contributions include the formulation and implementation into STAGS of the following:
 A finiteelementindependent corotational theory that has now become an important standard in computational mechanics the world over
 An arclength method that permits the traversal of limit points from prebuckling to postbuckling
 An algorithm to determine multiple bifurcation eigenvalues and eigenvectors from nonlinearly determined prebuckled equilibrium states
 A solution strategy that permits the successive introduction of a sequence of buckling modal imperfections into a nonlinear equilibrium analysis
 A strategy that permits successive smooth transitions from nonlinear static to transient and from nonlinear transient to static analyses
 A strategy that permits the simulation of unzipping of a through crack in a shell possibly with multiple crack tips and turning of a crack during loading
 A “sandwich” finite element that efficiently accounts for soft, sheardeformable cores and stiff face sheets
 Unique nonlinear material models in separate “material modules” which are independent of the rest of the software.
Charles was a great man, loved by many for his great sense of humor, respected by many for his towering intellect and superb accomplishments in the field of structural mechanics. He will remain forever in the minds and hearts of his many friends and colleagues.