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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

A variational approach to dynamics of flexible multibody systems

This paper presents a variational formulation of constrained dynamics of flexible multibody systems, using a vector-variational calculus approach. Body reference frames are used to define global position and orientation of individual bodies in the system, located and oriented by position of its origin and Euler parameters, respectively. Small strain linear elastic deformation of individual components, relative to their body references frames, is defined by linear combinations of deformation modes that are induced by constraint reaction forces and normal modes of vibration. A library of kinematic couplings between flexible and/or rigid bodies is defined and analyzed. Variational equations of motion for multibody systems are obtained and reduced to mixed differential-algebraic equations of motion. A space structure that must deform during deployment is analyzed, to illustrate use of the methods developed.

Wu, Shih-Chin↗
NASA NTRS1994

Late Cenozoic flexural deformation of the middle U.S. Atlantic passive margin

Despite the century-long recognition of regional epeirogeny along the middle Atlantic passive margin, relatively few studies have focused on understanding postrift uplift mechanisms. Here, we demonstrate that epeirogenic uplift of the central Appalachian Piedmont and subsidence of the Salisbury Embayment represent first-order, flexural isostatic processes driven by continental denudation and offshore deposition. Our results show that regional epeirogenic processes, present on all Atlantic-type passive margins, are best resolved by specific stratigraphic and geomorphic relationships, rather than topography. A simple one-dimensional geodynamic model, constrained by well-dated Baltimore Canyon trough, Coastal Plain, and lower Susquehanna River (piedmont) stratigraphy, simulates flexural deforamtion of the U.S. Atlantic margin. The model represents the passive margin lithosphree as a uniformly thick elastic plate, without horizontal compressive stresses, that deforms flexurally under the stress of strike-averaged, vertically applied line loads. Model results illustrate a complex interaction among margin stratigraphy and geomorphology, the isostatic repsonse to denudational and depositional processes, and the modulating influence of exogenic forces such as eustasy. The current elevation, with respect to modern sea level, of fluvial terraces and correlateive Coastal Plain deposits or unconformities is successfully predicted through the synthesis of paleotopography, eustatic change, and margin flexure. Results suggest that the middle U.S. Atlantic margin landward of East Coast Magnetic Anomaly is underlain by lithoshpere with an average elastic thickness of 40 km (flexural rigidity, D = 4 X 10(exp 23) N m), the margin experience an average, long-term denudation rate of approximately 10m/m.y., and the Piedmont has been flexurally upwaped between 35 and 130 meters in the last 15 m.y. Long term isostatic continental uplift resulting rom denudation and basin subsidence resulting rom sediment loading are accomodated primately by a convex-up flexural hinge, physiographically represented by the Fall Zone. Our results elucidate an inherent danger in using topography alone to constrain late-stage passive margin deformation mechanisms. Only through careful synthesis of field stratigraphic and geomorphic elements such as fluvial terraces, Coastal Plain deposits, and offshore stratigraphy can age control be extended from the offshore depositional setting to the erosionally dominated continent. This sudy demonstrates that despite a relatively subdued topography, the middle U.S. Atlantic margin experiences progressive flexural isostatic deformation similar to that proposed for high-relief margins characterized by great escarpments. Thus margin topographic diversity remains a function of other factors, such as lithospheric composition and/or structure, supracrustal stratigraphy and structure, degree of drainage integration, drainage divide migration and climate.

Pazzaglia, Frank J.↗
DOE OSTI2024

From anti-Arrhenius to Arrhenius behavior in a dislocation-obstacle bypass: Atomistic simulations and theoretical investigation

Dislocations are the primary carriers of plasticity in metallic materials. Understanding the basic mechanisms for dislocation movement is paramount to predicting the material mechanical response. Relying on atomistic simulations, we observe a transition from non-Arrhenius to Arrhenius behavior in the rate of an edge dislocation overcoming the long-range elastic interaction with a prismatic loop in tungsten. Close to the critical resolved shear stress, the process shows a non-Arrhenius behavior at low temperatures. However, as the temperature increases, the activation entropy starts to dominate, leading to a traditional Arrhenius-like behavior. Here, we have computed the activation entropy analytically along the minimum energy path following Schoeck’s method, which captures the cross-over between anti-Arrhenius and Arrhenius domains. Also, the Projected Average Force Integrator (PAFI), another simulation method to compute free energies along an initial transition path, exhibits considerable concurrence with Schoeck’s formalism. We conclude that entropic effects need to be considered to understand processes involving dislocations bypassing elastic barriers close to the critical resolved shear stress. More work needs to be performed to fully understand the discrepancies between Schoeck’s and PAFI compared to molecular dynamics.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Effective Elastic Modulus as a Function of Angular Leaf Span for Curved Leaves of Pyrolytic Boron Nitride

A theoretical equation was derived to predict the spring constant (load/deflection) for a simply supported cylindrical section with a line force applied at the center. Curved leaves of PBN were mechanically deformed and the force versus deflection data was recorded and compared to the derived theoretical equation to yield an effective modulus for each leaf. The effective modulus was found to vary from the pure shear modulus for a flat plate to a mixed mode for a half cylinder as a function of the sine of one half the angular leaf span. The spring constants of individual PBN leaves were usually predicted to within 30%.

Kaforey, M. L.↗
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.↗
NASA NTRS1975

Experimental transient and permanent deformation studies of steel-sphere-impacted or impulsively-loaded aluminum beams with clamped ends

The sheet explosive loading technique (SELT) was employed to obtain elastic-plastic, large-deflection transient and/or permanent strain data on simple well-defined structural specimens and materials: initially-flat 6061-T651 aluminum beams with both ends ideally clamped via integral construction. The SELT loading technique was chosen since it is both convenient and provides forcing function information of small uncertainty. These data will be useful for evaluating pertinent structural response prediction methods. A second objective was to obtain high-quality transient-strain data for a well-defined structural/material model subjected to impact by a rigid body of known mass, impact velocity, and geometry; large-deflection, elastic-plastic transient response conditions are of primary interest. The beam with both ends clamped and a steel sphere as the impacting body were chosen. The steel sphere was launched vertically by explosive propulsion to achieve various desired impact velocities. The sphere/beam impact tests resulted in producing a wide range of structural responses and permanent deformations, including rupture of the beam from excessive structural response in two cases. The transient and permanent strain data as well as the permanent deflection data obtained are of high quality and should be useful for checking and evaluating methods for predicting the responses of simple 2-d structures to fragment (sphere) impact. Transient strain data very close to the point of impact were not obtained over as long a time as desirable because the gage(s) in that region became detached during the transient response.

Witmer, E. A.↗
NASA NTRS1989

The effect of an unsteady three-dimensional wake on elastic blade-flapping eigenvalues in hover

This paper describes the formulation of a finite-state inflow model based on an acceleration potential and a helical wake geometry. The states of the model are coefficients of an inflow expansion in terms of a Fourier series (azimuthally) and of special polynomials (radially). The integrals over the wake are done in closed-form to obtain a set of ordinary differential equations for the inflow coefficients. The forcing functions for these equations are generalized forces which are integrals of the blade loading exactly as in structural dynamics. This model implicitly includes (for the hover case) Prandtl-Goldstein tip losses, dynamic inflow, and Theodorsen/Loewy lift deficiency. Thus, it is a fully three-dimensional unsteady wake model. Here, this model is coupled with elastic-blade equations in hover and eigenvalues are found. The results show that the three-dimensional wake has a large effect on the flap damping of all modes.

Peters, David A.↗
DOE OSTI2020

Analytical integration of the tractions induced by non-singular dislocations on an arbitrary shaped triangular quadratic element

Here, an analytical model is proposed to evaluate the nodal force induced by a segment of dislocation upon an arbitrary shaped triangular element. This calculation is required in hybrid methods that associate dislocation dynamics to boundary or finite element to solve simultaneously the evolution of large ensembles of dislocations with complex boundary conditions. Nodal forces are defined as the triple integration of the unbalanced traction field induced by a straight dislocation upon the surface of the element. Following our previous approach (Queyreau et al 2014 Modelling Simul. Mater. Sci. Eng. 22 035004) on a simpler geometry and in the case of linear isotropic elasticity, triple integrals are solved by sequences of integration by parts that exhibit recurrence relations. The traction field is defined and finite everywhere even at the core of dislocations, thanks to the use of the non-singular stress expression formulated by Cai et al (2006 J. Mech. Phys. Solids 54 561–587). The nodal force expressions can be used when considering both a single convolution or double convolution of the Green's function with the core distribution. A solution is also proposed for the case of a semi-infinite segment through the study of the asymptotic behavior of the analytical expressions. The proposed approach is exact and very computationally efficient. The choice of an arbitrary shaped triangular element and quadratic shape functions allow the consideration of complex geometries and comply with automatic meshing procedures. These analytical expressions could also be employed to estimate dislocation interactions with interfaces.

42 ENGINEERING↗
DOE OSTI2021

Deep learning of accurate force field of ferroelectric HfO 2

The discovery of ferroelectricity in HfO 2 -based thin films opens up new opportunities for using this silicon-compatible ferroelectric to realize low-power logic circuits and high-density nonvolatile memories. The functional performances of ferroelectrics are intimately related to their dynamic responses to external stimuli such as electric fields at finite temperatures. Molecular dynamics is an ideal technique for investigating dynamical processes on large length and time scales, though its applications to new materials are often hindered by the limited availability and accuracy of classical force fields. Here we present a deep neural network–based interatomic force field of HfO 2 learned from ab initio data using a concurrent learning procedure. The model potential is able to predict structural properties such as elastic constants, equation of states, phonon dispersion relationships, and phase transition barriers of various hafnia polymorphs with accuracy comparable with density functional theory calculations. Here, the validity of this model potential is further confirmed by the reproduction of experimental sequences of temperature-driven ferroelectric-paraelectric phase transitions of HfO 2 with isobaric-isothermal ensemble molecular dynamics simulations. We suggest a general approach to extend the model potential of HfO 2 to related material systems including dopants and defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Detailed study of quark-hadron duality in spin structure functions of the proton and neutron

Background: The response of hadrons, the bound states of the strong force (QCD), to external probes can be described in two different, complementary frameworks: as direct interactions with their fundamental constituents, quarks and gluons, or alternatively as elastic or inelastic coherent scattering that leaves the hadrons in their ground state or in one of their excited (resonance) states. The former picture emerges most clearly in hard processes with high momentum transfer, where the hadron response can be described by the perturbative expansion of QCD, while at lower energy and momentum transfers, the resonant excitations of the hadrons dominate the cross section. The overlap region between these two pictures, where both yield similar predictions, is referred to as quark-hadron duality and has been extensively studied in reactions involving unpolarized hadrons. Some limited information on this phenomenon also exists for polarized protons, deuterons, and 3 He nuclei, but not yet for neutrons. Purpose: In this paper, we present comprehensive and detailed results on the correspondence between the extrapolated deep inelastic structure function g 1 of both the proton and the neutron with the same quantity measured in the nucleon resonance region. Thanks to the fine binning and high precision of our data, and using a well-controlled perturbative QCD (pQCD) fit for the partonic prediction, we can make quantitative statements about the kinematic range of applicability of both local duality and global duality. Method: We use the most updated QCD global analysis results at high x from the Jefferson Lab Angular Momentum Collaboration to extrapolate the spin structure function g 1 into the nucleon resonance region and then integrate over various intervals in the scaling variable x. We compare the results with the large data set collected in the quark-hadron transition region by the CLAS Collaboration, including, for the first time, deconvoluted neutron data, integrated over the same intervals. We present this comparison as a function of the momentum transfer Q 2 . Results: We find that, depending on the integration interval and the minimum momentum transfer chosen, a clear transition to quark-hadron duality can be observed in both nucleon species. Furthermore, we show, for the first time, the approach to scaling behavior for g 1 measured in the resonance region at sufficiently high momentum transfer. Conclusions: Here, our results can be used to quantify the deviations from the applicability of pQCD for data taken at moderate energies and can help with extraction of quark distribution functions from such data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1972

Blade frequency program for nonuniform helicopter rotors, with automated frequency search

A computer program for determining the natural frequencies and normal modes of a lumped parameter model of a rotating, twisted beam, with nonuniform mass and elastic properties was developed. The program is used to solve the conditions existing in a helicopter rotor where the outboard end of the rotor has zero forces and moments. Three frequency search methods have been implemented. Including an automatic search technique, which allows the program to find up to the fifteen lowest natural frequencies without the necessity for input estimates of these frequencies.

Sadler, S. G.↗
NASA NTRS1971

Dielectric siphons.

The normally weak polarization force density, exerted on insulating dielectric liquids by a nonuniform electric field, is enhanced if high pressures are used. The nonuniform electric field acts as an elastic ?wall' to contain and guide the dielectric fluid. A general theory for these electrohydrodynamic (EHD) conduits has been developed. An illustrative example of the EHD conduits is the dielectric siphon consisting of two U-shaped electrodes held adjacent to each other by insulating nylon screws.

Jones, T. B.↗
NASA NTRS1973

Rotor dynamic considerations for large wind power generator systems

Successful large, reliable, low maintenance wind turbines must be designed with full consideration for minimizing dynamic response to aerodynamic, inertial, and gravitational forces. Much of existing helicopter rotor technology is applicable to this problem. Compared with helicopter rotors, large wind turbines are likely to be relatively less flexible with higher dimensionless natural frequencies. For very large wind turbines, low power output per unit weight and stresses due to gravitational forces are limiting factors. The need to reduce rotor complexity to a minimum favors the use of cantilevered (hingeless) rotor configurations where stresses are relieved by elastic deformations.

Ormiston, R. A.↗
NASA NTRS1989

Rift propagation

A model for rift propagation which treats the rift as a crack in an elastic plate which is filled from beneath by upwelling viscous asthenosphere as it lengthens and opens. Growth of the crack is driven by either remotely applied forces or the pressure of buoyant asthenosphere in the crack and is resisted by viscous stresses associated with filling the crack. The model predicts a time for a rift to form which depends primarily on the driving stress and asthenosphere viscosity. For a driving stress on the order of 10 MPa, as expected from the topography of rifted swells, the development of rifts over times of a few Myr requires an asthenosphere viscosity of 10 to the 16th Pa s (10 to the 17th poise). This viscosity, which is several orders of magnitude less than values determined by postglacial rebound and at least one order of magnitude less than that inferred for spreading center propagation, may reflect a high temperature or large amount of partial melting in the mantle beneath a rifted swell.

Parmentier, E. M.↗
NASA NTRS1920

Experience with Geared Propeller Drives for Aviation Engines

I. The development of the gear wheels: (a) bending stresses; (b) compressive stresses; (c) heating; (d) precision of manufacture. II. General arrangement of the gearing. III. Vibration in the shaft transmission. An overview is given of experience with geared propeller drives for aviation engines. The development of gear wheels is discussed with emphasis upon bending stresses, compressive stresses, heating, and precision in manufacturing. With respect to the general arrangement of gear drives for airplanes, some principal rules of mechanical engineering that apply with special force are noted. The primary vibrations in the shaft transmission are discussed. With respect to vibration, various methods for computing vibration frequency and the influence of elastic couplings are discussed.

Kutzbach, K↗
NASA NTRS1997

Force Model for Control of Tendon Driven Hands

Knowing the tendon forces generated for a given task such as grasping via a model, an artificial hand can be controlled. A two-dimensional force model for the index finger was developed. This system is assumed to be in static equilibrium, therefore, the equations of equilibrium were applied at each joint. Constraint equations describing the tendon branch connectivity were used. Gaussian elimination was used to solve for the unknowns of the Linear system. Results from initial work on estimating tendon forces in post-operative hands during active motion therapy were discussed. The results are important for understanding the effects of hand position on tendon tension, elastic effects on tendon tension, and overall functional anatomy of the hand.

Pena, Edward↗
DOE OSTI2020

Hydrokinetic Power Conversion Using Vortex-Induced Oscillation with Cubic Restoring Force

A cubic-spring restoring function with high-deformation stiffening is introduced to passively improve the harnessed marine hydrokinetic power by using flow-induced oscillations/vibrations (FIO/V) of a cylinder. In these FIO/V experiments, a smooth, rigid, single-cylinder on elastic end-supports is tested at Reynolds numbers ranging from 24,000 < Re < 120,000. The parameters of the tested current energy converter (CEC) are cubic stiffness and linear damping. Using the second generation of digital virtual spring-damping (Vck) controller developed by the Marine Renewable Energy Laboratory (MRELab), the cubic modeling of the oscillator stiffness is tested. Experimental results show the influence of the parameter variation on the amplitude, frequency, energy conversion, energy efficiency, and power of the converter. All experiments are conducted in the low turbulence-free surface water (LTFSW) channel of the MRELab of the University of Michigan. The main conclusions are: (1) The nonlinearity in the cubic oscillator is an effective way to extend the vortex-induced vibration (VIV) upper branch, which results in higher harnessing power and efficiency compared to the linear stiffness cylinder converter. (2) Compared to the linear converter, the overall power increase is substantial. The nonlinear power optimum, occurring at the end of the VIV upper branch, is 63% higher than its linear counterpart. (3) The cubic stiffness converter with low harnessing damping achieves consistently good performance in all the VIV regions because of the hardening restoring force, especially at higher flow velocity.

13 HYDRO ENERGY↗