# The Dynamical Behavior of a Rigid Body Relative Equilibrium Position

Authors: T. S. Amer
Journal: Advances in Mathematical Physics
 (2017)

---

## Abstract

In this paper, we will focus on the dynamical behavior of a rigid body suspended on an elastic spring as a pendulum model with three degrees of freedom. It is assumed that the body moves in a rotating vertical plane uniformly with an arbitrary angular velocity. The relative periodic motions of this model are considered. The governing equations of motion are obtained using Lagrange’s equations and represent a nonlinear system of second-order differential equations that can be solved in terms of generalized coordinates. The numerical solutions are investigated using the fourth-order Runge-Kutta algorithms through Matlab packages. These solutions are represented graphically in order to describe and discuss the behavior of the body at any instant for different values of the physical parameters of the body. The obtained results have been discussed and compared with some previous published works. Some concluding remarks have been presented at the end of this work. The importance of this work is due to its numerous applications in life such as the vibrations that occur in buildings and structures.


---

## Body

## 1. Introduction

The pendulum models have provided the researchers with a fertile source of examples in nonlinear dynamics and lately in nonlinear control. The most famous rigid pendulum consists of a mass particle that is attached to one end of a massless rigid arm and the other end of the arm is fixed to a pivot point that provides a rotational joint for the arm and mass particle. If the arm and mass particle are constrained to move within a fixed plane, the system is referred to as a planar one-dimension pendulum. If the arm and mass particle are unconstrained, the system is referred to as a spherical two-dimension pendulum. The three-dimensional motion of the swinging spring is studied in [1]. The resonance phenomenon that occurs during the motion is also investigated. Nonlinear normal vibration modes of the spring pendulum and the system containing a pendulum absorber are considered in [2].The published articles on such models are very large. Few researches view the pendulum as a rigid body. Slandered pendulum models are defined by a single rotational degree of freedom, referred to as a planar rigid rotational degree of freedom, or two rotational degrees of freedom, referred to as spherical rigid pendulum. Control problems for planar and spherical pendulum models have been studied by outstanding researchers; see [3–6].In [7], the motion of a variable length pendulum was studied to determine the characteristics of motion. In [8], the process analysis method is presented as the analytical method to obtain a second-order approximate solution for a simple pendulum. This method does not depend on small parameter and therefore can overcome the disadvantages and limitations of the perturbation expansion method. In [9], the authors studied a simple pendulum with a hinge exhibiting bilinear hysteretic moment-rotation characteristics subjected to periodic base motions. They have shown that the bilinear hysteretic nature of the system becomes an effective way to limit the growth of the response during parametric resonance. Various perturbation techniques [10] were employed to obtain the analytical solutions for many physical problems. In [11], the authors studied the motion of the supported point of a pendulum on an ellipse and the method of small parameter [10] was used to obtain the periodic solution of the equation of motion. The relative periodic motion of a pendulum with an elastic string was studied in [12] and generalized in [13] when the motion of the supported point of a rigid body pendulum with an elastic string moves on an elliptical path. The equations of motion were deduced using Lagrange’s equations and solved through the small parameter method to obtain their solutions up to the second order of approximation.In [14] Amer and Bek studied the chaotic responses of a harmonically excited spring pendulum which moves in a circular path under some conditions. The obtained equations of motion represent a nonautonomous system of two nonlinear differential equations of two degrees of freedom. The approximate solution was obtained up to the third order using the multiple scales method [10]. The parametric control of oscillations and rotations of a compound pendulum was studied in [15]. An approximate asymptotic approach of this problem, based on a combination of the averaging method [10] and the maximum principle, is proposed and applied. The limiting cases of small oscillations and rapid rotation of a pendulum are studied in [16].In [17], the authors studied the vibration and stability of the nonlinear spring pendulum to describe the motion of a ship. The effects of the longitudinal absorber on the system are described through the obtained results. This model is modified in [18] by connecting the spring pendulum to the transverse absorber. So, the motion has three degrees of freedom under multiparametric excitations. The approximate solution is obtained using the multiple scales method up to the second-order approximations.The nonlinear two-degree-of-freedom system has been examined in [19]. The analytical approximate solution up to the third order is obtained using the same previous method. All the possible resonances of this solution are examined.The aim of this work is to investigate the motion of a rigid body suspended on an elastic massless spring. The equations of motion are derived using Lagrange’s equation and are considered as a nonlinear system of second-order differential equations. Each equation of this system depends on all the body variables with their derivatives. So, it is not easy to separate these equations as explicit second-order differential equations of one variable in one side. In order to overcome this quandary, the Mathematica program is used. Consequently, the numerical solutions are achieved using the fourth-order Runge-Kutta procedure of ode45 solver [20] with the aid of more recent computer package, for example, Matlab program. Computer codes are carried out to obtain the graphical representations of the attained numerical solutions for the different parameters of the body. The stability of the solutions is checked during the time interval of motion. Discussion of the results is presented through the comparison between the different plots for different variables. The importance of this problem is due to its wide applications in many fields such as physics and engineering applications like swaying buildings.


## 2. Formulation of the Problem

This section is devoted to introduce the motion of a rigid body suspended on an elastic massless spring as a pendulum model. So, we considerOXY as a coordinates system, rotating with angular velocity ω with respect to the downward axis OY, relative to the motion of a rigid body of mass m. The elastic spring is suspended to a point O1, where OO1=hsin⁡ωt at any time t. Let us suppose that the rigid body is attached with the spring at the point O2 and has the point C as a center of mass, φ1 denotes the deformation angle between the spring and the vertical axis O1Y1, and φ2 refers to the angle between the straight line directed through O2 to C and the vertical. Choosing an orthogonal coordinates system Cξηζ of the body, in which Cη directed along O2C, Cξ is perpendicular to Cη and lying in the OXY plane, while Cζ is perpendicular to the OXY plane (see Figure 1). Assume, without loss of generality, that the axes Cξ,Cη, and Cζ are the principal axes of inertia of the body.Figure 1
The rigid body pendulum.The coordinates of the center of mass of the bodyxC and yC, relative to the system OXY, can be written as(1)xC=hsin⁡ωt+ρsin⁡φ1+asin⁡φ2,yC=ρcos⁡φ1+acos⁡φ2,a=O2C,where ρ is the length of elastic string after time t.The kinetic energyT and the potential energy V of the system have the form(2)T=12mhωcos⁡ωt2+ρ˙2+ρφ˙12+aφ˙22+2hωρ˙sin⁡φ1cos⁡ωt+2hωρφ˙1cos⁡ωtcos⁡φ1+2ahωφ˙2cos⁡ωtcos⁡φ2+2aρ˙φ˙2sin⁡φ1-φ2+2aφ˙1φ˙2cos⁡φ1-φ2+ω2hsin⁡ωt+ρsin⁡φ1+asin⁡φ22+J3mφ˙22+ω2mJ1sin2φ2+J2cos2φ2,V=12k2ρ-l2-mgρcos⁡φ1+acos⁡φ2,where J1,J2, and J3 are the principal moments of inertia with respect to the axes Cξ,Cη, and Cζ, respectively, k2 is the spring’s constant, l represents the unstretched length of the string, and g denotes the gravitational attraction.According to the above equations, one can obtain the Lagrangian of the system [21](3)L=T-V.An inspection of (1)–(3), we can observe that the Lagrangian L is expressed in terms of three generalized coordinates ρ,φ1,φ2 and three corresponding generalized velocities ρ˙,φ˙1,φ˙2.Use the following Lagrange’s equations(4)ddt∂L∂q˙i-∂L∂qi=0;qi≡ρ,φ1,φ2,q˙i≡ρ˙,φ˙1,φ˙2to obtain the equations of motion in the form(5)ρ¨+aφ¨2sin⁡φ1-φ2-aφ˙22cos⁡φ1-φ2-ρφ˙12-ω22hsin⁡ωt+ρsin⁡φ1+asin⁡φ2sin⁡φ1-gcos⁡φ1+K2ρ-l=0,ρφ¨1+2ρ˙φ˙1+aφ¨2cos⁡φ1-φ2+aφ˙22sin⁡φ1-φ2-ω22hsin⁡ωt+ρsin⁡φ1+asin⁡φ2cos⁡φ1+gsin⁡φ1=0,l1φ¨2+ρ¨-ρφ˙12sin⁡φ1-φ2+ρφ¨1+2ρ˙φ˙1cos⁡φ1-φ2-ω22hsin⁡ωt+ρsin⁡φ1+asin⁡φ2cos⁡φ1-φ2+ω22mJ2-J1sin⁡2φ2+gsin⁡φ2=0,where(6)K2=k2m,l1=1aa2-J3m.Here l1 is the derived length of the body relative to O2. Equations (5) are the governing equations of motion of our model that represent a nonlinear system of second-order differential equations.In order to study this problem we consider that the oscillations of our system are closing to the position of the relative equilibrium. So, we can assume(7)J1=J2.Hence, for the relative equilibrium state the angles φ10 and φ20 are equal, and then we can write(8)ρ=b+ξt,φ1=φ0+φt,φ2=φ0+ψt,where φ0 represents the value of φ10 and φ20 and b denotes the pendulum string’s length in the case of relative equilibrium. Moreover, the quantities φ0 and b can be determined from the following equations:(9)K2b-l=ω2a+bsin2φ0+gcos⁡φ0,g=ω2a+bcos⁡φ0.Substituting from (8) into (5), then using (7) and (9), we get(10)ξ¨+a11ξ+a12φ+a13ψ=f1,bφ¨+aψ¨+b11ξ+b12φ+b13ψ=f2,l1ψ¨+bφ¨+c11ξ+c12φ+c13ψ=f3,where(11)a11=K2-ω2sin2φ0+2hcos⁡φ0sin⁡ωt,a12=bc11,a13=ac11,b11=c11,b12=K2b-l-ω2bcos2φ0,b13=-ω2acos2φ0,c11=-ω2sin⁡φ0cos⁡φ0,c12=-ω2bcos2φ0,c13=K2b-l+ω22hsin⁡φ0sin⁡ωt-acos2φ0,(12)f1=ξ+bφ˙2+Aξφ+aψ¨ψ-φ+aψ˙2+Bφψ+C1,f2=-ξφ¨-2ξ˙φ˙+ψ-φaψ˙+c11ξ+bφ2+Dξ+2a13ψφ+C2,f3=ξ¨ψ-φ+bφ˙2+ξφ˙2φ-ψ-ξφ¨-2ξ˙φ˙+2ω2hsin⁡ωtcos⁡φ0-b+a1-ψ2-b+ξφψc11-gsin⁡φ0+ω2ξφcos2φ0-ψsin2φ0,(13)A=-2c11,B=aω2cos2φ0,D=ω2cos2φ0,C1=-bK2-ω2sin2φ0+ω22hsin⁡φ0sin⁡ωt+asin2ωt+sin2φ0+gcos⁡φ0+K2l,C2=2hω2sin⁡ωt-b+ac11-gsin⁡φ0.Our principle aim is to obtain the numerical solutions of system (10) which consists of three nonlinear differential equations of second-order. In view of the right hand sides of these equations, we found three functions f1,f2, and f3 given by (12). In fact, it is not easy to obtain the second derivatives of the generalized coordinates ξ,φ, and ψ such that each equation contains one of these derivatives only.


## 3. Numerical Solutions

This section is devoted to discuss the numerical solutions for the considered model in Section2. Computer programs are carried out to investigate the graphical representations for these solutions, to describe the motion, and to illustrate the behavior of the pendulum at any time.System (10) consists of three nonlinear differential equations of second-order in terms of ξ,φ, and ψ and is reconsidered to obtain the numerical solutions in framework of the fourth-order Runge-Kutta algorithms through Matlab packages [22]. Each equation of this system includes all variables ξ,φ,ψ and their derivatives from the first and second order; see systems of (10), (11), (12), and (13). So, the mentioned system is more complicated to deal with and to obtain another corresponding one consisting of second-order differential equations in terms of ξ¨,φ¨, and ψ¨ explicitly. Computer codes are utilized in order to overcome these difficulties and to separate each of ξ¨,φ¨, and ψ¨. Consequently, system (10) is transformed into the following system with the aid of (11), (12), and (13):(14)ξ¨=-1H-aψ-b+ξc11ξ+c12φ+c13ψ-f3+2ξ˙φ˙-b+ξφ-ψφ˙2+bb11ξ+b12φ+b13ψ+gsin⁡φ0+ω2·-ξφcos2φ0-2hsin⁡ωtcos⁡φ0+-b+b+ξφψ+aψ2-1·cos⁡φ0sin⁡φ0+ξφsin2φ0+2ξ˙φ˙-b+ξφ-ψφ˙2+ab-b+ξl1C1-a11ξ-a12φ-a13ψ+b+ξφ˙2+Aξ+Bψ+aψ˙2φ,φ¨=-1bf3-c11ξ-c12φ-c13ψ+l1aψ·C1-a11ξ-a12φ-a13ψ+b+ξφ˙2+Aξ+Bψ+aψ˙2φ+1aHψl1-aψb+ξ+2ξ˙φ˙-b+ξφ-ψφ˙2b+ξφ˙2+Aξ+Bψ+aψ˙2φb+ξφ˙2+Aξ+Bψ+aψ˙2φ·f3-c11ξ-c12φ-c13ψ+b·b11ξ+b12φ+b13ψ+gsin⁡φ0+ω2·-ξφcos2φ0-2hsin⁡ωtcos⁡φ0+b+ξφψ+aψ2-1-bcos⁡φ0sin⁡φ0+ξφsin2φ0+2ξ˙φ˙-b+ξφ-ψφ˙2+ab-b+ξl1C1-a11ξ-a12φ-a13ψAξ+Bψ+aψ˙2+b+ξφ˙2+Aξ+Bψ+aψ˙2φ,ψ¨=-1Hbφ-ψa11ξ+a12φ+a13ψ-Aξ+Bψφ2ξ˙φ˙+φω2ξsin2φ0+aφ-ψψ˙2-bb11ξ+b12φ+b13ψ+b+ξc11ξ+c12φ+c13ψ-bC1φ-ψ+bω2cos⁡φ0ξφcos⁡φ0+2hsin⁡ωt-b+ξf3-bgsin⁡φ0+bω2·a+b-ψb+ξφ+aψcos⁡φ0sin⁡φ0-b2ξ˙φ˙+φω2ξsin2φ0+aφ-ψψ˙2,where (15)H=abψ2-φψ-1+b+ξl1.It is clear that the left hand sides of the equations of the previous system are given explicitly in terms ofξ¨,φ¨, and ψ¨, respectively. On the other hand, the right hand sides are functions of ξ,φ,ψ,ξ˙,φ˙, and ψ˙.The ode45 solver is used in order to obtain the numerical solutions of the nonstiff ordinary differential equations of the previous system (14), in which this solver uses a variable step of Runge-Kutta technique [20]. So we can rewrite system (14) as a system of coupled first-order differential equations as follows.A choice of the state variables for this system is(16)X1=ξ,X2=φ,X3=ψ,X4=ξ˙,X5=φ˙,X6=ψ˙,which results in the following state-equations:(17)X˙1=X4,X˙2=X5,X˙3=X6.Use (16) and (17) into system (14) to get(18)X˙4=-1H-aX3-b+X1+2X˙1X˙2-b+X1X2-X3X˙22AX1+BX3+aX˙32·c11X1+c12X2+c13X3-f3+bb11X1+b12X2+b13X3+gsin⁡φ0+ω2·-X1X2cos2φ0-2hsin⁡ωtcos⁡φ0+-b+b+X1X2X3+aψ2-1cos⁡φ0sin⁡φ0+X1X2sin2φ0+2X˙1X˙2-b+X1X2-X3X˙22+ab-b+ξl1C1-a11X1-a12X2AX1+BX3+aX˙32-a13X3+b+X1X˙22+AX1+BX3+aX˙32X2,X˙5=-1bf3-c11X1-c12X2-c13X3+l1aX3C1-a11X1-a12X2-a13X3+b+X1X˙22+AX1+BX3+aX˙32X2+1aHX3·l1-aX3b+X1f3-c11X1-c12X2-c13X3+2X˙1X˙2-b+X1X2-X3X˙22+bb11X1+b12X2+b13X3+gsin⁡φ0+ω2-X1X2cos2φ0-2hsin⁡ωtcos⁡φ0+b+X1X2X3+aX32-1-bcos⁡φ0sin⁡φ0+X1X2sin2φ0+2X˙1X˙2-b+X1X2-X3X˙22+ab-b+X1l1·C1-a11ξ-a12X2-a13X3AX1+BX3+aX˙32+b+X1X˙22+AX1+BX3+aX˙32X2,X˙6=-1HbX2-X3ω2X1sin2φ0+aX2-X3X˙32·a11X1+a12X2+a13X3-AX1+BX3X2-bb11X1+b12X2+b13X3+b+X1c11X1+c12X2+c13X3-bC1X2-X3+bω2cos⁡φ0X1X2cos⁡φ0+2hsin⁡ωt-b+X1f3-bgsin⁡φ0+bω2a+b-X3b+X1X2+aX3·cos⁡φ0sin⁡φ0-b2X˙1X˙2+X2ω2X1sin2φ0+aX2-X3X˙32·ω2X1sin2φ0+aX2-X3X˙32.The following initial conditions are required to achieve the numerical solution of (18), by using the fourth-order Runge-Kutta method of ode45 solver in framework of Matlab program,(19)X10=0.001,X20=0.1,X30=0.01,X˙40=0,X˙50=0,X˙60=0in addition to the following physical parameters of the considered model:(20)m=50kg,g=9.8m·s-2,J1=8kg·m2,l=0.7m,ω=4rad·s-1,a=0.5m,b=0,3m,φ0=0,0.4radh=2.5,4.5,t=0⟶1.7min.Figure 2 shows the variation of the solutions ξ,φ,ψ and their derivatives ξ˙,φ˙,ψ˙ against time t when φ0=0 and φ0=0.4rad. This figure is drawn at b=3m, h=2.4, and a=0.5m. The variations of ξ,φ, and ψ with ξ˙,φ˙, and ψ˙, respectively, are illustrated in Figure 3, namely, the phase plane diagrams that are represented in Figures 3(a), 3(b), 3(c) and 3(d), 3(e), 3(f) when φ0=0 and φ0=0.4rad, respectively, with the same other parameters that are taken into consideration in Figure 2.Figure 2
Variation of the solutions and their derivatives versus timet when b=3m, h=4.5, and a=0.5m: (a, d) show the effect of t on the behavior of ξ and ξ˙ waves when φ0=0 and φ0=0.4rad, respectively, (b, e) show the effect of t on the behavior of φ and φ˙ waves when φ0=0 and φ0=0.4rad, respectively, and (c, f) show the effect of t on the wave that describes the behavior of ψ and ψ˙ when φ0=0 and φ0=0.4rad, respectively.
(a)
(b)
(c)
(d)
(e)
(f)Figure 3
The phase plane diagram whenb=3m, h=4.5, and a=0.5m: (a, d) represent the variation of the amplitude ξ with its velocity ξ˙ at φ0=0 and φ0=0.4rad, respectively, (b, e) represent the variation of the amplitude φ with its velocity φ˙ at φ0=0 and φ0=0.4rad, respectively, and (c, f) represent the variation of the amplitude ψ with its velocity ψ˙ at φ0=0 and φ0=0.4rad, respectively.
(a)
(b)
(c)
(d)
(e)
(f)In these figures, our principle aim is to investigate the effect of increasing time on the motion of pendulum.According to the calculations depicted in Figure2(a), we found that when φ0=0, the wave of the elongation ξ grows up with the increasing of time till t=0.9min. After that both of the elongation ξ and its derivative ξ˙ fluctuate between increasing and decreasing when time reaches t=1.43min. Thus the wave of the solution ξ is stable; see the phase plane Figure 3(a). With the passing of time, one can observe that ξ and ξ˙ are growing quickly, so the motion will be unstable after t=1.43min. The rage behavior of both ξ and ξ˙ is due to the weight of the rigid body and the values of the principal moments of inertia. Consequently, we expect that behavior of elongation becomes greater as observed in Figures 2(a) and 2(d). Moreover, the variation of the spring between stretching and contraction is consistent with the phase plane diagrams represented in Figures 3(a) and 3(d).It is worthwhile to notice from Figure2(b) that when time t increases from t=0 to t=0.4min, the behavior of the angle φ increases gradually to reach the value φ≃1rad≃57° and then decreases slowly to reach φ≃0.8rad≃46° during the time period t∈]0.4,0.9[min. After t=0.9min, the decline of the wave becomes quickly to reach φ≃-2.3rad≃-132° at the end of time period (minus sign indicates opposite direction). This is not possible because φ must belong to the interval -π/2,π/2. So, the motion of the wave is unstable as it is manifest from Figure 3(b). On the other hand, φ˙ increases till t≃0.2min and then fluctuates as indicated from Figure 2(b).The graphs displayed in Figures2(c) and 3(c) describe the variation of the (ψ and ψ˙) against time and the phase plane diagram (ψ˙ with ψ), respectively, when φ0=0. It is clear that when time belongs to the period 0,0.43min, the angle ψ remains stationary, and then its wave oscillates between decreasing and increasing till t=1.43min. After that time, the angle ψ increases up to the end of time interval and consequently the motion will be stable as seen from Figure 3(c) during the period 0<t≤1.43. It is obvious from Figure 2(c) that the behavior of ψ˙ remains stationary to some extent through the time interval 0,0.6min and then oscillates between increasing and decreasing till t=1.7min.It should be noticed that whenφ0=0.4rad, the stretching on the string ξ increases gradually till the time t becomes 0.9min, and then ξ and ξ˙ oscillate between increasing and decreasing when the time reaches the end of time interval; see Figure 2(d). Consequently, the wave of the solution is stable as seen from the phase plane Figure 3(d).An inspection of the graphs depicted in Figure2(e) shows that the wave describing the behavior of the angle φ increases gradually from φ=0 at t=0 to its maximum value φ≃0.9rad≃56° at t=0.4min and then decreases slowly at t≃1min to reach its minimum value φ≃-0.19rad≃-11° (minus sign indicates opposite direction) at t≃1.26min. With the increasing of time, the wave grows again to reach the value φ≃0.85rad≃49° at t≃1.39min. Thus, the motion is stable as it is manifest from Figure 3(e). On the other hand, φ˙ increases and decreases as indicated from Figure 2(e).Also, it is remarkable from Figure2(f) that the behavior of the angle ψ remains steady till t=0.5min; then its wave oscillates between decreasing and increasing. Consequently the motion will be stable as seen from Figure 3(f). It is worthwhile to notice also from Figure 2(f) that the behavior of ψ˙ oscillates between increasing and decreasing.From the above observations, we can conclude that the motion of our model is more stable whenφ0=0.4rad than when φ0=0. This highlights the importance of the effect of φ0 value on the motion. It is worthwhile to notice that the comparison between the solutions ξ,φ, and ψ included in Figures 2(a), 2(b), and 2(c) with the corresponding Figures 2(d), 2(e), and 2(f) reveals that the amplitude of the waves decreases when φ0 increases from 0 to 0.4rad. On the other hand, the comparison between their derivatives shows that the amplitude of the waves increases when φ0 increases. Figure 4 shows the variation of (ξ,ξ˙),(φ,φ˙), and (ψ,ψ˙) with time t when φ0 changes from 0 for Figures 4(a), 4(b), and 4(c) to 0.4rad for Figures 4(d), 4(e), and 4(f) at the same values of other parameters b=3m,h=4.5, and a=0.5m. According to the calculations depicted in these figures, we can consider these figures as a rotation of the corresponding parts of Figure 3 with time to observe the bending and crossing of the resulting curves.Figure 4
The 3D pattern whenb=3m, h=4.5, and a=0.5m: (a, d) indicate the variation of ξ and ξ˙ versus t when φ0=0 and φ0=0.4rad, respectively, (b, e) indicate the variation of φ and φ˙ versus t when φ0=0 and φ0=0.4rad, respectively, and (c, f) indicate the variation of ψ and ψ˙ versus t when φ0=0 and φ0=0.4rad, respectively.
(a)
(b)
(c)
(d)
(e)
(f)Figures5(a), 5(c) and 5(b), 5(d) represent 3D plots that illustrate the variation of the solutions ξ,φ via ψ and ξ˙,φ˙ via ψ˙, respectively, for different values of φ0 when b=3m,h=4.5, and a=0.5m. The graphs displayed in parts of Figure 6 show the variation of (ξ,ξ˙),(φ,φ˙), and (ψ,ψ˙) against time t when b=0 with consideration of the parameters φ0=0.4rad,h=4.5, and a=0.5m. The corresponding phase plane between the amplitudes ξ,φ,ψ and their derivatives ξ˙,φ˙,ψ˙ is represented in parts of Figure 7. Inspection of the graph depicted in Figure 6(a) shows that when time t increases from 0 to 0.45min the behavior of the solution ξ remains stationary and quickly growing during the time interval t∈0.45,1.05min and then oscillates till the end of time interval. This indicates that the motion is stable as seen from Figure 7(a). On the other side, the behavior of the derivative ξ˙ remains approximately stationary during the interval t∈[0,0.45]min and then fluctuates with the increasing of time; see Figure 6(a).Figure 5
The 3D diagrams whenb=3m, h=4.5, and a=0.5m: (a, c) elucidate the variation of ξ and φ versus ψ when φ0=0 and φ0=0.4rad, respectively, and (b, d) elucidate the variation of ξ˙ and φ˙ versus ψ˙ when φ0=0 and φ0=0.4rad, respectively.
(a)
(b)
(c)
(d)Figure 6
(a), (b), and (c) explain the variation of the solutionsξ,φ, and ψ with their derivatives ξ˙,φ˙, and ψ˙ via time t, respectively, when b=0, h=4.5, a=0.5m, and φ0=0.4rad.
(a)
(b)
(c)Figure 7
The phase plane diagrams between amplitudes and their velocities atb=0, h=4.5, a=0.5m, and φ0=0.4rad: (a) shows the influence of ξ on ξ˙, (b) shows the effect of φ on φ˙, and (c) shows the variation of ψ with ψ˙.
(a)
(b)
(c)By the same way, we can observe that the wave of the angleφ increases through a short time to reach its maximum value φ≃2.7rad≃155° at t≃0.23min, taking into consideration that -π/2<φ<π/2, and then decreases slowly to reach its minimum value φ≃-0.9rad≃-52° at the end of time interval; see Figure 6(b). This indicates that the motion is close to be stable as observed from the phase plane Figure 8(b). As seen from Figure 6(b), φ˙ increases and decreases quickly during the period t∈0,0.1min to reach its minimum value at the end of time interval.Figure 8
The 3D plots atb=0, h=4.5, a=0.5m, and φ0=0.4rad: (a) illustrates the variation of ξ and ξ˙ via t, (b) illustrates the variation of φ and φ˙ via t, and (c) illustrates the variation of ψ and ψ˙ via t.
(a)
(b)
(c)The variation ofψ and ψ˙ with time is illustrated in Figure 6(c). In this figure our main goal is to examine the influence of time on the motion of pendulum. It is clear that the behavior of ψ and ψ˙ remains stationary (to some extent) when t∈0,0.5min; then their waves fluctuate till the end of time interval. Consequently the motion is stable as seen from the phase plane diagram, Figure 7(c). The comparison between parts of Figure 6 with the corresponding Figures 2(d), 2(e), and 2(f) shows that when b changes from 0 to 3m the amplitude of the waves decreases. Also, the motion will be more stable when b=3m than when b=0 as seen from the corresponding phase plane diagrams, that is, Figures 3(d), 3(e), 3(f) and 7(a), 7(b), 7(c), respectively.On the other hand, parts of Figure8 show 3D plots that describe the variation of the solutions and their derivative via time when b=0, h=4.5, φ0=0.4rad, and a=0.5m. The plots displayed in the mentioned parts show bending and crossing of the resulting curves.Figures9(a), 9(b), and 9(c) show the variation of the solutions ξ,φ,ψ and their derivatives ξ˙,φ˙,ψ˙ with time t when h=2.5 for the given values of other parameters b=3m, φ0=0.4rad, and a=0.5m. In view of the first part, we can conclude that when time t increases, each of the waves ξ and ξ˙ oscillates between increasing and decreasing till t=1.46min and then increases gradually. So, the motion is stable as seen from Figure 10(a).Figure 9
(a), (b), and (c) demonstrate the variation of (ξ and ξ˙), (φ and φ˙), and (ψ and ψ˙) against time t, respectively, at b=3m, h=2.5, a=0.5m, and φ0=0.4rad.
(a)
(b)
(c)Figure 10
The phase plane diagrams which portray the relation between amplitudes and their velocities atb=3m, h=2.5, a=0.5m, and φ0=0.4rad: (a) describes the influence of ξ on ξ˙, (b) shows the effect of φ on φ˙, and (c) illustrates the variation of ψ with ψ˙.
(a)
(b)
(c)From a closer look on the second part of  Figure9(b) we can write, with the increasing of time, the behavior of φ wave increases to reach its maximum value φ≃0.9rad≃52° at t=0.43min and then decreases slowly through the period t∈0.43,1.19min. After that, its behavior has a sharp decline in a few seconds (about 2.4 s) and then increases till the end of time period and consequently the motion is stable; see Figure 10(b).According to the calculations depicted in Figure9(c), we can observe that the waves describing ψ and ψ˙ decrease slowly till t=0.9min and then increase and decline sharp. The phase plane Figure 10(c) shows that the behavior of ψ is not stable.When parts of Figure9 and their phase plane parts (of Figure 10) are generally compared with the corresponding Figures 2(d), 2(e), and 2(f) and their phase plane Figures 3(d), 3(e), and 3(f), we can observe that amplitude of the wave increases when h=4.5 compared to when h=2.5 and the motion is more stable when h=4.5. An inspection of parts of Figure 11 reveals the 3D plots when h=2.5 with the same other data considered in Figures 9 and 10. Figure 10 shows the variation of the solutions ξ,φ,ψ and their derivatives ξ˙,φ˙,ψ˙ with time t. It is worthwhile to notice that the comparison between Figures 4(d), 4(e), and 4(f) and Figures 11(a), 11(b), and 11(c) shows more bending and crossing of the curves in Figures 4(d), 4(e), and 4(f) when h=4.5 than the corresponding ones of Figure 11.Figure 11
The 3D patterns atb=3m, h=2.5, a=0.5m, and φ0=0.4rad: (a) illustrates the variation of ξ and ξ˙ versus t, (b) illustrates the variation of φ and φ˙ versus t, (c) illustrates the variation of ψ and ψ˙ versus t.
(a)
(b)
(c)Now, we study the last case whenh=0 with the same other data b=3m, φ0=0.4rad, and a=0.5m. The obtained results are represented graphically in Figures 12(a), 12(b), and 12(c), while their phase plane diagrams are given in Figures 12(d), 12(e), and 12(f). At the first glance, we can conclude that this case is not stable, so it is very important to notice that the dimensionless parameter h must take any value different from zero as it is pointed in Figure 2 (h=4.5) and Figure 9 (h=2.5). This elucidates the importance of h parameter on the motion.Figure 12
(a, b, and c) explain the variation of the solutionsξ,φ, and ψ with their derivatives ξ˙,φ˙, and ψ˙ via time t, respectively, when b=3m, h=0, a=0.5m, and φ0=0.4rad; (d, e, and f) illustrate the variation of the solutions against their first derivatives for the same values of the considered parameters.
(a)
(b)
(c)
(d)
(e)
(f)


## 4. Conclusion

A conclusion that may be made here is that the problem of the relative motion of a rigid body as a pendulum model is investigated. The governing deferential equations are obtained using Lagrange’s equations. Mathematica package was utilized, in order to overcome the difficulties that appear in the separation of the second derivatives of the generalized coordinatesξ,φ, and ψ for the nonlinear system (10). Computer codes are used to obtain the numerical solutions for system (14). These solutions are represented graphically using Matlab program, to study the influence of the different parameters on the motion. The good effect of the parameters h,b, and φ0 on the motion is obvious from the mentioned plots. The motion of our model is more stable when the parameters h,b, and φ0 take values run away from zero. This highlights the importance of the effect of these parameters on the motion. Such results have been confirmed by many works, such as Ismail [13] and Amer and Bek [14].



---

*Source: 8070525-2017-03-08.xml*