SV Shock and Vibration 1875-9203 1070-9622 Hindawi 10.1155/2019/9154272 9154272 Research Article Dynamic Analysis of Flexible Rotor Based on Transfer Symplectic Matrix Deng Hao 1 https://orcid.org/0000-0003-3849-5427 Fang Xi 1 https://orcid.org/0000-0001-5190-5465 Wu Huachun 2 Ding Yiming 1 Yu Jinghu 1 Zhang Xiaoyu 1 Li Hong 3 Liu Chenwei 4 Marchesiello Stefano 1 School of Science Wuhan University of Technology Wuhan 430070 China whut.edu.cn 2 School of Mechanical and Electronic Engineering Wuhan University of Technology Wuhan 430070 China whut.edu.cn 3 Ordnance NCO Academy Army Engineering University of PLA Wuhan 430075 China ust.com.cn 4 Midea Group Wuhan Refrigeration Equipment Co. Ltd. Wuhan Hubei 43005 China 2019 992019 2019 05 05 2019 16 07 2019 11 08 2019 992019 2019 Copyright © 2019 Hao Deng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

In view of the numerical instability and low accuracy of the traditional transfer matrix method in solving the high-order critical speed of the rotor system, a new idea of incorporating the finite element method into the transfer matrix is proposed. Based on the variational principle, the transfer symplectic matrix of gyro rotors suitable for all kinds of boundary conditions and supporting conditions under the Hamilton system is derived by introducing dual variables. To verify the proposed method in rotor critical speed, a numerical analysis is adopted. The simulation experiment results show that, in the calculation of high-order critical speed, especially when exceeding the sixth critical speed, the numerical accuracy of the transfer symplectic matrix method is obviously better than that of the reference method. The relative errors between the numerical solution and the exact solution are 0.0347% and 0.2228%, respectively, at the sixth critical speed. The numerical example indicates the feasibility and superiority of the method, which provides the basis for the optimal design of the rotor system.

National Natural Science Foundation of China 61573012 61702388 Fundamental Research Funds for Central University WUT2016A004
1. Introduction

With the development of modern industrial technology, the application prospect and running speed of rotary machinery have been greatly improved. Correspondingly, the performance and design requirements of rotary machinery have also been put forward with greater challenges . In order to make the rotor work smoothly, avoid resonance and other vibration faults, and ensure high rotation accuracy, the study of the natural frequency and mode of the rotor dynamics system, the critical speed, and the unbalanced response has become an indispensable important content of rotor dynamics. Critical speed and vibration mode, which are important dynamic characteristic parameters of rotating machinery, are closely related to various structural parameters of the rotor system itself, such as the influence of unbalance and support stiffness on critical speed. Therefore, it is of great significance to analyze the sensitivity of the dynamic characteristic parameters of the rotor system to the structural parameters of the system and to optimize the rotor structure and shorten the design period .

So far, the numerical analysis methods of rotor dynamics are mainly divided into two categories: transfer matrix method and finite element method. The finite element method combines the variational principle of mathematics with the energy principle of mechanics [3, 4]. According to the principle of minimum potential energy of elasticity, when a deformed elastomer is in an equilibrium state under the action of external forces, the curve that minimizes the total potential energy is the real deformation curve among many possible deformation permissible curves. The expression of the rotor finite element method is concise and standard, which has prominent advantages when solving the problems of the large complex system composed of rotor and surrounding structures, but this method will increase the system freedom degree and result in a lot of time and space spending .

For the multisupported multidisk rotor system, when the system has more degrees of freedom, whether by establishing the motion differential equation to solve the critical speed and the unbalanced response of the system or by solving the eigenvalue problem, the calculation workload is very large. Therefore, a new computational method, the classical transfer matrix method (Prohl transfer matrix method), was developed. This method is widely used because the order of the matrix will not increase with the increase of the degrees of freedom of the system and less storage units are needed . Nevertheless, with the increase of frequency of trial computation, the accuracy decreases, especially in solving the critical speeds and unbalance response of the rotor structure system, which is of large scale and high speed. The vibration mode of the structural system firstly appears numerical instability, which is mainly manifested as a sharp increase in the amplitude at the end of the shaft system, and then the root of the frequency equation begins to be errors and omissions, resulting in a larger error in the calculation results. The Riccati transfer matrix method is then proposed to transform the original two-point boundary value problem of the differential equation into a one-point initial value problem through the Riccati transformation, which improves the numerical stability and calculation accuracy while retaining the advantages of the classical transfer matrix method . However, when solving the critical speed, there may be some problems such as increasing roots or missing roots because the interval of residual quantity with different signs does not correspond to the interval of residual quantity root.

The transfer matrix method mentioned above is to analyze the dynamic characteristics of the rotor under the Lagrange system. Academician Zhong Wanxie introduced some important issues of classical mechanics into the Hamilton system; he studied the dynamics of the gyro rotor from a new perspective and proved the effectiveness of analyzing the gyro rotor system under the Hamilton system .

Consider the advantages and disadvantages of each method mentioned above, in order to make the rotor work smoothly, avoid resonance, and guarantee its high rotation accuracy; the idea of finite element is incorporated into the transfer matrix in this paper. Starting from the variational principle, dual variables are introduced to deduce the transfer matrix of the gyro rotor under the Hamilton system. Then, numerical examples are given to verify the feasibility and superiority of this method. Since the transfer matrix between state vectors is derived by symplectic method, it has significant advantages in long-term numerical stability. In addition, compared with the traditional transfer matrix method, this method avoids a series of complex operations such as solving the aggregation model of the beam, which is easy to understand, and can obtain relatively accurate results when calculating the high-order vortex frequency.

A brief outline of this paper is as follows. In the second section, we firstly review the classical transfer matrix methods such as Prohl transfer matrix method and Riccati transfer matrix method and then introduce the dual-rotor system and its transfer matrix method. In the third section, we derive the transfer symplectic matrix between state vectors of the rotor system under the Hamiltonian system from the variational principle and briefly explained that this type of transfer matrix is symplectic. In the fourth section, the numerical example is used for verification and different numerical algorithms are compared. In the fifth section, we mainly introduce the conclusions and innovations.

2. Theoretical Basis 2.1. Classical Transfer Matrix Method

As one of the analysis methods of rotor dynamics, the transfer matrix method is usually used to solve rotor dynamics problems in the frequency domain. As the transfer matrix method does not require the storage and the dimension does not increase with the number of elements (for it is replaced with more matrix multiplications), this kind of method is especially suitable for chain systems such as rotor systems.

Figure 1 depicts the rotor structure adopted in this section. Before analyzing the dynamics of the subsystem by the transfer matrix method, the rotor subsystem is essential to be discretized and simplified into the lumped mass model  (shown in Figure 2). The assembled parts on the rotor, such as the front and rear lamination stacks, are modeled as rigid disks, with consideration of the gyroscopic moment. Considering the xOz plane of the discrete rotor model, a typical unit consists of a massless shaft segment and a lumped mass disc .

Rotor model and structure drawing.

Shaft torsional vibration system.

Considering the ith element, defining Z as the stator vector containing deflection y, with θ as the slope, M as the bending moment, and Q as the shear force, we obtain(1)Zi=y,θ,M,QiT.

For any axial segment i, there is a certain relationship between the state vectors Zi and Zi+1, that is(2)Zi+1=TiZi,where Ti is called the transfer matrix of components between the two sections and the matrix elements can be obtained by force analysis and deformation relation of components.

In order to facilitate the study of the relationship between adjacent nodes and axis segments, this paper did the simplification as follows: The bending moment and shear force to the right side of node i1 represent the bending moment and shear force to the left side of axis segment i, respectively, while the bending moment and shear force to the left of node i represent the bending moment and shear force to the right side of axis segment i, respectively.

According to the bending deformation formula and displacement analysis of beam, the following equation is obtained :(3)Zi+1L=AiZiR,where the field transfer matrix Ai is(4)Ai=1lili22EiIili31υ6EiIi01liEiIili22EiIi001li0001,in which υ=ks6EiIi/GiAili2 represents the shear effect coefficient, ks is called the cross-sectional shape coefficient, and ksli/GiAi is the influence of shear deformation on bending deformation.

For isotropic rotor systems, when the system moves in a simple harmonic motion with the angular velocity ω, we get y¨i=ω2yi. Suppose that the jth node has a spring bearing and the spring stiffness is kj. According to rigid body dynamics and D’Alembert’s principle, the relationship between the state vectors at the left and right ends of the rotating disk supported by elasticity can be obtained as (5)ZiR=CiZiL,where the transfer matrix Ci is(6)Ci=100001000JdiJpiΩωω210miω2kj001,in which Jdi is the moment of inertia of diameter, Jpi is the polar moment of inertia, Ω is the angular velocity of axis rotation, and ω is the angular speed of revolution.

By combining the flexible rotor and rigid disk with a transfer element, the relation of the state vectors from ith to i+1th can be written as(7)Zi+1L=AiZiR=AiCiZiLTiZiL.

When the transfer matrix of each part of the rotor system is obtained from the above calculation, the whole system is constructed by successively multiplying the matrix of each part according to equation (7):(8)ZN+1=TNZN==TNTN1T2T1Z1=DZ1,where the assembled transferring matrix D is(9)D=d11d12d13d14d21d22d23d24d31d32d33d34d41d42d43d44.

It is found that the value of the element containing ω in the total transfer matrix D is only restricted by the boundary conditions on both sides, so the critical angular velocity of the system can be obtained by solving the ω, which satisfies the boundary conditions .

When the traditional transfer matrix method is used to calculate the response of high-frequency dynamics or systems with an excessive number of nodes, the numerical instability occurs. In order to solve the problem, the Riccati transfer matrix method appears. The Riccati transfer matrix method occurs which converts a numerically unstable two-point boundary value problem into a numerically stable value problem. Through the Riccati transform, the state vector Zi+1 is converted into(10)Zi+1=fei+1=u11u12u21u22fei,where f and e represent the force part of the state vector and the displacement part of the state vector, respectively. A generalized Riccati transformation at the ith section is given as follows:(11)fi=siei,where si is the so-called Riccati transfer matrix. And through the derivation, we get the recursion formula(12)si+1=u11s+u12iu21s+u22i1.

2.2. Dual-Rotor System Transfer Matrix Method

Consider a dual-rotor system with two intermediate support, as shown in Figure 3. Defining the state vectors of the ith element of the two rotors as Zi1 and Zi2, respectively:(13)Zi1=y1,θ1,M1,Q1iT,Zi2=y2,θ2,M2,Q2iT.

Schematic diagram of dual-rotor system.

As for the noncoupling element, according to the transfer matrix method of single rotor system, their respective transfer relations can be obtained:(14)Z1Z2i+1=Ti100Ti2Z1Z2i,denoted by(15)Zi+1=TiZi,Ti1=1+l36EIC1C2l+l22EIC3l22EIl36EIC1l22EIC21+lEIC3lEIl22EIlC2C31lC2001i1,Ti2=1+l36EIC1C2l+l22EIC3l22EIl36EIC1l22EIC21+lEIC3lEIl22EIlC2C31lC2001i2,in which Ti1 and Ti2 are the transfer matrix of the two rotors. Their form can be derived from the Prohl, Riccati, or other transfer matrix methods. C1=1v, C2=mω2k, and C3=JdJpΩ/ωω2+K, where k and K are defined as 0 when there is no spring bearing in the ith element.

As for the coupling element, assume that the coupling point is in the ith unit. According to the balance of forces, the coordination equations are obtained:(16)y1θ1M1Q1iR=G1G2y1θ1M1Q1y2θ2M2Q2iLT,y2θ2M2Q2iR=G2G1y1θ1M1Q1y2θ2M2Q2iLT,so that we can get state vectors on both sides of the coupling element:(17)ZiR=TGiZiL,where TGi represents the transfer matrix for the coupling element, denoted as(18)TGi=G1G2G2G1i,G1=11Kh1k1i,G2=Khki,where k stands for support stiffness of the spring and Kh is the torque stiffness.

2.3. Critical Speed and Mode of Vibration

Supposing that the assembled transferring matrix is denoted by ZN+1R=DZ1R, according to the boundary conditions of the model(19)M1=0,Q1=0,MN+1=0,QN+1=0.

Combine equations (17) and (19), there is(20)d31d32d42d42yθ0=00.

Equation (20) is a system of homogeneous algebraic equations, where each element is a function of vortex frequency, and its coefficient determinant Δω2 called surplus is(21)Δω2=d31d32d42d42.

The condition for the existence of nonzero solutions to the homogeneous algebraic equations is that the determinant of coefficients is equal to zero, thus obtaining the rotor frequency equation:(22)Δω2=Δd31d32d42d42=0.

The value ω satisfying equation (22) is not only a synchronous forward vortex angular velocity of the rotor but also the critical angular velocity of the rotor.

3. Construction of the Transfer Matrix in Symplectic Space

Different from the traditional transfer matrix construction, this paper deduces the transfer symplectic matrix of the linear rotor system in symplectic space from the variational principle. The idea of finite element is incorporated into the transfer matrix, and the transfer symplectic matrix of gyro rotor under the Hamilton system is derived by introducing dual variables from the variational principle. Using this method to calculate the critical speed has the advantages of good numerical stability, high calculation accuracy and easy to compile general calculation program, and it has been verified by numerical examples.

Consider the following types of rotor systems: isotropic Timoshenko beams are used at the connection beams. Take a section of beam, which considers the moment of inertia and shear deformation for analysis. Let the z-coordinate be along the axis of the beam before deformation. Let μ be the cross-section shear coefficient, G be the shear modulus, E be the elastic modulus, I be the moment of inertia, A be the cross-section area of the beam, and ux,uy and ψx,ψy be the linear and angular displacements of the cross section, respectively. The deformation energy of this beam can be expressed as(23)V=120lEIdψxdz2+dψydz2+μAGduxdzψx2+duydzψy2dz.

Take the xOz perpendicular plane, for example, the internal force of the beam can be expressed as(24)M=EIdψxdz,Q=μAGduxdzψx.

The dynamic equations of the beam can be deduced to(25)Qz=ρA2uxt2,QMz=ρI2ψxt2.

For the sake of analysis, the time domain is converted to the frequency domain and the variable is still represented by the previous notation(26)u=uz,ωeiωt,ψ=ψz,ωeiωt,where ω is the rotation angular velocity of the rotor, and then, the dynamic equation (25) can be converted into(27)ddzμAGduxdzψx+ρAω2ux=0,ddzEIdψxdz+μAGduxdzψx+ρIω2ψx=0.

By the same logic, we can get that of the yOz plane. Thus, we can obtain the kinetic energy of this beam:(28)Ek=120lρAduxdz2+duydz2+ρIdψxdz2+dψydz2+2ΩρIdψxdzψydψydzψxdz.

Therefore, the Lagrangian density function is(29)L=12ρAduxdz2+duydz2+ρIdψxdz2+dψydz2+2ΩρIdψxdzψydψydzψxEIdψxdz2+dψydz2μAGduxdzψx2+duydzψy2.

For the convenience of analysis and calculation, the following notation is introduced: q=u,ψT;q˙=u˙,ψ˙T. Take the matrix form(30)L=12qTρAω2ρIω2qqTΩρIωq12q˙TμAGEIq˙q˙TμAGq12qTμAGq.

Then, the potential energy variation principle can be given for the above equation as follows:(31)δ120lLq,q˙dz=0,Lq,q˙=12q˙TK22q˙+q˙TK21q+12qTK11q.

The corresponding variational principle is Sqz=0LLq,q˙dz,δS=0 and of which(32)K22=μAGEI,K21=μAG,K11=ρAω2ρIω2μAG2ΩρIω,K12=K21T.

Then, the dynamic equation can be re-expressed as(33)K22q¨+K21K12q˙K11q=0.

Considering that both the variational expression and the dynamic equation only include the displacement type of variables, the method of introducing dual variables is adopted, let p=L/q˙ and then it can be re-expressed as(34)p=K22q˙+K21q.

Substitute equation (34) into the dynamic equation, we can get(35)q˙=Aq+Cp,p˙=BqATp,where A=K221K21,B=K11K12K221K21, and EC=K221.

Then, the Hamiltonian Hp,q is introduced as(36)Hp,q=pTAq12qTBq+12pTCp.

The corresponding variational principle is Sqz,pz=z0zfpTq˙Hq,pdz,δS=0, and the dual vectors q and p are unified to form the full state vector vz=qT,pTT.

The equation above is for the general n-dimensional displacement vector; in this paper, there is n=2, then p=K22q˙+K21q=μAGu˙ψ,EIψ˙T=Q.M. Moreover, at the same time, we can get the Hamiltonian equation of state as(37)v˙=Hv,H=ACBAT,in which(38)A=11,B=ρω2Aρω2I,C=μAG1EI1.

It is observed that the mechanical meaning of dual vector p is composed of generalized internal force, namely, shear force Q and bending moment M, which shows the rationality and validity of the introduced dual variables.

When solving the ordinary differential equation (37), we let li=Zi+1Zi and(39)vi+1=eHilivi.

Along the direction of the Z-axis, divide the Z-axis into arbitrarily small intervals. The transformation matrix when transferring from section i to section i+1 is(40)Fi=eHili.

Now that the li=Zi+1Zi is small enough, we can expand eHili by Taylor series(41)Fi=eHili=I+liHi+li22!Hi2+li33!Hi3+.

Equation (41) indicates that the series converges so that the method put forward is available. Then, vi+1=Fivi, in which Fi is the transfer matrix Ai as mentioned in Section 2.

The symplectic algorithm maintains strictly the symplectic structure of the Hamilton system and has unique advantages in long-term numerical stability . Therefore, the transfer symplectic matrix method proposed in this paper has better numerical stability and accuracy in solving the critical speed of the rotor.

4. Numerical Results

Traditional single stator and rotor motors, whether dc, asynchronous, or synchronous, have only one mechanical port. In recent years, the concept of dual-rotor motor has been put forward, which enables the design of the motor to realize the dual mechanical shaft and realize the independent transfer of energy on the two mechanical shafts. The dual-rotor motor can greatly reduce the actual size and weight of the motor and greatly improve the working efficiency of the motor, which can well meet the specific requirements of energy saving and speed regulation of the motor, showing better motor performance. Therefore, dual-rotor motor has been widely used in many fields and this paper will take the dual-rotor system as an example for verification analysis.

In this section, in order to validate the accuracy and efficiency of the method presented in this paper, a numerical example is given and the critical speeds of each order are compared with those presented in [26, 27] and other numerical algorithms; the numerical analysis and symbolic computations are performed in MATLAB.

As an example, we choose the two-rotor system with two intermediate supports, which is shown in Figure 4. Some material and geometrical properties are given as E = 206 GPa, μ = 0.3, and ρ = 7850 kg/m3. The detailed original data of the calculation example are shown in Table 1.

Computational model of a dual-rotor system.

Original data table.

Rotor number, N Revolving speed, ω (rad/s1) Equatorial moment of inertia, J×108 (m4) Unit number, N Quality, m (kg) Rotational inertia, Ip(kgm2) Spring stiffness, k×107 (Nm1) Length, l (m)
1 994.88 3.9760 1 0.325 0 2.5 0.075
2 10.45 0.1125 0.075
3 0.935 0 0.225
4 1.3205 0 0.225
5 0.925 0 0.075
6 8.45 0.09 0.075
7 0.305 0 1.5 0

2 1361.36 5.2002 8 0.1932 0 1.2 0.075
9 7.7796 0.0735 0.15
10 0.3278 0 0.15
11 6.2796 0.0432 0.075
12 0.1732 0 1.0 0

In order to establish the transfer matrix of the dual-rotor system, the rotor model in Figure 4 is divided into elements as shown in Figure 5, where 10, 11, 17, and 18 are imaginary elements which have no length and mass.

Cell partition diagram.

The critical speed results of the dual-rotor system calculated by the method in this paper are listed in Table 2. The exact solutions are obtained by Lagrange’s equation of motion.

First six-order critical speed results of various algorithms.

Order Proposed Reference Theoretical
1 400.75 400.97 400.72
2 894.42 894.42 894.41
3 1433.01 1434.26 1432.09
4 1829.76 1829.76 1829.74
5 2119.37 2119.42 2119.30
6 3057.71 3063.46 3056.65
7 3423.41 3429.56 3422.19
8 3723.98 3730.15 3721.94
9 4755.99 4763.05 4752.94
10 12046.95 12077.82 12033.93

It can be seen from the results in Table 2 that when calculating the first five-order critical speeds, the results of the two methods are close to the exact solutions. However, when exceeding the sixth critical speed, the numerical errors of the algorithm proposed in the literature are obviously increased, while the transfer symplectic matrix method proposed in this paper retains high accuracy.

Through the above analysis, it can be found that the difference between the method used in the literature and the method proposed in this paper is very little, and it is difficult to find the advantages and disadvantages of the two methods from the results. Therefore, the relative error analysis is made as follows. The relevant results are shown in Figure 6.

From the figure, we can find that the error curves of the two methods show an upward trend overall. The results are consistent with our cognition: when the rotation speed is high, the values of matrix elements may vary greatly, which directly leads to the reduction of calculation accuracy and the increase of error.

When the rotor speed reaches the third-, sixth-, and ninth-order critical speed, the relative errors between the numerical and exact solutions obtained by the two methods are significantly increased. This phenomenon may be caused by the structure of rotor itself.

When the rotor rotates at a low speed (below the fourth critical speed), both the reference method and the method proposed in this paper have satisfactory numerical results in terms of accuracy. However, when the rotor speed exceeds the fifth critical speed, the numerical stability of the method proposed in this paper is obviously better than that adopted in the literature. The relative errors are 0.0347% and 0.2228, respectively, at the sixth critical speed.

Comparison diagram of relative error between numerical solution and exact solution.

The calculated natural modes are shown in Figure 7 (positive vortex).

5. Conclusion

In this paper, the idea of finite element is included into the transfer matrix. Starting from the variational principle, dual variables are introduced to deduce the transfer symplectic matrix of gyro rotor under the Hamilton system.

Compared with the traditional transfer matrix method, a series of complex operations such as solving the aggregation model of beam are avoided, which is easy to understand.

This method is applicable to all kinds of boundary conditions and supporting conditions, so it has universal significance. The state vector reflects the mechanical characteristics of bearing rotor, and the physical meaning is clear and intuitive.

The simulation experimental results show that, in the calculation of high-order critical speed (exceed the fifth critical speed), the numerical accuracy of the transfer symplectic matrix method is obviously better than that of the reference method. Meanwhile, the numerical stability is improved compared with the reference method to some extent, which is reflected by the trend of the curve. The numerical example indicates the feasibility and superiority of the method, which provides the basis for the optimal design of the rotor system.

Data Availability

The data used to support the findings of this study have not been made available because the data that have been used are confidential and the authors do not have permission to share data.

Conflicts of Interest

The authors declare no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant nos. 61573012 and 61702388) and Fundamental Research Funds for Central University (WUT2016A004).

Chen Y. S. Huabiao Z. Review and prospect on the research of dynamics of complete aero-engine systems Acta Aeronautica Et Astronautica Sinica 2011 32 8 1371 1391 Tatar A. Schwingshackl C. W. Friswell M. I. Dynamic behaviour of three-dimensional planetary geared rotor systems Mechanism and Machine Theory 2019 134 39 56 10.1016/j.mechmachtheory.2018.12.0232-s2.0-85059040122 Changliang T. Jinfu Y. Dongjiang H. Huan L. Rotor dynamics research of the composite flywheel spin test system: modeling and simulation Proceedings 2016 IEEE International Conference on Power and Renewable Energy October 2017 Shanghai, China IEEE 10.1109/ICPRE.2016.78711332-s2.0-85017209620 Zhang A. P. Zhao C. B. Research on the verification method of the rotor dynamics model of the locomotive turbocharger Proceedings 2015 International Conference on Computer Science and Mechanical Automation October 2016 Hangzhou, China IEEE 318 321 10.1109/CSMA.2015.702-s2.0-84963529405 Saimi A. Hadjoui A. An engineering application of the h-p version of the finite elements method to the dynamics analysis of a symmetrical on-board rotor European Journal of Computational Mechanics 2016 25 5 388 416 10.1080/17797179.2016.12455972-s2.0-84992529774 Knotek J. Novotný P. Maršálek O. Multibody based tool for simulation of the turbocharger rotor dynamics Applied Mechanics and Materials 2016 821 1 229 234 10.4028/www.scientific.net/amm.821.229 Brouwer M. D. Sadeghi F. Ashtekar A. Archer J. Lancaster C. Combined explicit finite and discrete element methods for rotor bearing dynamic modeling Tribology Transactions 2015 58 2 300 315 10.1080/10402004.2014.9686992-s2.0-85032096265 Lee J. W. Lee J. Y. In-plane bending vibration analysis of a rotating beam with multiple edge cracks by using the transfer matrix method Meccanica 2017 52 4-5 1143 1157 10.1007/s11012-016-0449-42-s2.0-84968611362 Dongyang C. Abbas L. K. Xiaoting R. Qing X. Marzocca P. Dynamic modeling of sail mounted hydroplanes system-part I: modal characteristics from a transfer matrix method Ocean Engineering 2017 130 629 644 10.1016/j.oceaneng.2016.12.0202-s2.0-85008698067 Al Rjoub Y. S. Hamad A. G. Free vibration of functionally Euler-Bernoulli and Timoshenko graded porous beams using the transfer matrix method KSCE Journal of Civil Engineering 2017 21 3 792 806 10.1007/s12205-016-0149-62-s2.0-84977071808 Gu J. J. Rui X. T. Zhang J. S. Riccati transfer matrix method for linear tree multibody systems Journal of Applied Mechanics 2016 84 1 011008 10.1115/1.40348662-s2.0-84993940385 Gu J. J. Rui X. T. Zhang J. S. Chen G. Computation of eigenvalues of linear tree multibody system based on riccati transfer matrix method Journal of Nanjing University of Science and Technology 2018 42 1 8 14 Mao W. G. Han X. Liu G. P. Riccati transfer matrix method combined with newmark acceleration formulation integration for analysing sliding bearings and rotor system Journal of Vibration and Shock 2015 10 20 14 20 Mao W. G. Han X. Liu G. P. Rolling bearing-rotor system Riccati-Newmark acceleration transfer matrix method Journal of Vibration and Shock 2015 34 20 80 84 Xu Z. X. Zhou D. Y. Deng Z. C. Numerical method based on Hamilton system and symplectic algorithm to differential games Applied Mathematics and Mechanics 2006 27 3 341 346 10.1007/s10483-006-0309-y2-s2.0-33744736667 Hongqing Z. Alatancang The Hamiltonian system and completeness of symglectic orthogonal system Applied Mathematics and Mechanics 1997 18 3 237 242 10.1007/BF02453366 Hua-jiang O. Wan-xie Z. Hamiltonian system and simpletic geometry in mechanics of materials (III)-flexure and free vibration of plates Applied Mathematics and Mechanics 1993 14 1 19 23 10.1007/BF024512172-s2.0-51649140540 Hsieh S. C. Chen J. H. Lee A. C. A modified transfer matrix method for the coupling lateral and torsional vibrations of symmetric rotor-bearing systems Journal of Sound and Vibration 2008 289 1-2 294 333 10.1016/j.jsv.2005.02.0042-s2.0-27344459464 Sun Z. He Y. Zhao J. Shi Z. Zhao L. Yu S. Identification of active magnetic bearing system with a flexible rotor Mechanical Systems and Signal Processing 2014 49 1-2 302 316 10.1016/j.ymssp.2014.05.0042-s2.0-84902548858 Jiang S. Zheng S. A modeling approach for analysis and improvement of spindle-drawbar-bearing assembly dynamics International Journal of Machine Tools and Manufacture 2010 50 1 131 142 10.1016/j.ijmachtools.2009.08.0102-s2.0-70449530524 Jiang S. Zheng S. Dynamic design of a high-speed motorized spindle-bearing system Journal of Mechanical Design 2010 132 3 034501 034505 10.1115/1.40011092-s2.0-77955195776 Jiang S. Mao H. Investigation of variable optimum preload for a machine tool spindle International Journal of Machine Tools and Manufacture 2010 50 1 19 28 10.1016/j.ijmachtools.2009.10.0012-s2.0-70449518007 Aleyaasin M. Ebrahimi M. Whalley R. Vibration analysis of distributed-lumped rotor systems Computer Methods in Applied Mechanics and Engineering 2000 189 2 545 558 10.1016/s0045-7825(99)00308-42-s2.0-0034282608 Zhang J. J. Cui S. Feng Y. X. A transfer matrix method for a multi-support rotor system in the symplectic space and its application Journal of Vibration and Shock 2017 36 16 32 36 Kang F. Qin M. Symplectic Geometry Algorithm of Hamilton System 2003 Hangzhou, China Zhejiang Science and Technology Press Pang H. Li X. L. Wang S. M. Study on the calculation of critical speed of double rotors by the integral transfer matrix method based on Riccati transform Journal of Mechanical Science and Technology 2005 24 7 832 834 Aviation Industry Press General Editorial Board of Aeroengine Design Manual, Rotor Dynamics and Whole Machine Vibration 2000 Beijing, China Aviation Industry Press