Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear stability analysis”

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 685 records · Page 38

Stability of flow over axisymmetric bodies with porous suction strips

Linear triple deck, closed form solutions for mean-flow quantities are developed for axisymmetric incompressible flow past a body with porous strips. The solutions account for upstream influence and are linear superpositions of the flow past the body without suction plus the perturbations due to the suction strips. Flow past the suctionless body is calculated using the Transition Analysis Program System, and a simple linear optimization scheme to determine number, spacing, and mass flow rate through the strips on an axisymmetric body is developed using the linear, triple-deck, closed-form solutions. The theory is demonstrated by predicting optimal strip distributions, and the effect of various adverse pressure-gradient situations on stability is studied.

Nayfeh, A. H.↗

Numerical analysis of stiffened shells of revolution. Volume 1: Theory manual for STARS-2S, 2B, 2V digital computer programs

The theoretical analysis background for the STARS-2 (shell theory automated for rotational structures) program is presented. The theory involved in the axisymmetric nonlinear and unsymmetric linear static analyses, and the stability and vibrations (including critical rotation speed) analyses involving axisymmetric prestress are discussed. The theory for nonlinear static, stability, and vibrations analyses, involving shells with unsymmetric loadings are included.

Svalbonas, V.↗

Stability characteristics of elastic airplane

System of computer programs uses linear theories to evaluate static and dynamic stability, trim state, inertial, and aerodynamic loading, and elastic deformations of aircraft configurations at subsonic and supersonic speeds. Primary emphasis is on analysis of stability and control characteristics of flexible aircraft. System also solves simple rigid aerodynamic problems.

Erickson, L. L.↗

An analysis of a nonlinear instability in the implementation of a VTOL control system

The contributions to nonlinear behavior and unstable response of the model following yaw control system of a VTOL aircraft during hover were determined. The system was designed as a state rate feedback implicit model follower that provided yaw rate command/heading hold capability and used combined full authority parallel and limited authority series servo actuators to generate an input to the yaw reaction control system of the aircraft. Both linear and nonlinear system models, as well as describing function linearization techniques were used to determine the influence on the control system instability of input magnitude and bandwidth, series servo authority, and system bandwidth. Results of the analysis describe stability boundaries as a function of these system design characteristics.

Weber, J. M.↗

Simultaneous stabilization and simultaneous pole placement by nonswitching dynamic compensation

The 'simultaneous stabilization problem' is defined and theorems are proposed for its solution. The problem consists in answering the question: given an r-tuple G sub 1(s), G sub r(s) of p x m proper transfer functions, does there exist a compensator K(s) such that the closed loop systems G sub 1(s) (I+K(s)G sub 1(s)) (-1), G sub r(s) (I+K(s) G sub r(s)) (-1) are (internally) stable. This question arises in reliability theory, where G sub 2(s), G sub r(s) represents a plant G sub 1(s) operating in various modes of failure and K(s) is a nonswitching stabilizing compensator. It is important in the stability analysis and design of a plant which can be switched into various operating modes. The simultaneous stabilization problem can also apply to the stabilization of a nonlinear system which is linearized at several equilibria. Conditions are defined for pole placement and the generalized Sylvestor matrix is discussed.

Ghosh, B. K.↗

Simultaneous stabilization and simultaneous pole-placement by nonswitching dynamic compensation

The 'simultaneous stabilization problem' is defined and theorems are proposed for its solution. The problem consists in answering the question: given an r-tuple G sub 1(s), G sub r(s) of p x m proper transfer functions, does there exist a compensator K(s) such that the closed loop systems G sub 1(s) (I+K(s)G sub 1(s)) (-1), G sub r(s) (I+K(s) G sub r(s)) (-1) are (internally) stable. This question arises in reliability theory, where G sub 2(s), G sub r(s) represents a plant G sub 1(s) operating in various modes of failure and K(s) is a nonswitching stabilizing compensator. It is important in the stability analysis and design of a plant which can be switched into various operating modes. The simultaneous stabilization problem can also apply to the stabilization of a nonlinear system which is linearized at several equilibria. Conditions are defined for pole placement and the generalized Sylvestor matrix is discussed. Previously announced in STAR as N82-31031

Ghosh, B. K.↗

A Measure-Theoretic Proof of the Markov Property for Hybrid Systems with Markovian Inputs

The behavior of a general hybrid system in discrete time can be represented by a non-linear difference equation x(k+1) = Fk(x(k), theta(k)), where theta(k) is assumed to be a finite state Markov chain. An important step in the stability analysis of these systems is to establish the Markov property of (x(k), theta(k)). There are, however, no complete proofs of this property which are simple to understand. This paper aims to correct this problem by presenting a complete and explicit proof, which uses only basic measure-theoretical concepts.

Tejada, Arturo↗

Modeling of the ECCD injection effect on the Heliotron J and LHD plasma stability

The aim of the study is to analyze the stability of the energetic particle modes (EPM) and Alfven Eigenmodes (AE) in Helitron J and LHD plasma if the electron cyclotron current drive (ECCD) is applied. Additionally, the analysis is performed using the code FAR3d that solves the reduced MHD equations describing the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moments for the energetic particle (EP) species, including the effect of the acoustic modes. The Landau damping and resonant destabilization effects are added via the closure relation. The simulation results show that the n = 1 EPM and n = 2 global AE (GAE) in Heliotron J plasma can be stabilized if the magnetic shear is enhanced at the plasma periphery by an increase (co-ECCD injection) or decrease (ctr-ECCD injection) of the rotational transform at the magnetic axis ($\rlap{-} \iota_{0}$). In the ctr-ECCD simulations, the EPM/AE growth rate decreases only below a given $\rlap{-} \iota_{0}$, similar to the ECCD intensity threshold observed in the experiments. In addition, ctr-ECCD simulations show an enhancement of the continuum damping. The simulations of the LHD discharges with ctr-ECCD injection indicate the stabilization of the n = 1 EPM, n = 2 toroidal AE (TAE) and n = 3 TAE, caused by an enhancement of the continuum damping in the inner plasma leading to a higher EP β threshold with respect to the co- and no-ECCD simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theoretical and experimental research in aeroelastic stability of an advanced bearingless rotor for future helicopters

The aeroelastic stability of a shaft-fixed bearingless rotor is analyzed in wind-tunnel tests for a wide range of operating conditions in order to determine whether such a system could be made aeroelastically stable without incorporating auxiliary dampers. The model rotor and blade properties are determined and used as an input to a bearingless-rotor analysis. Theoretical predictions are compared with experimental results in hover and forward flights. The analysis predicts the lag mode damping satisfactorily for collective pitch between 5 deg and 10 deg; however, the quasi-steady linear aerodynamic modeling overpredicts the damping values for higher collective pitch settings. It is noted that soft blade pitch links improve aeroelastic stability in hover and at low advance ratio.

Wang, James M.↗

Radiation effects on stress evolution and dimensional stability of large fusion energy structures

Here we assess the effects of neutron irradiation on the deformation and stress evolution of large-scale fusion energy structures. This is accomplished through non-linear finite element structural analysis of the coupled thermal and mechanical fields at the Beginning-Of-Life (BOL), at 45 dpa, and at 90 dpa. Radiation effects include volumetric swelling and the influence of radiation on the mechanical properties. The system studied here is a large section of a full inboard module of an integrated structure comprising the First Wall and Blanket (FW/B) of a Dual Cooled Lithium-Lead (DCLL) energy conversion unit in the Fusion Nuclear Science Facility (FNSF). The structural material is the ferritic/martensitic steel F82H. We analyze several radiation effects phenomena that can lead to significant impact on the mechanical design and lifetime of the structure. These include volumetric swelling effects on the assembly and disassembly of modules and the impact of spatial gradients in neutron damage and helium generation on deformation, stress evolution, and plastic strain accumulation. We show that the stress state is strongly influenced by the accumulation of swelling strains, radiation hardening and softening, and by spatial gradients in displacement damage and helium gas generation. Several key regions in the FW/B structure are identified where severe plastic strains accumulate and may be potential sites for failure, especially during cyclic reactor operations. We finally show that minimum gaps between Inboard (IB) blanket sectors must be present during assembly so as to avoid excessive stress and strain in FW/B modules critical regions and raising chances of failure. A range of structural displacements are shown, from as low as 4 mm at BOL for free side walls to as high as 46 mm at 90 dpa for fully constrained side walls.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of Curved-Plate Elements for the Exact Buckling Analysis of Composite Plate Assemblies Including Transverse-Shear Effects

The analytical formulation of curved-plate non-linear equilibrium equations including transverse-shear-deformation effects is presented. The formulation uses the principle of virtual work. A unified set of non-linear strains that contains terms from both physical and tensorial strain measures is used. Linearized, perturbed equilibrium equations (stability equations) that describe the response of the plate just after buckling occurs are then derived after the application of several simplifying assumptions. These equations are then modified to allow the reference surface of the plate to be located at a distance z(sub c) from the centroidal surface. The implementation of the new theory into the VICONOPT exact buckling and vibration analysis and optimum design computer program is described as well. The terms of the plate stiffness matrix using both Classical Plate Theory (CPT) and first-order Shear-Deformation Plate Theory (SDPT) are presented. The necessary steps to include the effects of in-plane transverse and in-plane shear loads in the in-plane stability equations are also outlined. Numerical results are presented using the newly implemented capability. Comparisons of results for several example problems with different loading states are made. Comparisons of analyses using both physical and tensorial strain measures as well as CPT and SDPF are also made. Results comparing the computational effort required by the new analysis to that of the analysis currently in the VICONOPT program are presented. The effects of including terms related to in-plane transverse and in-plane shear loadings in the in-plane stability equations are also examined. Finally, results of a design-optimization study of two different cylindrical shells subject to uniform axial compression are presented.

McGowan, David Michael↗

A unified perspective on robot control - The energy Lyapunov function approach

A unified framework for the stability analysis of robot tracking control is presented. By using an energy-motivated Lyapunov function candidate, the closed-loop stability is shown for a large family of control laws sharing a common structure of proportional and derivative feedback and a model-based feedforward. The feedforward can be zero, partial or complete linearized dynamics, partial or complete nonlinear dynamics, or linearized or nonlinear dynamics with parameter adaptation. As result, the dichotomous approaches to the robot control problem based on the open-loop linearization and nonlinear Lyapunov analysis are both included in this treatment. Furthermore, quantitative estimates of the trade-offs between different schemes in terms of the tracking performance, steady state error, domain of convergence, realtime computation load and required a prior model information are derived.

Wen, John T.↗

Three-Dimensional Nonequilibrium Boundary Layer Stability & Transition Analysis For Martian Entry

In this work, we apply the boundary layer stability code LASTRAC, modified to include chemical and thermal nonequilibrium capabilities, to a simulated test case based on literature, with three-dimensional effects. LASTRAC provides Parabolized Stability Equations (PSE)as well as Linear Stability Theory (LST) to predict the stability of a boundary layer and transition with semiempirical eN methods. Tests of the Mars Science Laboratory (MSL)indicated a range of laminar to transitional flow conditions. In order to reduce the weight of Thermal Protection Systems (TPS), reliable prediction of transition and an understanding of the environmental factors that contribute to transition are necessary. This work is a step towards addressing this need, including the effects of thermochemical nonequilibrium but neglecting radiation, ablation, and ionization effects. Our results indicate that an oblique first mode produces maximum N-Factors on the symmetry plane, however with N-Factors too small to explain transition. Although we are able to show nonequilibrium three-dimensional PSE results inCO2, further study with improved mean flow solutions, accounting for some of the uncertainties inherent in this analysis, or evaluation of alternative transition mechanisms is needed to explain transition in the wind tunnel cases examined here.

H L Kline↗

Stochastic Control Synthesis of Systems with Structured Uncertainty

This paper presents a study on the design of robust controllers by using random variables to model structured uncertainty for both SISO and MIMO feedback systems. Once the parameter uncertainty is prescribed with probability density functions, its effects are propagated through the analysis leading to stochastic metrics for the system's output. Control designs that aim for satisfactory performances while guaranteeing robust closed loop stability are attained by solving constrained non-linear optimization problems in the frequency domain. This approach permits not only to quantify the probability of having unstable and unfavorable responses for a particular control design but also to search for controls while favoring the values of the parameters with higher chance of occurrence. In this manner, robust optimality is achieved while the characteristic conservatism of conventional robust control methods is eliminated. Examples that admit closed form expressions for the probabilistic metrics of the output are used to elucidate the nature of the problem at hand and validate the proposed formulations.

Padula, Sharon L.↗

On the Markov Property for Nonlinear Discrete-Time Systems with Markovian Inputs

The behavior of a general hybrid system in discrete-time can be represented by a non-linear difference equation x(k + 1) = Fk(x(k), theta(k)), where theta(k) is assumed to be a finite-state Markov chain. An important step in the stability analysis of these systems is to establish the Markov property of (x(k), theta(k)). There are, however, no complete proofs of this property which are simple to understand. This paper aims to correct this problem by presenting a complete and explicit proof, which uses only fundamental measure-theoretical concepts.

Tejada, Arturo↗

Wave Driven Non-Linear Flow Oscillator for the 22-Year Solar Cycle

We propose that waves generate an oscillation in the Sun to account for the 22-year magnetic cycle. The mechanism we envision is analogous to that driving the Quasi Biennial Oscillation (QBO) observed in the terrestrial atmosphere, which is well understood in principal. Planetary waves and gravity waves deposit momentum in the background atmosphere and accelerate the flow under viscous dissipation. Analysis shows that such a momentum source represents a non-linearity of third or generally odd order, which generates also the fundamental frequency/period so that an oscillation is maintained without external time dependent forcing. For the Sun, we propose that the wave driven oscillation would occur just below the convection region, where the buoyancy frequency or convective stability becomes small to favor wave breaking and wave mean flow interaction. Using scale analysis to extrapolate from terrestrial to solar conditions, we present results from a simplified analytical model, applied to the equator, that incorporates Hines'Doppler Spread Parameterization for gravity waves (GW). Based on a parametric study, we conclude: (1) Depending on the adopted horizontal wavelengths of GW's, wave amplitudes < 10 m/s can be made to produce oscillating zonal winds of about 25 m/s that should be large enough to generate a corresponding oscillation in the main poloidal magnetic field; (2) The oscillation period can be made to be 22 years provided the buoyancy frequency (stability) is sufficiently small, which would place the oscillating wind field near the base of the convection region; (3) In this region, the turbulence associated with wave processes would be enhanced by low stability, and this also helps to produce the desired oscillation period and generate the dynamo currents that would produce the reversing magnetic field. We suggest that the above mechanism may also drive other long-period metronomes in planetary and stellar interiors.

Mayr, H. G.↗

Third-order resonance effects and the nonlinear stability of drop oscillations

The three-dimensional nonlinear oscillations of an isolated, inviscid drop with surface tension are studied by a multiple timescale analysis and pre-averaging applied to the variational principle for the appropriate Lagrangian. Amplitude equations are derived which describe the generic cubic resonance caused by the spatial degeneracy of the eigenfrequencies of the linear normal modes. This resonant coupling leads to the instability of the finite amplitude axisymmetric oscillations to small nonaxisymmetric perturbations, as is demonstrated here for the three- and four-lobed normal modes. Solutions to the interaction equations that describe finite amplitude, nonaxisymmetric traveling-wave solutions are also obtained and their stability is investigated. A nongeneric cubic resonance between the two-lobed and four-lobed oscillatory modes leads to quasi-periodic motions.

Natarajan, Ramesh↗

Analysis of limit cycles in control systems for joint dominated structures

An approach to the modeling of limit cycles due to joint nonlinearities in large joint dominated space structures is presented which makes it possible to predict limit cycles and determine their stability. An actively controlled truss structure with nonlinear joints is modeled as a linear system with nonlinear feedback by separating the joint load-displacement characteristics into a linear part, which prevails at large displacements, and a nonlinear part. By replacing the joints by their linear parts, it is possible to perform a standard model decomposition which yields a reduced order linear model. Linear control laws can be easily included into the linear part of the system; nonlinear control laws can be implemented but they must be fed back to the linear model. The modeling approach described here allows straightforward limit cycle analysis.

Mercadal, M.↗