SEARCH · Search NASA

Results for “Thick shell finite element formulation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 19 records

All sources NASA NTRS DOE OSTI

NASA NTRS1997

Implementation of Free-Formulation-Based Flat Shell Elements into NASA Comet Code and Development of Nonlinear Shallow Shell Element

This study presents a transient nonlinear finite element analysis within the realm of a multi-body dynamics formulation for determining the dynamic response of a moderately thick laminated shell undergoing a rapid and large rotational motion and nonlinear elastic deformations. Nonlinear strain measure and rotation, as well as 'the transverse shear deformation, are explicitly included in the formulation in order to capture the proper motion-induced stiffness of the laminate. The equations of motion are derived from the virtual work principle. The analysis utilizes a shear deformable shallow shell element along with the co-rotational form of the updated Lagrangian formulation. The shallow shell element formulation is based on the Reissner-Mindlin and Marguerre theory.

Barut, A.↗

NASA NTRS2001

Finite Rotation Analysis of Highly Thin and Flexible Structures

Deployable space structures such as sunshields and solar sails are extremely thin and highly flexible with limited bending rigidity. For analytical investigation of their responses during deployment and operation in space, these structures can be modeled as thin shells. The present work examines the applicability of the solid shell element formulation to modeling of deployable space structures. The solid shell element formulation that models a shell as a three-dimensional solid is convenient in that no rotational parameters are needed for the description of kinematics of deformation. However, shell elements may suffer from element locking as the thickness becomes smaller unless special care is taken. It is shown that, when combined with the assumed strain formulation, the solid shell element formulation results in finite element models that are free of locking even for extremely thin structures. Accordingly, they can be used for analysis of highly flexible space structures undergoing geometrically nonlinear finite rotations.

Clarke, Greg V.↗

NASA NTRS1990

EZDESIT: A Computer Program for Structural Element Sizing and Vehicle Weight Prediction

A method is presented for formulating stiffness terms and thermal coefficients of fiber reinforced composite, stiffened panels for input to finite element analysis (FEA). The method is robust enough to handle panels with general cross sectional shapes, including those which are unsymmetric and/or unbalanced. New thermal coefficients are introduced to quantify panel response from through-the-thickness temperature gradients. Equations are defined for stiffness, thermal expansion, and thermal bending that consider the full complement of membrane, bending, membrane-bending coupling, and transverse shear behavior. A technique of implementing this capability with a single plane of shell finite elements using the MSC/NASTRAN™ FEA program is revealed. Finally, an example of a composite, hat-stiffened panel is included to demonstrate errors that occur when an unsymmetric panel is symmetrically formulated as traditionally done. These erroneous values and the correct ones produced from the presented method are listed for the panel's stiffness matrices, thermal expansion and bending vectors, and the forces and moments that are thermally induced from in-plane and through-the-thickness temperature gradients.

Jeffrey A Cerro↗

NASA NTRS1992

Composite Stiffened Panels: Part I - Stiffness, Thermal Expansion, and Thermal Bending Formulation for Finite Element Analysis

A method is presented for formulating stiffness terms and thermal coefficients of fiber-reinforced composite, stiffened panels for input to finite element analysis (FEA). The method is robust enough to handle panels with general cross sectional shapes, including those which are unsymmetric and/or unbalanced. New thermal coefficients are introduced to quantify panel response from through-the-thickness temperature gradients. Equations are defined for stiffness, thermal expansion, and thermal bending that consider the full complement of membrane, bending, membrane-bending coupling, and transverse shear behavior. A technique of implementing this capability with a single plane of shell finite elements using the MSC/NASTRAN(TradeMark) FEA program is revealed. Finally, an example of a composite, hat-stiffened panel is included to demonstrate errors that occur when an unsymmetric panel is symmetrically formulated as traditionally done. These erroneous values and the correct ones produced from the presented method are listed for the panel's stiffness matrices, thermal expansion and bending vectors, and the forces and moments that are thermally induced from in-plane and through-the-thickness temperature gradients.

Craig S Collier↗

NASA NTRS1986

On the numerical performance of three-dimensional thick shell elements using a hybrid/mixed formulation

Three-dimensional thick shell elements with 8, 16, and 18 nodes are formulated by using the hybrid/mixed method. In bending applications, these elements are free from locking effect and give improved stress predictions. Finite element equations are derived from the Hellinger-Reissner variational principle in which both the displacement and stress fields are approximated by independent interpolation functions. For the assumption of stress parameters, three guidelines are followed: (1) suppression of kinematic deformation modes, (2) invariant element property, and (3) the constraint index exhibited by the element, when applied to constrained-media problems, must be greater than or equal to one. Numerical results are presented to show the element's behavior characteristics regarding sensitivity to locking, distortion effect (patch tests), mesh convergence and the accuracy of stress evaluation.

Graf, W.↗

NASA NTRS1987

A quadrilateral shell element using a mixed formulation

A simple quadrilateral shell element consisting of five nodes, four corner nodes and a central node, is developed for linear elastic analysis of thin as well as moderately thick shells. Based on a modified Hellinger-Reissner principle, finite element equations are derived from the assumed displacement and strain fields. By carefully choosing appropriate strain terms, all kinematic deformation modes are suppressed. Although the present element is similar to a displacement-based degenerated shell, no locking is experienced when it is applied to thin shell problems. Five examples are given to illustrate the analysis capability of the shell element. Numerical results indicate that the element shows fast mesh convergence and gives excellent stress predictions.

Saleeb, A. F.↗

NASA NTRS2001

Finite Element Analysis of Active and Sensory Thermopiezoelectric Composite Materials

Analytical formulations are developed to account for the coupled mechanical, electrical, and thermal response of piezoelectric composite materials. The coupled response is captured at the material level through the thermopiezoelectric constitutive equations and leads to the inherent capability to model both the sensory and active responses of piezoelectric materials. A layerwise laminate theory is incorporated to provide more accurate analysis of the displacements, strains, stresses, electric fields, and thermal fields through-the-thickness. Thermal effects which arise from coefficient of thermal expansion mismatch, pyroelectric effects, and temperature dependent material properties are explicitly accounted for in the formulation. Corresponding finite element formulations are developed for piezoelectric beam, plate, and shell elements to provide a more generalized capability for the analysis of arbitrary piezoelectric composite structures. The accuracy of the current formulation is verified with comparisons from published experimental data and other analytical models. Additional numerical studies are also conducted to demonstrate additional capabilities of the formulation to represent the sensory and active behaviors. A future plan of experimental studies is provided to characterize the high temperature dynamic response of piezoelectric composite materials.

Lee, Ho-Jun↗

NASA NTRS1993

Stiffness, thermal expansion, and thermal bending formulation of stiffened, fiber-reinforced composite panels

A method is presented for formulating stiffness terms and thermal coefficients of stiffened, fiber-reinforced composite panels. The method is robust enough to handle panels with general cross sectional shapes, including those which are unsymmetric and/or unbalanced. Nonlinear, temperature and load dependent constitutive material data of each laminate are used to 'build-up' the stiffened panel membrane, bending, and membrane-bending coupling stiffness terms and thermal coefficients. New thermal coefficients are introduced to quantify panel response from through-the-thickness temperature gradients. A technique of implementing this capability with a single plane of shell finite elements using the MSC/NASTRAN analysis program (FEA) is revealed that provides accurate solutions of entire airframes or engines with coarsely meshed models. An example of a composite, hat-stiffened panel is included to demonstrate errors that occur when an unsymmetric panel is symmetrically formulated as traditionally done. The erroneous results and the correct ones produced from this method are compared to analysis from discretely meshed three-dimensional FEA.

Collier, Craig S.↗

NASA NTRS1994

Discrete-Layer Piezoelectric Plate and Shell Models for Active Tip-Clearance Control

The objectives of this work were to develop computational tools for the analysis of active-sensory composite structures with added or embedded piezoelectric layers. The targeted application for this class of smart composite laminates and the analytical development is the accomplishment of active tip-clearance control in turbomachinery components. Two distinct theories and analytical models were developed and explored under this contract: (1) a discrete-layer plate theory and corresponding computational models, and (2) a three dimensional general discrete-layer element generated in curvilinear coordinates for modeling laminated composite piezoelectric shells. Both models were developed from the complete electromechanical constitutive relations of piezoelectric materials, and incorporate both displacements and potentials as state variables. This report describes the development and results of these models. The discrete-layer theories imply that the displacement field and electrostatic potential through-the-thickness of the laminate are described over an individual layer rather than as a smeared function over the thickness of the entire plate or shell thickness. This is especially crucial for composites with embedded piezoelectric layers, as the actuating and sensing elements within these layers are poorly represented by effective or smeared properties. Linear Lagrange interpolation polynomials were used to describe the through-thickness laminate behavior. Both analytic and finite element approximations were used in the plane or surface of the structure. In this context, theoretical developments are presented for the discrete-layer plate theory, the discrete-layer shell theory, and the formulation of an exact solution for simply-supported piezoelectric plates. Finally, evaluations and results from a number of separate examples are presented for the static and dynamic analysis of the plate geometry. Comparisons between the different approaches are provided when possible, and initial conclusions regarding the accuracy and limitations of these models are given.

Heyliger, P. R.↗

NASA NTRS2003

A Variational Principle for Reconstruction of Elastic Deformations in Shear Deformable Plates and Shells

A variational principle is formulated for the inverse problem of full-field reconstruction of three-dimensional plate/shell deformations from experimentally measured surface strains. The formulation is based upon the minimization of a least squares functional that uses the complete set of strain measures consistent with linear, first-order shear-deformation theory. The formulation, which accommodates for transverse shear deformation, is applicable for the analysis of thin and moderately thick plate and shell structures. The main benefit of the variational principle is that it is well suited for C(sup 0)-continuous displacement finite element discretizations, thus enabling the development of robust algorithms for application to complex civil and aeronautical structures. The methodology is especially aimed at the next generation of aerospace vehicles for use in real-time structural health monitoring systems.

Tessler, Alexander↗

NASA NTRS1986

On a third-order shear deformation theory for laminated composite shells

A higher-order theory based on an assumed displacement field in which the surface displacements are expanded in powers of the thickness coordinate up to the third order is presented. The theory allows parabolic description of the transverse shear stresses, and therefore the shear correction factors of the usual shear deformation theory are not required in the present theory. The theory also accounts for small strains but moderately large displacements (i.e., von Karman strains). A finite-element model based on independent approximations of the displacements and bending moments (i.e., mixed formulation) is developed. The element is used to analyze cross-ply and angle-ply laminated shells for bending.

Liu, C. F.↗

NASA NTRS1988

Hierarchic plate and shell models based on p-extension

This paper is concerned with formulations of finite element models for beams, arches, plates and shells based on the principle of virtual work. The focus is on computer implementation of hierarchic sequences of finite element models suitable for numerical solution of a large variety of practical problems which may concurrently contain thin and thick plates and shells, stiffeners, and regions where three dimensional representation is required. The approximate solutions corresponding to the hierarchic sequence of models converge to the exact solution of the fully three dimensional model. The stopping criterion is based on (1) estimation of the relative error in energy norm; (2) equilibrium tests; and (3) observation of the convergence of quantities of interest.

Szabo, Barna A.↗

NASA NTRS1988

Hierarchic plate and shell models based on p-extension

Formulations of finite element models for beams, arches, plates and shells based on the principle of virtual work was studied. The focus is on computer implementation of hierarchic sequences of finite element models suitable for numerical solution of a large variety of practical problems which may concurrently contain thin and thick plates and shells, stiffeners, and regions where three dimensional representation is required. The approximate solutions corresponding to the hierarchic sequence of models converge to the exact solution of the fully three dimensional model. The stopping criterion is based on: (1) estimation of the relative error in energy norm; (2) equilibrium tests, and (3) observation of the convergence of quantities of interest.

Szabo, Barna A.↗

NASA NTRS2018

An Efficient Inverse Finite Element Method for Shape and Stress Sensing of Laminated Composite and Sandwich Plates and Shells

The inverse Finite Element Method (iFEM) is a revolutionary methodology for real-time reconstruction of full-field structural displacements and stresses in structures that are instrumented with strain sensors. This inverse problem is commonly referred to as shape and stress sensing, which is well-recognized as an enabling technology for structural health monitoring systems. In this study, an improved iFEM formulation is proposed for shape and stress sensing of laminated composite and sandwich plates and shells. The formulation includes the kinematics of a shear deformation plate theory known as Refined Zigzag Theory (RZT) as its baseline. The present iFEM formulation is based upon the minimization of a weighted-least-squares functional that uses the complete set of section strains of RZT. The improved iFEM methodology is applicable for shape and stress sensing of thin and moderately thick plate and shell structures involving a relatively small number of strain gauges. The main advantage of the current formulation is that highly accurate through-the-thickness distributions of displacements, strains, and stresses are attainable using an element based on simple C0-continuous displacement interpolation functions. A three-node inverseshell element, named i3-RZT, is developed. Two example problems are examined in detail: (1) a simply supported rectangular laminated composite plate and (2) a wedge structure with a hole near one of the clamped ends. For both problems, the experimental strain data are generated numerically by the direct finite element analysis using high-fidelity discretizations. These strains are then regarded as the experimental strains obtained from surface mounted strain gauges or embedded fiber Bragg grating (FBG) sensors. The numerical results demonstrate the superior capability and potential applicability of the i3- RZT/iFEM methodology for performing accurate shape and stress sensing of complex composite structures.

Kefal, Adnan↗

NASA NTRS1988

An efficient Mindlin finite strip plate element based on assumed strain distribution

A simple two node, linear, finite strip plate bending element based on Mindlin-Reissner plate theory for the analysis of very thin to thick bridges, plates, and axisymmetric shells is presented. The new transverse shear strains are assumed for constant distribution in the two node linear strip. The important aspect is the choice of the points that relate the nodal displacements and rotations through the locking transverse shear strains. The element stiffness matrix is explicitly formulated for efficient computation and ease in computer implementation. Numerical results showing the efficiency and predictive capability of the element for analyzing plates with different supports, loading conditions, and a wide range of thicknesses are given. The results show no sign of the shear locking phenomenon.

Chulya, Abhisak↗

NASA NTRS1988

A preliminary investigation of finite-element modeling for composite rotor blades

The results from an initial phase of an in-house study aimed at improving the dynamic and aerodynamic characteristics of composite rotor blades through the use of elastic couplings are presented. Large degree of freedom shell finite element models of an extension twist coupled composite tube were developed and analyzed using MSC/NASTRAN. An analysis employing a simplified beam finite element representation of the specimen with the equivalent engineering stiffness was additionally performed. Results from the shell finite element normal modes and frequency analysis were compared to those obtained experimentally, showing an agreement within 13 percent. There was appreciable degradation in the frequency prediction for the torsional mode, which is elastically coupled. This was due to the absence of off-diagonal coupling terms in the formulation of the equivalent engineering stiffness. Parametric studies of frequency variation due to small changes in ply orientation angle and ply thickness were also performed. Results showed linear frequency variations less than 2 percent per 1 degree variation in the ply orientation angle, and 1 percent per 0.0001 inch variation in the ply thickness.

Lake, Renee C.↗

NASA NTRS1988

A preliminary investigation of finite-element modeling for composite rotor blades

Lake, Renee C.↗

NASA NTRS1991

Probabilistic Structures Analysis Methods (PSAM) for select space propulsion system components

The basic formulation for probabilistic finite element analysis is described and demonstrated on a few sample problems. This formulation is based on iterative perturbation that uses the factorized stiffness on the unperturbed system as the iteration preconditioner for obtaining the solution to the perturbed problem. This approach eliminates the need to compute, store and manipulate explicit partial derivatives of the element matrices and force vector, which not only reduces memory usage considerably, but also greatly simplifies the coding and validation tasks. All aspects for the proposed formulation were combined in a demonstration problem using a simplified model of a curved turbine blade discretized with 48 shell elements, and having random pressure and temperature fields with partial correlation, random uniform thickness, and random stiffness at the root.

Source record↗