Efficient and Robust Approaches to the Stability Analysis of Large Multibody Systems

Author(s):  
Olivier A. Bauchau ◽  
Jielong Wang

Linearized stability analysis methodologies that are applicable to large scale, multiphysics problems are presented in this paper. Two classes of closely related algorithms based on a partial Floquet and on an autoregressive approach, respectively, are presented in common framework that underlines their similarity and their relationship to other methods. The robustness of the proposed approach is improved by using optimized signals that are derived from the proper orthogonal modes of the system. Finally, a signal synthesis procedure based on the identified frequencies and damping rates is shown to be an important tool for assessing the accuracy of the identified parameters; furthermore, it provides a means of resolving the frequency indeterminacy associated with the eigenvalues of the transition matrix for periodic systems. The proposed approaches are computationally inexpensive and consist of purely post processing steps that can be used with any multiphysics computational tool or with experimental data. Unlike classical stability analysis methodologies, it does not require the linearization of the equations of motion of the system.

2011 ◽  
Vol 137 ◽  
pp. 72-76
Author(s):  
Wei Zhang ◽  
Xian Wen ◽  
Yan Qun Jiang

A proper orthogonal decomposition (POD) method is applied to study the global stability analysis for flow past a stationary circular cylinder. The flow database at Re=100 is obtained by CFD software, i.e. FLUENT, with which POD bases are constructed by a snapshot method. Based on the POD bases, a low-dimensional model is established for solving the two-dimensional incompressible NS equations. The stability of the flow solution is evaluated by a POD-Chiba method in the way of the eigensystem analysis for the velocity disturbance. The linear stability analysis shows that the first Hopf bifurcation takes place at Re=46.9, which is in good agreement with available results by other high-order accurate stability analysis methods. However, the calculated amount of POD is little, which shows the availability and advantage of the POD method.


Author(s):  
Kamal Kumar Basumatary ◽  
Gaurav Kumar ◽  
Karuna Kalita ◽  
Sashindra K Kakoty

Rotors supported on gas foil bearings have low damping characteristics, which limits its application. A possible solution could be an integration of a gas foil bearing with an electromagnetic actuator. This paper discusses the effect of electromagnetic actuators on the stability of a rotor supported on gas foil bearings. A coupled dynamic model combining the dynamics of gas foil bearing and electromagnetic actuator has been developed. The fluid film forces from the gas foil bearings and the electromagnetic forces from the electromagnetic actuators are integrated into the equations of motion of the rotor. The sub-synchronous vibration present in case of conventional gas foil bearings is reduced and the stability band of the rotor is increased due to the implementation of electromagnetic actuator.


Author(s):  
Ram Turaga

The influence of deterministic surface texture on the sub-synchronous whirl stability of a rigid rotor has been studied. Non-linear transient stability analysis has been performed to study the stability of a rigid rotor supported on two symmetric journal bearings with a rectangular dimple of large aspect ratio. The surface texture parameters considered are dimple depth to minimum film thickness ratio and the location of the dimple on the bearing surface. Journal bearings of different Length to diameter ratios have been studied. The governing Reynolds equation for finite journal bearings with incompressible fluid has been solved using the Finite Element Method under isothermal conditions. The trajectories of the journal center have been obtained by solving the equations of motion of the journal center by the fourth-order Runge-Kutta method. When the dimple is located in the raising part of the pressure curve the positive rectangular dimple is seen to decrease the stability whereas the negative rectangular dimple is seen to improve the stability of the rigid rotor.


Author(s):  
Zhanhong Wan ◽  
Saihua Huang ◽  
Zhilin Sun ◽  
Zhenjiang You

Purpose – The present work is devoted to the numerical study of the stability of shallow jet. The effects of important parameters on the stability behavior for large scale shallow jets are considered and investigated. Connections between the stability theory and observed features reported in the literature are emphasized. The paper aims to discuss these issues. Design/methodology/approach – A linear stability analysis of shallow jet incorporating the effects of bottom topography, bed friction and viscosity has been carried out by using the shallow water stability equation derived from the depth averaged shallow water equations in conjunction with both Chézy and Manning resistance formulae. Effects of the following main factors on the stability of shallow water jets are examined: Rossby number, bottom friction number, Reynolds number, topographic parameters, base velocity profile and resistance model. Special attention has been paid to the Coriolis effects on the jet stability by limiting the rotation number in the range of Ro∈[0, 1.0]. Findings – It is found that the Rossby number may either amplify or attenuate the growth of the flow instability depending on the values of the topographic parameters. There is a regime where the near cancellation of Coriolis effects due to other relevant parameters influences is responsible for enhancement of stability. The instability can be suppressed by the bottom friction when the bottom friction number is large enough. The amplification rate may become sensitive to the relatively small Reynolds number. The stability region using the Manning formula is larger than that using the Chézy formula. The combination of these effects may stabilize or destabilize the shallow jet flow. These results of the stability analysis are compared with those from the literature. Originality/value – Results of linear stability analysis on shallow jets along roughness bottom bed are presented. Different from the previous studies, this paper includes the effects of bottom topography, Rossby number, Reynolds number, resistance formula and bed friction. It is found that the influence of Reynolds number on the stability of the jet is notable for relative small value. Therefore, it is important to experimental investigators that the viscosity should be considered with comparison to the results from inviscid assumption. In contrast with the classical analysis, the use of multi-parameters of the base velocity and topographic profile gives an extension to the jet stability analysis. To characterize the large scale motion, besides the bottom friction as proposed in the related literature, the Reynolds number Re, Rossby number Ro, the topographic parameters and parameters controlling base velocity profile may also be important to the stability analysis of shallow jet flows.


2013 ◽  
Vol 392 ◽  
pp. 156-160
Author(s):  
Ju Seok Kang

Multibody dynamics analysis is advantageous in that it uses real dimensions and design parameters. In this study, the stability analysis of a railway vehicle based on multibody dynamics analysis is presented. The equations for the contact points and contact forces between the wheel and the rail are derived using a wheelset model. The dynamics equations of the wheelset are combined with the dynamics equations of the other parts of the railway vehicle, which are obtained by general multibody dynamics analysis. The equations of motion of the railway vehicle are linearized by using the perturbation method. The eigenvalues of these linear dynamics equations are calculated and the critical speed is found.


2005 ◽  
Vol 1 (1) ◽  
pp. 71-80 ◽  
Author(s):  
Olivier A. Bauchau ◽  
Jielong Wang

The linearized stability analysis of dynamical systems modeled using finite element-based multibody formulations is addressed in this paper. The use of classical methods for stability analysis of these systems, such as the characteristic exponent method or Floquet theory, results in computationally prohibitive costs. Since comprehensive multibody models are “virtual prototypes” of actual systems, the applicability to numerical models of the stability analysis tools that are used in experimental settings is investigated in this work. Various experimental tools for stability analysis are reviewed. It is proved that Prony’s method, generally regarded as a curve-fitting method, is equivalent, and sometimes identical, to Floquet theory and to the partial Floquet method. This observation gives Prony’s method a sound theoretical footing, and considerably improves the robustness of its predictions when applied to comprehensive models of complex multibody systems. Numerical and experimental applications are presented to demonstrate the efficiency of the proposed procedure.


1986 ◽  
Vol 39 (3) ◽  
pp. 331 ◽  
Author(s):  
B Eckhardt ◽  
JA Louw ◽  
W-H Steeb

We review two criteria which have been used to predict the onset of large scale stochasticity in Hamiltonian systems. We show that one of them, due to Toda and based on a local stability analysis of the equations of motion, is inconclusive. An approach based on the local Riemannian curvature K of trajectories correctly predicts chaos if K < 0 everywhere, but�no further conclusions can be drawn. New (counter-)examples are provided.


1977 ◽  
Vol 17 (01) ◽  
pp. 79-91 ◽  
Author(s):  
D.W. Peaceman

Abstract The usual linearized stability analysis of the finite-difference solution for two-phase flow in porous media is not delicate enough to distinguish porous media is not delicate enough to distinguish between the stability of equations using semi-implicit mobility and those using completely implicit mobility. A nonlinear stability analysis is developed and applied to finite-difference equations using an upstream mobility that is explicit, completely implicit, or semi-implicit. The nonlinear analysis yields a sufficient (though not necessary) condition for stability. The results for explicit and completely implicit mobilities agree with those obtained by the standard linearized analysis; in particular, use of completely implicit mobility particular, use of completely implicit mobility results in unconditional stability. For semi-implicit mobility, the analysis shows a mild restriction that generally will not be violated in practical reservoir simulations. Some numerical results that support the theoretical conclusions are presented. Introduction Early finite-difference, Multiphase reservoir simulators using explicit mobility were found to require exceedingly small time steps to solve certain types of problems, particularly coning and gas percolation. Both these problems are characterized percolation. Both these problems are characterized by regions of high flow velocity. Coats developed an ad hoc technique for dealing with gas percolation, but a more general and highly successful approach for dealing with high-velocity problems has been the use of implicit mobility. Blair and Weinaug developed a simulator using completely implicit mobility that greatly relaxed the time-step restriction. Their simulator involved iterative solution of nonlinear difference equations, which considerably increased the computational work per time step. Three more recent papers introduced the use of semi-implicit mobility, which proved to be greatly superior to the fully implicit method with respect to computational effort, ease of use, and maximum permissible time-step size. As a result, semi-implicit mobility has achieved wide use throughout the industry. However, this success has been pragmatic, with little or no theoretical work to justify its use. In this paper, we attempt to place the use of semi-implicit mobility on a sounder theoretical foundation by examining the stability of semi-implicit difference equations. The usual linearized stability analysis is not delicate enough to distinguish between the semi-implicit and completely implicit difference equation. A nonlinear stability analysis is developed that permits the detection of some differences between the stability of difference equations using implicit mobility and those using semi-implicit mobility. DIFFERENTIAL EQUATIONS The ideas to be developed may be adequately presented using the following simplified system: presented using the following simplified system: horizontal, one-dimensional, two-phase, incompressible flow in homogeneous porous media, with zero capillary pressure. A variable cross-section is included so that a variable flow velocity may be considered. The basic differential equations are (1) (2) The total volumetric flow rate is given by (3) Addition of Eqs. 1 and 2 yields =O SPEJ P. 79


Sign in / Sign up

Export Citation Format

Share Document