Search NASA⌕ Search

SEARCH · Search NASA

Results for “shell model”

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 253 records · Page 14

Numerical investigation of sound transmission through double wall cylinders with respect to active noise control

Market pressure for more fuel efficient air travel has led to increased use of turboprop and higher bypass turbofan engines. The low frequency components of propeller, jet and boundary layer noise are difficult to attenuate with conventional passive techniques. Weight and geometric restrictions for sound absorbing meterials limit the amount and type of treatment that may be applied. An active noise control (ANC) method is providing to be an attractive alternative. The approach taken in this paper uses a numerical finite/boundary element method (FEM/BEM) that may be easilty adapted to arbitrary geometries. A double walled cylinder is modeled using commercially available software. The outer shell is modeled as an aluminum cylinder, similar to that of aircraft skins. The inner shell is modeled as a composite material representative of a lightweight, stiff trim panel. Two different inner shell materials are used. The first is representative of current trim structure, the second a much stiffer composite. The primary source is generated by an exterior acoustic monopole. Control fields are generated using normal force inputs to the inner cylindrical shell. A linear least mean square (LMS) algorithm is used to determine amplitudes of control forces that minimize the interior acoustic field. Coupling of acoustic and structural modes and noise reductions are discussed for each of the inner shell materials.

Coats, T. J.↗

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↗

Planetary Geophysics and Tectonics

The broad objective of this work is to improve understanding of the internal structures and thermal and stress histories of the solid planets by combining results from analytical and computational modeling, and geophysical data analysis of gravity, topography and tectonic surface structures. During the past year we performed two quite independent studies in the attempt to explain the Mariner 10 magnetic observations of Mercury. In the first we revisited the possibility of crustal remanence by studying the conditions under which one could break symmetry inherent in Runcorn's model of a uniformly magnetized shell to produce a remanent signal with a dipolar form. In the second we applied a thin shell dynamo model to evaluate the range of intensity/structure for which such a planetary configuration can produce a dipole field consistent with Mariner 10 results. In the next full proposal cycle we will: (1) develop numerical and analytical and models of thin shell dynamos to address the possible nature of Mercury s present-day magnetic field and the demise of Mars magnetic field; (2) study the effect of degree-1 mantle convection on a core dynamo as relevant to the early magnetic field of Mars; (3) develop models of how the deep mantles of terrestrial planets are perturbed by large impacts and address the consequences for mantle evolution; (4) study the structure, compensation, state of stress, and viscous relaxation of lunar basins, and address implications for the Moon s state of stress and thermal history by modeling and gravity/topography analysis; and (5) use a three-dimensional viscous relaxation model for a planet with generalized vertical viscosity distribution to study the degree-two components of the Moon's topography and gravity fields to constrain the primordial stress state and spatial heterogeneity of the crust and mantle.

Zuber, Maria↗

Neotectonics of Asia: Thin-shell finite-element models with faults

As India pushed into and beneath the south margin of Asia in Cenozoic time, it added a great volume of crust, which may have been (1) emplaced locally beneath Tibet, (2) distributed as regional crustal thickening of Asia, (3) converted to mantle eclogite by high-pressure metamorphism, or (4) extruded eastward to increase the area of Asia. The amount of eastward extrusion is especially controversial: plane-stress computer models of finite strain in a continuum lithosphere show minimal escape, while laboratory and theoretical plane-strain models of finite strain in a faulted lithosphere show escape as the dominant mode. We suggest computing the present (or neo)tectonics by use of the known fault network and available data on fault activity, geodesy, and stress to select the best model. We apply a new thin-shell method which can represent a faulted lithosphere of realistic rheology on a sphere, and provided predictions of present velocities, fault slip rates, and stresses for various trial rheologies and boundary conditions. To minimize artificial boundaries, the models include all of Asia east of 40 deg E and span 100 deg on the globe. The primary unknowns are the friction coefficient of faults within Asia and the amounts of shear traction applied to Asia in the Himalayan and oceanic subduction zones at its margins. Data on Quaternary fault activity prove to be most useful in rating the models. Best results are obtained with a very low fault friction of 0.085. This major heterogeneity shows that unfaulted continum models cannot be expected to give accurate simulations of the orogeny. But, even with such weak faults, only a fraction of the internal deformation is expressed as fault slip; this means that rigid microplate models cannot represent the kinematics either. A universal feature of the better models is that eastern China and southeast Asia flow rapidly eastward with respect to Siberia. The rate of escape is very sensitive to the level of shear traction in the Pacific subduction zones, which is below 6 MPa. Because this flow occurs across a wide range of latitudes, the net eastward escape is greater than the rate of crustal addition in the Himalaya. The crustal budget is balanced by extension and thinning, primarily within the Tibetan plateau and the Baikal rift. The low level of deviation stresses in the best models suggests that topographic stress plays a major role in the orogeny; thus, we have to expect that different topography in the past may have been linked with fundamentally different modes of continental collision.

Kong, Xianghong↗

Finite Element Modeling of the Buckling Response of Sandwich Panels

A comparative study of different modeling approaches for predicting sandwich panel buckling response is described. The study considers sandwich panels with anisotropic face sheets and a very thick core. Results from conventional analytical solutions for sandwich panel overall buckling and face-sheet-wrinkling type modes are compared with solutions obtained using different finite element modeling approaches. Finite element solutions are obtained using layered shell element models, with and without transverse shear flexibility, layered shell/solid element models, with shell elements for the face sheets and solid elements for the core, and sandwich models using a recently developed specialty sandwich element. Convergence characteristics of the shell/solid and sandwich element modeling approaches with respect to in-plane and through-the-thickness discretization, are demonstrated. Results of the study indicate that the specialty sandwich element provides an accurate and effective modeling approach for predicting both overall and localized sandwich panel buckling response. Furthermore, results indicate that anisotropy of the face sheets, along with the ratio of principle elastic moduli, affect the buckling response and these effects may not be represented accurately by analytical solutions. Modeling recommendations are also provided.

Rose, Cheryl A.↗

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.↗

Gamma-rays from cosmic ray interactions in supernova shells

A model is presented for the transport and interaction of cosmic rays accelerated by a pulsar and confined inside an expanding supernova remnant. Assuming that protons are accelerated at the reverse shock in the confined pulsar wind and convected into the shell via the Rayleigh-Taylor instability, the diffusion and interaction of these protons in the expanding envelope is modeled. The resulting gamma-ray flux is lower than previous estimates due primarily to the inclusion of proton adiabatic losses in the expanding pulsar wind. Energy-dependent diffusion causes the higher energy gamma-ray light curves to decay faster than those at lower energy. The predicted flux from SN1987A, for proton luminosity less than 10 exp 40 erg/s, is below the present detector sensitivities at both GeV and TeV energies, although supernovae occurring within the Galaxy may be detectable.

Harding, A. K.↗

Application of the line-spring model to a cylindrical shell containing a circumferential or axial part-through crack

The line-spring model developed by Rice and Levy (1972) is used to obtain an approximate solution for a cylindrical shell containing a part-through surface crack. A Reissner type theory is used to account for the effects of the transverse shear deformations, and the stress intensity factor at the deepest penetration point of the crack is tabulated for bending and membrane loading by varying three-dimensionless length parameters of the problem formed from the shell radius, the shell thickness, the crack length, and the crack depth. The upper bounds of the stress intensity factors are provided, and qualitatively the line-spring model gives the expected results in comparison with elasticity solutions.

Delale, F.↗

Neptune's magnetic field - Calculation of field geometric invariants derived from the I8E1 GSFC/BRI model

L shell values along the Voyager 2 encounter trajectory and those associated with the N1 through N6 moons and N1R through N6R rings of Neptune are computed numerically on the basis of a simplified description of the Neptunian magnetic field derived from the Goddard Space Flight Center/Bartol Research Institute I8E1 model, which includes internal terms up to and including the octupole (but no external terms). Like Uranus, the large tilt between the dipole term and the rotation axis causes the moons and rings to sweep a very large range of L shells. Their orbital motion introduces additional periodicities, causing the maxima and minima in L space to vary systematically with time.

Acuna, M. H.↗

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↗

Finite Element Analysis of Geodesically Stiffened Cylindrical Composite Shells Using a Layerwise Theory

Layerwise finite element analyses of geodesically stiffened cylindrical shells are presented. The layerwise laminate theory of Reddy (LWTR) is developed and adapted to circular cylindrical shells. The Ritz variational method is used to develop an analytical approach for studying the buckling of simply supported geodesically stiffened shells with discrete stiffeners. This method utilizes a Lagrange multiplier technique to attach the stiffeners to the shell. The development of the layerwise shells couples a one-dimensional finite element through the thickness with a Navier solution that satisfies the boundary conditions. The buckling results from the Ritz discrete analytical method are compared with smeared buckling results and with NASA Testbed finite element results. The development of layerwise shell and beam finite elements is presented and these elements are used to perform the displacement field, stress, and first-ply failure analyses. The layerwise shell elements are used to model the shell skin and the layerwise beam elements are used to model the stiffeners. This arrangement allows the beam stiffeners to be assembled directly into the global stiffness matrix. A series of analytical studies are made to compare the response of geodesically stiffened shells as a function of loading, shell geometry, shell radii, shell laminate thickness, stiffener height, and geometric nonlinearity. Comparisons of the structural response of geodesically stiffened shells, axial and ring stiffened shells, and unstiffened shells are provided. In addition, interlaminar stress results near the stiffener intersection are presented. First-ply failure analyses for geodesically stiffened shells utilizing the Tsai-Wu failure criterion are presented for a few selected cases.

Gerhard, Craig Steven↗

The modeling of piezoceramic patch interactions with shells, plates and beams

General models describing the interactions between a pair of piezoceramic patches and elastic substructures consisting of a cylindrical shell, plate and beam are presented. In each case, the manner in which the patch loads enter both the strong and weak forms of the time-dependent structural equations of motion is described. Through force and moment balancing, these loads are then determined in terms of material properties of the patch and substructure (thickness, elastic properties, Poisson ratios), the geometry of the patch placement, and the voltages into the patches. In the case of the shell, the coupling between banding and inplane deformations, which is due to the curvature, is retained. These models are sufficiently general to allow for potentially different patch voltages which implies that they can be suitably employed when using piezoceramic patches for controlling system dynamics when both extensional and bending vibrations are present.

Banks, H. T.↗

Shear-flexible finite-element models of laminated composite plates and shells

Several finite-element models are applied to the linear static, stability, and vibration analysis of laminated composite plates and shells. The study is based on linear shallow-shell theory, with the effects of shear deformation, anisotropic material behavior, and bending-extensional coupling included. Both stiffness (displacement) and mixed finite-element models are considered. Discussion is focused on the effects of shear deformation and anisotropic material behavior on the accuracy and convergence of different finite-element models. Numerical studies are presented which show the effects of increasing the order of the approximating polynomials, adding internal degrees of freedom, and using derivatives of generalized displacements as nodal parameters.

Noor, A. K.↗

Nonlinear probabilistic finite element models of laminated composite shells

A probabilistic finite element analysis procedure for laminated composite shells has been developed. A total Lagrangian finite element formulation, employing a degenerated 3-D laminated composite shell with the full Green-Lagrange strains and first-order shear deformable kinematics, forms the modeling foundation. The first-order second-moment technique for probabilistic finite element analysis of random fields is employed and results are presented in the form of mean and variance of the structural response. The effects of material nonlinearity are included through the use of a rate-independent anisotropic plasticity formulation with the macroscopic point of view. Both ply-level and micromechanics-level random variables can be selected, the latter by means of the Aboudi micromechanics model. A number of sample problems are solved to verify the accuracy of the procedures developed and to quantify the variability of certain material type/structure combinations. Experimental data is compared in many cases, and the Monte Carlo simulation method is used to check the probabilistic results. In general, the procedure is quite effective in modeling the mean and variance response of the linear and nonlinear behavior of laminated composite shells.

Engelstad, S. P.↗

Characterization and Modeling of Cell Wall Imperfections in Aluminum Honeycomb Cores using X-ray CT Imaging

Honeycomb cell imperfections may significantly affect the failure of sandwich composites in core crushing and shear modes. Identification and quantification of the different types of imperfections and evaluation of their effects on the mechanical response of honeycomb cores in flatwise compression is presented. X-Ray computed tomography (CT) scans of honeycomb cores in co-cured composite sandwich panels were performed. Image analysis methods are applied to characterize and quantify core cell geometry and cell wall waviness. An approach is developed for extracting cell wall mid surface planes and waviness for constructing finite element models using shell elements. Analysis results from models with the as-manufactured imperfections are compared to results from models with idealized geometries to evaluate the effect of the as-manufactured imperfections on the initiation of cell wall buckling and subsequent core crushing.

aluminum honeycomb↗

Imperfection Insensitivity Analyses of Advanced Composite Tow-Steered Shells

Two advanced composite tow-steered shells, one with tow overlaps and another without overlaps, were previously designed, fabricated and tested in end compression, both without cutouts, and with small and large cutouts. In each case, good agreement was observed between experimental buckling loads and supporting linear bifurcation buckling analyses. However, previous buckling tests and analyses have shown historically poor correlation, perhaps due to the presence of geometric imperfections that serve as failure initiators. For the tow-steered shells, their circumferential variation in axial stiffness may have suppressed this sensitivity to imperfections, leading to the agreement noted between tests and analyses. To investigate this further, a numerical investigation was performed in this study using geometric imperfections measured from both shells. Finite element models of both shells were analyzed first without, and then, with measured imperfections that were then, superposed in different orientations around the shell longitudinal axis. Small variations in both the axial prebuckling stiffness and global buckling load were observed for the range of imperfections studied here, which suggests that the tow steering, and resulting circumferentially varying axial stiffness, may result in the test-analysis correlation observed for these shells.

Wu, K. Chauncey↗