Engineering PapersSearch

SEARCH · Engineering Papers

Results for “FEM”

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 91 records · Page 5

Comparison of boundary element and finite element methods in spur gear root stress analysis

The boundary element method (BEM) is used to compute fillet stress concentration in spur gear teeth. The results are shown to compare favorably with analogous results obtained using the finite element method (FEM). A partially supported thin rim gear is studied. The loading is applied at the pitch point. A three-dimensional analysis is conducted using both the BEM and FEM (NASTRAN). The results are also compared with those of a two-dimensional finite element model. An advantage of the BEM over the FEM is that fewer elements are needed with the BEM. Indeed, in the current study the BEM used 92 elements and 270 nodes whereas the FEM used 320 elements and 2037 nodes. Moreover, since the BEM is especially useful in problems with high stress gradients it is potentially a very useful tool for fillet stress analyses.

Sun, H.

Validation of finite element and boundary element methods for predicting structural vibration and radiated noise

Analytical and experimental validation of methods to predict structural vibration and radiated noise are presented. A rectangular box excited by a mechanical shaker was used as a vibrating structure. Combined finite element method (FEM) and boundary element method (BEM) models of the apparatus were used to predict the noise radiated from the box. The FEM was used to predict the vibration, and the surface vibration was used as input to the BEM to predict the sound intensity and sound power. Vibration predicted by the FEM model was validated by experimental modal analysis. Noise predicted by the BEM was validated by sound intensity measurements. Three types of results are presented for the total radiated sound power: (1) sound power predicted by the BEM modeling using vibration data measured on the surface of the box; (2) sound power predicted by the FEM/BEM model; and (3) sound power measured by a sound intensity scan. The sound power predicted from the BEM model using measured vibration data yields an excellent prediction of radiated noise. The sound power predicted by the combined FEM/BEM model also gives a good prediction of radiated noise except for a shift of the natural frequencies that are due to limitations in the FEM model.

Seybert, A. F.

Modal Test/Analysis Correlation of Space Station Structures Using Nonlinear Sensitivity

The modal correlation problem is formulated as a constrained optimization problem for validation of finite element models (FEM's). For large-scale structural applications, a pragmatic procedure for substructuring, model verification, and system integration is described to achieve effective modal correlation. The space station substructure FEM's are reduced using Lanczos vectors and integrated into a system FEM using Craig-Bampton component modal synthesis. The optimization code is interfaced with MSC/NASTRAN to solve the problem of modal test/analysis correlation; that is, the problem of validating FEM's for launch and on-orbit coupled loads analysis against experimentally observed frequencies and mode shapes. An iterative perturbation algorithm is derived and implemented to update nonlinear sensitivity (derivatives of eigenvalues and eigenvectors) during optimizer iterations, which reduced the number of finite element analyses.

Gupta, Viney K.

Modeling structural dynamic behavior of SSME components

FEM studies are presented of the nozzle and the low-pressure fuel-pump inducer designs for the Space Shuttle Main Engine (SSME) to analyze the effects of structural vibrations. FEM preprocessing software based on a CAD system is employed to develop a model of the component's sophisticated geometry. The nozzle geometry is also defined by means of the preprocessing technique and subsequently analyzed with respect to time-transient loading. The analysis is conducted with a Cray supercomputer using the SPAR/EAL FEM program. The investigation of the nozzle demonstrates the advantageous use of symmetry in the determination of nozzle response to SSME start-up transients. Plots of time vs strain are developed for gages on the nozzle wall and steerhorn tubing. The results of the inducer modeling are found to be adequate for investigating the component's principle modes, and the nozzle results indicate the suitability of the FEM techniques for optimizing the design of engine components.

Kiefling, Larry A.

Modal test/analysis correlation of Space Station structures using nonlinear sensitivity

The modal correlation problem is formulated as a constrained optimization problem for validation of finite element models (FEM's). For large-scale structural applications, a pragmatic procedure for substructuring, model verification, and system integration is described to achieve effective modal correlations. The space station substructure FEM's are reduced using Lanczos vectors and integrated into a system FEM using Craig-Bampton component modal synthesis. The optimization code is interfaced with MSC/NASTRAN to solve the problem of modal test/analysis correlation; that is, the problem of validating FEM's for launch and on-orbit coupled loads analysis against experimentally observed frequencies and mode shapes. An iterative perturbation algorithm is derived and implemented to update nonlinear sensitivity (derivatives of eigenvalues and eigenvectors) during optimizer iterations, which reduced the number of finite element analyses.

Gupta, Viney K.

Experimental validation of finite element and boundary element methods for predicting structural vibration and radiated noise

This research report is presented in three parts. In the first part, acoustical analyses were performed on modes of vibration of the housing of a transmission of a gear test rig developed by NASA. The modes of vibration of the transmission housing were measured using experimental modal analysis. The boundary element method (BEM) was used to calculate the sound pressure and sound intensity on the surface of the housing and the radiation efficiency of each mode. The radiation efficiency of each of the transmission housing modes was then compared to theoretical results for a finite baffled plate. In the second part, analytical and experimental validation of methods to predict structural vibration and radiated noise are presented. A rectangular box excited by a mechanical shaker was used as a vibrating structure. Combined finite element method (FEM) and boundary element method (BEM) models of the apparatus were used to predict the noise level radiated from the box. The FEM was used to predict the vibration, while the BEM was used to predict the sound intensity and total radiated sound power using surface vibration as the input data. Vibration predicted by the FEM model was validated by experimental modal analysis; noise predicted by the BEM was validated by measurements of sound intensity. Three types of results are presented for the total radiated sound power: sound power predicted by the BEM model using vibration data measured on the surface of the box; sound power predicted by the FEM/BEM model; and sound power measured by an acoustic intensity scan. In the third part, the structure used in part two was modified. A rib was attached to the top plate of the structure. The FEM and BEM were then used to predict structural vibration and radiated noise respectively. The predicted vibration and radiated noise were then validated through experimentation.

Seybert, A. F.

Finite element method for eigenvalue problems in electromagnetics

Finite element method (FEM) has been a very powerful tool to solve many complex problems in electromagnetics. The goal of the current research at the Langley Research Center is to develop a combined FEM/method of moments approach to three-dimensional scattering/radiation problem for objects with arbitrary shape and filled with complex materials. As a first step toward that goal, an exercise is taken to establish the power of FEM, through closed boundary problems. This paper demonstrates the developed of FEM tools for two- and three-dimensional eigenvalue problems in electromagnetics. In section 2, both the scalar and vector finite elements have been used for various waveguide problems to demonstrate the flexibility of FEM. In section 3, vector finite element method has been extended to three-dimensional eigenvalue problems.

Reddy, C. J.

Space station static and dynamic analyses using parallel methods

Algorithms for high-performance parallel computers are applied to perform static analyses of large-scale Space Station finite-element models (FEMs). Several parallel-vector algorithms under development at NASA Langley are assessed. Sparse matrix solvers were found to be more efficient than banded symmetric or iterative solvers for the static analysis of large-scale applications. In addition, new sparse and 'out-of-core' solvers were found superior to substructure (superelement) techniques which require significant additional cost and time to perform static condensation during global FEM matrix generation as well as the subsequent recovery and expansion. A method to extend the fast parallel static solution techniques to reduce the computation time for dynamic analysis is also described. The resulting static and dynamic algorithms offer design economy for preliminary multidisciplinary design optimization and FEM validation against test modes. The algorithms are being optimized for parallel computers to solve one-million degrees-of-freedom (DOF) FEMs. The high-performance computers at NASA afforded effective software development, testing, efficient and accurate solution with timely system response and graphical interpretation of results rarely found in industry. Based on the author's experience, similar cooperation between industry and government should be encouraged for similar large-scale projects in the future.

Gupta, V.

Analysis of Waveguide Junction Discontinuities Using Finite Element Method

A Finite Element Method (FEM) is presented to determine reflection and transmission coefficients of rectangular waveguide junction discontinuities. An H-plane discontinuity, an E-plane ridge discontinuity, and a step discontinuity in a concentric rectangular waveguide junction are analyzed using the FEM procedure. Also, reflection and transmission coefficients due to presence of a gap between two sections of a rectangular waveguide are determined using the FEM. The numerical results obtained by the present method are in excellent agreement with the earlier published results. The numerical results obtained by the FEM are compared with the numerical results obtained using the Mode Matching Method (MMM) and also with the measured data.

Deshpande, Manohar D.

True Concurrent Thermal Engineering Integrating CAD Model Building with Finite Element and Finite Difference Methods

Thermal engineering has long been left out of the concurrent engineering environment dominated by CAD (computer aided design) and FEM (finite element method) software. Current tools attempt to force the thermal design process into an environment primarily created to support structural analysis, which results in inappropriate thermal models. As a result, many thermal engineers either build models "by hand" or use geometric user interfaces that are separate from and have little useful connection, if any, to CAD and FEM systems. This paper describes the development of a new thermal design environment called the Thermal Desktop. This system, while fully integrated into a neutral, low cost CAD system, and which utilizes both FEM and FD methods, does not compromise the needs of the thermal engineer. Rather, the features needed for concurrent thermal analysis are specifically addressed by combining traditional parametric surface based radiation and FD based conduction modeling with CAD and FEM methods. The use of flexible and familiar temperature solvers such as SINDA/FLUINT (Systems Improved Numerical Differencing Analyzer/Fluid Integrator) is retained.

Panczak, Tim

An Investigation of Two Acoustic Propagation Codes for Three-Dimensional Geometries

The ability to predict fan noise within complex three-dimensional aircraft engine nacelle geometries is a valuable tool in studying low-noise designs. Recent years have seen the development of aeroacoustic propagation codes using various levels of approximation to obtain such a capability. In light of this, it is beneficial to pursue a design paradigm that incorporates the strengths of the various tools. The development of a quasi-3D methodology (Q3D-FEM) at NASA Langley has brought these ideas to mind in relation to the framework of the CDUCT-LaRC acoustic propagation and radiation tool. As more extensive three dimensional codes become available, it would seem appropriate to incorporate these tools into a framework similar to CDUCT-LaRC and use them in a complementary manner. This work focuses on such an approach in beginning the steps toward a systematic assessment of the errors, and hence the trade-offs, involved in the use of these codes. To illustrate this point, CDUCT-LaRC was used to study benchmark hardwall duct problems to quantify errors caused by wave propagation in directions far removed from that defined by the parabolic approximation. Configurations incorporating acoustic treatment were also studied with CDUCT-LaRC and Q3D-FEM. The cases presented show that acoustic treatment diminishes the effects of CDUCT-LaRC phase error as the solutions are attenuated. The results of the Q3D-FEM were very promising and matched the analytic solution very well. Overall, these tests were meant to serve as a step toward the systematic study of errors inherent in the propagation module of CDUCT-LaRC, as well as an initial test of the higher fidelity Q3D-FEM code.

Nark, D. M.

A Statistical Approach for the Concurrent Coupling of Molecular Dynamics and Finite Element Methods

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, increasing the size of the MD domain quickly presents intractable computational demands. A robust approach to surmount this computational limitation has been to unite continuum modeling procedures such as the finite element method (FEM) with MD analyses thereby reducing the region of atomic scale refinement. The challenging problem is to seamlessly connect the two inherently different simulation techniques at their interface. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the typical boundary value problem used to define a coupled domain. The method uses statistical averaging of the atomistic MD domain to provide displacement interface boundary conditions to the surrounding continuum FEM region, which, in return, generates interface reaction forces applied as piecewise constant traction boundary conditions to the MD domain. The two systems are computationally disconnected and communicate only through a continuous update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM) as opposed to a direct coupling method where interface atoms and FEM nodes are individually related. The methodology is inherently applicable to three-dimensional domains, avoids discretization of the continuum model down to atomic scales, and permits arbitrary temperatures to be applied.

Saether, E.

A New Concurrent Multiscale Methodology for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Yamakov, Vesselin

An Embedded Statistical Method for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Saether, E.

Elastic Model Transitions: a Hybrid Approach Utilizing Quadratic Inequality Constrained Least Squares (LSQI) and Direct Shape Mapping (DSM)

A method for transitioning linear time invariant (LTI) models in time varying simulation is proposed that utilizes both quadratically constrained least squares (LSQI) and Direct Shape Mapping (DSM) algorithms to determine physical displacements. This approach is applicable to the simulation of the elastic behavior of launch vehicles and other structures that utilize multiple LTI finite element model (FEM) derived mode sets that are propagated throughout time. The time invariant nature of the elastic data for discrete segments of the launch vehicle trajectory presents a problem of how to properly transition between models while preserving motion across the transition. In addition, energy may vary between flex models when using a truncated mode set. The LSQI-DSM algorithm can accommodate significant changes in energy between FEM models and carries elastic motion across FEM model transitions. Compared with previous approaches, the LSQI-DSM algorithm shows improvements ranging from a significant reduction to a complete removal of transients across FEM model transitions as well as maintaining elastic motion from the prior state.

Jurenko, Robert J.

Aeroelastic Airworthiness Assesment of the Adaptive Compliant Trailing Edge Flaps

The Adaptive Compliant Trailing Edge (ACTE) demonstrator is a joint task under the National Aeronautics and Space Administration Environmentally Responsible Aviation Project in partnership with the Air Force Research Laboratory and FlexSys, Inc. (Ann Arbor, Michigan). The project goal is to develop advanced technologies that enable environmentally friendly aircraft, such as adaptive compliant technologies. The ACTE demonstrator flight-test program encompassed replacing the Fowler flaps on the SubsoniC Aircraft Testbed, a modified Gulfstream III (Gulfstream Aerospace, Savannah, Georgia) aircraft, with control surfaces developed by FlexSys. The control surfaces developed by FlexSys are a pair of uniquely-designed unconventional flaps to be used as lifting surfaces during flight-testing to validate their structural effectiveness. The unconventional flaps required a multidisciplinary airworthiness assessment to prove they could withstand the prescribed flight envelope. Several challenges were posed due to the large deflections experienced by the structure, requiring non-linear analysis methods. The aeroelastic assessment necessitated both conventional and extensive testing and analysis methods. A series of ground vibration tests (GVTs) were conducted to provide modal characteristics to validate and update finite element models (FEMs) used for the flutter analyses for a subset of the various flight configurations. Numerous FEMs were developed using data from FlexSys and the ground tests. The flap FEMs were then attached to the aircraft model to generate a combined FEM that could be analyzed for aeroelastic instabilities. The aeroelastic analysis results showed the combined system of aircraft and flaps were predicted to have the required flutter margin to successfully demonstrate the adaptive compliant technology. This paper documents the details of the aeroelastic airworthiness assessment described, including the ground testing and analyses, and subsequent flight-testing performed on the unconventional ACTE flaps.

aeroelasticity

Stress Intensity Factors for Layered Pressure Vessel Inner Layer Through Cracks

A difficulty encountered when performing Fitness-for- Service assessments for layered pressure vessels (LPVs) is the lack of stress intensity factor solution in literature that produce accurate results for inner layer longitudinal through cracks. Using surrogate solutions such as a through crack in a plate or cylinder produce results that can be overly conservative especially for longer cracks. This is largely due to the ability of a layered pressure vessel to redistribute hoop load to other layers, the restricted radial movement of the cracked layer, and the friction forces applied in the cracked region. To understand this problem, a parametric finite element model (FEM) generator was developed that is capable of producing layered pressure vessel models with inner layer through cracks. The results from the FEMs were used to create a dataset of inner layer through crack stress intensity factors (KI) for layered pressure vessels corresponding to variations of internal pressure, radius, layer thicknesses, friction factor, and crack length. The elastic modulus of the material also has an effect on KI but, for this dataset, the elastic modulus was fixed at the typical value for steel - 29,500 ksi (203 GPa). Finally, a non-dimensional model was developed and calibrated using the dataset. This allows KI to be calculated without the need of a FEM using a closed-form equation. The results of the closed-form solution were then compared to FEM results showing accuracy was generally within 10%.

Hobbs, Joel R.

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