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,297 records · Page 72

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↗

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↗

Simulating the effect of grain structure and porosity on creep for powder bed fusion 316H

This report details results from modeling studies focused on identifying the potential factors influencing creep performance of 316H stainless steel manufactured using Laser Powder Bed Fusion (LPBF), with a focus on experimentally-observed and predicted differences between the response of the AM material and conventionally manufactured wrought 316H. The studies presented here systematically look at the role of grain structure, porosity and texture in the anisotropic creep behavior of the material under conditions expected in high temperature advanced nuclear reactors. A physics-based Crystal Plasticity Finite Element modeling approach has allowed us to look at these factors in isolation and study their role in the material's long-term deformation behavior. An update to the precipitation model used in the constitutive framework for 316H is also presented here, with the aim of improving the accuracy of our results from this modeling effort. This work is a step in the direction of gaining mechanistic understanding of the long-term deformation behavior of LPBF 316H, helping us model its long term performance and reducing the qualification time for the material.

36 MATERIALS SCIENCE↗

Advances in reduction techniques for tire contact problems

Some recent developments in reduction techniques, as applied to predicting the tire contact response and evaluating the sensitivity coefficients of the different response quantities, are reviewed. The sensitivity coefficients measure the sensitivity of the contact response to variations in the geometric and material parameters of the tire. The tire is modeled using a two-dimensional laminated anisotropic shell theory with the effects of variation in geometric and material parameters, transverse shear deformation, and geometric nonlinearities included. The contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the stress resultants, the generalized displacements, and the Lagrange multipliers associated with the contact conditions. The elemental arrays are obtained by using a modified two-field, mixed variational principle. For the application of reduction techniques, the tire finite element model is partitioned into two regions. The first region consists of the nodes that are likely to come in contact with the pavement, and the second region includes all the remaining nodes. The reduction technique is used to significantly reduce the degrees of freedom in the second region. The effectiveness of the computational procedure is demonstrated by a numerical example of the frictionless contact response of the space shuttle nose-gear tire, inflated and pressed against a rigid flat surface. Also, the research topics which have high potential for enhancing the effectiveness of reduction techniques are outlined.

Noor, Ahmed K.↗

Reduction technique for tire contact problems

A reduction technique and a computational procedure are presented for predicting the tire contact response and evaluating the sensitivity coefficients of the different response quantities. The sensitivity coefficients measure the sensitivity of the contact response to variations in the geometric and material parameters of the tire. The tire is modeled using a two-dimensional laminated anisotropic shell theory with the effects of variation in geometric and material parameters, transverse shear deformation, and geometric nonlinearities included. The contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the stress resultants, the generalized displacements, and the Lagrange multipliers associated with the contact conditions. The elemental arrays are obtained by using a modified two-field, mixed variational principle. For the application of the reduction technique, the tire finite element model is partitioned into two regions. The first region consists of the nodes that are likely to come in contact with the pavement, and the second region includes all the remaining nodes. The reduction technique is used to significantly reduce the degrees of freedom in the second region. The effectiveness of the computational procedure is demonstrated by a numerical example of the frictionless contact response of the space shuttle nose-gear tire, inflated and pressed against a rigid flat surface.

Noor, Ahmed K.↗

Reduced-Order Blade Mistuning Analysis Techniques Developed for the Robust Design of Engine Rotors

The primary objective of this research program is to develop vibration analysis tools, design tools, and design strategies to significantly improve the safety and robustness of turbine engine rotors. Bladed disks in turbine engines always feature small, random blade-to-blade differences, or mistuning. Mistuning can lead to a dramatic increase in blade forced-response amplitudes and stresses. Ultimately, this results in high-cycle fatigue, which is a major safety and cost concern. In this research program, the necessary steps will be taken to transform a state-of-the-art vibration analysis tool, the Turbo-Reduce forced-response prediction code, into an effective design tool by enhancing and extending the underlying modeling and analysis methods. Furthermore, novel techniques will be developed to assess the safety of a given design. In particular, a procedure will be established for using eigenfrequency curve veerings to identify "danger zones" in the operating conditions--ranges of rotational speeds and engine orders in which there is a great risk that the rotor blades will suffer high stresses. This work also will aid statistical studies of the forced response by reducing the necessary number of simulations. Finally, new strategies for improving the design of rotors will be pursued. Several methods will be investigated, including the use of intentional mistuning patterns to mitigate the harmful effects of random mistuning, and the modification of disk stiffness to avoid reaching critical values of interblade coupling in the desired operating range. Recent research progress is summarized in the following paragraphs. First, significant progress was made in the development of the component mode mistuning (CMM) and static mode compensation (SMC) methods for reduced-order modeling of mistuned bladed disks (see the following figure). The CMM method has been formalized and extended to allow a general treatment of mistuning. In addition, CMM allows individual mode mistuning, which accounts for the realistic effects of local variations in blade properties that lead to different mistuning values for different mode types (e.g., mistuning of the first torsion mode versus the second flexural mode). The accuracy and efficiency of the CMM method and the corresponding Turbo-Reduce code were validated for an example finite element model of a bladed disk.

Min, James B.↗

MSC/NASTRAN Stress Analysis of Complete Models Subjected to Random and Quasi-Static Loads

Space payloads, such as those which fly on the Space Shuttle in Spacelab, are designed to withstand dynamic loads which consist of combined acoustic random loads and quasi-static acceleration loads. Methods for computing the payload stresses due to these loads are well known and appear in texts and NASA documents, but typically involve approximations such as the Miles' equation, as well as possible adjustments based on "modal participation factors." Alternatively, an existing capability in MSC/NASTRAN may be used to output exact root mean square [rms] stresses due to the random loads for any specified elements in the Finite Element Model. However, it is time consuming to use this methodology to obtain the rms stresses for the complete structural model and then combine them with the quasi-static loading induced stresses. Special processing was developed as described here to perform the stress analysis of all elements in the model using existing MSC/NASTRAN and MSC/PATRAN and UNIX utilities. Fail-safe and buckling analyses applications are also described.

Hampton, Roy W.↗

Structural Characterization of a Toppling Rock Slab From Array-Based Ambient Vibration Measurements and Numerical Modal Analysis

Accurate assessments of the internal structure and boundary conditions of unstable rock slopes are imperative for evaluating landslide hazard scenarios. However, instability characterization at depth remains challenging and is often limited by costly or invasive subsurface investigations. Here, we develop a new approach coupling array-based ambient vibration modal analysis and numerical modeling to improve structural characterization of rock slope instabilities at depth. We used ambient noise cross-correlation on 4 hr of seismic data recorded by an array of 30 nodal geophones at a 500-m-long toppling rock slab in Utah, USA to identify modal frequencies between 0.8 and 3.5 Hz and derive modal displacements. We show that transverse and longitudinal bending modes span the length of the instability, indicating an interconnected slab. Further, statistical comparison of field results with outputs from >1,000 finite element models with varying boundary conditions showed that the instability depth varies between 40–70 and 10–20 m in the middle and lateral regions, respectively. Our approach yields new information on the structural conditions of rock cliff and column instabilities at depth, which is not easily obtained by other means but is imperative for change detection monitoring and improved hazard assessments.

58 GEOSCIENCES↗

Vertical Drop Testing and Analysis of the WASP Helicopter Skid Gear

Human occupant modeling and injury risk assessment have been identified as areas of research for improved prediction of rotorcraft crashworthiness within the NASA Aeronautics Program's Subsonic Rotary Wing Project. As part of this effort, an experimental program was conducted to assess the impact performance of a skid gear for use on the WASP kit-built helicopter, which is marketed by HeloWerks, Inc. of Hampton, Virginia. Test data from a drop test at an impact velocity of 8.4 feet-per-second were used to assess a finite element model of the skid gear test article. This assessment included human occupant analytic models developed for execution in LS-DYNA. The test article consisted of an aluminum skid gear mounted beneath a steel plate. A seating platform was attached to the upper surface of the steel plate, and two 95th percentile Hybrid III male Aerospace Anthropomorphic Test Devices (ATDs) were seated on the platform and secured using a four-point restraint system. The goal of the test-analysis correlation is to further the understanding of LS-DYNA ATD occupant models and responses in the vertical (or spinal) direction. By correlating human occupant experimental test data for a purely vertical impact with the LS-DYNA occupant responses, improved confidence in the use of these tools and better understanding of the limitations of the automotive-based occupant models for aerospace application can begin to be developed.

Fuchs, Yvonne T.↗