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 1,585 records · Page 88

A COMPACT model reduction methodology for articulated, multi-flexible bodies structures

To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one need appropriate reduced-order models for all the flexible components involved. The component modes projection and assembly model reduction methodology is one way to construct these reduced-order component models. However, the generated reduced-order system models typically contain high frequency extraneous modes. The presence of high frequency extraneous modes in a reduced-order system model forces us to use a small stepsize in the numerical integration of that model, resulting in an increase in the elapsed time of the integration process. To eliminate these high-frequency modes, we first identify and discard the component modes that 'contributed' to the generations of these modes. A sensitivity analysis is then used to determine how the remaining component modes should be perturbed in order to still closely capture the selected system modes. The effectiveness of the proposed methodology has been successfully verified using a finite-element model of the Galileo spacecraft.

Lee, Allan Y.↗

Crystal plasticity modeling of β phase deformation in Ti-6Al-4V

Ti-6Al-4V is an alloy of titanium that dominates titanium usage in applications ranging from mass-produced consumer goods to high-end aerospace parts. The material's structure on a microscale is known to affect its mechanical properties but these effects are not fully understood. Specifically, this work will address the effects of low volume fraction intergranular β phase on Ti-6Al-4V's mechanical response during the transition from elastic to plastic deformation. A crystal plasticity-based finite element model is used to fully resolve the deformation of the β phase for the first time. This high fidelity model captures mechanisms difficult to access via experiments or lower fidelity models. Lastly, the results are used to assess lower fidelity modeling assumptions and identify phenomena that have ramifications for failure of the material.

36 MATERIALS SCIENCE↗

Spectroscopic and Computational Analysis of Chemical Pollutants Interactions on Steel - 18409

There is limited understanding of the scientific interaction of the chemical contaminants to DOE facility materials which can alter the stability of the materials, affecting the general safety. The resulting studies will have a direct impact towards developing an enhanced cost-effective decontamination method and understanding the negative effects caused by the chemical waste on the structural integrity of the solid surface in contact. First order kinematic models provide a better fit to the experimental data, implying that one layer of Hg adsorbs to the steel surface. Energy dispersive X-ray analysis confirmed that stainless steel samples are mainly composed of Fe, Ni, and Cr, and have a low carbon content. Objective: To determine the physiochemical interactions of Mercury (Hg) and Cesium (Cs) to steel and concrete using spectroscopic, adsorption and computational studies. To develop a multiscale simulation that shows the effects of adsorption and diffusion of alkali metals on the mechanical performance of steel alloys used in nuclear power plants. Approach: Molecular Dynamics (MD) simulations were used to compute the stress gradient generated in steel due to Hg adsorption. Finite element models were developed for simulating the effects of temperature and pressure in a diffusion analysis. Based on experimental data, one layer of Hg was used for the MD simulations. Material properties at the atomic level were determined from the MD simulations The stress gradient from the MD simulations is used as an input parameter in the FE simulation. Conclusions and Future Work: Computational methods such as MD and FE simulations were used to determine the effects of mercury adsorption in stainless steel samples. The adsorption process induces a stress gradient in the facility materials, which can affect the mechanical performance of facility materials under service. Tensile stresses help the diffusion process while compressive have the opposite effect. Temperature has a positive effect in the mass diffusion process. Future work will evaluate the effects of temperature on the material properties at the molecular level using both experimental and MD methods. Future work will evaluate the mechanical performance of facility materials at a wide range of environmental variables.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A static analysis of three-dimensional functionally graded beams through hierarchical one-dimensional finite elements

In this paper, the static response of three-dimensional beams made of functionally graded materials is investigated through a family of hierarchical one-dimensional finite elements. A wide variety of elements is proposed differing by the kinematic formulation and the number of nodes per elements along the beam axis. Elements’ stiffness matrix and load vector are derived in a unified nuclear form that does not depend upon the a priori expansion order over the cross-section nor the finite element approximation along the beam axis. Results are validated towards three-dimensional finite element models as well as equivalent Navier-type analytical solutions. The numerical investigations show that accurate and efficient solutions (when compared with full three-dimensional FEM solutions) can be obtained by the proposed family of hierarchical one-dimensional elements’ family.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Loss Factor Estimation Using the Impulse Response Decay Method on a Stiffened Structure

High-frequency vibroacoustic modeling is typically performed using energy-based techniques such as Statistical Energy Analysis (SEA). Energy models require an estimate of the internal damping loss factor. Unfortunately, the loss factor is difficult to estimate analytically, and experimental methods such as the power injection method can require extensive measurements over the structure of interest. This paper discusses the implications of estimating damping loss factors using the impulse response decay method (IRDM) from a limited set of response measurements. An automated procedure for implementing IRDM is described and then evaluated using data from a finite element model of a stiffened, curved panel. Estimated loss factors are compared with loss factors computed using a power injection method and a manual curve fit. The paper discusses the sensitivity of the IRDM loss factor estimates to damping of connected subsystems and the number and location of points in the measurement ensemble.

Cabell, Randolph↗

Thermal mechanical analysis of sprag clutches

Work done at Case Western Reserve University on the Thermal Mechanical analysis of sprag helicopter clutches is reported. The report is presented in two parts. The first part is a description of a test rig for the measurement of the heat generated by high speed sprag clutch assemblies during cyclic torsional loading. The second part describes a finite element modeling procedure for sliding contact. The test rig provides a cyclic torsional load of 756 inch-pounds at 5000 rpm using a four-square arrangement. The sprag clutch test unit was placed between the high speed pinions of the circulating power loop. The test unit was designed to have replaceable inner ad outer races, which contain the instrumentation to monitor the sprag clutch. The torque loading device was chosen to be a water cooled magnetic clutch, which is controlled either manually or through a computer. In the second part, a Generalized Eulerian-Lagrangian formulation for non-linear dynamic problems is developed for solid materials. This formulation is derived from the basic laws and axioms of continuum mechanics. The novel aspect of this method is that we are able to investigate the physics in the spatial region of interest as material flows through it without having to follow material points. A finite element approximation to the governing equations is developed. Iterative Methods for the solution of the discrete finite element equations are explored. A FORTRAN program to implement this formulation is developed and a number of solutions to problems of sliding contact are presented.

Mullen, Robert L.↗

Thermal Residual Stress in Environmental Barrier Coated Silicon Nitride - Modeled

When exposed to combustion environments containing moisture both un-reinforced and fiber reinforced silicon based ceramic materials tend to undergo surface recession. To avoid surface recession environmental barrier coating systems are required. However, due to differences in the elastic and thermal properties of the substrate and the environmental barrier coating, thermal residual stresses can be generated in the coated substrate. Depending on their magnitude and nature thermal residual stresses can have significant influence on the strength and fracture behavior of coated substrates. To determine the maximum residual stresses developed during deposition of the coatings, a finite element model (FEM) was developed. Using this model, the thermal residual stresses were predicted in silicon nitride substrates coated with three environmental coating systems namely barium strontium aluminum silicate (BSAS), rare earth mono silicate (REMS) and earth mono di-silicate (REDS). A parametric study was also conducted to determine the influence of coating layer thickness and material parameters on thermal residual stress. Results indicate that z-direction stresses in all three systems are small and negligible, but maximum in-plane stresses can be significant depending on the composition of the constituent layer and the distance from the substrate. The BSAS and REDS systems show much lower thermal residual stresses than REMS system. Parametric analysis indicates that in each system, the thermal residual stresses can be decreased with decreasing the modulus and thickness of the coating.

Ali, Abdul-Aziz↗

European Service Module - Structural Test Article (E-STA) Building Block Test Approach and Model Correlation Observations

The Orion European Service Module - Structural Test Article (E-STA) underwent sine vibration testing in 2016 using the Mechanical Vibration Facility (MVF) multi-axis shaker system at NASA Glenn Research Center’s (GRC) Plum Brook Station (PBS) Space Power Facility (SPF). The main objective was to verify the structural integrity of the European Service Module (ESM) under sine sweep dynamic qualification vibration testing. A secondary objective was to perform a fixed-base modal survey, while E-STA was still mounted to MVF, in order to achieve a test correlate the finite element model (FEM). To facilitate the E-STA system level correlation effort, a building block test approach was implemented. Modal tests were performed on two major subassemblies, the crew module/launch abort structure (CM/LAS) and the crew module adapter (CMA) mass simulators. These subassembly FEMs were individually correlated and then integrated into the E-STA FEM prior to the start of the E-STA sine vibration test. This paper summarizes the modal testing and model correlation efforts of both of these subassemblies and how the building block approach assisted in the overall correlation of the E-STA FEM. This paper will also cover modeling practices that should be avoided, recommended instrumentation positioning on complex structures, and the importance of the FEM geometrically matching CAD in sufficient detail in order to adequately replicate internal load paths. The goal of this paper is to inform the reader of the hard earned lessons learned and pitfalls to avoid when applying a building block test approach.

Winkel, James P.↗

Novel application of the modal strain energy technique for state-of-the-art damping predictions

Structural damping, which measures the energy dissipation of a vibrating structure, is a key modeling input for lightweight structures but is notoriously hard to predict. This work utilizes vibration-based measurements of centimeter-scale coupons and the modal strain energy approach to predict structural damping of a lightly damped structure. The approach was originally validated with panels shorter than a meter in length. This work extends the validation to a 2.75 m beam made of unidirectional and biaxial glass fiber laminates bonded by adhesive. The comparison between three-dimensional finite element model predictions and full-scale experimental measurements of damping show an average error of 5.2% for the first five modes. Additionally, the modal strain energy approach is newly applied with a one-dimensional geometrically exact beam theory model and a two-dimensional sectional analysis solver. This beam approach accurately predicts the damping behavior of the first bending modes but loses accuracy for higher order modes that are dominated by three-dimensional effects. This novel approach provides faster simulations while allowing arbitrary beam cross sections. The paper also investigates traditional and high-force dynamic mechanical analysis to measure structural damping of coupons. Both alternatives show significant errors in attempted validation against the theoretical thermoelastic damping of aluminum coupons.

17 WIND ENERGY↗

Mechanical Design Studies of the MQXF Long Model Quadrupole for the HiLumi LHC

The Large Hadron Collider Luminosity upgrade (HiLumi) program requires new low-β triplet quadrupole magnets, called MQXF, in the Interaction Region (IR) to increase the LHC peak and integrated luminosity. The MQXF magnets, designed and fabricated in collaboration between CERN and the U.S. LARP, will all have the same cross section. The MQXF long model, referred as MQXFA, is a quadrupole using the Nb3Sn superconducting technology with 150 mm aperture and a 4.2 m magnetic length and is the first long prototype of the final MQXF design. The MQXFA magnet is based on the previous LARP HQ and MQXFS designs. In this paper we present the baseline design of the MQXFA structure with detailed 3D numerical analysis. A detailed tolerance analysis of the baseline case has been performed by using a 3D finite element model, which allows fast computation of structures modelled with actual tolerances. Tolerance sensitivity of each component is discussed to verify the actual tolerances to be achieved by vendors. In conclusion, tolerance stack-up analysis is presented in the end of this paper.

43 PARTICLE ACCELERATORS↗

Correlation of AH-1G airframe flight vibration data with a coupled rotor-fuselage analysis

The formulation and features of the Rotor-Airframe Comprehensive Analysis Program (RACAP) is described. The analysis employs a frequency domain, transfer matrix approach for the blade structural model, a time domain wake or momentum theory aerodynamic model, and impedance matching for rotor-fuselage coupling. The analysis is applied to the AH-1G helicopter, and a correlation study is conducted on fuselage vibration predictions. The purpose of the study is to evaluate the state-of-the-art in helicopter fuselage vibration prediction technology. The fuselage vibration predicted using RACAP are fairly good in the vertical direction and somewhat deficient in the lateral/longitudinal directions. Some of these deficiencies are traced to the fuselage finite element model.

Sangha, K.↗

Modeling and control for vibration suppression of a flexible smart structure

Theoretical and experimental results of the modeling and control of a flexible ribbed antenna are presented. The antenna consists of eight flexible ribs which constitutes a smart antenna in the sense that the actuator and sensors are an integral part of the structure. The antenna exhibits closely space and repeated modes, thus multi-input multi-output (MIMO) control is necessary for controllability and observability of the structure. The structure also exhibits mode localization phenomenon and contains post buckled members making an accurate finite element model of the structure difficult to obtain. An identified MIMO minimum order model of the antenna is synthesized from identified single-input single-output (SISO) transfer functions curve fit in the frequency domain. The identified model is used to design a positive position feedback (PPF) controller that increases damping in all of the modes in the targeted frequency range. Due to the accuracy of the open loop model of the antenna, the closed loop response predicted by the identified model correlates well wtih experimental results.

Dosch, J.↗

Acoustic Modeling and Analysis for the Space Shuttle Main Propulsion System Liner Crack Investigation

Cracks were found on bellows flow liners in the liquid hydrogen feedlines of several space shuttle orbiters in 2002. An effort to characterize the fluid environment upstream of the space shuttle main engine low-pressure fuel pump was undertaken to help identify the cause of the cracks and also provide quantitative environments and loads of the region. Part of this effort was to determine the duct acoustics several inches upstream of the low-pressure fuel pump in the region of a bellows joint. A finite element model of the complicated geometry was made using three-dimensional fluid elements. The model was used to describe acoustics in the complex geometry and played an important role in the investigation. Acoustic mode shapes and natural frequencies of the liquid hydrogen in the duct and in the cavity behind the flow liner were determined. Forced response results were generated also by applying an edgetone-like forcing to the liner slots. Studies were conducted for state conditions and also conditions assuming two-phase entrapment in the backing cavity. Highly instrumented single-engine hot fire data confirms the presence of some of the predicted acoustic modes.

Casiano, Matthew J.↗

Select Variables Affecting Thermal System Design of a Liquid-Cooled Stator

NASA Glenn Researchers have endeavored to create a high-power-density, efficient electric motor called the High Efficiency Megawatt Motor (HEMM) with a goal of exceeding 98% efficiency and 1.46 MW of power. These aggressive goals were set because increasing the power density and efficiency of electric machines (motors and generators) is integral to bringing electrified aircraft (EA) to commercial realization. Resistive losses in the stator windings are by far the largest (34%) contributor to this design’s total energy losses. Previous testing of a HEMM-like stator test article (statorette) showed good agreement between measured performance and model predictions. Therefore, a foundation was provided to test other engineering parameters that are encountered when designing an electric machine. Of particular interest are the effects of changes to the thickness of the epoxy from the potting process, changes in the fluid volume/flow path around the end windings, and changes in the viscosity of the cooling fluid on the overall thermal performance of the unit. The impact of those parameters as revealed through a third-generation build, test, and finite element model analysis of a statorette test section is the focus of this paper.

EAP↗

Nonlinear shell analyses of the SRB/ETA ring interface

Two-dimensional shell models of the portion of the solid rocket booster (SRB) with the external tank attachment (ETA) ring were developed using the STAGSC-1 computer code. These analyses were performed in support of the SRB recertification program underway at the NASA Marshall Space Flight Center. The purpose of these analyses is to calculate the overall stress and deflection distributions for the SRB/ETA ring interface. An overview of the SRB/ETA ring geometry is presented followed by a discussion of the analysis methods and finite element models. Finally, the nonlinear structural response of the SRB/ETA ring interface is discussed.

Knight, Norman F., Jr.↗

A geometric nonlinear degenerated shell element using a mixed formulation with independently assumed strain fields

A mixed formulation is chosen to overcome deficiencies of the standard displacement-based shell model. Element development is traced from the incremental variational principle on through to the final set of equilibrium equations. Particular attention is paid to developing specific guidelines for selecting the optimal set of strain parameters. A discussion of constraint index concepts and their predictive capability related to locking is included. Performance characteristics of the elements are assessed in a wide variety of linear and nonlinear plate/shell problems. Despite limiting the study to geometric nonlinear analysis, a substantial amount of additional insight concerning the finite element modeling of thin plate/shell structures is provided. For example, in nonlinear analysis, given the same mesh and load step size, mixed elements converge in fewer iterations than equivalent displacement-based models. It is also demonstrated that, in mixed formulations, lower order elements are preferred. Additionally, meshes used to obtain accurate linear solutions do not necessarily converge to the correct nonlinear solution. Finally, a new form of locking was identified associated with employing elements designed for biaxial bending in uniaxial bending applications.

Graf, Wiley E.↗

Meissner-London Susceptibility of Superconducting Right Circular Cylinders in an Axial Magnetic Field

Analysis of magnetic susceptibility of nonellipsoidal samples is a long-standing problem in experimental studies of magnetism and superconductivity. Here the quantitative description of the Meissner-London response (no Abrikosov vortices) of right circular cylinders in an axial magnetic field is given. Three-dimensional adaptive finite-element modeling is used to calculate the total magnetic moment ( m ) in a wide range of London penetration depth ( λ ) to sample size ( a ) ratios, λ /a. Furthermore, by fitting numerical solutions, the closed-form empirical magnetic susceptibility is obtained involving only sample dimensions and λ, thus providing a recipe for determining the London penetration depth from the precision measurements of the magnetic susceptibility. Detailed protocols of the experimental data analysis using the developed approach are given. The results can be readily extended to the most frequently used cuboid-shaped samples.

3-dimensional systems↗

Mechanical Modeling of a WIPP Drum Under Pressure

Mechanical modeling was undertaken to support the Waste Isolation Pilot Plant (WIPP) technical assessment team (TAT) investigating the February 14th 2014 event where there was a radiological release at the WIPP. The initial goal of the modeling was to examine if a mechanical model could inform the team about the event. The intention was to have a model that could test scenarios with respect to the rate of pressurization. It was expected that the deformation and failure (inability of the drum to contain any pressure) would vary according to the pressurization rate. As the work progressed there was also interest in using the mechanical analysis of the drum to investigate what would happen if a drum pressurized when it was located under a standard waste package. Specifically, would the deformation be detectable from camera views within the room. A finite element model of a WIPP 55-gallon drum was developed that used all hex elements. Analyses were conducted using the explicit transient dynamics module of Sierra/SM to explore potential pressurization scenarios of the drum. Theses analysis show similar deformation patterns to documented pressurization tests of drums in the literature. The calculated failure pressures from previous tests documented in the literature vary from as little as 16 psi to 320 psi. In addition, previous testing documented in the literature shows drums bulging but not failing at pressures ranging from 69 to 138 psi. The analyses performed for this study found the drums failing at pressures ranging from 35 psi to 75 psi. When the drums are pressurized quickly (in 0.01 seconds) there is significant deformation to the lid. At lower pressurization rates the deformation of the lid is considerably less, yet the lids will still open from the pressure. The analyses demonstrate the influence of pressurization rate on deformation and opening pressure of the drums. Analyses conducted with a substantial mass on top of the closed drum demonstrate that the drums will still open provided the pressure is high enough. Investigation teams should look for displaced drum lids when searching for drums that have pressurized and failed. The mechanical modeling study for this program is summarized in the following memo. Following a brief introduction, there is a summary of a brief literature review of previous pressure testing of drums, an explanation of the model, presentation of the key results, some discussion, and concluding with a summary and key points.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗