Dynamic characteristics of self-acting gas bearing flexible rotor coupling system based on the forecasting orbit method
|
|
- Melvyn Carroll
- 5 years ago
- Views:
Transcription
1 Nonlinear Dyn (2012) 69: DOI /s z Dynamic characteristics of self-acting gas bearing flexible rotor coupling system based on the forecasting orbit method Guang-hui Zhang Yi Sun Zhan-sheng Liu Min Zhang Jia-jia Yan Received: 26 November 2010 / Accepted: 4 November 2011 / Published online: 3 December 2011 Springer Science+Business Media B.V Abstract This paper studies the nonlinear dynamic characteristics of a flexible rotor supported by selfacting gas bearings theoretically. The multiple degree freedom model of flexible rotor is established by the finite element method and analyzed coupled with the transient gas lubricated Reynolds equation by employing the forecasting orbit method. The Reynolds equation is solved by the alternating direction implicit method and the dynamic response of the rotor is calculated by the Newmark integral method. To settle the problem that the two kinds of transient solving processes (transient Reynolds equation for bearing and transient equation of motion for rotor) cannot be solved simultaneously, which arises from the fact that they need each other s results as their initial values, the multi-field coupling algorithm based on the forecasting method is proposed and applied in this paper. By employing the numerical method, the rotor trajectory diagram, phase diagram, frequency spectrum, power spectrum, bifurcation diagram, and vibration mode diagram were obtained. It is to note that the dynamic characteristics of self-acting gas G.-h. Zhang ( ) Z.-s. Liu M. Zhang J.-j. Yan School of Energy Science and Engineering, Harbin Institute of Technology, Harbin, China qq @gmail.com Y. Sun School of Astronautic, Harbin Institute of Technology, Harbin, China bearing rotor system and whirling instability of the system could be depicted successfully. This would establish the foundation for contributing to a further understanding of the gas bearing flexible rotor system. Keywords Gas bearing rotor system Stability Nonlinear dynamic Bifurcation Forecasting orbit method Nomenclature c average air film thickness e eccentricity e x,e y eccentricity in x, y direction h air film thickness l T the length of the rotor element p air film pressure p a atmospheric pressure t time x,y,z coordinates A the section area of the rotor F x,f y the gas film forces act on the rotor F n,f n+1 the gas film forces at time n and n + 1 H dimensionless air film thickness H n gas film thickness function at time n H n+1 forecast gas film thickness function at time n + 1by forecasting I y moment of the inertia for rotor L the length of the bearing NX the mesh number of circumferential direction
2 342 G.-h. Zhang et al. NZ the mesh number of bearing length direction P dimensionless pressure P n+1 dimensionless pressure at time n + 1 Q the square of the P ; Q = P 2 R radius of the bearing W load capacity ε the eccentricity ratio η the dynamic viscosity of air θ the angular coordinate θ 0 attitude angle ξ the dimensionless coordinates in length direction ρ the air density under arbitrary pressure ρ a the air density under atmospheric pressure ρ r the density of the rotor τ the dimensionless time ϕ the attitude angle ω rotating angular speed of the shaft Λ the bearing number shear factor for the rotor Φ s 1 Introduction The gas bearings are more and more widely used in aerospace and mechanical engineering fields for the benefit of high rotating speeds, low friction, no wear, no contamination and smooth functioning, which are generated from the features of gas-lubricated film. However, the self-excited whirling phenomenon exists in the gas bearing which is similar to that of liquid bearing. The whirling speed of gas bearing is 0.45 times of rotating speeds [1], which will cause the instability of gas bearing rotor system and limit the application of gas bearing. So for the proper design of such systems, the dynamic of gas bearing rotor system should be identified accurately and effectively. Although abundant literature for predicting the nonlinear whirling behavior of gas bearing can be found, most of it has limitations: some ignore the partial derivative of pressure or film thickness with respect to time for gas film force; some simplify the rotor as the rigid kinematic model; and nearly none takes into account the coupling mechanism between the gas film forces and the rotor motion. Those factors should be considered simultaneously to obtain the dynamic characteristics of gas bearing rotor system accurately and effectively. Usually there are two ways to analyze the dynamic characteristics of gas bearing rotor system: one is the perturbation method based on the steady-state operations, the other one is the orbit method. In the perturbation method, the load on perturbation state is calculated by applying small displacement or velocity on the steady-state operation; the variation of load is decomposed by perturbation method and then the stiffness and damping coefficients are obtained. This method was employed by Sternlicht [2], Rentzepis and Sternlicht [3], Ausman [4], and Lund [5] and has the advantage of high maneuverability and accuracy on analyzing the smaller whirling orbit (for example, the lower rotating angular velocity and higher load); however, for the case of high rotating velocity and light rotor weight, where the whirling orbit of the rotor is relatively bigger, the difficult-to-estimate numerical error will be generated by applying the perturbation method. Another way to solve the stability problem of gas bearing rotor system is the orbit method, which was proposed by Castelli in Ref. [6]. In the orbit method, the complete nonlinear equations are included, the calculation procedures highly coincide with the assumed governing equations, and the trajectory orbit of the rotor can be obtained. So the orbit method enables both establishment of the stability threshold and prediction of the behavior of the bearing into the instability region. Because of the above-mentioned advantages, the orbit method is being more and more widely employed. Castelli discussed the numerical method for solving the transient Reynolds equation in Ref. [7]. The proposed method comprises the satisfactory numerical stability of implicit arithmetic and the fast convergence of the explicit method. In Ref. [8] Dimofte modified the alternating direction implicit method and developed the code for analyzing the steady and dynamic characteristics. Bou-Sa solved the nonlinear Reynolds equation with time terms to simulate the dynamic characteristics of the rotor in Ref. [9]. The proposed model was compared with conventional linear journal bearing model which employs the stiffness and damping coefficients. It is to note that the linear model would lead to the incorrect solutions; the nonlinear model could forecast the rotor dynamic characteristics better. In Refs. [10 12], the gas-lubricated Reynolds equation with time term was solved by employing the finite difference method and over-relaxation method. The dynamic behavior in horizontal and vertical directions of the journal was analyzed by locus diagram, Poincaré map, power spectrum and bifurcation diagram. The concept of whirling ratio was imported into
3 Dynamic characteristics of self-acting gas bearing flexible rotor coupling system 343 the solving procedure. Some limits existed in the research for journal gas bearing rotor system mentioned above. Most of the analysis was particular for bearing without considering the effect of the rotor; some assumed the rotor model as rigid rotor with little number of freedoms, but nearly no-one modeled the multiple degrees of freedom for gas bearing rotor system. Also most of the papers employed the last-step results to solve the next step, without considering the coupling between the transient Reynolds equation for bearing and transient equation of motion for rotor. The essence of the rotor gas bearing coupling dynamic system is the fluid structure interaction problem. Because of the characteristics of different physical fields, the numerical methods for solving them are distinct. Many different coupling algorithms are proposed. The full coupling algorithm which is called the monolithical approach is proposed to solve the fluid structure interaction problem in Refs. [13 16]. The full coupling algorithm has the advantages of good stabilization and introducing no energy error. But this approach cannot employ the originally developed flow domain and structure domain simulation code, which limits the application of the monolithical approach. With the status of the sophisticated flow domain and structure domain calculation, it is better to develop the proper multi-field algorithm based on the interactive relation. Piperno and Farhat in Refs. [17, 18] developed the partitioned procedures for the twodimensional and three-dimension aeroelastic problems. Farhat in Ref. [19] proposed the structure forecasting predictor that is based on the energy conservation law to simulate the aeroelastic problems. The choice of the structures predictor has to do with the calculation form of the flow domain. Usually it is difficult to choose the structures predictor. Also, for the realistic cases, the exactly energy conservation is impossible. Piperno [20] analyzed the predictor, and obtained the same conclusion. Since there are few references regarding the coupling between the gas bearing flow and the rotor motion this topic needs to be studied in depth. The self-acting gas bearing rotor system with multiple degrees of freedom is modeled by finite element method in this paper. The ADI (alternating direction implicit) is employ to solve the transient gas lubricated Reynolds equation and the Newmark method is used to solve the equations of motion for rotor. The forecasting orbit method is proposed to couple the gaslubricated Reynolds equation with the equation of motion for rotor. The nonlinear dynamic characteristics of journal gas bearing rotor system are analyzed, and the phenomenon and mechanism of the whirling in the gas bearing are studied. The effects of the imbalance mass on the dynamic behavior are investigated in this paper. 2 Mathematical modeling 2.1 Modeling for transient gas bearing forces The mathematical model of gas bearings is much more complicated to deal with than that of the liquid film bearing because of the compressibility of the gas. So we need the following assumptions: the variation of viscosity and density across the film direction is neglected, the flow is isothermal, the mass flow loss caused by the end leakage of the bearing is neglectable, and the air is an ideal gas. The well-known transient Reynolds equation of an isothermal gas journal bearing is as follows: ( ρh 3 ( ρh 3 ) z ) + 12 t (ρh) p p R θ η R θ η z + 6ω (ρh) = 0 (1) θ The film thickness is as follows: h = c ( 1 + ε cos(θ ϕ) ) (2) ε = 1 (e x ) c 2 + (e y ) 2 ; ϕ = arctg e y (3) e x By employing the dimensionless parameters, ξ = z R, P = p p a, H = h c, Λ = 6ηωR2 p a c 2, τ = ωt 2Λ The dimensionless form of transient Reynolds equation is derived by introducing the isothermal gas condition p/ρ = p a /ρ a into (1): ( PH 3 P ) ( PH 3 P ) θ θ ξ ξ + Λ (PH) (P H ) + = 0 (4) θ τ Equation (4) is the nonlinear partial differential equation, so the iteration method is employed to solve the equation.
4 344 G.-h. Zhang et al. Fig. 1 The model of self-acting gas journal bearings Fig. 2 Calculation domain of journal gas bearings (self-acting type) The structure schematic of self-acting gas journal bearing is shown in Fig. 1. The symmetry boundary condition, periodic boundary condition, mass flow boundary condition and atmospheric boundary condition are defined in Fig. 2. The computational domain is θ(0 : NX+1) and ξ(0 : NZ+1). The solution domain is θ(1 : NX) and ξ(1 : NZ). The atmospheric boundary condition is p(ξ = 0) = p(ξ = NZ)= p a ; the periodic boundary condition is p(θ = 0) = p(θ = NX) and p(θ = 1) = p(θ = NX+ 1). By submitting Q = P 2, the central difference method is employed to separate the variables in the θ and ξ direction, and the ADI (alternating direction implicit) method is used to separate the variables in the form of the time τ. According to the rules governing the scheme, the computations are carried out for three discrete time levels, n, n + 1 and n + 2, which are separated from each other by segments of the dimensionless time τ. One iterative cycle consists of solving procedure for the levels n + 1 and n + 2. In the beginning of the computation, the pressure distribution in the level n is assumed to be known. In the first stage of the iterative cycle, for the level n + 1, the derivatives with respect to θ are written as unknowns, and the derivatives with respect to ξ are considered to be the values known from the level n. In the second stage (level n + 2), the derivatives with respect to ξ are unknown, and the derivatives with respect to θ are treated as the known values from the level n + 1. On finishing the iterative cycle, the level n + 1 is treated as a new level n and the next cycle is started. In this scheme the following formulas are applied: h n p n Q n+1 Q n + 2 p n τ + 2Λ Qn+1 p n h n h n +1 h n 1 2 θ h n+1 h n τ Q n Λ Qn+1 1 p n 2 θ 3 ( h n ) 2 ( h n +1 h n 1 )(Qn+1 +1 Qn+1 1 ) 4 θ 2 3 ( h n ) 2 ( h n i+1,j h n i, 1j )(Qn i+1,j Qn i 1,j ) 4 ξ 2 ( h n ) 3 Q n Qn+1 + Q n+1 1 θ 2 ( h n ) 3 Q n i+1,j 2Qn + Qn i 1,j ξ 2 = 0 (5), Qn+1 1 are the unknown param- where Q n+1 eters for step 1. +1, Qn+1 h n+1 Q n+2 Q n+1 p n+1 τ + 2Λ Qn+1 p n p n+1 h n+1 +1 h n θ h n+2 h n+1 τ h n+1 Q n Λ Qn+1 1 p n+1 2 θ 3 ( h n+1 ) 2 ( h n+1 +1 h n+1 1 )(Qn+1 +1 Qn+1 1 ) 4 θ 2 3 ( h n+1 ) 2 ( h n+1 i+1,j h n+1 i, 1j )(Qn+2 i+1,j Qn+2 i 1,j ) 4 ξ 2 ( h n+1 ) 3 Q n Qn+1 + Q n+1 1 θ 2 ( h n+1 ) 3 Q n+2 i+1,j 2Qn+2 + Q n+2 i 1,j ξ 2 = 0 (6)
5 Dynamic characteristics of self-acting gas bearing flexible rotor coupling system 345 EI y K = lt 3 (1 + Φ s) 12 6l T 12 6l T (4 + Φ s )lt 2 6l T (2 Φ s )l 2 T 12 6l T symm. (4 + Φ s )lt 2 Fig. 3 A model of gas journal bearing system where Q n+2 i+1,j, Qn+2, Qn+2 i 1,j are the unknown parameters for step 2. Equations (5) and (6) are the nonlinear equations, so the ADI method is not unconditionally stable. To keep the stability of the calculation scheme, the condition following should be considered: the step time τ should be less than During the solving procedure, the over-relaxation method is employed to form and update the triangular matrix and the discrete equations are solved. By integrating the transient gas film pressure p(θ, ξ,τ) acting on the rotor, the transient gas film forces can be obtained as follows: 2π L/R F x = P a R 2 p(θ,ξ,τ)cos θdξdθ; (7) 0 F y = P a R 2 2π 0 0 L/R 0 p(θ,ξ,τ)sin θdξdθ (8) 2.2 Modeling for motion equation of elastic rotor The model of journal gas bearing system is presented in Fig. 3. The rotor system is formed by discrete disk, elastic shaft segment with distributed mass. The numbering scheme is shown in Fig. 3, where the journal gas bearing locates on nodes 2, 6 and the imbalance mass locates on node 4. The finite element method is employed to form the stiffness and mass matrixes. The Timoshenko beam element is used. Then the equation of motion is as follows: M{ü}+C{ u}+k{u}={f } (9) where {u} is the displacement and angular vector of each node for the rotor system is {x 1,θ y1,y 1, θ x1, x 2,θ y2,y 2 θ x2,...,x N,θ yn,y N, θ xn }. N is the total number of nodes of the rotor system. K is the stiffness matrix, and M is the mass matrix. By employing the strain energy and kinematic energy expression, the above matrix is as follows: M = M T + M R ρal T M T = 420(1 + Φ s ) 2 m 1 l T m 2 m 3 l T m 4 l T 2 m 5 l T m 4 lt 2 m 6 m 1 l T m 2 symm. lt 2 m 5 ρi y M R = 30l T (1 + Φ s ) 2 m 7 l T m 8 m 7 l T m 8 l T 2 m 9 l T m 8 lt 2 m 10 m 7 l T m 8 symm. lt 2 m 9 m 1 = Φ s + 140Φ 2 s ; m 2 = Φ s Φ 2 s ; m 3 = Φ s + 70Φ 2 s ; m 4 = Φ s Φ 2 s ; m 5 = 4 + 7Φ s + 3.5Φ 2 s ; m 6 = 3 + 7Φ s + 3.5Φ 2 s ; m 7 = 36; m 8 = 3 15Φ s ; m 9 = 4 + 5Φ s + 10Φ 2 s ; m 10 = 1 + 5Φ s 5Φ 2 s The gyroscopic matrix is equal to 2M R ; C is the damping matrix, which includes the gyroscopic matrix; {F } are the forces acting on the rotor. In this paper, {F } include the transient gas film forces F x and F y which act on the elastic rotor. Equation (9) is the nonlinear ordinary differential equation; the step-by-step integration method is applied to carry out the numerical simulation. The Newmark integration method is employed to calculate the displacement, velocity and acceleration for each node of the rotor.
6 346 G.-h. Zhang et al. Fig. 4 The flow chart for weak coupling orbit method 2.3 Solving process of forecasting orbit method The orbit method based on weak coupling process The weak coupling process between the rotor dynamic response algorithm and transient gas film forces algorithm is depicted as follows: as shown in Fig. 4, the initial condition for the loop is the rotor initial position in the bearing at time n, expressed by the initial gas film thickness function H n. In weak coupling algorithm, the gas film thickness function H n+1 forecast at time n + 1 is equal to the gas film thickness function at time n, and then it will be taken into the ADI algorithm as the initial condition to solve the transient gas lubricated Reynolds equation at time n + 1, where the transient pressure distribution will be obtained. By integrating the gas pressure distribution P n+1, the transient gas film forces acting on the rotor will be obtained. With the obtained forces, the Newmark method will be employed to calculate the dynamic response of the rotor system. The orbit trajectory of the rotor at time n + 1 expressed by the gas film thickness function H n+1 will be obtained and then taken as the initial value for the next loop. From the above weak coupling process, the rotor dynamic equation is the ordinary differential equation and the gas-lubricated Reynolds equation is the second-order partial differential equation, which comprise the two analysis steps for one step of the bearing rotor system solving procedure. Solving results of the two types of analyses are needed as each other s initial value, which cannot be realized simultaneously in the coupling process. For this problem, the dynamic orbit method based on the weak coupling process is being proposed. In this procedure, the gas film thickness function from the rotor dynamic response calculation at last time step will be used as initial value for transient gas lubricated Reynolds equation for the next time step. Because of the delay of the time step of the two kinds of equation, this method is called the weak coupling process The orbit method based on forecasting coupling process The forecasting coupling process between the rotor dynamic response algorithm and transient gas film forces algorithm is depicted as follows: as shown in Fig. 5, based on the weak coupling process, the initial condition for the loop is also the rotor initial position in the bearing at time n, expressed by the initial gas film thickness function H n. Differently from the weak coupling algorithm, the gas film thickness function H n+1 forecast at time n + 1 will be calculated by forecasting method in forecasting coupling procedure. Then it will be taken into the ADI algorithm as the initial condition to solve the transient gas lubricated Reynolds equation at time n + 1, where the transient pressure distribution will be obtained. By integrating the gas pressure distribution P n+1, the transient gas film forces acting on the rotor will be obtained. With the obtained forces, the Newmark method will be employed to calculate the dynamic response of the rotor system. The orbit trajectory of the rotor at time n + 1 expressed by the gas film thickness function H n+1 will be obtained and then taken as the initial value for the next loop.
7 Dynamic characteristics of self-acting gas bearing flexible rotor coupling system 347 Fig. 5 The flow chart for forecasting coupling orbit method There are many ways to predict the gas film thickness function H n+1 forecast at time n + 1. Different methods need a varied number of the gas film thickness at previous time. The linear interpolation method is employed in this study to forecast he gas film thickness function at the next time step. It assumes that the trend of the function varying with time is linear, and then the following equation will be listed: H n+1 H n τ n = H n H n 1 τ n 1 and then ( H H n+1 n H n 1 = τ n 1 + H n ) τ n τ n So it is obvious that the gas film thickness function of the two previous time steps is needed for the linear forecasting procedure. 3 Numerical simulation and dynamic analysis The simulation model for numerical analysis is as follows: the length of the bearing L = 0.05 m, the diameter of the bearing D = m, air density ρ a = kg/m 3, air viscosity μ = 1.82E 5 N s/m 2, specific heat ratio k = 1.401, atmospheric pressure p a = Pa, the mesh grid size: in z direction (axial direction) nz = 20, in x direction (circumferential direction) nx = 40, the length of the rotor: 0.3 m, the diameter of the rotor: 0.05 m, the elastic modulus of the rotor 2.0E11 Pa, the density of the rotor: 7825 kg/m 3. The convergence condition is Max p n+1 p n Table 1 Comparison of calculated results Λ ε W θ 0 FEM DQM FEM DQM Verification of the proposed method To verify the accuracy of the orbit method employed in this study, the comparison is presented with the results from Ref. [21]bymodelingthe bearingof length diameter equal to 1. By ignoring the effects of rotor and the transient term of the gas lubricated Reynolds equation, the load capacity W and attitude angle θ 0 are obtained as shown in Table 1. The columns with are the results of the proposed method. It indicates that the numerical results of this paper are in good agreement with the reference results, which makes the foundation for the analysis of the gas bearing rotor system. 3.2 Comparison between the weak coupling and forecasting coupling methods The comparison between the weak coupling orbit method and the forecasting coupling orbit method is
8 348 G.-h. Zhang et al. Fig. 6 Shaft locus of different rotating speeds for different coupling algorithms without imbalance mass presented in this paper. The time step in this analysis is s. The effect of the shaft locus obtained by weak coupling and forecasting coupling methods are observed in Fig. 6, where the imbalance mass is ignored. As shown in Fig. 6(a,b), when the rotating speed is 5000 r/min, the difference between the two method is neglectable. As the rotating speed increases to 11,000 r/min, it indicates that the shaft locus amplitude of the forecasting method is a bit smaller than that of the weak coupling method in Fig. 6(c,d), but the same trend exists. As the rotating speed reaches 14,000 r/min, the regularity is the same as for rotating speed of 11,000 r/min. It is observed that in the case of ignoring the imbalance mass, the tendency for the weak coupling and forecasting coupling methods is the same, but the amplitude of the shaft locus for weak coupling is bigger. The effect of the shaft locus obtained by weak coupling and forecasting coupling method is observed in Fig. 7, where the imbalance mass is considered. AsshowninFig.7(a,b), when the rotating speed is 5000 r/min, the difference between the two methods is small. As the rotating speed increases to 11,000 r/min, it indicates that the shaft locus amplitude of the forecasting method is a bit smaller than that of the weak coupling method in Fig. 7(c,d), but the same trend exists. As the rotating speed reaches 14,000 r/min, the regularity is the same as for rotating speed of 11,000 r/min. It is observed that in the case of considering the imbalance mass, the tendency for the weak coupling and forecasting coupling methods is the same, but the amplitude of the shaft locus for weak coupling is bigger. The case which considers the effect of gravity forces and ignores the imbalancing mass is presented in Fig. 8. It indicates that for different rotating speeds (5000, 11,000 and 14,000 r/min), the tendency of the shaft locus for the two coupling algorithms is same. The static balance positions for the two methods are also same. In is observed that the weak coupling method and forecasting coupling method have little influence on the dynamic response when the gravity forces are considered. The variation of iterative number for different time steps which reach the stable condition is shown in Fig. 9, where the step times , , , , and s are employed. The calculated number of steps de-
9 Dynamic characteristics of self-acting gas bearing flexible rotor coupling system 349 Fig. 7 Shaft locus of different rotating speeds for different coupling algorithms with imbalance mass Fig. 8 Shaft locus of different rotating speeds for different coupling algorithms with imbalance under gravity force
10 350 G.-h. Zhang et al. creases as the time step increases, but the total time of reaching the stable condition increases. The reason for this phenomenon is that the essential asynchronous characteristics of weak coupling method introduce the delay error. For the cases of forecasting coupling method, the tendency for decrement of calculated step is more intensive than that of weak coupling method. It is observed that as the time step increases, the iterative number decreases, and the prod- uct of the two nearly keeps unchanged. The total time of reaching the stable condition for forecasting coupling method does change with the variation of the time step. So the forecasting coupling method can eliminate the delay of the weak coupling algorithms effectively. In the following study, the forecasting coupling method will be used to analyze the dynamics of the gas bearing rotor system. 3.3 The dynamic analysis of self-acting gas bearing rotor system The dynamic analysis without imbalance mass Fig. 9 The variation of iterative number with time steps The shaft locus and time response spectra diagrams are presented in Fig. 10 by analyzing the rotor bearing system without imbalance mass. As shown in Fig. 10(a), when the rotating speed is 5000 r/min, the dynamic characteristics of the rotor bearing system is point stabilization, where the rotor is running near some stable position after initial motion. When the rotating speed reaches 11,500 r/min (Fig. 10(b)), the Fig. 10 The shaft locus, time response and frequency spectrum diagram for different rotating speeds with no imbalance mass
11 Dynamic characteristics of self-acting gas bearing flexible rotor coupling system 351 Fig. 11 The power spectrum, phase diagram and Poincaré maps for different rotating speeds with no imbalance mass stable limit cycle motion appears. With the increment of the rotating speed as in Fig. 10(c), the motion of the rotor does not attenuate, and keeps the stable limit cycle motion status. The phase diagram in X direction in Fig. 11 also indicates the existence of the limit cycle motion. The frequency analysis is performed and the results are presented in Fig. 10(d,e,f) and Fig. 11(a,b,c) by ignoring the imbalance mass effect. It indicates that as the rotating speed increases, the dynamic response amplitude for the cases of 11,500 and 14,000 r/min changes little. The lower frequency component that is smaller than 0.5 frequency multiplication appears generally. The range of lower frequency whirl for gas bearing is smaller than that of oil bearing because of the lower viscosity of the gas. For frequency analysis, the data of stable motion are employed, but for the shaft locus the initial data are used. The power spectral diagram is presented in Fig. 11(d,e,f). It is observed that the frequency components of the gas bearing rotor system are complicated, where two lower frequency components appear. The difference between the two lower frequencies is small for the point stable condition (Fig. 4(a)). For the cases of the limit cycle, the frequency components are complicated, and the amplitude of different fre-
12 352 G.-h. Zhang et al. quency components varies with the rotating speeds in power spectrum (Fig. 11(b,c)), which is single frequency in frequency spectrum diagram. The Poincaré maps for different rotating speeds are presented in Fig. 11(g,h,i). In non-autonomous system, a point on the Poincaré section is referred to as the return point of the time series at a constant interval of T, where T is the period of the rotating speed. As the rotating speed is 5000 r/min (Fig. 11(g)), the T period motion is stable. As the rotating speed increased to 11,500 r/min (Fig. 11(h)), the T period motion loses the stability and is replaced by 8T period motion. When the rotating speed reaches 14,000 r/min (Fig. 11(i)), the system continues to exhibit the 8T period motion. The bifurcation diagram of the gas bearing rotor system is shown in Fig. 12, which takes the rotating speed as the bifurcation parameter. It indicates that as the rotating speed increases to 8000 r/min, the Hopf bifurcation occurs. The motion of the system changes from the balance points of stabilization to the limit cycle condition The dynamic analysis with imbalance mass Fig. 12 The x direction displacement bifurcation with no imbalance mass The cases that consider the imbalance mass are carried out. The imbalance mass acts on the node 4. The mass weight is 0.2 kg, and the eccentricity is m. Fig. 13 The shaft locus, time response and frequency spectrum diagram for different rotating speeds with imbalance mass
13 Dynamic characteristics of self-acting gas bearing flexible rotor coupling system 353 Fig. 14 The power spectrum, phase diagram and Poincaré maps for different rotating speeds with imbalance mass When the rotating speeds are 5000 and 10,000 r/m, the rotor runs with the status of stable periodicity which is shown in Fig. 13(a,b). The amplitude of the motion increases as the rotating speed increases. When the rotating speed reaches 3000 r/min, the stable periodic motion destabilizes. From the phase diagram of Fig. 14(d,e,f), it indicates that the behavior of the motion becomes more complicated. Then the frequency analysis for the gas bearing rotor system with imbalance mass is carried out. When the rotating speed is 5000 r/min, the basic frequency generated from the imbalance mass is the main component, as shown in Fig. 13(d). As the rotating speed increases to 10,000 r/min, the 0.35 frequency multiplication that the amplitude is small appears in Fig. 13(e). As the rotating speed reaches 30,000 r/min, the am-
14 354 G.-h. Zhang et al. plitude of the 0.35 frequency multiplication increases and then exceeds the value of the basic frequency component. The lower frequency whirling phenomenon is obvious. As the rotating speed increases, the lower frequency components that are smaller than the 0.35 frequency multiplication appear based on the original 0.35 frequency multiplication in Fig. 13(e). But the amplitude of 0.35 frequency multiplication is the biggest in the system. The power spectral analysis is performed for different rotating speeds, which is shown in Fig. 14(a,b,c). It is observed that the frequency component of the power spectrum is more complicated than that of the no imbalance mass. With the increment of the rotating speed, the power of the basic frequency is becoming smaller. On the contrary, the lower frequency component is the main part of the whole frequency component. And also some certain frequency components that exist between the basic and second frequency components appear. The Poincaré maps for different rotating speeds with imbalance mass are presented in Fig. 14(g,h,i). As the rotating speed is 5000 r/min (Fig. 14(g)), the T period motion is stable. As the rotating speed increases to 10,000 r/min (Fig. 14(h)), the system continues to exhibit the T period motion. When the rotating speed reaches 30,000 r/min (Fig. 14(i)), the T period motion loses the stability and is replaced by 8T period motion. AsshowninFig.15, the rotating speed is taken as the bifurcation parameter. It is observed that as the rotating speed reaches 11,000 r/min, the stable periodic motion destabilizes. It indicates that the motion changes into the multiply periodic. Compared with the Fig. 12, the imbalance mass can be employed to increase the whirling speed by delaying the occurrence of whirling. Similar conclusion is given in Ref. [22]: The tests previously described indicated that large amplitude self-excited whirl cannot exist if the magnitude of the eccentricity ratio is kept above a certain critical level. By plotting the deformation of the nodes, the modal diagram of the rotor system with rotating speed of 47,000 r/min is presented in Fig. 16. It indicates that the rotor is running with the first bending natural frequency modes. 4 Conclusion This paper has studied the nonlinear dynamic characteristic of self-acting gas bearing flexible rotor sys- Fig. 15 The x direction displacement bifurcation with imbalance mass Fig. 16 The rotor modal diagram with rotating speeds of r/min with no imbalance mass tem by forecasting coupling orbit method. To settle the problem that the two kinds of transient solving processes (transient Reynolds equation for bearing and transient equation of motion for rotor) cannot be solved simultaneously, which arises from the fact that they need each other s results as their initial values, the multi-field coupling algorithm based on the forecasting method is proposed and applied in this paper. By employing the numerical method, the nonlinear gas film forces, the trajectory diagram, phase diagram, frequency spectrum, power spectrum, bifurcation diagram, vibration mode diagram and the equilibrium position diagram were obtained. It is to be noted that the dynamic characteristics of journal gas bearings rotor system and whirling instability of the gas film could be depicted successfully. The results reveal that the presence of a complex dynamic behavior is comprised of periodic, subharmonic responses of the rotor. By comparing the weak coupling and forecasting coupling method, the forecasting coupling method can
15 Dynamic characteristics of self-acting gas bearing flexible rotor coupling system 355 eliminate the delay of the weak coupling algorithms effectively. Regarding the behavior of the rotor bearing system with respect to variation in rotating speed and imbalance mass, the results indicate that the imbalance mass can be employed to increase the whirling speeds by delaying the occurrence of whirling. The method presented in this study provides the foundation for the design of gas bearing flexible rotor system. The results indicate suitable operating conditions that suppress subharmonic motion and prevent instable behavior. Acknowledgements The research presented here was supported by the National Natural Science Foundation of China (Grant No ) and the Fundamental Research Funds for the Central Universities (Grant No. HIT. NSRIF ). The authors are grateful for the support. References 1. Yang, P., Zhu, K.Q., et al.: On the nonlinear stability of self-acting gas journal bearings. Tribol. Int. 42(1), (2009) 2. Sternlicht, B.: Elastic and damping properties of cylindrical journal bearings. J. Basic Eng. 81, (1959) 3. Rentzepis, G.M., Sternlicht, B.: On the stability of rotors in cylindrical journal bearings. J. Basic Eng. 84(3), (1961) 4. Ausman, J.S.: An improved analytical solution for selfacting, gas-lubricated journal bearings of finite length. J. Basic Eng. 83(2), (1961) 5. Lund, J.W.: The stability of an elastic rotor in journal bearings with flexible, damped supports. J. Appl. Mech. 87, (1965) 6. Castelli, V., Elrod, H.G.: Solution of the stability problem for 360 degree self-acting, gas-lubricated bearing. J. Basic Eng. 87(1), (1961) 7. Castelli, V., Stevenson, C.H.: A Semi-implicit Numerical Method for Treating the Time Transient Gas Lubrication Equation. Mechanical Technology Inc, New York, TID (1967) 8. Dimofte, F.: Fast methods to numerically integrate the Reynolds equation for gas fluid films. New York, NASA Technical Memorandum, (1992) 9. Bou-sa, B., Grau, G., Iordanoff, I.: On nonlinear rotor dynamic effects of aerodynamic bearings with simple flexible rotors. J. Eng. Gas Turbine Power 130(1), (2008) 10. Wang, C., Jang, M., Yeh, Y.: Bifurcation and nonlinear dynamic analysis of a flexible rotor supported by relative short gas journal bearings. Chaos Solitons Fractals 32(2), (2007) 11. Wang, C.C.: Application of a hybrid method to the nonlinear dynamic analysis of a flexible rotor supported by a spherical gas-lubricated bearing system. Nonlinear Anal. 70(5), (2009) 12. Wang, C.C.: Bifurcation and nonlinear analysis of a flexible rotor supported by a relative short spherical gas bearing system. Commun. Nonlinear Sci. Numer. Simul. 15(9), (2010) 13. Frederic, J.B.: A monolithical fluid-structure interaction algorithm applied to the piston problem. Comput. Methods Appl. Mech. Eng. 167, (1998) 14. Michler, C., Hulshoff, S.J., van Brummelen, E.H., et al.: A monolithic approach to fluid-structure interaction. Comput. Fluids 33, (2004) 15. Longatte, E., Verreman, V., Bendjeddou, Z., et al.: Comparison of strong and partitioned fluid structure code coupling methods. In: ASME, Pressure Vessels and Piping Division (Publication) PVP, v 4, Proceedings of the ASME Pressure Vessels and Piping Conference 2005 Fluid-Structure Interaction, PVP2005, pp (2005) 16. Etienne, S., Pelletier, D., Garon, A.: An updated lagrangian monolithic formulation for steady-state fluidstructure interaction problems. In: 43rd AIAA Aerospace Sciences Meeting and Exhibit Meeting Papers, pp (2005) 17. Piperno, S., Farhat, C., Larrouturou, B.: Partitioned procedures for the transient solution of coupled aeroelastic problems Part I: Model problem, theory and twodimensional applications. Comput. Methods Appl. Mech. Eng. 124, (1995) 18. Piperno, S.: Explicit/implicit fluid/structure staggered procedures with a structural predictor and fluid subcycling for 2d inviscid aeroelastic simulations. Int. J. Numer. Methods Fluids 25, (1997) 19. Farhat, C., Lesoinne, C., Letallec, P.: Load and motion transfer algorithms for fluid/structure interaction problems with non-matching discrete interfaces: momentum and energy conservation, optimal discrimination and application to aeroelasticity. Comput. Methods Appl. Mech. Eng. 157(1 2), (1998) 20. Piperno, S.: Partitioned procedures for the transient solution of coupled aeroelastic problems Part II: Energy transfer analysis and three-dimensional applications. Comput. Methods Appl. Mech. Eng. 190(24 25), (2001) 21. Malik, M., Bert, C.W.: Differential quadrature solutions for steady-state incompressible and compressible lubrication problems. J. Tribol. 116(2), (1994) 22. Reynolds, D.B., Gross, W.A.: Experimental investigation of whirl in self-acting air-lubricated journal bearings. Tribol. Trans. 5(2), (1962)
Nonlinear Dynamic Analysis of a Hydrodynamic Journal Bearing Considering the Effect of a Rotating or Stationary Herringbone Groove
G. H. Jang e-mail: ghjang@hanyang.ac.kr J. W. Yoon PREM, Department of Mechanical Engineering, Hanyang University, Seoul, 133-791, Korea Nonlinear Dynamic Analysis of a Hydrodynamic Journal Bearing Considering
More informationSTATIC AND DYNAMIC CHARACTERISTICS OF HYDRODYNAMIC FOUR- LOBE JOURNAL BEARING WITH COUPLE STRESS LUBRICANTS
STATIC AND DYNAMIC CHARACTERISTICS OF HYDRODYNAMIC FOUR- LOBE JOURNAL BEARING WITH COUPLE STRESS LUBRICANTS B. Chetti, b.chetti@gmail.com, Institute of sciences and Technology, Center University of Khemis
More informationEffect of an hourglass shaped sleeve on the performance of the fluid dynamic bearings of a HDD spindle motor
DOI 10.1007/s00542-014-2136-5 Technical Paper Effect of an hourglass shaped sleeve on the performance of the fluid dynamic bearings of a HDD spindle motor Jihoon Lee Minho Lee Gunhee Jang Received: 14
More informationNatural frequency analysis of fluid-conveying pipes in the ADINA system
Journal of Physics: Conference Series OPEN ACCESS Natural frequency analysis of fluid-conveying pipes in the ADINA system To cite this article: L Wang et al 2013 J. Phys.: Conf. Ser. 448 012014 View the
More information1 Introduction. Minho Lee 1 Jihoon Lee 1 Gunhee Jang 1
DOI 10.1007/s005-015-5-5 TECHNICAL PAPER Stability analysis of a whirling rigid rotor supported by stationary grooved FDBs considering the five degrees of freedom of a general rotor bearing system Minho
More informationRobust shaft design to compensate deformation in the hub press fitting and disk clamping process of 2.5 HDDs
DOI 10.1007/s00542-016-2850-2 TECHNICAL PAPER Robust shaft design to compensate deformation in the hub press fitting and disk clamping process of 2.5 HDDs Bumcho Kim 1,2 Minho Lee 3 Gunhee Jang 3 Received:
More informationResearch Article Stability Analysis of Journal Bearing: Dynamic Characteristics
Research Journal of Applied Sciences, Engineering and Technology 9(1): 47-52, 2015 DOI:10.19026/rjaset.9.1375 ISSN: 2040-7459; e-issn: 2040-7467 2015 Maxwell Scientific Publication Corp. Submitted: July
More informationTheory and Practice of Rotor Dynamics Prof. Rajiv Tiwari Department of Mechanical Engineering Indian Institute of Technology Guwahati
Theory and Practice of Rotor Dynamics Prof. Rajiv Tiwari Department of Mechanical Engineering Indian Institute of Technology Guwahati Module - 7 Instability in rotor systems Lecture - 4 Steam Whirl and
More informationStability Analysis of a Hydrodynamic Journal Bearing With Rotating Herringbone Grooves
G. H. Jang e-mail: ghjang@hanyang.ac.kr J. W. Yoon PREM, Department of Mechanical Engineering, Hanyang University, Seoul, 33-79, Korea Stability Analysis of a Hydrodynamic Journal Bearing With Rotating
More informationAnalysis of Fluid Film Stiffness and Damping coefficient for A Circular Journal Bearing with Micropolar Fluid
et International Journal on Emerging Technologies 5(1): 206-211(2014) ISSN No. (Print) : 0975-8364 ISSN No. (Online) : 2249-3255 Analysis of Fluid Film Stiffness Damping coefficient for A Circular Journal
More informationAEROELASTIC ANALYSIS OF COMBINED CONICAL - CYLINDRICAL SHELLS
Proceedings of the 7th International Conference on Mechanics and Materials in Design Albufeira/Portugal 11-15 June 2017. Editors J.F. Silva Gomes and S.A. Meguid. Publ. INEGI/FEUP (2017) PAPER REF: 6642
More informationInfluence of friction coefficient on rubbing behavior of oil bearing rotor system
Influence of friction coefficient on rubbing behavior of oil bearing rotor system Changliang Tang 1, Jinfu ang 2, Dongjiang Han 3, Huan Lei 4, Long Hao 5, Tianyu Zhang 6 1, 2, 3, 4, 5 Institute of Engineering
More informationDynamics of Rotor Systems with Clearance and Weak Pedestals in Full Contact
Paper ID No: 23 Dynamics of Rotor Systems with Clearance and Weak Pedestals in Full Contact Dr. Magnus Karlberg 1, Dr. Martin Karlsson 2, Prof. Lennart Karlsson 3 and Ass. Prof. Mats Näsström 4 1 Department
More informationAEROELASTIC ANALYSIS OF SPHERICAL SHELLS
11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver
More informationSAMCEF For ROTORS. Chapter 1 : Physical Aspects of rotor dynamics. This document is the property of SAMTECH S.A. MEF A, Page 1
SAMCEF For ROTORS Chapter 1 : Physical Aspects of rotor dynamics This document is the property of SAMTECH S.A. MEF 101-01-A, Page 1 Table of Contents rotor dynamics Introduction Rotating parts Gyroscopic
More informationSimulation and Experimental Research on Dynamics of Low-Pressure Rotor System in Turbofan Engine
Simulation and Experimental Research on Dynamics of Low-Pressure Rotor System in Turbofan Engine Shengxiang Li 1, Chengxue Jin 2, Guang Zhao 1*, Zhiliang Xiong 1, Baopeng Xu 1 1. Collaborative Innovation
More information1415. Effects of different disc locations on oil-film instability in a rotor system
1415. Effects of different disc locations on oil-film instability in a rotor system Hui Ma 1, Xueling Wang 2, Heqiang Niu 3, Hui Li 4 School of Mechanical Engineering and Automation, Northeastern University,
More informationDynamic analysis of a HDD spindle system with FDBs due to the bearing width and asymmetric grooves of journal bearing
Microsystem Technologies Micro- and Nanosystems Information Storage and Pro Springer-Verlag 2005 10.1007/s00542-005-0606-5 Technical paper Dynamic analysis of a HDD spindle system with FDBs due to the
More informationCHAPTER 1 INTRODUCTION Hydrodynamic journal bearings are considered to be a vital component of all the rotating machinery. These are used to support
CHAPTER 1 INTRODUCTION Hydrodynamic journal bearings are considered to be a vital component of all the rotating machinery. These are used to support radial loads under high speed operating conditions.
More informationHakwoon Kim Gunhee Jang Sanghoon Lee. 1 Introduction
Microsyst Technol (2011) 17:749 759 DOI 10.1007/s00542-010-1188-4 TECHNICAL PAPER Complete determination of the dynamic coefficients of coupled journal and thrust bearings considering five degrees of freedom
More informationDETC EXPERIMENT OF OIL-FILM WHIRL IN ROTOR SYSTEM AND WAVELET FRACTAL ANALYSES
Proceedings of IDETC/CIE 2005 ASME 2005 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference September 24-28, 2005, Long Beach, California, USA DETC2005-85218
More informationVIBRATION ANALYSIS OF TIE-ROD/TIE-BOLT ROTORS USING FEM
VIBRATION ANALYSIS OF TIE-ROD/TIE-BOLT ROTORS USING FEM J. E. Jam, F. Meisami Composite Materials and Technology Center Tehran, IRAN jejaam@gmail.com N. G. Nia Iran Polymer & Petrochemical Institute, Tehran,
More informationNON-LINEAR ROTORDYNAMICS: COMPUTATIONAL STRATEGIES
The 9th International Symposium on Transport Phenomena and Dynamics of Rotating Machinery Honolulu, Hawaii, February 1-14, NON-LINEAR ROTORDNAMICS: COMPUTATIONAL STRATEGIES Tom J. Chalko Head of Rotordynamic
More informationStability analysis of a whirling disk-spindle system supported by FDBs with rotating grooves
Microsyst Technol (011) 17:787 797 DOI 10.1007/s0054-010-111-9 TECHNICAL PAPER Stability analysis of a whirling disk-spindle system supported by FDBs with rotating grooves Jihoon Lee Gunhee Jang Kyungmoon
More informationSelf-Excited Vibration
Wenjing Ding Self-Excited Vibration Theory, Paradigms, and Research Methods With 228 figures Ö Springer Contents Chapter 1 Introduction 1 1.1 Main Features of Self-Excited Vibration 1 1.1.1 Natural Vibration
More informationNumerical analysis of three-lobe journal bearing with CFD and FSI
Numerical analysis of three-lobe journal bearing with CFD and FSI Pankaj Khachane 1, Dinesh Dhande 2 1PG Student at Department of Mechanical Engineering, AISSMSCOE Pune, Maharashtra, India 2Assistant Professor
More informationFinite element analysis of rotating structures
Finite element analysis of rotating structures Dr. Louis Komzsik Chief Numerical Analyst Siemens PLM Software Why do rotor dynamics with FEM? Very complex structures with millions of degrees of freedom
More informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More informationCircular Bearing Performance Parameters with Isothermal and Thermo-Hydrodynamic Approach Using Computational Fluid Dynamics
Circular Bearing Performance Parameters with Isothermal and Thermo-Hydrodynamic Approach Using Computational Fluid Dynamics Amit Chauhan 1 Department of Mechanical Engineering, University Institute of
More informationAnalysis on propulsion shafting coupled torsional-longitudinal vibration under different applied loads
Analysis on propulsion shafting coupled torsional-longitudinal vibration under different applied loads Qianwen HUANG 1 ; Jia LIU 1 ; Cong ZHANG 1,2 ; inping YAN 1,2 1 Reliability Engineering Institute,
More informationImplementation of a Thermo- Hydrodynamic Model to Predict Morton Effect
Implementation of a Thermo- Hydrodynamic Model to Predict Morton Effect Antonini *, Fausti and Mor Polibrixia srl, Via A. Tadini 49, 25125 Brescia. *orresponding author: Via Branze 45, 25123 Brescia, massimo.antonini@polibrixia.it
More information1859. Forced transverse vibration analysis of a Rayleigh double-beam system with a Pasternak middle layer subjected to compressive axial load
1859. Forced transverse vibration analysis of a Rayleigh double-beam system with a Pasternak middle layer subjected to compressive axial load Nader Mohammadi 1, Mehrdad Nasirshoaibi 2 Department of Mechanical
More informationSTICK-SLIP WHIRL INTERACTION IN DRILLSTRING DYNAMICS
STICK-SLIP WHIRL INTERACTION IN DRILLSTRING DYNAMICS R. I. Leine, D. H. van Campen Department of Mechanical Engineering, Eindhoven University of Technology, P. O. Box 513, 5600 MB Eindhoven, The Netherlands
More informationAvailable online at ScienceDirect. Procedia Technology 23 (2016 ) 42 50
Available online at www.sciencedirect.com ScienceDirect Procedia Technology 23 (216 ) 42 5 3rd International Conference on Innovations in Automation and Mechatronics Engineering, ICIAME 216 On the stiffness
More informationLinear and Nonlinear Analysis of Plain Journal Bearings Lubricated With Couple Stress Fluid
ISSN 2395-1621 Linear and Nonlinear Analysis of Plain Journal Bearings Lubricated With Couple Stress Fluid #1 Deepali Kangude 1 deepalikangude94@gmail.com 1 P.G. student Mechanical Department, DYPIET Pimpri,
More informationDEVELOPMENT OF TOOL FOR THE IDENTIFICATION OF STIFFNESS AND DAMPING COEFFICIENTS OF JOURNAL BEARING
DEVELOPMENT OF TOOL FOR THE IDENTIFICATION OF STIFFNESS AND DAMPING COEFFICIENTS OF JOURNAL BEARING Sampath Kumar B 1 and G.Satish Babu 2 1 P.G student, Department of mechanical engineering, JNTUH College
More informationSimulation of fluid-structure interaction for an elastic cylinder in an axial flow
Advances in Fluid Mechanics IX 151 Simulation of fluid-structure interaction for an elastic cylinder in an axial flow Z. G. Liu 1, Y. Liu 1 & J. Lu 2 1 Department of Mechanical Engineering, Hong Kong Polytechnic
More informationContents. I Introduction 1. Preface. xiii
Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................
More informationJournal of Advanced Mechanical Design, Systems, and Manufacturing
Static Analysis of Viscoelastic Supported Gas Foil Thrust Bearing with Journal Inclination* Quan ZHOU**, Yu HOU***, Rugang CHEN***, Shuangtao CHEN*** and Chunzheng CHEN*** ** School of Environmental and
More informationHydrodynamic Lubrication
ME 383S Bryant February 15, 2005 1 Hydrodynamic Lubrication Fluid Lubricant: liquid or gas (gas bearing) Mechanism: Pressures separate surfaces o Normal loads on bodies o Convergent profile between surfaces
More information3 Mathematical modeling of the torsional dynamics of a drill string
3 Mathematical modeling of the torsional dynamics of a drill string 3.1 Introduction Many works about torsional vibrations on drilling systems [1, 12, 18, 24, 41] have been published using different numerical
More informationDynamics of assembled structures of rotor systems of aviation gas turbine engines of type two-rotor
Dynamics of assembled structures of rotor systems of aviation gas turbine engines of type two-rotor Anatoly А. Pykhalov 1, Mikhail А. Dudaev 2, Mikhail Ye. Kolotnikov 3, Paul V. Makarov 4 1 Irkutsk State
More informationAnalysis of flow characteristics of a cam rotor pump
IOP Conference Series: Materials Science and Engineering OPEN ACCESS Analysis of flow characteristics of a cam rotor pump To cite this article: Y Y Liu et al 2013 IOP Conf. Ser.: Mater. Sci. Eng. 52 032023
More informationDynamic behavior of turbine foundation considering full interaction among facility, structure and soil
Dynamic behavior of turbine foundation considering full interaction among facility, structure and soil Fang Ming Scholl of Civil Engineering, Harbin Institute of Technology, China Wang Tao Institute of
More informationApplication of CFD analysis for static and dynamic characteristics of hydrodynamic journal bearing
Application of CFD analysis for static and dynamic characteristics of hydrodynamic journal bearing Kenichiro Takenaka *, Tomoyuki Hayashi, Kazuyoshi Miyagawa Abstract ISROMAC 7 International Symposium
More informationDynamic Analysis of a Rotor-Journal Bearing System of Rotary Compressor
Purdue University Purdue e-pubs International Compressor Engineering Conference School of Mechanical Engineering 6 Dynamic Analysis of a Rotor-Journal Bearing System of Rotary Compressor Fei Xie Xi an
More informationStudy of the influence of the resonance changer on the longitudinal vibration of marine propulsion shafting system
Study of the influence of the resonance changer on the longitudinal vibration of marine propulsion shafting system Zhengmin Li 1, Lin He 2, Hanguo Cui 3, Jiangyang He 4, Wei Xu 5 1, 2, 4, 5 Institute of
More informationPROJECT 2 DYNAMICS OF MACHINES 41514
PROJECT 2 DYNAMICS OF MACHINES 41514 Dynamics of Rotor-Bearing System Lateral Vibrations and Stability Threshold of Rotors Supported On Hydrodynamic Bearing and Ball Bearing. Ilmar Ferreira Santos, Prof.
More informationDynamic Model of a Badminton Stroke
ISEA 28 CONFERENCE Dynamic Model of a Badminton Stroke M. Kwan* and J. Rasmussen Department of Mechanical Engineering, Aalborg University, 922 Aalborg East, Denmark Phone: +45 994 9317 / Fax: +45 9815
More informationAeroelastic effects of large blade deflections for wind turbines
Aeroelastic effects of large blade deflections for wind turbines Torben J. Larsen Anders M. Hansen Risoe, National Laboratory Risoe, National Laboratory P.O. Box 49, 4 Roskilde, Denmark P.O. Box 49, 4
More informationModel of the gas journal bearing dynamics with a flexibly supported foil
Model of the gas journal bearing dynamics with a flexibly supported foil (NUM209-15) Eliza Tkacz, Zbigniew Kozanecki, Dorota Kozanecka, Jakub Lagodziski Abstract: The work is devoted to an analysis of
More informationInvestigation into the Self-loosening Trend of Bolt Joints on the Tower Crane ShengChun Wang1,a, PeiWei Ni1,b, Ye Zhang1,c and MingXiao Dong1,2,d
2nd International Conference on Advances in Mechanical Engineering and Industrial Informatics (AMEII 2016) Investigation into the Self-loosening Trend of Bolt Joints on the Tower Crane ShengChun Wang1,a,
More informationGeometrically exact beam dynamics, with and without rotational degree of freedom
ICCM2014 28-30 th July, Cambridge, England Geometrically exact beam dynamics, with and without rotational degree of freedom *Tien Long Nguyen¹, Carlo Sansour 2, and Mohammed Hjiaj 1 1 Department of Civil
More informationEffects of Structural Forces on the Dynamic Performance of High Speed Rotating Impellers.
Effects of Structural Forces on the Dynamic Performance of High Speed Rotating Impellers. G Shenoy 1, B S Shenoy 1 and Raj C Thiagarajan 2 * 1 Dept. of Mechanical & Mfg. Engineering, Manipal Institute
More informationAnswers to questions in each section should be tied together and handed in separately.
EGT0 ENGINEERING TRIPOS PART IA Wednesday 4 June 014 9 to 1 Paper 1 MECHANICAL ENGINEERING Answer all questions. The approximate number of marks allocated to each part of a question is indicated in the
More informationModeling and Performance Analysis of a Flywheel Energy Storage System Prince Owusu-Ansah, 1, Hu Yefa, 1, Philip Agyeman, 1 Adam Misbawu 2
International Conference on Electromechanical Control Technology and Transportation (ICECTT 2015) Modeling and Performance Analysis of a Flywheel Energy Storage System Prince Owusu-Ansah, 1, Hu Yefa, 1,
More informationA novel fluid-structure interaction model for lubricating gaps of piston machines
Fluid Structure Interaction V 13 A novel fluid-structure interaction model for lubricating gaps of piston machines M. Pelosi & M. Ivantysynova Department of Agricultural and Biological Engineering and
More informationNew Representation of Bearings in LS-DYNA
13 th International LS-DYNA Users Conference Session: Aerospace New Representation of Bearings in LS-DYNA Kelly S. Carney Samuel A. Howard NASA Glenn Research Center, Cleveland, OH 44135 Brad A. Miller
More informationSYNTHESIS OF A FLUID JOURNAL BEARING USING A GENETIC ALGORITHM
SYNTHESIS OF A FLUID JOURNAL BEARING USING A GENETIC ALGORITHM A. MANFREDINI and P. VIGNI Dipartimento di Ingegneria Meccanica, Nucleare e della Produzione (DIMNP) - University of Pisa Via Diotisalvi,
More informationStudy on Nonlinear Dynamic Response of an Unbalanced Rotor Supported on Ball Bearing
G. Chen College of Civil Aviation, Nanjing University of Aeronautics and Astronautics, Nanjing, 210016, P.R.C. e-mail: cgzyx@263.net Study on Nonlinear Dynamic Response of an Unbalanced Rotor Supported
More informationThe Ultimate Load-Carrying Capacity of a Thin-Walled Shuttle Cylinder Structure with Cracks under Eccentric Compressive Force
The Ultimate Load-Carrying Capacity of a Thin-Walled Shuttle Cylinder Structure with Cracks under Eccentric Compressive Force Cai-qin Cao *, Kan Liu, Jun-zhe Dong School of Science, Xi an University of
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationVibrations in Mechanical Systems
Maurice Roseau Vibrations in Mechanical Systems Analytical Methods and Applications With 112 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Contents Chapter I. Forced Vibrations
More informationA consistent dynamic finite element formulation for a pipe using Euler parameters
111 A consistent dynamic finite element formulation for a pipe using Euler parameters Ara Arabyan and Yaqun Jiang Department of Aerospace and Mechanical Engineering, University of Arizona, Tucson, AZ 85721,
More information892 VIBROENGINEERING. JOURNAL OF VIBROENGINEERING. JUNE VOLUME 15, ISSUE 2. ISSN
1004. Study on dynamical characteristics of misalignrubbing coupling fault dual-disk rotor-bearing system Yang Liu, Xing-Yu Tai, Qian Zhao, Bang-Chun Wen 1004. STUDY ON DYNAMICAL CHARACTERISTICS OF MISALIGN-RUBBING
More informationJitter and Basic Requirements of the Reaction Wheel Assembly in the Attitude Control System
Jitter and Basic Requirements of the Reaction Wheel Assembly in the Attitude Control System Lulu Liu August, 7 1 Brief Introduction Photometric precision is a major concern in this space mission. A pointing
More informationStability of Water-Lubricated, Hydrostatic, Conical Bearings With Spiral Grooves for High-Speed Spindles
S. Yoshimoto Professor Science University of Tokyo, Department of Mechanical Engineering, 1-3 Kagurazaka Shinjuku-ku, Tokyo 16-8601 Japan S. Oshima Graduate Student Science University of Tokyo, Department
More informationStructural Dynamics Lecture 2. Outline of Lecture 2. Single-Degree-of-Freedom Systems (cont.)
Outline of Single-Degree-of-Freedom Systems (cont.) Linear Viscous Damped Eigenvibrations. Logarithmic decrement. Response to Harmonic and Periodic Loads. 1 Single-Degreee-of-Freedom Systems (cont.). Linear
More informationVibration Analysis of Hollow Profiled Shafts
International Journal of Current Engineering and echnology ISSN 77-406 04 INPRESSCO. All Rights Reserved. Available at http://inpressco.com/category/ijcet Research Article Vibration Analysis of Hollow
More informationRigid Flexible Contact Analysis of an Inflated Membrane Balloon with Various Contact Conditions
roceedings Rigid Flexible Contact Analysis of an Inflated Membrane Balloon with Various Contact Conditions MengXiong Liu and XiDe Li * Department of Engineering Mechanics, Tghua University, Beijing 100084,
More informationInvestigations of Slip Effect on the Performance of Micro Gas Bearings and Stability of Micro Rotor-Bearing Systems
Sensors 007, 7, 99-44 sensors ISSN 44-80 007 by MPI www.mdpi.org/sensors Full Research Paper Investigations of Slip Effect on the Performance of Micro Gas Bearings and Stability of Micro Rotor-Bearing
More informationShape Optimization of Oldham Coupling in Scroll Compressor
Purdue University Purdue e-pubs International Compressor Engineering Conference School of Mechanical Engineering 24 Shape Optimization of Oldham Coupling in Scroll Compressor In Hwe Koo LG Electronics
More informationDYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION
October 1-17,, Beijing, China DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION Mohammad M. Ahmadi 1 and Mahdi Ehsani 1 Assistant Professor, Dept. of Civil Engineering, Geotechnical Group,
More informationD : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown.
D : SOLID MECHANICS Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown. Q.2 Consider the forces of magnitude F acting on the sides of the regular hexagon having
More informationThermohydrodynamic analysis of a worn plain journal bearing
Tribology International 37 (2004) 129 136 www.elsevier.com/locate/triboint Thermohydrodynamic analysis of a worn plain journal bearing M. Fillon, J. Bouyer Université de Poitiers, Laboratoire de Mécanique
More information1544. Synchronous and subsynchronous vibration under the combined effect of bearings and seals: numerical simulation and its experimental validation
1544. Synchronous and subsynchronous vibration under the combined effect of bearings and seals: numerical simulation and its experimental validation Wanfu Zhang 1, Jiangang Yang 2, Chun Li 3, Ren Dai 4,
More informationA multiscale framework for lubrication analysis of bearings with textured surface
A multiscale framework for lubrication analysis of bearings with textured surface *Leiming Gao 1), Gregory de Boer 2) and Rob Hewson 3) 1), 3) Aeronautics Department, Imperial College London, London, SW7
More informationD : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each.
GTE 2016 Q. 1 Q. 9 carry one mark each. D : SOLID MECHNICS Q.1 single degree of freedom vibrating system has mass of 5 kg, stiffness of 500 N/m and damping coefficient of 100 N-s/m. To make the system
More informationStudy of coupling between bending and torsional vibration of cracked rotor system supported by radial active magnetic bearings
Applied and Computational Mechanics 1 (2007) 427-436 Study of coupling between bending and torsional vibration of cracked rotor system supported by radial active magnetic bearings P. Ferfecki a, * a Center
More informationPerturbation Method in the Analysis of Manipulator Inertial Vibrations
Mechanics and Mechanical Engineering Vol. 15, No. 2 (2011) 149 160 c Technical University of Lodz Perturbation Method in the Analysis of Manipulator Inertial Vibrations Przemys law Szumiński Division of
More informationDynamic behavior of thrust air bearings
Department of Mechanical Engineering Dynamics and Control Research Group Dynamic behavior of thrust air bearings R.H.M. Franssen D&C 6. Master s thesis Coaches: dr. ir. R.H.B. Fey Commitee: prof. dr. H.
More informationShape Optimization of Revolute Single Link Flexible Robotic Manipulator for Vibration Suppression
15 th National Conference on Machines and Mechanisms NaCoMM011-157 Shape Optimization of Revolute Single Link Flexible Robotic Manipulator for Vibration Suppression Sachindra Mahto Abstract In this work,
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationProgram System for Machine Dynamics. Abstract. Version 5.0 November 2017
Program System for Machine Dynamics Abstract Version 5.0 November 2017 Ingenieur-Büro Klement Lerchenweg 2 D 65428 Rüsselsheim Phone +49/6142/55951 hd.klement@t-online.de What is MADYN? The program system
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationα Cubic nonlinearity coefficient. ISSN: x DOI: : /JOEMS
Journal of the Egyptian Mathematical Society Volume (6) - Issue (1) - 018 ISSN: 1110-65x DOI: : 10.1608/JOEMS.018.9468 ENHANCING PD-CONTROLLER EFFICIENCY VIA TIME- DELAYS TO SUPPRESS NONLINEAR SYSTEM OSCILLATIONS
More informationExperimental Investigations of Whirl Speeds of a Rotor on Hydrodynamic Spiral Journal Bearings Under Flooded Lubrication
International Conference on Fluid Dynamics and Thermodynamics Technologies (FDTT ) IPCSIT vol.33() () IACSIT Press, Singapore Experimental Investigations of Whirl Speeds of a Rotor on Hydrodynamic Spiral
More informationOn The Finite Element Modeling Of Turbo Machinery Rotors In Rotor Dynamic Analysis
Proceedings of The Canadian Society for Mechanical Engineering International Congress 2018 CSME International Congress 2018 May 27-30, 2018, Toronto, On, Canada On The Finite Element Modeling Of Turbo
More informationAttempt ALL QUESTIONS IN SECTION A, ONE QUESTION FROM SECTION B and ONE QUESTION FROM SECTION C Linear graph paper will be provided.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 2015-2016 ENERGY ENGINEERING PRINCIPLES ENG-5001Y Time allowed: 3 Hours Attempt ALL QUESTIONS IN SECTION A, ONE QUESTION FROM
More informationCO-ROTATIONAL DYNAMIC FORMULATION FOR 2D BEAMS
COMPDYN 011 ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering M. Papadrakakis, M. Fragiadakis, V. Plevris (eds.) Corfu, Greece, 5-8 May 011 CO-ROTATIONAL
More informationImplementation of an advanced beam model in BHawC
Journal of Physics: Conference Series PAPER OPEN ACCESS Implementation of an advanced beam model in BHawC To cite this article: P J Couturier and P F Skjoldan 28 J. Phys.: Conf. Ser. 37 625 Related content
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More information1631. Dynamic analysis of offset press gear-cylinder-bearing system applying finite element method
1631. Dynamic analysis of offset press gear-cylinder-bearing system applying finite element method Tiancheng OuYang 1, Nan Chen 2, Jinxiang Wang 3, Hui Jing 4, Xiaofei Chen 5 School of Mechanical Engineering,
More informationFlow-Induced Vibration of Pipeline on Elastic Support
Available online at www.sciencedirect.com Procedia Engineering () 986 99 The Twelfth East Asia-Pacific Conference on Structural Engineering and Construction Flow-Induced Vibration of Pipeline on Elastic
More informationPRESSURE AND VELOCITY AMPLITUDES OF THE INCOMPRESSIBLE FLUID IN CONCENTRIC ANNULAR PASSAGE WITH OSCILLATORY BOUNDARY: TURBULENT FLOW
Journal of Engineering Science and Technology Vol. 9, No. 2 (2014) 220-232 School of Engineering, Taylor s University PRESSURE AND VELOCITY AMPLITUDES OF THE INCOMPRESSIBLE FLUID IN CONCENTRIC ANNULAR
More informationDynamic Response Of Laminated Composite Shells Subjected To Impulsive Loads
IOSR Journal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-1684,p-ISSN: 2320-334X, Volume 14, Issue 3 Ver. I (May. - June. 2017), PP 108-123 www.iosrjournals.org Dynamic Response Of Laminated
More informationWORK SHEET FOR MEP311
EXPERIMENT II-1A STUDY OF PRESSURE DISTRIBUTIONS IN LUBRICATING OIL FILMS USING MICHELL TILTING PAD APPARATUS OBJECTIVE To study generation of pressure profile along and across the thick fluid film (converging,
More informationLubrication and Journal Bearings
UNIVERSITY OF HAIL College of Engineering Department of Mechanical Engineering Chapter 12 Lubrication and Journal Bearings Text Book : Mechanical Engineering Design, 9th Edition Dr. Badreddine AYADI 2016
More informationTheory and Practice of Rotor Dynamics Prof. Dr. Rajiv Tiwari Department of Mechanical Engineering Indian Institute of Technology Guwahati
Theory and Practice of Rotor Dynamics Prof. Dr. Rajiv Tiwari Department of Mechanical Engineering Indian Institute of Technology Guwahati Module - 2 Simpul Rotors Lecture - 2 Jeffcott Rotor Model In the
More informationAnalytical Strip Method for Thin Isotropic Cylindrical Shells
IOSR Journal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-1684,p-ISSN: 2320-334X, Volume 14, Issue 4 Ver. III (Jul. Aug. 2017), PP 24-38 www.iosrjournals.org Analytical Strip Method for
More information