Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element methods”

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 1,405 records · Page 78

Nonlinear shell analyses of the SRB/ETA ring interface

Two-dimensional shell models of the portion of the solid rocket booster (SRB) with the external tank attachment (ETA) ring were developed using the STAGSC-1 computer code. These analyses were performed in support of the SRB recertification program underway at the NASA Marshall Space Flight Center. The purpose of these analyses is to calculate the overall stress and deflection distributions for the SRB/ETA ring interface. An overview of the SRB/ETA ring geometry is presented followed by a discussion of the analysis methods and finite element models. Finally, the nonlinear structural response of the SRB/ETA ring interface is discussed.

Knight, Norman F., Jr.↗

Dynamics of runaway tails with time-dependent sub-Dreicer dc fields in magnetized plasmas

The evolution of runaway tails driven by sub-Dreicer time-dependent dc fields in a magnetized plasma are studied numerically using a quasi-linear code based on the Ritz-Galerkin method and finite elements. It is found that the runaway tail maintained a negative slope during the dc field increase. Depending on the values of the dc electric field at t = 0 and the electron gyrofrequency to the plasma frequency ratio the runaway tail became unstable to the anomalous Doppler resonance or remained stable before the saturation of the dc field at some maximum value. The systems that remained stable during this stage became unstable to the anomalous Doppler or the Cerenkov resonances when the dc field was kept at the saturation level or decreased. Once the instability is triggered, the runaway tail is isotropized.

Moghaddam-Taaheri, E.↗

Effect of a shallow low-viscosity zone on the formation of midplate swells

A finite-element numerical method is used here to analyze the effect of a low-viscosity zone on convection driven by heating from below in the upper mantle, in particular on the formation of midplate swells. The convective temperature and velocity solutions are calculated for different combinations of the viscosity in the top layer, the fluid layer thicknesses, and the Rayleigh number based on the viscosity in the bottom layer. The temperature solutions are used to calculate the geoid, topography, and heat flow anomalies, the elastic plate thickness, the depth of compensation, and an upper bound on the uplift time that result from the flow. The results are compared to data at the Hawaii, Bermuda, Cape Verde, and Marquesas swells. The magnitudes and the trend with age are consistent with theoretical and other estimates of the viscosity variation in the shallow upper mantle. Convective models can therefore explain the uplift and observed anomalies at midplate swells.

Robinson, E. M.↗

Axisymmetric shapes and stability of charged drops in an external electric field

A highly conducting charged drop that is surrounded by a fluid insulator of another density can be levitated by suitably applying a uniform electric field. Axisymmetric equilibrium shapes and stability of the levitated drop are found by solving simultaneously the augmented Young-Laplace equation for surface shape and the Laplace equation for the elecric field, together with constraints of fixed drop volume, charge, and center of mass. The means are a method of subdomains, finite element basis functions, and Galerkin's method of weighted residuals, all facilitated by a large-scale computer. Shape families of fixed charge are treated systematically by first-order continuation. Previous analyses by Abbas et al. in 1967 and Abbas and Latham in 1969, in which the shapes of levitated drops are approximated as spheroids, are corrected. The new analysis shows that drops charged to less than the Rayleigh limit lose shape stability at turning points, with respect to external field strength, and that the instability seen in experiments of Doyle et al. in 1964 and others is not a bifurcation to a family of two-lobed shapes, but rather is a related imperfect bifurcation.

Basaran, O. A.↗

Multibody instantly centered moving Lagrangian observer schemes. I - Formulation. II - Application to vehicular simulations

A finite element modeling method employing an instantly centered moving Lagrangian observer is proposed to solve multibody problems involving several rotating components, each with its own rotational history. The present technique allows the steady-state behavior of the problem simulation to be rendered stationary and time independent. In the second part, the method is extended to model the steady and transient response of ground-based automotive-type vehicular systems, including the modeling of vehicular obstruction rollover events.

Padovan, Joseph↗

Innovative design of composite structures: The use of curvilinear fiber format in structural design of composites

The influence is investigated of a curvilinear fiber format on load carrying capacity of a layered fiber reinforced plate with a centrally located hole. A curvilinear fiber format is descriptive of layers in a laminate having fibers which are aligned with the principal stress directions in those layers. Laminates of five curvilinear fiber format designs and four straightline fiber format designs are considered. A quasi-isotropic laminate having a straightline fiber format is used to define a baseline design for comparison with the other laminate designs. Four different plate geometries are considered and differentiated by two values of hole diameter/plate width equal to 1/6 and 1/3, and two values of plate length/plate width equal to 2 and 1. With the plates under uniaxial tensile loading on two opposing edges, alignment of fibers in the curvilinear layers with the principal stress directions is determined analytically by an iteration procedure. In-plane tensile load capacity is computed for all of the laminate designs using a finite element analysis method. A maximum strain failure criterion and the Tsai-Wu failure criterion are applied to determine failure loads and failure modes. Resistance to buckling of the laminate designs to uniaxial compressive loading is analyzed using the commercial code Engineering Analysis Language. Results indicate that the curvilinear fiber format laminates have higher in-plane tensile load capacity and comparable buckling resistance relative to the straightline fiber format laminates.

Charette, R. F.↗

Calculation of flight vibration levels of the AH-1G helicopter and correlation with existing flight vibration measurements

NASA-Langley is sponsoring a rotorcraft structural dynamics program with the objective to establish in the U.S. a superior capability to utilize finite element analysis models for calculations to support industrial design of helicopter airframe structures. In the initial phase of the program, teams from the major U.S. manufacturers of helicopter airframes will apply extant finite element analysis methods to calculate loads and vibrations of helicopter airframes, and perform correlations between analysis and measurements. The aforementioned rotorcraft structural dynamics program was given the acronym DAMVIBS (Design Analysis Method for Vibrations). Sikorsky's RDYNE Rotorcraft Dynamics Analysis used for the correlation study, the specifics of the application of RDYNE to the AH-1G, and comparisons of the predictions of the method with flight data for loads and vibrations on the AH-1G are described. RDYNE was able to predict trends of variations of loads and vibrations with airspeed, but in some instances magnitudes differed from measured results by factors of two or three to one. Sensitivities were studied of predictions to rotor inflow modeling, effects of torsional modes, number of blade bending modes, fuselage structural damping, and hub modal content.

Sopher, R.↗

Acceleration of runaway electrons in solar flares

The dc electric field acceleration of electrons out of a thermal plasma and the evolution of the runaway tail are studied numerically, using a relativistic quasi-linear code based on the Ritz-Galerkin method and finite elements. A small field-aligned electric field is turned on at a certain time. The resulting distribution function from the runaway process is used to calculate the synchrotron emission during the evolution of the runaway tail. It is found that, during the runaway tail formation, which lasts a few tens of seconds for typical solar flare conditions, the synchrotron emission level is low, almost ot the same order as the emission from the thermal plasma, at the high-frequency end of the spectrum. However, the emission is enhanced explosively in a few microseconds by several orders of magnitude at the time the runaway tail stops growing along the magnetic field and tends toward isotropy due to the pitch-angle scattering of the fast particles. Results indicate that, in order to account for the observed synchrotron emission spectrum of a typical solar flare, the electric field acceleration phase must be accompanied or preceded by a heating phase which yields an enhanced electron temperature of about 2-15 keV in the flare region if the electric field is 0.1-0.2 times the Dreicer field and cyclotron-to-plasma frequency ratios are of order 1-2.

Moghaddam-Taaheri, E.↗

Aeroelastic effects on stability and control of hingeless rotor helicopters

The combined effect of torsional flexibility and offset of the aerodynamic center relative to the elastic axis on the stability and control of a hingeless rotor helicopter in forward flight is studied in this paper. The aeroelastic model of the blade is based on fully coupled flap-lag-torsion equations of motion, which include kinematic nonlinearities due to moderate deflections. The equations are discretized using a finite element Galerkin method in space, and a classical Galerkin method in time. The vehicle trim calculations are coupled to the blade aeroelastic response calculations. Quasilinearization is used to compute aeroelastic stability. Reducing torsional stiffness and moving the elastic axis ahead of the aerodynamic center strongly stabilizes the phugoid mode without reducing control sensitivity. Aeroelastic stability calculations must accompany flight dynamics calculations to prevent aeroelastic instabilities that might easily go undetected.

Celi, Roberto↗

Preliminary structural design of a lunar transfer vehicle aerobrake

An aerobrake concept for a Lunar transfer vehicle was weight optimized through the use of the Taguchi design method, structural finite element analyses and structural sizing routines. Six design parameters were chosen to represent the aerobrake structural configuration. The design parameters included honeycomb core thickness, diameter to depth ratio, shape, material, number of concentric ring frames, and number of radial frames. Each parameter was assigned three levels. The minimum weight aerobrake configuration resulting from the study was approx. half the weight of the average of all twenty seven experimental configurations. The parameters having the most significant impact on the aerobrake structural weight were identified.

Bush, Lance B.↗

Preliminary structural design of a lunar transfer vehicle aerobrake

An aerobrake for a lunar transfer vehicle was designed through the use of the Taguchi design method, structural finite element analyses, and structural sizing routines. A minimum weight aerobrake structural configuration was the objective of the study. Six design parameters were chosen to represent the aerobrake structural configuration. The design parameters included honeycomb core thickness, diameter-to-depth ratio, shape, material, number of concentric ring frames and number of radial frames. Each parameter had 3 levels. The optimum aerobrake configuration resulting from the study was approximately half the weight of the average of all the other experimental configurations. The parameters having the most significant impact on the aerobrake structural weight were identified.

Bush, Lance B.↗

Detailed analysis and test correlation of a stiffened composite wing panel

New finite element analysis methods are examined by application to a complicated composite wing panel from the V-22 rotorcraft. A detailed FEM model with a relatively coarse mesh of 9-node elements was generated, and linear and nonlinear stress analyses, first-ply failure analyses, and buckling analyses were conducted. At low values of applied load, i.e., up to the design ultimate load of the panel, the linear stress analysis accurately predicted the strains and structural response characteristics of the panel.

Davis, D. D., Jr.↗

Capacitive skin characterization

NASA is currently involved in research that utilizes a capacitive sensor that is used for proximity detection of objects. This sensor is sensitive to conductive and dielectric materials including metal objects and humans. The range of the sensor has been found to be about twelve inches. It is the goal of this research project to further characterize the sensor so that it can be tailored to specific requirements. The characterization of the sensor should be with respect to shield size, sensor size, object size, and object distance. The method of finite elements to calculate the capacitance of the sensor while varying different parameters was used. Each of the parameters was varied in turn, often by selecting data points from different runs. The plotted results are shown and an apparent functionality developed for each.

Mcconnell, Robert↗

Scattering of ultrasonic wave by cracks in a plate

A hybrid numerical method combining finite elements and the boundary integral representation is used to investigate the transient scattering of ultrasonic waves by a crack in a plate. The incident wave models the guided waves generated by a steel ball impact on the plate. Two surface-breaking cracks and one subsurface crack are studied here. The results show that the location and depth of cracks have measurable effects on the surface responses in time and frequency domains. Also, the scattered fields have distinct differences in the three cases.

Liu, S. W.↗

Recent advances in the structual analysis of scientific balloons

The analysis of balloon envelopes by the finite element (FE) method is plagued by convergence problems. A pratical FE analysis approach is based on the fact that in thin shells with non-zero Gaussian curvature the membrane solution component is essentially decoupled from the bending solution component. A proxy-problem is solved by using a small artificial bending stiffness that assures convergence without significantly affecting the membrane solution component. This approach has been previously validated on slightly overpressurized balloon envelopes. Extensions of this approach to more difficult problems in the structural analysis of balloon envelopes are presented. The convergence forcing modelling measures are discussed. Implications of the findings of the analysis results to future balloon designs are also discussed.

Schur, W. W.↗

Adaptive Finite-Element Computation In Fracture Mechanics

Report discusses recent progress in use of solution-adaptive finite-element computational methods to solve two-dimensional problems in linear elastic fracture mechanics. Method also shown extensible to three-dimensional problems.

Min, J. B.↗

Convergence of a Substructuring Method with LaGrange Multipliers

We analyze the convergence of a substructuring iterative method with Lagrange multipliers, proposed recently by Farhat and Roux. The method decomposes finite element discretization of an elliptic boundary value problem into Neumann problems on the subdomains and a coarse problem for the subdomain nullspace components. For linear conforming elements and preconditioning by the Dirichlet problems on the subdomains, we prove the asymptotic bound on the condition number C(1 + log(H/h))(sup gamma), gamma = 2 or 3, where h is the characteristic element size and H is the subdomain size.

Mandel, Jan↗

Design Oriented Structural Modeling for Airplane Conceptual Design Optimization

The main goal for research conducted with the support of this grant was to develop design oriented structural optimization methods for the conceptual design of airplanes. Traditionally in conceptual design airframe weight is estimated based on statistical equations developed over years of fitting airplane weight data in data bases of similar existing air- planes. Utilization of such regression equations for the design of new airplanes can be justified only if the new air-planes use structural technology similar to the technology on the airplanes in those weight data bases. If any new structural technology is to be pursued or any new unconventional configurations designed the statistical weight equations cannot be used. In such cases any structural weight estimation must be based on rigorous "physics based" structural analysis and optimization of the airframes under consideration. Work under this grant progressed to explore airframe design-oriented structural optimization techniques along two lines of research: methods based on "fast" design oriented finite element technology and methods based on equivalent plate / equivalent shell models of airframes, in which the vehicle is modelled as an assembly of plate and shell components, each simulating a lifting surface or nacelle / fuselage pieces. Since response to changes in geometry are essential in conceptual design of airplanes, as well as the capability to optimize the shape itself, research supported by this grant sought to develop efficient techniques for parametrization of airplane shape and sensitivity analysis with respect to shape design variables. Towards the end of the grant period a prototype automated structural analysis code designed to work with the NASA Aircraft Synthesis conceptual design code ACS= was delivered to NASA Ames.

Livne, Eli↗