Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuous finite element”

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 307 records · Page 17

Interlaminar Stresses by Refined Beam Theories and the Sinc Method Based on Interpolation of Highest Derivative

Computation of interlaminar stresses from the higher-order shear and normal deformable beam theory and the refined zigzag theory was performed using the Sinc method based on Interpolation of Highest Derivative. The Sinc method based on Interpolation of Highest Derivative was proposed as an efficient method for determining through-the-thickness variations of interlaminar stresses from one- and two-dimensional analysis by integration of the equilibrium equations of three-dimensional elasticity. However, the use of traditional equivalent single layer theories often results in inaccuracies near the boundaries and when the lamina have extremely large differences in material properties. Interlaminar stresses in symmetric cross-ply laminated beams were obtained by solving the higher-order shear and normal deformable beam theory and the refined zigzag theory with the Sinc method based on Interpolation of Highest Derivative. Interlaminar stresses and bending stresses from the present approach were compared with a detailed finite element solution obtained by ABAQUS/Standard. The results illustrate the ease with which the Sinc method based on Interpolation of Highest Derivative can be used to obtain the through-the-thickness distributions of interlaminar stresses from the beam theories. Moreover, the results indicate that the refined zigzag theory is a substantial improvement over the Timoshenko beam theory due to the piecewise continuous displacement field which more accurately represents interlaminar discontinuities in the strain field. The higher-order shear and normal deformable beam theory more accurately captures the interlaminar stresses at the ends of the beam because it allows transverse normal strain. However, the continuous nature of the displacement field requires a large number of monomial terms before the interlaminar stresses are computed as accurately as the refined zigzag theory.

Slemp, Wesley C. H.↗

The Space-Time Conservative Schemes for Large-Scale, Time-Accurate Flow Simulations with Tetrahedral Meshes

Despite decades of development of unstructured mesh methods, high-fidelity time-accurate simulations are still predominantly carried out on structured, or unstructured hexahedral meshes by using high-order finite-difference, weighted essentially non-oscillatory (WENO), or hybrid schemes formed by their combinations. In this work, the space-time conservation element solution element (CESE) method is used to simulate several flow problems including supersonic jet/shock interaction and its impact on launch vehicle acoustics, and direct numerical simulations of turbulent flows using tetrahedral meshes. This paper provides a status report for the continuing development of the space-time conservation element solution element (CESE) numerical and software framework under the Revolutionary Computational Aerosciences (RCA) project. Solution accuracy and large-scale parallel performance of the numerical framework is assessed with the goal of providing a viable paradigm for future high-fidelity flow physics simulations.

Venkatachari, Balaji Shankar↗

Development of a Computationally Efficient, High Fidelity, Finite Element Based Hall Thruster Model

This report documents the development of a two dimensional finite element based numerical model for efficient characterization of the Hall thruster plasma dynamics in the framework of multi-fluid model. Effect of the ionization and the recombination has been included in the present model. Based on the experimental data, a third order polynomial in electron temperature is used to calculate the ionization rate. The neutral dynamics is included only through the neutral continuity equation in the presence of a uniform neutral flow. The electrons are modeled as magnetized and hot, whereas ions are assumed magnetized and cold. The dynamics of Hall thruster is also investigated in the presence of plasma-wall interaction. The plasma-wall interaction is a function of wall potential, which in turn is determined by the secondary electron emission and sputtering yield. The effect of secondary electron emission and sputter yield has been considered simultaneously, Simulation results are interpreted in the light of experimental observations and available numerical solutions in the literature.

Jacobson, David↗

Modeling and Analysis of Wrinkled Membranes: An Overview

Thin-film membranes are basic elements of a variety of space inflatable/deployable structures. Wrinkling degrades the performance and reliability of these membrane structures, and hence has been a topic of continued interest. Wrinkling analysis of membranes for general geometry and arbitrary boundary conditions is quite challenging. The objective of this presentation is two-fold. Firstly, the existing models of wrinkled membranes and related numerical solution methods are reviewed. The important issues to be discussed are the capability of a membrane model to characterize taut, wrinkled and slack states of membranes in a consistent and physically reasonable manner; the ability of a wrinkling analysis method to predict the formation and growth of wrinkled regions, and to determine out-of-plane deformation and wrinkled waves; the convergence of a numerical solution method for wrinkling analysis; and the compatibility of a wrinkling analysis with general-purpose finite element codes. According to this review, several opening issues in modeling and analysis of wrinkled membranes that are to be addressed in future research are summarized, The second objective of this presentation is to discuss a newly developed membrane model of two viable parameters (2-VP model) and associated parametric finite element method (PFEM) for wrinkling analysis are introduced. The innovations and advantages of the proposed membrane model and PFEM-based wrinkling analysis are: (1) Via a unified stress-strain relation; the 2-VP model treat the taut, wrinkled, and slack states of membranes consistently; (2) The PFEM-based wrinkling analysis has guaranteed convergence; (3) The 2-VP model along with PFEM is capable of predicting membrane out-of-plane deformations; and (4) The PFEM can be integrated into any existing finite element code. Preliminary numerical examples are also included in this presentation to demonstrate the 2-VP model and PFEM-based wrinkling analysis approach.

Yang, B.↗

Integrated Modeling Activities for the James Webb Space Telescope: Structural-Thermal-Optical Analysis

The James Web Space Telescope (JWST) is a large, infrared-optimized space telescope scheduled for launch in 2011. This is a continuation of a series of papers on modeling activities for JWST. The structural-thermal-optical, often referred to as STOP, analysis process is used to predict the effect of thermal distortion on optical performance. The benchmark STOP analysis for JWST assesses the effect of an observatory slew on wavefront error. Temperatures predicted using geometric and thermal math models are mapped to a structural finite element model in order to predict thermally induced deformations. Motions and deformations at optical surfaces are then input to optical models, and optical performance is predicted using either an optical ray trace or a linear optical analysis tool. In addition to baseline performance predictions, a process for performing sensitivity studies to assess modeling uncertainties is described.

Johnston, John D.↗

Simulating Space Capsule Water Landing with Explicit Finite Element Method

A study of using an explicit nonlinear dynamic finite element code for simulating the water landing of a space capsule was performed. The finite element model contains Lagrangian shell elements for the space capsule and Eulerian solid elements for the water and air. An Arbitrary Lagrangian Eulerian (ALE) solver and a penalty coupling method were used for predicting the fluid and structure interaction forces. The space capsule was first assumed to be rigid, so the numerical results could be correlated with closed form solutions. The water and air meshes were continuously refined until the solution was converged. The converged maximum deceleration predicted is bounded by the classical von Karman and Wagner solutions and is considered to be an adequate solution. The refined water and air meshes were then used in the models for simulating the water landing of a capsule model that has a flexible bottom. For small pitch angle cases, the maximum deceleration from the flexible capsule model was found to be significantly greater than the maximum deceleration obtained from the corresponding rigid model. For large pitch angle cases, the difference between the maximum deceleration of the flexible model and that of its corresponding rigid model is smaller. Test data of Apollo space capsules with a flexible heat shield qualitatively support the findings presented in this paper.

Wang, John T.↗

Design and Evaluation of Modifications to the NASA Langley Flow Impedance Tube

The need to minimize fan noise radiation from commercial aircraft engine nacelles continues to provide an impetus for developing new acoustic liner concepts. If the full value of such concepts is to be attained, an understanding of grazing flow effects is crucial. Because of this need for improved understanding of grazing flow effects, the NASA Langley Research Center Liner Physics Group has invested a large effort over the past decade into the development of a 2-D finite element method that characterizes wave propagation through a lined duct. The original test section in the Langley Grazing IncidenceTube was used to acquire data needed for implementation of this finite element method. This test section employed a stepper motor-driven axial-traversing bar, embedded in the wall opposite the test liner, to position a flush-mounted microphone at pre-selected locations. Complex acoustic pressure data acquired with this traversing microphone were used to educe the acoustic impedance of test liners using this 2-D finite element method and a local optimization technique. Results acquired in this facility have been extensively reported, and were compared with corresponding results from various U.S. aeroacoustics laboratories in the late 1990 s. Impedance data comparisons acquired from this multi-laboratory study suggested that it would be valuable to incorporate more realistic 3-D aeroacoustic effects into the impedance eduction methodology. This paper provides a description of modifications that have been implemented to facilitate studies of 3-D effects. The two key features of the modified test section are (1) the replacement of the traversing bar and its flush-mounted microphone with an array of 95 fixed-location microphones that are flush-mounted in all four walls of the duct, and (2) the inclusion of a suction device to modify the boundary layer upstream of the lined portion of the duct. The initial results achieved with the modified test section are provided in this report, and a comparison of these results with those achieved using the original test section is used to demonstrate that the data acquisition and analysis with the new test section can be confidently used for impedance eduction.

Jones, Michael G.↗

Aeroelastic Analysis of a Flexible Wing Wind Tunnel Model with Variable Camber Continuous Trailing Edge Flap Design

This paper presents data analysis of a flexible wing wind tunnel model with a variable camber continuous trailing edge flap (VCCTEF) design for drag minimization tested at the University of Washington Aeronautical Laboratory (UWAL). The wind tunnel test was designed to explore the relative merit of the VCCTEF concept for improved cruise efficiency through the use of low-cost aeroelastic model test techniques. The flexible wing model is a 10%-scale model of a typical transport wing and is constructed of woven fabric composites and foam core. The wing structural stiffness in bending is tailored to be half of the stiffness of a Boeing 757-era transport wing while the torsional stiffness is about the same. This stiffness reduction results in a wing tip deflection of about 10% of the wing semi-span. The VCCTEF is a multi-segment flap design having three chordwise camber segments and five spanwise flap sections for a total of 15 individual flap elements. The three chordwise camber segments can be positioned appropriately to create a desired trailing edge camber. Elastomeric material is used to cover the gaps in between the spanwise flap sections, thereby creating a continuous trailing edge. Wind tunnel data analysis conducted previously shows that the VCCTEF can achieve a drag reduction of up to 6.31% and an improvement in the lift-to-drag ratio (L=D) of up to 4.85%. A method for estimating the bending and torsional stiffnesses of the flexible wingUWAL wind tunnel model from static load test data is presented. The resulting estimation indicates that the stiffness of the flexible wing is significantly stiffer in torsion than in bending by as much as 9 to 1. The lift prediction for the flexible wing is computed by a coupled aerodynamic-structural model. The coupled model is developed by coupling a conceptual aerodynamic tool Vorlax with a finite-element model of the flexible wing via an automated geometry deformation tool. Based on the comparison of the lift curve slope, the lift prediction for the rigid wing is in good agreement with the estimated lift coefficients derived from the wind tunnel test data. Due to the movement of the VCCTEF during the wind tunnel test, uncertainty in the lift prediction due to the indicated variations of the VCCTEF deflection is studied. The results show a significant spread in the lift prediction which contradicts the consistency in the aerodynamic measurements, thus suggesting that the indicated variations as measured by the VICON system may not be reliable. The lift prediction of the flexible wing agrees very well with the measured lift curve for the baseline configuration. The computed bending deflection and wash-out twist of the flexible wing also match reasonably well with the aeroelastic deflection measurements. The results demonstrate the validity of the aerodynamic-structural tool for use to analyze aerodynamic performance of flexible wings.

Wind Tunnel↗

A sequential linear optimization approach for controller design

A linear optimization approach with a simple real arithmetic algorithm is presented for reliable controller design and vibration suppression of flexible structures. Using first order sensitivity of the system eigenvalues with respect to the design parameters in conjunction with a continuation procedure, the method converts a nonlinear optimization problem into a maximization problem with linear inequality constraints. The method of linear programming is then applied to solve the converted linear optimization problem. The general efficiency of the linear programming approach allows the method to handle structural optimization problems with a large number of inequality constraints on the design vector. The method is demonstrated using a truss beam finite element model for the optimal sizing and placement of active/passive-structural members for damping augmentation. Results using both the sequential linear optimization approach and nonlinear optimization are presented and compared. The insensitivity to initial conditions of the linear optimization approach is also demonstrated.

Horta, L. G.↗

A Parametric Finite-Element Model for Evaluating Segmented Mirrors with Discrete, Edgewise Connectivity

Since future astrophysics missions require space telescopes with apertures of at least 10 meters, there is a need for on-orbit assembly methods that decouple the size of the primary mirror from the choice of launch vehicle. One option is to connect the segments edgewise using mechanisms analogous to damped springs. To evaluate the feasibility of this approach, a parametric ANSYS model that calculates the mode shapes, natural frequencies, and disturbance response of such a mirror, as well as of the equivalent monolithic mirror, has been developed. This model constructs a mirror using rings of hexagonal segments that are either connected continuously along the edges (to form a monolith) or at discrete locations corresponding to the mechanism locations (to form a segmented mirror). As an example, this paper presents the case of a mirror whose segments are connected edgewise by mechanisms analogous to a set of four collocated single-degree-of-freedom damped springs. The results of a set of parameter studies suggest that such mechanisms can be used to create a 15-m segmented mirror that behaves similarly to a monolith, although fully predicting the segmented mirror performance would require incorporating measured mechanism properties into the model. Keywords: segmented mirror, edgewise connectivity, space telescope

Gersh-Range, Jessica A.↗

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

The Use of the STAGS Finite Element Code in Stitched Structures Development

In the last 30 years NASA has worked in collaboration with industry to develop enabling technologies needed to make aircraft more fuel-efficient and more affordable. The focus on the airframe has been to reduce weight, improve damage tolerance and better understand structural behavior under realistic flight and ground loading conditions. Stitched structure is a technology that can address the weight savings, cost reduction, and damage tolerance goals, but only if it is supported by accurate analytical techniques. Development of stitched technology began in the 1990's as a partnership between NASA and Boeing (McDonnell Douglas at the time) under the Advanced Composites Technology Program and has continued under various titles and programs and into the Environmentally Responsible Aviation Project today. These programs contained development efforts involving manufacturing development, design, detailed analysis, and testing. Each phase of development, from coupons to large aircraft components was supported by detailed analysis to prove that the behavior of these structures was well-understood and predictable. The Structural Analysis of General Shells (STAGS) computer code was a critical tool used in the development of many stitched structures. As a key developer of STAGS, Charles Rankin's contribution to the programs was quite significant. Key features of STAGS used in these analyses and discussed in this paper include its accurate nonlinear and post-buckling capabilities, its ability to predict damage growth, and the use of Lagrange constraints and follower forces.

Jegley, Dawn C.↗

Analysis of delamination onset and growth in curved laminates

Delamination onset and growth from a matrix crack in a curved composite laminate are studied. Finite element analysis is used to determine values of G which can be compared to the interlaminar fracture toughness of the composite for delamination predictions. The delamination is assumed to grow from the matrix crack and to remain straight and perpendicular to the edge. G is larger at the free edge than in the center and increases with delamination growth. Thus, once delamination begins from the matrix crack it should continue around the curved section with no increase in applied load. In the local vicinity of the matrix crack the rate of change of G decreases initially and then begins to increase with a point of inflection in between. Measured failure loads and the interlaminar fracture toughness of the composite are used to show that the point of inflection corresponds to the critical value of G.

Martin, Roderick H.↗

Probabilistic system identification in the time domain

The objective of system identification is to determine reliable dynamical models of a structure by systematically using its measured excitation and response. It brings together in an integrated fashion, experimental, analytical, and computational techniques in structural dynamics. Areas of application for system identification include the following: (1) Model Evaluation--assessing assumptions (linearity and equivalent viscous damping) and techniques (finite-element modeling) used to construct theoretical models of a structure; (2) Model Improvement--updating of a theoretical model to enable more accurate response predictions for possible future loads on the structure, or for control of the structure; (3) Empirical Modelling--developing empirical relationships (nonlinear models) or empirical parameter values (modal damping) because the present state of the art does not provide theoretical results; and (4) Damage Detection and Assessment--continual or episodic updating of a structural model through vibration monitoring to detect and locate any structural damage. It can be argued that since the construction or modification of models using test data is subject to inherent uncertainties, the above problems should be properly treated within a Bayesian probabilistic framework. Such a methodology is presented which allows the precision of the estimates of the model parameters to be computed. It also leads to a guiding principle in applications. Namely, when selecting a single model from a given class of models, one should take the most probable model in the class based on the experimental data. Practical applications of this principle are given which are based on the utilization of measured seismic motions in large civil structures. Examples include the application of a computer program MODE-ID to identify modal properties directly from seismic excitation and response time histories from a nine-story steel-frame building at JPL and from a freeway overpass bridge.

Beck, James L.↗

High-Fidelity Generalization Method of Cells for Inelastic Periodic Multiphase Materials

An extension of a recently-developed linear thermoelastic theory for multiphase periodic materials is presented which admits inelastic behavior of the constituent phases. The extended theory is capable of accurately estimating both the effective inelastic response of a periodic multiphase composite and the local stress and strain fields in the individual phases. The model is presently limited to materials characterized by constituent phases that are continuous in one direction, but arbitrarily distributed within the repeating unit cell which characterizes the material's periodic microstructure. The model's analytical framework is based on the homogenization technique for periodic media, but the method of solution for the local displacement and stress fields borrows concepts previously employed by the authors in constructing the higher-order theory for functionally graded materials, in contrast with the standard finite-element solution method typically used in conjunction with the homogenization technique. The present approach produces a closed-form macroscopic constitutive equation for a periodic multiphase material valid for both uniaxial and multiaxial loading. The model's predictive accuracy in generating both the effective inelastic stress-strain response and the local stress said inelastic strain fields is demonstrated by comparison with the results of an analytical inelastic solution for the axisymmetric and axial shear response of a unidirectional composite based on the concentric cylinder model, and with finite-element results for transverse loading.

Aboudi, Jacob↗

Vibration-Based Data Used to Detect Cracks in Rotating Disks

Rotor health monitoring and online damage detection are increasingly gaining the interest of aircraft engine manufacturers. This is primarily due to the fact that there is a necessity for improved safety during operation as well as a need for lower maintenance costs. Applied techniques for the damage detection and health monitoring of rotors are essential for engine safety, reliability, and life prediction. Recently, the United States set the ambitious goal of reducing the fatal accident rate for commercial aviation by 80 percent within 10 years. In turn, NASA, in collaboration with the Federal Aviation Administration, other Federal agencies, universities, and the airline and aircraft industries, responded by developing the Aviation Safety Program. This program provides research and technology products needed to help the aerospace industry achieve their aviation safety goal. The Nondestructive Evaluation (NDE) Group of the Optical Instrumentation Technology Branch at the NASA Glenn Research Center is currently developing propulsion-system-specific technologies to detect damage prior to catastrophe under the propulsion health management task. Currently, the NDE group is assessing the feasibility of utilizing real-time vibration data to detect cracks in turbine disks. The data are obtained from radial blade-tip clearance and shaft-clearance measurements made using capacitive or eddy-current probes. The concept is based on the fact that disk cracks distort the strain field within the component. This, in turn, causes a small deformation in the disk's geometry as well as a possible change in the system's center of mass. The geometric change and the center of mass shift can be indirectly characterized by monitoring the amplitude and phase of the first harmonic (i.e., the 1 component) of the vibration data. Spin pit experiments and full-scale engine tests have been conducted while monitoring for crack growth with this detection methodology. Even so, published data are extremely limited, and the basic foundation of the methodology has not been fully studied. The NDE group is working on developing this foundation on the basis of theoretical modeling as well as experimental data by using the newly constructed subscale spin system shown in the preceding photograph. This, in turn, involved designing an optimal sub-scale disk that was meant to represent a full-scale turbine disk; conducting finite element analyses of undamaged and damaged disks to define the disk's deformation and the resulting shift in center of mass; and creating a rotordynamic model of the complete disk and shaft assembly to confirm operation beyond the first critical concerning the subscale experimental setup. The finite element analysis data, defining the center of mass shift due to disk damage, are shown. As an example, the change in the center of mass for a disk spinning at 8000 rpm with a 0.963-in. notch was 1.3 x 10(exp -4) in. The actual vibration response of an undamaged disk as well as the theoretical response of a cracked disk is shown. Experiments with cracked disks are continuing, and new approaches for analyzing the captured vibration data are being developed to better detect damage in a rotor. In addition, the subscale spin system is being used to test the durability and sensitivity of new NDE sensors that focus on detecting localized damage. This is designed to supplement the global response of the crack-detection methodology described here.

Gyekenyesi, Andrew L.↗

Dynamic Characterization of NASA Armstrong Test Facility Mechanical Vibration Facility

In preparation for Sierra Nevada Corporation’s (SNC) Dream Chaser spacecraft vibration test campaign at NASA Armstrong Test Facility (ATF) in Sandusky, Ohio, a dynamic characterization of the Mechanical Vibration Facility (MVF) is currently ongoing. The Mechanical Vibration Facility is comprised of an 18 ft diameter annulus table that is driven with sixteen hydraulic vertical actuator assemblies and four hydraulic horizontal actuator assemblies. During a test campaign this allows for single axis vibration testing in the vertical axis and in each of the two orthogonal horizontal axes without the need for reconfiguring the test article. The dynamic characterization of the facility was defined to follow a building-block approach requiring three vibration test configurations: (1) MVF bare table (2) MVF table and head expander (3) MVF table, head expander, and a dummy test article. Following the completion of the facility dynamic characterization, a facility finite element model (FEM) will be test-verified allowing for more accurate pretest analysis used for determining response limits and abort levels. To date two of the three test configurations have been completed. Test results have been used to verify the facility FEM that includes both the 18 ft diameter annulus table and the aluminum head expander that fills in the center opening of the annulus to provide a continuous flat mounting surface. Both the MVF table and the head expander models were previously correlated using test results from free-free boundary condition modal tests so model updates to the facility FEM focused on the uncorrelated vertical actuator assemblies and horizontal actuator assemblies. The accuracy of the test-verified facility FEM will be further determined based on comparison of analytical and test results from the vibration test of the dummy test article. Efforts will also be made to determine fixed base modes of the test article attached to the flexible vibration table using fixed base correction techniques. This paper will discuss test and analytical results from the MVF dynamic characterization test sequence and corresponding model updating effort.

Emma L Pierson↗

Interfacing modules for integrating discipline specific structural mechanics codes

An outline of the organization and capabilities of the Engine Structures Computational Simulator (Simulator) at NASA Lewis Research Center is given. One of the goals of the research at Lewis is to integrate various discipline specific structural mechanics codes into a software system which can be brought to bear effectively on a wide range of engineering problems. This system must possess the qualities of being effective and efficient while still remaining user friendly. The simulator was initially designed for the finite element simulation of gas jet engine components. Currently, the simulator has been restricted to only the analysis of high pressure turbine blades and the accompanying rotor assembly, although the current installation can be expanded for other applications. The simulator presently assists the user throughout its procedures by performing information management tasks, executing external support tasks, organizing analysis modules and executing these modules in the user defined order while maintaining processing continuity.

Endres, Ned M.↗