Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct coupling analysis”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Static aeroelastic analysis of wings using Euler/Navier-Stokes equations coupled with improved wing-box finite element structures

Strong interactions between flow about an aircraft wing and the wing structure can result in aeroelastic phenomena which significantly impact aircraft performance. Time-accurate methods for solving the unsteady Navier-Stokes equations have matured to the point where reliable results can be obtained with reasonable computational costs for complex non-linear flows with shock waves, vortices and separations. The ability to combine such a flow solver with a general finite element structural model is key to an aeroelastic analysis in these flows. Earlier work involved time-accurate integration of modal structural models based on plate elements. A finite element model was developed to handle three-dimensional wing boxes, and incorporated into the flow solver without the need for modal analysis. Static condensation is performed on the structural model to reduce the structural degrees of freedom for the aeroelastic analysis. Direct incorporation of the finite element wing-box structural model with the flow solver requires finding adequate methods for transferring aerodynamic pressures to the structural grid and returning deflections to the aerodynamic grid. Several schemes were explored for handling the grid-to-grid transfer of information. The complex, built-up nature of the wing-box complicated this transfer. Aeroelastic calculations for a sample wing in transonic flow comparing various simple transfer schemes are presented and discussed.

Guruswamy, Guru P.↗
NASA NTRS2006

Skin Temperature Analysis and Bias Correction in a Coupled Land-Atmosphere Data Assimilation System

In an initial investigation, remotely sensed surface temperature is assimilated into a coupled atmosphere/land global data assimilation system, with explicit accounting for biases in the model state. In this scheme, an incremental bias correction term is introduced in the model's surface energy budget. In its simplest form, the algorithm estimates and corrects a constant time mean bias for each gridpoint; additional benefits are attained with a refined version of the algorithm which allows for a correction of the mean diurnal cycle. The method is validated against the assimilated observations, as well as independent near-surface air temperature observations. In many regions, not accounting for the diurnal cycle of bias caused degradation of the diurnal amplitude of background model air temperature. Energy fluxes collected through the Coordinated Enhanced Observing Period (CEOP) are used to more closely inspect the surface energy budget. In general, sensible heat flux is improved with the surface temperature assimilation, and two stations show a reduction of bias by as much as 30 Wm(sup -2) Rondonia station in Amazonia, the Bowen ratio changes direction in an improvement related to the temperature assimilation. However, at many stations the monthly latent heat flux bias is slightly increased. These results show the impact of univariate assimilation of surface temperature observations on the surface energy budget, and suggest the need for multivariate land data assimilation. The results also show the need for independent validation data, especially flux stations in varied climate regimes.

Bosilovich, Michael G.↗
NASA NTRS2020

On Boundary-Value Problems for RANS Equations and Two-Equation Turbulence Models

Currently, in engineering computations for high Reynolds number turbulent flows, turbulence modeling continues to be the most frequently used approach to represent the effects of turbulence. Such models generally rely on solving either one or two transport equations along with the Reynolds-Averaged Navier–Stokes (RANS) equations. The solution of the boundary-value problem of any system of partial differential equations requires the complete delineation of the equations and the boundary conditions, including any special restrictions and conditions. In the literature, such a description is often incomplete, neglecting important details related to the boundary conditions and possible restrictive conditions, such as how to ensure satisfying prescribed values of the dependent variables of the transport equations in the far field of a finite domain. In this article, we discuss the possible influence of boundary values, as well as near-field and far-field behavior, on the solution of the RANS equations coupled with transport equations for turbulence modeling. In so doing, we defne the concept of a welldefined boundary-value problem. Additionally, a three-dimensional, rather than a simpler one-dimensional analysis is performed to analyze the near-wall and far-field behavior of the turbulence model variables. This allows an assessment of the decay rate of these variables required to realize the boundary conditions in the far field. This paper also addresses the impact of various transformations of two-equation models (e.g., the model of Wilcox) to remove the singular behavior of the dissipation rate (ω) at the surface boundary. Finally, the issue of well-posedness regarding the governing equations is considered. A compelling argument (although not a proof) for ill-posedness is made for both direct and inverse problems.

Navier–Stokes equations↗
NASA NTRS2006

Pointing Reference Scheme for Free-Space Optical Communications Systems

A scheme is proposed for referencing the propagation direction of the transmit laser signal in pointing a free-space optical communications terminal. This recently developed scheme enables the use of low-cost, commercial silicon-based sensors for tracking the direction of the transmit laser, regardless of the transmit wavelength. Compared with previous methods, the scheme offers some advantages of less mechanical and optical complexity and avoids expensive and exotic sensor technologies. In free-space optical communications, the transmit beam must be accurately pointed toward the receiver in order to maintain the communication link. The current approaches to achieve this function call for part of the transmit beam to be split off and projected onto an optical sensor used to infer the pointed direction. This requires that the optical sensor be sensitive to the wavelength of the transmit laser. If a different transmit wavelength is desired, for example to obtain a source capable of higher data rates, this can become quite impractical because of the unavailability or inefficiency of sensors at these wavelengths. The innovation proposed here decouples this requirement by allowing any transmit wavelength to be used with any sensor. We have applied this idea to a particular system that transmits at the standard telecommunication wavelength of 1,550 nm and uses a silicon-based sensor, sensitive from 0.5 to 1.0 micrometers, to determine the pointing direction. The scheme shown in the figure involves integrating a low-power 980-nm reference or boresight laser beam coupled to the 1,550-nm transmit beam via a wavelength-division-multiplexed fiber coupler. Both of these signals propagate through the optical fiber where they achieve an extremely high level of co-alignment before they are launched into the telescope. The telescope uses a dichroic beam splitter to reflect the 980- nm beam onto the silicon image sensor (a quad detector, charge-coupled device, or active-pixel-sensor array) while the 1,550- nm signal beam is transmitted through the optical assembly toward the remotely located receiver. Since the 980-nm reference signal originates from the same single-mode fiber-coupled source as the transmit signal, its position on the sensor is used to accurately determine the propagation direction of the transmit signal. The optics are considerably simpler in the proposed scheme due to the use of a single aperture for transmitting and receiving. Moreover, the issue of mechanical misalignment does not arise because the reference signal and transmitted laser beams are inherently co-aligned. The beam quality of the 980-nm reference signal used for tracking is required to be circularly symmetric and stable at the tracking-plane sensor array in order to minimize error in the centroiding algorithm of the pointing system. However, since the transmit signal is delivered through a fiber that supports a single mode at 1,550 nm, propagation of higher order 980-nm modes is possible. Preliminary analysis shows that the overall mode profile is dominated by the fundamental mode, giving a near symmetric profile. The instability of the mode was also measured and found to be negligible in comparison to the other error contributions in the centroid position on the sensor array.

Wright, Malcolm↗
NASA NTRS1987

Effect of dropped plies on the strength of graphite-epoxy laminates

The reduction in the compressive and tensile strengths of graphite-epoxy laminates with thickness discontinuities due to dropped plies was studied by experiment and analysis. The specimens were fabricated with all the dropped plies lumped together in the center of a sixteen-ply quasi-isotropic layup, such that one surface was flat and the slope of the opposite surface changed abruptly at the dropped ply location to accommodate the thickness change. Even though the thick and thin sections are symmetrically laminated, there exists bending-extension coupling due to the geometric eccentricity of the middle planes of the thick and thin sections. Experiments were conducted on fifty-four specimens that differed in the configuration of the dropped plies only. The strength of a laminate with dropped plies is less than the strength of its thin section, and the compressive strength of a laminate with dropped plies is less than its tensile strength. The reduction in strength is directly related to the axial stiffness change between the thick and thin sections. To examine the mechanism of failure, the three-dimensional state of stress in the dropped ply region was evaluated by the finite element method. A tensile interlaminar criterion predicted the correct location of failure, but underestimated the failure load.

Curry, James M.↗
NASA NTRS1995

Comparison of Two Multidisciplinary Optimization Strategies for Launch-Vehicle Design

The investigation focuses on development of a rapid multidisciplinary analysis and optimization capability for launch-vehicle design. Two multidisciplinary optimization strategies in which the analyses are integrated in different manners are implemented and evaluated for solution of a single-stage-to-orbit launch-vehicle design problem. Weights and sizing, propulsion, and trajectory issues are directly addressed in each optimization process. Additionally, the need to maintain a consistent vehicle model across the disciplines is discussed. Both solution strategies were shown to obtain similar solutions from two different starting points. These solutions suggests that a dual-fuel, single-stage-to-orbit vehicle with a dry weight of approximately 1.927 x 10(exp 5)lb, gross liftoff weight of 2.165 x 10(exp 6)lb, and length of 181 ft is attainable. A comparison of the two approaches demonstrates that treatment or disciplinary coupling has a direct effect on optimization convergence and the required computational effort. In comparison with the first solution strategy, which is of the general form typically used within the launch vehicle design community at present, the second optimization approach is shown to he 3-4 times more computationally efficient.

Braun, R. D.↗
NASA NTRS2010

Simultaneous Spectral Temporal Adaptive Raman Spectrometer - SSTARS

Raman spectroscopy is a prime candidate for the next generation of planetary instruments, as it addresses the primary goal of mineralogical analysis, which is structure and composition. However, large fluorescence return from many mineral samples under visible light excitation can render Raman spectra unattainable. Using the described approach, Raman and fluorescence, which occur on different time scales, can be simultaneously obtained from mineral samples using a compact instrument in a planetary environment. This new approach is taken based on the use of time-resolved spectroscopy for removing the fluorescence background from Raman spectra in the laboratory. In the SSTARS instrument, a visible excitation source (a green, pulsed laser) is used to generate Raman and fluorescence signals in a mineral sample. A spectral notch filter eliminates the directly reflected beam. A grating then disperses the signal spectrally, and a streak camera provides temporal resolution. The output of the streak camera is imaged on the CCD (charge-coupled device), and the data are read out electronically. By adjusting the sweep speed of the streak camera, anywhere from picoseconds to milliseconds, it is possible to resolve Raman spectra from numerous fluorescence spectra in the same sample. The key features of SSTARS include a compact streak tube capable of picosecond time resolution for collection of simultaneous spectral and temporal information, adaptive streak tube electronics that can rapidly change from one sweep rate to another over ranges of picoseconds to milliseconds, enabling collection of both Raman and fluorescence signatures versus time and wavelength, and Synchroscan integration that allows for a compact, low-power laser without compromising ultimate sensitivity.

Blacksberg, Jordana↗
NASA NTRSDate not supplied

Lunar South Pole Regenerative Fuel Cell System Efficiency Analysis

It is expected that future NASA mission energy storage needs can be adequately met by regenerative fuel cells (RFCs). Aerospace RFCs package proton exchange membrane (PEM) electrochemical conversion devices based on hydrogen, oxygen, and water. To optimize system designs and direct further development, the most useful system metric is round-trip efficiency (RTE). Improving RTE generally increases specific energy and reduces system mass. The following analysis models the impact of several factors on RTE for RFCs at the lunar south pole scaling from 0.1 kW to 50 kW. The analysis demonstrates that higher electrolyzer operational pressure most negatively impacts both RFC RTE and specific energy. Optimal electrolyzer membrane thickness is evaluated as a function of operational current density and pressure. Reactant storage heating and power management losses are significant for RFCs of any scale but are reduced in impact at the south pole compared to colder equator applications. Solenoid valve power remains a noticeable parasitic load for smaller RFCs. Larger electrolyzers enable higher peak theoretical efficiencies when recharging at higher rates, if the RFC can be coupled with a corresponding higher power source. RFC design is a compromise of many interrelated factors, but this analysis provides a starting point for difficult decisions that must be made when designing a lunar south pole RFC.

regenerative fuel cell↗
NASA NTRS1974

An analysis of the radiation from apertures in curved surfaces by the geometrical theory of diffraction

In this paper the geometrical theory of diffraction is extended to treat the radiation from apertures of slots in convex perfectly conducting surfaces. It is assumed that the tangential electric field in the aperture is known so that an equivalent infinitesimal source can be defined at each point in the aperture. Surface rays emanate from this source which is a caustic of the ray system. A launching coefficient is introduced to describe the excitation of the surface ray modes. If the field radiated from the surface is desired, the ordinary diffraction coefficients are used to determine the field of the rays shed tangentially from the surface rays. The field of the surface ray modes is not the field on the surface; hence if the mutual coupling between slots is of interest, a second coefficient related to the launching coefficient must be employed. In the region adjacent to the shadow boundary, the component of the field directly radiated from the source is represented by Fock-type functions. In the illuminated region the incident radiation from the source (this does not include the diffracted field components) is treated by geometrical optics. This extension of the geometrical theory of diffraction is applied to calculate the radiation from slots on elliptic cylinders, spheres, and spheroids.

Pathak, P. H.↗
NASA NTRS1993

Analysis of the tropospheric water distribution during FIRE 2

The Penn State/NCAR mesoscale model, as adapted for use at ARC, was used as a testbed for the development and validation of cloud models for use in General Circulation Models (GCM's). This modeling approach also allows us to intercompare the predictions of the various cloud schemes within the same dynamical framework. The use of the PSU/NCAR mesoscale model also allows us to compare our results with FIRE-II (First International Satellite Cloud Climatology Project Regional Experiment) observations, instead of climate statistics. Though a promising approach, our work to date revealed several difficulties. First, the model by design is limited in spatial coverage and is only run for 12 to 48 hours at a time. Hence the quality of the simulation will depend heavily on the initial conditions. The poor quality of upper-tropospheric measurements of water vapor is well known and the situation is particularly bad for mid-latitude winter since the coupling with the surface is less direct than in summer so that relying on the model to spin-up a reasonable moisture field is not always successful. Though one of the most common atmospheric constituents, water vapor is relatively difficult to measure accurately, especially operationally over large areas. The standard NWS sondes have little sensitivity at the low temperatures where cirrus form and the data from the GOES 6.7 micron channel is difficult to quantify. For this reason, the goals of FIRE Cirrus II included characterizing the three-dimensional distribution of water vapor and clouds. In studying the data from FIRE Cirrus II, it was found that no single special observation technique provides accurate regional distributions of water vapor. The Raman lidar provides accurate measurements, but only at the Hub, for levels up to 10 km, and during nighttime hours. The CLASS sondes are more sensitive to moisture at low temperatures than are the NWS sondes, but the four stations only cover an area of two hundred kilometers on a side. The aircraft give the most accurate measurements of water vapor, but are limited in spatial and temporal coverage. This problem is partly alleviated by the use of the MAPS analyses, a four-dimensional data assimilation system that combines the previous 3-hour forecast with the available observations, but its upper-level moisture analyses are sometimes deficient because of the vapor measurement problem. An attempt was made to create a consistent four-dimensional description of the water vapor distribution during the second IFO by subjectively combining data from a variety of sources, including MAPS analyses, CLASS sondes, SPECTRE sondes, NWS sondes, GOES satellite analyses, radars, lidars, and microwave radiometers.

Westphal, Douglas L.↗
NASA NTRS2000

Effects of Geometric Azimuthal Asymmetries of the PPM Stack on Electron Beam Characteristics

The effects of geometric azimuthally asymmetric properties of a periodic permanent magnet (PPM) focusing stack on electron beam characteristics obtained using a fully three dimensional (3D) particle-in-cell (PIC) code will be presented. The simulation model, using MAFIA (Solution of MAxwell's equations by the Finite-Integration-Algorithm), incorporates 3D behavior of the beam immersed in static fields calculated directly from the exact geometry and material properties of the 3D magnetic focusing structure. The Hughes 8916H, 18-40 GHz helical TWT for the millimeter-wave power module (MMPM) was used as a prototype. Firstly, the effects of C-magnets used at the input and output of the TWT to allow for coupling of the RF signal into and out of the tube are considered. The 8916H input and output C-magnets differ because coaxial couplers are used at the input and waveguide couplers are used at the output The repositioning of the beam from its central axis due to the inclusion of the output C-magnet was found to be most significant. The modeled output C-magnet and its orientation in the Cartesian coordinate system is shown, and a two-dimensional beam profile including the output C-magnet is also shown. A table presents the shift of the beam center off the central axis relative to the average radius of the beam at the longitudinal points A, B and C designated on an enclosed figure. Secondly, the addition of shunts, or rectangular iron pieces applied manually by a skilled technician in order to improve beam transmission, is considered. The shunts are applied to the top of the tube; thus, azimuthal symmetry of the focusing stack is interrupted. Although shunts are typically added during RF focusing, they are also typically added at the input section of the tube where RF forces are minimal, making an electron optics analysis meaningful. Because several shunts are usually applied to one pole piece, the simulations have been simplified by modeling a half washer with the same radius and longitudinal length as a shunt over the entire x, positive-y half of the transverse plane. A modeled pole piece and shunt as described are shown. Lastly, in order to study the effects of magnet misalignments, a magnet in the PPM stack was arbitrarily chosen and adjusted so that its central axis was shifted both 0.7 percent and 1.0 percent of the magnet outer diameter in the positive-y direction. In practice, positioning the magnets so that their central axis is accurately aligned with the central axis of the tube is challenging. Thus, it is a strong possibility that one or more magnets will be misaligned relative to the tube central axis.

Kory, Carol L.↗
NASA NTRS1976

Flutter and buckling of general laminated plates

An anisotropic buckling and flutter analysis is developed with allowance for both bending-extensional coupling and bending-twisting coupling within the framework of linear small deflection theory for simply supported general laminated plates. The extended Galerkin method is used to obtain approximate solutions to the coupled governing equations. The effects of various anisotropic stiffness parameters on the static and dynamic stability of laminated plates are evaluated, with particular emphasis on assessing the range of applicability of classical orthotropic plate theory. It is shown that bending-extensional coupling and bending-twisting stiffness terms have a destabilizing effect on buckling and flutter, the effect being more pronounced for a small number of layers. For symmetric plates, the number of layers required for orthotropic plate theory to be applicable is generally less for the buckling problem than for flutter. For square plates, aligning the fibers with the direction of airflow over the plate surface results in the highest flutter dynamic pressure.

Sawyer, J. W.↗
NASA NTRS1983

Non-planar MHD model for solar flare-generated disturbances in the heliospheric equatorial plane

An analysis, with a representative (canonical) example of solar-flare-generated equatorial disturbances, is made for the temporal and spatial changes in the solar wind plasma and magnetic field environment between the sun and 1 AU. The goal is to search for first-order global consequences rather than to make a parametric study. The analysis treats all three plasma velocity and magnetic field components in any convenient heliospheric plane of symmetry. The representative disturbance is examined for the canonical case in which the temporal and spatial changes in a homogeneous solar wind caused by a solar-flare-generated shock wave are described. All plasma and field parameters at three radial locations are examined. These are the central meridian and 33 deg W and 90 deg W of the flare's central meridian. It is found that the incorporation of a small meridional magnetic field in the ambient magnetic spiral field has negligible effect on the results. The magnetic field exhibits strong kinking within the interplanetary shocked flow, even reversed polarity that, coupled with low temperature and low density, suggests a plausible explanation for magnetic clouds' with accompanying double-streaming of electrons observed at directions approximately 90 deg to the heliocentric radius.

Wu, S. T.↗
NASA NTRS1992

PDF turbulence modeling and DNS

The problem of time discontinuity (or jump condition) in the coalescence/dispersion (C/D) mixing model is addressed in probability density function (pdf). A C/D mixing model continuous in time is introduced. With the continuous mixing model, the process of chemical reaction can be fully coupled with mixing. In the case of homogeneous turbulence decay, the new model predicts a pdf very close to a Gaussian distribution, with finite higher moments also close to that of a Gaussian distribution. Results from the continuous mixing model are compared with both experimental data and numerical results from conventional C/D models. The effect of Coriolis forces on compressible homogeneous turbulence is studied using direct numerical simulation (DNS). The numerical method used in this study is an eight order compact difference scheme. Contrary to the conclusions reached by previous DNS studies on incompressible isotropic turbulence, the present results show that the Coriolis force increases the dissipation rate of turbulent kinetic energy, and that anisotropy develops as the Coriolis force increases. The Taylor-Proudman theory does apply since the derivatives in the direction of the rotation axis vanishes rapidly. A closer analysis reveals that the dissipation rate of the incompressible component of the turbulent kinetic energy indeed decreases with a higher rotation rate, consistent with incompressible flow simulations (Bardina), while the dissipation rate of the compressible part increases; the net gain is positive. Inertial waves are observed in the simulation results.

Hsu, A. T.↗
NASA NTRS2000

An Improved Neutron Transport Algorithm for Space Radiation

A low-energy neutron transport algorithm for use in space radiation protection is developed. The algorithm is based upon a multigroup analysis of the straight-ahead Boltzmann equation by using a mean value theorem for integrals. This analysis is accomplished by solving a realistic but simplified neutron transport test problem. The test problem is analyzed by using numerical and analytical procedures to obtain an accurate solution within specified error bounds. Results from the test problem are then used for determining mean values associated with rescattering terms that are associated with a multigroup solution of the straight-ahead Boltzmann equation. The algorithm is then coupled to the Langley HZETRN code through the evaporation source term. Evaluation of the neutron fluence generated by the solar particle event of February 23, 1956, for a water and an aluminum-water shield-target configuration is then compared with LAHET and MCNPX Monte Carlo code calculations for the same shield-target configuration. The algorithm developed showed a great improvement in results over the unmodified HZETRN solution. In addition, a two-directional solution of the evaporation source showed even further improvement of the fluence near the front of the water target where diffusion from the front surface is important.

Heinbockel, John H.↗
NASA NTRS2014

Diagnosis of Middle Atmosphere Climate Sensitivity by the Climate Feedback Response Analysis Method

We present a new method to diagnose the middle atmosphere climate sensitivity by extending the Climate Feedback-Response Analysis Method (CFRAM) for the coupled atmosphere-surface system to the middle atmosphere. The Middle atmosphere CFRAM (MCFRAM) is built on the atmospheric energy equation per unit mass with radiative heating and cooling rates as its major thermal energy sources. MCFRAM preserves the CFRAM unique feature of an additive property for which the sum of all partial temperature changes due to variations in external forcing and feedback processes equals the observed temperature change. In addition, MCFRAM establishes a physical relationship of radiative damping between the energy perturbations associated with various feedback processes and temperature perturbations associated with thermal responses. MCFRAM is applied to both measurements and model output fields to diagnose the middle atmosphere climate sensitivity. It is found that the largest component of the middle atmosphere temperature response to the 11-year solar cycle (solar maximum vs. solar minimum) is directly from the partial temperature change due to the variation of the input solar flux. Increasing CO2 always cools the middle atmosphere with time whereas partial temperature change due to O3 variation could be either positive or negative. The partial temperature changes due to different feedbacks show distinctly different spatial patterns. The thermally driven globally averaged partial temperature change due to all radiative processes is approximately equal to the observed temperature change, ranging from 0.5 K near 70 km from the near solar maximum to the solar minimum.

MCFRAM↗
NASA NTRS2012

In Situ Geochemical Analysis and Age Dating of Rocks Using Laser Ablation-Miniature Mass Spectrometer

A miniaturized instrument for performing chemical and isotopic analysis of rocks has been developed. The rock sample is ablated by a laser and the neutral species produced are analyzed using the JPL-invented miniature mass spectrometer. The direct sampling of neutral ablated material and the simultaneous measurement of all the elemental and isotopic species are the novelties of this method. In this laser ablation-miniature mass spectrometer (LA-MMS) method, the ablated neutral atoms are led into the electron impact ionization source of the MMS, where they are ionized by a 70-eV electron beam. This results in a secondary ion pulse typically 10-100 microsecond wide, compared to the original 5-10-nanosecond laser pulse duration. Ions of different masses are then spatially dispersed along the focal plane of the magnetic sector of the miniature mass spectrometer and measured in parallel by a modified CCD (charge-coupled device) array detector capable of detecting ions directly. Compared to conventional scanning techniques, simultaneous measurement of the ion pulse along the focal plane effectively offers a 100% duty cycle over a wide mass range. LAMMS offers a more quantitative assessment of elemental composition than techniques that detect laser-ionized species produced directly in the ablation process because the latter can be strongly influenced by matrix effects that vary with the structure and geometry of the surface, the laser beam, and the ionization energies of the elements. The measurement of high-precision isotopic ratios and elemental composition of different rock minerals by LAMMS method has been demonstrated. The LA-MMS can be applied for the absolute age determination of rocks. There is no such instrument available presently in a miniaturized version that can be used for NASA space missions. Work is in progress in the laboratory for geochronology of rocks using LA-MMS that is based on K-Ar radiogenic dating technique.

Sinha, Mahadeva P.↗
NASA NTRS1996

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↗