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,027 records · Page 57

Vibration, performance, flutter and forced response characteristics of a large-scale propfan and its aeroelastic model

An investigation of the vibration, performance, flutter, and forced response of the large-scale propfan, SR7L, and its aeroelastic model, SR7A, has been performed by applying available structural and aeroelastic analytical codes and then correlating measured and calculated results. Finite element models of the blades were used to obtain modal frequencies, displacements, stresses and strains. These values were then used in conjunction with a 3-D, unsteady, lifting surface aerodynamic theory for the subsequent aeroelastic analyses of the blades. The agreement between measured and calculated frequencies and mode shapes for both models is very good. Calculated power coefficients correlate well with those measured for low advance ratios. Flutter results show that both propfans are stable at their respective design points. There is also good agreement between calculated and measured blade vibratory strains due to excitation resulting from yawed flow for the SR7A propfan. The similarity of structural and aeroelastic results show that the SR7A propfan simulates the SR7L characteristics.

August, Richard↗

Vibration, performance, flutter and forced response characteristics of a large-scale propfan and its aeroelastic model

An investigation of the vibration, performance, flutter, and forced response of the large-scale propfan, SR7L, and its aeroelastic model, SR7A, has been performed by applying available structural and aeroelastic analytical codes and then correlating measured and calculated results. Finite element models of the blades were used to obtain modal frequencies, displacements, stresses and strains. These values were then used in conjunction with a 3-D, unsteady, lifting surface aerodynamic theory for the subsequent aeroelastic analyses of the blades. The agreement between measured and calculated frequencies and mode shapes for both models is very good. Calculated power coefficients correlate well with those measured for low advance ratios. Flutter results show that both propfans are stable at their respective design points. There is also good agreement between calculated and measured blade vibratory strains due to excitation resulting from yawed flow for the SR7A propfan. The similarity of structural and aeroelastic results show that the SR7A propfan simulates the SR7L characteristics.

August, Richard↗

Model Of Bearing With Hydrostatic Damper

Improved mathematical model of rotational and vibrational dynamics of bearing package in turbopump incorporates effects of hydrostatic damper. Part of larger finite-element model representing rotational and vibrational dynamics of rotor and housing of pump. Includes representations of deadband and nonlinear stiffness and damping of ball bearings, nonlinear stiffness and damping of hydrostatic film, and stiffness of bearing support. Enables incorporation of effects of hydrostatic damper into overall rotor-dynamic mathematical model without addition of mathematical submodel of major substructure.

Goggin, David G.↗

The effect of friction in the hold down post spherical bearings on hold down post loads

The effect of friction at the connection of the Solid Rocket Booster (SRB) aft skirt and the mobile launch platform (MLP) hold down posts was analyzed. A simplified model of the shuttle response during the Space Shuttle Main Engine (SSME) buildup was constructed. The model included the effect of stick-slip friction for the rotation of the skirt about the spherical bearing. Current finite element models assume the joint is completely frictionless in rotation and therefore no moment is transferred between the skirt and the hold down posts. The model was partially verified against test data and preliminary parameter studies were performed. The parameter studies indicated that the coefficient of friction strongly influenced the moment on the hold down posts. The coefficient of friction had little effect on hold down post vertical loads, however. Further calibration of the model is necessary before the effect of friction on the hold down post horizontal loads can be analyzed.

Richardson, James A.↗

Enhancements to NURBS-Based FEA Airfoil Modeler: SABER

NURBS (Non-Uniform Rational B-Splines) have become a common way for CAD programs to fit a smooth surface to discrete geometric data. This concept has been extended to allow for the fitting of analysis data in a similar manner and "attaching" the analysis data to the geometric definition of the structure. The "attaching" of analysis data to the geometric definition allows for a more seamless sharing of data between analysis disciplines. NURBS have become a useful tool in the modeling of airfoils. The use of NURBS has allowed for the development of software that easily and consistently generates plate finite element models of the midcamber surface of a given airfoil. The resulting displacements can then be applied to the original airfoil surface and the deformed shape calculated.

Saleeb, A. F.↗

Orion Parachute Riser Cutter Development

This paper presents the tests and analytical approach used on the development of a steel riser cutter for the CEV Parachute Assembly System (CPAS) used on the Orion crew module. Figure 1 shows the riser cutter and the steel riser bundle which consists of six individual cables. Due to the highly compressed schedule, initial unavailability of the riser material and the Orion Forward Bay mechanical constraints, JSC primarily relied on a combination of internal ballistics analysis and LS-DYNA simulation for this project. Various one dimensional internal ballistics codes that use standard equation of state and conservation of energy have commonly used in the development of CAD devices for initial first order estimates and as an enhancement to the test program. While these codes are very accurate for propellant performance prediction, they usually lack a fully defined kinematic model for dynamic predictions. A simple piston device can easily and accurately be modeled using an equation of motion. However, the accuracy of analytical models is greatly reduced on more complicated devices with complex external loads, nonlinear trajectories or unique unlocking features. A 3D finite element model of CAD device with all critical features included can vastly improve the analytical ballistic predictions when it is used as a supplement to the ballistic code. During this project, LS-DYNA structural 3D model was used to predict the riser resisting load that was needed for the ballistic code. A Lagrangian model with eroding elements shown in Figure 2 was used for the blade, steel riser and the anvil. The riser material failure strain was fine tuned by matching the dent depth on the anvil with the actual test data. LS-DYNA model was also utilized to optimize the blade tip design for the most efficient cut. In parallel, the propellant type and the amount were determined by using CADPROG internal ballistics code. Initial test results showed a good match with LS-DYNA and CADPROG simulations. Final paper will present a detailed roadmap from initial ballistic modeling and LS-DYNA simulation to the performance testing. Blade shape optimization study will also be presented.

Oguz, Sirri↗

Inverse Force Determination on a Small Scale Launch Vehicle Model Using a Dynamic Balance

A launch vehicle can experience large unsteady aerodynamic forces in the transonic regime that, while usually only lasting for tens of seconds during launch, could be devastating if structural components and electronic hardware are not designed to account for them. These aerodynamic loads are difficult to experimentally measure and even harder to computationally estimate. The current method for estimating buffet loads is through the use of a few hundred unsteady pressure transducers and wind tunnel test. Even with a large number of point measurements, the computed integrated load is not an accurate enough representation of the total load caused by buffeting. This paper discusses an attempt at using a dynamic balance to experimentally determine buffet loads on a generic scale hammer head launch vehicle model tested at NASA Ames Research Center's 11' x 11' transonic wind tunnel. To use a dynamic balance, the structural characteristics of the model needed to be identified so that the natural modal response could be and removed from the aerodynamic forces. A finite element model was created on a simplified version of the model to evaluate the natural modes of the balance flexures, assist in model design, and to compare to experimental data. Several modal tests were conducted on the model in two different configurations to check for non-linearity, and to estimate the dynamic characteristics of the model. The experimental results were used in an inverse force determination technique with a psuedo inverse frequency response function. Due to the non linearity, the model not being axisymmetric, and inconsistent data between the two shake tests from different mounting configuration, it was difficult to create a frequency response matrix that satisfied all input and output conditions for wind tunnel configuration to accurately predict unsteady aerodynamic loads.

testing↗

Aeroelastic Modeling of X-56A Stiff-Wing Configuration Flight Test Data

Aeroelastic stability and control derivatives for the X-56A Multi-Utility Technology Testbed (MUTT), in the stiff-wing configuration, were estimated from flight test data using the output-error method. Practical aspects of the analysis are discussed. The orthogonal phase-optimized multisine inputs provided excellent data information for aeroelastic modeling. Consistent parameter estimates were determined using output error in both the frequency and time domains. The frequency domain analysis converged faster and was less sensitive to starting values for the model parameters, which was useful for determining the aeroelastic model structure and obtaining starting values for the time domain analysis. Including a modal description of the structure from a finite element model reduced the complexity of the estimation problem and improved the modeling results. Effects of reducing the model order on the short period stability and control derivatives were investigated.

Grauer, Jared A.↗

Relevence of 3D Simulations and Sandwich Core Topology for the Modeling of Honeycomb Core Sandwich Panels Undergoing Mode I Delamination

A recently developed cohesive zone traction-separation law, which includes the effects of fiber bridging in a novel way, is extended from 2D to 3D. The proposed cohesive model is applied to low fidelity (i.e. homogenized core representation) and high fidelity (i.e. directly accounting for the core topology) finite element models of a composite panel comprised of carbon fiber reinforced plastic facesheets and a honeycomb sandwich core. This enables the investigation of 2D to 3D parameter transferability, width-dependent effects such as thumbnail-shaped crack growth, and the verification of plane strain / plane stress assumptions. A pronounced curvature of the initial interface-related crack front is observed, while the bridging-related crack front is straight. Furthermore, it is found that the cohesive parameters can easily be transferred from 2D to 3D under plane stress assumptions, but not under plane strain assumptions. The numerical predictions are compared to experimental load-displacement and R-curves.

Composites↗

Multiple Shaker Placement for Ground Vibration Test of X-59 Aircraft Using Topology Optimization

A multiple shaker placement methodology is developed and tested using a topology optimization technique. Current multiple shaker placement methodology requires optimum accelerometer placement and optimum single-shaker placement techniques. The proposed methodology is tested using a finite element model of the X-59 Low Boom Flight Demonstrator aircraft. The effective independence and the driving point acceleration transfer function (DPATF) methods are used for the accelerometer placement study. In this study, four shakers are used to excite each mode more effectively during the ground vibration test; all the modes of interest thus are separated into four groups. Each shaker takes care of a separate group of modes. Grouping the modes of interest is performed utilizing topology optimization. The number of modes for each group therefore will be automatically decided during grouping. For each group of modes, perform the following two steps to determine optimal location of four shakers: 1) At each accelerometer location, compare the magnitude of DPATF values at natural frequencies, select the minimum value, and make a vector with these minimum values of the DPATF magnitudes for each group; and 2) Select the degrees of freedom corresponding to the maximum value of this vector. The objective function value is the maximum value of the vector with minimum value of the magnitude of the superposed acceleration transfer function. This objective function value is maximized by changing the modes for each group. Forty accelerometers are enough to have good correlation between mode shapes obtained from the reduced order model and the simulated ground vibration test.

Optimization Analysis↗

Finite Element Prediction of Acoustic Scattering and Radiation from Submerged Elastic Structures

A finite element formulation is derived for the scattering and radiation of acoustic waves from submerged elastic structures. The formulation uses as fundamental unknowns the displacement in the structure and a velocity potential in the field. Symmetric coefficient matrices result. The outer boundary of the fluid region is terminated with an approximate local wave-absorbing boundary condition which assumes that outgoing waves are locally planar. The finite element model is capable of predicting only the near-field acoustic pressures. Far-field sound pressure levels may be determined by integrating the surface pressures and velocities over the wet boundary of the structure using the Helmholtz integral. Comparison of finite element results with analytic results show excellent agreement. The coupled fluid-structure problem may be solved with general purpose finite element codes by using an analogy between the equations of elasticity and the wave equation of linear acoustics.

Everstine, G. C.↗

Elasto-plastic coupled temperature-displacement finite element analysis of two-dimensional rolling-sliding contact with a translating heat source

The present paper describes a transient translating elastoplastic thermomechanical finite element model to study two-dimensional frictional rolling contact. Frictional two-dimensional contact is simulated by repeatedly translating a nonuniform thermomechanical distribution across the surface of an elastoplastic half space. The half space is represented by a two-dimensional finite element mesh with appropriate boundaries. Calculations are for an elastic-perfectly plastic material and the selected thermophysical properties are assumed to be temperature independent. The paper presents temperature variations, stress and plastic strain distributions and deformations. Residual tensile stresses are observed. The magnitude and depth of these stresses depends on (1) the temperature gradients and (2) the magnitudes of the normal and tangential tractions.

Rubin, C. A.↗

Numerical and experimental analysis of mechanically induced failure in electric vehicle battery modules

Mitigating thermal runaway and cell-to-cell propagation is essential for improving the safety of electric and hybrid vehicles. Enhancing digital twin capabilities to predict battery mechanical abuse is particularly critical for automotive and aerospace applications, where crashworthiness is a key concern. Understanding failure conditions and propagation in battery modules during mechanical abuse is complex due to interactions between structural deformation, heat transfer, electrochemical processes, exothermic reactions and mechanical fracture. While prior studies have focused on modeling cell-level behavior, extending these models to module or pack level is necessary for a system level understating of electric vehicle safety. This study develops coupled large deformation finite element models that simultaneously solve for electrochemistry, material failure, internal short circuit and thermal runaway propagation. The models account for mechanical and thermal interactions between lithium-ion cells and other battery components while the contact interfaces are evolving with time. Model-predicted voltage, temperature and force responses are compared with experimental data for validation. The results demonstrate that the approach captures key failure mechanisms, including thermal propagation through heat transfer, electrical propagation from short circuits in parallel-connected cells, and mechanical propagation via penetration and crack formation. These findings show that computational models are valuable tools for understanding battery module failure and providing insight that can reduce the need for extensive experimental testing.

25 ENERGY STORAGE↗

Design and Analysis of the Rocket Experiment Demonstration of a Soft X-Ray (REDSoX) Polarimeter Mirror Module Assembly

The Rocket Experiment Demonstration of a Soft X-ray (REDSoX) Polarimeter is a NASA sounding rocket experiment that is designed to demonstrate the technology necessary for measuring linear X-ray polarization as a function of energy below 1 keV. In astrophysics, soft X-ray spectropolarimetry will be used to probe the nature of acceleration mechanisms in quasar jets and to test models of neutron star structure. NASA Marshall Space Flight Center (MSFC) has designed a grazing-incidence mirror module assembly (MMA) for the REDSoX payload and used Finite Element Modelling Software (ANSYS) to perform structural analysis of the design. In this paper we will describe the overall design of the REDSoX MMA, details of the analysis techniques used for predicting factor of safety in the mirror adhesive bonds and structural components, the buckling analysis of the outer housing, and the raytrace technique used to estimate the effect of gravity sag on optical performance during ground testing.

Structural Analysis↗

Micro-beam bending of FCC bicrystals: A comparison between defect dynamics simulations and experiments

To understand the role of the grain boundary (GB) in plasticity at small scale, a concurrently coupled mesoscale plasticity model was developed to simulate micro-bending of bicrystalline micron-sized beams. By coupling dislocation dynamics (DD) with a finite element model (FEM), a novel defect dynamics model provides the means to investigate intricate interactions between dislocations and GBs under various loading conditions. Our simulations of micro-bending agree well with corresponding micro-bending experiments, and they show that mechanical response of bicrystals could have not only hardening but also softening depending on the characters of the GB. In addition, changing the location of the GB in the microbeams results in different mechanical responses; GBs located at the neutral plane show softening compared to single crystals, while inclined GBs located halfway along the length of the beam show little effect. In conclusion, simulation results could provide a clear picture on detailed dislocation-GB interactions, and quantitative resolved shear stress analysis supplemented by dislocation density distribution is used to analyze the mechanical response of bicrystalline samples.

36 MATERIALS SCIENCE↗

Development and testing of a solar flux heating freeze recovery system for molten salt parabolic troughs

A major obstacle preventing commercialization of molten salt parabolic trough plants is freeze protection and recovery of the solar field. One of the most common nitrate salt compositions under consideration allows for solar field outlet temperatures up to 565°C and begins solidifying above 240°C. Using molten salt directly in the solar field can significantly reduce thermal energy storage costs and the higher temperature compared to oil heat transfer fluids can power a more efficient power cycle, but the very high freezing point is cause for concern. This work builds on a previous modeling effort investigating different methods for melting Solar Salt frozen in the solar field, which suggested the viability of a novel solar heating method to lower costs of the freeze protection system. In this work higher fidelity models are developed to better investigate the melting process with a more detailed 3D geometry which includes the insulated bellows. The thermal-fluid model of the melting process is updated to include temperature dependent density and data for solid phase properties. In addition, a finite element model is developed to resolve thermal stresses at the point with highest thermal gradients. Results with the new model confirm the viability of using controllable solar flux heating to thaw salt frozen in the solar field without damaging the receiver; however, the presence of non-illuminated sections significantly slows down the melting process. Solar heating simulations with pauses off-sun of 120 and 60 s require 18.5 and 11.9 h to melt salt from a night-time temperature of 10°C with 1000 W m2 DNI, which results in maximum thermal stresses of 48.2 and 53.5 MPa, respectively. Adding 150 W m−1 to the solar heating simulation with a 60 s pause off-sun decreases both the melt time and thermal stresses due to the more uniform heating, suggesting a combined heating method may be the best option.

14 SOLAR ENERGY↗

Effects of damping on mode shapes, volume 1

Displacement, velocity, and acceleration admittances were calculated for a realistic NASTRAN structural model of space shuttle for three conditions: liftoff, maximum dynamic pressure and end of solid rocket booster burn. The realistic model of the orbiter, external tank, and solid rocket motors included the representation of structural joint transmissibilities by finite stiffness and damping elements. Methods developed to incorporate structural joints and their damping characteristics into a finite element model of the space shuttle, to determine the point damping parameters required to produce realistic damping in the primary modes, and to calculate the effect of distributed damping on structural resonances through the calculation of admittances.

Gates, R. M.↗

Normal mode analysis of the IUS/TDRS payload in a payload canister/transporter environment

Special modeling techniques were developed to simulate an accurate mathematical model of the transporter/canister/payload system during ground transport of the Inertial Upper Stage/Tracking and Data Relay Satellite (IUS/TDRS) payload. The three finite element models - the transporter, the canister, and the IUS/TDRS payload - were merged into one model and used along with the NASTRAN normal mode analysis. Deficiencies were found in the NASTRAN program that make a total analysis using modal transient response impractical. It was also discovered that inaccuracies may exist for NASTRAN rigid body modes on large models when Given's method for eigenvalue extraction is employed. The deficiencies as well as recommendations for improving the NASTRAN program are discussed.

Meyer, K. A.↗