A Modified FEM for Transverse and Lateral Vibration Analysis of Thin Beams Under a Mass Moving with a Variable Acceleration

In this paper, a new modified finite element method that can be used in the analysis of transverse and lateral vibrations of the thin beams under a point mass moving with a variable acceleration and constant jerk is presented. Jerk is the change in acceleration over time. In this method, the classical finite element of the beam is modified by the inclusion of the inertial effects of the moving mass. This modification is made using the relations between nodal forces and nodal deflections and shape functions of six DOF beam element. The mass, stiffness, and damping matrices of the modified finite element are determined by forces caused by the corresponding transverse and lateral accelerations and jerks, and transverse Coriolis and centripetal accelerations and jerks, respectively. This method was first applied on a simply supported beam plate to provide a comparison with the previous studies in literature, and it was proved that the results were within acceptable limits. Secondly, it was applied on a CNC type box-framed beam to analyse the dynamic response of the beam in terms of variable acceleration and jerk as well as constant velocity and mass ratios.


INTRODUCTION
Dynamic behaviour of mechanical and structural systems with moving masses is a very important engineering problem.In the design of the running systems under very high-speed masses and with a very little working tolerance such as high speed precision metal processing, high-speed rail transportation, robotics applications in microchip production, high-speed projectile firing from a barrel, highspeed rocket launching systems, and etc., some new engineering problems arise together with speed increase.In such systems, the pre-engineering analysis must be accurate in order to prevent possible accidents and failures in the application.For example, the design and development studies of a missile and high-speed rail systems that cost millions of dollars should not only depend on experimental Latin American Journal of Solids and Structures 14 (2017) 485-511 studies.When an accident is happened, no one retrieve the situation in an experimental trial.In this case, realistic modelling and analysis is vital for such systems.The computation methods given in literature are commonly for low speeds and simple cases where the motion at variable speeds, acceleration and inertia effect of mass are neglected.In this study, a more realistic method of modelling and analysis of such systems is provided without neglecting the damping, inertia effects of the mass in motion, effect of a variable speed, effect of a variable acceleration, and effect of a constant jerk.Due to the variety of its applications, the moving load problem has been widely studied in literature by many researchers.For example, (Gerdemeli et al. 2011, Kahya 2012, Sharbati and Szyszkowski 2011, and Esen 2013) have studied the subject for a constant velocity of a moving mass.An analytical solution of the effects of moving load motion for Timoshenko beams has been provided by Lee 1996.Considering the effect of a variable velocity load on structures (Michaltsos 2002, Wang 2009, Esen 2011, and Dyniewicz and Bajer 2012 have studied the dynamic behaviour of beams under accelerated masses.Inertia effects of the moving mass continue to be a point of interest for bridge dynamics, railroad design, and other high-velocity delicate motion processes and some others (Michaltsos et al. 1996, Michaltsos and Kounadis 2001, and Dehestani et al. 2009) have studied and proposed some simplified analytical solutions of the inertial effects of moving masses.Wu 2005 has studied the vertical and horizontal displacements of an inclined simply supported beam under moving loads for the effects of the velocity and the mass of a moving-load including the Coriolis and centrifugal forces and inclination angle of the beam.Omolofe 2013 has presented a procedure involving spectral Galerkin and integral transformation methods for the problem of the dynamic deflections of beam structure resting on bi-parametric elastic subgrade and subjected to travelling loads.Oni andAwodola 2010, andAwodola 2014 have studied the dynamic response of an elastically supported non-prismatic Bernoulli-Euler beam and the flexural motions of elastically supported rectangular plates carrying moving masses and resting on variable Winkler elastic foundations.For the application of the moving load problem in defence field, a finite element model of a high-speed projectile and barrel interaction has provided by Esen andKoç 2015a, 2015b.Another Finite element application of the moving load on plate dynamics can be found in Esen 2015.The dynamics of a railcar bogie when moving over a flexible bridge has been studied by Mizrak and Esen 2015 as an application in high-speed rail transportation.
Generally, studies in the field of moving loads have been carried out by researchers working on vehicle-bridge interaction.In such problems, the mass of the bridge system is often too large when compared to the mass of the vehicle.For this reason, the mass of the vehicle does not significantly affect the natural frequency of the overall system.However, in mechanical engineering applications, the mass of the moving load can be much larger than the mass of the transport system.This is the main difference of the moving loads in the mechanical engineering applications where the mass of the moving load may significantly affect the natural frequencies of vibration modes of the transport system.In addition, for the vehicle-bridge-interaction, the moving loads on bridges can be modelled as a moving oscillator due to the vehicle's suspension system; thus, an amount of vibration energy dissipates in suspension system.However, the moving load is in tight contact with the beam especially in robotic and CNC systems.Due to the accuracy of the machined work-piece, any relative motion is not allowed between the mass and beam in the vertical direction.Moreover, the vibration-induced vertical position error should be less in the micron level; it is reasonably expected due to the accuracy Latin American Journal of Solids and Structures 14 (2017) 485-511 of the work performed.Therefore, the velocity, acceleration, and jerk of the mass are vital in terms of the accuracy of such mechanical engineering applications.The main purpose of this study is examining the issue in this perspective.

MATHEMATICAL MODELLING
In Figure 1, a simply supported Euler-Bernoulli beam with a lumped mass m p which is travelling on the beam from the left end to the right end with a variable velocity v m (t), a variable acceleration a m (t) and a constant jerk J, is given.
In the formulation the following assumptions will be adopted (Fig. 1)  The mass inertia is considered. The mass is always in contact with the beam. The beam is thin and small displacements in the beam are occurred according to thin beam theory. The beam can be of variable thickness but the material properties are constant trough length of the beam. In vibration, the longitudinal and flexural motions are independent of each other but the axial force due to the acceleration of the mass affects the flexural displacements. The trajectory of the mass is defined by time-dependent x p (t) Based on the above assumptions, the motion equation of the beam due to the effect of the mass located at the time-dependent point x p on the beam is provided by Eq. (1) Fryba 1999: The left hand side of Eq. ( 1) represents the resisting internal stiffness, inertia and damping forces due to the external forces on the right side with the parameters that: N(x,t) represents the time dependent axial load due to the acceleration of the moving mass; ρ represents the density, A(x) represents the none uniform cross-sectional area, σ represents equivalent viscous damping coefficient, E represents the Young's modulus of elasticity, I(x) represents the area moment of inertia, x represents the central coordinate of the beam, t represents time, w z (x, t) represents the vertical displacement of the beam, mp represents the mass of the load, ( The initial and boundary conditions of the beam are: In order to determine time dependent displacements of the beam, a rough analytical solution of the motion Eq. ( 1) shown above, can be obtained through simplifications that ignore the effects of inertia and damping; and accepting that the cross-section area is uniform and the mass moves with a constant velocity.For simplified cases, omitting geometric and dynamic nonlinearities such moving load problems have been extensively studied in the literature by numerous researchers.When the mass moves with a constant jerk on the beam, the exact solution of the motion equation using analytical methods will be very difficult; and no study of constant jerky motion of a mass has been reported in literature, yet.This kind of motion problems are faced frequently in today's engineering applications such as high-speed precision machining, precision grinding of hard cutting tools, some robotic applications of precision surgical robots and precision robots used in microchip production; where vibration of the linear run-ways are vital for the accuracy of the work done.For an exact or an admissible solution in the engineering sense, there are needs of some new methods that represent mass motion with all effects including jerk.This paper presents a solution method representing the mass as a time dependent finite element and combining this finite element with the classical finite elements of the beam.This method includes both inertial and damping effects and gives a solution for both transverse and longitudinal vibrations of the beam.From this perspective, the proposed FEM combined method can be very useful for design and control engineers if they pay attention to the requirements of modelling as described below.For any time t, the effective transverse (z) force at coordinate x, due to the interaction between the moving mass and the beam which is vibrating and deflected, can be determined as follows using he principle of D'Alembert :

Modified Mass, Damping and Stiffness Matrices of the Beam Element
( ) For the linear motion of a particle object, the instantaneous acceleration is derived using the constant jerk of the motion and integrating it with respect to time for an increment of time Δt from t-Δt to t, as given below: Considering continuous contact between the mass and beam, when the mass moves on the beam with a variable acceleration and a constant jerk, the instantaneous flexural acceleration of the contact point of the mass can also be determined using the definite integration of the differential acceleration and jerk equation over a small time increment Δt, in the same way as Eq. ( 4).In such case, the effective transverse force equation will be rearranged as follows: ( ) When the necessary total differentiations of deflection function w z (x, t) with respect to the time dependent contact point x p of the moving mass at time t, are determined and after making some mathematical arrangements, Eq.(5) becomes: where " ′ " and " • " are, respectively, spatial and time derivatives of deflection.Derivation of the acceleration and jerk terms are given in Appendix A. The terms ( ) (6), represent the inertia, the Coriolis and the centripetal acceleration components at time t-Δt, while ( )    represent jerks for above mentioned acceleration components of the interaction of the mass and beam at time t.Where z w  is inertia; is centripetal acceleration components, while is Centripetal jerks of the accelerations respectively.In vibration, the inertia jerk of the beam can be approximated using the average acceleration method proposed by Newmark in the Beta method of integration with in the assumptions of the method.In this case the jerk of the beam is: When the beam is in vibration, the longitudinal (x) force component, between the moving mass and the beam, induced by the vibration and curvature of the deflected beam is ( ) The equivalent nodal forces of the beam element s under a lumped moving mass that has a variable acceleration and a constant jerk are: Latin American Journal of Solids and Structures 14 (2017) 485-511 Where, ϕ i (i=1-6) are shape functions of the beam element given by Clough and Penzien, 2003; these are cubic Hermitian polynomials that are expressed as where l the length of s-th beam element.The relation between shape functions and transverse and longitudinal displacements of the s-th beam element at position x and time t, is (Clough and Penzien, 2003): where usi (i = 1-6) are the displacements for the nodes of the s-th beam element on which the moving mass m p locates at time t.
The property matrices of a beam element without any moving mass, in Fig. 2, having both transverse and longitudinal nodal forces and deflections, are derived from the usage of the principle of virtual works.Introducing the principle of virtual displacements, and applying unit displacements at the nodal points and then equating the work done by the external forces to the work done on the internal forces: WE = WI (Clough and Penzien, 2003) for a uniform beam segment using the interpolation functions of Eq. ( 11), a stiffness equation can be obtained for the stiffness coefficients.For the mass coefficients, an elemental balance equation can be obtained using the relation between nodal accelerations and resisting inertial forces with the help of the same interpolation functions.In the special case of a beam element with uniformly distributed mass and under a mass moving with a variable acceleration over the beam element, the coefficients of modified mass, stiffness and equivalent damping matrices of the beam element s can be obtained by taking into account the contribution of inertial, Coriolis and centripetal forces due to the related accelerations and jerks (7) induced by the moving mass.In such a case, the total force balance equation for the beam element s is obtained in matrix form as below: where m , c and k are, respectively, the modified mass, damping and stiffness matrices of the beam element s that includes all the effects of the moving mass, while m, c and k are consistent mass, equivalent viscous damping and consistent stiffness matrices of the beam element itself respectively, Latin American Journal of Solids and Structures 14 (2017) 485-511 and G k is the geometric stiffness matrix due to the axial force N(x,t).The matrices w m  and w m  , w c  and w c  , w k  and w k  are derived from the interaction of the moving mass and the beam element and they represent the moving mass with a constant jerk.The components of these matrices vary depending on the motion velocity, acceleration and jerk, and the position of the mass on the beam element.The coefficients of these matrices are obtained from the equivalent nodal force equation ( 10) of the beam element s.One can obtain these matrices using the relation between the nodal displacements and the deflection functions ( 12) with the interpolation functions (11), and then putting them in the nodal force equation ( 10).After obtaining the spatial and time derivatives of the deflection function, Eq. ( 10) is rearranged separately for acceleration and jerk in terms of nodal accelerations, nodal velocities and nodal displacements, thus these matrices are obtained, as given in Appendix B.
During the integration of motion equation of the overall system, the instantaneous values of these matrices are calculated at every time step Δt, an then added to the values of the matrices of the beam element s on which the moving mass locates at time t.Because, they are derived from the relation between nodal displacements and shape functions of the two-node beam element, the dimensions of these matrices are equal to the dimensions of the mass, damping, stiffness matrices of the beam element.Hence, the beam element has three displacements DOF at each end nodal point; the dimensions of these modification matrices will also be 6x6.
The consistent mass, equivalent damping and consistent stiffness matrices m, c, and k, of the beam element s are calculated using their coefficients m ij , c ij, and k ij , as given below ( , 2, 3, 5, 6), ( ) ( ) ,( , ( , 2, 3, 5, 6), Where y and t are the coefficients for the equivalent viscous damping.They can be determined by using Rayleigh damping theory in which the coefficients of damping matrix is proportional to the combination of coefficients of the mass and stiffness matrices, as given by where ζ i and ζ j are the damping ratios of the structural system for any corresponding natural frequencies ω i and ω j .
Latin American Journal of Solids and Structures 14 (2017) 485-511 In order to include the flexural effect of an axial load N(x,t) due to the acceleration of the mass, that is negative for positive acceleration and vice versa, the geometric stiffness matrix G k are obtained with the following coefficients: , 0 ( , ) ( ) ( ) , ( , 2, 3, 5, 6), where ,

Gi j
k are the coefficients of consistent geometric stiffness matrix and the Hermitian interpolation functions.Eqs. ( 11) are used in deriving the geometric stiffness coefficients.
In the time dependent elemental equation of motion of the beam element s, in Eq. ( 13a), the force t f is as given below: ( ) ( ) The nodal acceleration, velocity, and displacement vectors, , ,and

Equation of Motion of the entire System
The equation of motion for the multiple degree of freedom damped structural system, shown in Fig. 2, is given by where M, C and K are, respectively, the overall mass, damping and stiffness matrices, while u  , u  and u are respectively, the acceleration, velocity and displacement vectors.Besides, F is the overall external force vector of the system at time t.
In general, such a structural system shown in Fig. 2, one can obtain overall mass M, stiffness K and C matrices by assembling its element matrices and imposing given boundary conditions.If there is a mass moving with a variable acceleration over the structure, the mass, stiffness and damping matrices of entire system can be obtained by taking into account the contribution of inertial, Coriolis Latin American Journal of Solids and Structures 14 (2017) 485-511 and centripetal forces, and jerks induced by moving mass.This is done by modifying only the elemental matrices of the beam element s on which the moving mass locates at time t, and keeping the elemental matrices of the other elements the same as original.The instantaneous overall force vector is also time-depended.The coefficients of overall force vector are equal to zero except the nodal forces of the beam element s.Thus, the instantaneous overall force vector of entire system becomes as below: In evaluation of the matrices of the beam element s, for the instantaneous position of the mass, the term x in the Eq. ( 11) is replaced with x m (t).The instantaneous values of x m (t) and s can be determined as follows:

Solution of Equation of Motion
If the added matrices of w m  , w m  , w k  , w k  , G k and the equivalent force vector f are zero, the undamped natural frequencies and vibration mode-shapes of the beam are obtained from homogenous solution of ( 19).In such a case, the Equation ( 19) reduces to The solution of ( 22) can be assumed as Introducing the latter u and u  into (22) gives This is a set of homogeneous equations, for which a nontrivial solution exists if The above determinant equation is satisfied for a set of n values of frequency ω 1 , ω 2 ,.., ω n .The frequency ω i is called the i-th natural frequency.Substituting all the ω i into (24) yields the corresponding set of vectors 1 Ψ , 2 Ψ ,…, n Ψ that satisfy this equation.The ith vector i Ψ corresponding to the ith natural frequency is called the ith natural mode, or ith mode shape [6].
Latin American Journal of Solids and Structures 14 (2017) 485-511 For such a system, given in ( 19), one can obtain a solution by using a numerical integration method like Newmark's β method [5].Using the Newmark integration method, equation ( 19) is transformed to equivalent stiffness and force equations.Thus, for the n-th time step at time t n (=t n-1 +Δt): where are, respectively, the initial conditions for the accelerations, velocities and displacements of the structural system at time t = t n .Displacements at time t n 1 ( ) ( ), Accelerations and velocities at time t n ( ) where 0 1 2 3 4 5 6 7 , , , , , , and a a a a a a a a are the integration constants in the Newmark-method; that are where β and γ are the integration parameters.In this study, β=0.25 and γ=0.5 are used for integration accuracy and stability.
If the mass travels on the beam, the instantaneous mass and stiffness matrices should be used for the solution of instantaneous natural frequencies of the entire system.For the frequency solution of (32), one can use (33).That is: where i w is the ith forced vibration frequency of the entire system.If the mass and stiffness matrices, used in (33) are time-dependent, the frequency solution will also be time dependent.
For the calculation of the instantaneous overall mass and stiffness matrices of the entire system at every time step of ∆t, one may use following steps: 1. Determine the mass and stiffness matrices of each beam element.2. For time t, determine the element s on which the moving mass locates with (21).
3. Determine x m (t) which is the time dependent position of the moving mass on the s-th element with (21).4. Calculate the time dependent shape functions with (11) by substituting the value x=x m (t) which is defined in the previous step.5. Calculate the modified mass m , damping c and stiffness k matrices of the s-th element with the help of ( 14), and Appendix B. 6. Calculate the instantaneous overall mass and stiffness matrices of the entire system by combining the mass and stiffness matrices of each beam element.Then impose boundary conditions.Eigen solution of these matrices gives instantaneous natural frequency of the entire system at time t.For t+∆t go to step 2

RESULTS AND DISCUSSION WITH NUMERICAL SOLUTIONS
In this article, the Newmark β method of integration [27] with the time step size =0.00001s is used to obtain the solution of Eq. ( 19).The integration constants that manage the sensitiveness and stability of the Newmark procedure, β=0.25 and γ=0.5 are imposed.It has been reported that when β takes 0.25 value and γ 0.5, this numerical procedure is unconditionally stable Szilard 2004.
Example 1.In order to compare present method with the others at the same assumptions, excluding the inertia effects of the mass and accepting the mass as moving load that has been studied in literature, let us take a simple supported isotropic beam-plate under a F = 4.4 N moving load.The dimensional and material specifications of the plate are identical with those preferred in Meirovitch 1967, i.e. l x = 10.36 cm; l y = 0.635 cm, h = 0.635 cm; E = 206.8GPa, ρ = 10686.9kg/m 3 ; T f =8.149 s, where T f is the fundamental period.In Table 1, dynamic amplification factors (DAF), which are defined as the ratio of the maximum dynamic deflection to the maximum static deflection, are compared with several previous numerical, analytical, and experimental results available in literature.It is noted that T is the required time for moving load to travel the beam.It is seen that the results obtained by the current finite element (column 3) are very close to the analytical solution Meirovitch 1967, and  (1) Present method.
(5) From Michaltsos 2002.(6) From Mohebpour et al 2011 Example 2: For another example, a frame box-section beam that can be preferable for a CNC system due to the lightweight and rigidity, with length L = 1 m, height h=0.1 m and width b = h = 0.1, and a wall thickness of 7.5 mm, was chosen.Since h / L = 10, it can be assumed that it is a thin Euler Bernoulli beam.In order to satisfy the desired computational accuracy in the Finite Element Method, firstly the number of elements should be determined.For this reason in Tables 2-5, the first six natural frequencies depending on the number of elements are given.In this paper, the number of elements is preferred as 100 since its results in the frequency calculation are very close to the exact ones.However, one can prefer 10 elements for an acceptable result in the analysis, in that case, for the sixth frequency in Table 2 the relative error percentage is about 0.8 percent.Where, the relative error=100(69797.8459-69247.9106)/69247.9106=0.79,and the error will be very little and insignificant for the lower frequencies.In moving mass problem, the variation of the natural frequencies of the beam is an interesting phenomenon that every design engineer should take into account of this effect of the mass.When the mass moves on the beam depending on the position of the mass, the natural frequencies undergo a continuous change, thus the vibration of the beam has to change as can be seen in Figures 3-7 in which normalized frequencies fr versus normalized position x p (t)/L of the mass over the beam are depicted for different end support conditions.Where fr is ratio of the frequencies in cases of with and without mass on the beam.Figure 3 displays the variation of the first four natural frequencies depending on the mass position on the beam and the mass ratio of the moving mass and the mass of the beam.Because of this variation in frequencies, the vibration behaviour of the beams under moving loads are very complicated and depends upon various parameters of the system such as mass, velocity, acceleration, and jerk of the moving mass including end support conditions.The variation of a particular natural frequency is governed by its corresponding natural mode shape.For example, as can be seen from Figure 3 the variation of the first natural frequency resembles its first natural mode shape.This is also valid for the other natural frequencies, 2 nd , 3 rd , 4th and so on.Another interesting behaviour is that at the location of the nodes, there are no changes of the frequencies even for higher mass ratios.Between the nodal points of the corresponding mode shapes, the rate of change is increased by the increasing mass ratio.These variations are also nonlinear, and have a tendency to the left or right as seen from the graphics depending on the first and last nodal point that are near to the left and right ends.The other nodal points except for these there are change but there is no tendency to the right or to the left.This event is valid for other end support conditions that are presented through Figures 3-7.The determination of the change in the natural frequencies is easy by using the proposed method in this study.Because the added matrices of w m  , w m  , w c  , w c  , w k  , w k  and G k are always added to the left side of the motion equation ( 19) of the whole system.This is the most advantageous property of the proposed method that one can directly model a dynamic system by including the inertia, Coriolis and centripetal forces and jerks of the moving mass as internal forces, and imposing the gravitational force to the right as an external force.The above matrices represent the dynamic added mass, damping and stiffness properties of the beam under a moving mass.The main difference of the assumption of the moving mass is this modification of the unloaded beam matrices.In the moving load assumption, all these matrices are zero except for the stiffness matrix G k due to axial force effect in the flexural deflection of the beam if there is an accelerating mass.Including second and third derivations of the deflection function one can analyse the effects of mass moving with constant and variable velocity on a beam using the current method.For the same mass ratio, i.e., e=0.5 the effect of the velocity of the mass on the change of the natural frequencies is given in Figure 7, for a fixed-free support condition.While the velocity is small, the effect of the velocity in the added matrices will also be very small because the added terms from velocity are too small compared to the high flexural rigidity EI of the beam.For the case, that the beam flexural rigidity is higher and velocity of the moving mass is constant the effect of the velocity on the natural frequency can be omitted even for 100 m/s the effect will be at the level of one percent.If the beam is slender or for the case of cable transport the effect of the velocity should be taken into account even though the velocity is small.For barrel and projectile interaction in a weapon system, that has a longer barrel and high exit velocity of the projectile, and rocket launching structures the effect of the velocity is very important due to high rate of acceleration and jerk.
Latin American Journal of Solids and Structures 14 (2017) 485-511 Moving load problems have been studied widely in literature, but most of these studies are for vehicle bridge interaction; and the most of them have omitted the inertial effect of the mass and accepted only the gravitational force of the moving load.Some of them have been modelled as moving oscillator that the main mass of the load is modelled as a sprung mass, and the unsprung mass of the suspension system is omitted.However, for the application of CNC applications of the moving mass where the accuracy of the work piece is important all the mass of the moving load is firmly in contact with the beam, the dynamic amplification factor (DAF) for different masses of the moving load may be significant for design engineers.Where DAF= maximum dynamic displacement of the beam/static displacement of the beam when the mass is at the middle.Figure ( 8) displays the DAF's, and vibration acceleration and velocities of the midpoint for different constant travelling velocities of the mass and different mass ratios e, where e= mass of the moving load/mass of the beam.The mass of the given beam is 22.2 kg; and for e=0.2, 0.6, 1.0 and 1.4, the masses of the moving load are: 4.44, 13.32, 22.2 and 31.08 kg, such a mass ratio applies to CNC applications.Figure ( 8a) is for lower velocities while the others are for higher velocities.The figures (8e, 8f and 8g) are for moving load assumptions which the inertia, Coriolis and centripetal effects of the moving mass are omitted and only gravitational effect is considered, while Figures (8a to 8e) are for moving mass assumption that consider all the effects of the moving mass.From the figures, one can realize that they are very different in terms of the mentioned assumptions.The maximum DAF is about 1.7 for moving load assumption, but for the moving mass assumption it is about 6 that is higher about 3.53 times.This is due to the effect of mass with higher travelling velocities.For very low travelling velocities the difference is not considerably high when the rigidity of the beam is high in comparison with the velocity due to its little effect in the system mass, stiffness and damping matrices.The velocity increment in Figure (8a) is 0.001 m/s, while it is 5 m/s for the others.When Figure 8a, is examined closely, it is observed that there are many resonance peaks and cancellation points.In both cases the vibration acceleration and velocities are very different from each other.In the case of the moving mass in Fig. 8d, the maximum is ten times higher when compared to the moving load case in Fig. 8g.In addition, although the increase in the mass ratio is linear, the increase in the maximum of the moving mass is not linear, but the increase in the moving load is linear in acceleration and speed, but no increase is observed in the DAF that there is no difference with increasing mass ratio e.The effect of the mass on the frequency variation of the beam given in Figure 3-8 is valid for the moving mass assumption.At constant acceleration of the mass, the displacements and rotations at the midpoint of the beam are given in Fig. 9.The increase in displacement seems to be very high at high acceleration of motion.However, such high acceleration cannot be the case for CNC applications.In addition, from the standpoint of existing bridge engineering problems, high acceleration does not apply in existing applications.However, the acceleration effect can be accounted for in high speed rail transport applications.Moreover, if the projectile is considered that accelerations to be higher than the given examples, the impact of the exposure should be considered in such applications, in the case of projectile and barrel interaction.The jerk can be valid for a sudden acceleration increase takes place in braking or accelerating due to the applied force chances over time.Starting with zero acceleration at different jerks, displacements of the midpoint of the beam are illustrated in Figure 10.For the movement with little jerks, Figure 10 a, the displacements are not significantly changed but for higher jerks, Figure 10b, as it can be in projectile motion in a barrel it is observed that displacements are considerably increased.When the mass moves with an acceleration due to an axial force is applied at the contact point of the mass, the beam vibrates in axial direction.Figure 11 shows the axial displacements of the midpoint for e=1, and case of different lower accelerations in (a), and case of different higher accelerations in (b).Higher accelerations are valid for projectile and barrel interaction of a weapon system like a tank that the exit velocity of the projectile is about 1750 m/s.One can calculate the average acceleration of the projectile that starts from zero velocity to a 1750 m/s exit velocity in a 6 m long barrel, the average acceleration is 2.55e5 m/s 2 .For lower accelerations that can be applicable for CNC applications the midpoint displacements are increased by the increasing of the motion acceleration of the mass.The rate of displacement increase should be considered when the accuracy of the work is important.In terms of bridge engineering applications, the effect of axial forces may be significantly lower due to the motion acceleration of the vehicle, but this effect may be important for asphalt life.Another application of accelerated motion is rocket launching, where a massive mass can experience very high acceleration because of the ignition of a rocket.

CONCLUSIONS
In this study, the inertial, Coriolis and centripetal effects of the moving mass are determined in terms of variable acceleration and constant jerk.The classical beam element equations are modified from the terms that come from the moving mass.The variable acceleration and jerk of the mass are combined the accelerations and jerks of the beam element on which the mass locates at time t, at the variable contact point of the mass.This method, in which only the contact element is replaced and the other elements of the beam remain the same, is a very advantageous tool for calculating many motion patterns, including jerk effect.In the formulation of the time-dependent force equation of the contacted beam element, the effects of velocity, acceleration and jerk are separately determined in order to investigate separately the effects of them.Some interesting results of the different cases are presented in terms of the usage field of the moving load problem, including the possible application of the theory in design of runways of CNCs, vehicle-bridge-interaction, in design of the barrels and rocket launching structures under the effect higher acceleration and jerks, etc.The general approach in Bridge design, the DAF is criticized, in terms of omitting the mass inertia and accepting the moving mass as a moving load that accepts the only gravitational effect of the mass.It is shoved that, for lower velocity of travelling of the mass, this approach makes not very much difference, but for higher velocities and accelerations, including the inertial effect, the DAF can be three times greater depending on size of the velocity, acceleration, and jerk.One drawback of the proposed method is that it is not usable for hand calculations due to the modification of the system matrices in each time step of the integration.
Latin American Journal of Solids and Structures 14 (2017) 485-511 APPENDIX A Derivation of Eqs. ( 6) and ( 9): When the beam is in vibration, the longitudinal (x) force component, between the moving mass and the beam, induced by the vibration and curvature of the deflected beam is ( ) The axial acceleration of the beam is derived from the second order total differentiation of its deflection function ( , ) When the beam is in vibration, thetransverse (z) force component, between the moving mass and the beam, induced by the vibration and curvature of the deflected beam is ( ) The vertical acceleration of the beam is derived from the second order total differentiation of its deflection function ( , ) In Eq. (B.1b), the first three terms are for a two-node beam element that has three DOF at each nodal point with the axial, transverse displacements and slope.They represent the mass, damping and stiffness property matrices of the beam element; can be found in [30].The other terms stand for the interaction with the moving mass with a constant jerk.In such a case, the equivalent force equation for the beam element that includes only these terms can be obtained as below:           , 3

Figure 1 :
Figure 1: A simply supported beam and a moving mass system.
length of the beam by the moving mass, g and δ represent gravitational acceleration and the Dirac delta function, respectively, and d 2 w(x p ,t)/dt 2 represents the acceleration of the beam at the contact point of lumped mass.

Figure 2 Figure 2 :
Figure2displays mesh discretion of a thin beam under a lumped moving mass and the beam element s on which the mass mp applies, at time t.The beam element s has three equivalent nodal forces and displacements at each nodal point.The time-dependent global position of the moving mass in the span is xp (t), while local position on the length of the element s is x m (t).The beam has n elements and (n+1) nodes.

Figure 3 :
Figure 3: Variations in natural frequencies for pinned-pinned beam.

Figure 4 :
Figure 4: Variations in natural frequencies for fixed-pinned beam.

Figure 5 :
Figure 5: Variations in natural frequencies for fixed-fixed beam.

Figure 6 :
Figure 6: Variations in natural frequencies for fixed-free beam.

Figure 7 :
Figure 7: Variations in natural frequencies for fixed-free beam for e=0.5 and v=10, 20 and 30 m/s.

Figure 8 :
Figure 8: Comparisons of DAFs in moving mass (a, b, c and d) and moving load (e, f and g) assumptions.

Figure 9 :
Figure 9: Midpoint displacements and rotations for e=1, and different constant acceleration.

Figure 11 :
Figure 11: Axial displacements of midpoint for e=1, and different constant accelerations.
Lateral and transverse Force equation of the beam element s can be obtained as summing up all the forces including the accelerations and jerk terms of Eqs.(A.2 and A.4) . (B.2) is arranged in the matrix form, the following elemental equation is obtained with the property matrices given in Eqs.(B.4 -B10). 1

Table 2 :
Free vibration frequencies for pinned-pinned thin beam.

Table 3 :
Free vibration frequencies for fixed-pinned beam.

Table 4 :
Free vibration frequencies for fixed-fixed beam.

Table 5 :
Free vibration frequencies for fixed-free beam.