Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stabilized finite element method”

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 109 records · Page 6

Adaptive quadrilateral and triangular finite-element scheme for compressible flows

The development of an adaptive mesh refinement procedure for analyzing high-speed compressible flows using the finite-element method is described. This new adaptation procedure, which uses both quadrilateral and triangular elements, was implemented with two explicit finite-element algorithms - the two-step Taylor-Galerkin and the multistep Galerkin-Runge-Kutta schemes. A von Neumann stability analysis and a rotating 'cosine hill'problem demonstrate the instability of the Taylor-Galerkin scheme when coupled with the adaptation procedure. For the same adaptive refinement scheme, the Galerkin-Runge-Kutta procedure yields stable solutions within its explicit stability limit. The utility of this new adaptation procedure for the prediction of compressible flow features is illustrated for inviscid problems involving strong shock interactions at hypersonic speeds.

Ramakrishnan, R.↗

An Aeroelastic Analysis of a Thin Flexible Membrane

Studies have shown that significant vehicle mass and cost savings are possible with the use of ballutes for aero-capture. Through NASA's In-Space Propulsion program, a preliminary examination of ballute sensitivity to geometry and Reynolds number was conducted, and a single-pass coupling between an aero code and a finite element solver was used to assess the static aeroelastic effects. There remain, however, a variety of open questions regarding the dynamic aeroelastic stability of membrane structures for aero-capture, with the primary challenge being the prediction of the membrane flutter onset. The purpose of this paper is to describe and begin addressing these issues. The paper includes a review of the literature associated with the structural analysis of membranes and membrane utter. Flow/structure analysis coupling and hypersonic flow solver options are also discussed. An approach is proposed for tackling this problem that starts with a relatively simple geometry and develops and evaluates analysis methods and procedures. This preliminary study considers a computationally manageable 2-dimensional problem. The membrane structural models used in the paper include a nonlinear finite-difference model for static and dynamic analysis and a NASTRAN finite element membrane model for nonlinear static and linear normal modes analysis. Both structural models are coupled with a structured compressible flow solver for static aeroelastic analysis. For dynamic aeroelastic analyses, the NASTRAN normal modes are used in the structured compressible flow solver and 3rd order piston theories were used with the finite difference membrane model to simulate utter onset. Results from the various static and dynamic aeroelastic analyses are compared.

Scott, Robert C.↗

H Infinity Control of Magnetic Bearings to Ensure Both System and External Periodic Disturbance Robustness

Both self-excited and forced disturbances often lead to severe rotor vibrations in a magnetic bearing systems with long slender shafts. This problem has been studied using the H-infinity method, and stability with good robustness can be achieved for the linearized model of a magnetic bearing when small transient disturbances are applied. In this paper, the H-infinity control method for self-excited and forced disturbances is first reviewed. It is then applied to the control of a magnetic bearing rotor system. In modelling the system, the shaft is first discretized into 18 finite elements and then three levels of condensation are applied. This leads to a system with three masses and three compliant elements which can be described by six state variable coordinates. Simulation of the resultant system design has been performed at speeds up to 10,000 rpm. Disturbances in terms of different initial displacements, initial impulses, and external periodic inputs have been imposed. The simulation results show that good stability can be achieved under these different transient disturbances using the proposed controller while at the same time reducing the sensitivity to external periodic disturbances.

Jiang, Yuhong↗

Analysis of structures with rotating, flexible substructures applied to rotorcraft aeroelasticity in GRASP

Application to the General Rotorcraft Aeromechanical Stability Program (GRASP) of new methodology for structural dynamic analysis, including substructuring, frames of reference, nodes, finite elements, and constraints, is discussed. The structure is decomposed into a hierarchy of substructures, and discrete relative motion between substructures is analyzed exactly. The finite element method is used to treat deformation of continua, and the library of finite elements includes a nonlinear beam element incorporating aeroelastic effects. Analytical bases for the aeroelastic beam element and the screw constraint are considered, and the important role of geometric stiffness in the formulation is shown.

Hodges, Dewey H.↗

Numerical solution of the vertical structure equation in the normal mode method

In the present model of multilayered stability stratification, aimed at obtaining the analytic eigensolutions of the vertical structure equation, each layer is characterized by its own static stability value. By requiring continuity of pressure and vertical velocity across each interface level, and by imposing suitable upper and lower boundary conditions, matching eigensolutions are obtained in terms of the Bessel functions. Attention is given to an explicit example of a double-layered stratified atmosphere which demonstrates the mathematical manipulations involved; the resultant vertical structure functions are used to check the accuracy of the numerical solutions by the finite difference and finite element methods.

Sasaki, Y. K.↗

Profiles of electrified drops and bubbles

Axisymmetric equilibrium shapes of conducting drops and bubbles, (1) pendant or sessile on one face of a circular parallel-plate capacitor or (2) free and surface-charged, are found by solving simultaneously the free boundary problem consisting of the augmented Young-Laplace equation for surface shape and the Laplace equation for electrostatic field, given the surface potential. The problem is nonlinear and the method is a finite element algorithm employing Newton iteration, a modified frontal solver, and triangular as well as quadrilateral tessellations of the domain exterior to the drop in order to facilitate refined analysis of sharply curved drop tips seen in experiments. The stability limit predicted by this computer-aided theoretical analysis agrees well with experiments.

Basaran, O. A.↗

A Study of the Effect of Adhesive and Matrix Stiffnesses on the Axial, Normal, and Shear Stress Distributions of a Boron-epoxy Reinforced Composite Joint

The mechanical properties of a symmetrical, eight-step, titanium-boron-epoxy joint are discussed. A study of the effect of adhesive and matrix stiffnesses on the axial, normal, and shear stress distributions was made using the finite element method. The NASA Structural Analysis Program (NASTRAN) was used for the analysis. The elastic modulus of the adhesive was varied from 345 MPa to 3100 MPa with the nominal value of 1030 MPa as a standard. The nominal values were used to analyze the stability of the joint. The elastic moduli were varied to determine their effect on the stresses in the joint.

Howell, W. E.↗

Robust controller synthesis for a large flexible space antenna

The linear-quadratic Gaussian/loop-transfer-recovery method is used to synthesize a fine-pointing control system for a large space antenna. A finite-element model for the 122-m hoop/column antenna is employed, and a compensator, utilizing attitudes sensors and torque actuators, is designed which achieves pointing performance while maintaining stability robustness to unmodeled dynamics. Inclusion of the rigid-body modes plus the first three elastic modes is found to be necessary to achieve a 0.1-rad/s bandwidth. Results are obtained by employing a modification of the standard robustness recovery procedure, which reduces the conservative nature of the design methodology. Performance degradation is encountered due to the presence of unavoidable invariant zeros within the design bandwidth.

Sundararajan, N.↗

Structural Design Exploration of an Electric Powered Multi-Propulsor Wing Configuration

Advancements in aircraft electric propulsion may enable an expanded operational envelope for electrically powered vehicles compared to their internal combustion engine counterparts. High aspect ratio wings provide additional lift and drag reduction for a proposed multi-propulsor design, however, the challenge is to reduce the weight of wing structures while maintaining adequate structural and aeroelastic margins. Design exploration using a conventional design-and-build philosophy coupled with a finite element method (FEM)-based design of experiments (DOE) strategy are presented to examine high aspect ratio wing structures that have spanwise distributed electric motors. Multiple leading-edge-mounted engine masses presented a challenge to design a wing within acceptable limits for dynamic and aeroelastic stability. Because the first four primary bending eigenmodes of the proposed wing structure are very sensitive to outboard motor placement, safety-of-flight requirements drove the need for multiple spars, rib attachments, and outboard structural reinforcements in the design. Global aeroelasticity became an increasingly important design constraint during the on-going design process, with outboard motor pod flutter ultimately becoming a primary design constraint. Designers successively generated models to examine stress, dynamics, and aeroelasticity concurrently. This research specifically addressed satisfying multi-disciplinary design criteria to generate fluid-structure interaction solution sets, and produced high aspect ratio primary structure designs for the NASA Scalable Convergent Electric Propulsion Technology and Operations Research (SCEPTOR) project in the Aeronautic Research Mission Directorate at NASA. In this paper, a dynamics-driven, quasi-inverse design methodology is presented to address aerodynamic performance goals and structural challenges encountered for the SCEPTOR demonstrator vehicle. These results are compared with a traditional computer aided design based approach.

Moore, James B.↗

Stability of generally stiffened anisotropic noncircular cylinders

Continuous filament grid-stiffened structure is a stiffening concept that combines structural efficiency and damage tolerance. However, finite element design of such structures against buckling is expensive due to the complexities of the structure. An analytical model of such a structure is developed using a penalty method (artificial springs) with a first order shear deformation theory (FSDT). The buckling analysis under combined loadings is done using energy method with a penalty/Rayleigh-Ritz technique. The penalty/Rayleigh-Ritz approach is computationally less demanding when compared to the finite element solution and mesh generation. Apart from the published research works on buckling of stiffened plates and shells by finite element and finite strips, research works on buckling of stiffened plates and shells utilize three different approaches; smeared, column, and discrete approaches. The discrete approach considers the discrete effects of the stiffeners in the buckling behavior by modeling stiffeners as line of bending (EI) and torsion (GJ) stiffnesses on panel skin. Some local deformations are lost when stiffeners are modeled as (EI) and (GJ) stiffeners. This approach becomes difficult in the case of plate stiffened in more than two directions. Most of the work done using the discrete approach involved the Classical Plate Theory (CLPT) rather than the FSDT. We report on our formulation of a discrete approach coupled with a penalty formulation and FSDT.

Sobh, Nahil Atef↗

Modal forced response of propfans in yawed flow

A modal forced response method for propfans in yawed flow is presented. This capability exists in the Aeroelastic Stability and Response of Propfan (ASTROP3) code developed at the Lewis Research Center. The code uses three-dimensional steady and unsteady cascade aerodynamics by Williams and Hwang (1986) and a NASTRAN finite element model to represent the blade structure. In addition, many utility programs exist in ASTROP3 that help in both the preprocessing of the NASTRAN model and the postprocessing of modal response results. The postprocessing work that computes the blade vibratory displacements and stresses in yawed flow are highlighted here.

Narayanan, G. V.↗

Parallel aeroelastic computations for wing and wing-body configurations

The objective of this research is to develop computationally efficient methods for solving fluid-structural interaction problems by directly coupling finite difference Euler/Navier-Stokes equations for fluids and finite element dynamics equations for structures on parallel computers. This capability will significantly impact many aerospace projects of national importance such as Advanced Subsonic Civil Transport (ASCT), where the structural stability margin becomes very critical at the transonic region. This research effort will have direct impact on the High Performance Computing and Communication (HPCC) Program of NASA in the area of parallel computing.

Byun, Chansup↗

Finite element simulations of thermosolutal convection in vertical solidification of binary alloys

Dendritic vertical solidification of a binary alloy is modeled using the finite element method to assess the effect of thermosolutal convection in macrosegregation. The mathematical model assumes steady-state solidification with a planar, undeformable surface defined by the dendrite tips and the eutectic isotherm. The dendritic region is assumed to advance at a constant solidification velocity v. The stability of the modeled system has been investigated and nonlinear calculations performed that show finger-like convection when the system is unstable. Results for lead-tin alloys show that when the system is unstable, convection is only significant in the uppermost part of the mush and is entirely driven by convection in the bulk fluid.

Heinrich, J. C.↗

Aeroelasticity of wing and wing-body configurations on parallel computers

The objective of this research is to develop computationally efficient methods for solving aeroelasticity problems on parallel computers. Both uncoupled and coupled methods are studied in this research. For the uncoupled approach, the conventional U-g method is used to determine the flutter boundary. The generalized aerodynamic forces required are obtained by the pulse transfer-function analysis method. For the coupled approach, the fluid-structure interaction is obtained by directly coupling finite difference Euler/Navier-Stokes equations for fluids and finite element dynamics equations for structures. This capability will significantly impact many aerospace projects of national importance such as Advanced Subsonic Civil Transport (ASCT), where the structural stability margin becomes very critical at the transonic region. This research effort will have direct impact on the High Performance Computing and Communication (HPCC) Program of NASA in the area of parallel computing.

Byun, Chansup↗

Using Tabulated Experimental Data to Drive an Orthotropic Elasto-Plastic Three-Dimensional Model for Impact Analysis

An orthotropic elasto-plastic-damage three-dimensional model with tabulated input has been developed to analyze the impact response of composite materials. The theory has been implemented as MAT 213 into a tailored version of LS-DYNA being developed under a joint effort of the FAA and NASA and has the following features: (a) the theory addresses any composite architecture that can be experimentally characterized as an orthotropic material and includes rate and temperature sensitivities, (b) the formulation is applicable for solid as well as shell element implementations and utilizes input data in a tabulated form directly from processed experimental data, (c) deformation and damage mechanics are both accounted for within the material model, (d) failure criteria are established that are functions of strain and damage parameters, and mesh size dependence is included, and (e) the theory can be efficiently implemented into a commercial code for both sequential and parallel executions. The salient features of the theory as implemented in LS-DYNA are illustrated using a widely used composite - the T800S/3900-2B[P2352W-19] BMS8-276 Rev-H-Unitape fiber/resin unidirectional composite. First, the experimental tests to characterize the deformation, damage and failure parameters in the material behavior are discussed. Second, the MAT213 input model and implementation details are presented with particular attention given to procedures that have been incorporated to ensure that the yield surfaces in the rate and temperature dependent plasticity model are convex. Finally, the paper concludes with a validation test designed to test the stability, accuracy and efficiency of the implemented model.

Polymer Matrix Composites↗

A program to compute three-dimensional subsonic unsteady aerodynamic characteristics using the doublet lattice method, L216 (DUBFLEX). Volume 2: Supplemental system design and maintenance document

The information necessary for execution of the digital computer program L216 on the CDC 6600 is described. L216 characteristics are based on the doublet lattice method. Arbitrary aerodynamic configurations may be represented with combinations of nonplanar lifting surfaces composed of finite constant pressure panel elements, and axially summetric slender bodies composed of constant pressure line elements. Program input consists of configuration geometry, aerodynamic parameters, and modal data; output includes element geometry, pressure difference distributions, integrated aerodynamic coefficients, stability derivatives, generalized aerodynamic forces, and aerodynamic influence coefficient matrices. Optionally, modal data may be input on magnetic field (tape or disk), and certain geometric and aerodynamic output may be saved for subsequent use.

Harrison, B. A.↗

The dynamics and control of large flexible space structures, 3. Part A: Shape and orientation control of a platform in orbit using point actuators

The dynamics, attitude, and shape control of a large thin flexible square platform in orbit are studied. Attitude and shape control are assumed to result from actuators placed perpendicular to the main surface and one edge and their effect on the rigid body and elastic modes is modelled to first order. The equations of motion are linearized about three different nominal orientations: (1) the platform following the local vertical with its major surface perpendicular to the orbital plane; (2) the platform following the local horizontal with its major surface normal to the local vertical; and (3) the platform following the local vertical with its major surface perpendicular to the orbit normal. The stability of the uncontrolled system is investigated analytically. Once controllability is established for a set of actuator locations, control law development is based on decoupling, pole placement, and linear optimal control theory. Frequencies and elastic modal shape functions are obtained using a finite element computer algorithm, two different approximate analytical methods, and the results of the three methods compared.

Bainum, P. M.↗

Development of methodology for horizontal axis wind turbine dynamic analysis

Horizontal axis wind turbine dynamics were studied. The following findings are summarized: (1) review of the MOSTAS computer programs for dynamic analysis of horizontal axis wind turbines; (2) review of various analysis methods for rotating systems with periodic coefficients; (3) review of structural dynamics analysis tools for large wind turbine; (4) experiments for yaw characteristics of a rotating rotor; (5) development of a finite element model for rotors; (6) development of simple models for aeroelastics; and (7) development of simple models for stability and response of wind turbines on flexible towers.

Dugundji, J.↗