The impact of the computer on shell analysis
Thin walled shells structural vibration, buckling and stress analysis by computer methods, discussing shell theory and design
SEARCH · Engineering Papers
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.
Thin walled shells structural vibration, buckling and stress analysis by computer methods, discussing shell theory and design
Flutter boundaries, as well as flutter limit cycle amplitudes, frequencies and stresses were computed for a panel of length-width ratio 4.48 exposed to applied in-plane and transverse loads. The Mach number range was 1.1 to 1.4. The method used involved direct numerical integration of modal equations of motion derived from the nonlinear plate equations of von Karman, coupled with linearized potential flow aerodynamic theory. The flutter boundaries agreed reasonably well with experiment, except when the in-plane loading approached the buckling load. Structural damping had to be introduced, to produce frequencies comparable to the experimental values. Attempts to compute panel deflections or stress at a given point met with limited success. There is some evidence, however, that deflection and stress maxima can be estimated with somewhat greater accuracy.
This paper develops an algorithmic solution strategy which allows the handling of positive/indefinite stiffness characteristics associated with the pre- and post-buckling of structures subject to complex thermomechanical loading fields. The flexibility of the procedure is such that it can be applied to both finite difference and element-type simulations. Due to the generality of the algorithmic approach developed, both kinematic and thermal/mechanical type material nonlinearity including inelastic effects can be treated. This includes the possibility of handling completely general thermomechanical boundary conditions. To demonstrate the scheme, the results of several benchmark problems is presented.
Explore the source record for details and available documents.
The problem of calculating detailed stress distributions around discontinuities in buckled, composite structural components for use with the various analytical failure prediction techniques has not been thoroughly explored. The purpose here is the application of computational methods to the detailed stress analysis problem which is the focus of this session of the workshop. One approach to uncovering the difficulties of this type of analysis and to providing specific directions for future research in this area is a direct attack on the problem using currently available analysis tools. A candidate problem has been selected and experiences from calculating its structural response are described.
The development of the Delft Interactive Shell Design Code (DISDECO) is described. The purpose is to make the accumulated theoretical, numerical, and practical knowledge of the last years readily accessible to users interested in the analysis of buckling sensitive structures. With this open ended, hierarchical, interactive computer code, the user can access workstation programs of increasing complexity. As a first update of DISDECO, the computational modules of the initial level were modified so as to handle combined loads consisting of axial compression, internal and external pressure, and clockwise or counter clockwise torque. The applied load is assumed to have a uniform spatial distribution and is divided into a fixed and a variable part. The magnitude of the variable part is allowed to vary in proportion to a load parameter Lambda. This leads to an eigenvalue problem for the critical lambda sub C.
Multidisciplinary design optimization (MDO) is expected to play a major role in the competitive transportation industries of tomorrow, i.e., in the design of aircraft and spacecraft, of high speed trains, boats, and automobiles. All of these vehicles require maximum performance at minimum weight to keep fuel consumption low and conserve resources. Here, MDO can deliver mathematically based design tools to create systems with optimum performance subject to the constraints of disciplines such as structures, aerodynamics, controls, etc. Although some applications of MDO are beginning to surface, the key to a widespread use of this technology lies in the improvement of its efficiency. This aspect is investigated here for the MDO subset of structural optimization, i.e., for the weight minimization of a given structure under size, strength, and displacement constraints. Specifically, finite element based multilevel optimization of structures (here, statically indeterminate trusses and beams for proof of concept) is performed. In the system level optimization, the design variables are the coefficients of assumed displacement functions, and the load unbalance resulting from the solution of the stiffness equations is minimized. Constraints are placed on the deflection amplitudes and the weight of the structure. In the subsystems level optimizations, the weight of each element is minimized under the action of stress constraints, with the cross sectional dimensions as design variables. This approach is expected to prove very efficient, especially for complex structures, since the design task is broken down into a large number of small and efficiently handled subtasks, each with only a small number of variables. This partitioning will also allow for the use of parallel computing, first, by sending the system and subsystems level computations to two different processors, ultimately, by performing all subsystems level optimizations in a massively parallel manner on separate processors. It is expected that the subsystems level optimizations can be further improved through the use of controlled growth, a method which reduces an optimization to a more efficient analysis with only a slight degradation in accuracy. The efficiency of all proposed techniques is being evaluated relative to the performance of the standard single level optimization approach where the complete structure is weight minimized under the action of all given constraints by one processor and to the performance of simultaneous analysis and design which combines analysis and optimization into a single step. It is expected that the present approach can be expanded to include additional structural constraints (buckling, free and forced vibration, etc.) or other disciplines (passive and active controls, aerodynamics, etc.) for true MDO.
The previous review on stochastic buckling of structures was written by Amazigo in 1976. The present review summarizes some of the developments which took place in recent two decades. A brief overview is given of the effect on uncertainty in the initial geometric imperfections, elastic moduli, applied forces, and thickness variation. For the benefit of the thinking reader, the review has a critical nature. Present essay should be viewed as a direct continuation of our previous paper (1983) with the same title. In order not to repeat what was covered there, it appears instructive to read it although not necessarily prior to dwelling on this article. Accordingly the title is appended with the serial number. It is not promised that the third review will follow since the university science, both fortunately and unfortunately, stands on three things: relevance, interest, and grants.
The results of bending tests on six multiweb beams of optimum weight-strength design are presented. The internal structure of the beams consisted of various combinations of two types of full-depth solid webs and a post-stringer web. The observed structural behavior, buckling load, and failing load of the beams are compared with results obtained by the use of existing methods of analysis and found to be quite predictable.
In this work, an all-bonded out-of-autoclave (OoA) curved longitudinal composite joint concept, intended for use in the next generation of composite heavy lift launch vehicles, was evaluated and verified through finite element (FE) analysis, fabrication, testing, and post-test inspection. The joint was used to connect two curved, segmented, honeycomb sandwich panels representative of a Space Launch System (SLS) fairing design. The overall size of the resultant panel was 1.37 m by 0.74 m (54 in by 29 in), of which the joint comprised a 10.2 cm (4 in) wide longitudinal strip at the center. NASTRAN and ABAQUS were used to perform linear and non-linear analyses of the buckling and strength performance of the jointed panel. Geometric non-uniformities (i.e., surface contour imperfections) were measured and incorporated into the FE model and analysis. In addition, a sensitivity study of the specimens end condition showed that bonding face-sheet doublers to the panel's end, coupled with some stress relief features at corner-edges, can significantly reduce the stress concentrations near the load application points. Ultimately, the jointed panel was subjected to a compressive load. Load application was interrupted at the onset of buckling (at 356 kN 80 kips). A post-test non-destructive evaluation (NDE) showed that, as designed, buckling occurred without introducing any damage into the panel or the joint. The jointed panel was further capable of tolerating an impact damage to the same buckling load with no evidence of damage propagation. The OoA cured all-composite joint shows promise as a low mass factory joint for segmented barrels.
The validation of progressive damage analysis methods requires experimental data that captures the initiation of damage, its evolution with continued loading, and the morphology of damage throughout its evolution. In support of the NASA Advanced Composites Project, several validation specimens of increasing complexity have been tested and characterized in order to provide the data required to assess the capabilities of several progressive damage analysis methods for fiber-reinforced composite materials. The analytical capabilities of interest for this validation program include networks of interacting matrix cracks and delaminations in post-buckled stiffened structures. A series of single stringer compression specimens was tested and characterized at NASA Langley Research Center. The experimental data from these tests are compared with progressive damage finite element analysis results. Test/analysis comparisons are made in terms of load/displacement histories, damage morphology, and damage size versus load metrics.
A parametric study was conducted to gage the sensitivity of simulations of low-velocity impacts on laminated composite panels to various finite element model and material model features. The material model used in this parametric study was previously verified and validated for conducting progressive damage analyses of post-buckled composite structures subjected to quasi-static loading conditions. Based on the results of this parametric study, updates have been made to a progressive damage analysis material model and the best practices for its application. Experimental results are compared to computational results obtained with and without the updates to the models. Both thermoset and thermoplastic fiber-reinforced materials are considered. The resulting enhancements to the material model are presented in detail.
A parametric study was conducted to gage the sensitivity of simulations of low-velocity impacts on laminated composite panels to various finite element model and material model features. The material model used in this parametric study was previously verified and validated for conducting progressive damage analyses of post-buckled composite structures subjected to quasi-static loading conditions. Based on the results of this parametric study, updates have been made to a progressive damage analysis material model and the best practices for its application. Experimental results are compared to computational results obtained with and without the updates to the models. Both thermoset and thermoplastic fiber-reinforced materials are considered. The resulting enhancements to the material model are presented in detail.
DESAP 2 is described as a finite element program for computer-automated, minimum weight design of elastic structures with constraints on stresses (including local instability criteria) and buckling loads. No limits are placed on the number of load conditions for stress-constrained design, but only one of these load conditions can be chosen as the potential buckling load. A substantial portion of DESAP 2, particularly the analysis of the prebuckling state, is derived from the SOLID SAP finite element program. The stress-constrained design is based on the classical stress ratio method, which drives the design towards a fully stressed state. The constraints on the buckling load are handled by solving the appropriate optimality criterion by successive iterations. During each iteration, the element sizes determined by the stress ratio method are used as the minimum size constraints. The element subroutines are organized in a manner that permits the user to make additions and changes with a minimal programming effort. Consequently, DESAP 2 can readily be changed into a special-purpose program to handle the user's specific design requirements and failure criteria.
The feasibility of using the rapidly growing technology of the shape memory alloys actuators in actively controlling the buckling of large flexible structures is investigated. The need for such buckling control systems is becoming inevitable as the design trends of large space structures have resulted in the use of structural members that are long, slender, and very flexible. In addition, as these truss members are subjected mainly to longitudinal loading they become susceptible to structural instabilities due to buckling. Proper control of such instabilities is essential to the effective performance of the structures as stable platforms for communication and observation. Mathematical models are presented that simulate the dynamic characteristics of the shape memory actuator, the compressive structural members, and the associated active control system. A closed-loop computer-controlled system is designed, based on the developed mathematical models, and implemented to control the buckling of simple beams. The performance of the computer-controlled system is evaluated experimentally and compared with the theoretical predictions to validate the developed models. The obtained results emphasize the importance of buckling control and suggest the potential of the shape memory actuators as attractive means for controlling structural deformation in a simple and reliable way.
Structural alloys embody internal mechanisms that allow recovery of state with varying stress and elevated temperature, i.e., they can return to a softer state following periods of hardening. Such material behavior is known to strongly influence structural response under some important thermomechanical loadings, for example, that involving thermal ratchetting. The influence of dynamic and thermal recovery on the creep buckling of a column under variable loading is investigated. The column is taken as the idealized (Shanley) sandwich column. The constitutive model, unlike the commonly employed Norton creep model, incorporates a representation of both dynamic and thermal (state) recovery. The material parameters of the constitutive model are chosen to characterize Narloy Z, a representative copper alloy used in thrust nozzle liners of reusable rocket engines. Variable loading histories include rapid cyclic unloading/reloading sequences and intermittent reductions of load for extended periods of time; these are superimposed on a constant load. The calculated results show that state recovery significantly affects creep buckling under variable loading. Structural alloys embody internal mechanisms that allow recovery of state with varying stress and time.
In the last 20 years NASA has worked in collaboration with industry to develop enabling technologies needed to make aircraft safer and more affordable, extend their lifetime, improve their reliability, better understand their behavior, and reduce their weight. To support these efforts, research programs starting with ideas and culminating in full-scale structural testing were conducted at the NASA Langley Research Center. Each program contained development efforts that (a) started with selecting the material system and manufacturing approach; (b) moved on to experimentation and analysis of small samples to characterize the system and quantify behavior in the presence of defects like damage and imperfections; (c) progressed on to examining larger structures to examine buckling behavior, combined loadings, and built-up structures; and (d) finally moved to complicated subcomponents and full-scale components. Each step along the way was supported by detailed analysis, including tool development, to prove that the behavior of these structures was well-understood and predictable. This approach for developing technology became known as the "building-block" approach. In the Advanced Composites Technology Program and the High Speed Research Program the building-block approach was used to develop a true understanding of the response of the structures involved through experimentation and analysis. The philosophy that if the structural response couldn't be accurately predicted, it wasn't really understood, was critical to the progression of these programs. To this end, analytical techniques including closed-form and finite elements were employed and experimentation used to verify assumptions at each step along the way. This paper presents a discussion of the utilization of the building-block approach described previously in structural mechanics research and development programs at NASA Langley Research Center. Specific examples that illustrate the use of this approach are included from recent research and development programs for both subsonic and supersonic transports.
On December 9th 2009, the International Space Station (ISS) 2A solar array mast experienced prolonged longeron shadowing during a Soyuz undocking. Analytical reconstruction of induced thermal and dynamic structural loads showed an exceedance of the mast buckling limit. Possible structural damage to the solar array mast could have occurred during this event. A Low fidelity video survey of the 2A mast showed no obvious damage of the mast longerons or battens. The decision was made to conduct an on-orbit dynamic test of the 2A array on December 18th, 2009. The test included thruster pluming on the array while photogrammetry data was recorded. The test was similar to other Dedicated Thruster Firings (DTFs) that were performed to measure structural frequency and damping of a solar array. Results of the DTF indicated lower frequency mast modes than model predictions, thus leading to speculation of mast damage. A detailed nonlinear analysis was performed on the 2A array model to assess possible solutions to modal differences. The setup of the parametric nonlinear trade study included the use of a detailed array model and the reduced mass and stiffness matrices of the entire ISS being applied to the array interface. The study revealed that the array attachment structure is nonlinear and thus was the source of error in the model prediction of mast modes. In addition, a detailed study was performed to determine mast mode sensitivity to mast longeron damage. This sensitivity study was performed to assess if the ISS program has sufficient instrumentation for mast damage detection.