Abstract
The vibration of gun barrels would result in the change of impact point, which would further reduce the firing accuracy of weapons. In the past, the calculation model based on the Euler beam theory could not satisfy the accuracy requirements. Based on the Timoshenko beam theory, the vibration equation of the stepped beam is established by invoking continuum transfer matrix method. The forced vibration of the stepped beam under the inertial moving load is solved. The model has better precision than the Euler beam model. The endpoint of the cantilever beam is analyzed. It is shown that the endpoint response increases with the increasing mass and acceleration of moving load, so does the inertial coefficient. With the increase of moving load speed, the endpoint response decreases, and the inertia coefficient increases. Among the three parameters, the mass of moving load is the main factor affecting the inertia coefficient. Furthermore, both free and forced vibrations of other stepped beam shaped structures with arbitrary segments and boundary conditions can be explored by using the proposed method.
1. Introduction
In the process of artillery launching, the gravity and eccentricity which is caused by high speed projectile will induce barrel vibration. The coupling vibration between projectile and barrel is an important factor which affects the firing accuracy of artillery [1]. With the development of new artillery, it is necessary to study the barrel vibration caused by the motion of projectile.
Zhou Ding [2] and others obtained a general solution for the lateral vibration of the barrel during single and continuous firing, and established a complete theory and analysis method for the vibration of the barrel caused by the movement in the projectile bore. Peng Xian [3] established the dynamic equations of the artillery system under the motion of the base, and solved them separately using the perturbation method and the numerical method, and proposed an analytical expression for calculating the deflection of the barrel end. Wei Shechun [4] and others fully considered the influence of the barrel on the movement and force of the projectile during the movement of the projectile, used VC ++ to carry out the parametric secondary development of ANSYS / LSDYNA, and established a projectile launch finite element integration platform to Precision and efficiency of projectile launch dynamic analysis.
Kang Xinzhong [5] treated the barrel as a conical cantilever beam, and derived the finite element equation of motion of the beam element based on the Dalambert principle and the deformed body virtual work principle, and gave the overall mass matrix, damping matrix, stiffness matrix and the set rule of nodal force array. Liu Ning [6] simplified the barrel into a cantilever beam of equal crosssection. The vibration equation of the barrel was established according to the BernoulliEuler elementary beam theory, and the vibration equation of the barrel was solved by modal analysis. The built model basically describes the dynamic characteristics of the system, and provides a simple and effective method for studying the coupling problem of projectiles. Su Zhongting [7] applied the finite element method and cantilever beam theory based on the idea of moving load, simplified the barrel into a flexible variablesection cantilever beam, and applied the cylindrical coordinate system to discretize the barrel structure into standard fanshaped elements. Timevarying dynamic equations of the radial vibration of the barrel under the action of moving projectiles. Ma Jisheng [8] treated the barrel as a cantilever beam on a kinematic support and the projectile as a rigid body. Based on the Kane equation Huston method, a relatively complete dynamic model for studying the coupling problem between projectile and barrel was established. Using this model, the projectile launch process can be analyzed more accurately.
As shown in Fig. 1, the gun barrel can be simplified to a stepped beam. Many researchers have worked on vibration of stepped beams. Mao [9] studied the free vibration of a stepped EulerBernoulli beam composed of two uniform sections using the Adomian decomposition method (ADM). Each section is considered a substructure which can be modeled using ADM. Based on the variation method, Su [10] presented an effective formulation for vibration analysis of multiplestepped functionally graded beams with general boundary conditions. Giunta [11] and Cicirello presented an approach to analyze jointed Euler beams with step changes in material and crosssection under static and dynamic loads. ElSayed [12] and Farghaly used the normalized transfer matrix (NTM) to derive the exact solution of stepped Timoshenko beam. Yokoyamat [13] considers the effect of axial force, transverse shear deformation and rotational inertia on the Timoshenko beam of equal crosssection. The free vibration equation of Timoshenko beam of equal crosssection is derived from the Hamilton principle, and the natural frequency and mode of the Timoshenko beam are obtained. EISayed [14] solves the natural frequency and mode of a singlespan Timoshenko beam of equal crosssection under axial force and proposes a mathematical model for structural engineering. Hu [15] and Wang formed the Hamilton system for solving the dynamic characteristics of the cantilever Timoshenko beam on the basis of considering the effects of shear deformation and rotational inertia.
Fig. 1Simplified diagram of gun barrel
The dynamic problem of beams under moving loads is also a research hot spot. Yang and Wang [16] studied the dynamic response and stability of an inclined beam under a moving vertical concentrated load. They presented the governing equation for transverse motion of the beam. It contained the effect of the compressive axial component of the vertical load on the bending stiffness capacity of the beam. Sudheesh [17, 18] analysed the forcedfree responses of nonuniform beams under moving point loads. They presented simple approximate analytical formula for the forced responses of undamped nonuniform beams, which was derived using the fundamental mode by the RayleighRitz (RR) method. Saif and David [19] presented a rigorous study on the behaviour of buried pipes under static and moving traffic loads using a robust finite element analysis. Olga [20] analysed the dynamic behaviour of a Rayleigh multispan uniform continuous beam system that is traversed by a constant moving force. You and Yan [21] concerned with the mechanical responses of anisotropic multilayered medium under harmonic moving load. Fang [22] established a new trackmultilayer ground model to investigate railway subgrade dynamic responses induced by moving train load.
Most articles on gun barrel vibration is based on Euler beam theory. This method has a small scope of application which is only suitable for slender structures [23]. Based on Timoshenko beam theory, the transfer matrix of Timoshenko beam is derived in this paper. The equation of inertia moving load is obtained by analysing the motion of the moving load. The forced vibration of cantilever beam is solved by Newmark method. The influence of mass, speed and acceleration on inertia coefficient is analysed in detail. The modelling and analysis methods adopted in this paper are also applicable to other engineering vibration problems caused by moving loads.
2. The vibration equation of stepped cantilever beam
The primary suppose based on the Euler beam theory is simple. The calculated results are close to the true values when the beam is slender or the loworder frequency is analysed. The Timoshenko beam theory considers the shear deformation and the moment of inertia caused by the bending deformation. It guarantees the accuracy of dynamic parameters of beams when analyse highorder mode or deep beam model. Fig. 2 shows the force and deformation state of the beam microunit.
Fig. 2Timoshenko beam microunit force and deformation diagram
In the Fig. 2, ${Q}_{y}$ is the shear, ${M}_{z}$ is the moment, ${I}_{q}$ is the inertial force acting on the beam microunit, ${I}_{0}$ is the inertia moment acting on the beam microunit, $\gamma $ is the shear angle produced by shear deformation on the neutral axis, $\theta $ is the angle produced by bending, $y$ is the transverse deformation of the beam microunit. The inertial force ${I}_{q}$ and the inertia moment ${I}_{0}$ acting on the beam microunit are as shown in the following formulas:
where $\mu $ is the linear density, $I$ is the moment of inertia, and $\rho $ is the density. According to the theoretical mechanics, the dynamic equilibrium equations of beam microunit are listed:
According to the material mechanics:
Bring Eq. (3a) and Eq. (3b) into Eq. (2a) and Eq. (2b), then obtain the transverse vibration equation of the Timoshenko beam:
where $E$ is the elastic modulus, $G$ is the shear modulus, $A$ is the crosssectional area, $k$ is the correction coefficient considering the uneven distribution of shear strain on the section.
Separating $y$ and $\theta $ in Eq. (4a) and Eq. (4b) to obtain a new system of equations:
Supposing the system performs the samefrequency harmonic motion $y(x,t)=y\left(x\right)\mathrm{s}\mathrm{i}\mathrm{n}\omega t$. Where $\omega $ is the free vibration angular frequency of the system, $y\left(x\right)$ is mode shape function. Bring $y\left(x,t\right)$ into Eq. (5a):
where $g$ and $a$ are intermediate variables:
Supposing the general solution of $y\left(x\right)$:
Bringing Eq. (8) into Eq. (6):
Solving the roots of Eq. (9) are $\pm {\stackrel{}{\lambda}}_{1}$ and $\pm i{\stackrel{}{\lambda}}_{2}$:
Bringing the obtained root into Eq. (8) and get the general solution of $y\left(x\right)$ as the following equation:
where ${C}_{1}$ to ${C}_{4}$ are undetermined constants. The fourthorder partial differential equation of $\theta \left(x,t\right)$ is the same as $y\left(x,t\right)$. Therefore, the solution of the fourthorder partial differential equation of $\theta \left(x,t\right)$ is as follow:
where ${B}_{1}$ to ${B}_{4}$ are undetermined constants. Eq. (11) and Eq. (12) are not independent. Bringing Eq. (11) and Eq. (12) into Eq. (4b), the relationship between ${B}_{1}$~${B}_{4}$ and ${C}_{1}$~${C}_{4}$ can be obtained by using the condition that the coefficients of $\mathrm{s}\mathrm{i}\mathrm{n}{\stackrel{}{\lambda}}_{2}x$, $\mathrm{c}\mathrm{o}\mathrm{s}{\stackrel{}{\lambda}}_{2}x$, $\mathrm{s}\mathrm{i}\mathrm{n}\mathrm{h}{\stackrel{}{\lambda}}_{1}x$ and $\mathrm{c}\mathrm{o}\mathrm{s}\mathrm{h}{\stackrel{}{\lambda}}_{1}x$ on both sides of the equation are equal:
Assuming intermediate variables $\alpha $ and $\beta $:
Bringing Eq. (11) and Eq. (12) into Eq. (3a) and Eq. (3b), then can get the following equations:
Combining Eq. (11), (12), (15a), (15b) to establish equations:
Listing the boundary conditions at the left end of the beam:
Bringing the boundary into the Eq. (16):
By introducing Eq. (13a) to Eq. (13d) into the above equation. ${C}_{1}$ to ${C}_{4}$ can be obtained:
By taking the obtained ${C}_{1}$~${C}_{4}$ into the Eq. (16), the entire equations for $y\left(x\right)$, $\theta \left(x\right)$, $Q\left(x\right)$, $M\left(x\right)$ can be obtained. Bringing the boundary conditions at the right end of the beam:
By introducing Eq. (20) into Eq. (16), the transverse vibration transfer matrix of the continuous Timoshenko beam can be obtained:
where [$T$] is the transfer matrix:
Assuming intermediate variables $\sigma $, $\tau $, ${\beta}^{4}$ and ${c}_{0}$${c}_{5}$:
This gives the transfer matrix [$T$] of a single beam. The transfer matrix of stepped beam can be obtained by multiplying multiple transfer matrices:
The transfer matrix of the step beam is represented by [$H$]:
Bringing cantilever beam boundary conditions:
The transfer matrix of cantilever beam is shown:
The coefficient determinant must be zero as the equation having nonzero solution:
Solving the determinant, the natural frequency $\omega $ of the stepped beam can be obtained. The mode shape functions are obtained by bringing the natural frequency into the Eq. (11) and Eq. (12). Then normalize the mode function, the normal modes $\varphi \left(x\right)$ and $\phi \left(x\right)$ can be obtained. Where $\varphi \left(x\right)$ is transverse displacement mode shape function and $\phi \left(x\right)$ is bending angle mode shape function.
Fig. 3Moving load acceleration analysis diagram
As shown in the Fig. 3, since the beam itself is vibrating, the moving load has a transverse inertial force. When the barrel coupling is not considered, the moving load is only gravity mg. When considering the barrel coupling, the inertia term of the moving load needs to be increased. In the Fig. 3, θ is the bending angle; $v$ is the velocity of the moving load; ${a}_{y}$ is the transverse acceleration of the microunit due to the bending deformation of the beam; ${a}_{\tau}$ is the tangential acceleration of the moving load; ${a}_{n}$ is the normal acceleration of the moving load; ${a}_{k}$ is the Coriolis acceleration caused by the interaction between rotational angular velocity and axial motion velocity.
According to the knowledge of material mechanics and theoretical mechanics:
The actual vertical acceleration of moving load is as follows:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\frac{{\partial}^{2}y}{\partial {t}^{2}}+\frac{dv}{dt}\theta +{v}^{2}\frac{{\partial}^{2}y}{\partial {x}^{2}}+2v\frac{\partial \theta}{\partial t}.$
The force acting on the beam by the moving load considering the inertial effect is $P\left(x,t\right)$:
where, ${m}_{y}$ represents the mass of the moving load; $s\left(t\right)$ represents the displacement of the moving load. In the Eq. (33), the first term represents the gravitational acceleration. The second term represents the transverse acceleration. The third term represents the Coriolis acceleration. The fourth and fifth terms represent the normal acceleration component and the tangential acceleration component, respectively.
$\delta (xs(t\left)\right)$ is the Dirac Delta function:
Bring $y(x,t)=\varphi \left(x\right)q\left(t\right)$ and $\theta (x,t)=\varphi \left(x\right)q\left(t\right)$ into the Eq. (33):
In Eq. (35), $i$ represents the $i$th order national frequency. By using the method of separating variables and the orthogonality of modal functions, the integral of the whole beam can be obtained:
By separating and sorting out the terms of zero, primary and secondary derivatives of generalized coordinates in the formula, the [$M$], [$K$], [$C$] matrices of the system are obtained, where [$M$], [$K$], [$C$] are symmetric matrices. The vibration differential equation of the whole system is expressed as follows:
In the above Eq. (37), ${\varphi}_{i}\left(s\right)$ and ${\phi}_{i}\left(s\right)$ represent the mode function of the beam where the moving load displacement is $s$. For such timevarying coefficient differential equations, it can only be solved by numerical method of stepbystep integration. Newmark$\beta $ method is the most commonly used integration method. By calculating the generalized coordinate $q\left(t\right)$ at any time, the dynamic response $y\left(x,t\right)$ of the stepped beam at any time can be obtained.
3. Example analysis and discussion
3.1. Numerical validation
As shown in Fig. 4, a fivestepped beam with a moving load is analysed. In this section, the Timoshenko beam transfer matrix is compared with the Euler beam transfer matrix and FEM calculation results to verify the correctness of the model.
The parameters of the cantilever beam are given in Table 1.
Numerical solutions of Timoshenko beam model and Euler beam model are calculated by programming in MATLAB software. The APDL language of ANSYS software is used in the finite element method for parametric programming. The stepped beam uses the BEAM188 element. The moving load uses the birthdeath element. The finite element model is shown in Fig. 5.
The first four natural frequencies are extracted for comparison. The results are shown in the following Table 2.
Fig. 4Fivestepped cantilever beam model diagram
Fig. 5Finite element simulation diagram
Table 1Model parameter table
$E=$2.1×10^{5 }MPa, ${\mu}_{P}=$ 0.3, $\rho =$7800 kg/m^{3}, $k=$3/4  
Section number  Length (m)  Inner diameter (m)  Outer diameter(m) 
${L}_{1}$  1.0  0.20  0.4 
${L}_{2}$  1.0  0.18  0.35 
${L}_{3}$  1.0  0.15  0.30 
${L}_{4}$  0.5  0.12  0.25 
${L}_{5}$  0.5  0.10  0.20 
Table 2Frequency comparison table
FEM / Hz  Timoshenko beam / Hz  Euler beam / Hz  Timoshenko deviation (%)  Euler deviation (%)  
${\omega}_{1}$  28.45  28.54  28.78  0.32  1.16 
${\omega}_{2}$  118.03  119.13  123.68  0.93  4.79 
${\omega}_{3}$  272.17  276.88  299.63  1.73  10.09 
${\omega}_{4}$  467.72  479.57  544.75  2.53  16.47 
It is shown in Table 2 that the difference of the first natural frequency is small. For the following three natural frequencies, the difference between Euler beam and FEM is big. As mentioned in the introduction, the poor calculation accuracy of Euler beams is because each simplified stepped beam belongs to a deep beam structure, and does not consider the section shear deformation and moment of inertia. Therefore, Timoshenko beam has higher accuracy than Euler beam.
The mass of moving load is 10 kg, the speed is 100 m/s. It runs in the beam at a constant speed. The forced vibration of beam is analyzed by Newmark$\beta $ method. The step length is set to 2×10^{4} s. The Newmark coefficient ${\gamma}_{N}$ is 0.5, ${\beta}_{N}$ is 0.25. The response of the free endpoint is extracted as the research object. The calculation results are shown in the following figures.
Fig. 6 shows the time domain response of the free endpoint on the beam. With the movement of the moving load, the displacement of the free endpoint increases. The FEM ’s response is the biggest. The Euler beam’s response is the smallest. Fig. 7 shows the percentage of the difference between the Timoshenko beam model and the Euler beam model for FEM results. Through the error comparison, we can know that Timoshenko beam is more accurate than Euler beam.
The 4 m and 8 m length beams are selected for studying the influence of the beam size on natural frequencies. The dimensions are shown in Table 3.
Fig. 6Response contrast diagram
Fig. 7Error percentage contrast diagram
Table 3Beam size parameter table
Beam length  Segment length (m)  Size  Inner and outer diameter parameters (m) 
4 m  ${L}_{1}=$1, ${L}_{2}=$ 1, ${L}_{3}=$1, ${L}_{4}=$ 0.5, ${L}_{5}=$ 0.5  Thick  ${D}_{1}=$0.8, ${D}_{2}=$0.7, ${D}_{3}=$0.6, ${D}_{4}=$0.5, ${D}_{5}=$0.4, ${d}_{1}=$ 0.4, ${d}_{2}=$0.36, ${d}_{3}=$ 0.3, ${d}_{4}=$ 0.24, ${d}_{5}=$0.2 
Middling  ${D}_{1}=$0.4, ${D}_{2}=$0.35, ${D}_{3}=$0.3, ${D}_{4}=$0.25, ${D}_{5}=$0.2 ${d}_{1}=$0.2, ${d}_{2}=$0.18, ${d}_{3}=$ 0.15, ${d}_{4}=$ 0.12, ${d}_{5}=$0.1  
8 m  ${L}_{1}=$ 2, ${L}_{2}=$2, ${L}_{3}=$2, ${L}_{4}=$1, ${L}_{5}=$1  
Slender  ${D}_{1}=$0.2, ${D}_{2}=$0.175, ${D}_{3}=$ 0.15, ${D}_{4}=$ 0.125, ${D}_{5}=$0.1 ${d}_{1}=$0.1, ${d}_{2}=$0.09, ${d}_{3}=$0.075, ${d}_{4}=$0.06, ${d}_{5}=$0.05 
The first four natural frequencies are analyzed. The result is shown in Table 4. It can be seen that if the length of the beam is fixed, the smaller the diameter, the smaller the deviation; if the diameter is fixed, the longer the length of the model, the smaller the deviation. The slender the size of the beam, the higher the accuracy of the calculation. The accuracy of the Timoshenko beam model is 35 times higher than Euler beam which is greater for higher order natural frequencies. For the slender beam, both theoretical models can solve the problem well. However, if it is a short thick beam or under the highorder natural frequency, the Euler beam model is no longer applicable. For a 4m thick beam, the maximum lengthtodiameter ratio of a single segment reached 0.8, the error of the fourthorder frequency of Timoshenko beam is more than 5 %. Therefore, in actual calculation, the diameter of the beam should not be simplified too thick, which makes the calculation accuracy worse.
3.2. Study the effect of moving load single variable on beam vibration
The stepped cantilever beam model in Table 1 is selected for analysis. Moving loads can be divided into consideringinertiaeffect and consideringgravityonly. The response of the free endpoint is extracted for analysis. The inertia coefficient ${\lambda}_{I}$_{}is defined:
In the Eq. (38), ${y}_{I}$ represents the response of the free endpoint when the moving load considering inertia effect moving to the free endpoint; ${y}_{g}$ is the response of the free endpoint when the moving load only considering gravity moving to the free endpoint.
3.2.1. The impact of moving load mass on response results
The influence of mass on the response of the cantilever beam is studied by changing the mass of moving load. The moving load moves uniformly on the beam at the speed of 100 m/s. The mass of the moving load ranges from 5 kg to 500 kg, which the increment is 5 kg. When the moving load moves to the free endpoint, the response of the free endpoint is selected as the research object. The curves obtained are as follows.
Table 4Size effect on frequency table
Beam length  Size  FEM (Hz)  Timoshenko beam (Hz)  Euler beam (Hz)  Timoshenko deviation (%)  Euler deviation (%) 
4 m  Thick  55.16  55.62  57.57  0.83  4.37 
210.2  215.94  247.36  2.73  17.68  
444.56  464.18  599.25  4.41  34.80  
709.51  749.94  1089.49  5.70  53.56  
Middling  28.45  28.54  28.78  0.32  1.16  
118.03  119.13  123.68  0.93  4.79  
272.17  276.88  299.63  1.73  10.09  
467.72  479.57  544.75  2.53  16.47  
Slender  14.34  14.36  14.39  0.14  0.35  
61.07  61.24  61.84  0.28  1.26  
145.9  146.65  149.81  0.51  2.68  
260.73  262.77  272.37  0.78  4.46  
8 m  Thick  14.23  14.26  14.39  0.21  1.12 
59.02  59.54  61.84  0.88  4.78  
136.09  138.37  149.81  1.68  10.08  
233.86  239.66  272.37  2.48  16.47  
Middling  7.17  7.18  7.2  0.14  0.42  
30.53  30.61  30.92  0.26  1.28  
72.95  73.33  74.91  0.52  2.69  
130.36  131.39  136.19  0.79  4.47  
Slender  3.59  3.59  3.6  0.00  0.28  
15.4  15.42  15.46  0.13  0.39  
37.19  37.25  37.45  0.16  0.70  
67.31  67.46  68.09  0.22  1.16 
It can be seen from Fig. 8, the forced vibration response of the free endpoint increases with the increase of moving load mass. The vibration response of the load with inertia effect is smaller than the response only considering gravity, which accords with the result of Eq. (35). As can be seen in Fig. 9, the inertia coefficient increases with the increase of mass. It means that if the mass of moving load is large, the effect of inertia must be considered.
Fig. 8Massresponse contrast diagram
Fig. 9Massinertia coefficient diagram
3.2.2. The impact of moving load speed on response results
The influence of speed on the response of the cantilever beam is studied by changing the speed of moving load. The moving load moves uniformly on the beam. The mass of the moving load is 10 kg. The speed of the moving load ranges from 150 m/s to 750 m/s, which the increment is 5 m/s. When the moving load moves to the free endpoint, the response of the free endpoint is selected as the research object. The curves obtained are as follows.
It can be seen from Fig. 10, the forced vibration response of the free endpoint decreases gradually with the increase of speed. The vibration response of the load with inertia effect is smaller than the response only considering gravity. It can be known from Eq. (33) that as the speed increases, the amplitude of the excitation force decreases, so the forced response at the free end decreases. As can be seen in Fig. 11, the inertia coefficient increases with the increase of velocity. This means that if the speed of moving load is fast, the effect of inertia must be considered.
Fig. 10Speedresponse contrast diagram
Fig. 11Speedinertia coefficient diagram
3.2.3. The impact of moving load acceleration on response results
The influence of acceleration on the response of the cantilever beam is studied by changing the acceleration of moving load. The mass of the moving load is 10 kg. The speed is 100 m/s. The acceleration of the moving load ranges from 0 m/s^{2} to 2000 m/s^{2}, which the increment is 10 m/s^{2}. When the moving load moves to the free endpoint, the response of the free endpoint is selected as the research object. The curves obtained are as follows.
Fig. 12Accelerationresponse contrast diagram
Fig. 13Accelerationinertia coefficient diagram
It can be seen from Fig. 12, the forced vibration response of the free endpoint increases with the increase of moving load acceleration. The vibration response of load with inertia effect is smaller than that without inertia effect. As can be seen from Fig. 13, with the increase of the acceleration, the inertia effect coefficient increases accordingly, which means the effect of inertia must be considered if the acceleration of moving load is large. In the range of 01000 m/s^{2}, the slope of inertia coefficient curve is close to horizontal. At this time, the acceleration has little influence on the inertia effect. When the acceleration exceeds 1000 m/s^{2}, the slope of the curve suddenly increases. The acceleration has a great influence on the inertia effect.
3.3. Comparison of double variables of moving load
3.3.1. The impact of the mass and speed of moving load on response results
Assume that the moving load performs a uniform motion on the beam. Change the mass and the speed of the moving load. Record the inertia coefficient when the moving load moves to the free endpoint. The results are shown in Fig. 14.
It can be seen from Fig. 14 that both mass and speed affect the value of inertia coefficient. With the mass increases, the value of inertia coefficient increases accordingly. And the same is true of speed to the value of inertia coefficient. The change of mass can significantly change the inertia coefficient. It can be seen from the figure that when the speed is 150 m/s, the inertia coefficient of 100 kg mass is more than twice that of 50 kg mass.
Also, in Fig. 14, when the mass of moving load is large, the influence of inertial effect on vibration response cannot be neglected regardless of the speed. When the mass of moving load is small, the effect of inertia on vibration response at low speed is negligible. However, when the moving load whose mass is 50 kg and speed is 120 m/s, the inertia coefficient exceed 5 %. The inertia effect to the vibration response influence can’t be neglected in this situation.
Fig. 14Inertial coefficient of moving load with different mass and speed
3.3.2. The impact of the mass and acceleration of moving load on response results
The acceleration analyzed in this paper is uniform acceleration. The influence of inertia effect on vibration response is considered when the moving load performs uniform acceleration motion. Consider the uniform acceleration movement with an initial velocity of 100 m/s. Record the inertia coefficient of the free endpoint. Take the acceleration and mass as variables. Collect multiple sets of data for analysis. The results of uniform acceleration motion are shown in Fig. 15.
As can be seen from Fig. 15, when the mass of the moving load is small, the curve is close to horizontal. At this time, the effect of acceleration on inertia coefficient is negligible. When the mass of the moving load is 10 kg and 20 kg, the inertia coefficient does not exceed 2 %. However, when the mass of the moving load is large, the slope of the curve increases observably. At this time, the influence of acceleration on inertia coefficient cannot be ignored. The inertia coefficient increases consistently with the increase of the speed and the mass. Acceleration has little effect on inertia coefficient in this example. It is only effective at changing the inertia influence coefficient when the mass is large.
Fig. 15Inertial coefficient of moving load with different mass and acceleration
4. Conclusions
Based on the stepped Timoshenko beam model, the transverse vibration response of the cantilever beam under inertial moving loads is analyzed. The conclusions are followed below.
1) Based on Timoshenko beam theory, the transverse vibration transfer matrix is derived in this paper. The equation of the inertia moving load is obtained by analyzing the motion of the moving load. The accuracy of this method is verified by comparing with the Euler beam model and the FEM model. This method can provide theoretical basis for gun barrel vibration.
2) The mass, velocity and acceleration of the moving load all affect the inertia coefficient. The inertia coefficient increases with the increasing mass, velocity and acceleration. This also means that the inertia effect of the moving load must be considered when these three parameters exceed the setting limits. The mass of the moving load is a major factor affecting the inertia coefficient relative to speed and acceleration.
Acknowledgements
The research work is supported by Marine LowSpeed Engine ProjectPhase I.
References

Chen M. M. Projectile balloting attributable to gun tube curvature. Shock and Vibration, Vol. 17, Issue 1, 2010, p. 3953.

Zhou D., Xie Y. S. Small parameter solution of barrel vibration caused by motion in projectile chamber. Shock and Vibration, Vol. 1, Issue 1, 1999, p. 7883.

Peng X., Zhang P., Liu Z. J. Vibration analysis of the artillery system under the movement of the base. Shock and Vibration, Vol. 28, Issue 3, 2009, p. 2326+58+196.

Wei S. C., Bai Z., Yang L. Establishment of a finite element simulation platform for projectile launch dynamics of full gun model based on ANSYS/LSDYNA. Journal of Bullets and Guidance, Vol. 34, Issue 5, 2014, p. 148150+155.

Kang X. Z., Wang B. Y. Finite element analysis of barrel vibration. Ordnance Journal, Vol. 1, Issue 1, 1990, p. 1622.

Liu N., Yang G. L. Effect of lateral impact of bullet tube on barrel dynamic response. Journal of Ballistics, Vol. 22, Issue 2, 2010, p. 6770.

Su Z. T., Xu D., Li X. L., Han Z. F. Finite element timehistory analysis of smallcaliber artillery projectile coupling dynamic response. Shock and Vibration, Vol. 31, Issue 1, 2012, p. 104108+147.

Ma J. S., Wang R. L. Theoretical model of the artillerycoupling problem. Ordnance Journal, Vol. 1, Issue 1, 2004, p. 7377.

Mao Q. B., Pietrzko S. Free vibration analysis of multiplestepped beams by using Adomian decomposition method. Mathematical and Computer Modelling of Dynamical Systems, Vol. 54, Issue 1, 2011, p. 756764.

Su Z., Jin G. Y., Ye T. G. Vibration analysis of multiplestepped functionally graded beams with general boundary conditions. Composite Structures, Vol. 186, Issue 3, 2018, p. 315323.

Giunta F., Cicirello A. On the analysis of jointed EulerBernoulli beams with step changes in material and crosssection under static and dynamic loads. Engineering Structures, Vol. 179, Issue 1, 2019, p. 6678.

ElSayed T.A., Farghaly S. H. A normalized transfer matrix method for the free vibration of stepped beams: comparison with experimental and FE(3D) methods. Shock and Vibration, Vol. 46, Issue 1, 2017, p. 123.

Yokoyama T. Vibration analysis of Timoshenko beamcolumns on twoparameter elastic foundations. Computers and Structures, Vol. 61, Issue 1, 1996, p. 9951007.

EiSayed T., Farghaly S. H. Exact vibration of Timoshenko beam combined with multiple mass spring subsystems. Structural Engineering and Mechanics, Vol. 57, Issue 1, 2016, p. 9891014.

Hu Q. P., Wang Q. Q. A new method for analyzing the natural vibration characteristics of Timoshenko cantilever beam. Journal of Hebei University of Engineering, Vol. 30, Issue 1, 2013, p. 13.

Yang D. S., Wang C. M. Dynamic response and stability of an inclined Euler beam under a moving vertical concentrated load. Engineering Structures, Vol. 186, Issue 3, 2019, p. 243254.

Sudheesh K. C. P., Sujatha C., Shankar K. Vibration of nonuniform beams under moving point loads: an approximate analytical solution in time domain. International Journal of Structural Stability and Dynamics, Vol. 17, Issue 3, 2016, p. 175186.

Sudheesh K. C. P., Sujatha C., Shankar K. Nonuniform EulerBernoulli beams under a single moving oscillator: an approximate analytical solution in time domain. Journal of Mechanical Science and Technology, Vol. 30, Issue 10, 2016, p. 44794487.

Saif A., David C., Faramarzi A. A comparative study of the response of buried pipes under static and moving loads. Transportation Geotechnics, Vol. 15, Issue 2, 2018, p. 3946.

Olga S. B., Paweł S., Filip Z. Application of Volterra integral equations in the dynamics of a multispan Rayleigh beam subjected to a moving load. Mechanical Systems and Signal Processing, Vol. 121, Issue 2, 2019, p. 777790.

You L. Y., Yan K. Z., Liu N. Y., Shi T. W., Lv S. T. Assessing the mechanical responses for anisotropic multilayered medium under harmonic moving load by Spectral Element Method (SEM). Applied Mathematical Modelling, Vol. 67, Issue 3, 2019, p. 2237.

Fang R., Lu Z., Yao H. L., Luo X. W., Yang M. L. Study on dynamic responses of unsaturated railway subgrade subjected to moving train load. Soil Dynamics and Earthquake Engineering, Vol. 115, Issue 6, 2018, p. 319323.

Li H. L., Zhang B., Gong Y. X., Wang D. H. Response analysis of ladder beams under inertial moving load. International Journal of Acoustics and Vibration, Vol. 23, Issue 3, 2018, p. 402410.