Abstract
In this paper, we consider the dynamical description of a pendulum model consists of a heavy solid connection to a nonelastic string which suspended on an elliptic path in a vertical plane. We suppose that the dimensions of the solid are large enough to the length of the suspended string, in contrast to previous works which considered that the dimensions of the body are sufficiently small to the length of the string. According to this new assumption, we define a large parameter and apply Lagrange’s equation to construct the equations of motion for this case in terms of this large parameter. These equations give a quasi-linear system of second order with two degrees of freedom. The obtained system will be solved in terms of the generalized coordinates and using the large parameter procedure. This procedure has an advantage over the other methods because it solves the problem in a new domain when fails all other methods for solving the problem in such a domain under these conditions. It is one of the most important applications, when we study the slow spin motion of a rigid body in a Newtonian field of force under an external moment or the rotational motion of a heavy solid in a uniform gravity field or the gyroscopic motions with a sufficiently small angular velocity component about the major or the minor axis of the ellipsoid of inertia. There are many applications of this technique in aerospace science, satellites, navigations, antennas, and solar collectors. This technique is also useful in all perturbed problems in physics and mechanics, for example, the perturbed pendulum motions and the perturbed mechanical systems. The results of this paper also are useful in moving bridges and the swings. For satisfying the validation of the obtained solutions, we consider numerical considerations by one of the numerical methods and compare the obtained analytical and numerical solutions.
1. Introduction
The pendulum models have attracted scientists and researchers with many descriptions of motions and their analysis as important examples in physics and theoretical and applied dynamics. The most important pendulum motions come from the moving of a heavy particle suspended a light rod which is jointed pivotally at a point on the -axis which rotates by an angular velocity about the horizontal fixed axis. In [1], the authors presented the elastic pendulum problem. They derived the equations of motion and gave real-life examples for elastic pendulum motions. Equilibrium states and trivial cases are studied. In [2], a harmonically excited spring pendulum motion with 3 degrees of freedom is considered. The authors were connected to a direct transversely tuned absorber. The authors used the multiple scale method for solving the system of equations of motion. They used the phase-plane method for applying the stability of the motion. They studied the influence of the tuned absorber and the parameters of the model on the change of the motion numerically. In [3], the elastic pendulum with resonance appearing at cubic approximation in the Lagrangian function was investigated. The modulated equations of motion were obtained. The authors used the Hamiltonian reduction technique to describe the dynamical behavior of this system. In [4], the approximated damping elastic pendulum model was considered. The authors derived the governing equation with six degrees of freedom. The natural frequencies are derived from the linearized system of equations. The stability diagrams are used to show regular periodicity. Due to damping, there are unstable waves (irregular) depending on the damping level and the amount of parametric excitation for the natural frequency. The article [5] presented the rotational motion in a plane for a mathematical pendulum when its axis swinging vertically and horizontally, so the axis path is an ellipse nears a circle. The study is considered as precise spin solutions for a circular axis path case with zero gravity. The existence of the stability conditions of these solutions is considered. Supposing that the excitation amplitude is not small and the coaxial path contains a small ellipsoid. Approximated solutions for both high and low linear damping are found. A comparison of approximate and digital solutions was performed for different values of the damping factor. In [6], the governing system of two generalized coordinates of a nonlinear dynamic model with the spring pendulum damped motion is considered in presence of the flow of an inviscid fluid. The motion control system is attained by applying the equations of Lagrange. The multiple scale procedure is used to solve the equations of motion for this system. The authors obtained approximated solutions of the 2nd approximation. The resonance cases and the steady-state ones are studied. The stability procedure is given using the phase plane diagrams. In [7], the dynamics of nonlinear multiple degrees of freedom for a spring pendulum model moving in an auxiliary circle for the elliptical path are considered. Using the method of multiple scales, the authors solved the system of equations of motion. The temporal history of the solutions obtained and the phase plane projections are given to explain the dynamic behavior of the mentioned system.
In [8] the authors studied the elastic pendulum motion with resonance appearing at the third order in its Lagrange’s approximation. The governing dynamical equations of motion and Hamiltonian reduction are obtained for this system. In [9], the frequencies of the elastic pendulum oscillations are in the ratio 2 : 1 which is named the pulsation. The dynamics of this problem and the modulation equations for the resonant motion with small amplitudes are obtained as 3 wave equations. The complete solutions are obtained by using the Hamiltonian reduction. The phases, amplitudes, and precession angles of the solutions are given. The validity of the obtained results is given in high accuracy by numerical experiments. The article [10] presented a dynamical behavior of an elastic pendulum model. The author assumed that the pendulum moves in a vertical plane which rotates with uniform angular velocity . The author used Lagrange’s equations to find the differential equations of motion. The obtained equations are reduced to a nonlinear system of 2nd order which is solved by the small parameter perturbed technique. Numerical considerations are obtained through the Matlab program. These considerations show the validity of both analytical and numerical solutions.
In [11], the nonlinear multiple degrees of freedom response is investigated for a normal dynamic system given by the spring pendulum moving in an elliptical path. The Lagrange equations were used to construct the governing system of equations for the motion. The multiple scale method is used for deriving the approximated solutions of the resonance system. The steady-state motions and the stabilities of the solutions are considered. The authors in [12] studied the problem of spring pendulum dynamics in the presence of pendulum absorber using the theory of nonlinear normal patterns and asymptotic numerical procedures. They investigated the dynamics of the pendulum for both low and high vibration amplitudes. Using different methods, the stability of the vibration patterns is analyzed. The authors in [13] studied the relative periodic solutions for an elastic pendulum motion of a rigid body with a corresponding suspended point moves on the auxiliary sphere of an ellipse. The Lagrange’s function is applied to construct the equations of motion which are solved using the perturbed small parameter method. Numerical considerations are given using the Runge-Kutta method to prove the validity of the solutions.
In this article, we consider a heavy solid suspended in a string which suspended to an ellipse in a uniform vertical plane. The equations of motion for this mathematical pendulum in terms of the two degrees of freedom are obtained, and a large parameter depends on the model properties is assumed. By the definition of the large parameter, the approximated periodic solutions are obtained using the large parameter perturbed procedure. The accuracy of these solutions is investigated through a numerical technique and computerized programs.
2. Formulation of the Problem
Let us consider the coordinate system on which -axis is horizontal and -axis is downward vertical. A rigid body with mass is suspended by a string of length which attached to a point moving on the ellipse. The point on the circle, which has the same center of the ellipse and its radius equals such that the line is normal to the largest axis of the ellipse, is called the corresponding point of the point on that circle. When the point moves on the ellipse, the point moves on the circle with a constant angular velocity (i.e., the angle between the line and the -axis must depend on time only). Let us assume that the angle between the straight line passing through and the center of gravity of the body and the vertical will be denoted by , the angle of deflection of the string from the vertical straight line passing through will be denoted by , and the distance between the center of gravity of the body and the point will be denoted by (i.e., ) (see Figure 1). Choosing the body system of coordinates such that passes through , is orthogonal to and lying in the plane , while is orthogonal to .
We assume that the axes of the system represent the principal axes of inertia of the body. The principal moment of inertia concerning will be denoted by . The motion of the body starts at the instant in the anticlockwise direction; after time , the point is expected to make an angle with the vertical axis.
The coordinate of the point is () which takes the form
The kinetic energy of the system is which takes the form
The potential energy of the system can be expressed as follows:
The Lagrangian of the system is [14]
It can be written as follows:
Let us define the following parameters. where , , and is a large parameter.
We introduce the variables so that and .
The new Lagrangian function of the system can be written as follows:
For a small value of , we can expand the sines and cosines in the form of power series expansions to obtain
3. Equations of Motion
According to Lagrange’s equations [15], where are the generalized coordinates and are the generalized velocities, and we take the coordinates and .
Introducing the parameters where is the acceleration of gravity; the system of differential equations of motion can be written as follows:
The expressions and are expected to be functions of a large parameter and dependent on the values of , , and .
4. Approximate Periodic Solutions
Now, to find the perturbed solution of the nonresonance case up to the second approximation, we apply the method of the large parameter [16] in the following form.
Substituting from (15) into (13) and (14) and equating to zero each coefficient of , we find that
Coefficient of
Coefficient of
Coefficient of
Coefficient of
Making use of the equations (15)-(20), the required periodic solutions take the form
Neglecting the secular terms [17] in the formulas (21) and (22), the constants and will vanish; then,
5. Discussion of the Results
In this subsection, we give a parametric analysis of the obtained results for the behavior of the obtained analytical solutions and as functions dependent on the time , the values, , , , and , and the large parameter ; that is,
We note that the domain of the obtained solutions under the assumed conditions is as follows: while the domain of the previous corresponding cases is as follows:
This means that the obtained solutions are treated in a new domain which is considered as a complement space for the previous work domains.
In what follow, we study the validity of both analytical and numerical solutions. Using one of the numerical methods, we obtain numerical solutions and give more analysis of the results. The graphical representations for both solutions are obtained through computer programming.
6. Computerized Data
This section is devoted to ascertaining the accuracy of the solutions being considered in Sections 3 and 4. Computer programs are developed for the representation of the obtained analytical solutions and and their derivatives in a definite period. On the other hand, the fourth-order Runge-Kutta method [18] is applied to solve the autonomous system (13) and (14) for obtaining the numerical solutions and their derivatives and investigating their graphical representations. The characteristic curves for both the analytical and numerical solutions and their derivatives appear in Figures 2–5. From Figures 2 and 3, we deduce that when increases, the number of the waves and their amplitudes of the solutions and decrease and vice versa, while from Figures 4 and 5, we deduce that when increases, the number of the waves is the same as the solutions and but their amplitudes decrease and vice versa. The phasing diagram procedures [19] for the stability of the solutions are given in Figures 6–11. From these figures, we deduce that the solutions and their derivatives are stable but their amplitudes decrease with increasing . In all stability figures, the curves take the cardioid form except Figure 9. We deduce from Figures 2–5 that the analytical and numerical solutions are in full agreement, and Figures 6–11 show that the stability diagrams are largely identical. That is, the numerical results are largely correspondence with the analytical ones which prove the validation of both solutions and the used techniques.
7. Conclusion
This motion is very important as a mathematical pendulum model for many problems in fluid materials and gas translation in big and small tanks and is considered as a generalization for the problem studied in [11]. The description of the motion of this pendulum model is given. The equations of motion are obtained in terms of two degrees of freedom. New conditions for the motion are considered when the length of the string and radii of the ellipse are sufficiently small comparing to the dimensions of the solid body. A large parameter is constructed such that . Using the large parameter technique, we solve this problem in terms of two degrees of freedom and . The large parameter procedure gives us the chance to study the motion in new conditions and a new domain of the solutions. The obtained analytical solutions are represented graphically through a computer program. On the other hand, the fourth-order Runge-Kutta numerical technique is used to show the validity of the solutions. The obtained phasing diagrams prove full agreement of the obtained analytical and numerical solutions. We deduce also that the change of the parameters of the solid body affects the behavior of the obtained solutions. For example, the increase in for the solutions of and gives a decrease in the number and amplitude of the waves and vice versa, while the increase of for the solutions and gives a decrease in the amplitudes of these solutions with the same number of the waves. When , we have a straightened wave of the motion and the number of the waves is infinite.
Data Availability
Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.
Conflicts of Interest
The author declares that he has no competing interests.