Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic forces”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

An approximation method for configuration optimization of trusses

Two- and three-dimensional elastic trusses are designed for minimum weight by varying the areas of the members and the location of the joints. Constraints on member stresses and Euler buckling are imposed and multiple static loading conditions are considered. The method presented here utilizes an approximate structural analysis based on first order Taylor series expansions of the member forces. A numerical optimizer minimizes the weight of the truss using information from the approximate structural analysis. Comparisons with results from other methods are made. It is shown that the method of forming an approximate structural analysis based on linearized member forces leads to a highly efficient method of truss configuration optimization.

Hansen, Scott R.↗
NASA NTRS2004

Characterization of Composite Fan Case Resins

The majority of commercial turbine engines that power today s aircraft use a large fan driven by the engine core to generate thrust which dramatically increases the engine s efficiency. However, if one of these fan blades fails during flight, it becomes high energy shrapnel, potentially impacting the engine or puncturing the aircraft itself and thus risking the lives of passengers. To solve this problem, the fan case must be capable of containing a fan blade should it break off during flight. Currently, all commercial fan cases are made of either just a thick metal barrier or a thinner metal wall surrounded by Kevlar-an ultra strong fiber that elastically catches the blade. My summer 2004 project was to characterize the resins for a composite fan case that will be lighter and more efficient than the current metal. The composite fan case is created by braiding carbon fibers and injecting a polymer resin into the braid. The resin holds the fibers together, so at first using the strongest polymer appears to logically lead to the strongest fan case. Unfortunately, the stronger polymers are too viscous when melted. This makes the manufacturing process more difficult because the polymer does not flow as freely through the braid, and the final product is less dense. With all of this in mind, it is important to remember that the strength of the polymer is still imperative; the case must still contain blades with high impact energy. The research identified which polymer had the right balance of properties, including ease of fabrication, toughness, and ability to transfer the load to the carbon fibers. Resin deformation was studied to better understand the composite response during high speed impact. My role in this research was the testing of polymers using dynamic mechanical analysis and tensile, compression, and torsion testing. Dynamic mechanical analysis examines the response of materials under cyclic loading. Two techniques were used for dynamic mechanical analysis. The ARES Instrument analyzed the material through torsion. The second machine, TA Instruments apparatus, applied a bending force to the specimen. These experiments were used to explore the effects of temperature and strain rate on the stiffness and strength of the resins. The two different types of loading allowed us to verify our results. An axial-torsional load frame, manufactured by MTS Systems, Inc., was used to conduct the tensile, compression, and torsional testing. These tests were used to determine the stress-strain curves for the resins. The elastic and plastic deformation data was provided to another team member for characterization of high fidelity material property predictions. This information was useful in having a better understanding of the polymers so that the fan cases could be as sturdy as possible. Deformation studies are the foundation for the computational modeling that provides the structural design of a composite engine case as well as detailed analysis of the blade impact event.

Dvoracek, Charlene M.↗
NASA NTRS1974

Transverse vibrations of non-uniform beams

The free vibrations of elastic beams with nonuniform characteristics are investigated theoretically by a new method. The new method is seen to combine the advantages of a finite element approach and of a Rayleigh-Ritz analysis. Comparison with the known analytical results for uniform beams shows good convergence of the method for natural frequencies and modes. For internal shear forces and bending moments, the rate of convergence is less rapid. Results from experiments conducted with a cantilevered helicopter blade with strong nonuniformities and also from alternative theoretical methods, indicate that the theory adequately predicts natural frequencies and mode shapes. General guidelines for efficient use of the method are presented.

Klein, L.↗
NASA NTRS2003

Aeroelastic System Development Using Proper Orthogonal Decomposition and Volterra Theory

This research combines Volterra theory and proper orthogonal decomposition (POD) into a hybrid methodology for reduced-order modeling of aeroelastic systems. The out-come of the method is a set of linear ordinary differential equations (ODEs) describing the modal amplitudes associated with both the structural modes and the POD basis functions for the uid. For this research, the structural modes are sine waves of varying frequency, and the Volterra-POD approach is applied to the fluid dynamics equations. The structural modes are treated as forcing terms which are impulsed as part of the uid model realization. Using this approach, structural and uid operators are coupled into a single aeroelastic operator. This coupling converts a free boundary uid problem into an initial value problem, while preserving the parameter (or parameters) of interest for sensitivity analysis. The approach is applied to an elastic panel in supersonic cross ow. The hybrid Volterra-POD approach provides a low-order uid model in state-space form. The linear uid model is tightly coupled with a nonlinear panel model using an implicit integration scheme. The resulting aeroelastic model provides correct limit-cycle oscillation prediction over a wide range of panel dynamic pressure values. Time integration of the reduced-order aeroelastic model is four orders of magnitude faster than the high-order solution procedure developed for this research using traditional uid and structural solvers.

Lucia, David J.↗
NASA NTRS2002

Investigation of Oxidation Profile in PMR-15 Polyimide using Atomic Microscope (AFM)

Nanoindentation measurements are made on thermosetting materials using cantiever deflection vs. piezoelectric scanner position behavior determined by AFM. The spring model is used to determine mechanical properties of materials. The generalized Sneddon's equation is utilized to calculate Young's moduli for thermosetting materials at ambient conditions. Our investigations show that the force-penetration depth curves during unloading in these materials can be described accurately by a power law relationship. The results show that the accuracy of the measurements can be controlled within 7%. The above method is used to study oxidation profiles in Pl\1R-15 polyimide. The thermo-mechanical profiles ofPNIR-15 indicate that the elastic modulus at the surface portion of the specimen is different from that at the interior of the material. It is also shown that there are two zones within the oxidized portion of the samples. Results confirm that the surface layer and the core material have substantially different properties.

Meador, Mary Ann B.↗
DOE OSTI2025

Assessment of critical flaw sizes and crack driving forces during additive manufacturing of metallic materials

Additive manufacturing (AM) of complex engineering components is often plagued by a high susceptibility to cracking, particularly in high-strength metallic materials. While alloy design efforts have made progress in mitigating solidification defects, there remains a need for mechanistic guidelines to predict susceptibility to solid-state cracking. To address this gap, driving forces for the growth of melt pool cracks are calculated across a wide range of alloys using an efficient computational framework. Calculations are coupled with rapid single track laser experiments to elucidate trends in cracking from laser melting. The analyses conducted here highlight the important role of material properties in susceptibility to cracking, notably fracture toughness and elastic modulus. An important finding is that residual stresses that are limited in magnitude to the yield stress of the material are likely insufficient to drive cracking during cooling. Furthermore, the implications of these results are discussed in the context of alloy design for AM and residual stress accumulation during AM.

36 MATERIALS SCIENCE↗
NASA NTRS2005

Surface Collisions Involving Particles and Moisture (SCIP'M)

Experiments were performed on the collision of a solid sphere with a nearly horizontal flat surface covered with a thin layer of viscous liquid. High-speed collisions were obtained by dropping the ball onto the surface from various heights, using gravitational acceleration. Low-speed collisions were obtained using pendulums with long strings or by launching the balls at low velocities in the reduced-gravity environment of parabolic flight. The sphere bounces only when the impact velocity exceeds a critical value. The coefficient of restitution (ratio of rebound velocity to impact velocity) increases with increasing impact velocity above the critical value, indicating the increasing relative importance of elastic deformation to viscous dissipation. The critical impact velocity increases, and the coefficient of restitution decreases, with increasing viscosity or thickness of the liquid layer and with decreasing density or size of the sphere. The ratio of the wet and dry coefficients is expressed as a function of the Stokes number (ratio of particle inertia and viscous forces), showing good agreement between theory and experiment. Similar experiments were performed with the flat surface inclined at various angles to the approaching sphere. A modified Stokes number, which is a measure of the ratio of inertia of the sphere in the normal direction to the viscous forces exerted by the fluid layer, was used for the analysis of oblique collisions. Even for these oblique collisions, it was found that no rebound of the ball was observed below a certain critical Stokes number. The coefficient of normal restitution, defined as a ratio of normal rebound velocity to normal approach velocity, was found to increase beyond the critical Stokes number and even out as it approaches the value for dry restitution at high Stokes numbers. It was also found that, for smooth spheres like steel, the normal restitution at the same modified Stokes number is independent of the angle of impact. The tangential coefficient of restitution, defined as the ratio of tangential rebound velocity to tangential approach velocity, is found to be nearly unity, except for very low approach velocities. Thus, as a first approximation, the theories that predict the coefficient of restitution for head-on wet collisions can be extended to predict the coefficient of normal restitution for oblique wet collisions. Additional experiments were performed with soft surfaces in which a porous cloth or sponge layer was placed over the hard, flat surface. In these experiments, the coefficient of restitution was found to decrease with increasing impact velocity, due to inelastic losses in the soft material. A model combining inelastic deformation and flow through porous media was developed to describe these findings.

Davis, Robert H.↗
NASA NTRS2015

An Efficient Solution Method for Multibody Systems with Loops Using Multiple Processors

This paper describes a multibody dynamics algorithm formulated for parallel implementation on multiprocessor computing platforms using the divide-and-conquer approach. The system of interest is a general topology of rigid and elastic articulated bodies with or without loops. The algorithm divides the multibody system into a number of smaller sets of bodies in chain or tree structures, called "branches" at convenient joints called "connection points", and uses an Order-N (O (N)) approach to formulate the dynamics of each branch in terms of the unknown spatial connection forces. The equations of motion for the branches, leaving the connection forces as unknowns, are implemented in separate processors in parallel for computational efficiency, and the equations for all the unknown connection forces are synthesized and solved in one or several processors. The performances of two implementations of this divide-and-conquer algorithm in multiple processors are compared with an existing method implemented on a single processor.

Multibody dynamics↗
DOE OSTI2021

A process to verify numerical models for seismic fluid-structure interaction in advanced reactor vessels

Seismic design and qualification of a liquid-filled advanced nuclear reactor will have to account for fluid-structure interaction (FSI). Interaction between the tank, internal components, and contained liquid will rely on analysis of numerical models that must be verified and validated. Here this study demonstrates a verification process for models of a base-supported cylindrical tank by comparing numerical predictions and analytical solutions. The numerical models are consistent with the assumptions made to derive analytical solutions, namely, either a rigid or a linear elastic tank, ideal fluid, and small-amplitude, unidirectional, horizontal inputs. One software platform is used to illustrate the process. Seismic FSI analysis is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Reported responses are those used for design, including hydrodynamic pressures on the tank wall, shear forces and moments at the tank base, and wave heights of the contained liquid. The accuracy of the numerical results is discussed. The numerical models are verified for calculating the pressures on the tank wall and reactions at its base. Accurate simulation of wave action is challenging for both solvers. Recommendations for modeling, code development, and steps for verification are provided. Although focused on reactor vessels and one software platform, the verification process described herein is broadly applicable to liquid-filled vessels and other finite element codes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

In Situ Nano-Indentation of a Gold Sub-Micrometric Particle Imaged by Multi-Wavelength Bragg Coherent X-ray Diffraction

The microstructure of a sub-micrometric gold crystal during nanoindentation is visualized by in situ multi-wavelength Bragg coherent X-ray diffraction imaging. The gold crystal is indented using a custom-built atomic force microscope. A band of deformation attributed to a shear band oriented along the (221) lattice plane is nucleated at the lower left corner of the crystal and propagates towards the crystal center with increasing applied mechanical load. After complete unloading, an almost strain-free and defect-free crystal is left behind, demonstrating a pseudo-elastic behavior that can only be studied by in situ imaging while it is invisible to ex situ examinations. The recovery is probably associated with reversible dislocations nucleation/annihilation at the side surface of the particle and at the particle-substrate interface, a behavior that has been predicted by atomistic simulations. The full recovery of the particle upon unloading sheds new light on extraordinary mechanical properties of metal nanoparticles obtained by solid-state dewetting.

36 MATERIALS SCIENCE↗
NASA NTRS1991

LAGEOS geodetic analysis-SL7.1

Laser ranging measurements to the LAGEOS satellite from 1976 through 1989 are related via geodetic and orbital theories to a variety of geodetic and geodynamic parameters. The SL7.1 analyses are explained of this data set including the estimation process for geodetic parameters such as Earth's gravitational constant (GM), those describing the Earth's elasticity properties (Love numbers), and the temporally varying geodetic parameters such as Earth's orientation (polar motion and Delta UT1) and tracking site horizontal tectonic motions. Descriptions of the reference systems, tectonic models, and adopted geodetic constants are provided; these are the framework within which the SL7.1 solution takes place. Estimates of temporal variations in non-conservative force parameters are included in these SL7.1 analyses as well as parameters describing the orbital states at monthly epochs. This information is useful in further refining models used to describe close-Earth satellite behavior. Estimates of intersite motions and individual tracking site motions computed through the network adjustment scheme are given. Tabulations of tracking site eccentricities, data summaries, estimated monthly orbital and force model parameters, polar motion, Earth rotation, and tracking station coordinate results are also provided.

Smith, D. E.↗
NASA NTRS2015

Gear Tooth Root Stresses of a Very Heavily Loaded Gear Pair-Case Study: Orbiter Body Flap Actuator Pinion and Ring Gear

The space shuttle orbiter's body flap actuator gearing was assessed as a case study of the stresses for very heavily loaded external-internal gear pairs (meshing pinion and ring gear). For many applications, using the high point of single tooth contact (HPSTC) to locate the position of the tooth force is adequate for assessing the maximum tooth root stress condition. But for aerospace gearing such an approach may be inadequate for assessing the stress condition while also simultaneously minimizing mass. In this work specialized contact analyses and finite element methods were used to study gear tooth stresses of body flap actuator gears. The analytical solutions considered the elastic deformations as an inherent part of the solutions. The ratio for the maximum tooth stresses using the HPSTC approach solutions relative to the contact analysis and finite element solutions were 1.40 for the ring gear and 1.28 for the pinion gear.

finite element method↗
NASA NTRS1985

Nonlinear quasi-fracture behavior of polymers

This paper deals with the nonlinear quasi-fracture behavior of polymers. Using the finite element method, the quasi-fracture model composed of an isolated craze opening with a nonlinear boundary envelope forces represented by oriented fibril domains is analyzed in simply-stressed rectangular plate. The stress distribution, the displacement profile, as well as the yield or plastic or flow zone adjacent to the craze tips are obtained. It is found that the variation of the stress distribution near the craze tip is not as sharp as that for an elastic crack. The initial load required for yielding or flowing is greater for the quasifracture as compared with that for a classical fracture. Under the same load the yield or flow zone is somewhat smaller as compared with that obtained on the basis of linear elastic-plastic fracture mechanics theory.

Sun, B. N.↗
NASA NTRS1993

Space Station/Orbiter berthing dynamics during an assembly flight

A large-angle, multi-body, dynamic modeling capability was developed to help validate numerical simulations of the dynamic motion and control forces which occur while berthing Space Station Freedom to the Shuttle Orbiter during early assembly flights. The paper describes the dynamics and control of the station, the attached Shuttle Remote Manipulator System, and the Orbiter during a maneuver from a gravity-gradient attitude to a torque equilibrium attitude using the station reaction control jets. The influence of the elastic behavior of the station and of the remote manipulator system on the attitude control of the station/Orbiter system during the maneuver is investigated. The flexibility of the station and the arm had only a minor influence on the attitude control of the system during the maneuver.

Cooper, Paul A.↗
NASA NTRS2001

Structural and Acoustic Numerical Modeling of a Curved Composite Honeycomb Panel

The finite and boundary element modeling of the curved section of a composite honeycomb aircraft fuselage sidewall was validated for both structural response and acoustic radiation. The curved panel was modeled in the pre-processor MSC/PATRAN. Geometry models of the curved panel were constructed based on the physical dimensions of the test article. Material properties were obtained from the panel manufacturer. Finite element models were developed to predict the modal parameters for free and supported panel boundary conditions up to a frequency of 600 Hz. Free boundary conditions were simulated by providing soft foam support under the four comers of the panel or by suspending the panel from elastic bands. Supported boundary conditions were obtained by clamping the panel between plastic tubing seated in grooves along the perimeter of a stiff and heavy frame. The frame was installed in the transmission loss window of the Structural Acoustic Loads and Transmission (SALT) facility at NASA Langley Research Center. The structural response of the curved panel due to point force excitation was predicted using MSC/NASTRAN and the radiated sound was computed with COMET/Acoustics. The predictions were compared with the results from experimental modal surveys and forced response tests on the fuselage panel. The finite element models were refined and updated to provide optimum comparison with the measured modal data. Excellent agreement was obtained between the numerical and experimental modal data for the free as well as for the supported boundary conditions. Frequency response functions (FRF) were computed relating the input force excitation at one panel location to the surface acceleration response at five panel locations. Frequency response functions were measured at the same locations on the test specimen and were compared with the calculated FRF values. Good agreement was obtained for the real and imaginary parts of the transfer functions when modal participation was allowed up to 3000 Hz. The validated finite element model was used to predict the surface velocities due to the point force excitation. Good agreement was obtained between the spatial characteristics of the predicted and measured surface velocities. The measured velocity data were input into the acoustic boundary element code to compute the sound radiated by the panel. The predicted sound pressure levels in the far-field of the panel agreed well with the sound pressure levels measured at the same location.

Grosveld, Ferdinand W.↗
DOE OSTI2021

Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics

Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.

Hu, Jing↗
NASA NTRS1989

Computing Viscoplastic Behavior Of A Material

Finite-element implementation developed for Robinson's unified model of viscoplasticity. Implemented via MARC general-purpose finite-element computer program by incorporating all of nonlinearity of material into initial load vector and treating vector as pseudo body force in finite-element equilibrium equations. In technique, each finite increment split into several equal subincrements, constitutive equations integrated by explicit Euler forward-difference method. HYPELA subroutine in MARC includes adaptive integration scheme selecting optimal sizes of subincrements. Spline function smooths discontinuous boundaries in Robinson's model to facilitate finite-element numerical calculations. Method used to calculate elastic, plastic, and creep deformations in solid materials at various temperatures; also used to study thermomechanical behavior of such things as hot pressure vessels and highly stressed components in engines.

Arya, V. K.↗
DOE OSTI2022

Local conformations and heterogeneities in structures and dynamics of isotactic polypropylene adsorbed onto carbon fiber

Carbon fiber (CF) reinforced polymers (CFRPs) have experienced widespread use in various industries. One of the most important parameters that controls the macroscopic property of CFRPs is the interface between a polymer matrix and CF. There is growing evidence to suggest the formation of a bound polymer layer (BPL), i.e., polymer chains that physically adsorb on a filler surface. However, this interface is always in contact with the thicker part of a polymer matrix, rendering its understanding a difficult task. Therein, we use CF-reinforced isotactic polypropylene (iPP) as a rational CFRP. To characterize the BPL on the CF surface, we extracted it from the CFRP using solvent-rinsing with p-xylene. The physical and thermal properties of the BPL were characterized by differential scanning calorimetry and thermogravimetric analysis, while its microscopic structures and dynamics were probed by small-angle neutron scattering and quasi-elastic neutron scattering (QENS) techniques. Subsequently, we employed atomistic molecular dynamics (MD) simulations to complement the QENS results above the bulk melting temperature and reveal details that were experimentally inaccessible. In this work, we observed that the degree of crystallinity of the BPL was quite lower than the bulk, while the melting temperature of the BPL remained the same as the bulks. Within the given length and time scales probed by QENS and MD, we also observed that most of the bound chains were mobile, with the formation of a high-density region (less than 1 nm in thickness) near the CF surface. The segmental dynamics of the bound chains probed by both QENS and MD were also much faster than those of the free chains, possibly due to the presence of a free surface region at the topmost surface of the BPL. Furthermore, the MD results demonstrated that the backbone chains and side groups lie nearly flat on the CF surface, which is the driving force for the flattening process of the iPP BPL to overcome the conformational entropy loss in the total free energy.

36 MATERIALS SCIENCE↗