Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element model”

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,063 records · Page 59

A model of the wall boundary layer for ducted propellers

The objective of the present study is to include a representation of a wall boundary layer in an existing finite element model of the propeller in the wind tunnel environment. The major consideration is that the new formulation should introduce only modest alterations in the numerical model and should still be capable of producing economical predictions of the radiated acoustic field. This is accomplished by using a stepped approximation in which the velocity profile is piecewise constant in layers. In the limit of infinitesimally thin layers, the velocity profile of the stepped approximation coincides with that of the continuous profile. The approach described here could also be useful in modeling the boundary layer in other duct applications, particularly in the computation of the radiated acoustic field for sources contained in a duct.

Eversman, Walter↗

Combustion Dynamics and Stability Modeling of a Liquid Oxygen/RP-2 Oxygen-Rich Staged Combustion Preburner and Thrust Chamber Assembly with Gas-Centered Swirl Coaxial Injector Elements

The Combustion Stability Tool Development (CSTD) project, funded by the Air Force Space and Missile Systems Center, began in March 2015 supporting a renewed interest in the development of a liquid oxygen/hydrocarbon, oxygen-rich combustion engine. The project encompasses the design, assembly, and hot-fire testing of the NASA Marshall Space Flight Center 40-klbf Integrated Test Rig (MITR). The test rig models a staged-combustion configuration by combining an oxygen-rich preburner (ORPB), to generate hot gas, with a thrust chamber assembly (TCA) using gas-centered swirl coaxial injector elements. There are five separately designed interchangeable injectors in the TCA that each contain 19- or 27- injector elements. A companion paper in this JANNAF conference describes the design characteristics, rationale, and fabrication issues for all the injectors. The data acquired from a heavily instrumented rig encompasses several injectors, several operating points, and stability bomb tests. Another companion paper in this JANNAF conference describes this test program in detail. In this paper, dynamic data from the hot-fire testing is characterized and used to identify the responses in the ORPB and TCA. A brief review of damping metrics are discussed and applied as a measure of stability margin for damped acoustic modes. Chug and longitudinal combustion stability models and predictions are described which includes new dynamic models for compressible flow through an orifice and a modification to incorporate a third feed line for inclusion of the fuel-film coolant. Flow-acoustics finite element modeling is used to investigate the anticipated TCA acoustics, the effects of injector element length on stability margin, and the potential use of an ORPB orifice trip ring for improving longitudinal stability margin.

Casiano, M. J.↗

Multifidelity uncertainty quantification with models based on dissimilar parameters

Multifidelity uncertainty quantification (MF UQ) sampling approaches have been shown to significantly reduce the variance of statistical estimators while preserving the bias of the highest-fidelity model, provided that the low-fidelity models are well correlated. However, maintaining a high level of correlation can be challenging, especially when models depend on different input uncertain parameters, which drastically reduces the correlation. Existing MF UQ approaches do not adequately address this issue. In this work, we propose a new sampling strategy that exploits a shared space to improve the correlation among models with dissimilar parameterization. We achieve this by transforming the original coordinates onto an auxiliary manifold using the adaptive basis (AB) method (Tipireddy and Ghanem, 2014). The AB method has two main benefits: (1) it provides an effective tool to identify the low-dimensional manifold on which each model can be represented, and (2) it enables easy transformation of polynomial chaos representations from high- to low-dimensional spaces. This latter feature is used to identify a shared manifold among models without requiring additional evaluations. Here we present two algorithmic flavors of the new estimator to cover different analysis scenarios, including those with legacy and non-legacy high-fidelity (HF) data. We provide numerical results for analytical examples, a direct field acoustic test, and a finite element model of a nuclear fuel assembly. For all examples, we compare the proposed strategy against both single-fidelity and MF estimators based on the original model parameterization.

42 ENGINEERING↗

Development of Computer Models for the Assessment of Foreign Body Impact Events on Composite Structures

The objective of this project was to model the 5-3/4 inch pressure vessels used on the NASA RTOP program in an attempt to learn more about how impact damage forms and what are the residual effects of the resulting damage. A global-local finite element model was developed for the bottle and the states of stress in the bottles were determined down to the constituent level. The experimental data that was generated on the NASA RTOP program was not in a form that enabled the model developed under this grant to be correlated with the experimental data. As a result of this exercise it is recommended that an experimental program be designed using statistical design of experiment techniques to generate data that can be used to isolate the phenomenon that control the formation of impact damage. This data should include residual property determinations so that models for post impact structural integrity can be developed. It is also recommended that the global-local methodology be integrated directly into the finite element code. This will require considerable code development.

Bucinell, Ronald B.↗

Control systems design for large flexible space structures

Several approaches for the design of reduced-order linear-quadratic-Gaussian type controllers for large space structures were proposed and evaluated using a continuous model of a long free-free beam. Sufficient conditions were derived for the asymptotic stability with this type of controller. A finite-element model of a free-free-free-free square plate was obtained for use in control systems studies. A method was developed for optimal damping enhancement in large space structures.

Joshi, S. M.↗

Reduced basis technique for calculating sensitivity coefficients of nonlinear structural response

An efficient reduced basis technique is presented for calculating the sensitivity of nonlinear structural response to variations in the design variables. The structure is discretized by using two-field mixed finite element models. The vector of structural response and its sensitivity coefficients (derivatives with respect to design variables) are each expressed as a linear combination of a small number of basis (or global approximation) vectors. The Bubnov-Galerkin technique is then used to approximate each of the finite element equations governing the response and the sensitivity coefficients by a small number of algebraic equations in the amplitudes of these vectors. The path derivatives (derivatives of the response vector with respect to path parameters, e.g., load parameters) are used as basis vectors for approximating the response. A combination of the path derivatives and their derivatives with respect to the design variables is used for approximating the sensitivity coefficients. The potential of the proposed technique is discussed and its effectiveness is demonstrated by means of numerical examples of laminated composite plates subjected to mechanical and thermal loads.

Noor, Ahmed K.↗

Application of maximum entropy optimal projection design synthesis to the NASA Spacecraft Control Laboratory Experiment (SCOLE)

The scope of this study covered steady-state, continuous-time vibration control under disturbances applied to the Space Shuttle and continuous-time models of actuators, sensors, and disturbances. Focus was on a clear illustration of the methodology, therefore sensor/actuator dynamics were initially ignored, and a finite element model of the NASA Spacecraft Control Laboratory Experiment (SCOLE) was conducted, including products of inertia and offset of reflector CM from the mast tip.

Hyland, Dave↗

Modeling and Analysis of Structural Dynamics for a One-Tenth Scale Model NGST Sunshield

New modeling and analysis techniques have been developed for predicting the dynamic behavior of the Next Generation Space Telescope (NGST) sunshield. The sunshield consists of multiple layers of pretensioned, thin-film membranes supported by deployable booms. Modeling the structural dynamic behavior of the sunshield is a challenging aspect of the problem due to the effects of membrane wrinkling. A finite element model of the sunshield was developed using an approximate engineering approach, the cable network method, to account for membrane wrinkling effects. Ground testing of a one-tenth scale model of the NGST sunshield were carried out to provide data for validating the analytical model. A series of analyses were performed to predict the behavior of the sunshield under the ground test conditions. Modal analyses were performed to predict the frequencies and mode shapes of the test article and transient response analyses were completed to simulate impulse excitation tests. Comparison was made between analytical predictions and test measurements for the dynamic behavior of the sunshield. In general, the results show good agreement with the analytical model correctly predicting the approximate frequency and mode shapes for the significant structural modes.

Johnston, John↗

General purpose program to generate compatibility matrix for the integrated force method

An efficient procedure for obtaining the compatibility conditions of finite-element models involves the generation of both field and compatibility conditions from deformation-displacement relations, using (1) the compatibility bandwidth, and (2) the node-determinacy concept. A computer program thus structured will generate sparse and banded compatibility conditions for a structure that is idealized by the finite elements.

Nagabhusanam, J.↗

Structural Modeling of a Five-Meter Thin Film Inflatable Antenna/Concentrator With Rigidized Support Struts

Dynamic characterization of a non-rigidized thin film inflatable antenna/solar concentrator structure with rigidized composite support struts is described in detail. A two-step finite element modeling approach in MSC/NASTRAN is utilized, consisting of: (1) a nonlinear static pressurization procedure used to obtain the updated stiffness matrix, and (2) a modal "restart" eigen solution that uses the modified stiffness matrix. Unique problems encountered in modeling of this large 5-m lightweight inflatable are identified, including considerable difficulty in obtaining convergence in the nonlinear pressurization solution. It was found that the extremely thin polyimide film material (.001 in or I mil) presents tremendous problems in obtaining a converged solution when internal pressure loading is applied. It was concluded that the ratios of film thickness to other geometric dimensions such as torus cross-sectional and ring diameter and lenticular diameter are the critical parameters for convergence of the pressurization procedure. Comparison of finite element predictions for frequency and mode shapes with experimental results indicated reasonable agreement considering the complexity of the structure, the film-to-air interaction, and the nonlinear material properties of the film. It was also concluded that analysis should be done using different finite element to codes to determine if a more robust and stable solution can be obtained.

Smalley, Kurt B.↗

Predicting the Dynamic Crushing Response of a Composite Honeycomb Energy Absorber Using Solid-Element-Based Models in LS-DYNA

This paper describes an analytical study that was performed as part of the development of an externally deployable energy absorber (DEA) concept. The concept consists of a composite honeycomb structure that can be stowed until needed to provide energy attenuation during a crash event, much like an external airbag system. One goal of the DEA development project was to generate a robust and reliable Finite Element Model (FEM) of the DEA that could be used to accurately predict its crush response under dynamic loading. The results of dynamic crush tests of 50-, 104-, and 68-cell DEA components are presented, and compared with simulation results from a solid-element FEM. Simulations of the FEM were performed in LS-DYNA(Registered TradeMark) to compare the capabilities of three different material models: MAT 63 (crushable foam), MAT 26 (honeycomb), and MAT 126 (modified honeycomb). These material models are evaluated to determine if they can be used to accurately predict both the uniform crushing and final compaction phases of the DEA for normal and off-axis loading conditions

Jackson, Karen E.↗

Observation and Modeling of Dynamic Fracture Behaviors of Battery Cell Under Impact Loading Using Enhanced Representative Volume Element Concept

The burgeoning electric automobile industry has increased interest in battery safety. Battery cells experience significant mechanical stress during operation, including the impact of accidents and vibrations from driving. The potential for thermal runaway reactions in battery cells raises safety concerns. Although numerous researchers have defined the dynamic behavior of battery cells and proposed numerical models to describe it, few studies have focused on the high-strain rate mechanical impact phase correlated with the onset of fracture. In this study, we describe the dynamic behavior of pouch battery cells and propose a modeling method to study their mechanical failure under impact situations. Impact tests are conducted at various velocities and heights. To overcome numerical issues commonly encountered under rapid deformation scenarios, a new finite element model is developed based on the representative volume element model. The proposed approach efficiently simulates continuous crack propagation and brittleness behavior during impact by permitting the individual behavior of the cell components. Therefore, engineers can reliably design safer electric vehicle battery cells by measuring the properties of the cell components.

ENERGY STORAGE↗

Measurement of the Static Nonlinear Third-Order Elastic Moduli of Rocks: Problems and Applicability

The third-order elastic (TOE) model has been used to describe the widely observed nonlinear mechanical behaviors of earth materials. In addition to linear elastic constants ( λ , μ ), three nonlinear elastic moduli ( A , B , C ) are required for isotropic rocks. Contrary to previous research on dynamic TOE moduli, this study followed the protocol to measure strain and stress under uniaxial and hydrostatic compressive tests statically, which were later used to invert for the full set of TOE moduli for four standard rock types with differing pore structures of a porous oolitic limestone, quartz-rich sandstones, and a dense crystalline basalt. The applicability of the TOE model to characterize nonlinearity depends on the fulfillment of path-independence and small-strain assumptions. Using the measured static TOE moduli, the finite element model demonstrates that the stress in the vicinity of the wellbore is more amplified than the stress in the linear elastic case, which leads to a wider zone of rock failure around the wellbore. Due to the long-standing discrepancy between static and dynamic moduli, the rarely reported full set of static TOE moduli is necessary and will benefit future research in understanding the effect of rock nonlinearity on geophysical and geomechanical applications, such as long-term safe storage of CO 2 and generating process of geohazards.

58 GEOSCIENCES↗

A critical assessment of viscous models of trench topography and corner flow

Stresses for Newtonian viscous flow in a simple geometry (e.g., corner flow, bending flow) are obtained in order to study the effect of imposed velocity boundary conditions. Stress for a delta function velocity boundary condition decays as 1/R(2); for a step function velocity, stress goes as 1/R; for a discontinuity in curvature, the stress singularity is logarithmic. For corner flow, which has a discontinuity of velocity at a certain point, the corresponding stress has a 1/R singularity. However, for a more realistic circular-slab model, the stress singularity becomes logarithmic. Thus the stress distribution is very sensitive to the boundary conditions, and in evaluating the applicability of viscous models of trench topography it is essential to use realistic geometries. Topography and seismicity data from northern Hoshu, Japan, were used to construct a finite element model, with flow assumed tangent to the top of the grid, for both Newtonian and non-Newtonian flow (power law 3 rheology). Normal stresses at the top of the grid are compared to the observed trench topography and gravity anomalies. There is poor agreement. Purely viscous models of subducting slables with specified velocity boundary conditions do not predict normal stress patterns compatible with observed topography and gravity. Elasticity and plasticity appear to be important for the subduction process.

Zhang, J.↗

A critical assessment of viscous models of trench topography and corner flow

Stresses for Newtonian viscous flow in a simple geometry (e.g., corner flow, bending flow) are obtained in order to study the effect of imposed velocity boundary conditions. Stress for a delta function velocity boundary condition decays as 1/R(2); for a step function velocity, stress goes as 1/R; for a discontinuity in curvature, the stress singularity is logarithmic. For corner flow, which has a discontinuity of velocity at a certain point, the corresponding stress has a 1/R singularity. However, for a more realistic circular-slab model, the stress singularity becomes logarithmic. Thus the stress distribution is very sensitive to the boundary conditions, and in evaluating the applicability of viscous models of trench topography it is essential to use realistic geometries. Topography and seismicity data from northern Hoshu, Japan, were used to construct a finite element model, with flow assumed tangent to the top of the grid, for both Newtonian and non-Newtonian flow (power law 3 rheology). Normal stresses at the top of the grid are compared to the observed trench topography and gravity anomalies. There is poor agreement. Purely viscous models of subducting slables with specified velocity boundary conditions do not predict normal stress patterns compatible with observed topography and gravity. Elasticity and plasticity appear to be important for the subduction process.

Zhang, J.↗

Incorporation of a progressive failure analysis method in the CSM testbed software system

Analysis of the postbuckling behavior of composite shell structures pose many difficult and challenging problems in the field of structural mechanics. Current analysis methods perform well for most cases in predicting the postbuckling response of undamaged components. To predict component behavior accurately at higher load levels, the analysis must include the effects of local material failures. The CSM testbed software system is a highly modular structural analysis system currently under development at Langley Research Center. One of the primary goals of the CSM testbed is to provide a software environment for the development of advanced structural analysis methods and modern numerical methods which will exploit advanced computer architecture such as parallel-vector processors. Development of a progressive failure analysis method consists of the design and implementation of a processor which will perform the ply-level progressive failure analysis and the development of a geometrically nonlinear analysis procedure which incorporates the progressive failure processor. Regarding the development of the progressive failure processor, two components are required: failure criteria and a degradation model. For the initial implementation, the failure criteria of Hashin will be used. For a matrix failure which typically indicates the development of transverse matrix cracks, the ply properties will be degraded. Work to date includes the design of the progressive failure analysis processor and initial plans for the controlling geometrically nonlinear analysis procedure. The implementation of the progressive failure analysis has begun. Access to the model database and the Hashin failure criteria are completed. Work is in progress on the input/output operations for the processor related data and the finite element model updating procedures. In total the progressive failure processor is approximately one-third complete.

Arenburg, Robert T.↗

Time simulation of flutter with large stiffness changes

Time simulation of flutter, involving large local structural changes, is formulated with a state-space model that is based on a relatively small number of generalized coordinates. Free-free vibration modes are first calculated for a nominal finite-element model with relatively large fictitious masses located at the area of structural changes. A low-frequency subset of these modes is then transformed into a set of structural modal coordinates with which the entire simulation is performed. These generalized coordinates and the associated oscillatory aerodynamic force coefficient matrices are used to construct an efficient time-domain, state-space model for a basic aeroelastic case. The time simulation can then be performed by simply changing the mass, stiffness, and damping coupling terms when structural changes occur. It is shown that the size of the aeroelastic model required for time simulation with large structural changes at a few apriori known locations is similar to that required for direct analysis of a single structural case. The method is applied to the simulation of an aeroelastic wind-tunnel model. The diverging oscillations are followed by the activation of a tip-ballast decoupling mechanism that stabilizes the system but may cause significant transient overshoots.

Karpel, Mordechay↗

Time simulation of flutter with large stiffness changes

Time simulation of flutter, involving large local structural changes, is formulated with a state-space model that is based on a relatively small number of generalized coordinates. Free-free vibration modes are first calculated for a nominal finite-element model with relatively large fictitious masses located at the area of structural changes. A low-frequency subset of these modes is then transformed into a set of structural modal coordinates with which the entire simulation is performed. These generalized coordinates and the associated oscillatory aerodynamic force coefficient matrices are used to construct an efficient time-domain, state-space model for basic aeroelastic case. The time simulation can then be performed by simply changing the mass, stiffness and damping coupling terms when structural changes occur. It is shown that the size of the aeroelastic model required for time simulation with large structural changes at a few a priori known locations is similar to that required for direct analysis of a single structural case. The method is applied to the simulation of an aeroelastic wind-tunnel model. The diverging oscillations are followed by the activation of a tip-ballast decoupling mechanism that stabilizes the system but may cause significant transient overshoots.

Karpel, M.↗