Dynamic of a cylindrical column on a vibrating base

The solution of the problem of cylindrical column on a rigid horizontal foundation motion is presented. It is assumed that the base moves in a horizontal direction, there is no slippage between the column and the base, and the column height is significantly greater than the base diameter. The modes of oscillations of the column, corresponding to the main resonance for two laws of the base motion are found. The results are presented in the form of amplitude-frequency characteristics. It is shown that the solutions corresponding to the column motion in the plane of vibration excitation are unstable.


Introduction
The study of the motion of a rigid body on a horizontal supporting surface is of great interest for applications in technology. The problem of vibrations of a rigid body in the form of a parallelepiped on a vibrating base was first considered in [1]. In this and a number of other works [2][3][4], this problem was considered in relation to the problems of seismic stability. There are other technical areas of application of this problem, for example, the problem of the movement of parts on a conveyor. The exact solution to the problem of vibrations of a parallelepiped on a vibrating base was obtained in [5], and in [6,7] the averaging method was used, which made it possible to investigate the stability of the solutions obtained. A detailed list of publications on the issue of seismic stability of elongated structures is given in [8].
In his article in 1963, G. Housner [1], analyzing reports on a number of earthquakes, drew attention to the fact that structures that are highly elongated in height can safely withstand earthquakes, despite their unstable appearance. The explanation for this phenomenon lies in the existence of stable vibration modes, in which the body consistently rests on the ribs of the base. In this work and a large number of subsequent ones, a body in the form of a parallelepiped was considered as a model, although most of the structures (water tanks, stone columns) have a base in the form of a circle.
In this paper, the dynamics of similar systems, but with a base in the form of a circle, is investigated.

Formulation of the problem
Consider the motion of a rigid body in the form of a column with circular cross-section on a horizontal supporting surface, performing harmonic oscillations in the horizontal direction. A rigid body on a support surface in the absence of slippage is a system with two degrees of freedom. In the equilibrium position, the central axis of the column coincides with the normal to the supporting plane, and the bottom base rests completely on the plane. When motions, the axis of the column deviates from the vertical by a certain angle (Fig. 1). We assume that there is no slippage at the point of contact between the column and the plane, and, in addition, the column cannot rotate about the vertical axis. This formulation of the problem corresponds to the case of taking into account the forces of interaction between the base and the column according to the Coulomb friction model, with a finite area of the contact patch and a sufficiently large value of the friction coefficient.
We will also assume that the column is strongly elongated in the vertical direction and the ratio of the base radius to the half-height of the column is small ℎ ⁄ ≪ 1. We take , -the angles of rotation of the column axis relative to the axes , associated with the movable base as a generalized coordinates. From the last assumption , ≪ 1 and we can put . The kinetic energy of the system is 2 ⁄ . For a column extended in height, we put . The potential energy of the system Π √ℎ cos − , where is the mass, is the gravitational acceleration, is the angle between the axis of the column and the line connecting the center of mass of the column C with the point contact. Generalized forces: where is the vibration frequency of the supporting surface, , , are the amplitudes and phase of the base motions in two orthogonal directions.
Substituting the obtained formulas into the Lagrange equations of the second kind, we obtain the equations of motion of the system: where, the dimensionless time is introduced, where ℎ/ , the dimensionless amplitude and frequency of oscillation excitation are ⁄ , ⁄ , 1,2. The dot denotes the derivative with respect to the dimensionless time. Dissipation according to the model of viscous friction with a dimensionless friction parameter is introduced into the equations of motion.
The system of Eq. (1) describes the motion of a strongly nonlinear system. With 0 and 0, we obtain the equation of column vibrations on a fixed plane, which has two solutions. The first solution corresponds to the vibrations of the column in the plane located at an angle to the axis. For one oscillation period, this solution has the form [6]: where , are the amplitude and frequency of natural vibrations, interconnected by the dependence: In this case, the base of the column rolls in a circle with center on the reference plane (Fig. 1). After substituting Eq. (4) into Eq. (1), we obtain the dependence of the oscillation amplitude on frequency for this solution: Next, we consider the problem of forced vibrations of a column on a moving base in the vicinity of the generating solution Eq. (4). In new variables, this solution = 0, = , and are constants related by relation Eq. (5). We restrict ourselves to solving the system of Eq. (6) for two types of base motion.

Base motion in a circle
Consider the special case = = , = − /2, which corresponds to the motion of all points of the supporting plane along a circle. For steady-state forced vibrations, there is an exact solution to the problem. Setting = + , we obtain a system of equations that does not explicitly contain time on the right-hand side: The solution of the system Eq. (7) is: = 0, = , In this case, the column has a constant angle of deviation from the vertical, and the point of contact moves in a circle with a constant angular velocity equal to the angular velocity of movement of the base. The last two formulas Eq. (8) are amplitude-frequency and phase-frequency characteristics, respectively. To study the stability of the obtained solution, we consider a matrix consisting of partial derivatives of Eqs. (7) with respect to , , and . After substituting = 0, = from solution Eq. (8), this matrix has the form: The asymptotic stability was investigated using the Routh-Hurwitz criterion.

Base motion along a straight line
Consider one more special case = 0, in which all points of the support plane move along a straight line. As in the previous case, we put = + . The system of equations has the form: The system of differential Eq. (9) explicitly contains the fast variable = 2 on the right-hand side. We will use the averaging method [9,10], for which we assume that the value of the viscous friction coefficient 2 ⁄ , the excitation amplitude 4 ⁄ , the angular velocity 2 ⁄ and the frequency detuning − 2 ⁄ are small. In the vicinity of the generating solution, the term 1 + − in the first equation is small. The functions , μ, , are slow variables of the dimensionless time . Applying the averaging procedure and keeping the previous notation for the averaged variables, we obtain: For steady-state oscillations, we assume the left-hand side of equations Eq. (10) to be equal to zero. The solution of the resulting system coincides with Eq. (8) taking into account the replacement = . The amplitude-frequency and phase-frequency characteristics coincide with those shown in Fig. 2  has the form: ( ) has period 2 ⁄ : The second equation in Eq. (11) has a solution = ⁄ and system can be written in the form: The solution corresponding to the motion of the column in the plane is realized only when = 0 . When the initial velocity is nonzero, the modulus of the angle increases indefinitely. This means that the motion of a column with a round base in one plane is unstable.

Conclusions
At small values of the acceleration amplitude of the base < or < and zero initial conditions, the column remains in a vertical position, the movement of the base only leads to a redistribution of the pressure of the lower edge of the column to the support plane. Under arbitrary initial conditions, or large amplitudes of excitation, the column can overturn. From the analysis of the results obtained, it follows that at low frequencies of vibration excitation, there are only these two variants of the motion of the system. When the excitation frequency exceeds a certain critical value, the value of which depends on the parameters of the problem, a stable solution exists for the system under consideration (Figs. 2 and 3). In this case, the base of the column rolls in a circle on the reference plane. Thus, even with large amplitudes of base motion, the column may not tip over.
The form of the amplitude-frequency dependences for a column with a round base and a parallelepiped [6] qualitatively coincide, but these are completely different curves adjacent to dependences Eqs. (5) and (3), respectively.
The assumption introduced at the beginning of the article about the absence of slippage at the point of contact of the column with the base reduces the scope of the results obtained. The problem requires additional research for various models of forces and moments of friction at the point of contact.