Next Article in Journal
Review of Steady-State Electric Power Distribution System Datasets
Previous Article in Journal
Shared Solar and Battery Storage Configuration Effectiveness for Reducing the Grid Reliance of Apartment Complexes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Energy-Based Vibration Modeling and Solution of High-Speed Elevators Considering the Multi-Direction Coupling Property

1
State Key Laboratory of Fluid Power and Mechatronic Systems, Zhejiang University, Hangzhou 310027, China
2
Canny Elevator Co., Ltd., Suzhou 215000, China
*
Authors to whom correspondence should be addressed.
Energies 2020, 13(18), 4821; https://doi.org/10.3390/en13184821
Submission received: 9 August 2020 / Revised: 9 September 2020 / Accepted: 13 September 2020 / Published: 15 September 2020
(This article belongs to the Section A1: Smart Grids and Microgrids)

Abstract

:
Vibration of a high-speed elevator is one of the major factors contributing to a safe, stable and comfortable ride. Vibration control and systematic optimization rely on the establishment of an accurate and reliable model. An energy-based vibration model (EVM) is proposed in this paper to describe the multi-directional coupling properties of high-speed elevators. After analyzing the structural and kinematic characteristics, a combination of the kinematic energy, elastic potential energy and virtual works is implemented to describe the vibrational characteristics. The Gaussian precise integration method is implemented to obtain accurate solutions because of the infinite degrees of freedom and the time-varying parameters. An example study is conducted using the KLK2 high-speed elevator. EVM is established to obtain the M a x ( P k ) / P k and A95 indices of the X, Y, Z axis vibrational acceleration through simulations. A comparison is proposed between the proposed EVM method, conventional methods based on differential-equation vibration model (DVM) and real prototype experimental data. Results show that the proposed EVM method could achieve more accurate results than the conventional DVM method. Deviations of these indicators are less than 5%, which suggest the efficiency, precision and reliability of the proposed EVM.

Graphical Abstract

1. Introduction

Operating speed and lifting elevation have been rapidly increasing with the high-rise buildings springing up. Generally, many causes, including installation error of the guide rail, stretch of the traction rope, airflow disturbance in the lift well and systematic intrinsic frequency, have impacts on the vibration of the lift car, which would significantly affect riding comfortability and safety. Vibration analysis, control, reduction and optimization of high-speed elevators demand a comprehensive solution urgently, which have attracted attention in the past few decades.
Vibration modeling and simulation techniques of high-speed elevators have constructed a solid theoretic basis in the literature. Through analysis on the non-linear vibration property of a moving elevator string with a weighted end, a time-varying differential equation set was constructed in [1], and further approximately solved by Hamilton’s principle. The natural frequency and modes of the elevator traction system were analyzed in [2] after applying external excitation. The free vibration equation is constructed based on the obtained response characteristics, and were applied in normal operation and emergency braking tests. To avoid the potential danger caused by the natural frequency coincidence between the high-rise building and the traction rope, a comprehensive mathematical model is presented in [3] considering the combined lateral stiffness of the guiding system. Also targeting at the influence of the lateral vibration on the dynamic response of traction rope, a vibration model was derived in [4] with quadratic and cubic non-linear terms. Further, the influence of cage motion and rope stretch were also considered in [5] to construct a dynamic model for both vibration analysis and prediction. Modal analysis in the operating condition is conducted in [6] for stiffness coefficients estimation through minimizing the summed mode shape and frequency errors between the experiments and the analytical model. A forced vibration differential equation was proposed in [7,8] to model the non-linear characteristics of the rolling guide shoes in consideration of the random perturbation of the guide rails and uncertainty of the car parameters. Then, horizontal vibration response of the car is proved to be influenced by the irregularity of the guide rails. The nonlinear behavior between the guide roller and the guiding rails was observed and investigated in [9]. Subsequently, a mathematical method along with an efficient parameter estimation method was proposed. Considering the combined effect of the guide rail perturbation and air flow disturbance, a transverse vibration model with four degrees of freedom is proposed in [10] based on computational fluid dynamic analysis under different working conditions. An adaptive vibration model based on neural network was proposed in [11] for vibration estimation, evaluation and prediction. An active roller guide system was developed and fully researched in [12] to prove its outstanding performance of controling elevator vibration and reducing electric power consumption. The nonlinear oscillation and stability of the elevator was investigated in [13], and an analytical solution based on parameter expansion was provided after estimating the influence of different parameters on the systematic response.
Some existing researches have focused on vibration optimization and control techniques. Vibration cannot be eliminated totally and thoroughly, but could be controlled and reduced to increase riding comfortability for passengers. A reduced order control model of the ultra high speed elevator was generated in [14], and the performance was proved through numerical simulations. Lateral vibrations of the traction rope caused by wind perturbation were studied in [15] considering the rope stretch, with the assumption that the axial, vertical velocity was treated constant. A specially designed input device with gaps installed near the upper end of the traction rope was proved to be efficient and robust in vibration control of the elevator [16]. A new control suppressor presented in [17] made use of the control force and actuator stroke to suppress vibration. The efficiency of controlling the first and second mode lift car vibration has been proved. Position of the installed vibration suppressor was researched in [18] to obtain the actual resonant position. Finite difference analysis with or without vibration suppressor were conducted for comparison. Based on test results of acceleration responses performed on an experimental elevator model, a non-stationary mathematical model was generated in [19]. The parameters were obtained through the Finite Element Method. A vibration control strategy based on the state-dependent Ricatti equation was proposed in [20] to suppress the horizontal nonlinear behavior after investigating the tilting motion of the lift car. A control device composed of two rotary electric motors was proposed in [21] to limit the irregular horizontal motion of the lift car on the guide rail, while not introducing any additional vertical vibration. An active suspension system was developed and verified in [22,23] which is composed of six actuators to suppress the horizontal vibrations. To suppress the lateral vibration, an active roller guide system was proposed in [24,25] with significant advantages of controlling performance and electric consumption. It was observed in [26], that the fundamental vibration frequency of the front-back direction is lower than the lateral direction. The necessity of considering the coupling effect from both directions was emphasized. An innovative method of vibration control was proposed in [27] through optimizing the number of poles slots of a gearless elevator system. An extended full-order observer configuration was introduced in [28] to obtain a systematic observation simultaneously, based on which an acceleration feed compensation was implemented. Sway and stretch of both ropes were controlled in [29] while treating the systematic natural frequency and damping ratio as time varying parameters. Based on investigations on the nonlinearities of the horizontal response under rail perturbations, a linear quadratic regulator control strategy along with a magnetorheological damper was developed in [30] to adjust the control force while considering parametric errors and measurement noises. Numerical simulations were then implemented to verify the effectiveness and robustness in improving driving comfort.
Based on the kinematic differential equations, researches in the literature have reduced the impact of vibration on riding comfortability through optimizing design parameters and developing vibration reduction devices. Establishment of a comprehensive vibrational model would construct the theoretic basis of subsequent analysis for vibration patterns, influencing factors, and vibration-control-oriented optimization. Traditionally, vertical [28,31,32] and horizontal [20,33,34] vibrations are treated separately. The coupling effect has not been taken into account properly. The mass of the traction rope is usually neglected, which would cause damage to the modeling reliability because of its significant length. In addition, the commonly used partial differential equation sets have unlimited degrees of freedom, which makes it difficult to obtain accurate solutions. The uncertainty and complexity of external perturbances would lower the accuracy and reliability of vibration models formulated by kinematic differential equations.
As a solution, an energy-based vibration model (EVM) is established in this paper to describe the multi-direction coupling vibration property of high-speed traction system. Through a combination of the kinematic energy, elastic potential energy and damping virtual work, the mass of the traction rope is considered and the uncertain perturbances are avoided. Computational accuracy is promoted through the implementation of Gaussian precise integration method. A real example is analyzed using KLK2 high-speed elevator of a known manufacturer. The EVM is constructed to compute key indicators M a x ( P k ) / P k and A 95 of vibrational acceleration along X, Y, and Z axis. Finally, a comparison is shown between results obtained by EVM, traditional differential-equation vibration model (DVM) and a prototype. Solution results and experimental data on KLK2 prototype have shown the established EVM is more effective and efficient than the conventional DVM.

2. Structural and Vibrational Characteristics

2.1. Basic Structure of a High-Speed Traction System

The mechanical structures of high-speed traction systems are varied according to their designed speed and load capacity. One of the most commonly used structure is shown in Figure 1, in which some principle parts including the traction sheave, lift car, traction rope, rope tension device, guiding pulley and compensation rope are presented. The lift car is driven up and down by the friction effect between the traction sheave grooves and the traction rope. Connecting the lift car and the counterweight, the traction rope bears a position and time-varying axial tension and is extended dynamically. The high-speed traction system is a typical slender flexible system.

2.2. Analysis of Multi-Direction Coupling Vibration

In consideration of the traction system, axial stretch and sidewise swing of the traction rope are the direct causes of systematic vibration. To provide a precise description of the multi-direction coupling vibration while reducing computation, the following assumptions are made:
  • The impact of the rope stiffness on vibration could be neglected [1,35,36]. The traction rope is simplified to flexible strings.
  • Radial deformation scale could be neglected compared with axial stretch. Small elements on the traction rope share the same scalar of velocity, and physical parameters (such as elasticity modulus, linear density) are generally constant.
  • The guide rail and the elevator well do not belong to the traction system. Fabrication uncertainties, installation deviation of the guide rail, and airflow of the elevator well are not the direct causes of traction system vibration, which could be simplified as constant perturbations at the contact point between the traction sheave and the rope.
  • The traction sheave and the tension pulley are treated as mass points considering the significant dimensional difference between the length of the traction rope and the radius of pulleys. Therefore, the segmented strings are connected at the mass centers of the traction sheave and the tension pulley.
Structural simplification of the high-speed traction system is conducted based on those aforementioned rules. The systematic multi-direction coupling vibration characteristic is shown in Figure 2. Establish the Cartesian coordinate system at the mass center of the traction sheave. The vertical X axis describes the positive traction stretch; the Y axis is horizontal; the direction of the Z axis follows the right-hand rule. The traction rope is segmented to two strings, both of which bear axial tension. The lift car is considered as a rigid object centered m e , which is connected in both ends by the two strings. The impact of hydraulic guiding shoe, which limits the kinematic degree of freedom and provides lubrication, is simulated using a spring-damper system with a stiffness coefficient k and a damper coefficient c . s ( t ) represents an external perturbation along Y axis to simulate one of the actual vibrational sources in consideration of modeling reliability and computational convenience.
Assume the linear density of the traction rope is ρ , cross-section area is S , and elasticity modulus is E . l ( t ) is used to denote the length of the string at time t , and L is used to represent the projection distance along X axis between the traction sheave and the tension pulley. Consider an element P ( x , 0 , 0 ) at distance x away from the spatial origin. The projections of the actual displacements along X, Y and Z axis are respectively represented as u ( x , t ) , v ( x , t ) and w ( x , t ) . Therefore, the vertical velocity and acceleration of the lift car could be formulated as v ( t ) = l ˙ ( t ) and a ( t ) = v ˙ ( t ) using the differential operation.

3. Vibration Modeling Based on Energy Methods

3.1. Formulations of the Vibration Model Based on Energy Analysis

The vibrational models and the solutions of the traction system are similar regarding the acceleration and deacceleration stages in both height-lifting and height-reducing process. Only the displacement, velocity and acceleration vectors could be opposite in signs. Therefore, in consideration of limited space available, only the deacceleration stage in the height-lifting process is discussed and elaborated in this paper. To construct the vibration model, translation vector of an element P ( x , 0 , 0 ) is formulated as [1,37]:
R = [ x ( t ) + u ( x , t ) ] i + v ( x , t ) j + w ( x , t ) k ,
where i , j and k are used for translational unit vector along X, Y and Z axis respectively. Velocity vector V is then calculated through the deviation operation D on both sides:
V = [ v ( t ) + D u ( x , t ) D t ] i + D v ( x , t ) D t j + D w ( x , t ) D t k
The deviation operation D is defined by:
D D t = t + v ( t ) x
u , v and w are used to simply record u ( x , t ) , v ( x , t ) , w ( x , t ) , and subscript e is used to mark physical quantities regarding properties of the lift car. The translation vector and velocity vector of the lift car are then obtained:
R e = [ l ( t ) + u ( l , t ) ] i + v ( l , t ) j + w ( l , t ) k
V e = [ v ( t ) + D u ( l , t ) D t ] i + D v ( l , t ) D t j + D w ( l , t ) D t k
The kinematic energy of the traction system is formulated in:
E k = 1 2 m e V e 2 + 1 2 ρ 0 l V 2 d x + 1 2 ρ l H V 2 d x
where the kinematic energy of the lift car, traction string and tension string at time t are included respectively as the first, second and third term in Equation (6).
The traction string bears the gravity and pre-tension of the tension string, as well as the gravity of itself and the lift car. The tension string is only influenced by its gravity and the pre-tension force. The gravity acceleration g is a positive constant scalar quantity along the positive direction of X axis. The static tension of an element P ( x , 0 , 0 ) in the traction string at time t could be formulated as:
T ( x , t ) = { [ m + ρ ( l ( t ) x ) ] g + F ( 0 < x < l ( t ) ) ρ ( L x ) g + F ( l + ( t ) < x < L )
Then, the elastic potential energy U E is defined as [38]:
U E ( t ) = 1 2 k v 2 ( l , t ) + 1 2 k w 2 ( l , t ) + 0 l ( T 0 < x < l ( t ) ε + 1 2 E S ε 2 + 1 2 E l κ 2 ) 2 d x + l L ( T l + ( t ) < x < L ε + 1 2 E S ε 2 + 1 2 E l κ 2 )
where E S and E l represent the tensile strength and bending strength respectively. The deformation energy terms of the guiding shoe are the first and the second terms. The elastic potential energy of the traction string and the tension string caused by the static tension, vertical stretch and bending deformation are included as the third and the fourth terms. The parameters ε and κ respectively represent the positive strain and curvature defined in Equations (9) and (10), where subscript x represent the partial differential operation.
ε = u x + 1 2 v x 2 + 1 2 w x 2
κ = ( v x x 2 + w x x 2 )
The systematic gravity potential energy U g is formulated as:
U g = 0 l ρ g u ( x , t ) d x l L ρ g u ( x , t ) d x m e g u ( l , t )
The three terms in Equation (11) respectively represent the gravity potential energy of the traction string, tension string and the lift car. The signs of terms in Equation (11) are worth noticing, and they should be carefully handled to guarantee the correct sign of U g . Basically, U g and u ( x , t ) are position-dependent signed scalar quantities, i.e., the signs of U g and u ( x , t ) should be determined by the actual position of the lift car comparing with the datum position at time 0. The sign of g is constantly positive because the direction of the gravity acceleration is always along the positive direction of X axis. In the discussed height-lifting process, the function value of u ( x , t ) is always negative, but the U g is positive because the lift car is being elevated. Therefore, minus signs are applied to terms in Equation (11). Virtual work of the damping force along Y and Z axis are represented as the first and the second term respectively in Equation (12):
δ W ( t ) = c D v ( l , t ) D t δ v ( l , t ) c D w ( l , t ) D t δ w ( l , t )

3.2. Basic Energy-Based Vibration Model (EVM)

The basic energy-based vibration model without external perturbations is constructed firstly. Because of the systematic non-conservative property, the general Hamilton’s principle [1,5] is introduced for kinematic modeling. After combining Equations (6), (8), (11), and (12), formulation in Equation (13) could be obtained.
t 1 t 2 ( δ E k δ U E δ U g + δ W ) d t = 0
In consideration of the Leibnitz’s rule and the partial integration method, Equations (14) and (15) could be achieved.
{ 0 l D v D t δ v t d x + v ( D v D t δ v ) l + 0 l t ( D v D t ) δ v d x = t 0 l D v D t δ v d x 0 l D w D t δ w t d x + v ( D w D t δ w ) l + 0 l t ( D w D t ) δ w d x = t 0 l D w D t δ w d x 0 l ( v + D u D t ) δ w t d x + v [ ( v + D u D t ) δ u ] + 0 l t ( v + D u D t ) δ u d x = t + 0 l ( v + D u D t ) δ u d x ( 0 < x < l ( t ) )
{ l + L D v D t δ v t d x + l + L t ( D v D t ) δ v d x = t l + L D v D t δ v d x + v ( D v D t δ v ) l + l + L D w D t δ w t d x + l + L t ( D w D t ) δ w d x = t l + L D w D t δ w d x + w ( D w D t δ w ) l + l + L ( v + D u D t ) δ u t d x + v [ ( v + D u D t ) δ u ] + l + L t ( v + D u D t ) δ u d x = t l + L ( v + D u D t ) δ u d x ( l + ( t ) < x < L )
The temporal and geometric boundary conditions of Equations (14) and (15) are:
{ u ( 0 , t ) = 0 , δ u ( 0 , t ) = δ u ( x , t 1 ) = δ u ( x , t 2 ) = 0 u ( l , t ) = u ( l , t ) = u ( l + , t )
{ v ( l , t ) = v ( L , t ) = 0 , δ v ( 0 , t ) = δ v ( x , t 1 ) = δ v ( x , t 2 ) = 0 v ( l , t ) = v ( l , t ) = v ( l + , t )
{ w ( l , t ) = w ( L , t ) = 0 , δ w ( 0 , t ) = δ w ( x , t 1 ) = δ w ( x , t 2 ) = 0 w ( l , t ) = w ( l , t ) = w ( l + , t )
The δ u , δ v and δ w are non-zeros because they are individual differentials. Therefore, the basic EVM and boundary conditions are calculated and re-organized in Equations (19)–(23).
{ T 2 + E S ( u x + 1 2 v x 2 ) + c t ( u t + v u x ) + k t u = 0 , x = L ( T 1 + E S u x + 1 2 E S v 2 x ) l 1 + m ( u t t + 2 v u x x + a u x + v 2 u x x + a ) = ( T 2 + E S u x + 1 2 E S v 2 x ) + m g , x = l ( t )
m ( v t t + 2 v v x t + a v x + v 2 v x x ) + c ( v t + v v x ) + k v + ( T v x ) l ( T 2 v x ) l + + ( E S u x v x + 1 2 E S v x 3 E I v x x x ) l ( E S u x v x + 1 2 E S v x 3 E I v x x x ) l + = 0 , x = l ( t )
{ ρ ( u t t + 2 v u x t + a u x + v 2 u x x + a ) E S ( u x x + v x v x x ) T l ρ g = 0 , 0 < x < l ( t ) ρ ( u t t + 2 v u x t + a u x + v 2 u x x + a ) E S ( u x x + v x v x x ) T l ρ g = 0 , l + ( t ) < x < L
{ ρ ( v t t + 2 v v x t + a v x + v 2 v x x ) + x ( T l v x ) E S x ( u x v x + 1 2 v x 3 ) + E I v x x x = 0 , 0 < x < l ( t ) ρ ( v t t + 2 v v x t + a v x + v 2 v x x ) + x ( T l + v x ) E S x ( u x v x + 1 2 v x 3 ) + E I v x x x = 0 , l ( t ) < x < L
{ ρ ( v t t + 2 v w x t + a w x + v 2 w x x ) + x ( T l w x ) E S x ( u x w x + 1 2 w x 3 ) + E I v w x x x = 0 , 0 < x < l ( t ) ρ ( v t t + 2 v w x t + a w x + v 2 w x x ) + x ( T l + w x ) E S x ( u x w x + 1 2 w x 3 ) + E I w x x x = 0 , l ( t ) < x < L

3.3. EVM with External Perturbations

Generally, main causes of the Y axis vibration are: Y axis perturbations induced by resonance between traction motor and lift car, axial runout induced by the wear of traction motor, elastic deformation and stretch of the traction rope, loose joint nuts. Elastic deformation of the traction rope is already considered in Section 3.1, while the impact of other causes could be reduced by adjustment of intrinsic frequency, periodic inspection of traction motor and guiding system. Main causes of the Z axis vibration are unsmooth joints of the guide rail, installation deviations, straightness errors and deformation, circular error of the guiding pulley, and the squeeze effect of lift well airflow. Some of these aspects could only be slightly relieved through improving the fabrication and installation techniques, adjusting airflow characteristics. Considering the types and properties of the remaining non-negligible and inevitable causes, their impacts are simplified to a Y axis perturbation s ( t ) in this paper. Boundary condition at x = 0 could be obtained as:
v ( 0 , t ) = s ( t )
As a result, vibration boundary conditions described in Equation (17) are transformed into a non-homogeneous form v ( x , t ) , which is further decomposed to a combination of term v 1 ( x , t ) with homogeneous boundary conditions, and term v 2 ( x , t ) with non-homogeneous boundary conditions:
v ( x , t ) = v 1 ( x , t ) + v 2 ( x , t )
Substitute Equation (25) to Equation (23), EVM with external perturbations could be formulated as:
δ W ( t ) = c D v ( l , t ) D t δ v ( l , t ) c D w ( l , t ) D t δ w ( l , t )
{ ρ ( v 1 . t t + 2 v v 1 . x t + a v 1 . x + v 2 v 1 . x x ) + T l v 1 . x E S ( u x v 1 . x + u x v 1 . x x + 1 2 v 1 . x 2 ) + E I v 1 . x x x x = ρ ( v 2 . t t + 2 v v 2 . x t + a v 2 . x + v 2 v 2 . x x ) + E S ( u x x v 2 . x + 3 v 1 . x v 1 . x x v 2 . x + 3 v 1 . x v 2 . x v 2 . x x + u x v 2 . x x ) + E S ( 3 2 v 1 . x x v 2 x 2 + 3 2 v 1 . x v 2 x x + 3 2 v 2 x 2 v 2 . x x ) + T l v 2 . x + T l v 2 . x x E I v 2 . x x x x , 0 < x < l ( t ) ρ ( v 1 . t t + 2 v v 1 . x t + a v 1 . x + v 2 v 1 . x x ) + T l + v 1 . x E S ( u x v 1 . x + u x v 1 . x x + 1 2 v 1 x 2 ) + E I v 1 . x x x = ρ ( v 2 . t t + 2 v v 2 x t + a v 2 x + v 2 v 2 x x ) + E S ( u x x v 2 x + 3 v 1 x v 1 . x x v 2 x + 3 v 1 x v 2 x v 2 x x + u x v 2 x x ) + E S ( 3 2 v 1 . x x v 2 x 2 + 3 2 v 1 x 2 v 2 x x + 3 2 v 2 x 2 v 2 . x x ) + T l + v 2 x + T l + v 2 . x x E I v 2 , x x x x , l + ( t ) < x < L
{ ρ ( v 1 . t t + 2 v w 1 . x t + a w 1 . x + v 2 w 1 . x x ) + T l w 1 . x E S ( u x w 1 . x + u x w 1 . x x + 1 2 w 1 . x 2 ) + E I w 1 . x x x x = ρ ( w 2 . t t + 2 v w 2 . x t + a w 2 . x + v 2 w 2 . x x ) + E S ( u x x w 2 . x + 3 w 1 . x w 1 . x x w 2 . x + 3 w 1 . x w 2 . . x w 2 . x x + u x w 2 . x x ) + E S ( 3 2 w 1 . x x w 2 . x 2 + 3 2 w 1 . x 2 w 2 . x x + 3 2 w 2 . x 2 w 2 . x x ) + T l w 2 . x + T l w 2 . x x E I w 2 . x x x x , 0 < x < l ( t ) ρ ( w 1 . t t + 2 v w 1 . x t + a w 1 . x + v 2 w 1 . x x ) + T l + w 1 . x E S ( u x w 1 . x + u x w 1 . x x + 1 2 w 1 . x 2 ) + E I w 1 . x x x x = ρ ( w 2 . t t + 2 v w 2 . x t + a w 2 . x + v 2 w 2 . x x ) + E S ( u x x w 2 . x + 3 w 1 . x w 1 . x x w 2 . x + 3 w 1 . x w 2 . x w 2 . x x + u x w 2 . x x ) + E S ( 3 2 w 1 . x x w 2 . x 2 + 3 2 w 1 . x 2 w 2 . x x + 3 2 w 2 . x 2 w 2 . x x ) + T l + w 2 . x + T l + w 2 . x x E I w 2 . x x x x , l + ( t ) < x < L
Boundary conditions when x = l ( t ) is presented as:
{ m ( u t t + 2 v u x x + a u x + v 2 u x x + a ) m g + ( T 1 + E S u x + 1 2 E S v 2 x ) l 1 ( T 2 + E S u x + 1 2 E S v 2 x ) = 0 , x = l ( t ) T 2 + E S ( u x + 1 2 v x 2 ) + c t ( u t + v u x ) + k t u = 0 , x = L
Equations (26)–(28) are mathematical expressions of the established EVM with homogeneous boundary conditions, and Equation (29) is the boundary condition. Vibration response of the system under Y axis external perturbations is obtained after defining v 2 ( x , t ) using a linear interpolation:
v 2 ( x , t ) = a 0 ( t ) x x 1 x 0 x 1 + a 1 ( t ) x x 0 x 1 x 0
Based on the boundary conditions, we could see:
{ v 2 ( 0 , t ) = s ( t ) , v 2 ( l ( t ) , t ) = 0 .
Substitute Equation (31) to Equation (30), and the coefficients in Equation (30) are determined as a 0 ( t ) = s ( t ) , a 1 ( t ) = s ( t ) . Also, v 2 ( x , t ) could be obtained by:
v 2 ( x , t ) = s ( t ) [ 1 x l ( t ) ]
Comprehensively considering Equations (27)–(29) and (32), the EVM of high-speed elevator traction system under external perturbations is finally formulated.

4. Solution of EVM

4.1. Discretization of the Vibration Model

Assume there are n small elements in the traction rope based on time-varying element method at time t, and the length is l ( t ) . As stated in Section 2.2, projective displacements along X, Y and Z axis are denoted as v , u , w . Actual displacement of the j -th element could be estimated by multiplying shape function N v , j ( x , l ( t ) ) , N u , j ( x , l ( t ) ) , N w , j ( x , l ( t ) ) and displacement vector q v , j ( t ) , q u , j ( t ) , q w , j ( t ) :
{ v ( x , t ) = N v . j ( x , l ( t ) ) q v . j ( t ) x j < x < x j + 1 , 1 j n u ( x , t ) = N u . j ( x , l ( t ) ) q u . j ( t ) x j < x < x j + 1 , 1 j n w ( x , t ) = N w . j ( x , l ( t ) ) q w . j ( t ) x j < x < x j + 1 , 1 j n
Then, transform the shape functions to the form of Lagrange interpolating polynomials:
N v . j ( x , l ( t ) ) = N u . j ( x , l ( t ) ) = N w . j ( x , l ( t ) ) = [ x j + 1 x x j + 1 x j , x x j x j + 1 x j ] = [ j n x l ( t ) , n x l ( t ) j + 1 ]
The shape functions are recorded as N j   because N v , j ( x , l ( t ) ) = N u , j ( x , l ( t ) ) = N w , j ( x , l ( t ) ) . Basically, function N j ( x , l ( t ) ) is time-varying and position-dependent because the length of traction rope l ( t ) changes over time. Through partial derivatives on x and t , Equation (34) is approximately expanded in Equations (35)–(37) before the established EVM equations could be solved:
{ u x = N j q u . j , u x x = N j q u . j , u t = N j q ˙ u . j + v N j q u , j u x t = N j q ˙ u . j + v N j q u . j u t t = N j q ¨ u . j + v N j q ˙ u . j + v N j q ˙ u , j + v 2 N j q u , j + a N j q u , j
{ v x = N j q v . j , v x x = N j q v . j , v x x x = N j q v . j , v x x x N ( 4 ) j q v , j v x t = N j q ˙ v . j + v N j q v . j v x t = N j q ˙ v . j + v N j q v . j v t t = N j q ¨ v . j + v N j q ˙ v . j + a N j q ˙ v , j + v N j q v , j + v 2 N j q v , j
{ w x = N j q w . j , w x x = N j q w . j , w x x x = N j q w . j , w x x x N ( 4 ) j q w , j w t = N j q ˙ w . j + v N j q w . j , w x t = N j q ˙ w . j + w N j q w . j w t t = N j q ¨ w . j + v N j q ˙ w . j + a N j q ˙ w , j + v N j q w , j + v 2 N j q w , j
Substitute Equations (35)–(37) to the established EVM equations (Equations (19)–(23)), and we could get an equivalent form in Equation (38):
M ( t ) q ¨ + C ( t ) q ˙ + K ( t ) q + P ( t ) + F ( t ) = 0
where the M ( t ) , C ( t ) ,   K ( t ) , F ( t ) respectively denote the mass matrix, damping matrix, stiffness matrix and generalized force matrix. P ( t ) is a non-linear generalized coordinates coupling term. M ( t ) and K ( t ) are symmetric matrices. C ( t ) is a combination of systematic gyroscopic matrix and stamping matrix. Use superscript ‘to represent differential operation, and the expressions are shown in:
M ( t ) = [ M u 0 0 0 M v 0 0 0 M w ] ,   C ( t ) = [ C u 0 0 0 C v 0 0 0 C w ] ,   K ( t ) = [ K u 0 0 0 K v 0 0 0 K w ] P ( t ) = [ P u P v P w ] ,   F ( t ) = [ F u F v F w ] ,   q = [ q u q v q w ]
where
M u . i j = M v . i j = M w . i j = ρ δ i j + N i N j | x = l
C u . i j = 2 ρ v 0 l N i N j d x + 2 m v N i N j + c N i N j | x = L
C v . i j = C w . i j = 2 ρ v 0 l N i N j d x + 2 m v N i N j + c N i N j | x = l
K u . i j = ρ a 0 l N i N j d x + ρ v 2 0 l N i N j d x E S 0 l N i N j d x + m a N i N j | x = l + m v 2 N i N j | x = l + k 1 N i N j | x = L
K v . i j = K w . i j = ρ a 0 l N i N j d x + ρ v 2 0 l N i N j d x + ρ ( g a ) 0 l N i N j d x + m a N i N j | x = l E I 0 l N i N j d x   + F 0 l N i N j d x + m a N i N j | x = l + m v 2 N i N j | x = l + c v N i N j | x = l + k N i N j | x = l
P u . j = E S 2 0 l ( N i q u . i ) 2 N j d x
P v . j = E S 0 l ( N i q u . i ) ( N i q v . i ) N j d x + E S 0 l ( N i q v . i ) 3 N j d x
P w . j = E S 0 l ( N i q u . i ) ( N i q w . i ) N j d x + E S 0 l ( N i q w . i ) 3 N j d x
F u . j = ρ a 0 l N j d x m a N j | x = L F N j | x = L
F v . j = F w . j = 0
δ i j is the delta function, which is defined in Equation (39):
δ i j = { 0 , i j 1 , i = j

4.2. Solution of the Discretized EVM

4.2.1. Solving EVM through Gaussian Precise Integration

Basic equations of a discretized kinematic model are shown in Equations (40) and (41):
M X ¨ + C X ˙ + K X = F ( t )
X ( 0 ) = X 0 , X ˙ ( 0 ) = X ( 0 )
The mass matrix M is commonly definite. A dual variable method in reference to the Hamilton’s principle:
P = M X ˙ + C X / 2
which is further transformed to:
X ˙ = M 1 P M 1 C X / 2
Equation (40) could be transformed to:
P ˙ = ( K C M 1 C / 4 ) X C M 1 P / 2 + F
Combine Equations (43) and (44), and we could obtain:
[ X ˙ P ˙ ] = [ A D B G ] [ X P ] + [ 0 F ]
Equation (45) could be written as Equation (46), where v = [ X P ] , H = [ A D B G ] , r = [ 0 F ] , and A = M 1 C / 2 , B = ( K C M 1 C / 4 ) , G = C M 1 / 2 , D = M 1 :
v ˙ = H v + r
Equation (46) is the state equation of structural kinematic response, whose common solutions are in the form of Equation (47):
v ( t ) = e H t v 0 + 0 t e H ( t τ ) r ( τ ) d τ
The load bearing time are equally discretized to several time intervals Δ t . Any time could be represented as t k = k Δ t . Therefore, the state equation is presented as:
v ( t k ) = e H t k v 0 + 0 t k e H ( t k τ ) r ( τ ) d τ
v ( t k + 1 ) = e H t k + 1 v 0 + 0 t k + 1 e H ( t k + 1 τ ) r ( τ ) d τ
Combine Equations (48) and (49), the relation between v ( t k ) and v ( t k + 1 ) could be represented in Equation (50):
v ( t k + 1 ) = e H Δ t v ( t k ) + t k t k + 1 e H ( t k + 1 τ ) r ( τ ) d τ
Let:
T ( Δ t ) = e H Δ t
and Equation (49) turns into:
v ( t k + 1 ) = T ( Δ t ) v ( t k ) + t k t k + 1 e H ( t k + 1 τ ) r ( τ ) d τ
Integration term in Equation (52) could be represented in Equation (53) after implementing Gaussian integration method:
t k t k + 1 e H ( t k + 1 τ ) r ( τ ) d τ = Δ t 2 i = 1 n ω i e Δ t 2 ( 1 ξ i ) H r ( t k + Δ t 2 ( 1 + ξ i ) ) + O ( Δ t 2 n 1 ) = Δ t 2 i = 1 n ω i T ( Δ t 2 ( 1 ξ i ) ) r ( t k + Δ t 2 ( 1 + ξ i ) ) + O ( Δ t 2 n 1 )
where n represents the quantity of integration point, ξ i is their coordinates, ω i is the weight coefficient, and T ( Δ t ( 1 + ξ i ) / 2 ) is the index matrix. Substitute Equation (53) to Equation (52) would result in Equation (54), which is in the form of Gaussian precise integration [39]:
v ( t k + 1 ) = T ( Δ t ) v ( t k ) + Δ t 2 i = 1 n ω i T ( Δ t 2 ( 1 ξ i ) ) r ( t k + Δ t 2 ( 1 + ξ i ) ) + O ( Δ t 2 n 1 )
Computational precision of the obtained state vector v ( t k + 1 ) is determined by the quantity of integration point n and time step Δ t . It is obvious that there is a recurrence relation between v ( t k ) and v ( t k + 1 ) . Therefore, the remaining problem is the computation of the physical property index matrix (PPIM) T .

4.2.2. Computation of the physical property index matrix (PPIM)

The key step of computing the PPIM lies in the natural of the exponential functions:
e H τ = ( e ( H τ ) / m ) m
Assume m = 2 N and N = 20 . Δ t = τ / m is a minor temporal interval, in which Equations (56) and (57) are satisfied:
e H Δ t I + H Δ t + ( H Δ t ) 2 / 2 = I + T a
T a = H Δ t ( I + ( H Δ t ) / 2 )
The scale of T a is so small that a precision loss would be encountered because of computer storing roundoff if being added with the unit matrix. Therefore, T a obtained from Equation (57) rather than T a + I could be stored in the computing device. Reserved terms in Equations (56) and (57) are determined by the computational precision requirement.
Equation (55) is decomposed to obtain the T matrix:
T = [ T a + I ] 2 N = [ T a + I ] 2 N 1 [ T a + I ] 2 N 1
Because ( I + T b ) ( I + T c ) = I + T b + T c + T b T c , Equation (58) could be treated as a loop: for   ( i t e r = 0 ;   i t e r < N ;   i t e r + + ) :   T a = 2 T a + T a T a . After the loop terminates, Equation (59) is implemented to obtain the PPIM:
T = I + T a

5. Case Study

A case study is presented in this section taking the KLK2 high-speed elevator as an example, which is analyzed to constructed the EVM for simulation and computation. A comparison is then conducted between the obtained results and experimental data for elaboration and verification.

5.1. Operating Parameters of the Example Elevator

There are six steel traction ropes in KLK2 high-speed elevation system. The linear density of the steel ropes is ρ = 0.575   kg / m 2 , and the elastic modulus is E = 8 × 10 10   N / m 2 . The nominal diameter of the traction rope and tension rope is d = 13   mm . The tensile strength is E A = 4.72 × 10 6   N , and the bending strength is E I = 22.16   N · m 2 . Mass of the lift car is m e = 1250   kg . The pre-tension load is F = 750   N . The stiffness coefficient of pre-tension device k t = 2 × 10 3   N / m . The damping coefficient is c t = 25   N · s / m . Length of the lift well is l = 225   m , and the lifting travel is L = 200   m . Nominal operating time from the bottom level to the top level is t t o t a l = 30   s . Maximum operating speed and acceleration is v m a x = 10   m / s and a m a x = 1   m / s 2 . Take g = 9.81   m / s 2 .

5.2. Operating State Equation of the Example Elevator

There are basically seven phases between the activation and termination according to the variations of lifting speed and acceleration, as shown in Table 1. The activation time is recorded as t 0 , and the end time of the i -th phase is denoted as t i . The length, velocity, acceleration and jerk are represented as l ( t i ) , v e ( t i ) , a e ( t i ) and j e ( t i ) .
Change of jerk during acceleration phases shows a quadratic parabolic property. Because jerk is the third derivative of elevation, length of the traction rope could be represented as a quintic polynomial, which could be represented as:
l e ( i ) ( t ) = C 0 ( i ) + C 1 ( i ) ( t t i 1 ) + C 2 ( i ) ( t t i 1 ) 2 + C 3 ( i ) ( t t i 1 ) 3 + C 4 ( i ) ( t t i 1 ) 4 + C 5 ( i ) ( t t i 1 ) 5 .
Velocity, acceleration and jerk could be represented as the first, second and third derivatives of l e ( i ) ( t ) in Equation (60):
{ a e ( i ) ( t ) = 2 C 2 ( i ) + 6 C 3 ( i ) ( t t i 1 ) 1 + 12 C 4 ( i ) ( t t i 1 ) 2 + 20 C 5 ( i ) ( t t i 1 ) 3 j e ( i ) ( t ) = 6 C 3 ( i ) + 24 C 4 ( i ) ( t t i 1 ) 1 + 60 C 5 ( i ) ( t t i 1 ) 2 v e ( i ) ( t ) = C 1 ( i ) + 2 C 2 ( i ) ( t t i 1 ) 1 + 3 C 3 ( i ) ( t t i 1 ) 2 + 4 C 4 ( i ) ( t t i 1 ) 3 + 5 C 5 ( i ) ( t t i 1 ) 4
where t i 1 t t i , t [ 1 , 7 ] . C n ( i ) ( n = 0 , 1 , , 5 ) are undetermined coefficients.
Jerk in phase 1 could be described by j ( t ) = α ( t t 0 ) + β ( t t 0 ) 2 . β = α / ( t 1 t 0 ) because the jerk is zero at time t 1 . Assume the maximum acceleration is a m a x , the maximum velocity is v m a x , the initial length of the traction rope is l ( 0 ) , the elevation height is H, and the total operating time is t t o t a l . Operation state equation of all these seven phases could be presented:
{ l ( 1 ) ( t ) = l ( 0 ) + α ( t t 0 ) 4 24 α ( t t 0 ) 5 60 t j l ( 2 ) ( t ) = l ( t 1 ) + v ( t 1 ) ( t t 1 ) + a ( t 1 ) 2 ( t t 1 ) 2 l ( 3 ) ( t ) = l ( t 2 ) + v ( t 2 ) ( t t 2 ) + a ( t 2 ) 2 ( t t 2 ) 2 α 24 ( t t 2 ) 4 + α 60 t j ( t t 2 ) 5 l ( 4 ) ( t ) = l ( t 3 ) + v m a x ( t t 3 ) l ( 5 ) ( t ) = l ( t 4 ) + v ( t 4 ) ( t t 4 ) α 24 ( t t 4 ) 4 + α 60 t j ( t t 4 ) 5 l ( 6 ) ( t ) = l ( t 5 ) + v ( t 5 ) ( t t 5 ) + a ( t 5 ) 2 ( t t 5 ) 2 l ( 7 ) ( t ) = l ( t 6 ) + v ( t 6 ) ( t t 6 ) + a ( t 6 ) 2 ( t t 6 ) 2 + α 24 ( t t 6 ) 4 α 60 t j ( t t 6 ) 5
where a = ( 6 a m a x 3 v m a x 2 ) / ( t t o t a l v m a x a m a x a m a x l t o t a l v m a x 2 ) 2 , and t j = ( 6 a m a x ) / α .

5.3. EVM of KLK2 High-Speed Elevator

Substitute these parameters and operating state equations to Equations (28) and (29), and EVM of KLK2 high-speed elevator is obtained. Properties and types of EVM between phases are similar. EVM equations and boundary conditions of the first phase are listed in Equations (63)–(66) as an example:
{ 0.575 ( u t t + 2 v u x t + ( 0.074 t 3 0.33 t 2 ) u x + v 2 u x x + ( 0.074 t 3 0.33 t 2 ) ) 4.72 × 10 6 ( u x x + v 1 x v x x ) [ 24500 + 5.635 ( 225 + 0.0037 t 5 0.027 t 4 x ) ] 5.635 = 4.72 × 10 6 ( v 2 x v 1 . x x + v 1 . x v 2 x x + v 2 . x v 2 . x x ) , 0 < x < 225 + 0.0037 t 5 0.027 t 4 0.575 ( u t t + 2 v u x t + ( 0.074 t 3 0.33 t 2 ) u x + v 2 u x x + ( 0.074 t 3 0.33 t 2 ) ) 4.72 × 10 6 ( u x x + v 1 . x v x x ) [ 24500 + 5.635 ( 200 x ) ] 5.635 = 4.72 × 10 6 ( v 2 x v 1 . x x + v 1 . x v 2 . x x + v 2 x v 2 . x x ) , 225 + 0.0037 t 5 0.027 t 4 < x < 200
{ 0.575 ( v 1 . t + 2 v v 1 . x t + ( 0.074 t 3 0.33 t 2 ) v 1 . x + v 2 v 1 . x x ) + [ 24500 + 5.635 ( 225 + 0.0037 t 5 0.027 t 4 x ) ] v 1 x 4.72 × 10 6 ( u x v 1 . x + u x v 1 . x x + 1 2 v 1 . x 2 ) + 22.16 v 1 . x 0 x x = 0.575 ( v 2 . t t + 2 v v 2 . x t + a v 2 . x v 2 v 2 . x x ) + 4.72 × 10 6 ( u x x v 2 x + 3 v 1 x v 1 . x x v 2 x + 3 v 1 . x v 2 x v 2 . x x + u x v 2 x x ) + 4.72 × 10 6 ( 3 2 v 1 . x x v 2 . x 2 + 3 2 v 1 . x 2 v 2 . x x + 3 2 v 2 . x 2 v 2 . x x ) + [ 24500 + 5.635 ( 225 + 0.0037 t 2 0.027 t 4 x ) ] v 2 . x + [ 24500 + 5.635 ( 225 + 0.0037 t 2 0.027 t 4 x ) ] v 2 . x x 22.16 v 2 × 000 , 0 < x < 225 + 0.0037 t 5 0.027 t 4 0.575 ( v 1 . t + 2 v v 1 . x t + ( 0.074 t 3 0.33 t 2 ) v 1 . x + v 2 v 1 . x x ) + [ 24500 + 5.635 ( 200 x ) ] v 1 . x 4.72 × 10 6 ( u x v 1 . x + u x v 1 . x x + 1 2 v 1 . x 2 ) + 22.16 v 1 . x 0 x x = 0.575 ( v 2 , t + 2 v v 2 . x t + ( 0.074 t 3 0.33 t 2 ) v 2 . x + v 2 v 2 . x x + 4.72 × 10 6 ( u x x v 2 . x + 3 v 1 . x v 1 . x x v 2 . x + 3 v 1 x v 2 x v 2 x x + u x v 2 x x ) + 4.72 × 10 6 ( 3 2 v 1.2 x v 2 x 2 + 3 2 v 1 . x 2 v 2 x x 3 2 v 2 x 2 v 2 . x x ) + [ 24500 + 5.635 ( 200 x ) ] v 2 x + [ 24500 + 5.635 ( 200 x ) ] v 2 x x 22.16 v 2 , x 00 x , 225 + 0.0037 t 5 0.027 t 4 < x < 200
{ 0.575 ( w 1 . t + 2 v w 1 . x t + ( 0.074 t 3 0.33 t 2 ) w 1 . x + v 2 w 1 . x x ) + [ 24500 + 5.635 ( 225 + 0.0037 t 5 0.027 t 4 x ) ] w 1 . x 4.72 × 10 6 ( u x w 1 . x + u x w 1 . x x + 1 2 w 1 . x 2 ) + 22.16 w 1 . x 00 x = 0.575 ( w 2 t t + 2 v w 2 . x t + ( 0.074 t 3 0.33 t 2 ) w 2 . x v 2 w 2 . x x + 4.72 × 10 6 ( u x x w 2 x + 3 w 1 . x w 1 . x x w 2 . x + 3 w 1 . x w 2 . x w 2 . x x + u x w 2 . x x ) + 4.72 × 10 6 ( 3 2 w 1 . x x w 2 . x 2 + 3 2 w 1 . x 2 w 2 . x x + 3 2 w 2 . x 2 w 2 . x x ) + [ 24500 + 5.635 ( 225 + 0.0037 t 2 0.027 t 4 x ) ] w 2 . x + [ 24500 + 5.635 ( 225 + 0.0037 t 2 0.027 t 4 x ) ] w 2 . x x 22.16 w 2 . x 00 x , 0 < x < 225 + 0.0037 t 5 0.027 t 4 0.575 w ( v 1 . t + 2 v w 2 . x t + ( 0.074 t 3 0.33 t 2 ) w 1 . x + v 2 w 1 . x x ) + [ 24500 + 5.635 ( 200 x ) ] w 1 . x 4.72 × 10 6 ( u x w 1 . x + u x w 1 . x x + 1 2 w 1 . x 2 ) + 22.16 w 1 . x 0 x = 0.575 ( w 2 t t + 2 v w 2 . x t + ( 0.074 t 3 0.33 t 2 ) w 2 x + v 2 w 2 . x x + 4.72 × 10 6 ( u x w 2 x + 3 v 1 . x w 1 . x x w 2 x + 3 w 1 . x w 2 . x w 2 . x x + u x w 2 . x x ) + 4.72 × 10 6 ( 3 2 w 1 . x x w 2 x 2 + 3 2 w 1 . x 2 w 2 . x x 3 2 w 2 . x 2 w 2 . x x ) + [ 24500 + 5.635 ( 200 x ) ] w 2 . x + [ 24500 + 5.635 ( 200 x ) ] w 2 , x x 22.16 w 2 , x 0 x x , 225 + 0.0037 t 5 0.027 t 4 < x < 200
{ 2500 u t t + 2 v u x x + ( 0.074 t 3 0.33 t 2 ) u x + v 2 u x x + ( 0.074 t 3 0.33 t 2 ) 5.635 + [ 24500 + 5.635 ( 225 + 0.0037 t 5 0.027 t 4 x ) ] + 4.72 × 10 6 u x + 2.36 × 10 6 v 2 x ) r 1 [ 24500 + 5.635 ( 200 x ) ] + E S u x + 1 2 E S v x 2 = 0 , x = 225 + 0.0037 t 5 0.027 t 4 [ 24500 + 5.635 ( 200 x ) ] + 4.72 × 10 6 ( u x + 1 2 v x 2 ) + 25 ( u t + v u x ) + 2 × 10 5 u = 0 , x = 200

5.4. Solution of the Established EVM

Use the solution method proposed in Section 4.2 and segment the traction rope to 10 parts. Therefore, initial displacements of the discrete nodes along Y and Z axis are q v , j ( 0 ) = 0.1 sin ( π ( j 1 ) / n ) and q w , j ( 0 ) = 0.1 sin ( π ( j 1 ) / n ) respectively, when t = 0 and 1 j n . There are structural flexibility and remaining degrees of freedom because of loose joints in the fabrication and installation of the lift car, the guide rail, the guiding sheave and the traction rope. Therefore, after the simulation begins and before the elevation process is started, the lift car would initially tend to gain stability under the effect of self-gravity and free-hanging. The kinematic conditions would appear unstable for a short period of time before self-stability is reached. As the real starting point of time, the systematic response would be recorded 10 s after the simulation is started. The main reason of the 10-s gap is to leave some time for automatic balancing before valuable data could be produced and collected. To ensure the accuracy and reliability of experimental simulation and analysis, all those curves are recorded and presented after the 10-s gap. Take the time step as 0.01 s, and the vibrational acceleration at the joint point of the traction string and the lift car along X, Y, Z axis is obtained and presented in Figure 3, Figure 4 and Figure 5 respectively. Key indicators including M a x ( P k ) / P k and A95 could estimate the operating vibration [40]. As a result, M a x ( P k ) / P k along X, Y, Z axis are 0.376   m / s 2 , 0.317   m / s 2 , 0.503   m / s 2 respectively. A95 are 0.251   m / s 2 , 0.118   m / s 2 , 0.429   m / s 2 along X, Y, Z axis.

5.5. Comparison and Verification between Vibration Models

To verify the effectiveness and efficiency of the established EVM model, the DVM model [41] is built and simulated to obtain the response curve of vibrational acceleration at the joint point between the traction string and the lift car. As shown in Figure 6, Figure 7 and Figure 8, M a x ( P k ) / P k along X, Y, Z axis are 0.363   m / s 2 , 0.305   m / s 2 , 0.495   m / s 2 respectively. A95, along X, Y, Z axis are 0.245   m / s 2 , 0.114   m / s 2 , 0.417   m / s 2 .
Prototype experiments and measurements for the KLK2 high-speed elevator are conducted in a 288-meter elevator experimental tower to obtain the vibration along X, Y, Z axis. The experimental devices are shown in Figure 9. Figure 9a shows the elevator experimental tower. Figure 9b shows the PMT EVA-625 elevator vibration analysis system of physical measurement technologies (PMT) [42]. The testing and measuring site is shown in Figure 9c.
Vibrational acceleration curve obtained from PMT EVA-625 is shown in Figure 10, in which the four curves from top to bottom are (I) variation curve of the absolute velocity during the height lifting process, (II) variation curve of the vibrational acceleration along X axis, (III) variation curve of the vibrational acceleration along Y axis, (IV) variation curve of the vibrational acceleration along Z axis. M a x ( P k ) / P k along X, Y, Z axis are 0.384   m / s 2 , 0.324   m / s 2 , 0.512   m / s 2 respectively. A95, along X, Y, Z axis are 0.264 m / s 2 , 0.120   m / s 2 , 0.440   m / s 2 .
A comparison between the obtained indicators through DVM method, prototype experiments, and the proposed EVM method is shown in Table 2. M a x ( P k ) / P k obtained from the DVM method are 0.363   m / s 2 , 0.305   m / s 2 , 0.495   m / s 2 along X, Y, Z axis respectively, whose deviations from the experimental results are 5.47%, 5.86%, 3.32%. A95 are 0.245   m / s 2 , 0.114   m / s 2 , 0.417   m / s 2 along X, Y, Z axis, whose deviations from the experimental results are 7.20%, 5.00%, 5.23% respectively. As a contrary, M a x ( P k ) / P k obtained from the proposed EVM method are 0.376   m / s 2 , 0.317   m / s 2 , 0.503   m / s 2 along X, Y, Z axis, whose deviations from the experimental results are 2.08%, 2.16%, 1.76%. A95 are 0.251   m / s 2 , 0.118   m / s 2 , 0.429   m / s 2 , whose deviations from the experimental results are 4.92%, 1.67% and 2.50% respectively. Obviously, the proposed EVM method could achieve more precise results than the conventional DVM method. It is worth noting that the deviations from experimental data are less than 5%. The manufacturing, fabrication, installation error and the operating environment would significantly enhance systematic vibration. However, as detailed in Section 3.3, these influencing factors are not carefully considered when establishing EVM. The proposed EVM would probably appear more accurate and the simulation results would be closer to the experimental data when precision engineering and periodic inspection could relieve their impacts. Therefore, the example study has shown the efficiency, potential and reliability of the proposed EVM method.

6. Conclusions

An EVM method is proposed in this paper to describe the multi-direction coupling vibration properties of high-speed traction system, which would further benefit vibration control and optimization design. The main contributions are:
  • The characteristics of the traction system is analyzed and simulated using masses, springs and dampers. A combination of kinematic energy, elastic potential energy and virtual work provides a description of the multi-direction coupling vibration properties, based on which EVM method is established.
  • EVM contains a complicated partial differential equation set with unlimited degrees of freedom and time-varying parameters. This paper discretizes the EVM to an ordinary differential equation set through the time-varying element method, and obtain the solutions using the Gaussian precise integration method.
  • An example study is conducted using KLK2 high-speed elevator. EVM is then established to obtain the simulation results and the solutions of vibrational acceleration indices, M a x ( P k ) / P k and A95, along X, Y, Z axis. The obtained results are then compared with the conventional DVM method and experimental data of KLK2 prototype. The proposed EVM method could achieve more accurate results than the conventional DVM method, and deviations of these indicators are less than 5%. In this way, efficiency, potential and reliability of the proposed EVM method have been elaborated and verified.

Author Contributions

Conceptualization, L.Q.; data curation, Y.W.; formal analysis, L.Q., C.H. and Y.W.; funding acquisition, L.Q.; investigation, L.Q. and C.H.; methodology, L.Q.; project administration, G.Y. and S.Z.; resources, S.Z. and L.Z.; software, C.H. and Y.R.; validation, L.Z.; visualization, C.H. and Y.R.; writing—original draft, L.Q. and C.H.; writing—review & editing, G.Y. and Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been funded by the National Natural Science Foundation of China (51875516), the Jiangsu Province Science and Technology Achievement Transforming Fund Project (BA2018083).

Conflicts of Interest

No conflict of interest exits in the submission of this manuscript.

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon request.

References

  1. Fung, R.F.; Lin, J.H.; Yao, C.M. Vibration analysis and suppression control of an elevator string actuated by a PM synchronous servo motor. J. Sound Vib. 1997, 206, 399–423. [Google Scholar] [CrossRef]
  2. Peng, Q.; Jiang, A.; Yuan, H.; Huang, G.; He, S.; Li, S. Study on Theoretical Model and Test Method of Vertical Vibration of Elevator Traction System. Math. Probl. Eng. 2020, 2020. [Google Scholar] [CrossRef] [Green Version]
  3. Crespo, R.S.; Kaczmarczyk, S.; Picton, P.; Su, H. Modelling and simulation of a stationary high-rise elevator system to predict the dynamic interactions between its components. Int. J. Mech. Sci. 2018, 137, 24–45. [Google Scholar] [CrossRef]
  4. Kaczmarczyk, S.; Andrew, J.P.; Adams, J.P. The Modelling and Prediction of the Influence of Building Vibration on the Dynamic Response of Elevator Ropes. Mater. Sci. Forum 2003, 440–441, 489–496. [Google Scholar] [CrossRef]
  5. Yang, D.H.; Kim, K.Y.; Kwak, M.K.; Lee, S. Dynamic modeling and experiments on the coupled vibrations of building and elevator ropes. J. Sound Vib. 2017, 390, 164–191. [Google Scholar] [CrossRef]
  6. Kobayashi, S.; Yoshimura, T.; Noguchi, N.; Omiya, A. Estimation of dynamic characteristics of an elevator car using operational modal analysis. Nihon Kikai Gakkai Ronbunshu C Hen/Trans. Jpn. Soc. Mech. Eng. Part C 2008, 74, 548–553. [Google Scholar] [CrossRef] [Green Version]
  7. Zhang, R.; Wang, C.; Zhang, Q. Response analysis of the composite random vibration of a high-speed elevator considering the nonlinearity of guide shoe. J. Braz. Soc. Mech. Sci. Eng. 2018, 40. [Google Scholar] [CrossRef]
  8. Zhang, R.; Wang, C.; Zhang, Q.; Liu, J. Response analysis of non-linear compound random vibration of a high-speed elevator. J. Mech. Sci. Technol. 2019, 33, 51–63. [Google Scholar] [CrossRef]
  9. Wee, H.; Kim, Y.Y.; Jung, H.; Lee, G.N. Nonlinear rate-dependent stick-slip phenomena: Modeling and parameter estimation. Int. J. Solids Struct. 2001, 38, 1415–1431. [Google Scholar] [CrossRef]
  10. Yang, Z.; Zhang, Q.; Zhang, R.; Zhang, L. Transverse Vibration Response of a Super High-Speed Elevator under Air Disturbance. Int. J. Struct. Stab. Dyn. 2019, 19, 1950103. [Google Scholar] [CrossRef]
  11. Taplak, H.; Erkaya, S.; Yildirim, Ş.; Uzmay, I. The Use of Neural Network Predictors for Analyzing the Elevator Vibrations. Arab. J. Sci. Eng. 2014, 39, 1157–1170. [Google Scholar] [CrossRef]
  12. Utsunomiya, K.; Okamoto, K.I.; Yumura, T.; Funai, K.; Kuraoka, H. Active roller guide system for high-speed elevators. Elev. World 2002, 50, 198–205. [Google Scholar]
  13. Mirabdollah Yani, R.; Ghodsi, R.; Darabi, E. A closed-form solution for nonlinear oscillation and stability analyses of the elevator cable in a drum drive elevator system experiencing free vibration. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 4467–4484. [Google Scholar] [CrossRef]
  14. Teshima, N.; Kamimura, K.; Nagai, M.; Kou, S.; Kamada, T. Vibration Control of Ultra High Speed Elevator by Active Mass Damper: 1st Report, Study by Optimal Control Theory. Nihon Kikai Gakkai Ronbunshu, C Hen/Trans. Jpn. Soc. Mech. Eng. Part C 1999, 65, 3479–3485. [Google Scholar] [CrossRef]
  15. Gaiko, N.V.; van Horssen, W.T. Resonances and vibrations in an elevator cable system due to boundary sway. J. Sound Vib. 2018, 424, 272–292. [Google Scholar] [CrossRef]
  16. Otsuki, M.; Ushijima, Y.; Yoshida, K.; Kimura, H.; Nakagawa, T. Application of nonstationary sliding mode control to suppression of transverse vibration of elevator rope using input device with gaps. JSME Int. J. Ser. C Mech. Syst. Mach. Elem. Manuf. 2006, 49, 385–394. [Google Scholar] [CrossRef] [Green Version]
  17. Yamazaki, Y.; Tomisawa, M.; Okada, K.; Sugiyama, Y. Vibration control of super-high-speed elevators: (Car vibration control method based on computer simulation and experiment using a full-size car model). JSME Int. J. Ser. C Dyn. Control Robot. Des. Manuf. 1997, 40, 74–81. [Google Scholar] [CrossRef]
  18. Kimura, H.; Nakagawa, T. Vibration analysis of elevator rope with vibration suppressor. Nihon Kikai Gakkai Ronbunshu C Hen/Trans. Jpn. Soc. Mech. Eng. Part C 2005, 71, 442–447. [Google Scholar] [CrossRef]
  19. Arrasate, X.; Kaczmarczyk, S.; Almandoz, G.; Abete, J.M.; Isasa, I. The modelling, simulation and experimental testing of the dynamic responses of an elevator system. Mech. Syst. Signal Process. 2014, 42, 258–282. [Google Scholar] [CrossRef]
  20. Santo, D.R.; Balthazar, J.M.; Tusset, A.M.; Piccirilo, V.; Brasil, R.M.L.R.F.; Silveira, M. On nonlinear horizontal dynamics and vibrations control for high-speed elevators. JVC J. Vib. Control 2018, 24, 825–838. [Google Scholar] [CrossRef]
  21. Nakano, K.; Hayashi, R.; Suda, Y.; Noguchi, N.; Arakawa, A. Active Vibration Control of an Elevator Car Using Two Rotary Actuators. J. Syst. Des. Dyn. 2011, 5, 155–163. [Google Scholar] [CrossRef]
  22. Noguchi, N.; Arakawa, A.; Miyoshi, K.; Kawamura, Y.; Yoshimura, T. Active vibration control technology for elevator cars considering controllability. Nihon Kikai Gakkai Ronbunshu C Hen/Trans. Jpn. Soc. Mech. Eng. Part C 2013, 79, 3192–3205. [Google Scholar] [CrossRef] [Green Version]
  23. Noguchi, N.; Arakawa, A.; Miyata, K.; Yoshimura, T.; Shin, S. Study on Active Vibration Control for High-Speed Elevators. J. Syst. Des. Dyn. 2011, 5, 164–179. [Google Scholar] [CrossRef]
  24. Utsunomiya, K.; Okamoto, K.; Yumura, T.; Sakuma, Y. Vibration control of high-speed elevators taking account of electricity consumption reduction. Nihon Kikai Gakkai Ronbunshu C Hen/Trans. Jpn. Soc. Mech. Eng. Part C 2006, 72, 2048–2055. [Google Scholar] [CrossRef]
  25. Utsunomiya, K. Vibration Damping Device for an Elevator. J. Acoust. Soc. Am. 2011, 129, 3419. [Google Scholar] [CrossRef]
  26. Baek, K.-H.; Kim, K.-Y.; Kwak, M.-K. Dynamic Modeling and Controller Design for Active Control of High-speed Elevator Front-back Vibrations. Trans. Korean Soc. Noise Vib. Eng. 2011, 21, 74–80. [Google Scholar] [CrossRef] [Green Version]
  27. Kim, D.Y.; Park, M.R.; Sim, J.H.; Hong, J.P. Advanced Method of Selecting Number of Poles and Slots for Low-Frequency Vibration Reduction of Traction Motor for Elevator. IEEE/ASME Trans. Mechatron. 2017, 22, 1554–1562. [Google Scholar] [CrossRef]
  28. Kang, J.K.; Sul, S.K. Vertical-vibration control of elevator using estimated car acceleration feedback compensation. IEEE Trans. Ind. Electron. 2000, 47, 91–99. [Google Scholar] [CrossRef]
  29. Nguyen, T.X.; Miura, N.; Sone, A. Analysis and control of vibration of ropes in a high-rise elevator under earthquake excitation. Earthq. Eng. Eng. Vib. 2019, 18, 447–460. [Google Scholar] [CrossRef]
  30. Dos Santos, L.C.C.; Brasil, R.M.L.R.F.; Balthazar, J.M.; Piccirillo, V.; Tusset, A.M.; Santo, D.R. Active vibration control of an elevator system using magnetorheological damper actuator. Int. J. Nonlinear Dyn. Control 2017, 1, 114. [Google Scholar] [CrossRef]
  31. Fan, W.; Zhu, W.D. Dynamic analysis of an elevator traveling cable using a singularity-free beam formulation. In Proceedings of the ASME Design Engineering Technical Conference, Cleveland, OH, USA, 6–9 August 2017; Volume 6. [Google Scholar]
  32. Rijanto, E.; Muramatsu, T.; Tagawa, Y. Control of elevator having parametric vibration using LPV control method: Simulation study in the case of constant vertical velocity. In Proceedings of the IEEE Conference on Control Applications, Kohala Coast, HI, USA, 22–27 August 1999. [Google Scholar]
  33. Mutoh, N.; Kagomiya, K.; Kurosawa, T.; Konya, M.; Andoh, T. Horizontal Vibration Suppression Method Suitable for Super-High-Speed Elevators. IEEJ Trans. Ind. Appl. 1998, 118, 353–362. [Google Scholar] [CrossRef] [Green Version]
  34. Feng, Y.; Zhang, J.; Zhao, Y. Modeling and robust control of horizontal vibrations for high-speed elevator. JVC/J. Vib. Control 2009, 15, 1375–1396. [Google Scholar] [CrossRef]
  35. Zhang, P.; Bao, J.H.; Zhu, C.M. Dynamic analysis of hoisting viscous damping string with time-varying length. J. Phys. Conf. Ser. 2013, 448. [Google Scholar] [CrossRef]
  36. Kimura, H.; Ito, H.; Fujita, Y.; Nakagawa, T. Forced vibration analysis of an elevator rope with both ends moving. J. Vib. Acoust. Trans. ASME 2007, 129, 471–477. [Google Scholar] [CrossRef]
  37. Lu, Y. Dynamics of the Flexble Multibody Syestems; Higher Education Press: Beijing, China, 1996. [Google Scholar]
  38. Irgens, F. Continuum Mechanics; Springer: Berlin/Heidelberg, Germany, 2008; ISBN 978-3-540-74298-2. [Google Scholar]
  39. Wan-Xie, Z. On precise integration method. J. Comput. Appl. Math. 2004. [Google Scholar] [CrossRef] [Green Version]
  40. ISO 18738-1:2012, Lifts (Elevators)-Measurement of Ride Quality[S]. 2012. Available online: https://www.iso.org/standard/54395.html (accessed on 7 July 2020).
  41. Shabana, A.A. Dynamics of Multibody Systems, 2nd ed.; Cambridge University Press: Cambridge, UK, 1998; ISBN 13: 978-0521594462. [Google Scholar]
  42. The World’s Standard of Measurement Elevator Ride Quality, Vibration & Sound. Available online: https://www.pmtvib.com/eva-625 (accessed on 7 July 2020).
Figure 1. Basic structure of a high-speed traction system.
Figure 1. Basic structure of a high-speed traction system.
Energies 13 04821 g001
Figure 2. Multi-direction vibration characteristic of the high-speed traction system.
Figure 2. Multi-direction vibration characteristic of the high-speed traction system.
Energies 13 04821 g002
Figure 3. Vibrational acceleration response along X axis using EVM.
Figure 3. Vibrational acceleration response along X axis using EVM.
Energies 13 04821 g003
Figure 4. Vibrational acceleration response along Y axis using EVM.
Figure 4. Vibrational acceleration response along Y axis using EVM.
Energies 13 04821 g004
Figure 5. Vibrational acceleration response along Z axis using EVM.
Figure 5. Vibrational acceleration response along Z axis using EVM.
Energies 13 04821 g005
Figure 6. Vibrational acceleration response along X axis using DVM.
Figure 6. Vibrational acceleration response along X axis using DVM.
Energies 13 04821 g006
Figure 7. Vibrational acceleration response along Y axis using DVM.
Figure 7. Vibrational acceleration response along Y axis using DVM.
Energies 13 04821 g007
Figure 8. Vibrational acceleration response along Z axis using DVM.
Figure 8. Vibrational acceleration response along Z axis using DVM.
Energies 13 04821 g008
Figure 9. Prototype experiments and measurements for the KLK2 high-speed elevator.
Figure 9. Prototype experiments and measurements for the KLK2 high-speed elevator.
Energies 13 04821 g009
Figure 10. Curves of lifting velocity and vibrational acceleration obtained from PMT EVA-625.
Figure 10. Curves of lifting velocity and vibrational acceleration obtained from PMT EVA-625.
Energies 13 04821 g010
Table 1. Seven phases of high-speed elevation operation state.
Table 1. Seven phases of high-speed elevation operation state.
PhaseTime UsedDescription
1 t 1 t 0 Acceleration increases to a e = a m a x > 0
2 t 2 t 1 Acceleration remains constant.
3 t 3 t 2 Acceleration reduces to a e = 0 . Velocity increases to v e = v m a x > 0 .
4 t 4 t 3 Velocity remains constant.
5 t 5 t 4 Deacceleration decreases to a e = a m a x < 0
6 t 6 t 5 Deacceleration remains constant.
7 t 6 t 7 Deacceleration increases to a e = 0 . Velocity increases to v e = 0 .
Table 2. A comparison of the obtained indicators through the DVM method, prototype experiments and the proposed EVM method.
Table 2. A comparison of the obtained indicators through the DVM method, prototype experiments and the proposed EVM method.
DVM Method ( m / s 2 ) EVM Method ( m / s 2 ) Prototype Experiments ( m / s 2 )
M a x ( P k ) / P k X axis0.3630.3760.384
Y axis0.3050.3170.324
Z axis0.4950.5030.512
A95X axis0.2450.2510.264
Y axis0.1140.1180.120
Z axis0.4170.4290.440

Share and Cite

MDPI and ACS Style

Qiu, L.; He, C.; Yi, G.; Zhang, S.; Wang, Y.; Rao, Y.; Zhang, L. Energy-Based Vibration Modeling and Solution of High-Speed Elevators Considering the Multi-Direction Coupling Property. Energies 2020, 13, 4821. https://doi.org/10.3390/en13184821

AMA Style

Qiu L, He C, Yi G, Zhang S, Wang Y, Rao Y, Zhang L. Energy-Based Vibration Modeling and Solution of High-Speed Elevators Considering the Multi-Direction Coupling Property. Energies. 2020; 13(18):4821. https://doi.org/10.3390/en13184821

Chicago/Turabian Style

Qiu, Lemiao, Ci He, Guodong Yi, Shuyou Zhang, Yang Wang, Yong Rao, and Lichun Zhang. 2020. "Energy-Based Vibration Modeling and Solution of High-Speed Elevators Considering the Multi-Direction Coupling Property" Energies 13, no. 18: 4821. https://doi.org/10.3390/en13184821

APA Style

Qiu, L., He, C., Yi, G., Zhang, S., Wang, Y., Rao, Y., & Zhang, L. (2020). Energy-Based Vibration Modeling and Solution of High-Speed Elevators Considering the Multi-Direction Coupling Property. Energies, 13(18), 4821. https://doi.org/10.3390/en13184821

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop