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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Three-dimensional reconstruction of implosion stagnation in laser direct drive on OMEGA

Multidimensional effects on hot-spot formation must be considered to better understand the current limits on the performance of direct-drive inertial confinement fusion experiments on OMEGA with cryogenically layered solid deuterium–tritium targets. A comprehensive reconstruction effort has been established at the Laboratory for Laser Energetics to infer hot-spot and shell conditions at stagnation from a large collection of x-ray, neutron, and particle detectors along multiple lines of sight. Several time-gated and time-integrated x-ray imagers are being used to record the shape of the hot-spot plasma. A 3D hot-spot x-ray emission tomography technique has been developed to infer low-mode drive asymmetries from the hot-spot shape. A suite of neutron diagnostics is used to provide measurements of hot-spot flow velocity, ion temperature, and areal density. Here, the information obtained from the x-ray and neutron detectors will be combined into a coherent model of the shape of the hot spot and shell assembly.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1992

Formation of shells in major mergers

Numerical simulations are used to study the fate of disks in mergers between equal-mass galaxies. Contrary to popular belief, mergers between equal-mass galaxies can form shells, loops, and ripples. Material that was originally in the outer disk of the premerger spirals falls into the remnant late in the merger event, long after the inner region of the remnant has relaxed. Thus, its evolution is similar to that of an accreted dwarf satellite that forms shells through 'phase wrapping'. However, the mechanism described in this letter avoids a number of difficulties with the accretion model; specifically, it can explain the observed correlation between properties of the host galaxies and shell alignments and luminosities. In view of this, we argue that many shell systems may have originated through 'major' mergers of comparable-mass galaxies rather than via 'minor' mergers or accretions. In particular, the shell elliptical NGC 3923 appears to be a good candidate for shell formation by a major merger.

Hernquist, Lars↗
NASA NTRS2015

In-Flight Aeroelastic Stability of the Thermal Protection System on the NASA HIAD, Part II: Nonlinear Theory and Extended Aerodynamics

Conical shell theory and a supersonic potential flow aerodynamic theory are used to study the nonlinear pressure buckling and aeroelastic limit cycle behavior of the thermal protection system for NASA's Hypersonic Inflatable Aerodynamic Decelerator. The structural model of the thermal protection system consists of an orthotropic conical shell of the Donnell type, resting on several circumferential elastic supports. Classical Piston Theory is used initially for the aerodynamic pressure, but was found to be insufficient at low supersonic Mach numbers. Transform methods are applied to the convected wave equation for potential flow, and a time-dependent aerodynamic pressure correction factor is obtained. The Lagrangian of the shell system is formulated in terms of the generalized coordinates for all displacements and the Rayleigh-Ritz method is used to derive the governing differential-algebraic equations of motion. Aeroelastic limit cycle oscillations and buckling deformations are calculated in the time domain using a Runge-Kutta method in MATLAB. Three conical shell geometries were considered in the present analysis: a 3-meter diameter 70 deg. cone, a 3.7-meter 70 deg. cone, and a 6-meter diameter 70 deg. cone. The 6-meter configuration was loaded statically and the results were compared with an experimental load test of a 6-meter HIAD. Though agreement between theoretical and experimental strains was poor, the circumferential wrinkling phenomena observed during the experiments was captured by the theory and axial deformations were qualitatively similar in shape. With Piston Theory aerodynamics, the nonlinear flutter dynamic pressures of the 3-meter configuration were in agreement with the values calculated using linear theory, and the limit cycle amplitudes were generally on the order of the shell thickness. The effect of axial tension was studied for this configuration, and increasing tension was found to decrease the limit cycle amplitudes when the circumferential elastic supports were neglected, but resulted in more complex behavior when the supports were included. The nominal flutter dynamic pressure of the 3.7-meter configuration was significantly lower than that of the 3-meter, and it was found that two sets of natural modes coalesce to flutter modes near the same dynamic pressure. This resulted in a significant drop in the limit cycle frequencies at higher dynamic pressures, where the flutter mode with the lower frequency becomes more critical. Pre-buckling pressure loads and the aerodynamic pressure correction factor were studied for all geometries, and these effects resulted in significantly lower flutter boundaries compared with Piston Theory alone. The maximum dynamic pressure predicted by aerodynamic simulations of a proposed 3.7-meter HIAD vehicle was still lower than any of the calculated flutter dynamic pressures, suggesting that aeroelastic effects for this vehicle are of little concern.

Goldman, Benjamin D.↗
NASA NTRS2011

Saturn's Magnetosphere Interaction with Titan for T9 Encounter: 3D Hybrid Modeling and Comparison with CAPS Observations

Global dynamics of ionized and neutral gases in the environment of Titan plays an important role in the interaction of Saturn s magnetosphere with Titan. Several hybrid simulations of this problem have already been done (Brecht et al., 2000; Kallio et al., 2004; Modolo et al., 2007a; Simon et al., 2007a, 2007b; Modolo and Chanteur, 2008). Observational data from CAPS for the T9 encounter (Sittler et al., 2009) indicates an absence of O(+) heavy ions in the upstream that change the models of interaction which were discussed in current publications (Kallio et al., 2004; Modolo et al., 2007a; Simon et al., 2007a, 2007b; Ma et al., 2007; Szego et al., 2007). Further analysis of the CAPS data shows very low density or even an absence of H(+) ions in upstream. In this paper we discuss two models of the interaction of Saturn s magnetosphere with Titan: (A) high density of H(+) ions in the upstream flow (0.1/cu cm), and (B) low density of H(+) ions in the upstream flow (0.02/cu cm). The hybrid model employs a fluid description for electrons and neutrals, whereas a particle approach is used for ions. We also take into account charge-exchange and photoionization processes and solve self-consistently for electric and magnetic fields. The model atmosphere includes exospheric H(+), H(2+), N(2+)and CH(4+) pickup ion production as well as an immobile background ionosphere and a shell distribution for active ionospheric ions (M(sub i)=28 amu). The hybrid model allows us to account for the realistic anisotropic ion velocity distribution that cannot be done in fluid simulations with isotropic temperatures. Our simulation shows an asymmetry of the ion density distribution and the magnetic field, including the formation of Alfven wing-like structures. The results of the ion dynamics in Titan s environment are compared with Cassini T9 encounter data (CAPS).

Lipatov, A. S.↗
NASA NTRS1999

A Mixed Multi-Field Finite Element Formulation for Thermopiezoelectric Composite Shells

Analytical formulations are presented which account for the coupled mechanical, electrical, and thermal response of piezoelectric composite shell structures. A new mixed multi-field laminate theory is developed which combines "single layer" assumptions for the displacements along with layerwise fields for the electric potential and temperature. This laminate theory is formulated using curvilinear coordinates and is based on the principles of linear thermopiezoelectricity. The mechanics have the inherent capability to explicitly model both the active and sensory responses of piezoelectric composite shells in thermal environment. Finite element equations are derived and implemented for an eight-noded shell element. Numerical studies are conducted to investigate both the sensory and active responses of piezoelectric composite shell structures subjected to thermal loads. Results for a cantilevered plate with an attached piezoelectric layer are com- pared with corresponding results from a commercial finite element code and a previously developed program. Additional studies are conducted on a cylindrical shell with an attached piezoelectric layer to demonstrate capabilities to achieve thermal shape control on curved piezoelectric structures.

Lee, Ho-Jun↗
NASA NTRS1998

Analysis and Design of Variable Stiffness Composite Cylinders

An investigation of the possible performance improvements of thin circular cylindrical shells through the use of the variable stiffness concept is presented. The variable stiffness concept implies that the stiffness parameters change spatially throughout the structure. This situation is achieved mainly through the use of curvilinear fibers within a fiber-reinforced composite laminate, though the possibility of thickness variations and discrete stiffening elements is also allowed. These three mechanisms are incorporated into the constitutive laws for thin shells through the use of Classical Lamination Theory. The existence of stiffness variation within the structure warrants a formulation of the static equilibrium equations from the most basic principles. The governing equations include sufficient detail to correctly model several types of nonlinearity, including the formation of a nonlinear shell boundary layer as well as the Brazier effect due to nonlinear bending of long cylinders. Stress analysis and initial buckling estimates are formulated for a general variable stiffness cylinder. Results and comparisons for several simplifications of these highly complex governing equations are presented so that the ensuing numerical solutions are considered reliable and efficient enough for in-depth optimization studies. Four distinct cases of loading and stiffness variation are chosen to investigate possible areas of improvement that the variable stiffness concept may offer over traditional constant stiffness and/or stiffened structures. The initial investigation deals with the simplest solution for cylindrical shells in which all quantities are constant around the circumference of the cylinder. This axisymmetric case includes a stiffness variation exclusively in the axial direction, and the only pertinent loading scenarios include constant loads of axial compression, pressure, and torsion. The results for these cases indicate that little improvement over traditional laminates exists through the use of curvilinear fibers, mainly due to the presence of a weak link area within the stiffness variation that limits the ultimate load that the structure can withstand. Rigorous optimization studies reveal that even though slight increases in the critical loads can be produced for designs with an arbitrary variation of the fiber orientation angle, the improvements are not significant when compared to traditional design techniques that utilize ring stiffeners and frames. The second problem that is studied involves arbitrary loading of a cylinder with a stiffness variation that changes only in the circumferential direction. The end effects of the cylinder are ignored, so that the problem takes the form of an analysis of a cross-section for a short cylinder segment. Various load cases including axial compression, pressure, torsion, bending, and transverse shear forces are investigated. It is found that the most significant improvements in load-carrying capability exist for cases which involve loads that also vary around the circumference of the shell, namely bending and shear forces. The stiffness variation of the optimal designs contribute to the increased performance in two ways: lowering the stresses in the critical areas through redistribution of the stresses; and providing a relatively stiff region that alters the buckling behavior of the structure. These results lead to an in-depth optimization study involving weight optimization of a fuselage structure subjected to typical design constraints. Comparisons of the curvilinear fiber format to traditional stiffened structures constructed of isotropic and composite materials are included. It is found that standard variable stiffness designs are quite comparable in terms of weight and load-carrying capability yet offer the added advantage of tailorability of distinct regions of the structure that experience drastically different loading conditions. The last two problems presented in this work involve the nonlinear phenomenon of long tubes under bending. Though this scenario is not as applicable to fuselage structures as the previous problems, the mechanisms that produce the nonlinear effect are ideally suited to be controlled by the variable stiffness concept. This is due to the fact that the dominating influence for long cylinders under bending is the ovalization of the cross-section, which is governed mainly by the stiffness parameters of the cylindrical shell. Possible improvement of the critical buckling moments for these structures is investigated using either a circumferential or axial stiffness variation. For the circumferential case involving infinite length cylinders, it is found that slight improvements can be observed by designing structures that resist the cross-sectional deformation yet do not detract from the buckling resistance at the critical location. The results also indicate that buckling behavior is extremely dependent on cylinder length. This effect is most easily seen in the solution of finite length cylinders under bending that contain an axial stiffness variation. For these structures, the only mechanism that exhibits improved response are those that effectively shorten the length of the cylinder, thus reducing the cross-sectional deformation due to the forced restraint at the ends. It was found that the use of curvilinear fibers was not able to achieve this effect in sufficient degree to resist the deformation, but that ring stiffeners produced the desired response admirably. Thus, it is shown that the variable stiffness concept is most effective at improving the bending response of long cylinders through the use of a circumferential stiffness variation.

Tatting, Brian F.↗
NASA NTRS2010

Finite Element Model Calibration Approach for Ares I-X

Ares I-X is a pathfinder vehicle concept under development by NASA to demonstrate a new class of launch vehicles. Although this vehicle is essentially a shell of what the Ares I vehicle will be, efforts are underway to model and calibrate the analytical models before its maiden flight. Work reported in this document will summarize the model calibration approach used including uncertainty quantification of vehicle responses and the use of nonconventional boundary conditions during component testing. Since finite element modeling is the primary modeling tool, the calibration process uses these models, often developed by different groups, to assess model deficiencies and to update parameters to reconcile test with predictions. Data for two major component tests and the flight vehicle are presented along with the calibration results. For calibration, sensitivity analysis is conducted using Analysis of Variance (ANOVA). To reduce the computational burden associated with ANOVA calculations, response surface models are used in lieu of computationally intensive finite element solutions. From the sensitivity studies, parameter importance is assessed as a function of frequency. In addition, the work presents an approach to evaluate the probability that a parameter set exists to reconcile test with analysis. Comparisons of pre-test predictions of frequency response uncertainty bounds with measured data, results from the variance-based sensitivity analysis, and results from component test models with calibrated boundary stiffness models are all presented.

Horta, Lucas G.↗
NASA NTRS2010

Finite Element Model Calibration Approach for Area I-X

Ares I-X is a pathfinder vehicle concept under development by NASA to demonstrate a new class of launch vehicles. Although this vehicle is essentially a shell of what the Ares I vehicle will be, efforts are underway to model and calibrate the analytical models before its maiden flight. Work reported in this document will summarize the model calibration approach used including uncertainty quantification of vehicle responses and the use of non-conventional boundary conditions during component testing. Since finite element modeling is the primary modeling tool, the calibration process uses these models, often developed by different groups, to assess model deficiencies and to update parameters to reconcile test with predictions. Data for two major component tests and the flight vehicle are presented along with the calibration results. For calibration, sensitivity analysis is conducted using Analysis of Variance (ANOVA). To reduce the computational burden associated with ANOVA calculations, response surface models are used in lieu of computationally intensive finite element solutions. From the sensitivity studies, parameter importance is assessed as a function of frequency. In addition, the work presents an approach to evaluate the probability that a parameter set exists to reconcile test with analysis. Comparisons of pretest predictions of frequency response uncertainty bounds with measured data, results from the variance-based sensitivity analysis, and results from component test models with calibrated boundary stiffness models are all presented.

Horta, Lucas G.↗
NASA NTRS1993

Circumstellar shells resolved in the IRAS survey data. I - Data processing procedure, results, and confidence tests

We have examined the IRAS 60 and 100 micron survey data covering 512 evolved stars and young planetary nebulae for evidence of spatially resolved structure. A simple model, consisting of a central unresolved source surrounded by a resolved isothermal shell, was fitted to the data for each star. Seventy-six stars were found to be resolved in the 60 micron data. Tests have been performed to verify that the extended structure seen is not an artifact of the data-processing algorithm.

Young, K.↗
NASA NTRS1972

Plasticity and the crack opening displacement in shells.

The plastic deformations and the crack opening displacement, delta, in cylindrical shells with an axial or circumferential crack and spherical shells with a meridional crack are considered. It is assumed that outside the perturbation zone of the crack the shell is subjected to uniform membrane loads perpendicular to the crack. The plastic strip model is used to calculate the plastic zone size. The crack opening displacement is calculated as the crack surface displacement at the crack tips by using the conventional superposition technique. In cylindrical shells with an axial crack the crack surface displacement perpendicular to the shell surface (i.e., the bulging) is also evaluated. The results are applied to a set of existing experimental data on the fracture of cylindrical shells. The tentative conclusion is that in dealing with the fracture of thin-walled structures, among various fracture models delta = constant appears to be the most satisfactory criterion.

Erdogan, F.↗
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 NTRS1986

Convection electric field effects on outer radiation belt electron precipitation

A model is presented for the possible diurnal modulation of outer radiation belt electron precipitation by considering the effect of the convection electric field on geomagnetically trapped electrons. The modulation flux is the flux due to electrons in the drift loss cone, i.e., those which drift into the bounce loss cone. The electron flux in the drift loss cone is related to the time allowable for diffusion from the stably trapped population to the drift loss cone for precipitation at a specific geographic location. This time, which is termed the maximum L-shell lifetime, is obtained by computing electron trajectories, using a realistic magnetic field model and a simple model for the electric field. The maximum L-shell lifetimes are taken to be the times between successive entries into the bounce loss cone. Conservation of the first two adiabatic invariants, as electrons are slowly energized by the convection electric field, leads to variations in pitch angle, maximum L-shell lifetimes, and, consequently, to changes in the electron flux in the drift loss cone. These results are compared with observations of precipitating electrons made with sounding rocket payloads.

Gelpi, C.↗
NASA NTRS1994

Thin shell formation in radiative shocks. 1: Supernova remnants in low-density media

This paper explores the onset of thin-shell formation in interstellar shocks associated with supernova explosions. We outline a simple but useful scheme that indicates the time at which thin shell formation begins for supernova remnants (SNRs) evolving in a range of interstellar environments, extending the previous analytical models to arbitrary power-law density media. The result depends on the gas cooling properties and the shock velocity and radius. This is then applied to the specific case of SNRs in low-density media. The procedure for defining the time for the onset of shell formation, t(sub sf), equates the value of the adiabat, kappa = p/rho(exp gamma), to zero using the known time dependence of the shock radius and velocity. For the case of a power-law density ambient medium of the form rho(r) = Br(exp -omega), it is found that shell formation can be prevented when the ambient density drops faster than a critical rate. For a cooling function of the form Lambda = Lambda(sub 0) tau(exp beta), with beta = -0.5 (appropriate for line cooling), shell formation never occurs for omega greater than or equal to 9/5. The shell formation time is then computed for spherical shocks in a power-law density medium. For omega = 0, the onset of shell formation is found to be at t(sub sf) approx. equal to 2.87 x 10(exp 4) E51(exp 3/14) n(sub 0 exp -4.7) yr, which agrees well with previous estimates derived by other means. We compare the analytical shell formation time with the results of detailed numerical models for omega = 0 and three different ambient densities and find good agreement. The extension of the criterion for the onset of thin shell formation using the ratio of cooling to swept-up column density is also described. This method provides a useful approximation for cases when the exact solution is not known.

Franco, Jose↗
NASA NTRS2003

Multi-Temperature Emission and Abundances in the Hot Gaseous Halo

We present the results of XMM-Newton observations of NGC 507, a dominant elliptical galaxy in a small group of galaxies. After carefully considering various systematic effects on abundance measurements, we report 'super-solar' metal abundances (both Fe and a-elements) present in the hot ISM: ZFe = 2-3 times solar with an observational limit of as high as 4 times solar inside the D25 ellipse of NGC 507. This is the highest ZFe reported so far, and fully consistent with those expected by the stellar evolution models where heavy elements are enriched by both type II and Ia supernovae ejecta. No unusual constraint either on the SNe rate or IMF is required. Among various factors affecting the accurate abundance measurement, we find that selecting a proper emission model is most important. As opposed to the X-ray spectral data with limited s/n and poor spatial/spectral resolution obtained in the previous missions, the spatially resolved XMM spectra provide enough statistics to untie the model-Z degeneracy and statistically require at least 3 emission components in each concentric shell (2 thermal components representing a finite range of kT in the hot ISM + 1 hard LMXB component). We show that a simpler model (such as a two-component model) produce a much lower best-fit ZFe. The abundances of a-elements (most accurately determined by Si) is also found to be super-solar and its radio to Fe is close to the solar ratio, suggesting a considerably contribution of heavy elements from Type Ia SNe. We estimate approx. 70% of MFe in the hot ISM originate from Type Ia.

Kim, Dong-Woo↗
NASA NTRS1976

On the elastic stability of shells

A synopsis of a series of investigations into the instability of axially compressed cylindrical shells is given. The objective of the research, which was made with models, was to devise a technique of nondestructive evaluation. The results show that, with models at any rate, success was achieved. Probing methods which can be used to determine the locations of weakness and the pertinent instability load levels were devised. The research on large scale shells was undertaken to determine the critical loads under as uniform a circumferential distribution of axial compressive force as possible. It is clear from the results presented that this objective was met.

Horton, W. H.↗
NASA NTRS2004

Topographic variations in chaos on Europa: Implications for diapiric formation

Disrupted terrain, or chaos, on Europa, might have formed through melting of a floating ice shell from a subsurface ocean [Cam et al., 1998; Greenberg et al., 19991, or breakup by diapirs rising from the warm lower portion of the ice shell [Head and Pappalardo, 1999; Collins et al., 20001. Each model makes specific and testable predictions for topographic expression within chaos and relative to surrounding terrains on local and regional scales. High-resolution stereo-controlled photoclinometric topography indicates that chaos topography, including the archetypal Conamara Chaos region, is uneven and commonly higher than surrounding plains by up to 250 m. Elevated and undulating topography is more consistent with diapiric uplift of deep material in a relatively thick ice shell, rather than melt-through and refreezing of regionally or globally thin ice by a subsurface ocean. Vertical and horizontal scales of topographic doming in Conamara Chaos are consistent with a total ice shell thickness >15 km. Contact between Europa's ocean and surface may most likely be indirectly via diapirism or convection.

Schenk, Paul M.↗
NASA NTRS2001

The Large-Scale Plasmaspheric Density Trough Associated With the 24 May 2000 Geomagnetic Storm: IMAGE EUV Observations and Global Core Plasma Modeling

The IMAGE EUV imager observed a plasmaspheric density, trough in association with a geomagnetically active period on 24 May 2000. In EUV, this density, trough appears as an Archimedes spiral extending from Earth's shadow to approximately 1900 MLT. We present an analysis of this density trough using simulated EUV images, Observational EUV images are subjected to edge analysis to establish the plasmapause L-shell and the location of the density trough in terms of L-shell, MLT extent, and radial width. The plasmaspheric density distribution is modeled using both static and dynamic models for the plasmasphere. The background plasmasphere is then numerically simulated using the 4-parameter plasmaspheric density model contained within the Global Core Plasma Model (GCPM) and the Dynamic Global Core Plasma Model (DGCPM). Simulated EUV images of the model plasmasphere are produced once an artificial density, depletion, matching the observed MLT extent and width, has been removed. Once the azimuthal extent and width of the trough have been simulated, the depth of the artificial density depletion is iteratively adjusted to produce simulated EUV images that approximate observation. The results of this analysis and discussion of possible origins for this density trough will be presented.

Adrian, M. L.↗
NASA NTRS1989

Assessment of computational models for multilayered composite cylinders

A study is made of the effects of variation in the lamination and geometric parameters of multilayered composite cylinders on the accuracy of the static and vibrational responses predicted by eight modeling approaches, based on two-dimensional shear-deformation shell theories. The standard of comparison is taken to be the exact three-dimensional elasticity solutions, and the quantities compared include both the gross response characteristics and detailed, through-the-thickness, distributions of displacements, stresses and strain energy densities. Based on the numerical studies conducted, a predictor-corrector approach, used in conjunction with the first-order shear-deformation theory (with five displacement parameters in the predictor phase), appears to be the most effective among the eight modeling approaches considered. For multilayered orthotropic cylinders the response quantities obtained by the predictor-corrector approach are shown to be in close agreement with the exact three-dimensional elasticity solutions for a wide range of lamination and geometric parameters.

Noor, A. K.↗