A finite element approach to the structural instability of beam columns, frames, and arches
Stiffness matrix for structural stability analyses on beam columns, frames, and arches
SEARCH · Search NASA
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.
Stiffness matrix for structural stability analyses on beam columns, frames, and arches
The results of research in applying statistical methods to the problem of structural dynamic system identification are presented. The study is in three parts: a review of previous approaches by other researchers, a development of various linear estimators which might find application, and the design and development of a computer program which uses a Bayesian estimator. The method is tried on two models and is successful where the predicted stiffness matrix is a proper model, e.g., a bending beam is represented by a bending model. Difficulties are encountered when the model concept varies. There is also evidence that nonlinearity must be handled properly to speed the convergence.
An algorithm is presented to relabel automatically the nodes of an arbitrary finite-element mesh. The purpose of such relabeling is to reduce the bandwidth of the master stiffness matrix produced by the finite-element method. The algorithm uses a random process for the relabeling. Computing time is reduced substantially, compared to systematic methods.
A self-correcting, incremental-solution procedure has been developed which seems to be promising in solving geometrically highly nonlinear problems. The procedure is also applicable to systems with many degrees of freedom where the bandwidth of the stiffness matrix is large compared with the narrow bandwidth for the shell of revolution.
A unique and relatively simple approach is presented for obtaining the linearized equations of motion. A conservative system consisting of two rigid bodies connected by any number of massless cables with linear axial stiffness is assumed. The cable forces and torques are expanded in a Taylor series about the equilibrium values of the system coordinates which results in a cable stiffness matrix. The method of obtaining the equilibrium values of the coordinates is discussed and results are presented. The range of validity of the linear model is determined by comparing results with a digital simulation of the nonlinear system.
Using stiffness matrix and lumped-mass representation specified number of natural frequencies are obtained using inverse iteration method. Mode shapes for each frequency are also obtained. These frequencies and mode shapes can be found in reasonable periods of computer time utilizing this code.
The elastostatic interaction problem between a penny-shaped crack and a slender inclusion or filament in an elastic matrix was formulated. For a single filament as well as multiple identical filaments located symmetrically around the crack the problem is shown to reduce to a singular integral equation. The solution of the problem is obtained for various geometries and filament-to-matrix stiffness ratios, and the results relating to the angular variation of the stress intensity factor and the maximum filament stress are presented.
The formulation of the linear stiffness matrix for a doubly-curved triangular thin shell element, using a modified potential energy principle, is described. The strain energy component of the potential energy is expressed in terms of displacements and displacement gradients by use of consistent Koiter strain-displacement equations. The element inplane and normal displacement fields are approximated by complete cubic polynomials. The interelement displacement admissibility conditions are met in the global representation by imposition of constraint conditions on the interelement boundaries; the constraints represent the modification of the potential energy. Errors due to the nonzero strains under rigid body motion are shown to be of small importance for practical grid refinements through performance of extensive comparison analyses.
Aspects of the formulation of a triangular thin shell finite element which pertain to geometrically nonlinear (small strain, finite displacement) behavior are described. The procedure for solution of the resulting nonlinear algebraic equations combines a one-step incremental (tangent stiffness) approach with one iteration in the Newton-Raphson mode. A method is presented which permits a rational estimation of step size in this procedure. Limit points are calculated by means of a superposition scheme coupled to the incremental side of the solution procedure while bifurcation points are calculated through a process of interpolation of the determinants of the tangent-stiffness matrix. Numerical results are obtained for a flat plate and two curved shell problems and are compared with alternative solutions.
A capability for the nonlinear vibration analysis of beam and frame structures suitable for use with NASTRAN level 15.5 is described. The nonlinearity considered is due to the presence of axial loads induced by longitudinal end restraints and lateral displacements that are large compared to the beam height. A brief discussion is included of the mathematical analysis and the geometrical stiffness matrix for a prismatic beam (BAR) element. Also included are a brief discussion of the equivalent linearization iterative process used to determine the nonlinear frequency, the required modifications to subroutines DBAR and XMPLBD of the NASTRAN code, and the appropriate vibration capability, four example problems are presented. Comparisons with existing experimental and analytical results show that excellent accuracy is achieved with NASTRAN in all cases.
The problem of case neck separation malfunction in the field of ammunition structural analysis is investigated. The axi-symmetric solid of revolution RING element was utilized in the manual piecewise linear analysis to obtain the expansion of the wall of the cartridge case and barrel chamber by the pressure of propellant gases and the stresses in the structure. The analysis included the varying material properties along the wall of the case and the chamber. Additional instructions were provided to change the element material ID's without recomputing the entire stiffness matrix.
A quadrilateral shell element, CQUAD4, was added to level 15.5 and subsequently to level 16.0 of NASTRAN. The element exhibited doubly curved surfaces and used biquadratic interpolation functions. Reduced integration techniques were used to improve the performance of the element in thin shell problems. The creation of several new bulk data items is discussed, along with a special module, GPNORM, to process SHLNORM bulk data cards. In addition to the theoretical basis for the element stiffness matrix, consistent mass and load matrices are presented. Several potential sources of degenerate behavior of the element were investigated. Guidelines for proper use of the element were suggested. Performance of the element on several widely published classical examples was demonstrated. The results showed a significant improvement over presently available NASTRAN shell elements for even the coarsest meshes. Potential applications to two classes of practical problems are discussed.
A modal residual-flexibility approach due to Schwendler and MacNeal (1962) is adapted to account for the 'static' contribution of higher-frequency modes without requiring their integration. It is assumed that each rotor mode acts as a lightly damped second-order system. The additional accuracy provided by residual-flexibility corrections becomes progressively more important as a modal model's actual boundary conditions are forced to deviate from the boundary conditions used to define the rotor's original structural model (stiffness matrix) and its associated eigendata input to the transient modal model. An analysis of the high-pressure-oxygen turbopump of the Space Shuttle main engine shows that the residual-flexibility approach ensures a substantial improvement in accuracy for a relatively moderate increase in computer-time requirements.
The effect of debonds on composite structures under compressive loads is investigated in this paper. Because of the inclusion of debonds, unsymmetric laminated rectangular plates are considered. The resulting bending-extension coupling is included and is shown to be quite significant for thin laminates typical of debonded regions. The solution technique is based on the Theorem of Minimum Potential Energy; and both simply supported and clamped boundaries are considered. The need to include transverse shear deformation in the analysis is documented. Also, the range of applicability of the reduced stiffness matrix technique for such plates is stated. The consideration of non-rectangular boundaries, as would normally occur for debonded regions, has been included through a finite element study. Using the results of this study, an effective laminated plate theory for irregularly shaped debonded regions is proposed.
The nonlinear behavior of a wave like deformed metal strip during the levelling process were calculated. Elastic-plastic material behavior as well as nonlinearities due to large deformations were considered. The considered problem lead to a combined stability and contact problem. It is shown that, despite the initially concentrated loading, neglecting the change of loading conditions due to altered contact domains may lead to a significant error in the evaluation of the nonlinear behavior and particularly to an underestimation of the stability limit load. The stability was examined by considering the load deflection path and the behavior of a load-dependent current stiffness parameter in combination with the determinant of the current stiffness matrix.
Equations are developed for the buckling of a general lattice structure that has repetitive geometry. Equilibrium at a typical node is expressed using finite element techniques, and the only assumption is that the response is periodic. By basing the stiffness matrix on the exact solution of the beam column equation, accurate results are obtained for complex buckling behavior that would require a very large system of equations using conventional techniques. The present method requires the eigenvalues of only a 6x6 determinant. The results are used to study the buckling of isogrid cylinders, three-element truss columns and polygonal rings. Details of the analysis including expressions for all terms in the governing stability determinant are given.
The effects of solar radiation pressure as the main environmental disturbance torque were incorporated into the model of the rigid orbiting shallow shell and computer simulation results indicate that within the linear range the rigid modal amplitudes are excited in proportion to the area to mass ratio. The effect of higher order terms in the gravity-gradient torque expressions previously neglected was evaluated and found to be negligible for the size structures under consideration. A graph theory approach was employed for calculating the eigenvalues of a large flexible system by reducing the system (stiffness) matrix to lower ordered submatrices. The related reachability matrix and term rank concepts are used to verify controllability and can be more effective than the alternate numerical rank tests. Control laws were developed for the shape and orientation control of the orbiting flexible shallow shell and numerical results presented.
Equations are developed for vibration of general lattice structures that have repetitive geometry. The method of solution is an extension of a previous paper for buckling of similar structures. The theory is based on representing each member of the structure with the exact dynamic stiffness matrix and taking advantage of the repetitive geometry to obtain an eigenvalue problem involving the degrees-of-freedom at a single node in the lattice. Results are given for shell-and beam-like lattice structures and for rings stiffened with tension cables and a central mast. The variation of frequency with external loading and the effect of local member vibration on overall modes is shown.