# Dynamic Response of a Gantry Crane’s Beam Subjected to a Two-Axle Moving Trolley

Authors: Qingrong Chen; Wenming Cheng; Lingchong Gao; Run Du
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

In this paper, a novel moving load model for 2D crane systems of which the trolley has two axles is proposed. Based on this model, the dynamics of a 2D gantry crane, which is modelled as a simply supported Euler–Bernoulli beam carrying a two-axle trolley from which a single-pendulum payload is suspended, is studied. The proposed model was verified by comparing with two models in existing papers and can be considered as an extended version of comparative models. Then, the effect of the trolley’s axle base on the dynamic responses of the beam is studied. It can be observed that increasing the length of trolley’s axle base will decrease the deflection of the beam, and a larger initial swing angle will cause a larger deflection of the beam without controlling the swing of the payload.


---

## Body

## 1. Introduction

Industrial crane systems, which are a class of representative nonlinear systems, play an important role in modern industries. They can be modelled as a coupled system of moving payload and flexible beam. Vibration is a serious problem in the applications of crane systems because the structural flexibility will certainly affect the behaviour of the controller leading to precise positioning problems of the trolley and payload [1]. For instance, in the metallurgical industry, it has not only a high demand for horizontal positioning accuracy of the crane transport but also a high demand for height positioning accuracy of the payload. It is well known that a larger span will cause a larger deflection of the main beam of crane systems while the cross section and the payload mass keep unchanged. Therefore, with the increase of crane span, the dynamic responses of the main beam cannot be neglected.Many researchers used the assumed mode method to establish the dynamic model of overhead crane systems [2–5] and quayside container crane systems [6]. It is found that the value of the beam deflection increases by increasing the total mass of the trolley and payload [2–4]. In [2], the authors also found that, for given masses of the trolley and payload, the location and the value of the maximum beam deflection are dependent upon the trolley speed. In [4], it is observed that the mass of the driver’s cabin has little effect on the deflection of the beam, but its position has obvious effect. In [5], it is found that the swing response of the payload is not sensitive to the stiffness of the beam. Some researchers applied the finite element method to study the transverse and longitudinal vibrations of the flexible beam of crane systems [7, 8]. In [7], the cable is considered as an elastic body, and the authors found that increase of trolley speed and acceleration/deceleration do not have significant influence for vertical displacements, but have on horizontal displacements. In [8], the author presented a technique to replace the moving load by an equivalent moving finite element. They demonstrated the effect of the cable length, the trolley mass, and the velocity of the trolley on the deflection of the beam.In the aforementioned papers, the trolley is modelled as a particle/mass point. That means there is only one contact point between the trolley and the flexible beam. However, a trolley has at least four wheels in contact with the flexible beam. If we consider it as a planar model, it can be reduced to a two-axle model. In the field of vehicle-bridge coupling system, many researchers studied the dynamic responses of bridges/beams subjected to moving vehicles with two axles [9–11]. In these papers, vehicles are in contact with the beams through two wheels, and the assumed mode method is adopted to derive the deflections of the beams. Inspired by [9–11], we proposed a novel moving load model for 2D crane systems of which the trolley has two axles. This new model can provide more accurate results of the dynamic responses of the crane beam.The rest parts of the paper are structured as follows. Section2 presents the new dynamic model of the crane system. Section 3 validates the accuracy of the mathematical model by comparing it with several existing models and analyses the effect of the trolley’s axle base. Section 4 draws the conclusions of the entire paper.


## 2. Dynamic Modelling

### 2.1. Crane System Description

Figure1 shows the overview of a L-type single-beam gantry crane, and Figure 2 shows the 2D schematic of the L-type single-beam gantry crane. In the 2D schematic, the trolley has two wheels in contacting with the beam. The full length of the beam is L, the distance between two legs is L0, and the length of the cantilever is L1. When the trolley reaches the limit position on the cantilever, the distance between the centre of the trolley and the centre line of the supporting leg is L2. The mass per unit length of the beam is ρ, the cross-sectional moment of inertia is I, and Young’s modulus is E.Figure 1
The overview of a L-type single-beam gantry crane.Figure 2
The 2D schematic diagram of the L-type single-beam gantry crane.Figure3 shows the schematic diagram of the coupling dynamics modelling of the crane system. O1 and O2 are support points simplifying for the supported legs. A1 and A2 are the end points of the cantilever. B1 and B2 are the limit positions of the centre of the trolley. The mass of the trolley is denoted as mc, the axle base of the trolley is d, the contact forces between the trolley and the beam are P1 and P2, the length of the cable is l, the mass of the payload is mp, and the sway angle of the payload is θ.Figure 3
The schematic diagram of the coupling dynamics modelling of the crane system.


### 2.2. Modelling

The following assumptions are made before the modelling:(1)
Compared with the payload’s mass, the mass of the cable connecting the payload and trolley can be ignored, and the cable is always in a tight state due to the effect of the payload. Therefore, the cable is considered massless and rigid, and the payload sways as a single pendulum.(2)
The dynamic responds of the support legs are not considered, and support points of the main beam are simplified as pin-pin supports.(3)
The cantilevers (O1−A1 and O2−A2) are considered as rigid beams because their displacements subjected to the moving trolley are small. The part of the main beam between two legs (O1−O2) is modelled as a simply supported Euler–Bernoulli beam.(4)
The cross section of the beam is constant.(5)
It is assumed that the trolley is always in contact with the top beam. The centre of the trolley has a constant distanceh from the beam top surface. The rotating movement around the centre of mass of the trolley is not considered.In this paper, based on these assumptions, the part of the main beamO1−O2 is considered as simply supported, and the general solution of its shape function is(1)ϕx=C1cosβx+C2sinβx+C3coshβx+C4sinhβx,with pined-pined boundary conditions(2)ϕ0=0,ϕ″0=0,ϕL0=0,ϕ″L0=0.Then, theith model shape function of the simply supported beam is derived as(3)ϕix=siniπxL0.The vertical displacement of the beam can be expressed as(4)yb=wx,t=∑i=1nϕixqit.The coordinates of the centre of the trolley are(5)xc=12x1+x2,yc=12wx1,t+wx2,t+h=12∑i=1nϕix1+ϕix2qit+h.The horizontal velocity of the centre of the trolley is(6)x˙c=x˙1=x˙2.And the vertical velocity is(7)y˙c=12dwx1,tdt+dwx2,tdt=12∑i=1nx˙1ϕi′x1qi+ϕix1q˙i+x˙2ϕi′x2qi+ϕix2q˙i=12∑i=1nx˙1ϕi′x1+x˙2ϕi′x2qi+12∑i=1nϕix1+ϕix2q˙i.The acceleration of trolley is known and expressed as(8)x¨c=x¨1=x¨2=a.The kinetic energy of the flexible beam, trolley, and payload are denoted byTb, Tc, and Tp, respectively, which can be calculated as follows:(9)Tb=12∫0Lρ∂yb∂t2dx=12∫0Lρ∑i=1nϕixq˙it2dx=12∑i=1n∑j=1nq˙itq˙jt∫0Lρϕixϕjxdx,(10)Tc=12mcx˙c2+y˙c2,(11)Tp=12mpx˙c+lθ˙cosθ2+y˙c−lθ˙sinθ2.Then, the total kinetic energy of the coupled system can be obtained:(12)T=Tb+Tc+Tp.The potential energy of the flexible beam, trolley, and payload are denoted byVb, Vc, and Vp, respectively, which can be calculated as follows:(13)Vb=12∫0LEI∂2yb∂x22dx=12∫0LEI∑i=1nϕi″xqit2dx=12∑i=1n∑j−1nqitqjt∫0LEIϕi″xϕj″xdx,(14)Vc=−mcgyc,(15)Vp=−mpgyc+lcosθ.Then, the total potential energy of the coupled system can be obtained:(16)V=Vb+Vc+Vp.By applying the Lagrange equation,(17)ddt∂T∂q˙i−∂T∂qi+∂V∂qi=0,i=1,2,…,n,one can derive the dynamic equation of the generalized coordinate qi:(18)14mc+mpϕix1+ϕix2∑j=1nx¨1ϕj′x1+x˙12ϕj″x1+x¨2ϕj′x2+x˙22ϕj″x2qj+14mc+mpϕix1+ϕix2∑j=1n2x˙1ϕj′x1+x˙2ϕj′x2q˙j+ϕjx1+ϕjx2q¨j+∑j=1n∫0Lρϕixϕjxdxq¨jt−12mplϕix1+ϕix2θ¨sinθ+θ˙2cosθ+∑j=1n∫0LEIϕi″xϕj″xdxqjt−12mc+mpgϕix1+ϕix2=0.By applying the Lagrange equation,(19)ddt∂T∂θ˙−∂T∂θ+∂V∂θ=0,one can derive the dynamic equation of the generalized coordinate θ.(20)−12mplsinθ·∑j=1nx¨1ϕj′x1+x˙12ϕj″x1+x¨2ϕj′x2+x˙22ϕj″x2qj−12mplsinθ·∑j=1n2x˙1ϕj′x1+x˙2ϕj′x2q˙j+ϕjx1+ϕjx2q¨j+mplx¨ccosθ+lθ¨+mpgl1sinθ=0.Equations (18) and (20) compose the dynamic equation of the coupled system.


## 2.1. Crane System Description

Figure1 shows the overview of a L-type single-beam gantry crane, and Figure 2 shows the 2D schematic of the L-type single-beam gantry crane. In the 2D schematic, the trolley has two wheels in contacting with the beam. The full length of the beam is L, the distance between two legs is L0, and the length of the cantilever is L1. When the trolley reaches the limit position on the cantilever, the distance between the centre of the trolley and the centre line of the supporting leg is L2. The mass per unit length of the beam is ρ, the cross-sectional moment of inertia is I, and Young’s modulus is E.Figure 1
The overview of a L-type single-beam gantry crane.Figure 2
The 2D schematic diagram of the L-type single-beam gantry crane.Figure3 shows the schematic diagram of the coupling dynamics modelling of the crane system. O1 and O2 are support points simplifying for the supported legs. A1 and A2 are the end points of the cantilever. B1 and B2 are the limit positions of the centre of the trolley. The mass of the trolley is denoted as mc, the axle base of the trolley is d, the contact forces between the trolley and the beam are P1 and P2, the length of the cable is l, the mass of the payload is mp, and the sway angle of the payload is θ.Figure 3
The schematic diagram of the coupling dynamics modelling of the crane system.


## 2.2. Modelling

The following assumptions are made before the modelling:(1)
Compared with the payload’s mass, the mass of the cable connecting the payload and trolley can be ignored, and the cable is always in a tight state due to the effect of the payload. Therefore, the cable is considered massless and rigid, and the payload sways as a single pendulum.(2)
The dynamic responds of the support legs are not considered, and support points of the main beam are simplified as pin-pin supports.(3)
The cantilevers (O1−A1 and O2−A2) are considered as rigid beams because their displacements subjected to the moving trolley are small. The part of the main beam between two legs (O1−O2) is modelled as a simply supported Euler–Bernoulli beam.(4)
The cross section of the beam is constant.(5)
It is assumed that the trolley is always in contact with the top beam. The centre of the trolley has a constant distanceh from the beam top surface. The rotating movement around the centre of mass of the trolley is not considered.In this paper, based on these assumptions, the part of the main beamO1−O2 is considered as simply supported, and the general solution of its shape function is(1)ϕx=C1cosβx+C2sinβx+C3coshβx+C4sinhβx,with pined-pined boundary conditions(2)ϕ0=0,ϕ″0=0,ϕL0=0,ϕ″L0=0.Then, theith model shape function of the simply supported beam is derived as(3)ϕix=siniπxL0.The vertical displacement of the beam can be expressed as(4)yb=wx,t=∑i=1nϕixqit.The coordinates of the centre of the trolley are(5)xc=12x1+x2,yc=12wx1,t+wx2,t+h=12∑i=1nϕix1+ϕix2qit+h.The horizontal velocity of the centre of the trolley is(6)x˙c=x˙1=x˙2.And the vertical velocity is(7)y˙c=12dwx1,tdt+dwx2,tdt=12∑i=1nx˙1ϕi′x1qi+ϕix1q˙i+x˙2ϕi′x2qi+ϕix2q˙i=12∑i=1nx˙1ϕi′x1+x˙2ϕi′x2qi+12∑i=1nϕix1+ϕix2q˙i.The acceleration of trolley is known and expressed as(8)x¨c=x¨1=x¨2=a.The kinetic energy of the flexible beam, trolley, and payload are denoted byTb, Tc, and Tp, respectively, which can be calculated as follows:(9)Tb=12∫0Lρ∂yb∂t2dx=12∫0Lρ∑i=1nϕixq˙it2dx=12∑i=1n∑j=1nq˙itq˙jt∫0Lρϕixϕjxdx,(10)Tc=12mcx˙c2+y˙c2,(11)Tp=12mpx˙c+lθ˙cosθ2+y˙c−lθ˙sinθ2.Then, the total kinetic energy of the coupled system can be obtained:(12)T=Tb+Tc+Tp.The potential energy of the flexible beam, trolley, and payload are denoted byVb, Vc, and Vp, respectively, which can be calculated as follows:(13)Vb=12∫0LEI∂2yb∂x22dx=12∫0LEI∑i=1nϕi″xqit2dx=12∑i=1n∑j−1nqitqjt∫0LEIϕi″xϕj″xdx,(14)Vc=−mcgyc,(15)Vp=−mpgyc+lcosθ.Then, the total potential energy of the coupled system can be obtained:(16)V=Vb+Vc+Vp.By applying the Lagrange equation,(17)ddt∂T∂q˙i−∂T∂qi+∂V∂qi=0,i=1,2,…,n,one can derive the dynamic equation of the generalized coordinate qi:(18)14mc+mpϕix1+ϕix2∑j=1nx¨1ϕj′x1+x˙12ϕj″x1+x¨2ϕj′x2+x˙22ϕj″x2qj+14mc+mpϕix1+ϕix2∑j=1n2x˙1ϕj′x1+x˙2ϕj′x2q˙j+ϕjx1+ϕjx2q¨j+∑j=1n∫0Lρϕixϕjxdxq¨jt−12mplϕix1+ϕix2θ¨sinθ+θ˙2cosθ+∑j=1n∫0LEIϕi″xϕj″xdxqjt−12mc+mpgϕix1+ϕix2=0.By applying the Lagrange equation,(19)ddt∂T∂θ˙−∂T∂θ+∂V∂θ=0,one can derive the dynamic equation of the generalized coordinate θ.(20)−12mplsinθ·∑j=1nx¨1ϕj′x1+x˙12ϕj″x1+x¨2ϕj′x2+x˙22ϕj″x2qj−12mplsinθ·∑j=1n2x˙1ϕj′x1+x˙2ϕj′x2q˙j+ϕjx1+ϕjx2q¨j+mplx¨ccosθ+lθ¨+mpgl1sinθ=0.Equations (18) and (20) compose the dynamic equation of the coupled system.


## 3. Simulation Results and Discussion

In this paper, Newmark-β method was used to solve the numerical solution of the dynamic models (18) and (20). The integration time step was 0.001 s. The initial position of the trolley is x1=0 and x2=−d (see Figure 3), and the trolley stops at position x1=L0+d and x2=L0. All parameters are listed in Table 1.Table 1
System parameters.
ParametersCase ACase BCase CEN/m22.07×10112.15×10112.06×1011Im41.04×10−60.80.0643ρkg/m7.041.53×103599.64L0m1010050L1m0010L2m000mckg703.06×10413.165×103mpkg03.06×10436×103lm055gm/s29.89.89.8


### 3.1. Validation of the Proposed Model

To validate the proposed model, a moving mass model (without payload) [4] and a moving trolley model (the trolley is simplified as a particle/mass point) [12] were adopted.


#### 3.1.1. Results Compared with [4]

If the payload swing and the axle base of the trolley are equal to 0, i.e.,θ,d=0, equation (18) describes a simply supported beam subjected to a moving mass. For the system parameters, refer to Case A in Table 1. The velocity of the moving load is constant and its value is 3.34 m/s. The simulation results are shown in Figure 4.Figure 4
The midspan deflections of the beam compared with [4].Figure4 shows the time histories of the midspan deflections of the simply supported beam. The assumed mode method is used to solve the deflections of the beam in both [4] and this paper. As can be seen from Figure 4, the maximum deflection of the midspan is 0.0694 m, and the calculation results of the proposed model are consistent with those of [4].


#### 3.1.2. Results Compared with [12]

If the axle base of the trolley is equal to 0, i.e.,d=0, equations (18) and (20) describe a simply supported beam subjected to a moving trolley suspended with a single-pendulum payload. In this case, the trolley is simplified as a mass point. For the system parameters, refer to Case B in Table 1. The simulation results are shown in Figure 5. For Figures 5(a)–5(c), the initial velocities of the trolley are 0, while the accelerations are 1, 2, and 3 m/s2, respectively.Figure 5
The midspan deflections of the beam compared with [12].
(a)(b)(c)It can be observed that the values of the maximum midspan deflection of the beam in Figures5(a)–5(c) are 0.0751, 0.0762, and 0.0862 m, respectively. That means the maximum midspan deflection of the beam will increase as the acceleration of the trolley increases. It can be also found that the calculation results of the proposed model are consistent with those of [12].In summary, it can be easily found from Figures4 and 5 that the simulation results of this paper are quite consistent with that of [4, 12]. Therefore, the dynamic model developed in this paper is considered as an acceptable model for calculating the dynamic responses of a crane system. Besides, the trolley’s axle base is equal to 0 in [4, 12] and the payload swing angle is also equal to 0 in [4], so the proposed model can be considered as an extended version of the models in [4, 12]. In other words, the proposed model is downward compatible with those in [4, 12].


### 3.2. Effect of the Trolley’s Axle Base

In this section, we study the effect of the trolley’s axle base under zero initial swing angle condition and nonzero initial swing angle condition. At the beginning of the simulation, the initial velocity of the trolley is 0. The trolley accelerates at a constant acceleration of 0.15 m/s2. Then, it runs at a constant speed at the maximum speed. Finally, it uniformly decelerates to a standstill. The acceleration and deceleration time are both 6.67 s. For the system parameters, refer to Case C in Table 1.


#### 3.2.1. Zero Initial Swing Angle Condition

For comparison study, we choose three control groups as follows:Case C-I:d = 4 m and θ = 0 radCase C-II:d = 2 m and θ = 0 radCase C-III:d = 0 m and θ = 0 radThe simulation results are shown in Figures6 and 7. Figure 6 shows the midspan deflections with regard to time. From Figure 6, it can be easily found that the maximum midspan deflection in Case C-II (0.0946 m) is slightly smaller than that in Case C-III (0.0947 m), and the maximum midspan deflection in Case C-I (0.0940 m) is the minimum among three cases. Figure 7 shows the relationship between the vertical displacements of the centre of the trolley and the position of the centre of the trolley. Suppose that the initial displacement of the trolley is 0. Due to the different axle bases, the midspan deflections and displacements of the trolley were quite different among three cases at the initial and final phases of the simulation. It can be easily found that the larger the axle base, the smaller the displacement of the trolley.Figure 6
The midspan deflections of the beam.Figure 7
The displacement of the centre of the trolley.


#### 3.2.2. Nonzero Initial Swing Angle Condition

In previous research studies, the effect of the initial swing angle is not considered, i.e., at the beginning of the simulation, the swing angle is set to 0. In practice, a large swing angle should be avoided as much as possible because there is a risk of dangerous accidents.For comprehensive understanding of initial swing angle effects, we choose the following three nonzero initial swing angle cases for comparison study:Case C-IV:d = 4 m and θ = 0.1 radCase C-V:d = 2 m and θ = 0.1 radCase C-VI:d = 0 m and θ = 0.1 radThe simulation results are shown in Figures8 and 9. From Figure 8, it can be easily found that, regarding the maximum midspan deflection, Case C-VI has the maximum value (0.0958 m), following Case C-V (0.0954 m), and Case C-IV has the minimum value (0.0949 m). From Figure 9, one can obtain that the larger the axle base, the smaller the vertical displacement of the trolley, too.Figure 8
The midspan deflections of the beam.Figure 9
The displacement of the centre of the trolley.To further study the effect of the initial swing angle of the payload on the midspan deflection of the beam, the following two cases were added:Case C-VII:d = 4 m and θ = 0.2 radCase C-VIII:d = 4 m and θ = 0.3 radFrom Figure10, one can observe that Case C-I has the minimum value of the maximum midspan deflection of the beam (0.0940 m), following Case C-IV (0.0949 m) and Case VII (0.0977 m), and Case C-VIII has the maximum value (0.1018 m).Figure 10
The midspan deflections of the beam.Thus, as the initial swing angle increased, the maximum midspan deflection of the beam increased. As shown in Figure11, the vertical displacement of the centre of the trolley increased as the initial swing angle increased.Figure 11
The displacement of the centre of the trolley.


## 3.1. Validation of the Proposed Model

To validate the proposed model, a moving mass model (without payload) [4] and a moving trolley model (the trolley is simplified as a particle/mass point) [12] were adopted.


### 3.1.1. Results Compared with [4]

If the payload swing and the axle base of the trolley are equal to 0, i.e.,θ,d=0, equation (18) describes a simply supported beam subjected to a moving mass. For the system parameters, refer to Case A in Table 1. The velocity of the moving load is constant and its value is 3.34 m/s. The simulation results are shown in Figure 4.Figure 4
The midspan deflections of the beam compared with [4].Figure4 shows the time histories of the midspan deflections of the simply supported beam. The assumed mode method is used to solve the deflections of the beam in both [4] and this paper. As can be seen from Figure 4, the maximum deflection of the midspan is 0.0694 m, and the calculation results of the proposed model are consistent with those of [4].


### 3.1.2. Results Compared with [12]

If the axle base of the trolley is equal to 0, i.e.,d=0, equations (18) and (20) describe a simply supported beam subjected to a moving trolley suspended with a single-pendulum payload. In this case, the trolley is simplified as a mass point. For the system parameters, refer to Case B in Table 1. The simulation results are shown in Figure 5. For Figures 5(a)–5(c), the initial velocities of the trolley are 0, while the accelerations are 1, 2, and 3 m/s2, respectively.Figure 5
The midspan deflections of the beam compared with [12].
(a)(b)(c)It can be observed that the values of the maximum midspan deflection of the beam in Figures5(a)–5(c) are 0.0751, 0.0762, and 0.0862 m, respectively. That means the maximum midspan deflection of the beam will increase as the acceleration of the trolley increases. It can be also found that the calculation results of the proposed model are consistent with those of [12].In summary, it can be easily found from Figures4 and 5 that the simulation results of this paper are quite consistent with that of [4, 12]. Therefore, the dynamic model developed in this paper is considered as an acceptable model for calculating the dynamic responses of a crane system. Besides, the trolley’s axle base is equal to 0 in [4, 12] and the payload swing angle is also equal to 0 in [4], so the proposed model can be considered as an extended version of the models in [4, 12]. In other words, the proposed model is downward compatible with those in [4, 12].


## 3.1.1. Results Compared with [4]

If the payload swing and the axle base of the trolley are equal to 0, i.e.,θ,d=0, equation (18) describes a simply supported beam subjected to a moving mass. For the system parameters, refer to Case A in Table 1. The velocity of the moving load is constant and its value is 3.34 m/s. The simulation results are shown in Figure 4.Figure 4
The midspan deflections of the beam compared with [4].Figure4 shows the time histories of the midspan deflections of the simply supported beam. The assumed mode method is used to solve the deflections of the beam in both [4] and this paper. As can be seen from Figure 4, the maximum deflection of the midspan is 0.0694 m, and the calculation results of the proposed model are consistent with those of [4].


## 3.1.2. Results Compared with [12]

If the axle base of the trolley is equal to 0, i.e.,d=0, equations (18) and (20) describe a simply supported beam subjected to a moving trolley suspended with a single-pendulum payload. In this case, the trolley is simplified as a mass point. For the system parameters, refer to Case B in Table 1. The simulation results are shown in Figure 5. For Figures 5(a)–5(c), the initial velocities of the trolley are 0, while the accelerations are 1, 2, and 3 m/s2, respectively.Figure 5
The midspan deflections of the beam compared with [12].
(a)(b)(c)It can be observed that the values of the maximum midspan deflection of the beam in Figures5(a)–5(c) are 0.0751, 0.0762, and 0.0862 m, respectively. That means the maximum midspan deflection of the beam will increase as the acceleration of the trolley increases. It can be also found that the calculation results of the proposed model are consistent with those of [12].In summary, it can be easily found from Figures4 and 5 that the simulation results of this paper are quite consistent with that of [4, 12]. Therefore, the dynamic model developed in this paper is considered as an acceptable model for calculating the dynamic responses of a crane system. Besides, the trolley’s axle base is equal to 0 in [4, 12] and the payload swing angle is also equal to 0 in [4], so the proposed model can be considered as an extended version of the models in [4, 12]. In other words, the proposed model is downward compatible with those in [4, 12].


## 3.2. Effect of the Trolley’s Axle Base

In this section, we study the effect of the trolley’s axle base under zero initial swing angle condition and nonzero initial swing angle condition. At the beginning of the simulation, the initial velocity of the trolley is 0. The trolley accelerates at a constant acceleration of 0.15 m/s2. Then, it runs at a constant speed at the maximum speed. Finally, it uniformly decelerates to a standstill. The acceleration and deceleration time are both 6.67 s. For the system parameters, refer to Case C in Table 1.


### 3.2.1. Zero Initial Swing Angle Condition

For comparison study, we choose three control groups as follows:Case C-I:d = 4 m and θ = 0 radCase C-II:d = 2 m and θ = 0 radCase C-III:d = 0 m and θ = 0 radThe simulation results are shown in Figures6 and 7. Figure 6 shows the midspan deflections with regard to time. From Figure 6, it can be easily found that the maximum midspan deflection in Case C-II (0.0946 m) is slightly smaller than that in Case C-III (0.0947 m), and the maximum midspan deflection in Case C-I (0.0940 m) is the minimum among three cases. Figure 7 shows the relationship between the vertical displacements of the centre of the trolley and the position of the centre of the trolley. Suppose that the initial displacement of the trolley is 0. Due to the different axle bases, the midspan deflections and displacements of the trolley were quite different among three cases at the initial and final phases of the simulation. It can be easily found that the larger the axle base, the smaller the displacement of the trolley.Figure 6
The midspan deflections of the beam.Figure 7
The displacement of the centre of the trolley.


### 3.2.2. Nonzero Initial Swing Angle Condition

In previous research studies, the effect of the initial swing angle is not considered, i.e., at the beginning of the simulation, the swing angle is set to 0. In practice, a large swing angle should be avoided as much as possible because there is a risk of dangerous accidents.For comprehensive understanding of initial swing angle effects, we choose the following three nonzero initial swing angle cases for comparison study:Case C-IV:d = 4 m and θ = 0.1 radCase C-V:d = 2 m and θ = 0.1 radCase C-VI:d = 0 m and θ = 0.1 radThe simulation results are shown in Figures8 and 9. From Figure 8, it can be easily found that, regarding the maximum midspan deflection, Case C-VI has the maximum value (0.0958 m), following Case C-V (0.0954 m), and Case C-IV has the minimum value (0.0949 m). From Figure 9, one can obtain that the larger the axle base, the smaller the vertical displacement of the trolley, too.Figure 8
The midspan deflections of the beam.Figure 9
The displacement of the centre of the trolley.To further study the effect of the initial swing angle of the payload on the midspan deflection of the beam, the following two cases were added:Case C-VII:d = 4 m and θ = 0.2 radCase C-VIII:d = 4 m and θ = 0.3 radFrom Figure10, one can observe that Case C-I has the minimum value of the maximum midspan deflection of the beam (0.0940 m), following Case C-IV (0.0949 m) and Case VII (0.0977 m), and Case C-VIII has the maximum value (0.1018 m).Figure 10
The midspan deflections of the beam.Thus, as the initial swing angle increased, the maximum midspan deflection of the beam increased. As shown in Figure11, the vertical displacement of the centre of the trolley increased as the initial swing angle increased.Figure 11
The displacement of the centre of the trolley.


## 3.2.1. Zero Initial Swing Angle Condition

For comparison study, we choose three control groups as follows:Case C-I:d = 4 m and θ = 0 radCase C-II:d = 2 m and θ = 0 radCase C-III:d = 0 m and θ = 0 radThe simulation results are shown in Figures6 and 7. Figure 6 shows the midspan deflections with regard to time. From Figure 6, it can be easily found that the maximum midspan deflection in Case C-II (0.0946 m) is slightly smaller than that in Case C-III (0.0947 m), and the maximum midspan deflection in Case C-I (0.0940 m) is the minimum among three cases. Figure 7 shows the relationship between the vertical displacements of the centre of the trolley and the position of the centre of the trolley. Suppose that the initial displacement of the trolley is 0. Due to the different axle bases, the midspan deflections and displacements of the trolley were quite different among three cases at the initial and final phases of the simulation. It can be easily found that the larger the axle base, the smaller the displacement of the trolley.Figure 6
The midspan deflections of the beam.Figure 7
The displacement of the centre of the trolley.


## 3.2.2. Nonzero Initial Swing Angle Condition

In previous research studies, the effect of the initial swing angle is not considered, i.e., at the beginning of the simulation, the swing angle is set to 0. In practice, a large swing angle should be avoided as much as possible because there is a risk of dangerous accidents.For comprehensive understanding of initial swing angle effects, we choose the following three nonzero initial swing angle cases for comparison study:Case C-IV:d = 4 m and θ = 0.1 radCase C-V:d = 2 m and θ = 0.1 radCase C-VI:d = 0 m and θ = 0.1 radThe simulation results are shown in Figures8 and 9. From Figure 8, it can be easily found that, regarding the maximum midspan deflection, Case C-VI has the maximum value (0.0958 m), following Case C-V (0.0954 m), and Case C-IV has the minimum value (0.0949 m). From Figure 9, one can obtain that the larger the axle base, the smaller the vertical displacement of the trolley, too.Figure 8
The midspan deflections of the beam.Figure 9
The displacement of the centre of the trolley.To further study the effect of the initial swing angle of the payload on the midspan deflection of the beam, the following two cases were added:Case C-VII:d = 4 m and θ = 0.2 radCase C-VIII:d = 4 m and θ = 0.3 radFrom Figure10, one can observe that Case C-I has the minimum value of the maximum midspan deflection of the beam (0.0940 m), following Case C-IV (0.0949 m) and Case VII (0.0977 m), and Case C-VIII has the maximum value (0.1018 m).Figure 10
The midspan deflections of the beam.Thus, as the initial swing angle increased, the maximum midspan deflection of the beam increased. As shown in Figure11, the vertical displacement of the centre of the trolley increased as the initial swing angle increased.Figure 11
The displacement of the centre of the trolley.


## 4. Conclusions

A novel moving load model for 2D crane systems of which the trolley has two axles was proposed in this paper. By comparing with existing papers, the proposed model was verified and can be considered as an extended version of comparative models. After the validation, the dynamics of a gantry crane system with a two-axle trolley were studied. On the basis of the forgoing analyses, one can draw the conclusions:(1)
There is small amplitude vibration of the beam under the action of the two-axle trolley. Thus, the dynamic responses of the beam subjected to a two-axle trolley is more accurate than that to a single-axle trolley.(2)
If the system parameters keep unchanged, the larger the axle base of the trolley, the smaller the maximum midspan deflection. Therefore, with a large total mass of the trolley and payload, increasing the length of trolley’s axle base will decrease the maximum midspan deflection.(3)
The greater swing angle coincides with the greater midspan deflection. In order to avoid impact on the flexible beam, the swing angle should be controlled as small as possible.



---

*Source: 3096213-2020-07-09.xml*