Research Article  Open Access
Nonlinear Vibration of Ladle Crane due to a Moving Trolley
Abstract
The structural vibration of the main beam of a crane causes fatigue damage and discomfort to the driver. The swing of the payload has an effect on positioning precision, especially for a ladle crane, and this directly affects production safety. To study the influence of system parameters on the vibration of a craneās main beam and the angle of the payload, a system consisting of the main beam, trolley, payload, and cabin was constructed. A rigidflexible coupling dynamic model of a moving trolley with a hanging payload that moves on the flexible main beam with a concentrated cabin mass is established, and the direct integration method is used to solve the nonlinear differential equations of system vibration, which are obtained through Lagrangeās equation. Then, the time domain responses of the flexible main beam, payload angle, and cabin vibration are obtained. The influences of the trolley running speed, quality of the payload, and quality and position of the cabin on the vibration of the main beam and payload angle are analyzed. The results indicate that the amplitude of the main beam is directly proportional to the quality of the trolley, payload, and cab; the position of the cabin is closer to the midspan; the amplitude of the main beam is larger; the structural damping has some influence on the vibration of the main beam; and the swing angle of the payload is related to the maximum running speed of the trolley, acceleration time, and length of the wire rope. In order to reduce the vibration of the main beam and cabin, the connection stiffness of the cabin should be ensured during installation.
1. Introduction
The ladle crane is a type of overhead crane and carries out onerous work under extremely hostile environments. Its mechanism and structure can withstand the strong impact vibration of the trolley on the main beam. Many scholars have simplified the vibration system model to the moving load beam coupling system and have carried out a significant amount of research on the dynamic response of beams under moving mass.
Michaltsos et al. [1] dealt with the linear dynamic response of a simply supported uniform beam under a moving load of constant magnitude and velocity by including the effect of its mass. By using a series solution for the dynamic deflection in terms of normal modes, the individual and coupling effect of moving load mass, its velocity, and other parameters were fully assessed.
ÅimÅek [2] investigated the dynamic characteristics of the main beam when the moving mass was running on a simply supported beam, established the equilibrium equations of the system by using Lagrangeās equations, and discussed the effects of shear deformation, various material distributions, velocity of the moving mass, inertia, Coriolis and the centripetal effects of moving mass on the dynamic displacements, and the stresses of the beam, in detail. Kiani et al. [3] investigated the maximum deflection and bending moment of the beams under various boundary conditions due to a moving mass by employing Hamiltonās principle and using EulerāBernoulli, Timoshenko, and thirdorder beam theories.
A comprehensive parametric investigation on the effects of moving mass weight and velocity on the dynamic behavior of a simply supported EulerāBernoulli was carried out by Nikkhoo et al. [4], by employing the eigenfunction expansion method. They introduced a concept of critical velocity in terms of beam fundamental period, midspan, and moving mass weight, in which the effect of convective accelerations in moving mass formulation was not negligible for masses moving with velocities greater than this critical one.
Esmailzadeh and Ghorashi [5] analyzed the dynamic behavior of beams with simply supported boundary conditions, which carried either uniform partially distributed moving masses or forces, and evaluated the critical speeds of the moving load, by considering the maximum deflection for the midspan of the beam. Moreover, they verified that the length of the distributed moving mass affects the dynamic response considerably.
Lee [6] investigated the dynamic response of a Timoshenko beam, which is acted upon by a moving mass, by using the Lagrange approach and the assumed mode method. He verified the results of the model by comparing them with the results of an equivalent moving load for the midspan deflection of a simply supported beam for a small number of moving mass weights and velocities as well as for different slenderness ratios of the base beam. Lou et al. [7] presented a finite element formulation of the Timoshenko beam subjected to a moving mass. Their results were in good agreement with those obtained by using the assumed mode method employed by Lee [6]. Many scholars have applied the model of mass moving on the beam to engineering practice. Lou et al. [8] and Cheng et al. [9] established a bridgetrackvehicle model in order to analyze the vibration of railway bridges under a moving train by taking into account the response of the track structure. The abovementioned studies have only investigated the dynamic characteristics of the simply supported beam under a moving mass, or the quality of moving mass, running speed, acceleration, or other factors. There is no model of moving mass that accounts for swing quality.
The dynamic response of cranes has been studied by many scholars. Niu et al. [10] presented a general mathematical modeling approach for electric crane system dynamics during operation, which could be used to analyze the dynamic responses of electromagnetism, mechanism, and structure, during the operation of electric cranes. Wu et al. [11ā14] established a mathematical model of the horizontal and vertical vibrations of a crane structure under the trolley running and lifting processes by using finite element theory and calculated the vibration response of the main beam. The correctness of the theoretical solution was verified by experimental verification, and the parameters were modified. Oguamanam et al. [15, 16] established an EulerāBernoulli equation in order to investigate the vibration response of a fixed crane beam. ZrniÄ et al. [17] constructed multidegreeoffreedom vibration models for crane structures in order to analyze the effects of wire rope length, damping ratio, and tilt angle of the sling load on structural vibration. Xin et al. [18] constructed a ninedegreeoffreedom mathematical model of a āhumanācraneārailā system and an annoyance rate model for use in the optimization of the structural parameters of overhead traveling cranes and used the particle swarm optimization algorithm to optimize the structural design of overhead cranes. Some scholars predicted crane fatigue life by analyzing dynamics of vibration system [19, 20].
Although the dynamic characteristics of the main beam and the swing angle of the crane have been studied; these two aspects have not been discussed simultaneously. The models reported in the literature contain only the main beam, trolley, and crane quality, without considering factors such as cabin quality, position, and structural damping of the main beam, which also affect the vibration of the main beam. The vibration time response of the cabin can be used to analyze the vibration of the human body during operation and thereby reduce the possibility of occupational disease. Therefore, this paper considers these factors in order to establish a dynamic model of the ladle crane vibration system. Based on the principle of energy conservation, the system vibration equation was established through Lagrangeās equation. In combination with a direct integral method for solving the approximate solution of nonlinear vibration, the influence of factors such as the quality (trolley, payload, and cabin), trolley running speed, length of wire rope, position and connection stiffness of cabin on the vibration of the main beam, and the influence of the swing of payload was analyzed. In this study, the vibration acceleration and amplitude time domain responses of the cabin were obtained, which provide a theoretical reference for the optimization of crane design.
2. Crane Vibration System Modeling
2.1. Description of Vibration System
The ladle floating crane includes a frame, trolley movement organization, cabin, and suspension system, and the ladle crane system vibration dynamics model is established according to the size and structural characteristics, based on the following assumptions:(1)The main beam is regarded as a flexible body, and only the vibration in the vertical direction is considered. If all the damping in the elastic body is viscous damping, the damping of the beam has very little influence on the vibration of the structure [21]; the damping ratio is usually 0.1%ā0.7%.(2)The cabin is simplified as a lumped mass, and only the vertical vibration is considered. The connection between the cabin and the main beam is simplified as a spring and damping model with greater stiffness.(3)According to the characteristics of the ladle crane, the rope is regarded as rigid without quality [22].(4)The payload of the crane is reduced to a swinging mass, which is suspended in mass blocks by a nonmass rigid rope, which moves with the mass in the plane.(5)The deflection of main beam is considered as linear, while the nonlinear effect of main beam deflection in the dynamic response analysis of the ladle crane is ignored.
The design method of the main beam consists of the allowable stress method and limit state method. The allowable stress method assumes that the deflection of the main beam is linear; however, the limit state method takes into consideration the nonlinear deflection of the main beam, while the fatigue analysis and the calculation of the crane beam structureās dynamic characteristics by these two methods have been reported in the literature [23, 24], where it was found that the results were basically the same with the exception of jib crane. Therefore, the nonlinear deflection of the beam had little effect on the dynamic response analysis of ladle cranes. In addition, after the completion of design and manufacturing, the beam of the bridge crane was formed on the precambering curve; therefore, the crane beam vibration amplitude was relatively small. Thus, this study adopted the mentioned hypothesis .
A schematic of the problem is depicted in Figure 1. The beam is simply supported and is assumed to be adequately modeled by using the EulerāBernoulli beam theory. Only the bending stiffness was calculated, and the shear stiffness and torsion stiffness of the beam were not considered. The properties of the beam are Youngās modulus , volume density , crosssectional area , length , and second moment of area . is the quality of the trolley, is the quality of the payload, is the quality of unit length for the main beam, is the length of the wire rope, is the connection stiffness between the cabin and the main beam, is the damping between the cabin and the main beam, is the quality of the cabin, is the mounting position coordinates of the cabin, is the gravitational acceleration, is the maximum running speed of the trolley, is the th modal damping ratio of the main beam structure, is the extreme left position of the trolley, is the extreme right position of the trolley, and is the swinging angle of the payload.
2.2. Dynamic Equations of the Vibration System
The trolley runs at the maximum speed on the main beam, and the elastic displacement curve of at moment in the direction of the main beam can be expressed as follows:where is the th modal of the simply supported beam, and and are the generalized coordinates and coordinate numbers of the elastic displacement of the main beam, respectively.
The coordinate vectors can be defined according to the simplified physical model shown in Figure 1. Therefore, the position vector of an elemental mass of beam , the position vector of carriage , the position vector of the payload , and the position vector of crane cabin can be expressed as
If the total kinetic energy of the system is , the energy in the two main beams is , the kinetic energy of the moving trolley is , and the kinetic energy of payload is , then the kinetic energy of each part of the system can be expressed as
In the coupled system model, the system potential energy includes the elastic strain energy of beams and the elastic potential energy of the spring between the cab and beam . In the process of movement, the work done by external forces (the gravity of cart, trolley, and payload) will also change the system potential.
The total potential energy for the system is
According to the principle of elastic mechanics, the main beam elastic strain energy is
The elastic potential energy of the spring between the cab and beam is
The work done by external forces (the gravity of cart, trolley, and payload) is
Equations (5), (6), and (7) are substituted into (4), and the potential energy in the system is expressed as
The dissipated energy in the system can be expressed as
The main beam damping in (9) can be obtained by (10) and (11). The th modal circular frequency of a simply supported beam with crosssection is
Thus, the structural damping of the th modal can be expressed aswhere is the damping ratio of th modal.
Lagrangeās equations of nonconservative systems are expressed as follows:where is the kinetic energy of the system, is the potential energy of the system, is the energy dissipation function of the system, is the damping force caused by the energy dissipation function , is the generalized excitation force of external action, is the generalized coordinate, and is the generalized velocity.
According to (12), the differential equations for the motion of the system can be derived as follows:
For a generic linear dynamic structural system, the equation of motion (see (13)) can be expressed in the following form:
In (16), , , and are mass, damping, and stiffness matrices of the system; , and are the acceleration, velocity, and displacement vectors; and is a timedependent loading vector. More details are included in Appendix.
2.3. Solving Method of Dynamic Response
The time domain response of the dynamic system can be obtained by using the direct integration method for a linear dynamic system with constant coefficients [18]. The Newmark method [25ā27] is a typical implicit solution approach to structural dynamics problems and has the advantages of good convergence, precision, and stability. Not only is the method condition stable, but its stability condition is also easily satisfied.
The vibration equations of the system ((13)ā(15)) constitute an dimensional twoorder nonlinear timevarying differential equation system, which cannot obtain the exact solution directly. Therefore, a direct integration method and a matrix iteration method are used in this paper to approximate the vibration response of the system.
It is assumed that the acceleration varies linearly within the interval , and that the Newmark scheme can be given in the following form:where and are parameters determined by the stability and integration accuracy, respectively. The Newmark family of the method is unconditionally stable if
For any timestep, if (17) are substituted into (16), the system equation can be obtained as
In (16), the equivalent stiffness and equivalent load can be expressed in the following form:
In (20), ; ; ; ;āā; . In this study, the parameters and were chosen for all examples.
3. Procedure of Solution
The procedures for calculating the dynamic responses of the main beam and a swinging object undergoing a moving trolley () are as follows:(1)The initial values of the vibration responses of the main beam , , and are taken as 0 in the calculation process, similar to the initial values of the payload swing angles , , and , and the initial values of the cabin vibration response , , and . The dynamic response , , and of (16) at moment is calculated by (17)ā(20).(2)The results obtained in step are substituted into (14) and (15). The payload swing angle responses , , and are calculated according to (17)ā(20). The cabin vibration response , , and is the same.(3)The results obtained by step are substituted into step in order to calculate the vibration response of the main beam at the next moment, and then step is repeated in order to iteratively calculate steps and , until a preset stop time is reached.(4)Determine whether the calculation results , , and are convergent. If not, return to step for further analysis and calculation.
The Newmark method used in this study should be unconditionally convergent, when calculating the numerical solution of system vibration [25, 26]. Therefore, it is necessary to determine the convergence of the calculation results through the time domain response of the vibration of the main beam. In engineering practice, when the trolley is running close to the midspan, the vertical central displacement of the main beam is larger, and when the midspan distance is greater, the vertical central displacement of the main beam is smaller. If the calculation results align with this situation, the result of the calculation converges [5, 14].
4. Numerical Results and Discussion
This study considered a 100/40āt and 28.5ām ladle crane made by Taiyuan Heavy Industry Limited Company as an example. During the service period of the crane, the rails were defective due to the influence of the foundation's differential settlement. In addition, there were vibrations of the human body and damage of equipment. The initial parameters of the crane vibration system are shown in Table 1.

4.1. System Model Validation
4.1.1. Beam Subjected to Moving Mass
To further confirm the reliability of the presented formulae and the developed computer programs, a pinnedpinned beam subjected to a moving mass and constant speed of was also investigated. The length of the beam was and the crosssectional area was , and second moment of area was . The mass density was and Youngās modulus was . The above information of the beam was almost the same as references [5, 14].
Figure 2 shows the time histories for the vertical central displacements of the pinnedpinned beam. In [14], the influence of the moving mass on the midpoint shape of the simply supported beam was solved by using finite element theory. Additionally, in reference [5], the dynamic response of the middle point of the main beam was solved under a modal assumption. However, in this study, the time domain response of the midpoint vibration of the main beam was solved by combining the modal hypothesis with the energy conservation of the system. As can be seen from Figure 2, three calculation results were consistent with the underlying trend. Since the differences between the last three curves were small, the formulations and computer programs developed in this study should be available for calculating the dynamic responses of a structure due to a moving load.
4.1.2. Swing Angle Verification
The payload swing is affected by the carriage traverse acceleration and the length of the wire rope. The reliability of the presented formulae and the developed computer programs were further confirmed with regard to the swing angle. A pinnedpinned beam was subjected to a moving mass of and a payload quality of . The beam had a length of crosssectional area of , and second moment of area . The mass density was and Youngās modulus was . In practice, the carriage starts with an initial velocity of zero and accelerates to a particular speed, which could be held constant for some time before decelerating to rest. In all cases, the terminal times of the intervals were , , and , the speed during the constant speed phase was , and the deceleration phase occurred over the last meter of the motion. The initial displacements for each swing angle were .
The quality of the cabin, the connection stiffness, and damping were all set to zero. By ignoring the structural damping effect, the factors were considered to be consistent with those reported in [16]. From Figure 3, it can be seen that the results of this study are similar to those reported in [16], which proves the validity of the calculation principle and procedure.
4.2. Dynamic Response of Crane Main Beam and Payload Swing Angles
The main beam on the side of the cabin was selected as the object of analysis. Therefore, the effect of trolley on the main beam was half of the trolleyās quality, and only half of the wire rope was applied to the main beam. According to the system vibration equation model, the factors affecting the dynamic response of the craneās main beam are related not only to the crosssection properties and dimension parameters of the main beam, but also to the quality of the trolley, running speed, and quality and position of the cabin. This study uses these aspects to carry out an analysis on the vibration of the main beam and payload swing angle response.
4.2.1. The Effect of Trolley Quality and Running Speed on Dynamic Response
In comparison to other cranes, the ladle crane lifting steel ladle is a highrisk item. Running too fast may directly lead to an excessively large swing angle, which greatly reduces the safety of the production process. Therefore, the trolley of the ladle crane is run slower than that of other cranes by up to approximately 80/min. In fact, for the trolley running speed patterns shown in Figure 4, the time for trolley acceleration and deceleration was set to 5ās. The track of the trolley runs from the limit position on one side to the limit position on the other side. The maximum running speed of the trolley was . This ensured the safety of the trolley during the operation. The responses in Figure 4 are based on the three operating patterns of the trolley.
Figure 5 shows the vertical central displacement response of the craneās main beam vibration for the three speed patterns of the trolley. When the trolley was running near the center of the main beam, the amplitude of the central displacement was the largest. The maximum running speeds of the trolley were 0.4ām/s, 0.8ām/s, and 1.2ām/s, while the maximum vertical central displacements were 3.7āmm, 3.8āmm, and 3.9āmm, respectively. Therefore, the faster the maximum speed of trolley in this range, the greater the vertical central displacement. However, the change was not obvious.
Figure 6 shows the time domain response of the payload swing angle at different running speeds of the trolley. In the pattern of speed , the swing angle first increased and then reduced in the 0ā5ās acceleration phase. In the 5ā20ās uniform running phase, the payload presented a reciprocating periodic swing. The swing period of the payload was . By substituting the initial parameters, the period was calculated as 4.7894ās, which is the same as that of the uniform running phase, thus proving that the results are correct. In 20ā25ās, the deceleration phase, payload, and swing angle first increased and then decreased. This is due to the slowdown in the uniform motion effect to increase the amplitude of the payload. The swing amplitude of the acceleration phase was , while the amplitude of the deceleration phase was . The swing amplitude of the deceleration phase was slightly larger than that of the acceleration phase. This is because the payload still had a swing speed at the end of the uniform motion, and the direction of the velocity coincided with the direction of the deceleration run.
When the maximum running speed of the trolley was 0.4ām/s, 0.8ām/s, and 1.2ām/s, the maximum swing angle in the acceleration phase was 0.17ārad, 0.32ārad, and 0.49ārad, respectively, and the deceleration phase was 0.18ārad, 0.38ārad, and 0.55ārad, respectively, as shown in Figure 6. Since the maximum speed of the trolley was faster, the amplitude of the payload increased significantly. This is because at the same acceleration or deceleration time (5ās), the faster the acceleration, the greater the swing angle of the payload. By changing the quality of the trolley and payload or the cabin stiffness value, the calculated results remain essentially unchanged. Thus, the swing angle of payload is related to the speed and acceleration of the trolley, regardless of the trolley quality and payload, or the cabin stiffness value.
It can be seen from Figure 7 that the trend of central displacement of the main beam for different trolley qualities remains essentially unchanged. For trolley qualities of 10āt, 20āt, and 30āt, the maximum deflection was 2.7āmm, 3.6āmm, and 4.4āmm, respectively. With the increase in trolley quality, there was a corresponding increase in vertical central displacement.
4.2.2. Influence of Payload Quality and Rope Length on Crane Dynamic Response
As shown in Figure 8, with the increase in payload quality, the deflection of central displacement of main beam increases, and the nonlinearity of vibration becomes more and more obvious.
As can be seen from Figure 9, the results are in agreement with the different lengths of the wire rope. Therefore, the length of the wire rope has little effect on the vibration of central displacement of the main beam. This is because the swing angle was very small, as was the change in energy of the wire rope length. The quality of the payload in the swing model plays a major role in the deformation of the main beam.
According to the swing period of the payload , the length of the wire rope affects the swing period of the payload. As can be seen in Figure 10, the payload swing angles of different rope lengths differ slightly in the 0ā5ās acceleration phase; the amplitude of swing angle increases with the rope length in the 5ā60ās uniform running phase; there is no necessary linear relationship between the swing amplitude and rope length in the 60ā65ās deceleration phase. The main reason for these phenomena is that rope length determines the swing period of the payload. In this example project, the acceleration and deceleration times of the trolley were designed as 5ās, while the swing period of the payload was 4.79ās, 5.35, 6.18ās, and 7.57ās. In the acceleration phase, when the incentive effect of the trolley on the payload disappeared, the initial position of the payload swing was very different. Therefore, there was a big difference between the swing amplitude in the uniform running phase. Similarly, due to the differences in the starting phase angle and the period in the deceleration phase, there was no linear change in the swing amplitude. In the design of the crane, it is necessary to investigate the swing period of the payload and the acceleration time of the trolley, so that the payload of the swing angle can be fundamentally reduced. The angle of the payload can be reduced by controlling the length of the wire rope.
4.2.3. Influence of Associated Parameters of Cabin on Dynamic Response of Crane
In this model, the cabin is reduced to a lumped mass. The quality of the cabin has some influence on the deformation of the main beam (see Figure 11). The amplitude of the vertical direction of the main beam increases as the quality of the cabin becomes heavier. However, the quality of the cabin was very light, in comparison to the quality of the main beam and trolley; therefore, the impact on the vibration of the main beam was smaller.
The cabin position is usually designed to be the cross end of the main beam of the crane, which reduces the effect of the vertical direction of the main beam. This helps to improve the comfort of the driver during operation and reduces the possibility of occupational diseases. The position of the cabin can be changed within a small range in the design process. The closer the cabin is to the midspan, the larger the amplitude in the vertical direction, and the higher the sensitivity (see Figure 12).
The influence of the connection stiffness on the central displacement of the main beam is shown in Figure 13. When the values of connection stiffness are and , the nonlinearity of the midspan vibration is not obvious; however, when the values of connection stiffness are and , it is obvious that the amplitude is inversely proportional to the connection stiffness. When the connection stiffness is small, the vibration of the cabin is strong and has a great influence on the vibration of the main beam. Therefore, in the installation of the cabin, it should be kept close to the main beam, and the space for activity should be avoided.
4.2.4. Structural Damping of Simply Supported Beam
As has been reported in the literature, although the structural damping of the main beam is small, it has a subtle effect on the vibration of the craneās main beam system. These subtle effects may cause problems in the crane design process. As can be seen in Figure 14, when , the vibration is nonlinear and the maximum amplitude is 3.8āmm. The vibration of the trolley, when it runs to the end of the main beam, is still nonlinear. When , the fluctuation of the vibration is reduced and the maximum amplitude is 3.55āmm. The damping dissipates some energy, which reduces the vibration energy.
4.3. Dynamic Response Analysis of Cabin
The vibration of the cabin can reflect the comfort of the driver during operation. In the model, the connection stiffness of the cabin is a direct cause of vibration. Figures 15ā17 show the vibration response under different values of connection stiffness. As can be seen from Figure 15, the stronger the connection stiffness of the cabin, the smaller the vibration amplitude of the cabin. When , the vibration amplitude of the cabin is basically zero. In addition, when , the amplitude of the vibration of the cabin fluctuates violently; the nonlinearity is obvious and the maximum amplitude is 7.7āmm. When the values of connection stiffness are and , the relative displacement of the cabin and the main beam is basically zero, and when , it is 5.2āmm (see Figure 16). As can be seen in Figure 17, the stronger the cabin connection stiffness, the slower the acceleration of the cabin. In summary, in order to reduce the vibration intensity of the cabin and to improve the operating comfort of the driver, the cabin must be well fixed at the end of the main beam.
4.4. Influence of Nonlinear Analysis on the Main Beam and Cabin
Nonlinearity is a common phenomenon in nature, and linearity is a special case of nonlinearity. The nonlinear theory for interpreting the physical phenomena of crane system vibration has universal significance. The nonlinear vibration equations of the rigidflexible coupling vibration system of cranes obtained in Section 2 include the nonlinear elements of the system. The motivation force in the Appendix is related not only to the mass of trolley and hoisting but also to the swing angle and vibration of the cabin, which indicates nonlinear characteristics that are more complex. In this section, the dynamic response effect of the crane beam and cabin in a linear system are calculated by ignoring the nonlinearities in the equation. By comparing the calculation results of the linear system and nonlinear system, the influence of nonlinear factors on the vibration of main beam and cabin can be seen.
For analyzing the influence of nonlinear factors on the vibration of a crane system, we selected the following parameters: trolley mass , payload quality , maximum running speed of trolley , and connection stiffness of cabin .
The vibration time domain responses of the main beam and cabin, in both linear and nonlinear systems, are shown in Figures 18 and 19, respectively. As can be seen from the figures, the curve of displacement response of the main beam and cabin for a nonlinear system is below that for a linear system. Therefore, nonlinear factors affect the vibration of the main beam and cabin. The maximum amplitudes of vibration of the main beam and cabin in a linear system are 10.1āmm and 5.1āmm, respectively, whereas the maximum amplitudes of vibration of the main beam and cabin in a nonlinear system are 10.7āmm and 5.5āmm, respectively. In a nonlinear system, the vibration responses of the main beam and cabin are increased by 5.9% and 7.8%, respectively. Ladle cranes are mainly used for lifting molten steel. When ladle cranes are fully loaded, the deformation of the main beam will reach its maximum. These effects of nonlinear factors may directly cause accidents and result in significant loss of personnel and property. Therefore, in the analysis and calculation process, it is very important to consider the nonlinear factors.
5. Conclusion
Based on the principle of system dynamics, a system vibration model of the flexible main beam, cabin, trolley, and payload is established, which can provide a theoretical basis for crane design. The vibration function of the main beam is obtained by the mode superposition method, and the system vibration equation is established by Lagrangeās equation. The dynamic responses of the main beam and payload are obtained from the approximate solution of nonlinear system vibration.
The vertical amplitude of the main beam is directly proportional to the quality of the trolley and payload. In the acceleration time of 5ās, when the maximum running speed of the trolley is 0.4ām/s, 0.8ām/s, and 1.2ām/s, the maximum amplitude of the midspan of the main beam is not significant. The quality of the cabin has little influence on the vibration of the main beam; however, the influence of the position is obvious. The closer the cabin is to the main beam, the larger the amplitude is. At the same time, when the connection stiffness of the cabin is weak, the amplitude fluctuation of the main beam is high, and the nonlinearity of the vibration is more obvious. In contrast, structural damping can only change the vibration amplitude of the main beam and the speed of energy attenuation within a certain range, which has little effect on the vibration of the main beam.
For the same acceleration and deceleration times, the swinging angle of the payload increases with the speed of the trolley. Since the length of the wire rope affects the swing period of the payload, the influence of acceleration and deceleration times on the payload swing angle is more complex and requires further research. As the next step, research on antiswing control, based on the flexibility of the wire rope, can be considered.
Appendix
In these matrices: ; ; ; ; ; ; ; .
Conflicts of Interest
All the authors do not have any possible conflicts of interest.
Acknowledgments
This paper is sponsored by the Fund for Shanxi ā1331Projectā Key Subjects Construction and Postgraduate Innovation Project in Shanxi Province (2017BY118). The data of this work is provided by Taiyuan Heavy Industry Co., Ltd.
References
 G. Michaltsos, D. Sophianopoulos, and A. N. Kounadis, āThe effect of a moving mass and other parameters on the dynamic response of a simply supported beam,ā Journal of Sound and Vibration, vol. 191, no. 3, pp. 357ā362, 1996. View at: Publisher Site  Google Scholar
 M. ÅimÅek, āVibration analysis of a functionally graded beam under a moving mass by using different beam theories,ā Composite Structures, vol. 92, no. 4, pp. 904ā917, 2010. View at: Publisher Site  Google Scholar
 K. Kiani, A. Nikkhoo, and B. Mehri, āPrediction capabilities of classical and shear deformable beam models excited by a moving mass,ā Journal of Sound and Vibration, vol. 320, no. 3, pp. 632ā648, 2009. View at: Publisher Site  Google Scholar
 A. Nikkhoo, F. R. Rofooei, and M. R. Shadnam, āDynamic behavior and modal control of beams under moving mass,ā Journal of Sound and Vibration, vol. 306, no. 35, pp. 712ā724, 2007. View at: Publisher Site  Google Scholar
 E. Esmailzadeh and M. Ghorashi, āVibration analysis of beams traversed by uniform partially distributed moving masses,ā Journal of Sound and Vibration, vol. 184, no. 1, pp. 9ā17, 1995. View at: Publisher Site  Google Scholar
 H. P. Lee, āThe dynamic response of a timoshenko beam subjected to a moving mass,ā Journal of Sound and Vibration, vol. 198, no. 2, pp. 249ā256, 1996. View at: Publisher Site  Google Scholar
 P. Lou, G.L. Dai, and Q.Y. Zeng, āFiniteelement analysis for a Timoshenko beam subjected to a moving mass,ā Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 220, no. 5, pp. 669ā678, 2006. View at: Publisher Site  Google Scholar
 P. Lou, Z. Yu, and F. T. K. Au, āRailbridge coupling element of unequal lengths for analysing traintrackbridge interaction systems,ā Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems, vol. 36, no. 4, pp. 1395ā1414, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 Y. S. Cheng, F. T. K. Au, and Y. K. Cheung, āVibration of railway bridges under a moving train by using bridgetrackvehicle element,ā Engineering Structures, vol. 23, no. 12, pp. 1597ā1606, 2001. View at: Publisher Site  Google Scholar
 C. M. Niu, H. W. Zhang, and H. Ouyang, āA comprehensive dynamic model of electric overhead cranes and the lifting operations,ā ARCHIVE Proceedings of the Institution of Mechanical Engineers Part C Journal of Mechanical Engineering Science, vol. 203210, pp. 1484ā1503, 19891996. View at: Google Scholar
 J.J. Wu, āFinite element modelling and experimental modal testing of a threedimensional framework,ā International Journal of Mechanical Sciences, vol. 46, no. 8, pp. 1245ā1266, 2004. View at: Publisher Site  Google Scholar
 J.J. Wu, āFinite element analysis and vibration testing of a threedimensional crane structure,ā Measurement, vol. 39, no. 8, pp. 740ā749, 2006. View at: Publisher Site  Google Scholar
 J. J. Wu, A. R. Whittaker, and M. P. Cartmell, āDynamic responses of structures to moving bodies using combined finite element and analytical methods,ā International Journal of Mechanical Sciences, vol. 43, no. 11, pp. 2555ā2579, 2001. View at: Publisher Site  Google Scholar
 J.J. Wu, āTransverse and longitudinal vibrations of a frame structure due to a moving trolley and the hoisted object using moving finite element,ā International Journal of Mechanical Sciences, vol. 50, no. 4, pp. 613ā625, 2008. View at: Publisher Site  Google Scholar
 D. C. D. Oguamanam, J. S. Hansen, and G. R. Heppler, āDynamic response of an overhead crane system,ā Journal of Sound and Vibration, vol. 213, no. 5, pp. 889ā906, 1998. View at: Publisher Site  Google Scholar
 D. C. D. Oguamanam, J. S. Hansen, and G. R. Heppler, āDynamics of a threedimensional overhead crane system,ā Journal of Sound and Vibration, vol. 242, no. 3, pp. 411ā426, 2001. View at: Publisher Site  Google Scholar
 N. D. ZrniÄ, V. M. GaÅ”iÄ, and S. M. BoÅ”njak, āDynamic responses of a gantry crane system due to a moving body considered as moving oscillator,ā Archives of Civil and Mechanical Engineering, vol. 15, no. 1, pp. 243ā250, 2015. View at: Publisher Site  Google Scholar
 Y. Xin, G. Xu, and N. Su, āDynamic optimization design of cranes based on humancranerail system dynamics and annoyance rate,ā Shock and Vibration, vol. 1, Article ID 8376058, pp. 1ā19, 2017. View at: Publisher Site  Google Scholar
 Q. Dong, G. N. Xu, and H. L. Ren, āResearch on risk of remanufacturing telescopic job structure system under potential multifailure mode correlation,ā Journal of Advanced Mechanical Design, Systems, and Manufacturing, vol. 11, no. 2, pp. 1ā18, 2017. View at: Publisher Site  Google Scholar
 Q. Dong, G. N. Xu, H. L. Ren, and A. H. Wang, āFatigue remaining life estimation for remanufacturing truck crane Jib structure based on random load spectrum,ā Fatigue & Fracture of Engineering Materials & Structures, vol. 40, no. 5, pp. 706ā731, 2017. View at: Publisher Site  Google Scholar
 G. T. Michaltsos, āDynamic behaviour of a singlespan beam subjected to loads moving with variable speeds,ā Journal of Sound and Vibration, vol. 258, no. 2, pp. 359ā372, 2002. View at: Publisher Site  Google Scholar
 Z. Jun and C. Zhijian, āThe research about horizontal moving rule of bridge crane and swing rule of grab,ā Journal of Shandong University (Natural Science Edition), vol. 33, no. 4, pp. 393ā397, 1998. View at: Google Scholar
 G. Xu, X. Fan, F. Lu, and H. Yang, āMonte Carlo simulation of fatigue reliabilities and initial cracks for welded box girders of crane,ā Journal of Mechanical Engineering, vol. 47, no. 20, pp. 41ā44, 2011. View at: Publisher Site  Google Scholar
 G. Xu and H. Li, āGeneral bridge crane structure CAD based on working stress method and limit state design method,ā Hoisting and Conveying Machinery, vol. 10, pp. 72ā78, 2008 (Chinese). View at: Google Scholar
 B. Zheng, X. Gao, and C. Zhang, āRadial integration BEM for vibration analysis of two and threedimensional elasticity structures,ā Applied Mathematics and Computation, vol. 277, pp. 111ā126, 2016. View at: Publisher Site  Google Scholar
 T. C. Chen, C. Su, Z. Q. Hu, and H. T. Ma, āAn explicit timedomain method in sensitivity analysis of nonstationary random responses,ā Journal of Vibration Engineering, vol. 28, no. 4, pp. 503ā509, 2015 (Chinese). View at: Google Scholar
 Z. Hu, C. Su, T. Chen, and H. Ma, āAn explicit timedomain approach for sensitivity analysis of nonstationary random vibration problems,ā Journal of Sound and Vibration, vol. 382, pp. 122ā139, 2016. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Yunsheng Xin 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.