Search NASA⌕ Search

SEARCH · Search NASA

Results for “shell elements”

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.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Developing Procedures to Implement Geometric Imperfections Beyond Right Circular Cylindrical Shells in Finite Element Method Models

Analysis of aerospace structures is frequently conducted using nominal dimensions and frequently assumes ideal conditions in loading, contact, constraints, et cetera. Off-nominal dimensions and nonideal conditions, however, are present in all structures. These are the result of widely ranging causes from coefficient of thermal expansion mismatches, manufacturing tooling anomalies, to assembly procedures that inadvertently alter the structure. Specifically, geometric imperfections can have potentially significant influence on the response of a structural test article observed in an experiment versus the response given by a numerical simulation. The Python Tool for Implementing Geometric Imperfections in Reduced Structures (Py_TIGIRS) was previously presented as a set of Python scripts to calculate and implement as-manufactured geometric midsurface and thickness imperfections into finite element method (FEM) shell models of nominally right circular cylinders. By taking advantage of the simple shape of a right circular cylinder, interpolations of the measured data points were able to be performed along directions that aligned to the cylindrical coordinate system axes of the entire structure. By taking advantage of the shell representation of the real structure as opposed to modeling using a continuum representation, the thickness variation was able to be implemented by shell section definitions instead of having to modify the position of multiple nodes in the thickness direction. Py_TIGIRS is a useful tool that established a procedural example on how to implement geometric imperfections in right circular cylindrical shell structures. Three new procedures, each expanded from concepts established in Py_TIGIRS, are proposed for various test-article designs and are intended to broaden the range of structures that can be modeled with measured geometric imperfections in the structural analysis community. Each test-article design introduces new challenges to successfully implement geometric imperfections into a FEM model. The first test-article design consists of a carbon fiber reinforced polymer square plate with a hat-shaped stiffener co-cured on one side. This test-article design was for a novel seven-point bend test that was also previously presented. Manufacturing and cure-cycle imperfections are observed using digital image correlation (DIC) techniques. As thermal expansion coefficient mismatches between the plate and stiffener materials were anticipated, a thermal analysis study with continuum shell and solid elements was conducted to capture the global shape observed prior to testing. The second test-article design is of a similar hat-stiffened plate configuration, but with a side length ratio near 3:1 with elongation in the stiffener direction. The test article was used to characterize the response to uniaxial compressive loading in the direction of the stiffener. Due to differing manufacturing steps, a thermal analysis like the one developed for the seven-point bend configuration was unable to mimic the observed geometric imperfections. Instead, a strategy based on applying deformations directly to the structure during analysis was developed for continuum shell and solid element representation of a stiffened panel.

Geometric imperfections↗
NASA NTRS2012

An Efficient Analysis Methodology for Fluted-Core Composite Structures

The primary loading condition in launch-vehicle barrel sections is axial compression, and it is therefore important to understand the compression behavior of any structures, structural concepts, and materials considered in launch-vehicle designs. This understanding will necessarily come from a combination of test and analysis. However, certain potentially beneficial structures and structural concepts do not lend themselves to commonly used simplified analysis methods, and therefore innovative analysis methodologies must be developed if these structures and structural concepts are to be considered. This paper discusses such an analysis technique for the fluted-core sandwich composite structural concept. The presented technique is based on commercially available finite-element codes, and uses shell elements to capture behavior that would normally require solid elements to capture the detailed mechanical response of the structure. The shell thicknesses and offsets using this analysis technique are parameterized, and the parameters are adjusted through a heuristic procedure until this model matches the mechanical behavior of a more detailed shell-and-solid model. Additionally, the detailed shell-and-solid model can be strategically placed in a larger, global shell-only model to capture important local behavior. Comparisons between shell-only models, experiments, and more detailed shell-and-solid models show excellent agreement. The discussed analysis methodology, though only discussed in the context of fluted-core composites, is widely applicable to other concepts.

Oremont, Leonard↗
NASA NTRS1999

Analytical Prediction of Damage Growth in Notched Composite Panels Loaded in Axial Compression

A progressive failure analysis method based on shell elements is developed for the computation of damage initiation and growth in stiffened thick-skin stitched graphite-epoxy panels loaded in axial compression. The analysis method involves a step-by-step simulation of material degradation based on ply-level failure mechanisms. High computational efficiency is derived from the use of superposed layers of shell elements to model each ply orientation in the laminate. Multiple integration points through the thickness are used to obtain the correct bending effects through the thickness without the need for ply-by-ply evaluations of the state of the material. The analysis results are compared with experimental results for three stiffened panels with notches oriented at 0, 15 and 30 degrees to the panel width dimension. A parametric study is performed to investigate the damage growth retardation characteristics of the Kevlar stitch lines in the pan

Ambur, Damodar R.↗
NASA NTRS2004

Analysis of Composite Skin-Stiffener Debond Specimens Using a Shell/3D Modeling Technique and Submodeling

The application of a shell/3D modeling technique for the simulation of skin/stringer debond in a specimen subjected to tension and three-point bending was studied. The global structure was modeled with shell elements. A local three-dimensional model, extending to about three specimen thicknesses on either side of the delamination front was used to model the details of the damaged section. Computed total strain energy release rates and mixed-mode ratios obtained from shell/3D simulations were in good agreement with results obtained from full solid models. The good correlation of the results demonstrated the effectiveness of the shell/3D modeling technique for the investigation of skin/stiffener separation due to delamination in the adherents. In addition, the application of the submodeling technique for the simulation of skin/stringer debond was also studied. Global models made of shell elements and solid elements were studied. Solid elements were used for local submodels, which extended between three and six specimen thicknesses on either side of the delamination front to model the details of the damaged section. Computed total strain energy release rates and mixed-mode ratios obtained from the simulations using the submodeling technique were not in agreement with results obtained from full solid models.

OBrien, T. Kevin↗
NASA NTRS1985

Nonlinear analysis of laminated shells including transverse shear strains

Numerical results obtained using a doubly curved, shear deformable shell element are presented for geometrically nonlinear analysis of laminated composite shells. The element is based on an extension of Sanders' shell theory and accounts for the von Karman strains and transverse shear strains. The sample numerical results presented here for the geometrically nonlinear analysis of laminated composite shells should serve as references for future investigations.

Reddy, J. N.↗
NASA NTRS1985

Application of the finite element method in the calculation of transmission loss of flat and curved panels

This investigation represents an extension of a study of Roussos (1985) who considered the noise transmission loss of a rectangular plate in an infinite baffle. Roussos, who employed an analytical formulation, considered an unstiffened plate. While it is difficult to consider stiffeners by means of analytical methods, the difficulties can be avoided by employing a finite element procedure. For this reason, the present study is concerned with the implementation of a finite element method. The representation of the panel transmission loss is discussed, and the determination of the panel motion by means of the finite element technique is described, taking into account an isotropic flat panel, the exciting force, an eigenvalue problem, the radiation pressure, a plate element, and a cylindrical shell element. Numerical results are considered for a flat panel, a curved panel, and a stiffened flat panel.

Koval, L. R.↗
NASA NTRS1992

Beam, plate, and shell finite elements for CSI design code

Efficient triangular plate bending finite elements (with 6 degree-of-freedom per node) are surveyed. The most promising element will be selected, coded, and implemented into the current Control Structures Interactions (CSI) Design Code. Numerical example, v.i.a. CSI Design Code will be conducted to confirm the efficiency and accuracy of the selected plate bending element. Procedures to further enhance the computations of the CSI design code will be recommended.

Nguyen, Duc T.↗
NASA NTRS1993

Fluid-structural interactions using Navier-Stokes flow equations coupled with shell finite element structures

A computational procedure is presented to study fluid-structural interaction problems for three-dimensional aerospace structures. The flow is modeled using the three-dimensional unsteady Euler/Navier-Stokes equations and solved using the finite-difference approach. The three dimensional structure is modeled using shell/plate finite-element formulation. The two disciplines are coupled using a domain decomposition approach. Accurate procedures both in time and space are developed to combine the solutions from the flow equations with those of the structural equations. Time accuracy is maintained using aeroelastic configuration-adaptive moving grids that are computed every time step. The work done by aerodynamic forces due to structural deformations is preserved using consistent loads. The present procedure is validated by computing the aeroelastic response of a wing and comparing with experiment. Results are illustrated for a typical wing-body configuration.

Guruswamy, Guru P.↗
NASA NTRS2007

Simulating Space Capsule Water Landing with Explicit Finite Element Method

A study of using an explicit nonlinear dynamic finite element code for simulating the water landing of a space capsule was performed. The finite element model contains Lagrangian shell elements for the space capsule and Eulerian solid elements for the water and air. An Arbitrary Lagrangian Eulerian (ALE) solver and a penalty coupling method were used for predicting the fluid and structure interaction forces. The space capsule was first assumed to be rigid, so the numerical results could be correlated with closed form solutions. The water and air meshes were continuously refined until the solution was converged. The converged maximum deceleration predicted is bounded by the classical von Karman and Wagner solutions and is considered to be an adequate solution. The refined water and air meshes were then used in the models for simulating the water landing of a capsule model that has a flexible bottom. For small pitch angle cases, the maximum deceleration from the flexible capsule model was found to be significantly greater than the maximum deceleration obtained from the corresponding rigid model. For large pitch angle cases, the difference between the maximum deceleration of the flexible model and that of its corresponding rigid model is smaller. Test data of Apollo space capsules with a flexible heat shield qualitatively support the findings presented in this paper.

Wang, John T.↗
NASA NTRS2015

Development and Characterization of a Rate-Dependent Three-Dimensional Macroscopic Plasticity Model Suitable for Use in Composite Impact Problems

Several key capabilities have been identified by the aerospace community as lacking in the material/models for composite materials currently available within commercial transient dynamic finite element codes such as LS-DYNA. Some of the specific desired features that have been identified include the incorporation of both plasticity and damage within the material model, the capability of using the material model to analyze the response of both three-dimensional solid elements and two dimensional shell elements, and the ability to simulate the response of composites composed with a variety of composite architectures, including laminates, weaves and braids. In addition, a need has been expressed to have a material model that utilizes tabulated experimentally based input to define the evolution of plasticity and damage as opposed to utilizing discrete input parameters (such as modulus and strength) and analytical functions based on curve fitting. To begin to address these needs, an orthotropic macroscopic plasticity based model suitable for implementation within LS-DYNA has been developed. Specifically, the Tsai-Wu composite failure model has been generalized and extended to a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The coefficients in the yield function are determined based on tabulated stress-strain curves in the various normal and shear directions, along with selected off-axis curves. Incorporating rate dependence into the yield function is achieved by using a series of tabluated input curves, each at a different constant strain rate. The non-associative flow-rule is used to compute the evolution of the effective plastic strain. Systematic procedures have been developed to determine the values of the various coefficients in the yield function and the flow rule based on the tabulated input data. An algorithm based on the radial return method has been developed to facilitate the numerical implementation of the material model. The presented paper will present in detail the development of the orthotropic plasticity model and the procedures used to obtain the required material parameters. Methods in which a combination of actual testing and selective numerical testing can be combined to yield the appropriate input data for the model will be described. A specific laminated polymer matrix composite will be examined to demonstrate the application of the model.

Impact↗
NASA NTRS1973

The finite element method in shell stability analysis.

A development of the finite element method for thin shell instability analysis is presented, covering three principal aspects: (1) representation of shell geometry, (2) representation of element behavior, and (3) algorithmic tools for solution of the large-order systems of nonlinear algebraic equations which characterize various phases of shell instability. Two shell elements are described, an arbitrary quadrilateral and a triangle, and numerical results are presented for two widely-employed comparison problems for linear (stable) analysis. Two shell problems which include instability effects are also solved.-

Gallagher, R. H.↗
NASA NTRS1988

A finite element study of the EIDI system

This paper presents a method for modeling the structural dynamics of an Electro-Impulse De-Icing System, using finite element analyses procedures. A guideline for building a representative finite element model is discussed. Modeling was done initially using four noded cubic elements, four noded isoparametric plate elements and eight noded isoparametric shell elements. Due to the size of the problem and due to the underestimation of shear stress results when compared to previous analytical work an approximate model was created to predict possible areas of shedding of ice. There appears to be good agreement with the test data provided by The Boeing Commercial Airplane Company. Thus these initial results of this method were found to be encouraging. Additional analytical work and comparison with experiment is needed in order to completely evaluate this approach.

Khatkhate, A. A.↗
NASA NTRS1993

Current capabilities for simulating the extreme distortion of thin structures subjected to severe impacts

The explicit transient dynamics technology in use today for simulating the impact and subsequent transient dynamic response of a structure has its origins in the 'hydrocodes' dating back to the late 1940's. The growth in capability in explicit transient dynamics technology parallels the growth in speed and size of digital computers. Computer software for simulating the explicit transient dynamic response of a structure is characterized by algorithms that use a large number of small steps. In explicit transient dynamics software there is a significant emphasis on speed and simplicity. The finite element technology used to generate the spatial discretization of a structure is based on a compromise between completeness of the representation for the physical processes modelled and speed in execution. That is, since it is expected in every calculation that the deformation will be finite and the material will be strained beyond the elastic range, the geometry and the associated gradient operators must be reconstructed, as well as complex stress-strain models evaluated at every time step. As a result, finite elements derived for explicit transient dynamics software use the simplest and barest constructions possible for computational efficiency while retaining an essential representation of the physical behavior. The best example of this technology is the four-node bending quadrilateral derived by Belytschko, Lin and Tsay. Today, the speed, memory capacity and availability of computer hardware allows a number of the previously used algorithms to be 'improved.' That is, it is possible with today's computing hardware to modify many of the standard algorithms to improve their representation of the physical process at the expense of added complexity and computational effort. The purpose is to review a number of these algorithms and identify the improvements possible. In many instances, both the older, faster version of the algorithm and the improved and somewhat slower version of the algorithm are found implemented together in software. Specifically, the following seven algorithmic items are examined: the invariant time derivatives of stress used in material models expressed in rate form; incremental objectivity and strain used in the numerical integration of the material models; the use of one-point element integration versus mean quadrature; shell elements used to represent the behavior of thin structural components; beam elements based on stress-resultant plasticity versus cross-section integration; the fidelity of elastic-plastic material models in their representation of ductile metals; and the use of Courant subcycling to reduce computational effort.

Key, Samuel W.↗
NASA NTRS2001

Notch Sensitivity of Woven Ceramic Matrix Composites Under Tensile Loading: An Experimental, Analytical, and Finite Element Study

The stress concentrations associated with circular notches and subjected to uniform tensile loading in woven ceramic matrix composites (CMCs) have been investigated for high-efficient turbine engine applications. The CMC's were composed of Nicalon silicon carbide woven fabric in SiNC matrix manufactured through polymer impregnation process (PIP). Several combinations of hole diameter/plate width ratios and ply orientations were considered in this study. In the first part, the stress concentrations were calculated measuring strain distributions surrounding the hole using strain gages at different locations of the specimens during the initial portion of the stress-strain curve before any microdamage developed. The stress concentration was also calculated analytically using Lekhnitskii's solution for orthotropic plates. A finite-width correction factor for anisotropic and orthotropic composite plate was considered. The stress distributions surrounding the circular hole of a CMC's plate were further studied using finite element analysis. Both solid and shell elements were considered. The experimental results were compared with both the analytical and finite element solutions. Extensive optical and scanning electron microscopic examinations were carried out for identifying the fracture behavior and failure mechanisms of both the notched and notched specimens. The stress concentration factors (SCF) determined by analytical method overpredicted the experimental results. But the numerical solution underpredicted the experimental SCF. Stress concentration factors are shown to increase with enlarged hole size and the effects of ply orientations on stress concentration factors are observed to be negligible. In all the cases, the crack initiated at the notch edge and propagated along the width towards the edge of the specimens.

Haque, A.↗
NASA NTRS2021

Application of a Developmental Composite Material Model to Predict the Crush Response of Two Energy Absorbers

In 2012, a consortium was formed with the goal of creating a new composite material model capable of predicting the wide range of properties, accumulated damage behavior, and the many different types of failure in composites under impact loading. The material model was developed for execution in the commercially available nonlinear, explicit transient dynamic finite element code, LS-DYNA®. This material model incorporates three sub-models: deformation, damage, and failure. In addition, the model accounts for strain rate and temperature effects and relies heavily on the input of tabulated material response data. The model is designated *MAT_COMPOSITE_TABULATED_PLASTICITY_DAMAGE, or *MAT_213. Initially, *MAT_213 was developed for use with solid elements only; however, a thin shell element formulation for *MAT_213 has been adapted recently. The objective of this project was to find a suitable modeling example to investigate the capabilities and performance of *MAT_213. In 2012, two composite energy absorbers were designed and evaluated at NASA Langley Research Center through multi-level testing and simulation. The first was a conical-shaped energy absorber, designated the conusoid, which consisted of four layers of hybrid carbon-Kevlar® plain-weave fabric oriented at [+45°/-45°/-45°/+45°] with respect to the vertical direction. The second was a sinusoidal-shaped energy absorber, designated the sinusoid, which consisted of hybrid carbon-Kevlar® plain-weave fabric face sheets, two layers for each face sheet oriented at ±45° with respect to the vertical direction, and a closed-cell ELFOAM® P200 polyisocyanurate foam core. Finite element models were developed of the energy absorbers and simulations were performed using LS-DYNA®. In this paper, the development of a *MAT_213 model of a hybrid carbon-Kevlar® plain-weave fabric is presented. Next, comparisons with material characterization tests are presented. Then, test-analysis results are documented for each energy absorber as comparisons of time-history responses, as well as predicted and experimental structural deformations and progressive damage under impact loading using the *MAT_213 material model. Since a prior *MAT_58, or *MAT_LAMINATED_COMPOSITE_FABRIC material model was used in previous simulations of the energy absorbers, comparisons are made between *MAT_58 and *MAT_213 model predictions with test data. Finally, the paper includes a comprehensive list of “lessons learned,” in which a series of parametric studies are documented that were performed to investigate specific issues related to the material model. These “lessons learned” are included in hopes that they may help future *MAT_213 users.

crashworthiness, LS-DYNA transient dynamic simulat↗
NASA NTRS1987

Comparison of composite rotor blade models: A coupled-beam analysis and an MSC/NASTRAN finite-element model

A methodology was developed for the structural analysis of composite rotor blades. This coupled-beam analysis is relatively simple to use compared with alternative analysis techniques. The beam analysis was developed for thin-wall single-cell rotor structures and includes the effects of elastic coupling. This paper demonstrates the effectiveness of the new composite-beam analysis method through comparison of its results with those of an established baseline analysis technique. The baseline analysis is an MSC/NASTRAN finite-element model built up from anisotropic shell elements. Deformations are compared for three linear static load cases of centrifugal force at design rotor speed, applied torque, and lift for an ideal rotor in hover. A D-spar designed to twist under axial loading is the subject of the analysis. Results indicate the coupled-beam analysis is well within engineering accuracy.

Hodges, Robert V.↗