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 433 records · Page 24

Development of a Continuum Damage Mechanics Material Model of a Graphite-Kevlar(Registered Trademark) Hybrid Fabric for Simulating the Impact Response of Energy Absorbing Kevlar(Registered Trademark) Hybrid Fabric for Simulating the Impact Response of Energy Absorbing

This paper describes the development of input properties for a continuum damage mechanics based material model, Mat 58, within LS-DYNA(Registered Trademark) to simulate the response of a graphite-Kevlar(Registered Trademark) hybrid plain weave fabric. A limited set of material characterization tests were performed on the hybrid graphite-Kevlar(Registered Trademark) fabric. Simple finite element models were executed in LS-DYNA(Registered Trademark) to simulate the material characterization tests and to verify the Mat 58 material model. Once verified, the Mat 58 model was used in finite element models of two composite energy absorbers: a conical-shaped design, designated the "conusoid," fabricated of four layers of hybrid graphite-Kevlar(Registered Trademark) fabric; and, a sinusoidal-shaped foam sandwich design, designated the "sinusoid," fabricated of the same hybrid fabric face sheets with a foam core. Dynamic crush tests were performed on components of the two energy absorbers, which were designed to limit average vertical accelerations to 25- to 40-g, to minimize peak crush loads, and to generate relatively long crush stroke values under dynamic loading conditions. Finite element models of the two energy absorbers utilized the Mat 58 model that had been verified through material characterization testing. Excellent predictions of the dynamic crushing response were obtained.

Jackson, Karen E.↗

Comparison of Structural Concepts for Transport Aircraft with a Tail Cone Turbine

An evaluation of three structural concepts for an advanced aircraft design with a tail cone turbine is presented. Structural models were developed using an innovative rapid finite element modeling tool called Conceptual Design Shop (CDS). CDS is an attempt to fill a gap in current finite element modeling software to automatically connect wings and tails to the fuselage with sufficient support structure in airframe models. For comparison with actual aircraft structures, a model of a transport aircraft design developed using CDS is compared with published structural weight data for a Boeing 737-200. A method for computing aerodynamic stability parameters from an MSC/NASTRAN finite element analysis output deck is also discussed. Finally, the weight effects for using a tail cone turbine in an advanced transport aircraft design are evaluated by comparing three composite structural models, for which component thicknesses and cross-sections are sized by the Hypersizer software.

Mason, Brian H.↗

Nonlinear transient analysis by energy minimization: A theoretical basis for the ACTION computer code

The formulation basis for establishing the static or dynamic equilibrium configurations of finite element models of structures which may behave in the nonlinear range are provided. With both geometric and time independent material nonlinearities included, the development is restricted to simple one and two dimensional finite elements which are regarded as being the basic elements for modeling full aircraft-like structures under crash conditions. Representations of a rigid link and an impenetrable contact plane are added to the deformation model so that any number of nodes of the finite element model may be connected by a rigid link or may contact the plane. Equilibrium configurations are derived as the stationary conditions of a potential function of the generalized nodal variables of the model. Minimization of the nonlinear potential function is achieved by using the best current variable metric update formula for use in unconstrained minimization. Powell's conjugate gradient algorithm, which offers very low storage requirements at some slight increase in the total number of calculations, is the other alternative algorithm to be used for extremely large scale problems.

Kamat, M. P.↗

Structural Model Tuning Tool: User's Reference Manual

This report presents an efficient approach for tuning finite element models to match the measured ground vibration test and static test data. Frequencies, mode shapes, total weight, location of the center of gravity, and static deformation computed from the finite element model are matched to the measured data. The model tuning procedure used in this work is based on solving an optimization problem in which the errors for the considered metrics, between the finite element prediction and the measured data, are minimized. Analytical sensitivity values of performance indices are computed using the NASTRAN-generated sensitivity values together with the in-house computer codes, which allow for faster computational time and the use of gradient-based optimizers. The method is applied to the Aerostructures Test Wing IV model. The study shows that the military standard and the National Aeronautics and Space Administration standard for comparing analytical and experimental modal data are all satisfied. The final finite element model correlates well with the test data. The flutter speed decreases by 8.91 percent after model tuning compared with the original Aerostructures Test Wing IV design.

Chan-gi Pak↗

A finite element propagation model for extracting normal incidence impedance in nonprogressive acoustic wave fields

A propagation model method for extracting the normal incidence impedance of an acoustic material installed as a finite length segment in a wall of a duct carrying a nonprogressive wave field is presented. The method recasts the determination of the unknown impedance as the minimization of the normalized wall pressure error function. A finite element propagation model is combined with a coarse/fine grid impedance plane search technique to extract the impedance of the material. Results are presented for three different materials for which the impedance is known. For each material, the input data required for the prediction scheme was computed from modal theory and then contaminated by random error. The finite element method reproduces the known impedance of each material almost exactly for random errors typical of those found in many measurement environments. Thus, the method developed here provides a means for determining the impedance of materials in a nonprogressirve wave environment such as that usually encountered in a commercial aircraft engine and most laboratory settings.

Watson, Willie R.↗

Vascular mechanics of the coronary artery

This paper describes our research into the vascular mechanics of the coronary artery and plaque. The three sections describe the determination of arterial mechanical properties using intravascular ultrasound (IVUS), a constitutive relation for the arterial wall, and finite element method (FEM) models of the arterial wall and atheroma. METHODS: Inflation testing of porcine left anterior descending coronary arteries was conducted. The changes in the vessel geometry were monitored using IVUS, and intracoronary pressure was recorded using a pressure transducer. The creep and quasistatic stress/strain responses were determined. A Standard Linear Solid (SLS) was modified to reproduce the non-linear elastic behavior of the arterial wall. This Standard Non-linear Solid (SNS) was implemented into an axisymetric thick-walled cylinder numerical model. Finite element analysis models were created for five age groups and four levels of stenosis using the Pathobiological Determinants of Atherosclerosis Youth (PDAY) database. RESULTS: The arteries exhibited non-linear elastic behavior. The total tissue creep strain was epsilon creep = 0.082 +/- 0.018 mm/mm. The numerical model could reproduce both the non-linearity of the porcine data and time dependent behavior of the arterial wall found in the literature with a correlation coefficient of 0.985. Increasing age had a strong positive correlation with the shoulder stress level, (r = 0.95). The 30% stenosis had the highest shoulder stress due to the combination of a fully formed lipid pool and a thin cap. CONCLUSIONS: Studying the solid mechanics of the arterial wall and the atheroma provide important insights into the mechanisms involved in plaque rupture.

NASA Program Biomedical Research and Countermeasur↗

Calibration of Airframe and Occupant Models for Two Full-Scale Rotorcraft Crash Tests

Two full-scale crash tests of an MD-500 helicopter were conducted in 2009 and 2010 at NASA Langley's Landing and Impact Research Facility in support of NASA s Subsonic Rotary Wing Crashworthiness Project. The first crash test was conducted to evaluate the performance of an externally mounted composite deployable energy absorber under combined impact conditions. In the second crash test, the energy absorber was removed to establish baseline loads that are regarded as severe but survivable. Accelerations and kinematic data collected from the crash tests were compared to a system integrated finite element model of the test article. Results from 19 accelerometers placed throughout the airframe were compared to finite element model responses. The model developed for the purposes of predicting acceleration responses from the first crash test was inadequate when evaluating more severe conditions seen in the second crash test. A newly developed model calibration approach that includes uncertainty estimation, parameter sensitivity, impact shape orthogonality, and numerical optimization was used to calibrate model results for the second full-scale crash test. This combination of heuristic and quantitative methods was used to identify modeling deficiencies, evaluate parameter importance, and propose required model changes. It is shown that the multi-dimensional calibration techniques presented here are particularly effective in identifying model adequacy. Acceleration results for the calibrated model were compared to test results and the original model results. There was a noticeable improvement in the pilot and co-pilot region, a slight improvement in the occupant model response, and an over-stiffening effect in the passenger region. This approach should be adopted early on, in combination with the building-block approaches that are customarily used, for model development and test planning guidance. Complete crash simulations with validated finite element models can be used to satisfy crash certification requirements, thereby reducing overall development costs.

Annett, Martin S.↗

Application of Single Crystal Failure Criteria: Theory and Turbine Blade Case Study

The orientation of the single crystal material within a structural component is known to affect the strength and life of the part. The first stage blade of the High Pressure Fuel Turbopump (HPFTP)/ Alternative Turbopump Development (ATD), of the Space Shuttle Main Engine (SSME) was used to study the effects of secondary axis'orientation angles on the failure rate of the blade. A new failure criterion was developed based on normal and shear strains on the primary crystallographic planes. The criterion was verified using low cycle fatigue (LCF) specimen data and a finite element model of the test specimens. The criterion was then used to study ATD/HPFTP first stage blade failure events. A detailed ANSYS finite element model of the blade was used to calculate the failure parameter for the different crystallographic orientations. A total of 297 cases were run to cover a wide range of acceptable orientations within the blade. Those orientations are related to the base crystallographic coordinate system that was created in the ANSYS finite element model. Contour plots of the criterion as a function of orientation for the blade tip and attachment were obtained. Results of the analysis revealed a 40% increase in the failure parameter due to changing of the primary and secondary axes of material orientations. A comparison between failure criterion predictions and actual engine test data was then conducted. The engine test data comes from two ATD/HPFTP builds (units F3- 4B and F6-5D), which were ground tested on the SSME at the Stennis Space Center in Mississippi. Both units experienced cracking of the airfoil tips in multiple blades, but only a few cracks grew all the way across the wall of the hollow core airfoil.

Sayyah, Tarek↗

Controls for space structures

Assembly and operation of large space structures (LSS) in orbit will require robot-assisted docking and berthing of partially-assembled structures. These operations require new solutions to the problems of controls. This is true because of large transient and persistent disturbances, controller-structure interaction with unmodeled modes, poorly known structure parameters, slow actuator/sensor dynamical behavior, and excitation of nonlinear structure vibrations during control and assembly. For on-orbit assembly, controllers must start with finite element models of LSS and adapt on line to the best operating points, without compromising stability. This is not easy to do, since there are often unmodeled dynamic interactions between the controller and the structure. The indirect adaptive controllers are based on parameter estimation. Due to the large number of modes in LSS, this approach leads to very high-order control schemes with consequent poor stability and performance. In contrast, direct model reference adaptive controllers operate to force the LSS to track the desirable behavior of a chosen model. These schemes produce simple control algorithms which are easy to implement on line. One problem with their use for LSS has been that the model must be the same dimension as the LSS - i.e., quite large. A control theory based on the command generator tracker (CGT) ideas of Sobel, Mabins, Kaufman and Wen, Balas to obtain very low-order models based on adaptive algorithms was developed. Closed-loop stability for both finite element models and distributed parameter models of LSS was proved. In addition, successful numerical simulations on several LSS databases were obtained. An adaptive controller based on our theory was also implemented on a flexible robotic manipulator at Martin Marietta Astronautics. Computation schemes for controller-structure interaction with unmodeled modes, the residual mode filters or RMF, were developed. The RMF theory was modified to compensate slow actuator/sensor dynamics. These new ideas are being applied to LSS simulations to demonstrate the ease with which one can incorporate slow actuator/sensor effects into our design. It was also shown that residual mode filter compensation can be modified for small nonlinearities to produce exponentially stable closed-loop control. A theory for disturbance accommodating controllers based on reduced order models of structures was developed, and stability results for these controllers in closed-loop with large-scale finite element models of structures were obtained.

Balas, Mark↗

Modes of interconnected lattice trusses using continuum models, part 1

This represents a continuing systematic attempt to explore the use of continuum models--in contrast to the Finite Element Models currently universally in use--to develop feedback control laws for stability enhancement of structures, particularly large structures, for deployment in space. We shall show that for the control objective, continuum models do offer unique advantages. It must be admitted of course that developing continuum models for arbitrary structures is no easy task. In this paper we take advantage of the special nature of current Large Space Structures--typified by the NASA-LaRC Evolutionary Model which will be our main concern--which consists of interconnected orthogonal lattice trusses each with identical bays. Using an equivalent one-dimensional Timoshenko beam model, we develop an almost complete continuum model for the evolutionary structure. We do this in stages, beginning only with the main bus as flexible and then going on to make all the appendages also flexible-except for the antenna structure. Based on these models we proceed to develop formulas for mode frequencies and shapes. These are shown to be the roots of the determinant of a matrix of small dimension compared with mode calculations using Finite Element Models, even though the matrix involves transcendental functions. The formulas allow us to study asymptotic properties of the modes and how they evolve as we increase the number of bodies which are treated as flexible. The asymptotics, in fact, become simpler.

Balakrishnan, A. V.↗

Adaptive Multi-Layer LMS Controller Design and Application to Active Vibration Suppression on a Truss and Proposed Impact Analysis Technique

This research develops an adaptive controller that actively suppresses a single frequency disturbance source at a remote position and tests the system on the NPS Space Truss. The experimental results are then compared to those predicted by an ANSYS finite element model. The NPS space truss is a 3.7-meter long truss that simulates a space-borne appendage with sensitive equipment mounted at its extremities. One of two installed piezoelectric actuators and an Adaptive Multi-Layer LMS control law were used to effectively eliminate an axial component of the vibrations induced by a linear proof mass actuator mounted at one end of the truss. Experimental and analytical results both demonstrate reductions to the level of system noise. Vibration reductions in excess of 50dB were obtained through experimentation and over 100dB using ANSYS, demonstrating the ability to model this system with a finite element model. This report also proposes a method to use distributed quartz accelerometers to evaluate the location, direction, and energy of impacts on the NPS space truss using the dSPACE data acquisition and processing system to capture the structural response and compare it to known reference Signals.

Barney, Timothy A.↗

Calculation of Thermally-Induced Displacements in Spherically Domed Ion Engine Grids

An analytical method for predicting the thermally-induced normal and tangential displacements of spherically domed ion optics grids under an axisymmetric thermal loading is presented. A fixed edge support that could be thermally expanded is used for this analysis. Equations for the displacements both normal and tangential to the surface of the spherical shell are derived. A simplified equation for the displacement at the center of the spherical dome is also derived. The effects of plate perforation on displacements and stresses are determined by modeling the perforated plate as an equivalent solid plate with modified, or effective, material properties. Analytical model results are compared to the results from a finite element model. For the solid shell, comparisons showed that the analytical model produces results that closely match the finite element model results. The simplified equation for the normal displacement of the spherical dome center is also found to accurately predict this displacement. For the perforated shells, the analytical solution and simplified equation produce accurate results for materials with low thermal expansion coefficients.

Soulas, George C.↗

Uncertainty Quantification Framework for Predicting Material Response with Large Number of Parameters: Application to Creep Prediction in Ferritic-Martensitic Steels Using Combined Crystal Plasticity and Grain Boundary Models

This paper presents an uncertainty quantification (UQ) framework for the physics-based model prediction of material response with a large number of parameters. The application problem presented in this work is that of predicting creep in Grade 91 steel at 600°C. The material response is defined with a physically based microstructural model with constitutive equations emulating several observed phenomena in Grade 91 and embodied into an explicit geometry mesoscale finite element model for prior austenite grains and grain boundaries. Creep within the grains and in grain boundaries are represented by crystal plasticity for dislocation motion and a physics-based model for cavity growth and nucleation, respectively. The creep behavior of this material is influenced by several parameters, some of which have a wide range of variation based on experimental data. UQ combined with microstructural modeling can discover the core microstructural causes of experimental variability, leading to improved materials with lower variability in critical long-term material properties. In this study, we investigate the model's uncertainty to identify material properties that may be modified during production to increase creep life and analyze different components of the crystal plasticity model for improvements. For this purpose, a quantity of interest is defined as time to minimum creep rate, which correlates well to the creep failure of the material. A deep neural network model was trained and validated to be used as a surrogate for the finite element model. Then, a variance-based sensitivity analysis is performed on the surrogate model to find the Sobol indices of the input parameters in respect to the output quantity of interest. The Sobol indices are used to reduce the dimensionality of the model. Generalized polynomial chaos expansion is used on the reduced basis models to propagate the uncertainty from the input parameters to the quantity of interest using the deep neural network surrogate model. These results are benchmarked against uncertainty propagation using Monte Carlo simulations. In conclusion, the UQ performed through the reduced basis model captures almost all the uncertainty in the model with significantly fewer simulations, making it possible to perform the UQ directly via simulations with the finite element model rather than surrogate machine-learned models.

36 MATERIALS SCIENCE↗

Hierarchic plate and shell models based on p-extension

Formulations of finite element models for beams, arches, plates and shells based on the principle of virtual work was studied. The focus is on computer implementation of hierarchic sequences of finite element models suitable for numerical solution of a large variety of practical problems which may concurrently contain thin and thick plates and shells, stiffeners, and regions where three dimensional representation is required. The approximate solutions corresponding to the hierarchic sequence of models converge to the exact solution of the fully three dimensional model. The stopping criterion is based on: (1) estimation of the relative error in energy norm; (2) equilibrium tests, and (3) observation of the convergence of quantities of interest.

Szabo, Barna A.↗

Vibration Propagation of Gear Dynamics in a Gear-Bearing-Housing System Using Mathematical Modeling and Finite Element Analysis

Vibration and noise caused by gear dynamics at the meshing teeth propagate through power transmission components to the surrounding environment. This study is devoted to developing computational tools to investigate the vibro-acoustic propagation of gear dynamics through a gearbox using different bearings. Detailed finite element/contact mechanics and boundary element models of the gear/bearing/housing system are established to compute the system vibration and noise propagation. Both vibration and acoustic models are validated by experiments including the vibration modal testing and sound field measurements. The effectiveness of each bearing type to disrupt vibration propagation is speed-dependent. Housing plays an important role in noise radiation .It, however, has limited effects on gear dynamics. Bearings are critical components in drivetrains. Accurate modeling of rolling element bearings is essential to assess vibration and noise of drivetrain systems. This study also seeks to fully describe the vibro-acoustic propagation of gear dynamics through a power-transmission system using rolling element and fluid film wave bearings. Fluid film wave bearings, which have higher damping than rolling element bearings, could offer an energy dissipation mechanism that reduces the gearbox noise. The effectiveness of each bearing type to disrupt vibration propagation in explored using multi-body computational models. These models include gears, shafts, rolling element and fluid film wave bearings, and the housing. Radiated noise is mapped from the gearbox surface to surrounding environment. The effectiveness of rolling element and fluid film wave bearings in breaking the vibro-acoustic propagation path from the gear to the housing is investigated.

Parker, Robert G.↗

All-Electric X-Plane, X-57 Mod II Ground Vibration Test

As part of the National Aeronautics and Space Administration New Aviation Horizons initiative to demonstrate and validate future high-impact concepts and technologies, the X-57 Maxwell airplane--the first all-electric X-plane--was conceived to advance research in the area of electric propulsion to show the feasibility of minimizing fuel use, reducing emissions, and lowering noise during flight. Through several configuration modifications to the X-57 airplane, validation of electrical-powered flight with increasing efficiency between each modification when compared to the baseline original airplane is anticipated. In the case of the X-57 Modification II airplane, a ground vibration test was needed to identify the airplane structural modes and use them to update and validate the finite element model. To determine the airworthiness of the airplane the updated finite element model will be utilized to investigate both classical and whirl flutter. The X-57ModificationII ground vibration test was performed by the National Aeronautics and Space Administration Armstrong Flight Research Center Flight Loads Laboratory. This paper will highlight the testing performed to acquire the modal data as well as the results.

Natalie Dawn Spivey↗

Adaptively remeshed multiphysical modeling of resistance forge welding with experimental validation of residual stress fields and measurement processes

Welding processes used in the production of pressure vessels impart residual stresses in the manufactured component. Computational modeling is critical to predicting these residual stress fields and understanding how they interact with notches and flaws to impact pressure vessel durability. Here, in this work, we present a finite element model for a resistance forge weld and validate it using laboratory measurements. Extensive microstructural changes, near-melt temperatures, and large localized deformations along the weld interface pose significant challenges to Lagrangian finite element modeling. The proposed modeling approach overcomes these roadblocks in order to provide a high-fidelity simulation that can predict the residual stress state in the manufactured pressure vessel; a rich microstructural constitutive model accounts for material recrystallization dynamics, a frictional-to-tied contact model is coordinated with the constitutive model to represent interfacial bonding, and adaptive remeshing is employed to alleviate severe mesh distortion. An interrupted-weld approach is applied to the simulation to facilitate comparison to displacement measures. Several techniques are employed for residual stress measurement in order to validate the finite element model: neutron diffraction, the contour method, and the slitting method. Model-measurement comparisons are supplemented with detailed simulations that reflect the configurations of the residual-stress measurement processes themselves. The model results show general agreement with experimental measurements, and we observe some similarities in the features around the weld region. Factors that contribute to model-measurement differences are identified. Finally, we conclude with some discussion of the model development and residual stress measurement strategies, including how to best leverage the efforts put forth here for other weld problems.

36 MATERIALS SCIENCE↗

A Shell/3D Modeling Technique for Delaminations in Composite Laminates

A shell/3D modeling technique was developed for which a local solid finite element model is used only in the immediate vicinity of the delamination front. The goal was to combine the accuracy of the full three-dimensional solution with the computational efficiency of a plate or shell finite element model. Multi-point constraints provide a kinematically compatible interface between the local 3D model and the global structural model which has been meshed with plate or shell finite elements. For simple double cantilever beam (DCB), end notched flexure (ENF), and single leg bending (SLB) specimens, mixed mode energy release rate distributions were computed across the width from nonlinear finite element analyses using the virtual crack closure technique. The analyses served to test the accuracy of the shell/3D technique for the pure mode I case (DCB), mode II case (ENF) and a mixed mode I/II case (SLB). Specimens with a unidirectional layup where the delamination is located between two 0 plies, as well as a multidirectional layup where the delamination is located between two non-zero degree plies, were simulated. For a local 3D model extending to a minimum of about three specimen thicknesses in front of and behind the delamination front, the results were in good agreement with mixed mode strain energy release rates obtained from computations where the entire specimen had been modeled with solid elements. For large built-up composite structures modeled with plate elements, the shell/3D modeling technique offers a great potential, since only a relatively small section in the vicinity of the delamination front needs to be modeled with solid elements.

Krueger, Ronald↗