Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic instability”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

Some conclusions regarding the aeroelastic stability of hingeless helicopter blades in hover and in forward flight

In this paper results and conclusions obtained from the study of the aeroelastic instability of hingeless helicopter blades are presented. First, the large amplitude coupled flap-lag equations of motion of a hingeless elastic helicopter blade are solved using an asymptotic expansion procedure in multiple time scales. Both hover and forward flight cases are considered. Stability boundaries and amplitudes of nonlinear response are obtained. From these, the importance of the nonlinear coupling and the effect of the periodic coefficients is determined. Next, using a system of linearized coupled flap-lag-pitch equations in hover, various divergence mechanisms for hingeless blades are shown. Finally, the flutter boundaries for coupled flap-lag-pitch are obtained. The effect of the torsional degree of freedom on the flap-lag type of instability is investigated. Similarly the effect of lag on the flap-pitch type of instability is considered. In addition, the effect of various blade parameters on the stability boundaries is shown.

Friedmann, P.↗
NASA NTRS1999

Rippling Instability of a Collapsing Bubble

The rippling instability of a liquid sheet was first observed by Debregeas, de Gennes, an Brochard-Wyart [Science 279, 1704 (1998)] on a hemispherical bubble resting on a free surface. Unlike a soap bubble, it collapses under its own weight while bursting, and folds into a wavy structure which breaks the original axisymmetry. In fact, this effect occurs for both purely elastic and purely viscous (liquid) sheets, and an analogy can be made between the two mechanisms. We present a theory for the onset of the instability in both cases, in which the growth of the corrugation out of an inextensible initial condition is governed by the competition between gravitational and bending (shearing) forces. The instability occurs for a range of densities, stiffnesses (viscosities), and sizes, a result which arises less from dynamics than from geometry, suggesting a wide validity. We further obtain a quantitative expression for the number of ripples. Finally, we present the results of experiments, which are consistent with our predictions.

daSilveira, Rava↗
DOE OSTI2021

R -mode Stability of GW190814's Secondary Component as a Supermassive and Superfast Pulsar

The nature of GW190814's secondary component m 2 of mass (2.50–2.67)M ⊙ in the mass gap between the currently known maximum mass of neutron stars and the minimum mass of black holes is currently under hot debate. Among the many possibilities proposed in the literature, the m 2 was suggested as a superfast pulsar while its r-mode stability against the run-away gravitational radiation through the Chandrasekhar-Friedman-Schutz mechanism is still unknown. Previously, Fortin et al. constructed a sample of 33 unified equations of state using the same nuclear interactions from the crust to the core consistently; from that sample we use those equations that fulfill all currently known astrophysical and nuclear physics constraints to compare the minimum frequency required for m 2 to rotationally sustain a mass greater than 2.50 M ⊙ with the critical frequency above which the r-mode instability occurs. Here, we use two extreme damping models assuming that the crust is either perfectly rigid or elastic. Using the stability of 19 observed low-mass X-ray binaries as an indication that the rigid crust damping of the r-mode dominates within the models studied, we find that m 2 is r-mode-stable while rotating with a frequency higher than 870.2 Hz (0.744 times its Kepler frequency of 1169.6 Hz) as long as its temperature is lower than about 3.9 × 10 7 K, further supporting the proposal that GW190814's secondary component is a supermassive and superfast pulsar.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

Influence of third-degree geometric nonlinearities on the vibration and stability of pretwisted, preconed, rotating blades

The governing coupled flapwise bending, edgewise bending, and torsional equations are derived including third-degree geometric nonlinear elastic terms by making use of the geometric nonlinear theory of elasticity in which the elongations and shears are negligible compared to unity. These equations are specialized for blades of doubly symmetric cross section with linear variation of pretwist over the blade length. The nonlinear steady state equations and the linearized perturbation equations are solved by using the Galerkin method, and by utilizing the nonrotating normal modes for the shape functions. Parametric results obtained for various cases of rotating blades from the present theoretical formulation are compared to those produced from the finite element code MSC/NASTRAN, and also to those produced from an in-house experimental test rig. It is shown that the spurious instabilities, observed for thin, rotating blades when second degree geometric nonlinearities are used, can be eliminated by including the third-degree elastic nonlinear terms. Furthermore, inclusion of third degree terms improves the correlation between the theory and experiment.

Subrahmanyam, K. B.↗
NASA NTRS1987

Influence of third-degree geometric nonlinearities on the vibration and stability of pretwisted, preconed, rotating blades

The governing coupled flapwise bending, edgewise bending, and torsional equations are derived including third-degree goemetric nonlinear elastic terms by making use of the geometric nonlinear theory of elasticity in which the elongations and shears are negligible compared to unity. These equations are specialized for blades of doubly symmetric cross section with linear variation of Pretwist over the blade length. The nonlinear steady state equations and the linearized perturbation equations are solved by using the Galerkin method, and by utilizing the nonrotating normal modes for the shape functions. Parametric results obtained for various cases of rotating blades from the present theoretical formulation are compared to those produced from the finite element code MSC/NASTRAN, and also to those produced from an in-house experimental test rig. It is shown that the spurious instabilities, observed for thin, rotating blades when second degree geometric nonlinearities are used, can be eliminated by including the third-degree elastic nonlinear terms. Furthermore, inclusion of third degree terms improves the correlation between the theory and experiment.

Subrahmanyam, K. B.↗
DOE OSTI2022

Design optimization for Richtmyer–Meshkov instability suppression at shock-compressed material interfaces

We report the Richtmyer–Meshkov instability (RMI) is a phenomenon that occurs at the interface of two substances of different densities due to an impulsive acceleration, such as a shock wave passing through this interface. Under these conditions, the instability can be seen as interface perturbations begin to grow into narrow jets or spikes of one substance that propagate into the other. In some cases, this interface may involve an elastic–plastic material, which can play a significant role in the development and behavior of the RMI. The ability to effectively control RMI jetting and spike growth is one major limiting factor in technological challenges, such as inertial confinement fusion, that involve using high-pressure shock waves to implode a fuel target. The propagation of RMI growth can lead to increased asymmetry in this implosion process and significantly reduce the obtained energy yield. We use hydrodynamics simulations of impactor shock-compression experiments and methods based in design optimization to suppress RMI spike growth by altering the geometry and other properties of a shock-compressed elastic–plastic material target that shares an interface with atmospheric air. These hydrodynamics simulations use an arbitrary Lagrangian–Eulerian method with a high-order finite element approach. Our results demonstrate that RMI suppression can be achieved by intentionally creating a separate upstream interface instability to counteract the growth of long narrow RMI spikes at an interface with initial perturbations.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Numerical simulation of the non-Newtonian mixing layer

This work is a continuing effort to advance our understanding of the effects of polymer additives on the structures of the mixing layer. In anticipation of full nonlinear simulations of the non-Newtonian mixing layer, we examined in a first stage the linear stability of the non-Newtonian mixing layer. The results of this study show that, for a fluid described by the Oldroyd-B model, viscoelasticity reduces the instability of the inviscid mixing layer in a special limit where the ratio (We/Re) is of order 1 where We is the Weissenberg number, a measure of the elasticity of the flow, and Re is the Reynolds number. In the present study, we pursue this project with numerical simulations of the non-Newtonian mixing layer. Our primary objective is to determine the effects of viscoelasticity on the roll-up structure. We also examine the origin of the numerical instabilities usually encountered in the simulations of non-Newtonian fluids.

Azaiez, Jalel↗
NASA NTRS1999

Development of a Practical Methodology for Elastic-Plastic and Fully Plastic Fatigue Crack Growth

A practical engineering methodology has been developed to analyze and predict fatigue crack growth rates under elastic-plastic and fully plastic conditions. The methodology employs the closure-corrected effective range of the J-integral, delta J(sub eff) as the governing parameter. The methodology contains original and literature J and delta J solutions for specific geometries, along with general methods for estimating J for other geometries and other loading conditions, including combined mechanical loading and combined primary and secondary loading. The methodology also contains specific practical algorithms that translate a J solution into a prediction of fatigue crack growth rate or life, including methods for determining crack opening levels, crack instability conditions, and material properties. A critical core subset of the J solutions and the practical algorithms has been implemented into independent elastic-plastic NASGRO modules. All components of the entire methodology, including the NASGRO modules, have been verified through analysis and experiment, and limits of applicability have been identified.

McClung, R. C.↗
NASA NTRS1999

Development of a Practical Methodology for Elastic-Plastic and Fully Plastic Fatigue Crack Growth

A practical engineering methodology has been developed to analyze and predict fatigue crack growth rates under elastic-plastic and fully plastic conditions. The methodology employs the closure-corrected effective range of the J-integral, (Delta)J(sub eff), as the governing parameter. The methodology contains original and literature J and (Delta)J solutions for specific geometries, along with general methods for estimating J for other geometries and other loading conditions, including combined mechanical loading and combined primary and secondary loading. The methodology also contains specific practical algorithms that translate a J solution into a prediction of fatigue crack growth rate or life, including methods for determining crack opening levels, crack instability conditions, and material properties. A critical core subset of the J solutions and the practical algorithms has been implemented into independent elastic-plastic NASGRO modules. All components of the entire methodology, including the NASGRO modules, have been verified through analysis and experiment, and limits of applicability have been identified.

McClung, R. C.↗
NASA NTRS1976

Variational theorems for superimposed motions in elasticity, with application to beams

Variational theorems are presented for a theory of small motions superimposed on large static deformations and governing equations for prestressed beams on the basis of 3-D theory of elastodynamics. First, the principle of virtual work is modified through Friedrichs's transformation so as to describe the initial stress problem of elastodynamics. Next, the modified principle together with a chosen displacement field is used to derive a set of 1-D macroscopic governing equations of prestressed beams. The resulting equations describe all the types of superimposed motions in elastic beams, and they include all the effects of transverse shear and normal strains, and the rotatory inertia. The instability of the governing equations is discussed briefly.

Doekmeci, M. C.↗
NASA NTRS1978

The dynamics and control of large flexible space structures. Part B: Development of continuum model and computer simulation

The equations of motion of an arbitrary flexible body in orbit were derived. The model includes the effects of gravity with all its higher harmonics. As a specific example, the motion of a long, slender, uniform beam in circular orbit was modelled. The example considers both the inplane and three dimensional motion of the beam in orbit. In the case of planar motion with only flexible vibrations, the pitch motion is not influenced by the elastic motion of the beam. For large values of the square of the ratio of the structural modal frequency to the orbital angular rate the elastic motion was decoupled from the pitch motion. However, for small values of the ratio and small amplitude pitch motion, the elastic motion was governed by a Hill's 3 term equation. Numerical simulation of the equation indicates the possibilities of instability for very low values of the square of the ratio of the modal frequency to the orbit angular rate. Also numerical simulations of the first order nonlinear equations of motion for a long flexible beam in orbit were performed. The effect of varying the initial conditions and the number of modes was demonstrated.

Bainum, P. M.↗
NASA NTRS1978

Dynamics of a flexible body in orbit

The equations of motion of an arbitrary flexible body in orbit are derived. The model includes the effects of gravity with all its higher harmonics. As a specific example, the motion of a long, slender, uniform beam in circular orbit is modeled. The example considers the inplane motion of the beam in orbit. In the case of planar motion with only flexural vibrations, the pitch motion is not influenced by the elastic motion of the beam. For large values of the square of the ratio of the structural modal frequency to the orbital angular rate the elastic motion is decoupled from the pitch motion. However, for small values of this ratio and small amplitude pitch motion, the elastic motion is governed by a Hill's 3-term equation. Numerical simulation of this equation indicates the possibilities of instability for very low values of the square of the ratio of the modal frequency to the orbit angular rate.

Kumar, V. K.↗
NASA NTRS1981

Interactive aircraft flight control and aeroelastic stabilization

Several examples are presented in which flutter involving interaction between flight mechanics modes and elastic wind bending occurs for a forward swept wing flight vehicle. These results show the basic mechanism by which the instability occurs and form the basis for attempts to actively control such a vehicle.

Weisshaar, T. A.↗
NASA NTRS1992

Parametric studies for tiltrotor aeroelastic stability in high-speed flight

The influence of several system design parameters on tiltrotor aeroelastic stability is examined for the high-speed (axial) flight mode. Coupling of the rotor flapping modes with the wing elastic modes produces a whirl motion, typical of tiltrotors, that can become unstable at high speeds. The sensitivity of this instability with respect to rotor frequencies, wing stiffness, forward wing sweep, and rotor thrust level is examined. Some important new trends are identified regarding the role of blade lag dynamics and forward wing sweep in tiltrotor aeroelastic stability. The blade lag frequency may be tuned to improve tiltrotor stability, and forward wing sweep is destabilizing because of changes in rotor force components associated with the sweep.

Nixon, Mark W.↗
DOE OSTI2023

Order-disorder charge density wave instability in the kagome metal (Cs,Rb)V 3 Sb 5

The origin of the charge density wave phases in the kagome metal compound AV 3 Sb 5 is still under great scrutiny. Here, we combine diffuse and inelastic x-ray scattering to identify a 3-dimensional precursor of the charge order at the L point that condenses into a CDW through a first order phase transition. The quasi-elastic critical scattering indicates that the dominant contribution to the diffuse precursor is the elastic central peak without phonon softening. However, the inelastic spectra show a small broadening of the Einstein-type phonon mode on approaching T CDW . Our results point to the situation where the Fermi surface instability at the L point is of order-disorder type with critical growth of quasi-static domains. The experimental data indicate that the CDW consists on an alternating Star of David and trihexagonal distortions and its dynamics goes beyond the classical weak-coupling scenario and is discussed within strong-electron phonon coupling and non-adiabatic models.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Machine Learning Predicts the Timing and Shear Stress Evolution of Lab Earthquakes Using Active Seismic Monitoring of Fault Zone Processes

Abstract Machine learning (ML) techniques have become increasingly important in seismology and earthquake science. Lab‐based studies have used acoustic emission data to predict time‐to‐failure and stress state, and in a few cases, the same approach has been used for field data. However, the underlying physical mechanisms that allow lab earthquake prediction and seismic forecasting remain poorly resolved. Here, we address this knowledge gap by coupling active‐source seismic data, which probe asperity‐scale processes, with ML methods. We show that elastic waves passing through the lab fault zone contain information that can predict the full spectrum of labquakes from slow slip instabilities to highly aperiodic events. The ML methods utilize systematic changes in P‐wave amplitude and velocity to accurately predict the timing and shear stress during labquakes. The ML predictions improve in accuracy closer to fault failure, demonstrating that the predictive power of the ultrasonic signals improves as the fault approaches failure. Our results demonstrate that the relationship between the ultrasonic parameters and fault slip rate, and in turn, the systematically evolving real area of contact and asperity stiffness allow the gradient boosting algorithm to “learn” about the state of the fault and its proximity to failure. Broadly, our results demonstrate the utility of physics‐informed ML in forecasting the imminence of fault slip at the laboratory scale, which may have important implications for earthquake mechanics in nature.

58 GEOSCIENCES↗
DOE OSTI2020

A fluid-structure interaction solver for investigating torsional galloping in solar-tracking photovoltaic panel arrays

Solar-tracking photovoltaic arrays are susceptible to aeroelastic fluttering during high-wind events. This dynamic fluttering behavior can grow in amplitude until the panels enter an unstable mode known as torsional galloping which can lead to panel failure or total array destruction. To better understand the physics of the torsional galloping phenomenon and to inform the discussion around panel design and recommended panel stow positions during high wind events, a fluid-structure interaction solver composed of a simulated atmospheric boundary layer with simplified panel structural responses was designed. The simulation choices and features of this solver were informed by the geometry and physical properties of an experimental panel array known to exhibit torsional galloping behavior during hind-wind events. Furthermore, these simulations revealed that the torsional galloping instability is driven by a combination of cyclic vortex shedding from the sun-facing side of the panel and the elastic properties of the torque tube linking the panel assemblies. Testing different stow angles across a range of wind speeds indicates that panels are generally more stable when stowed at negative angles where the leading edge is closer to the ground, hypothesized to be due to ground-blocking effects. These results are supplemented by a discussion of stability trends noted during testing and possible implications when considering multi-row array interactions.

14 SOLAR ENERGY↗
DOE OSTI2024

Edgewise Structural Damping of a 2.8-MW Land-Based Wind Turbine Rotor Blade

Modern wind turbines push the predictive capabilities of state-of-the-art aero-servo-elastic tools. The existing limitations hide across the numerical tool chain and can result in serious issues, such as missing aeroelastic instabilities during the design phase. Structural damping is an input that is usually hard to estimate, but also has a major impact on the turbine behavior. In this paper, we discuss an experiment that aims to accurately quantify the structural damping characterizing the edgewise modes of modern wind turbine blades. The experiment is carried out on a 2.8-MW land-based wind turbine and features a fast yaw actuation that induces an edgewise motion on one of the three blades. The Covariant-subspace system identification (Cov-SSI) method is then used to post-process the blade root moment to estimate the short-term edgewise structural damping. Despite limitations of the Cov-SSI method, which consistently under-predicts the absolute values of damping, we observe that structural damping decreases across the first three blade edgewise modes, which is different from the stiffness-proportional damping model that assumes that structural damping increases with the modes. This paper argues that a stiffness-proportional damping model, which is implemented in most aeroelastic tools, is therefore not conservative and might hide aeroelastic instabilities that can instead appear in the field.

17 WIND ENERGY↗