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 829 records · Page 46

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.↗

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↗

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.↗

Three-dimensional numerical modeling of indium phosphide Point-Contact Solar Cells

The Point-Contact Solar Cell (PCSC) geometry has proven very effective for silicon cells. To date, it has not been implemented in III-V materials. In addition, modeling such a geometry is very difficult because of its three-dimensional nature. We have developed a three-dimensional finite element modeling code (FIESTA ROC). In this paper, we present results from a three-dimensional modeling study of InP point-contact solar cells.

Clark, Ralph O.↗

Influence of Back-Up Bearings and Support Structure Dynamics on the Behavior of Rotors With Active Supports

This report presents a synopsis of the research work. Specific accomplishments are itemized below: (1) Experimental facilities have been developed. This includes a magnetic bearing test rig and an auxiliary bearing test rig. In addition, components have been designed, constructed, and tested for use with a rotordynamics test rig located at NASA Lewis Research Center. (2) A study of the rotordynamics of an auxiliary bearing supported T-501 engine model was performed. (3) An experimental/simulation study of auxiliary bearing rotordynamics has been performed. (4) A rotordynamical model for a magnetic bearing supported rotor system, including auxiliary bearing effects has been developed and simulation studies performed.(5) A finite element model for a foil bearing has been developed and studies of a rotor supported by foil bearings have been performed. (6) Two students affiliated with this project have graduated with M.S. degrees.

Flowers, George T.↗

Simulation of Hypervelocity Impact on Aluminum-Nextel-Kevlar Orbital Debris Shields

An improved hybrid particle-finite element method has been developed for hypervelocity impact simulation. The method combines the general contact-impact capabilities of particle codes with the true Lagrangian kinematics of large strain finite element formulations. Unlike some alternative schemes which couple Lagrangian finite element models with smooth particle hydrodynamics, the present formulation makes no use of slidelines or penalty forces. The method has been implemented in a parallel, three dimensional computer code. Simulations of three dimensional orbital debris impact problems using this parallel hybrid particle-finite element code, show good agreement with experiment and good speedup in parallel computation. The simulations included single and multi-plate shields as well as aluminum and composite shielding materials. at an impact velocity of eleven kilometers per second.

Fahrenthold, Eric P.↗

Development and Application of Benchmark Examples for Mixed-Mode I/II Quasi-Static Delamination Propagation Predictions

The development of benchmark examples for quasi-static delamination propagation prediction is presented. The example is based on a finite element model of the Mixed-Mode Bending (MMB) specimen for 50% mode II. The benchmarking is demonstrated for Abaqus/Standard, however, the example is independent of the analysis software used and allows the assessment of the automated delamination propagation prediction capability in commercial finite element codes based on the virtual crack closure technique (VCCT). First, a quasi-static benchmark example was created for the specimen. Second, starting from an initially straight front, the delamination was allowed to propagate under quasi-static loading. Third, the load-displacement as well as delamination length versus applied load/displacement relationships from a propagation analysis and the benchmark results were compared, and good agreement could be achieved by selecting the appropriate input parameters. The benchmarking procedure proved valuable by highlighting the issues associated with choosing the input parameters of the particular implementation. Overall, the results are encouraging, but further assessment for mixed-mode delamination fatigue onset and growth is required.

Krueger, Ronald↗

Three-Dimensional High Fidelity Progressive Failure Damage Modeling of NCF Composites

Performance prediction of off-axis laminates is of significant interest in designing composite structures for energy absorption. Phenomenological models available in most of the commercial programs, where the fiber and resin properties are smeared, are very efficient for large scale structural analysis, but lack the ability to model the complex nonlinear behavior of the resin and fail to capture the complex load transfer mechanisms between the fiber and the resin matrix. On the other hand, high fidelity mesoscale models, where the fiber tows and matrix regions are explicitly modeled, have the ability to account for the complex behavior in each of the constituents of the composite. However, creating a finite element model of a larger scale composite component could be very time consuming and computationally very expensive. In the present study, a three-dimensional mesoscale model of non-crimp composite laminates was developed for various laminate schemes. The resin material was modeled as an elastic-plastic material with nonlinear hardening. The fiber tows were modeled with an orthotropic material model with brittle failure. In parallel, new stress based failure criteria combined with several damage evolution laws for matrix stresses were proposed for a phenomenological model. The results from both the mesoscale and phenomenological models were compared with the experiments for a variety of off-axis laminates.

Aitharaju, Venkat↗

Computationally Efficient Beam Elements for Accurate Stresses in Sandwich Laminates and Laminated Composites with Delaminations

Laminated composites are prone to delamination failure due to the lack of reinforcement through the thickness. Therefore, during the design process the initiation and propagation of delaminations should be accounted for as early as possible. This paper presents computationally efficient nine degree-of-freedom (dof) and eight-dof shear locking-free beam elements using the mixed form of the refined zigzag theory (RZT(sup(m))). The corresponding nine-dof and eight-dof elements use the anisoparametric and constrained anisoparametric interpolation schemes, respectively, to eliminate shear locking in slender beams. The advantage of the present element over previous RZT beam elements is that no post-processing is required to accurately model the transverse shear stress while maintaining the computational efficiency of a low-order beam element. Comparisons with high-fidelity finite element models and three-dimensional elasticity solutions show that the elements can robustly and accurately predict the displacement field, axial stress and transverse shear stress through the thickness of a sandwich beam or a composite laminate with an embedded delamination. In fact, the accuracy and computational efficiency of predicting stresses in laminates with embedded delaminations make the present elements attractive choices for RZT-based delamination initiation and propagation methodologies available in the literature.

R. M. J. Groh↗

Transient Response Dynamic Module Modifications to Include Static and Kinetic Friction Effects

A methodology that supports forced transient response dynamic solutions when both static and kinetic friction effects are included in a structural system model is described. Modifications that support this type of nonlinear transient response solution are summarized for the transient response dynamics (TRD) NASTRAN module. An overview of specific modifications for the NASTRAN processing subroutines, INITL, TRD1C, and TRD1D, are described with further details regarding inspection of nonlinear input definitions to define the type of nonlinear solution required, along with additional initialization requirements and specific calculation subroutines to successfully solve the transient response problem. The extension of the basic NASTRAN nonlinear methodology is presented through several stages of development to the point where constraint equations and residual flexibility effects are introduced into the finite difference Newmark-Beta recurrsion formulas. Particular emphasis is placed on cost effective solutions for large finite element models such as the Space Shuttle with friction degrees of freedom between the orbiter and payloads mounted in the cargo bay. An alteration to the dynamic finite difference equations of motion is discussed, which allows one to include friction effects at reasonable cost for large structural systems such as the Space Shuttle. Data are presented to indicate the possible impact of transient friction loads to the payload designer for the Space Shuttle. Transient response solution data are also included, which compare solutions without friction forces and those with friction forces for payloads mounted in the Space Shuttle cargo bay. These data indicate that payload components can be sensitive to friction induced loads.

Misel, J. E.↗

Structural Damage Detection using Modal Test Data

This study presents a methodology for updating the finite element model of a structure for damage detection purposes using an incomplete set of experimentally obtained modal frequencies and modeshapes.

Damage↗