Torsional Vibrations of Synchronous Motor Driven Trains Using p-Method

1995 ◽  
Vol 117 (1) ◽  
pp. 152-160 ◽  
Author(s):  
W. J. Chen

Synchronous motors produce damaging oscillating torques during startup. If the system is not properly analyzed and designed, the torsional excitation can be very destructive. This paper presents a systematic approach to the dynamic analysis of synchronous motor driven rotating machinery. The p-version of the finite element method is used in the formulation of the equations of motion which provides a great deal of simplicity in the modeling process. The convergence is achieved by increasing the polynomial order of the basis functions of the geometric elements. The system damping matrix can be constructed from the element level or can be calculated by specifying the critical damping factors for a given number of modes of interest. A modified Newmark integration method is employed in the nonlinear transient response calculation. The nonlinearity of the flexible resilient couplings can be easily implemented into this direct numerical integration algorithm. The dynamic stiffness and damping of the resilient couplings are updated at each time step to ensure the dynamic equilibrium. Two examples have been employed to illustrate the validity of the proposed algorithm. The effectiveness, accuracy, and simplicity of the use of p-method on the torsional vibration of synchronous motor driven trains are demonstrated in this paper.

2012 ◽  
Vol 134 (1) ◽  
Author(s):  
M. Amine Hassini ◽  
Mihai Arghir

The traditional small perturbation method is successfully used for linear dynamic analysis of gas bearings but excludes any nonlinear study. Investigating large displacements requires the evaluation of the nonlinear aerodynamic forces in the thin film. To avoid solving the unsteady compressible thin film fluid equations, we propose a method based on the use of frequency dependent dynamic coefficients and on the rational function approximation of the resulting impedances. Calculating impedances for several eccentricities enables mapping the full dynamic behavior of the bearing. A set of ordinary differential equations is then developed by using the inverse of Laplace transform. The equations of motion of the rotor are subsequently solved numerically with local linearization at each time step. The numerical results obtained by using impedances are in good agreement with the reaction forces obtained by solving the full nonlinear transient Reynolds equation.


Author(s):  
Paulo Flores

The main purpose of this work is to present a general and comprehensive approach to automatically adjust the time step for the contact and non contact periods in multibody dynamics. The basic idea of the described methodology is to ensure that the first impact within a multibody system does not occur with a large value for relative bodies’ penetration in order to avoid the artificially large contact forces associated. The detection of the instant of contact takes place when the distance between two bodies change the sign between two discrete moments in time. In fact, in theory, the contact starts when this distance is zero, or a very small value to prevent the round-off errors. Thus, during the numerical solution of the system equations of motion if the first penetration is below this small value previously specified, then the current time is taken as the impact time. On the other hand, if the first penetration is larger than the specified tolerance, then the current time step is beyond the impact time. In this case, integration algorithm is forced to go back and take a smaller time step until a step can be taken within the acceptable tolerance. The main features of this approach are the easiness to implement and the good computational efficiency. In addition, it can easily deal with the transitions between non contact and contact cases in multibody dynamics. Finally, results obtained from dynamic simulations are presented and discussed to study the validity of the methodology proposed in this work.


Author(s):  
Jang-Der Jeng ◽  
Yuan Kang ◽  
Yeon-Pun Chang ◽  
Shyh-Shyong Shyr

The Duffing oscillator is well-known models of nonlinear system, with applications in many fields of applied sciences and engineering. In this paper, a response integration algorithm is proposed to analyze high-order harmonic and chaotic motions in this oscillator for modeling rotor excitations. This method numerically integrates the distance between state trajectory and the origin in the phase plane during a specific period and predicted intervals with excitation periods. It provides a quantitative characterization of system responses and can replace the role of the traditional stroboscopic technique (Poincare´ section method) to observe bifurcations and chaos of the nonlinear oscillators. Due to the signal response contamination of system, thus it is difficult to identify the high-order responses of the subharmonic motion because of the sampling points on Poincare´ map too near each other. Even the system responses will be made misjudgments. Combining the capability of precisely identifying period and constructing bifurcation diagrams, the advantages of the proposed response integration method are shown by case studies. Applying this method, the effects of the change in the stiffness and the damping coefficients on the vibration features of a Duffing oscillator are investigated in this paper. From simulation results, it is concluded that the stiffness and damping of the system can effectively suppress chaotic vibration and reduce vibration amplitude.


Author(s):  
SD Yu ◽  
BC Wen

This article presents a simple procedure for predicting time-domain vibrational behaviors of a multiple degrees of freedom mechanical system with dry friction. The system equations of motion are discretized by means of the implicit Bozzak–Newmark integration scheme. At each time step, the discontinuous frictional force problem involving both the equality and inequality constraints is successfully reduced to a quadratic mathematical problem or the linear complementary problem with the introduction of non-negative and complementary variable pairs (supremum velocities and slack forces). The so-obtained complementary equations in the complementary pairs can be solved efficiently using the Lemke algorithm. Results for several single degree of freedom and multiple degrees of freedom problems with one-dimensional frictional constraints and the classical Coulomb frictional model are obtained using the proposed procedure and compared with those obtained using other approaches. The proposed procedure is found to be accurate, efficient, and robust in solving non-smooth vibration problems of multiple degrees of freedom systems with dry friction. The proposed procedure can also be applied to systems with two-dimensional frictional constraints and more sophisticated frictional models.


Author(s):  
Kai Feng ◽  
Xueyuan Zhao ◽  
Zhiyang Guo

With increasing need for high-speed, high-temperature, and oil-free turbomachinery, gas foil bearings (GFBs) have been considered to be the best substitutes for traditional oil-lubricated bearings. A multi-cantilever foil bearing (MCFB), a novel GFB with multi-cantilever foil strips serving as the compliant underlying structure, was designed, fabricated, and tested. A series of static and dynamic load tests were conducted to measure the structural stiffness and equivalent viscous damping of the prototype MCFB. Experiments of static load versus deflection showed that the proposed bearing has a large mechanical energy dissipation capability and a pronounced nonlinear static stiffness that can prevents overly large motion amplitude of journal. Dynamic load tests evaluated the influence of motion amplitude, loading orientation and misalignment on the dynamic stiffness and equivalent viscous damping with respect to excitation frequency. The test results demonstrated that the dynamic stiffness and damping are strongly dependent on the excitation frequency. Three motion amplitudes were applied to the bearing housing to investigate the effects of motion amplitude on the dynamic characteristics. It is noted that the bearing dynamic stiffness and damping decreases with incrementally increasing motion amplitudes. A high level of misalignment can lead to larger static and dynamic bearing stiffness as well as to larger equivalent viscous damping. With dynamic loads applied to two orientations in the bearing midplane separately, the dynamic stiffness increases rapidly and the equivalent viscous damping declines slightly. These results indicate that the loading orientation is a non-negligible factor on the dynamic characteristics of MCFBs.


1997 ◽  
Vol 119 (1) ◽  
pp. 57-63 ◽  
Author(s):  
M. J. Goodwin ◽  
P. J. Ogrodnik ◽  
M. P. Roach ◽  
Y. Fang

This paper describes a combined theoretical and experimental investigation of the eight oil film stiffness and damping coefficients for a novel low impedance hydrodynamic bearing. The novel design incorporates a recess in the bearing surface which is connected to a standard commercial gas bag accumulator; this arrangement reduces the oil film dynamic stiffness and leads to improved machine response and stability. A finite difference method was used to solve Reynolds equation and yield the pressure distribution in the bearing oil film. Integration of the pressure profile then enabled the fluid film forces to be evaluated. A perturbation technique was used to determine the dynamic pressure components, and hence to determine the eight oil film stiffness and damping coefficients. Experimental data was obtained from a laboratory test rig in which a test bearing, floating on a rotating shaft, was excited by a multi-frequency force signal. Measurements of the resulting relative movement between bearing and journal enabled the oil film coefficients to be measured. The results of the work show good agreement between theoretical and experimental data, and indicate that the oil film impedance of the novel design is considerably lower than that of a conventional bearing.


Author(s):  
Jesús Cardenal ◽  
Javier Cuadrado ◽  
Eduardo Bayo

Abstract This paper presents a multi-index variable time step method for the integration of the equations of motion of constrained multibody systems in descriptor form. The basis of the method is the augmented Lagrangian formulation with projections in index-3 and index-1. The method takes advantage of the better performance of the index-3 formulation for large time steps and of the stability of the index-1 for low time steps, and automatically switches from one method to the other depending on the required accuracy and values of the time step. The variable time stepping is accomplished through the use of an integral of motion, which in the case of conservative systems becomes the total energy. The error introduced by the numerical integrator in the integral of motion during consecutive time steps provides a good measure of the local integration error, and permits a simple and reliable strategy for varying the time step. Overall, the method is efficient and powerful; it is suitable for stiff and non-stiff systems, robust for all time step sizes, and it works for singular configurations, redundant constraints and topology changes. Also, the constraints in positions, velocities and accelerations are satisfied during the simulation process. The method is robust in the sense that becomes more accurate as the time step size decreases.


Author(s):  
Sung-Soo Kim ◽  
Jeffrey S. Freeman

Abstract This paper details a constant stepsize, multirate integration scheme which has been proposed for multibody dynamic analysis. An Adams-Bashforth Moulton integration algorithm has been implemented, using the Nordsieck form to store internal integrator information, for multirate integration. A multibody system has been decomposed into several subsystems, treating inertia coupling effects of subsystem equations of motion as the inertia forces. To each subsystem, different rate Nordsieck form of Adams integrator has been applied to solve subsystem equations of motion. Higher order derivative information from the integrator provides approximation of inertia force computation in the decomposed subsystem equations of motion. To show the effectiveness of the scheme, simulations of a vehicle multibody system that consists of high frequency suspension motion and low frequency chassis motion have been carried out with different tire excitation forces. Efficiency of the proposed scheme has been also investigated.


Author(s):  
Ahmad W. Yacout

This study has theoretically analyzed the surface roughness, centripetal inertia and recess volume fluid compressibility effects on the dynamic behavior of a restrictor compensated hydrostatic thrust spherical clearance type of bearing. The stochastic Reynolds equation, with centripetal inertia effect, and the recess flow continuity equation with recess volume fluid compressibility effect have been derived to take into account the presence of roughness on the bearing surfaces. On the basis of a small perturbations method, the dynamic stiffness and damping coefficients have been evaluated. In addition to the usual bearing design parameters the results for the dynamic stiffness and damping coefficients have been calculated for various frequencies of vibrations or squeeze parameter (frequency parameter) and recess volume fluid compressibility parameter. The study shows that both of the surface roughness and the centripetal inertia have slight effects on the stiffness coefficient and remarkable effects on the damping coefficient while the recess volume fluid compressibility parameter has the major effect on the bearing dynamic characteristics. The cross dynamic stiffness showed the bearing self-aligning property and the ability to oppose whirl movements. The orifice restrictor showed better dynamic performance than that of the capillary tube.


Sign in / Sign up

Export Citation Format

Share Document