Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite Element Modeling”

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 613 records · Page 34

Formulation of an experimental substructure model using a Craig-Bampton based transmission simulator

An experimental–analytical substructuring is attractive when there is motivation to replace one or more system subcomponents with an experimental model. This experimentally derived substructure can then be coupled to finite element models of the rest of the structure to predict the system response. The transmission simulator method couples a fixture to the component of interest during a vibration test in order to improve the experimental model for the component. The transmission simulator is then subtracted from the tested system to produce the experimental component. This method reduces ill-conditioning by imposing a least squares fit of constraints between substructure modal coordinates to connect substructures, instead of directly connecting physical interface degrees of freedom. This paper presents an alternative means of deriving the experimental substructure model, in which a Craig–Bampton representation of the transmission simulator is created and subtracted from the experimental measurements. The corresponding modal basis of the transmission simulator is described by the fixed-interface modes, rather than free modes that were used in the original approach. Moreover, these modes do a better job of representing the shape of the transmission simulator as it responds within the experimental system, leading to more accurate results using fewer modes. The new approach is demonstrated using a simple finite element model based example with a redundant interface.

42 ENGINEERING↗

Advances and trends in the development of computational models for tires

Status and some recent developments of computational models for tires are summarized. Discussion focuses on a number of aspects of tire modeling and analysis including: tire materials and their characterization; evolution of tire models; characteristics of effective finite element models for analyzing tires; analysis needs for tires; and impact of the advances made in finite element technology, computational algorithms, and new computing systems on tire modeling and analysis. An initial set of benchmark problems has been proposed in concert with the U.S. tire industry. Extensive sets of experimental data will be collected for these problems and used for evaluating and validating different tire models. Also, the new Aircraft Landing Dynamics Facility (ALDF) at NASA Langley Research Center is described.

Noor, A. K.↗

Non-linear dynamic analysis of geared systems

Under driving conditions, a typical geared system may be subjected to large dynamic loads. Also, the vibration level of the geared system is directly related to the noise radiated from the gear box. The steady state dynamic behavior of the system is examined in order to design reliable and quiet transmissions. The scope is limited to a system containing a spur gear pair with backlash and periodically time varying mesh stiffness, and rolling element bearings with clearance type nonlinearities. The internal static transmission error at the gear mesh, which is of importance from high frequency noise and vibration control view point, is considered in the formulation in sinusoidal or periodic form. A dynamic finite element model of the linear time invariant (LTI) system is developed. Effects of several system parameters, such as torsional and transverse flexibilities of the shafts and prime mover/load inertias, on free and forced vibration characteristics are investigated. Several reduced order LTI models are developed and validated by comparing their eigen solutions with the finite element model results. Using the reduced order formulations, a three degree of freedom dynamic model is developed which includes nonlinearities associated with radical clearances in the radial rolling element bearings, backlash between a spur gear pair and periodically varying gear mesh stiffness. As a limiting case, a single degree of freedom model of the spur gear pair with backlash is considered and mathematical conditions for tooth separation and back collision are defined. Both digital simulation technique and analytical models such as method of harmonic balance and the method of multiple scales were used to develop the steady state frequency response characteristics for various nonlinear and/or time varying cases.

Singh, Rajendra↗

Non-linear dynamic analysis of geared systems, part 2

A good understanding of the steady state dynamic behavior of a geared system is required in order to design reliable and quiet transmissions. This study focuses on a system containing a spur gear pair with backlash and periodically time-varying mesh stiffness, and rolling element bearings with clearance type non-linearities. A dynamic finite element model of the linear time-invariant (LTI) system is developed. Effects of several system parameters, such as torsional and transverse flexibilities of the shafts and prime mover/load inertias, on free and force vibration characteristics are investigated. Several reduced order LTI models are developed and validated by comparing their eigen solution with the finite element model results. Several key system parameters such as mean load and damping ratio are identified and their effects on the non-linear frequency response are evaluated quantitatively. Other fundamental issues such as the dynamic coupling between non-linear modes, dynamic interactions between component non-linearities and time-varying mesh stiffness, and the existence of subharmonic and chaotic solutions including routes to chaos have also been examined in depth.

Singh, Rajendra↗

Combined near- and far-field high-energy diffraction microscopy dataset for Ti-7Al tensile specimen elastically loaded in situ

High-energy diffraction microscopy (HEDM) constitutes a suite of combined X-ray characterization methods, which hold the unique advantage of illuminating the microstructure and micromechanical state of a material during concurrent in situ mechanical deformation. The data generated from HEDM experiments provides a heretofore unrealized opportunity to validate meso-scale modeling techniques, such as crystal plasticity finite element modeling (CPFEM), by explicitly testing the accuracy of these models at the length scales where the models predict their response. Combining HEDM methods with in situ loading under known and controlled boundary conditions represents a significant challenge, inspiring the recent development of a new high-precision rotation and axial motion system for simultaneously rotating and axially loading a sample. In this paper, we describe the initial HEDM dataset collected using this hardware on an alpha-titanium alloy (Ti-7Al) under in situ tensile deformation at the Advanced Photon Source, Argonne National Laboratory. We present both near-field HEDM data that maps out the grain morphology and intragranular crystallographic orientations and far-field HEDM data that provides the grain centroid, grain average crystallographic orientation, and grain average elastic strain tensor for each grain. Finally, we provide a finite element mesh that can be utilized to simulate deformation in the volume of this Ti-7Al specimen.

36 MATERIALS SCIENCE↗

Elastic deformation effects on aerodynamic characteristics for a high-aspect-ratio supercritical-wing model

The results of an investigation of the deformations of a high-aspect-ratio, force/pressure, supercritical-wing model during wind tunnel tests and the effects these deformations have on the wing aerodynamics are presented. A finite element model of the wing was developed, and then, for conditions corresponding to wind tunnel test points, experimental aerodynamic loads and theoretical aerodynamic loads were applied to the finite element model. Comparisons were made between the results of these load conditions for changes in structural deflections and for changes in aerodynamic characteristics. The results show that the deformations are quite small and that the pressure data are not significantly affected by model deformation.

Watson, J. J.↗

Transient thermal analysis of a titanium multiwall thermal protection system

The application of the SPAR thermal analyzer to the thermal analysis of a thermal protection system concept is discussed. The titanium multiwall thermal protection system concept consists of alternate flat and dimpled sheets which are joined together at the crests of the dimples and formed into 30 cm by 30 cm (12 in. by 12 in.) tiles. The tiles are mechanically attached to the structure. The complex tile geometry complicates thermal analysis. Three modes of heat transfer were considered: conduction through the gas inside the tile, conduction through the metal, and radiation between the various layers. The voids between the dimpled and flat sheets were designed to be small enough so that natural convection is insignificant (e.g., Grashof number 1000). A two step approach was used in the thermal analysis of the multiwall thermal protection system. First, an effective normal (through-the-thickness) thermal conductivity was obtained from a steady state analysis using a detailed SPAR finite element model of a small symmetric section of the multiwall tile. This effective conductivity was then used in simple one dimensional finite element models for preliminary analysis of several transient heat transfer problems.

Blosser, M. L.↗

Initial dynamic load estimates during configuration design

This analysis includes the structural response to shock and vibration and evaluates the maximum deflections and material stresses and the potential for the occurrence of elastic instability, fatigue and fracture. The required computations are often performed by means of finite element analysis (FEA) computer programs in which the structure is simulated by a finite element model which may contain thousands of elements. The formulation of a finite element model can be time consuming, and substantial additional modeling effort may be necessary if the structure requires significant changes after initial analysis. Rapid methods for obtaining rough estimates of the structural response to shock and vibration are presented for the purpose of providing guidance during the initial mechanical design configuration stage.

Schiff, Daniel↗

The Simulation of Precipitation Evolutions and Mechanical Properties in Friction Stir Welding with Post-Weld Heat Treatments

A finite element model of friction stir welding capable of re-meshing is used to simulate the temperature variations. Re-meshing of the finite element model is used to maintain a fine mesh resolving the gradients of the solution. The Kampmann–Wagner numerical model for precipitation is then used to study the relation between friction stir welds with post-weld heat treatment (PWHT) and the changes in mechanical properties. Results indicate that the PWHT holding time and PWHT holding temperature need to be optimally designed to obtain FSW with better mechanical properties. Higher precipitate number with lower precipitate sizes gives higher strength in the stirring zone after PWHT. The coarsening of precipitates in HAZ are the main reason to hinder the improvement of mechanical property when PWHT is used.

36 MATERIALS SCIENCE↗

Low-cost fin-tube heat exchanger design for building thermal energy storage using phase change material

Heat transfer in phase change materials (PCMs) is complex because the melting and freezing fronts change as functions of stored or released heat. In prior attempts to optimize heat exchangers (HXs) in one or two dimensions, complex geometry has often been used to maximize the melt and freeze front area. This complex geometry is difficult and hence expensive to construct. This paper proposes a multiple-scale 3D finite element modeling approach to design fin-tube HXs for low-cost latent thermal energy storage applications. Further, the optimal fin and tube designs were determined at three scales (unit-scale, medium-scale, and large-scale) by modeling the melt and freeze front in three dimensions and using measured bulk thermal properties. The finite element model was validated by comparing it with the experimental data for a referenced design of a similar type. The results indicate that commercially available organic PCMs with low conductivity (<0.3 W/m·K) can have charge and discharge times appropriate for building thermal energy storage (i.e., 4–5 h) with fin-tube HX designs at costs <$26/kWh, even when the temperature difference (5.56°C) between the heat transfer fluid and the PCM phase change temperature is small. However, as the HX increases in length, the temperature reduction along the tube limits some larger-scale designs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Development of a Numerical Modeling Approach for Buckling Analysis of Sandwich Composite Cylindrical Shells with Selected Results

The buckling response of geometrically perfect and imperfect cylindrical sandwich shells can be investigated using nonlinear finite element analyses with two-dimensional general-purpose shell elements. Such analyses are used in the NASA Engineering and Safety Center Shell Buckling Knockdown Factor Project, which has the goal of developing new analysis-based buckling design recommendations for select classes of sandwich composite cylindrical structures under uniaxial compressive load. As such, finite element models of sandwich composite cylinders were developed and analyses were performed to predict the buckling responses of geometrically perfect and imperfect sandwich composite cylinders. The development of the selected finite-element modeling approach for a sandwich composite cylinder is discussed. Buckling-response sensitivity of geometrically imperfect sandwich cylinders for various shell element types were investigated as part of this study. Preliminary results of geometric imperfections influence on buckling response of sandwich cylinders are also presented.

Structural Modeling↗

Testing Lithium-Ion Battery with the Internal Reference Electrode: An Insight into the Blocking Effect

Long-term stability of reference electrode is a major issue in lithium-ion batteries. A composite Ni-Li reference electrode (RE) with a longer lifetime than traditional Cu-wire Li reference electrodes is proposed in this study. However, the larger size of the Ni-Li RE creates a blocking effect when positioned internally in the battery. The blocking effect manifests as an abnormal feature in the measured electrode potential. A two-dimensional (2D) finite element model of the cell with an internal RE is constructed to investigate the mechanism of the blocking effect. The internal RE blocks the path of local current flow in the electrolyte, and decreases the electrochemical reaction rate of particles in the anode and cathode, leading to measurement errors. Characteristic lengths of RE and applied currents were examined as parameters in the finite element model to determine their effect on the magnitude of error. A useful guide is provided for proper selection of size of RE and applied current in dynamic tests of lithium-ion batteries.

2-D finite element model↗

Modal Test and Analysis Correlation of Wind Tunnel Blades for the National Transonic Facility

NASA has a significant maintenance challenge with the aging wind tunnel infrastructure across the agency . Wind tunnel blades, such as the ones in the National Transonic Facility at NASA Langley Research Center, have been subjected to damage from impacts and fatigue. Accurate analytical models of these blades are necessary to either recertify replacement blades to be built or for redesign of these blades. In this paper an upgraded finite element model of these wind tunnel blades with the capability of representing damage is presented. An eigenvalue analysis of the structural model is performed, and the boundary conditions are adjusted to simulate the natural frequencies obtained by experimental testing of a set of twenty five blades in 1992. Additional numerical studies of the effect of delamination damage on the natural frequencies of the blade are conducted using the tuned finite element model.

wind tunnel↗

Modal Test and Analysis Correlation of Wind Tunnel Blades for the National Transonic Facility

NASA has a significant maintenance challenge with the aging wind tunnel infrastructure across the agency . Wind tunnel blades, such as the ones in the National Transonic Facility at NASA Langley Research Center, have been subjected to damage from impacts and fatigue. Accurate analytical models of these blades are necessary to either recertify replacement blades to be built or for redesign of these blades. In this paper an upgraded finite element model of these wind tunnel blades with the capability of representing damage is presented. An eigenvalue analysis of the structural model is performed, and the boundary conditions are adjusted to simulate the natural frequencies obtained by experimental testing of a set of twenty five blades in 1992. Additional numerical studies of the effect of delamination damage on the natural frequencies of the blade are conducted using the tuned finite element model.

wind tunnel↗

Postseismic deformation and stress diffusion due to viscoelasticity and comments on the modified Elsasser model

Finite element techniques have been used to investigate the far-field deformation and stress changes due to asthenospheric viscoelastic relaxation following a dip-slip earthquake. The diffusion of extensional stress toward the subduction zone following a thrust earthquake on land is qualitatively consistent with the modified Elsasser model as proposed by Rydelek and Sacks (1988, 1990) to explain the coupled occurrence of land and sea earthquakes near Japan. However, the magnitude of the diffusing tensional signal is significantly smaller. The nominal model consists of a partially faulted elastic lithosphere overlying a viscoelastic substrate. Other models consider thin channel flow, rupturing of the entire elastic lithosphere, and changes in the depth of faulting. While some of these changes have significant impact on the magnitude and spatial features of the stress and deformation field, the far-field stress remains small. Numerical experiments demonstrate that the assumptions of the modified Elsasser model accentuate uniaxial deformation. When these assumptions are replicated in the finite element calculations, reasonable agreement between the models is achieved.

Cohen, Steven C.↗

Predicting Stress vs. Strain Behaviors of Thin-Walled High Pressure Die Cast Magnesium Alloy with Actual Pore Distribution

In this paper, a three-dimensional (3D) microstructure-based finite element modeling method (i.e., extrinsic modeling method) is developed, which can be used in examining the effects of porosity on the ductility/fracture of Mg castings. For this purpose, AM60 Mg tensile samples were generated under high-pressure die-casting in a specially-designed mold. Before the tensile test, the samples were CT-scanned to obtain the pore distributions within the samples. 3D microstructure-based finite element models were then developed based on the obtained actual pore distributions of the gauge area. The input properties for the matrix material were determined by fitting the simulation result to the experimental result of a selected sample, and then used for all the other samples’ simulation. The results show that the ductility and fracture locations predicted from simulations agree well with the experimental results. This indicates that the developed 3D extrinsic modeling method may be used to examine the influence of various aspects of pore sizes/distributions as well as intrinsic properties (i.e., matrix properties) on the ductility/fracture of Mg castings.

Ductility↗

Three-dimensional finite element progressive failure analysis of composite laminates under axial extension

A three-dimensional (3D) progressive failure algorithm is developed, where the layerwise laminate theory (LWLT) of Reddy is used for kinematic description. The finite element model based on the layerwise theory predicts both inplane and interlaminar stresses with the same accuracy as that of a conventional 3D finite element model and provides a convenient format for modeling the 3D stress fields in composite laminates. A parametric study is conducted to investigate the effect of out-of-plane material properties, 3D stiffness reduction methods, and boundary conditions on the failure loads and strains of a composite laminate under axial extension. The results indicate that different parameters have a different degree of influence on the failure loads and strains. The predictive ability of various phenomenological failure criteria is evaluated in the light of experimental results available in the literature, and the predictions of the LWLT are compared with those of the first-order shear deformation theory. It is concluded that a 3D stress analysis is necessary to predict accurately the failure behavior of composite laminates.

Reddy, Yeruva S.↗

Developing Procedures to Implement Geometric Imperfections Beyond Right Circular Cylindrical Shells in Finite Element Method Models

Analysis of aerospace structures is frequently conducted using nominal dimensions and frequently assumes ideal conditions in loading, contact, constraints, et cetera. Off-nominal dimensions and nonideal conditions, however, are present in all structures. These are the result of widely ranging causes from coefficient of thermal expansion mismatches, manufacturing tooling anomalies, to assembly procedures that inadvertently alter the structure. Specifically, geometric imperfections can have potentially significant influence on the response of a structural test article observed in an experiment versus the response given by a numerical simulation. The Python Tool for Implementing Geometric Imperfections in Reduced Structures (Py_TIGIRS) was previously presented as a set of Python scripts to calculate and implement as-manufactured geometric midsurface and thickness imperfections into finite element method (FEM) shell models of nominally right circular cylinders. By taking advantage of the simple shape of a right circular cylinder, interpolations of the measured data points were able to be performed along directions that aligned to the cylindrical coordinate system axes of the entire structure. By taking advantage of the shell representation of the real structure as opposed to modeling using a continuum representation, the thickness variation was able to be implemented by shell section definitions instead of having to modify the position of multiple nodes in the thickness direction. Py_TIGIRS is a useful tool that established a procedural example on how to implement geometric imperfections in right circular cylindrical shell structures. Three new procedures, each expanded from concepts established in Py_TIGIRS, are proposed for various test-article designs and are intended to broaden the range of structures that can be modeled with measured geometric imperfections in the structural analysis community. Each test-article design introduces new challenges to successfully implement geometric imperfections into a FEM model. The first test-article design consists of a carbon fiber reinforced polymer square plate with a hat-shaped stiffener co-cured on one side. This test-article design was for a novel seven-point bend test that was also previously presented. Manufacturing and cure-cycle imperfections are observed using digital image correlation (DIC) techniques. As thermal expansion coefficient mismatches between the plate and stiffener materials were anticipated, a thermal analysis study with continuum shell and solid elements was conducted to capture the global shape observed prior to testing. The second test-article design is of a similar hat-stiffened plate configuration, but with a side length ratio near 3:1 with elongation in the stiffener direction. The test article was used to characterize the response to uniaxial compressive loading in the direction of the stiffener. Due to differing manufacturing steps, a thermal analysis like the one developed for the seven-point bend configuration was unable to mimic the observed geometric imperfections. Instead, a strategy based on applying deformations directly to the structure during analysis was developed for continuum shell and solid element representation of a stiffened panel.

Geometric imperfections↗