Study on dynamic characteristics of double cylinder double acting bilge pump transmission

Jiangmin Xia1 , Zhuying Li2 , Guodong Wei4 , Zhenhai Zhang3

1, 2Department of Naval Architect Engineering, Naval University of Engineering, Wuhan, Hubei Province, China

4, 3College of Power Engineering, Naval University of Engineering, Wuhan, Hubei Province, China

1Corresponding author

Vibroengineering PROCEDIA, Vol. 14, 2017, p. 1-6.
Received 1 September 2017; accepted 9 September 2017; published 21 October 2017

Copyright © 2017 - JVE International Ltd.


In view of the complex structure of the ship bilge pump, the excitation source in the process of work is difficult to accurately determine the fault condition of the bilge pump through the vibration test data. Based on the theory of multi - body dynamics, rotor dynamics and electrical mechanics, the kinematics and kinetic equations of the bilge pump drive are established, and the dynamics, kinematics and vibration characteristics of bottom Pump drive analysis of the cabin is carried out. The frequency of the excitation force corresponding to the transmission mechanism is deduced, and the frequency of the excitation force of the main transmission is calculated.

Keywords: bilge pump, drive mechanism, kinematics dynamics, vibration.

1. Introduction

Bilge pump is an important device to maintain the vitality of the ship, vibration abnormalities will seriously affect the stealth performance of ships [1]. It is helpful to judge the performance status of the equipment in time and to accurately locate the faulty parts of the vibration and noise spectrum of the normal and abnormal state, and the main factors that affect the vibration and noise performance of the pump [2].

At present, the main technical problem is that the vibration characteristics of the bilge pump are not accurate, and it is difficult to accurately judge the technical status of the equipment. It is necessary to study the kinematics, dynamics and excitation characteristics of the equipment itself, including the motor, rolling bearing and the centrifugal force generated by the rotation of the crankshaft, by analyzing the characteristic frequency and the fault condition of the components of the equipment by the characteristic frequency. And then to establish the dynamic model and numerical calculation model of the mechanical parts of the bilge pump mechanical system. The research results can be used as the important basis for the monitoring of the bilge pump state and the vibration and noise control during the whole life period [3].

2. Double cylinder double acting bilge pump structure

A kind of bilge pump transmission system diagram shown in Fig. 1, the prime mover armature through the bearing installed in the motor body to drive the worm rotation, driven by the worm gear rotation, and the worm wheel and crankshaft meshing drive crankshaft rotation. The excitation source of the device itself is the centrifugal force of the motor, rolling bearing, gear transmission, piston inertia force, piston imbalance leading to crankshaft rotation [4].

2.1. Research on mathematical model of bilge pump drive

The pump transmission structure is divided into four parts: armature, crankshaft, connecting rod and piston. The mechanical model and equation of each part are established according to the interconnection of each part [5-7]. The armature rotor is simplified as a multi-stage continuous mass ladder axis, the establishment of armature structure model shown in Fig. 2. The reference base and the conjoined base are established based on the center of the bearing of the lower end of the motor. It is written as er=xr,yr,zr and eb=xb,yb,zb. q1 of the armature coil in the magnetic field, the inertia force q2 caused by the rotor center of mass from the rotating shaft; the binding force includes the restraining force of the lower bearing to the shaft FO=Fox,Foy,Foz and FB=FBx,FBy,FBz [8].

Fig. 1. The transmission system diagram

Fig. 2. The simplified model of armature rotor system

The armature rotor transmits the interference force to the body through the bearing at both ends. Therefore, the analysis of the dynamic characteristics of the armature rotor is to determine the binding force FO and FB. The differential equations of motion of the armature rotor are established by using the Timoshenko beam – axis model under the influence of the moment of inertia of the rotor and the shear deformation of the shaft on the dynamic characteristics of the motor armature. The motor armature rotor and its output shaft is approximately divided into five sections, namely the lower end of the shaft to the lower bearing section, the lower bearing to the lower end surface of the armature rotor, the armature rotor, the armature rotor upper surface to the upper bearing section, the upper end of the shaft is numbered j= 1, 2, 3, 4, 5.

The deformation curves of the neutral axis are xz,t, yz,t, the cross section is θxz,t, θyz,t. To express in a plural number:

w = x + i y θ = θ x + i θ y .

The external force and shear stress Q are also expressed in the plural form:

M = M x + i M y ,           q = q x + i q y ,             Q = Q x + i Q y .

Use Timoshenko beam – axis model, the first – order differential equation of the axis is:

4 w 1 z 4 - ρ E 1 + E μ G 4 w 1 z 2 t 2 + i ρ Ω E 3 w 1 z 2 t + ρ A i E I 2 w 1 t 2 - i ρ 2 Ω E μ G ,

where ρ is the axis density, E is the axial modulus, G is the axial shear modulus, μ is the shear deformation coefficient, and Ω is the axis rotation angle. The lower boundary condition is:

θ 1 - w 1 z = 0 ,           θ 1 z = 0 .

Since the lower bearing is simplified to elastic support and its stiffness is kA, the boundary condition at the junction is:

μ G A 1 θ 1 l 1 - w 1 l 1 z - μ G A 2 θ 2 0 - w 2 0 z - k A w 2 0 = 0 , E I 2 θ 2 0 z - E I 1 θ 1 l 1 z = 0 , w 1 l 1 = w 2 0 = w , θ 1 l 1 = θ 2 0 ,

where A1and A2 are the cross-sectional areas of the first and second segments.

The armature is simplified as a clustering parameter model. As shown in Fig. 3, the translational displacement of the armature centroid C is xc, yc, and the deviation angle around the centroid c is θx, θy. The axial center-of-mass displacement of the both ends of the armature is x1, y1 and x2, y2, respectively:

  w 2 l 2 , t = x 1 + i y 1 ,           w 4 0 , t = x 2 + i y 2 .

Fig. 3. The armature rotor model

Fig. 4. The crank – slider body composition diagram

From the second law of rigid body Newton, we get the differential equation of motion of centroid c:

m x ¨ c + k 11 x D + k 12 θ y c = F x , m y ¨ c + k 21 y D + k 22 θ x c = F y ,

where k11, k12, k21, k22 is the bending stiffness, Fx, Fy is the electromagnetic force caused by the electromagnetic field. From the momentum theorem, the differential equation of centroid rotation is obtained:

J d θ x ¨ + L θ y ˙ - k 32 y + k 33 θ x = 0 , J d θ y ¨ - L θ x ˙ + k 41 x + k 44 θ y = 0 ,

where k32, k33, k41, k44 for the axis bending stiffness. The relationship between the coordinates xD, yD of the point D in Fig. 3 and the coordinates of the two ends is:

x D θ y y D θ x = b / l a / l - 1 / l 1 / l b / l - 1 / l a / l 1 / l x 1 x 2 y 1 y 2 .

Let the eccentricity of the disc be e the initial phase is φ, the relationship between the coordinate xD, yD and the center of mass C is:

x c = x D + e c o s Ω t + φ , y c = y D + e s i n Ω t + φ .

Substituting Eq. (10) into Eq. (7):

m x ¨ D + k 11 x D + k 12 θ y c = F x + m e Ω 2 c o s Ω t + φ , m y ¨ D + k 21 y D + k 22 θ x c = F y + m e Ω 2 s i n Ω t + φ .

Eqs. (3)-(11) are the differential equations and boundary conditions of the motor rotor steady-state motion. From Fig. 4, the crankshaft, connecting rod and slider as a whole, by the base, crank, connecting rod, slider, guide rod and piston composed of six rigid body, and the centroid Ci of the conjoined group is established. The vector diameters of the connected base points Ci are ei=xi,yiT. The angle between the base vector xi and the base vector xr of the reference base is the attitude angle φi of the rigid body. The vector of the rigid body Bi is formed by the vector ri and the posture angle φi of each rigid body, and the coordinate array of the rigid body Bi is:

q i = r i T , φ i T = x i , y i , φ i T ,           i = 1 ,   2 , , 6 .

The coordinate array of the system consists of 18 coordinates:

q = q 1 T , q 2 T               q 6 T T = x 1 , y 1 , φ 1 , x 2 , y 2 , φ 2   ,       x 6 , y 6 , φ 6 T .

The system has a total of four hinges (H1~H5) and two sliding hinges (H6~H7). Due to the slider B5, the guide rod B4 is in the same form as the slider B3, thus adding two constraints. The total number of system constrained equations is 17, and the degree of freedom δ= 18-17 = 1. As the crank B1 by the worm, to a certain angular velocity ω movement, for the additional drive constraints, so the system for the whole system. In the system, the dots on the base B6 and the crankshaft B1 are constrained by the rotation hinges, as shown in Fig. 5. Let the coordinates of the point on the conjoined body e1 is -a1,0T.

Is the direction cosine matrix of the conjoined base e1 relative to the reference coordinate base er:

A r 1 = c o s φ 1 - s i n φ 1 s i n φ 1 c o s φ 1 .

The coordinates of the crank B1 on the reference coordinate base er are:

r O = r 1 + A r 1 ρ O 1 .

The coordinate relationship is explicit:

x O r y O r = x 1 y 1 + c o s φ 1 - s i n φ 1 s i n φ 1 c o s φ 1 - a 1 0 = x 1 - a 1 c o s φ 1 y 1 - a 1 s i n φ 1 .

Thus the H1 constraint equation is obtained:

x 1 - a 1 c o s φ 1 y 1 - a 1 s i n φ 1 = 0 0 .

In the same way, we can get the constrain equation of H2~ H8.From the function Φq,t, the acceleration constraint equation is:

Φ q q , t q ¨ = - Φ q q ˙ q q ˙ - 2 Φ q t q ˙ - Φ t t .

Fig. 5. The base and crank restraints diagram

Fig. 6. The crank – slider mechanical model diagram

2.2. Research on the dynamic characteristics of crank – slider mechanism

Bilge pump crank – slider mechanism, including slider, connecting rod, crankshaft and other three main components. Crank movement of the slider mechanism generated by the interference force transmitted through the bearing to the body, causing the pump vibration, so the purpose of the study is to grasp the bearing on the crankshaft of the reaction force. The components are used as rigid bodies, both ends of the rolling bearing is simplified as a hinge, and the crank-slider movement is simplified as a plane motion [9].

First, the crank-slider mechanical system as the object, the force shown in Fig. 6, the main force of the system, including the worm on the worm force F1a, the hydraulic pressure on the piston force F5a, the Rigid body mass mig, rigid body motion process inertia force miai. The binding force includes the force of the bearing on the crankshaft F1n, the force of the body on the slider F3n, the body of the piston force F5n and force even M5n.

The crank and the slider are connected by a connecting rod, and the connecting rod can be simplified as a two-pole structure. The force of the worm on the worm wheel F1a along the tangential direction of the crank, the connecting rod on the crank force F1c along the connecting rod direction. The momentum theorem and the momentum moment of the moving point are applied to the crank centroid C1 to obtain the crank kinetic equation:

m 1 x ¨ 1 = F 1 x c - F 1 x a + F 1 x n ,                                                     m 1 y ¨ 1 = - F 1 y c + F 1 y a + F 1 y n - m 1 g ,                 J 1 φ ¨ 1 = M 1 a F 1 a + M 1 c F 1 c + M 1 n F 1 n ,


Z 1 = m 1 m 1 J 1 ,             q 1 = x 1 y 1 φ 1 ,             F ^ 1 a = F 1 x a F 1 y a M 1 a ,             F ^ 1 n = F 1 x n F 1 y n M 1 n ,           F ^ 1 C = F 1 x C F 1 y C M 1 C .

The Eq. (21) is written in the form of a matrix:

Z 1 q ¨ 1 F ^ 1 a + F ^ 1 n + F ^ 1 C .

3. Conclusions

1) Considering the influence of the moment of inertia of the rotor and the shear deformation of the shaft on the dynamic characteristics of the motor armature, the kinematic and kinetic differential equations of the armature rotor can be established by using the Timoshenko beam-axis model.

2) By moving the moving parts as rigid body, the two ends of the rolling bearing are simplified for the hinge, and the crank-slider movement is simplified to the plane movement. The seven parts of the rolling bearing, front and rear ends, cranks, connecting rod journals, worm gears and so on are modeled. The dynamics of the mechanism and the constraints are in line with the actual situation.


  1. Yuan Donghong, Hua Suobaom, Gu Zhehong, Pan Zhengguang, Lu Fei Analysis of vibration and noise mechanism of bilge pump and measures for vibration and noise reduction. Chinese Ship Research, Vol. 4, Issue 5, 2009, p. 75-80.
  2. Cui Hongfeng, Li Guanghua Vibration analysis and precautionary measures for boat bilge pump. China Ship Repair, Vol. 6, 2005, p. 26-27.
  3. Zhang Yunxia Design and Structural Design of Linear Motor Bilge Pump. Master’s Degree Thesis of China University of Petroleum, 2007.
  4. Zhang Junhong Multi-Body Dynamics Analysis of Vibration and Noise of Internal Combustion Engine. Technology of Tianjin University, 2012.
  5. Liu Yanzhu, Hong Jiazhen, Yang Haixing Multi-Rigid Body System Dynamics. Higher Education Press, Beijing, 1989.
  6. Haug E. J. Computer Aided Kinematics and Dynamics of Mechanical Systems: Basic Methods. Needham Heights, Allyn and Bacon, 1989, p. 199-237.
  7. Banerjee A. K. Lemak M E. Multi-flexible body dynamics capturing motion induced stiffness. Journal of Applied Mechanicac, Vol. 58, 1991, p. 766-75.
  8. Tan Wenjie Multi-Body System Dynamics Analysis Method for Vibration Response of Internal Combustion Engine Crankshaft. Master’s Thesis of Vehicle and Traffic Engineering, Beijing Institute of Technology, 2012.
  9. Ho J. Y. L. Direct path method for flexible multibody spacecraft dynamics. Journal of Spacecraft and Rocket, Vol. 14, Issue 2, 1977, p. 102-110.