Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic strain engineering”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Anomalous elastic behavior of tantalum at high pressures: Experimental and theoretical studies

Elastic wave velocities of polycrystalline tantalum have been measured up to 13.6 GPa at room temperature using ultrasonic interferometry technique in a multi-anvil apparatus. The bulk and shear moduli, as well as their pressure derivatives are obtained using Eulerian 3rd order finite strain equations, yielding K S0 = 193.9(29) GPa, G 0 = 69.3(10) GPa, K S0' = 3.18(5) and G 0' = 0.88(1). Here the Steinberg-Guinan yield stress model is examined using the current shear modulus and its pressure derivative, where a slightly lower A 0 = 0.0127 GPa –1 compared to previous studies is obtained. First-principles calculations on the elasticity of tantalum have also been performed up to 200 GPa. Anomalous softening of the elastic shear constant C 44 and shear wave velocity V S is observed between 80– 180 GPa, whereas the anisotropy shows a softening-stiffening behavior with increasing pressure and reaches a minimum at around 150 GPa. Such anomalies might commonly originate from the nesting of Fermi surface for Group VB metals.

36 MATERIALS SCIENCE↗
NASA NTRS1998

Experimental Techniques Verified for Determining Yield and Flow Surfaces

Structural components in aircraft engines are subjected to multiaxial loads when in service. For such components, life prediction methodologies are dependent on the accuracy of the constitutive models that determine the elastic and inelastic portions of a loading cycle. A threshold surface (such as a yield surface) is customarily used to differentiate between reversible and irreversible flow. For elastoplastic materials, a yield surface can be used to delimit the elastic region in a given stress space. The concept of a yield surface is central to the mathematical formulation of a classical plasticity theory, but at elevated temperatures, material response can be highly time dependent. Thus, viscoplastic theories have been developed to account for this time dependency. Since the key to many of these theories is experimental validation, the objective of this work (refs. 1 and 2) at the NASA Lewis Research Center was to verify that current laboratory techniques and equipment are sufficient to determine flow surfaces at elevated temperatures. By probing many times in the axial-torsional stress space, we could define the yield and flow surfaces. A small offset definition of yield (10 me) was used to delineate the boundary between reversible and irreversible behavior so that the material state remained essentially unchanged and multiple probes could be done on the same specimen. The strain was measured with an off-the-shelf multiaxial extensometer that could measure the axial and torsional strains over a wide range of temperatures. The accuracy and resolution of this extensometer was verified by comparing its data with strain gauge data at room temperature. The extensometer was found to have sufficient resolution for these experiments. In addition, the amount of crosstalk (i.e., the accumulation of apparent strain in one direction when strain in the other direction is applied) was found to be negligible. Tubular specimens were induction heated to determine the flow surfaces at elevated temperatures. The heating system induced a large amount of noise in the data. By reducing thermal fluctuations and using appropriate data averaging schemes, we could render the noise inconsequential. Thus, accurate and reproducible flow surfaces (see the figure) could be obtained.

Lerch, Brad A.↗
DOE OSTI2020

Primal interface debonding formulation for finite strain isotropic plasticity

In this work, a framework is developed for modeling ductile damage of nonlinear materials whose plastic deformation is characterized using rate independent classical plasticity. This method relies on the assumption that the free energy can be decomposed into elastic, plastic and damage parts. A thermodynamically consistent method is derived which satisfies the second law of thermodynamics in the Clausius–Duhem inequality form. The dissipation associated with plasticity takes place in the domain only, while damage dissipation is localized to the interface. The method is developed using Variational Multiscale ideas to obtain definitions of the interface fluxes within a primal formulation analogous to the Discontinuous Galerkin method, which ensures weakly vanishing interface gap prior to reaching a damage initiation criterion. The local nonlinear problem to calculate both plastic deformation gradient and damage variable follows an incremental approach similar to classical plasticity return mapping algorithm. This elastoplastic damage formulation is developed for material undergoing finite strain, and it naturally accommodates a trapezoidal traction separation law (TSL) whose shape can be varied to model either ductile interface behavior or brittle interface behavior. The formulation's performance is assessed through modeling a patch test and a compact tension specimen.

42 ENGINEERING↗
DOE OSTI2020

Localized Strain Analysis of Ce- and Mg-Treated Cast Iron under Uniaxial Compression

Cast iron exhibits a wide range of mechanical properties, depending on its microstructural features. The microstructure of cast iron consists of several microconstituents with different elastic-plastic behavior, making the strain non-uniform across the bulk material. To understand the individual effects of these microconstituents on the overall mechanical behavior, local strain analysis using digital image correlation analysis was carried out. Samples with two different compositions (varying cerium, magnesium and silicon) were processed at different solidification velocities in a Bridgman furnace. Sections of the directionally solidified samples were loaded under uniaxial compression to measure global and local strain behavior. Despite the variability of the microstructure, the stress–strain curves obtained by digital image correlation (DIC) were found to react in a well-controlled way to changes in solidification velocity. It was observed that high-strain failure (greater than 15%) was accompanied by local straining of the softer ferritic phase, but during low-strain failure, local straining was not prominent. Higher nodularities, due to higher solidification velocities, raised the compressive strength without affecting the toughness significantly. Higher percentages of carbides led to higher compressive strengths with corresponding losses in ductility. The continuity of the matrix was also found to play an important role in the behavior during compression.

ductile failure↗
NASA NTRS2024

Scalability in SLA lattice through lattice orientation and hybrid frame and plate architectures

Lightweighting has been a key goal for engineers and designers, with lattice structures widely explored as the building blocks for structural components. Cellular structure-inspired lattice truss frame designs made of struts, have been extensively studied. All plate-based lattice designs have superior mechanical performance. Limits in scalability occur from formation of closed pockets limiting uv curing in processes like stereolithography (SLA). We examine hybrid frame and plate body-centered cubic-simply cubic (BCC-SC) lattices under compression. Unit cells and scaled up lattice exhibit an increase in yield stress and modulus with the addition of plates. The loading direction on the hybrid frame and plate unit cells affected the magnitude of improvement. Simulations and measurements indicated that the optimal lightweight lattice was determined to be when two plates were placed opposite each other with plates buttressing the struts, inhibiting buckling of the struts aligned with the loading direction. This lattice resulted in 63% improvement in specific modulus, a 137% improvement in specific yield point, and a 360% improvement in specific energy absorption (SEA) and the scaled up 4 × 4 × 4 scaled-up structures, showed a 107%, 148%, and 297% in specific modulus, specific yield point, and SEA, respectively A combined stretching-bending behavior was identified in optimal orientations reflecting the delayed buckling mechanism paired to a rising stress-strain curve past the elastic yield indicating bending resistance. The mass moment of inertia was found to be a key parameter correlating optimum orientation for the same number of plates added to the BCC-SC frame.

Mahan Ghosh↗
DOE OSTI2026

Compressive Response and Energy Absorption of Additively Manufactured Elastomers with Varied Simple Cubic Architectures

Additive manufacturing, and particularly the vat photopolymerization process, enables the fabrication of complex geometries at high resolution and small length scales, making it well-suited for fabricating cellular structures (e.g., foams and lattices). Among these, elastomeric cellular structures are of growing interest due to their tunable compliance and energy dissipation. However, comprehensive data on the compressive behavior of these structures remains limited, especially for investigating the structure-property effects from changing the density and distribution of material within the cellular structure. This study explores how the mechanical response of polyurethane-based simple cubic structures changes when varying volume fraction, unit cell length, and unit cell patterning, which have not been systematically investigated previously in additively manufactured elastomers. Increasing volume fraction from 10% to 50% yielded significant changes in compressive stress–strain performance (decreasing strain at 0.5 MPa by 41.6% and increasing energy absorption density by 3962.5%). Although changing the unit cell length between 2.5 and 7 mm in ~30 mm parts did not result in statistically different stress–strain responses, modifying the configuration of struts of different thicknesses across designs with 30% volume fraction altered the stress–strain behavior (differences of 12.5% in strain at 0.5 MPa and 109.4% for energy absorption density). Power law relationships were developed to understand the interactions between volume fraction, unit cell length, and elastic modulus, and experimental data showed strong fits (R 2 > 0.91). These findings enhance the understanding of how multiple structural design aspects influence the performance of elastomeric cellular materials, providing a foundation for informing strategic design of tailorable materials for diverse mechanical applications.

36 MATERIALS SCIENCE↗
NASA NTRS2004

Validation of Design and Analysis Techniques of Tailored Composite Structures

Aeroelasticity is the relationship between the elasticity of an aircraft structure and its aerodynamics. This relationship can cause instabilities such as flutter in a wing. Engineers have long studied aeroelasticity to ensure such instabilities do not become a problem within normal operating conditions. In recent decades structural tailoring has been used to take advantage of aeroelasticity. It is possible to tailor an aircraft structure to respond favorably to multiple different flight regimes such as takeoff, landing, cruise, 2-g pull up, etc. Structures can be designed so that these responses provide an aerodynamic advantage. This research investigates the ability to design and analyze tailored structures made from filamentary composites. Specifically the accuracy of tailored composite analysis must be verified if this design technique is to become feasible. To pursue this idea, a validation experiment has been performed on a small-scale filamentary composite wing box. The box is tailored such that its cover panels induce a global bend-twist coupling under an applied load. Two types of analysis were chosen for the experiment. The first is a closed form analysis based on a theoretical model of a single cell tailored box beam and the second is a finite element analysis. The predicted results are compared with the measured data to validate the analyses. The comparison of results show that the finite element analysis is capable of predicting displacements and strains to within 10% on the small-scale structure. The closed form code is consistently able to predict the wing box bending to 25% of the measured value. This error is expected due to simplifying assumptions in the closed form analysis. Differences between the closed form code representation and the wing box specimen caused large errors in the twist prediction. The closed form analysis prediction of twist has not been validated from this test.

Jegley, Dawn C.↗
DOE OSTI2023

Flexoelectricity and New Phenomena

It is easy to miss the scientific implications of our recent work on Triboelectricity. Everyone knows that rubbing and contact can produce static electricity; less appreciated is that the thermodynamic driver has been an open question since static electricity was first observed by Thales of Miletus around 585 BC. People think they understand it, for instance one common explanation that can still be found in the current literature is that differences in the work function drives charge transfer, often called the Volta-Helmholtz hypothesis. As summarized in 1967 by Harper, this fails to explain many experimental observations, for instance that charging can occur when two pieces of the same material are rubbed against each other. We are the first to place triboelectricity on a solid foundation rooted in quantum mechanics – the flexoelectric effect. We were able to explain a significant number of previously unexplained phenomena: 1) Bipolar tribocurrents associated with stick-slip, due to the change in sign of the strain gradients. 2) A one-third power scaling of tribocurrents with indentation force. 3) Tribocharging when two identical materials are used, these being due to local variations in the asperities so there are usually local potential differences. 4) Inhomogeneous charging of insulators, related to the statistical nature of asperities. 5) An experimentally observed reversal in the sign of charge transfer for negative and positive curvature, which is related to a change in the sign of the strain gradient. Exploiting our DOE prior funded work on flexoelectricity, we obtained semi-quantitative matching to existing experimental measurements of the surface charge in triboelectric experiments. The work has been well received in the literature. The work has been the focus of a number of popular science press articles, and also formed the basis for a Podcast for children 6-10 “The Rise and Fall of Static Man” posted in December 2019 by NPR as part of their “Wow in the World” series. I was also briefly interviewed by the Chicago PBS station in January 2020. This work are significant for a wide range of energy applications; to quote from an independent source: "Triboelectric power has plenty of potential, says Wenzhuo Wu, an assistant professor of engineering at Purdue. If the basics of static electricity are better understood, we could maximize the efficiency of wind or wave power generators, Wu says. The body's own movement could be used to power internal medical devices. Imagine being able to create a roof shaped to harness the power of a raindrop — the friction of the rain passing over the surface — to generate triboelectricity, powering the building below it." This is the start of new science, some of which we already partially understand such as the role of band bending in charge transfer. We need to understand charge transfer combining elasticity, quantum mechanics, band bending and defect states. These directly involve several of the DOE Grand Challenges: "How do we control material processes at the level of electrons? How do remarkable properties of matter emerge from complex correlations of the atomic or electronic constituents and how can we control these properties? How do we characterize and control matter away—especially very far away— from equilibrium?" I will argue that this work truly falls into the class of disruptive science; it is not just a simple extension, linear science. Not everyone will accept the approach. Since we explain far more about triboelectricity than anyone before, the preponderance of evidence supports the model. The feedback I have received is that many agree with the work, to quote: "The model makes sense, says Michael McAlpine, a professor of engineering at the University of Minnesota. "It's such a simple explanation, I was surprised I didn't put my finger on that," McAlpine says." The proposal received strong reviews. It was also publicized on the Department of Energy Web Site.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Bending-induced isostructural transitions in ultrathin layers of van der Waals ferrielectrics

Using Landau-Ginzburg-Devonshire (LGD) phenomenological approach we analyze the bending-induced re-distribution of electric polarization and field, elastic stresses and strains inside ultrathin layers of van der Waals ferrielectrics. We consider a CuInP2S6 (CIPS) thin layer with fixed edges and suspended central part, the bending of which is induced by external forces. The unique aspect of CIPS is the existence of two ferrielectric states, FI1 and FI2, corresponding to big and small polarization values, which arise due to the specific four-well potential of the eighth-order LGD functional. When the CIPS layer is flat, the single-domain FI1 state is stable in the central part of the layer, and the FI2 states are stable near the fixed edges. With an increase of the layer bending below the critical value, the sizes of the FI2 states near the fixed edges decreases, and the size of the FI1 region increases. When the bending exceeds the critical value, the edge FI2 states disappear being substituted by the FI1 state, but they appear abruptly near the inflection regions and expand as the bending increases. The bending-induced isostructural FI1-FI2 transition is specific for the bended van der Waals ferrielectrics described by the eighth (or higher) order LGD functional with consideration of linear and nonlinear electrostriction couplings. Furthermore, the isostructural transition, which is revealed in the vicinity of room temperature, can significantly reduce the coercive voltage of ferroelectric polarization reversal in CIPS nanoflakes, allowing for the curvature-engineering control of various flexible nanodevices.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Revealing the unusual rate-dependent mechanical behaviors of nematic liquid crystal elastomers

Liquid crystal elastomers (LCEs) exhibit unique mechanical properties of soft elasticity and enhanced energy dissipation with rate dependency. They are potentially transformative materials for applications in mechanical impact mitigation and vibration isolation. However, previous studies have primarily focused on the mechanics of LCEs under equilibrium and quasistatic loading conditions. Critical knowledge gaps exist in understanding their rate-dependent behaviors, which are a complex mixture of traditional network viscoelasticity and the soft elastic behaviors with changes in the mesogen orientation and order parameter. Together, these inelastic mechanisms lead to unusual rate-dependent energy absorption responses of LCEs. In this work, we developed a viscoelastic constitutive theory for monodomain nematic LCEs to investigate how multiple underlying sources of inelasticity manifest in the rate-dependent and dissipative behaviors of monodomain LCEs. Here, the theoretical modeling framework combines the neo-classical network theory with evolution rules for the mesogen orientation and order parameter with conventional viscoelasticity. The model is calibrated with uniaxial tension and compression data spanning six decades of strain rates. The established 3D constitutive model enables general loading predictions taking the initial mesogen orientation and order parameter as inputs. Additionally, parametric studies were performed to further understand the rate dependence of monodomain LCEs in relation to their energy absorption characteristics. Based on the parametric studies, particularly loading scenarios are identified as conditions where LCEs outperform conventional elastomers regarding energy absorption.

42 ENGINEERING↗
DOE OSTI2021

Nonlinear convergence in contact mechanics: Immersed boundary finite volume

In this report we present an immersed boundary finite volume (IBM) method for simulating quasistatic contact mechanics of linearly elastic domains at small strains. In IBM, all external boundaries and internal contacts of an object are represented by embedded surfaces inside a Cartesian mesh, which need not conform to the grid lines. The contact constraints consist of the non-penetrability condition and Coulomb’s friction law, which are discretized using special interpolation stencils and enforced via penalty parameters. The resulting nonlinear system depends on displacement unknowns only. To solve it, we use the Newton method but find that it diverges frequently. To understand the divergence pattern, we analyze a simplified 2-cell problem and show that the global convergence of Newton cannot be ensured for any choice of penalty parameters. We thus propose a modified Newton solver, which guarantees convergence for the 2-cell problem and is numerically verified to converge for all the challenging simulations considered herein. While both 1 st - and 2 nd -order variants of IBM, in displacement unknowns, are proposed, the modified Newton solver applies only to the 1 st -order variant.

42 ENGINEERING↗
DOE OSTI2022

Non-destructive depth-resolved characterization of residual strain fields in high electron mobility transistors using differential aperture x-ray microscopy

Localized residual stress and elastic strain concentrations in microelectronic devices often affect the electronic performance, resistance to thermomechanical damage, and, likely, radiation tolerance. A primary challenge for the characterization of these concentrations is that they exist over sub-μm length-scales, precluding their characterization by more traditional residual stress measurement techniques. Here, in this study, we demonstrate the use of synchrotron x-ray-based differential aperture x-ray microscopy (DAXM) as a viable, non-destructive means to characterize these stress and strain concentrations in a depth-resolved manner. DAXM is used to map two-dimensional strain fields between the source and the drain in a gallium nitride (GaN) layer within high electron mobility transistors (HEMTs) with sub-μm spatial resolution. Strain fields at various positions in both pristine and irradiated HEMT specimens are presented in addition to a preliminary stress analysis to estimate the distribution of various stress components within the GaN layer. γ-irradiation is found to significantly reduce the lattice plane spacing in the GaN along the sample normal direction, which is attributed to radiation damage in transistor components bonded to the GaN during irradiation.

42 ENGINEERING↗
DOE OSTI2023

Type 2 diabetes impairs annulus fibrosus fiber deformation and rotation under disc compression in the University of California Davis type 2 diabetes mellitus (UCD-T2DM) rat model

Abstract Understanding the biomechanical behavior of the intervertebral disc is crucial for studying disease mechanisms and developing tissue engineering strategies for managing disc degeneration. We used synchrotron small-angle X-ray scattering to investigate how changes to collagen behavior contribute to alterations in the disc’s ability to resist compression. Coccygeal motion segments from 6-month-old lean Sprague-Dawley rats ( n=7) and diabetic obese University of California Davis type 2 diabetes mellitus (UCD-T2DM) rats ( n=6, diabetic for 68±7 days) were compressed during simultaneous synchrotron scanning to measure collagen strain at the nanoscale (beamline 7.3.3 of the Advanced Light Source). After compression, the annulus fibrosus was assayed for nonenzymatic cross-links. In discs from lean rats, resistance to compression involved two main energy-dissipation mechanisms at the nanoscale: (1) rotation of the two groups of collagen fibrils forming the annulus fibrosus and (2) straightening (uncrimping) and stretching of the collagen fibrils. In discs from diabetic rats, both mechanisms were significantly impaired. Specifically, diabetes reduced fibril rotation by 31% and reduced collagen fibril strain by 30% (compared to lean discs). The stiffening of collagen fibrils in the discs from diabetic rats was consistent with a 31% higher concentration of nonenzymatic cross-links and with evidence of earlier onset plastic deformations such as fibril sliding and fibril–matrix delamination. These findings suggest that fibril reorientation, stretching, and straightening are key deformation mechanisms that facilitate whole-disc compression, and that type 2 diabetes impairs these efficient and low-energy elastic deformation mechanisms, thereby altering whole-disc behavior and inducing the earlier onset of plastic deformation.

Rosenberg, James L.↗
DOE OSTI2022

Modeling wave propagation in elastic solids via high-order accurate implicit-mesh discontinuous Galerkin methods

Here, a high-order accurate implicit-mesh discontinuous Galerkin framework for wave propagation in single-phase and bi-phase solids is presented. The framework belongs to the embedded-boundary techniques and its novelty regards the spatial discretization, which enables boundary and interface conditions to be enforced with high-order accuracy on curved embedded geometries. High-order accuracy is achieved via high-order quadrature rules for implicitly-defined domains and boundaries, whilst a cell-merging strategy addresses the presence of small cut cells. The framework is used to discretize the governing equations of elastodynamics, written using a first-order hyperbolic momentum-strain formulation, and an exact Riemann solver is employed to compute the numerical flux at the interface between dissimilar materials with general anisotropic properties. The space-discretized equations are then advanced in time using explicit high-order Runge–Kutta algorithms. Several two- and three-dimensional numerical tests including dynamic adaptive mesh refinement are presented to demonstrate the high-order accuracy and the capability of the method in the elastodynamic analysis of single- and bi-phases solids containing complex geometries.

42 ENGINEERING↗
DOE OSTI2021

Anisotropic elastoplastic response of double-porosity media

Here we present a continuum framework to simulate fluid flow through anisotropic elastoplastic media with double porosity. Two effective stress measures σ' and δ" emerge from the thermodynamic formulation, which are energy-conjugate to the elastic and plastic components of strain, respectively. Both effective stress measures can be expressed as a combination of the total Cauchy stress σ and the average pore pressure $\bar{p}$ in the two pore scales. In the effective stress for elasticity, $\bar{p}$ is scaled with a rank-2 Biot tensor, whereas the effective stress for plasticity follows the Terzaghi form in which $\bar{p}$ is scaled by the Kronecker delta. The Biot tensor and storage coefficients are derived as functions of elasticity parameters and porosities. A mixed finite element formulation is introduced to discretize the domain and solve initial boundary-value problems. A stabilization scheme is employed on equal-order interpolation for both displacement and pressure fields throughout the entire range of drainage responses. Numerical simulations reproduce the hydromechanical response of Opalinus shale in one-dimensional consolidation tests throughout the range of primary and secondary consolidation under different external loads. Numerical simulations of the consolidation of a rectangular domain subjected to a strip load demonstrate the efficacy of the proposed stabilization scheme, as well as illustrate the impacts of stress history, mass transfer, and different pore systems on the hydromechanical response.

42 ENGINEERING↗
NASA NTRS2009

Modeling and Simulation of Variable Mass, Flexible Structures

The advent of the new Ares I launch vehicle has highlighted the need for advanced dynamic analysis tools for variable mass, flexible structures. This system is composed of interconnected flexible stages or components undergoing rapid mass depletion through the consumption of solid or liquid propellant. In addition to large rigid body configuration changes, the system simultaneously experiences elastic deformations. In most applications, the elastic deformations are compatible with linear strain-displacement relationships and are typically modeled using the assumed modes technique. The deformation of the system is approximated through the linear combination of the products of spatial shape functions and generalized time coordinates. Spatial shape functions are traditionally composed of normal mode shapes of the system or even constraint modes and static deformations derived from finite element models of the system. Equations of motion for systems undergoing coupled large rigid body motion and elastic deformation have previously been derived through a number of techniques [1]. However, in these derivations, the mode shapes or spatial shape functions of the system components were considered constant. But with the Ares I vehicle, the structural characteristics of the system are changing with the mass of the system. Previous approaches to solving this problem involve periodic updates to the spatial shape functions or interpolation between shape functions based on system mass or elapsed mission time. These solutions often introduce misleading or even unstable numerical transients into the system. Plus, interpolation on a shape function is not intuitive. This paper presents an approach in which the shape functions are held constant and operate on the changing mass and stiffness matrices of the vehicle components. Each vehicle stage or component finite element model is broken into dry structure and propellant models. A library of propellant models is used to describe the distribution of mass in the fuel tank or Solid Rocket Booster (SRB) case for various propellant levels. Based on the mass consumed by the liquid engine or SRB, the appropriate propellant model is coupled with the dry structure model for the stage. Then using vehicle configuration data, the integrated vehicle model is assembled and operated on by the constant system shape functions. The system mode shapes and frequencies can then be computed from the resulting generalized mass and stiffness matrices for that mass configuration. The rigid body mass properties of the vehicle are derived from the integrated vehicle model. The coupling terms between the vehicle rigid body motion and elastic deformation are also updated from the constant system shape functions and the integrated vehicle model. This approach was first used to analyze variable mass spinning beams and then prototyped into a generic dynamics simulation engine. The resulting code was tested against Crew Launch Vehicle (CLV-)class problems worked in the TREETOPS simulation package and by Wilson [2]. The Ares I System Integration Laboratory (SIL) is currently being developed at the Marshall Space Flight Center (MSFC) to test vehicle avionics hardware and software in a hardware-in-the-loop (HWIL) environment and certify that the integrated system is prepared for flight. The Ares I SIL utilizes the Ares Real-Time Environment for Modeling, Integration, and Simulation (ARTEMIS) tool to simulate the launch vehicle and stimulate avionics hardware. Due to the presence of vehicle control system filters and the thrust oscillation suppression system, which are tuned to the structural characteristics of the vehicle, ARTEMIS must incorporate accurate structural models of the Ares I launch vehicle. The ARTEMIS core dynamics simulation models the highly coupled nature of the vehicle flexible body dynamics, propellant slosh, and vehicle nozzle inertia effects combined with mass and flexible body properties that vary significant with time during the flight. All forces that act on the vehicle during flight must be simulated, including deflected engine thrust force, spatially distributed aerodynamic forces, gravity, and reaction control jet thrust forces. These forces are used to excite an integrated flexible vehicle, slosh, and nozzle dynamics model for the vehicle stack that simulates large rigid body translations and rotations along with small elastic deformations. Highly effective matrix math operations on a distributed, threaded high-performance simulation node allow ARTEMIS to retain up to 30 modes of flex for real-time simulation. Stage elements that separate from the stack during flight are propagated as independent rigid six degrees of freedom (6DOF) bodies. This paper will present the formulation of the resulting equations of motion, solutions to example problems, and describe the resulting dynamics simulation engine within ARTEMIS.

Tobbe, Patrick A.↗