Next Article in Journal
Reaction Center of Rhodobacter Sphaeroides, a Photoactive Protein for pH Sensing: A Theoretical Investigation of Charge Transport Properties
Next Article in Special Issue
ALIPPF-Controller to Stabilize the Unstable Motion and Eliminate the Non-Linear Oscillations of the Rotor Electro-Magnetic Suspension System
Previous Article in Journal
Alternative Methods of Bioactive Compounds and Oils Extraction from Berry Fruit By-Products—A Review
Previous Article in Special Issue
Nonlinear Analyses of Porous Functionally Graded Sandwich Piezoelectric Nano-Energy Harvesters under Compressive Axial Loading
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Stability Analysis of a Vibrating Auto-Parametric Dynamical System Near Resonance

1
Mathematics Department, Faculty of Science, Tanta University, Tanta 31527, Egypt
2
Institute of Applied Mechanics, Poznan University of Technology, 60-965 Poznan, Poland
3
Department of Physics and Engineering Mathematics, Faculty of Engineering, Tanta University, Tanta 31734, Egypt
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(3), 1737; https://doi.org/10.3390/app12031737
Submission received: 19 January 2022 / Revised: 2 February 2022 / Accepted: 4 February 2022 / Published: 8 February 2022
(This article belongs to the Special Issue Nonlinear Vibrations)

Abstract

:
This paper examines a new vibrating dynamical motion of a novel auto-parametric system with three degrees of freedom. It consists of a damped Duffing oscillator as a primary system attached to a damped spring pendulum as a secondary system. Lagrange’s equations are utilized to acquire the equations of motion according to the number of the system’s generalized coordinates. The perturbation technique of multiple scales is applied to provide the solutions to these equations up to a higher order of approximations, with the aim of obtaining more accurate novel results. The categorizations of resonance cases are presented, in which the case of primary external resonance is examined to demonstrate the conditions of solvability of the steady-state solutions and the equations of modulation. The time histories of the achieved solutions, the resonance curves in terms of the modified amplitudes and phases, and the regions of stability are outlined for various parameters of the considered system. The non-linear stability, in view of both the attained stable fixed points and the criterion of Routh–Hurwitz, is investigated. The results of this paper will be of interest for specialized research that deals with the vibration of swaying buildings and the reduction in the vibration of rotor dynamics, as well as studies in the fields of mechanics and space engineering.

1. Introduction

The dynamical motion of a pendulum in generality or the motion of any system containing a pendulum is considered as one of the oldest scientific subjects in non-linear dynamics. The planar motion of a simple pendulum has been investigated in several research works that reveal its complicated behavior, e.g., [1,2,3,4,5]. The periodic motion of a simple pendulum with a pivot point that moved on a vertical ellipse was investigated in [1], in which some special cases were examined to show the different motions of this point. On the other hand, the chaotic motion of a parametrically excited pendulum as a mathematical model has been investigated numerically, experimentally, and analytically in [2,3,4], respectively. The harmonic motion of the non-linear vibrations of an excited auto-parametric system was studied in [5]. The harmonic balance method (HBM) [6] was used to evaluate the response of this system in which its performance was studied.
The periodic solutions of the equations of motion (EOM) of a greatly damped pendulum were investigated in [7] under some certain constraints, in which the behavior of the pendulum led to chaos through a symmetrical pitchfork bifurcation. A global analytical inspection of the pendulum is presented in [8,9,10], where the pendulum was subjected to various applied external forces and moments. The dynamics of two planar flexible pendulums attached to an excited horizontal platform were studied in [11]. The authors asserted that the pendulum exhibited rotational and oscillatory motion simultaneously, in which the stability of the motion existed for the fixed in-phase and anti-phase states simultaneously.
The HBM was used in [12] to obtain the analytic results of an auto-parametric system connected to a non-linear elastic pendulum close to one of the arising parametric resonances. The authors verified the attained results by the comparison with the numerical results. On the other hand, the perturbation technique of multiple scales (PTMS) [6] has been used in several scientific works, such as [13,14,15,16,17,18,19,20,21], to obtain the approximate solutions of the controlling EOM of the corresponding dynamical systems. In refs. [13,14,15], the authors investigated the motion of an elastic pendulum under the influence of excitation force and moment, when its supported point moved in circular, Lissajous, and elliptic paths, respectively. The possible fixed points and the corresponding steady-state solutions of the stable points are gained. Some practical examples are presented in [15] to reveal the possible movements of the pendulum’s pivot point. The rigid body (RB) motion as a pendulum was investigated in [16] with a fixed pivot point, while its generalizations were examined in [17,18] for the motion of a damped spring RB with linear and non-linear stiffness, respectively, and when the movement of the supported pendulum’s point was constrained to elliptic trajectories. The obtained solutions were sketched with the possible resonances curves for different parameters of the investigated system. Another generalization is found in [19], where the path of the suspension point had the form of a Lissajous curve. The analytical results were compared with the numerical results, revealing a high consistency between them and showing the accuracy of the analytical technique.
Recently, the motion of a three degrees of freedom (DOF) double pendulum was examined in [20] under the existence of a harmonic excitation force and two moments. The first pendulum was considered to be rigid, in which its first point was constrained to move only in an elliptical path, while its second end was connected with a damped spring pendulum. The attained solutions were plotted and verified with the numerical solutions. On the other hand, the resonance curves corresponding to the stability and instability areas of a two DOF system were investigated in [21]. It was considered that the motion was constrained to a plane under the influence of two harmonic forces at the free end of a damped elastic pendulum and one harmonic moment that acted at the suspension point.
The motion of a spring-damper pendulum in the presence of a resistance force besides torque at an immovable suspension point and the force at the free end of the pendulum with harmonic magnitude were investigated in [22]. The effectiveness of the fluid flow properties, the buoyancy and drag forces, and an excitation force on the 2DOF spring-damper were examined in [23]. All resonance cases were characterized, in which the solutions at the steady state were verified in view of the conditions of solvability. The dynamical motion of a tuned absorber was studied in [24], in which the PTMS was applied to obtain the approximate solutions up to the third order for the case of a movable pivot point on an ellipse with a stationary angular velocity. The modulation equations were derived and reduced to two algebraic equations in order to investigate the stability of the fixed points according to selected data of the considered system.
Several studies [25,26,27,28,29,30,31] have investigated the Duffing pendulum systems about the primary and secondary resonances using analytical approaches or numerical analysis. The responses of an auto-parametric system containing a spring, mass, and damper was studied in [25] using the HBM. The approximate results of another auto-parametric system consisting of a damped pendulum connected to a non-linear vibrating oscillator near the region of principal resonance were investigated in [26]. In ref. [27], the author examined the oscillations of a parametrically self-excited 2DOF dynamical system consisting of coupled oscillators identified by self-excitation, with the non-linearity of a Duffing type. The vibrational analysis of a mechanical system similar to an auto-parametric system subjected to an external harmonic excitation was investigated in [28]. The authors reported that a semi-active magnetic damper and a non-linear spring acted on the system to improve its dynamics and control movements. The analytic solutions of this system were obtained in [29] using the HBM, and they were verified according to comparison with the numerical results. The discussion of the vibration mitigation and energy harvesting of the same system is found in [30], in which the concept of energy harvesting is closer to the vibration absorber concept. In ref. [31], the dynamics of a pendulum suspended on a forced Duffing oscillator were investigated using numerical analysis. The bifurcation analysis of the pendulum was also performed.
This paper focuses on the motion of a novel three DOF auto-parametric dynamical system consisting of a damped spring pendulum as the secondary system suspended with a forced non-linear Duffing damping oscillator as the primary system. The governing EOM are derived using the second kind of Lagrange’s equations according to the generalized coordinates of the considered system. The PTMS is used to solve the EOM analytically up to the third order of approximation and to obtain more accurate results. All cases of resonance are categorized, in which the case of primary external resonance is examined to obtain the conditions of solvability besides the modulation equations. The time histories of the modified amplitudes, the modified phases, the achieved results, resonance curves, and the path’s projections on selected space planes are graphed to describe the impact of the system’s parameters on the motion. The stability and instability regions in view of the investigated resonance case are presented using the Routh–Hurwitz criterion [32]. The results of this work are generalized as the previous results in [5,31] for the case of the motion of a rigid pendulum without any effectiveness of force and moment. These results have a significant impact in special dampers mounted in buildings, applied in the opposite directions of river vortices or earthquakes, and in the applications of pendulums that are mounted on bridge towers.

2. Formulation of the Dynamical System

Let us consider the dynamical vibrations of an auto-parametric system consisting of an excited Duffing oscillator with normal length 0 , static elongation Y s , stiffness k 1 , and mass M that is periodically forced to move vertically as a primary system. The secondary system is a damped spring pendulum with normal length 1 , static elongation X s , stiffness k 2 , and mass m . This spring is connected to the Duffing oscillator at a point O . Let Y ( t ) denote the elongation of the Duffing oscillator and ϕ represent the rotation angle of the spring that lies between the downward vertical Y 1 -axis and the direction of the spring’s elongation X ( t ) (see Figure 1).
Moreover, we consider that the motion is influenced by external linear excitation forces, F 1 ( t ) = F 1 cos ( Ω 1 t ) and F 3 ( t ) = F 3 cos ( Ω 3 t ) , and an excitation moment M ϕ ( t ) = F 2 cos ( Ω 2 t ) that are acting on the masses M ,   m and at the point O , respectively. Here, Ω 1 , Ω 2 , and Ω 3 are the frequencies of F 1 ( t ) , M ϕ ( t ) , and F 3 ( t ) , respectively. The investigated auto-parametric system is damped by viscous forces and moments C 1 Y ˙ , C 2 ϕ ˙ , and C 3 X ˙ , in which C 1 , C 2 , and C 3 are the respective coefficients of the viscous damping.
Based on the above description, we can write the potential energy V and kinetic energy T as follows:
V = 1 2 k 1 Y 2 + 1 2 k 2 X 2 M g ( l 0 + Y ) m g ( l 0 + Y ) m g ( l 1 + X ) cos ϕ , T = 1 2 ( M + m ) Y ˙ 2 + 1 2 m ( l 1 + X ) 2 ϕ ˙ 2 + 1 2 m X ˙ 2 m ( l 1 + X ) Y ˙ ϕ ˙ sin ϕ + m X ˙ Y ˙ cos ϕ ,
where g is the acceleration of Earth’s gravitational force and the dots are time derivatives.
The governing EOM can be derived using the following Lagrange’s Equations:
d d t ( L Y ˙ ) ( L Y ) = Q Y , d d t ( L ϕ ˙ ) ( L ϕ ) = Q ϕ , d d t ( L X ˙ ) ( L X ) = Q X .
Here, L = T V is the Lagrange’s function, ( Y , ϕ , X ) and ( Y ˙ , ϕ ˙ , X ˙ ) are the system’s generalized coordinates and their corresponding velocities, respectively, and ( Q Y , Q ϕ , Q X ) are the generalized forces which are formed as follows:
Q Y = F 1 cos ( Ω 1 t ) C 1 Y ˙ , Q ϕ = F 2 cos ( Ω 2 t ) C 2 ϕ ˙ , Q X = F 3 cos ( Ω 3 t ) C 3 X ˙ .
Consider the following definitions of the dimensionless parameters:
ω ¯ 1 2 = k 1 m + M , ω ¯ 2 2 = g l 1 , ω ¯ 3 2 = k 2 m , ω 1 2 = ω ¯ 1 2 ω ¯ 2 2 , ω 2 2 = ω ¯ 3 2 ω ¯ 2 2 , c 1 = C 1 ( m + M ) ω ¯ 2 , c 2 = C 2 m l 1 2 ω ¯ 2 , c 3 = C 3 m ω ¯ 2 , f 1 = F 1 ( m + M ) l 1 ω ¯ 2 2 , f 2 = F 2 m l 1 2 ω ¯ 2 2 , f 3 = F 3 m l 1 ω ¯ 2 2 , p 1 = Ω 1 ω ¯ 2 , p 2 = Ω 2 ω ¯ 2 , p 3 = Ω 3 ω ¯ 2 , β = m m + M , τ = ω 2 t , x = X X s l 1 , y = Y Y s l 1 , Y s = ( m + M ) g k 1 , X s = m g k 2 , l 0 = 0 + Y s , l 1 = 1 + X s .
The substitution of Equations (1), (3) and (4) into the system of Equation (2) produces the following dimensionless form of the EOM:
y ¨ ( τ ) + ω 1 2 y 2 β x ˙ ( τ ) ϕ ˙ ( τ ) sin ϕ β ( 1 + x ) ϕ ˙ 2 ( τ ) cos ϕ + β x ¨ ( τ ) cos ϕ β ( 1 + x ) ϕ ¨ ( τ ) sin ϕ + c 1 y ˙ ( τ ) 2 = f 1 cos ( Ω 1 τ ) ,
( 1 + x ) 2 ϕ ¨ ( τ ) ( 1 + x ) y ¨ ( τ ) sin ϕ + 2 x ˙ ( τ ) ϕ ˙ ( τ ) + 2 x x ˙ ( τ ) ϕ ˙ ( τ ) + ( 1 + x ) sin ϕ + c 2 ϕ ˙ ( τ ) = f 2 cos ( Ω 2 τ ) ,
x ¨ ( τ ) + ω 2 2 x + ( 1 + x ) ϕ ˙ 2 ( τ ) + y ¨ ( τ ) cos ϕ cos ϕ + c 3 x ˙ ( τ ) + 1 = f 3 cos ( Ω 3 τ ) .
The inspection of the system of Equations (5)–(7) shows that the governing EOM consists of three non-linear ordinary differential equations (ODE) in terms of y , ϕ , and x .

3. The Perturbation Technique

The main objective of the present section is to obtain the approximate analytic solutions of the EOM (5)–(7) using the PTMS to address the cases of resonance, obtain the conditions of solvability, and achieve the equations of the modulations. Therefore, we can approximate the functions sin ϕ and cos ϕ in Equations (5)–(7) using a Taylor series according to their first terms, for which these approximations are valid in a small neighborhood of the position of static equilibrium [33]. Then, Equations (5)–(7) become:
y ¨ ( τ ) + ω 1 2 y 2 β x ˙ ( τ ) ϕ ˙ ( τ ) ( ϕ ϕ 3 6 ) β ( 1 + x ) ϕ ˙ 2 ( τ ) ( 1 ϕ 2 2 ) + β x ¨ ( τ ) ( 1 ϕ 2 2 ) β ( 1 + x ) ϕ ¨ ( τ ) ( ϕ ϕ 3 6 ) + c 1 y ˙ ( τ ) 2 = f 1 cos ( Ω 1 τ ) ,
( 1 + x ) 2 ϕ ¨ ( τ ) ( 1 + x ) y ¨ ( τ ) ( ϕ ϕ 3 6 ) + 2 x ˙ ( τ ) ϕ ˙ ( τ ) + 2 x x ˙ ( τ ) ϕ ˙ ( τ ) + ( 1 + x ) ( ϕ ϕ 3 6 ) + c 2 ϕ ˙ ( τ ) = f 2 cos ( Ω 2 τ ) ,
x ¨ ( τ ) + ω 2 2 x + ( 1 + x ) ϕ ˙ 2 ( τ ) + ( 1 ϕ 2 2 ) y ¨ ( τ ) + ϕ 2 2 + c 3 x ˙ ( τ ) = f 3 cos ( Ω 3 τ ) .
Based on the motion in a small neighborhood of the static equilibrium position, we consider a small parameter 0 < ε < < 1 , in which the parameters y , ϕ , and x can be rewritten in terms of new variables, ξ , φ , and χ , respectively, according to the following forms:
y ( τ ) = ε ξ ( τ ; ε ) , ϕ ( τ ) = ε φ ( τ ; ε ) , x ( τ ) = ε χ ( τ ; ε ) .
According to PTMS, we can seek the approximate solutions of the variables ξ , φ , and χ in the following forms [34]:
ξ = k = 0 2 ε k ξ k + 1 ( τ 0 , τ 1 , τ 2 ) + O ( ε 3 ) , φ = k = 0 2 ε k φ k + 1 ( τ 0 , τ 1 , τ 2 ) + O ( ε 3 ) , χ = k = 0 2 ε k χ k + 1 ( τ 0 , τ 1 , τ 2 ) + O ( ε 3 ) ,
where τ n = ε n τ ( n = 0 , 1 , 2 ) are various time scales, in which τ 0 and τ 1 , τ 2 are defined as fast and slow time scales, respectively. Therefore, the time derivatives regarding τ can be transformed according to the time scales τ n as the following chain rule:
d d τ = τ 0 + ε τ 1 + ε 2 τ 2 + O ( ε 3 ) , d 2 d τ 2 = 2 τ 0 2 + 2 ε 2 τ 0 τ 1 + ε 2 ( 2 τ 1 2 + 2 2 τ 0 τ 2 ) + O ( ε 3 ) ,
where terms of O ( ε 3 ) and higher are neglected due to their minuteness.
As before, we consider the minuteness of the amplitudes of external forces and moments f j ( j = 1 , 2 , 3 ) , as well as the coefficients of damping c i . Therefore, we can write:
f j = ε 3 f ˜ j , c j = ε 2 c ˜ j ( j = 1 , 2 , 3 ) .
Furthermore, we consider the following:
β = ε β ˜ ,
where the amounts f ˜ j , c ˜ j , a ˜ , b ˜ , and β ˜ are of order unity.
Substituting the new variables into Equations of system (11), the assumed solutions (12), the transformations (13), and the use of the parameters of Equations (14) and (15) into the EOM (8)–(10) and equating the coefficients of like powers of ε in each side obtains the following groups of partial differential equations (PDEs):
Equations of order ε :
2 ξ 1 τ 0 2 + ω 1 2 ξ 1 = 0 ,
2 φ 1 τ 0 2 + φ 1 = 0 ,
2 χ 1 τ 0 2 + ω 2 2 χ 1 + 2 ξ 1 τ 0 2 = 0 .
Equations of order ε 2 :
2 ξ 2 τ 0 2 + ω 1 2 ξ 2 = 2 2 ξ 1 τ 0 τ 1 β ˜ 2 χ 1 τ 0 2 ,
2 ϕ 2 τ 0 2 + ϕ 2 = 2 2 ϕ 1 τ 0 τ 1 2 χ 1 2 ϕ 1 τ 0 2 + ϕ 1 2 ξ 1 τ 0 2 2 ϕ 1 τ 0 χ 1 τ 0 ϕ 1 χ 1 ,
2 χ 2 τ 0 2 + ω 2 2 χ 2 + 2 ξ 2 τ 0 2 = 2 2 χ 1 τ 0 τ 1 3 2 ϕ 1 2 2 2 ξ 1 τ 0 τ 1 .
Equations of order ε 3 :
2 ξ 3 τ 0 2 + ω 1 2 ξ 3 = 2 ξ 1 τ 0 2 2 2 ξ 1 τ 0 τ 2 2 2 ξ 2 τ 0 τ 1 + β ˜ ( ϕ 1 τ 0 ) 2 2 β ˜ 2 χ 1 τ 0 τ 1 β ˜ 2 χ 2 τ 0 2 + β ˜ ϕ 1 2 ϕ 1 τ 0 2 c ˜ 1 ξ 1 τ 0 + f ˜ 1 cos ( p 1 τ 0 ) ,
2 ϕ 3 τ 0 2 + ϕ 3 = 2 ϕ 1 τ 1 2 2 2 ϕ 1 τ 0 τ 2 2 2 ϕ 2 τ 0 τ 1 4 χ 1 2 ϕ 1 τ 0 τ 1 2 χ 1 2 ϕ 2 τ 0 2 2 χ 2 2 ϕ 1 τ 0 2 χ 1 2 2 ϕ 1 τ 0 2 + 2 ϕ 1 2 ξ 1 τ 0 τ 1 + ϕ 2 2 ξ 1 τ 0 2 + ϕ 1 2 ξ 2 τ 0 2 2 ϕ 1 τ 1 χ 1 τ 0 2 ϕ 1 τ 0 χ 1 τ 1 2 ϕ 2 τ 1 χ 1 τ 0 2 ϕ 1 τ 0 χ 2 τ 0 2 χ 1 ϕ 1 τ 0 χ 1 τ 0 χ 1 ϕ 2 χ 2 ϕ 2 + ϕ 1 3 6 c ˜ 2 ϕ 1 τ 0 + f ˜ 2 cos ( p 2 τ 0 ) ,
2 χ 3 τ 0 2 + ω 2 2 χ 3 + 2 ξ 3 τ 0 2 = 2 χ 1 τ 1 2 2 2 χ 1 τ 0 τ 2 2 2 χ 2 τ 0 τ 1 2 ϕ 1 τ 0 ϕ 2 τ 0 χ 1 ( ϕ 1 τ 0 ) 2 2 ξ 1 τ 1 2 2 2 ξ 1 τ 0 τ 2 2 2 ξ 2 τ 0 τ 1 + 3 2 ϕ 1 2 2 ξ 1 τ 0 2 ϕ 1 ϕ 2 c ˜ 3 χ 1 τ 0 + f ˜ 3 cos ( p 3 τ 0 ) .
The above nine PDE (16)–(24) can be solved sequentially one by one. Consequently, the general solutions of Equations (16) and (17) are:
ξ 1 = A 1 e i ω 1 τ 0 + A ¯ 1 e i ω 1 τ 0 ,
ϕ 1 = A 2 e i τ 0 + A ¯ 2 e i τ 0 .
The substitution of Equation (25) into Equation (18) yields a nonhomogeneous equation, and its solution can be written in the form:
χ 1 = A 3 e i ω 2 τ 0 + A ¯ 3 e i ω 2 τ 0 + ω 1 2 ω 1 2 ω 2 2 ( A 1 e i ω 1 τ 0 + A ¯ 1 e i ω 1 τ 0 ) ,
where A j ( j = 1 , 2 , 3 ) represent unknown complex functions of the slow time scales τ 1 , τ 2 , while A ¯ j are the complex conjugates of these functions.
The substitution of the first-order solutions of Equations (25)–(27) into Equations (19)–(21) yields secular terms. The elimination of these terms demands that:
2 i ω 1 A 1 τ 1 β ˜ ω 1 4 ( ω 1 2 ω 2 2 ) A 1 = 0 ,
A 2 τ 1 = 0 .
According to these conditions, the second-order solutions have the following forms:
ξ 2 = β ˜ ω 2 2 ( ω 1 2 ω 2 2 ) A 3 e i ω 2 τ 0 + C C . ,
ϕ 2 = ω 1 A 2 ( ω 1 2 ω 2 2 ) { [ ( 1 + 2 ω 1 ) ( ω 1 2 ω 2 2 ) ] ( ω 1 + 2 ) A 1 e i ( 1 + ω 1 ) τ 0 + [ ( 1 2 ω 1 ) ( ω 1 2 ω 2 2 ) ] ( ω 1 2 ) × A ¯ 1 e i ( 1 ω 1 ) τ 0 } A 2 ω 2 [ ( 1 + 2 ω 2 ) ( ω 2 + 2 ) A 3 e i ( 1 + ω 2 ) τ 0 + ( 1 2 ω 2 ) ( ω 2 2 ) A ¯ 3 e i ( 1 ω 2 ) τ 0 ] + C C .
Here, C C . represents the complex conjugates of the preceding terms.
After grasping the above procedure, we can obtain another condition for eliminating the secular terms by substituting the solutions from Equations (25)–(27), (30) and (31) into Equation (21), in the following form:
2 i ω 2 A 3 τ 1 β ˜ ω 2 4 ω 1 2 ω 2 2 A 3 = 0 .
Therefore, the solution χ 2 has the form:
χ 2 = 3 A 2 A ¯ 2 ω 2 2 3 A 2 2 2 ( ω 2 4 4 ) e 2 i τ 0 + β ˜ ω 1 4 ω 2 2 ( ω 1 2 ω 2 2 ) 3 A 1 e i ω 1 τ 0 + C C .
Making use of the previous solutions to Equations (25)–(27), (30), (31) and (33) and the third-order approximate Equations (22)–(24), secular terms are produced. The eliminating conditions of these terms required that:
2 i ω 1 A 1 τ 2 + i ω 1 c ˜ 1 A 1 β ˜ 2 ω 1 6 A 1 4 ( ω 1 2 ω 2 2 ) 3 [ ω 1 2 ( 4 ω 1 2 + 11 ) + ω 2 2 ( 7 4 ω 2 2 ) + 8 ω 1 2 ω 2 2 ] = 0 ,
2 i A 2 τ 2 + i c ˜ 2 A 2 6 ω 1 4 ( ω 2 2 2 ω 1 2 + 1 ) ( ω 1 2 4 ) ( ω 1 2 ω 2 2 ) 2 A 1 A 2 A ¯ 1 A 2 2 A ¯ 2 2 ω 2 2 ( ω 2 2 4 ) ( ω 2 4 25 ω 2 2 + 24 ) A 2 A 3 A ¯ 3 ( ω 2 2 4 ) ( ω 2 2 ω 1 2 ) { 6 ω 1 2 ( ω 2 2 1 ) + ω 2 2 [ 6 ( 1 ω 2 2 ) + β ˜ ω 2 ( 2 ω 2 2 + 3 ω 2 2 ) ] } = 0 ,
2 i ω 2 A 3 τ 2 + i c ˜ 3 ω 2 A 3 + β ˜ 2 ω 2 6 A 3 2 ( ω 1 2 ω 2 2 ) 2 ( 1 2 ω 1 2 ) 3 + 2 ω 2 ( ω 2 + 2 ) ω 2 ( ω 2 + 2 ) A 2 A 3 A ¯ 2 = 0 .
Based on the previous conditions of Equations (34)–(36), the approximate third-order solutions turn into:
ξ 3 = f ˜ 1 e i p 1 τ 0 2 ( ω 1 2 p 1 2 ) [ 2 β ˜ A 2 2 ( ω 2 2 4 ) + 6 β ˜ A 2 2 ] ( ω 1 2 4 ) ( ω 2 2 4 ) e 2 i τ 0 [ 1 ( ω 1 2 ω 2 2 ) 2 + ω 2 2 ( ω 1 2 ω 2 2 ) 3 ] β ˜ 2 ω 2 4 A 3 e i ω 2 τ 0 + C C ,
ϕ 3 = f ˜ 2 e i p 2 τ 0 2 ( 1 p 2 2 ) + A ¯ 3 A 2 [ d 1 d 2 A ¯ 1 e i ( 1 ω 1 ω 2 ) τ 0 + d 3 d 4 A 1 e i ( 1 + ω 1 ω 2 ) τ 0 ] + A 2 { d 5 d 6 A 1 2 e i ( 1 + 2 ω 1 ) τ 0 + A 3 [ d 7 d 8 A ¯ 1 e i ( 1 ω 1 + ω 2 ) τ 0 + d 9 d 10 A 1 e i ( 1 + ω 1 + ω 2 ) τ 0 ] } + C C ,
χ 3 = f ˜ 1 p 1 2 e i p 1 τ 0 2 ( ω 1 2 p 1 2 ) ( ω 2 2 p 1 2 ) + f ˜ 3 e i p 3 τ 0 2 ( ω 2 2 p 3 2 ) + d 11 d 12 A 1 A 2 A ¯ 2 e i ω 1 τ 0 + A 2 2 [ d 13 d 14 A ¯ 1 e i ( 2 ω 1 ) τ 0 + d 15 d 16 A 1 e i ( 2 + ω 1 ) τ 0 + d 17 d 18 A ¯ 3 e i ( 2 ω 2 ) τ 0 + d 19 d 20 A 3 e i ( 2 + ω 2 ) τ 0 8 β ˜ ( ω 2 2 1 ) ( ω 1 2 4 ) ( ω 2 2 4 ) e 2 i τ 0 ] + C C . ,
where d s ( s = 1 , 2 , 3 , , 20 ) are functions of ω 1 and ω 2 (see Appendix A).
According to the removal conditions of secular terms, Equations (34)–(36), and the following initial conditions, we can estimate the functions A j ( j = 1 , 2 , 3 ) :
y ( 0 ) = B 01 , y ˙ ( 0 ) = B 02 , ϕ ( 0 ) = B 03 , ϕ ˙ ( 0 ) = B 04 , x ( 0 ) = B 05 , x ˙ ( 0 ) = B 06 ,
where B 0 k ( k = 1 , 2 , , 6 ) are known values.
Now, we will investigate the arising cases of resonance according to the previous approximate solutions. These cases can be categorized into primary (main) external resonance, which takes place at p 1 = ω 1 , p 2 = 1 , p 3 = ω 2 , and internal resonance, which occurs at ω 1 = ω 2 , ω 1 = ω 2 , ω 1 = 2 , ω 1 = 2 , ω 2 = 2 , ω 2 = 2 .
In this situation, we investigate the three primary external resonances, p 1 ω 1 , p 2 1 , and p 3 ω 2 , simultaneously. These resonances depict the closeness of p 1 , p 2 and p 3 to ω 1 , 1 and ω 2 , respectively. Therefore, detuning parameters σ j ( j = 1 , 2 , 3 ) can be introduced according to:
p 1 = ω 1 + σ 1 , p 2 = 1 + σ 2 , p 3 = ω 2 + σ 3 .
Since the detuning parameters are known as the distances from the oscillations to the strict resonance [35,36], we can use the small parameter ε to express these parameters as follows:
σ j = ε σ ˜ j ( j = 1 , 2 , 3 ) .
Substituting parameters (41) and (42) into Equations (8)–(10) and removing the produced secular terms obtains the following conditions of solvability:
2 i ω 1 A 1 τ 2 + i ω 1 c ˜ 1 A 1 β ˜ 2 ω 1 6 A 1 4 ( ω 1 2 ω 2 2 ) 3 [ ω 1 2 ( 4 ω 1 2 + 11 ) + ω 2 2 ( 7 4 ω 2 2 ) + 8 ω 1 2 ω 2 2 ] f ˜ 1 2 e i σ ˜ 1 τ 1 = 0 , 2 i A 2 τ 2 + i c ˜ 2 A 2 6 ω 1 4 ( ω 2 2 2 ω 1 2 + 1 ) ( ω 1 2 4 ) ( ω 1 2 ω 2 2 ) 2 A 1 A 2 A ¯ 1 A 2 A 3 A ¯ 3 ( ω 2 2 4 ) ( ω 2 2 ω 1 2 ) × { 6 ω 1 2 ( ω 2 2 1 ) + ω 2 2 [ 6 ( 1 ω 2 2 ) + β ˜ ω 2 ( 2 ω 2 2 + 3 ω 2 2 ) ] } 2 2 A ¯ 2 2 ω 2 2 ( ω 2 2 4 ) × ( ω 2 4 25 ω 2 2 + 24 ) f ˜ 2 2 e i σ ˜ 2 τ 1 = 0 , 2 i ω 2 3 τ 2 + i c ˜ 3 ω 2 A 3 ( 2 ω 1 2 1 ) 2 ( ω 1 2 ω 2 2 ) 2 A 3 β ˜ 2 ω 2 6 [ 3 + 2 ω 2 ( ω 2 + 2 ) ] ω 2 ( ω 2 + 2 ) × A 2 A 3 A ¯ 2 f ˜ 3 2 e i σ ˜ 3 τ 1 = 0 .
It is obvious that these conditions constitute a system of three non-linear PDEs in terms of the functions A j ( j = 1 , 2 , 3 ) , in which A j = A j ( τ 1 , τ 2 ) . Therefore, we can write them in polar forms as follows:
A j = h ˜ j ( τ 1 , τ 2 ) 2 e i ψ j ( τ 1 , τ 2 ) , h j = ε h ˜ j ( j = 1 , 2 , 3 ) .
Here, h ˜ j and ψ j are real functions of the amplitudes and phases, respectively, for the functions ξ , φ and χ , while h j are the amplitudes of y , ϕ , and x , as presented in the assumptions of Equations (11) and (12).
Making use of the following modified phases:
θ j ( τ 1 , τ 2 ) = τ j σ ˜ j ψ j ( τ 1 , τ 2 )
and the time derivative operators (13), we can transform the PDE system (43) into another ODE system. The separation of the real and imaginary parts of the resulting system yields the following six ODEs from the first order that describe the modulation equations in terms of modified phases θ j and amplitudes h j in the following forms:
h 1 d θ 1 d τ = h 1 σ 1 + β 2 ω 1 5 h 1 4 ( ω 1 2 ω 2 2 ) 3 [ ω 1 2 ( 4 ω 1 2 + 11 ) + ω 2 2 ( 7 4 ω 2 2 ) + 8 ω 1 2 ω 2 2 ] + f 1 2 ω 1 cos θ 1 , d h 1 d τ = 1 2 h 1 c 1 + f 1 2 ω 1 sin θ 1 , h 2 d θ 2 d τ = h 2 σ 2 6 ω 1 2 h 1 2 h 2 8 ( ω 1 2 4 ) ( ω 1 2 ω 2 2 ) 2 ( ω 2 2 2 ω 1 2 + 1 ) + h 2 3 16 ω 2 2 ( ω 2 2 4 ) × ( ω 2 4 25 ω 2 2 + 24 ) 3 h 3 2 h 2 4 ( ω 2 2 4 ) ( ω 2 2 1 ) + f 2 2 cos θ 2 , d h 2 d τ = 1 2 h 2 c 2 + f 2 2 sin θ 2 , h 3 d θ 3 d τ = h 3 σ 3 + β 2 ω 2 5 h 3 4 ( ω 1 2 ω 2 2 ) 3 ( ω 2 2 3 ω 1 2 ) + h 2 2 h 3 8 ω 2 2 ( ω 2 + 2 ) [ 3 + 2 ω 2 ( ω 2 + 2 ) ] + f 3 2 ω 2 cos θ 3 , d h 3 d τ = 1 2 h 3 c 3 + f 3 2 ω 2 sin θ 3 .
According to the conditions (40), we can solve the above modulation in Equation (46). The time histories of the amplitudes h j ( j = 1 , 2 , 3 ) and the modified phases θ j are portrayed in Figure 2, Figure 3 and Figure 4, according to the uses of the following data:
σ 1 = 0.01 , σ 2 = 0.02 , σ 3 = 0.03 , f 1 = 0.00002 , f 2 = 0.002 , f 3 = 0.00005 , β = 0.167 , c 1 = 0.00188 , c 2 = 0.0135 , c 3 = 0.0034 , p 1 = ω 1 + σ 1 , p 2 = 1 + σ 2 , p 3 = ω 2 + σ 3 .
These figures were calculated for different values of ω ¯ 1 = ( 3.02 , 3.43 , 4.08 ) ,   ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) , and ω ¯ 3 = ( 7.07 , 7.75 , 8.37 ) as in the parts [(a), (b)], [(c), (d)], and [(e), (f)], respectively. The oscillations of the waves describing h 1 and h 3 decrease with increasing time, behaving in a decay manner as seen in Figure 2 and Figure 4, while the variation of the amplitude h 2 behaves in the form of steady-state motion as shown in Figure 3. On the other hand, the variation of the modified phases θ j increases with time to a specified value, then it tends to exhibit a stationary behavior, as portrayed in Figure 2, Figure 3 and Figure 4. The conclusion that can be made here is that the modified amplitudes and phases behave in a stable manner.
Figure 5, Figure 6 and Figure 7 illustrate the time histories of the attained approximate analytic solutions of Y , X , and ϕ till the third order of approximation. These figures were drawn when the previously used data with various values of ω ¯ j ( j = 1 , 2 , 3 ) were used. According to the plotted curves in these figures, we can conclude that these curves have periodic or quasi-periodic forms of wave packets, in which their amplitudes increase (see Figure 5 and Figure 7) and decrease, as shown in Figure 6, with the passage of time until the waves behave in the form of stationary waves at the end of the examined period of time.
To follow up on the vibrations of the examined system to investigate whether it is stable and to understand the specific transference of motion, the system of Equation (46) permits us to focus on this knowledge. The plotted trajectories of the modified amplitudes and phases are regarded as the best way to depict the system’s behavior.
The track projections on specific space planes are plotted in Figure 8, Figure 9 and Figure 10, when ω ¯ 1 = 4.08 , ω ¯ 2 = 5.72 , and ω ¯ 3 = 8.37 , respectively. The curves of these figures have spiral forms toward one point, in which the green, black, and red solid points in these figures refer to a stable state. After the transient state, we find that all pathways tend to have steady-state behavior.

4. Solutions in the Steady State

The main objective of this section is to examine the oscillations of the considered system in the steady state. It is known that this case is generated when the transient processes disappear due to the presence of damping [37,38]. Then, we consider the zero value of the left hand side of Equation (46), i.e., d θ j d t = 0 , d h j d t = 0 ( j = 1 , 2 , 3 ) . Therefore, we can easily obtain the following system of six algebraic equations:
h 1 σ 1 + β 2 ω 1 5 h 1 4 ( ω 1 2 ω 2 2 ) 3 [ ω 1 2 ( 4 ω 1 2 + 11 ) + ω 2 2 ( 7 4 ω 2 2 ) + 8 ω 1 2 ω 2 2 ] + f 1 2 ω 1 cos θ 1 = 0 , 1 2 h 1 c 1 + f 1 2 ω 1 sin θ 1 = 0 , h 2 σ 2 6 ω 1 2 h 1 2 h 2 8 ( ω 1 2 4 ) ( ω 1 2 ω 2 2 ) 2 ( ω 2 2 2 ω 1 2 + 1 ) + h 2 3 16 ω 2 2 ( ω 2 2 4 ) ( ω 2 4 25 ω 2 2 + 24 ) 6 h 3 2 h 2 8 ( ω 2 2 4 ) ( ω 2 2 1 ) + f 2 2 cos θ 2 = 0 , 1 2 h 2 c 2 + f 2 2 sin θ 2 = 0 , h 3 σ 3 + β 2 ω 2 5 h 3 4 ( ω 1 2 ω 2 2 ) 3 ( ω 2 2 3 ω 1 2 ) + h 2 2 h 3 8 ω 2 2 ( ω 2 + 2 ) [ 3 + 2 ω 2 ( ω 2 + 2 ) ] + f 3 2 ω 2 cos θ 3 = 0 , 1 2 h 3 c 3 + f 3 2 ω 2 sin θ 3 = 0 .
A closer look into this system shows that if we remove the phases θ j ( j = 1 , 2 , 3 ) , the following non-linear algebraic equations of amplitudes and frequencies are obtained:
f 1 2 = 4 ω 1 2 { [ h 1 σ 1 + β 2 ω 1 5 h 1 4 ( ω 1 2 ω 2 2 ) 3 ( ω 1 2 ( 4 ω 1 2 + 11 ) + ω 2 2 ( 7 4 ω 2 2 ) + 8 ω 1 2 ω 2 2 ) ] 2 + 1 4 h 1 2 c 1 2 } , f 2 2 = 4 { [ h 2 σ 2 6 ω 1 2 h 1 2 h 2 8 ( ω 1 2 4 ) ( ω 1 2 ω 2 2 ) 2 ( ω 2 2 2 ω 1 2 + 1 ) + h 2 3 16 ω 2 2 ( ω 2 2 4 ) ( ω 2 4 25 ω 2 2 + 24 ) 6 h 3 2 h 2 8 ( ω 2 2 4 ) ( ω 2 2 1 ) ] 2 + 1 4 h 2 2 c 2 2 } , f 3 2 = 4 ω 2 2 { [ h 3 σ 3 + β 2 ω 2 5 h 3 4 ( ω 1 2 ω 2 2 ) 3 ( ω 2 2 3 ω 1 2 ) + h 2 2 h 3 8 ω 2 2 ( ω 2 + 2 ) × ( 3 + 2 ω 2 ( ω 2 + 2 ) ) ] 2 + 1 4 h 3 2 c 3 2 } .
For further investigation of the stability, we will examine the behavior of the dynamical system in a region close to the fixed points. To complete this task, we will consider the following forms of amplitudes and phases:
h 1 = h 10 + h 11 , h 2 = h 20 + h 21 , h 3 = h 30 + h 31 , θ 1 = θ 10 + θ 11 , θ 2 = θ 20 + θ 21 , θ 3 = θ 30 + θ 31 ,
where h 10 ,   h 20 ,   h 30 ,   θ 10 ,   θ 20 , and θ 30 are the solutions in the steady state, while h 11 ,   h 21 ,   h 31 ,   θ 11 ,   θ 21 , and θ 31 are the corresponding perturbations, which are assumed to be very small.
The substitution of Equation (49) into Equation (46) yields the following equations:
h 10 d θ 11 d τ = h 11 σ 1 f 1 2 ω 1 θ 11 sin θ 10 + G 1 , d h 11 d τ = 1 2 h 11 c 1 + f 1 2 ω 1 θ 11 cos θ 10 , h 20 d θ 21 d τ = h 21 σ 2 f 2 2 θ 21 sin θ 20 + G 2 , d h 21 d τ = 1 2 h 21 c 2 + f 2 2 θ 21 cos θ 20 , h 30 d θ 31 d τ = h 31 σ 3 f 3 2 ω 2 θ 31 sin θ 30 + G 3 , d h 31 d τ = 1 2 h 31 c 3 + f 3 2 ω 2 θ 31 cos θ 30 ,
where:
G 1 = β 2 ω 1 5 h 11 4 ( ω 1 2 ω 2 2 ) 3 [ ω 1 2 ( 4 ω 1 2 + 11 ) + ω 2 2 ( 7 4 ω 2 2 ) + 8 ω 1 2 ω 2 2 ] , G 2 = 6 ( ω 2 2 1 ) 8 ( ω 2 2 4 ) ( 2 h 20 h 30 h 31 + h 30 2 h 21 ) 3 h 20 2 h 21 16 ω 2 2 ( ω 2 2 4 ) ( ω 2 4 25 ω 2 2 + 24 ) 6 ω 1 2 ( ω 2 2 2 ω 1 2 + 1 ) 8 ( ω 1 2 4 ) ( ω 1 2 ω 2 2 ) 2 ( 2 h 20 h 10 h 11 + h 10 2 h 21 ) , G 3 = [ 3 + 2 ω 2 ( ω 2 + 2 ) ] 8 ω 2 2 ( ω 2 + 2 ) ( 2 h 30 h 20 h 21 + h 20 2 h 31 ) + β 2 ω 2 5 h 31 4 ( ω 1 2 ω 2 2 ) 3 ( ω 2 2 3 ω 1 2 ) .
Based on the minuteness of the perturbation functions h 11 ,   h 21 ,   h 31 ,   θ 11 ,   θ 21 , and θ 31 , we can express each solution in the form of a linear combination of S k e λ τ   ( k = 1 , 2 , , 6 ) , where S k and λ are constants and the eigenvalue of the unknown perturbation, respectively.
If the solutions in the steady state (fixed points) h 10 ,   h 20 ,   h 30 ,   θ 10 ,   θ 20 , and θ 31 are stable asymptotically, then the real parts of the roots of the following characteristic equation must be negative:
λ 6 + Γ 1 λ 5 + Γ 2 λ 4 + Γ 3 λ 3 + Γ 4 λ 2 + Γ 5 λ + Γ 6 = 0 ,
where Γ k are functions depending on the parameters h 10 ,   h 20 ,   h 30 ,   θ 10 ,   θ 20 ,   θ 30 ,   f 1 ,   f 2 , and f 3 (see Appendix B).
Referring to the above, the substantial conditions for the stability of the solutions in the steady state can be written according to the criterion of Routh–Hurwitz [32] in the following forms:
Γ 1 > 0 ,   Γ 1 Γ 2 Γ 3 > 0 , Γ 3 ( Γ 1 Γ 2 Γ 3 ) + Γ 1 ( Γ 5 Γ 1 Γ 4 ) > 0 , Γ 1 ( Γ 2 Γ 3 Γ 4 Γ 2 2 Γ 5 + 2 Γ 4 Γ 5 Γ 3 Γ 6 ) + Γ 3 ( Γ 2 Γ 5 Γ 3 Γ 4 ) + Γ 1 2 ( Γ 2 Γ 6 Γ 4 2 ) Γ 5 2 > 0 , Γ 5 ( Γ 4 ( Γ 1 Γ 2 Γ 3 + Γ 3 2 + Γ 1 2 Γ 4 ) + Γ 1 ( Γ 2 2 2 Γ 4 ) Γ 5 Γ 2 Γ 3 Γ 5 + Γ 5 2 ) + ( Γ 3 3 Γ 1 Γ 3 ( Γ 2 Γ 3 + 3 Γ 5 ) + Γ 1 2 ( Γ 3 Γ 4 + 2 Γ 2 Γ 5 ) ) Γ 6 Γ 1 3 Γ 6 2 > 0 , Γ 6 ( Γ 5 ( Γ 4 ( Γ 1 Γ 2 Γ 3 + Γ 3 2 + Γ 1 2 Γ 4 ) + Γ 1 ( Γ 2 2 2 Γ 4 ) Γ 5 Γ 2 Γ 3 Γ 5 + Γ 5 2 ) + ( Γ 3 3 Γ 1 Γ 3 ( Γ 2 Γ 3 + 3 Γ 5 ) + Γ 1 2 ( Γ 3 Γ 4 + 2 Γ 2 Γ 5 ) ) Γ 6 Γ 1 3 Γ 6 2 ) > 0 .

5. Non-Linear Analysis of Stability

The main aim of this section is to investigate the stability analysis of the considered system using the non-linear stability approach of Routh–Hurwitz. The examined system is exposed to various forces and moments, such as the excited harmonic forces F 1 , F 3 and the moment M ϕ . It is found that the natural frequencies ω j ¯ ( j = 1 , 2 , 3 ) , damping coefficients c j , and the parameters of detuning perform a crucial destabilizing role in the criteria’s stability. To plot the stability diagrams of the system of Equation (46), the previous values of the different parameters are used. The modified amplitudes h j ( j = 1 , 2 , 3 ) are shown with the detuning parameter σ 1 in Figure 11, Figure 12, Figure 13 and Figure 14 for different values of ω j ¯ when σ 2 = σ 3 = 0 . The outlined curves of these figures contain the possible fixed points. Resonance curves are plotted in Figure 11, Figure 12, Figure 13 and Figure 14 according to the previous data of parameters.
Each of the resonance curves has different stable and unstable regions with the various values of the natural frequencies ω ¯ 1 and ω ¯ 2 (see Figure 11, Figure 12, Figure 13 and Figure 14). More precisely, the inspection of parts of Figure 11 shows that the instability and stability ranges are 0.1 < σ 1 < 0.009 and 0.009 σ 1 < 0.1 , respectively, as indicated in Figure 11 in blue. On the other hand, when ω ¯ 1 = 3.43 , the red curves represent the instability and stability ranges of 0.1 < σ 1 < 0.012 and 0.012 σ 1 < 0.1 , respectively, while the black curves represent the areas of instability and stability when ω ¯ 1 = 3.02 , which are 0.1 < σ 1 < 0.014 and 0.014 σ 1 < 0.1 , respectively. The included curves of each part have different fixed points associated with the different values of ω ¯ 1 , as indicated in Figure 11. Moreover, the instability regions decrease with an increase in the value of ω ¯ 1 , due to the impact of the natural frequency ω ¯ 1 .
Figure 12 illustrates the different ranges of stability criterion, where the instability and stability areas, indicated in blue, are determined according to the ranges 0.1 < σ 1 < 0.009 and 0.009 σ 1 < 0.1 , respectively. At ω ¯ 2 = 5.72 , the red curves show the instability and stability areas with ranges 0.1 < σ 1 < 0.021 and 0.021 σ 1 < 0.1 , respectively, while the black curves show the instability and stability areas at ω ¯ 2 = 6.68 , which are 0.1 < σ 1 < 0.033 and 0.033 σ 1 < 0.1 , respectively.
According to the calculations of the graphed curves in Figure 13 and Figure 14, we can state that three different peak fixed points are produced when ω ¯ 2 has different values for the resonance curves of h 2 and h 3 , which are varied with σ 1 . The instability and the stability areas become 0.1 < σ 1 < 0.009 and 0.009 σ 1 < 0.1 , respectively, denoted in blue (see Figure 13 and Figure 14), while at ω ¯ 2 = 5.72 , the red curves indicate the instability and stability areas 0.1 < σ 1 < 0.021 and 0.021 σ 1 < 0.1 , respectively. On the other hand, the black curves represent the areas of instability and stability when ω ¯ 2 = 6.68 , which are 0.1 < σ 1 < 0.033 and 0.033 σ 1 < 0.1 , respectively.
To clarify the properties of the non-linear amplitudes of the system of Equations (46) and to examine their stabilities, we consider the following transformations [39]:
A j = [ U j ( τ 1 , τ 2 ) + i V j ( τ 1 , τ 2 ) ] e i σ ˜ j τ 1 ; ( j = 1 , 2 , 3 )
where:
U j = ε u j , V j = ε v j .
Substituting the operators (13) and the previous transformation (53) into the system of the ODE (46), we can equate the real parts and the imaginary parts with zero to obtain of the following system:
d u 1 d τ = v 1 { σ 1 + G 4 [ ω 1 2 ( 4 ω 1 2 + 11 ) + ω 2 2 ( 7 4 ω 2 2 ) + 8 ω 1 2 ω 2 2 ] } c 1 2 u 1 , d v 1 d τ = u 1 { σ 1 + G 4 [ ω 1 2 ( 4 ω 1 2 + 11 ) + ω 2 2 ( 7 4 ω 2 2 ) + 8 ω 1 2 ω 2 2 ] } + c 1 2 v 1 f 1 4 ω 1 , d u 2 d τ = v 2 [ σ 2 G 5 ( u 1 2 + v 1 2 ) + G 6 ( u 2 2 + v 2 2 ) + G 7 ( u 3 2 + v 3 2 ) ] c 2 2 u 2 , d v 2 d τ = u 2 [ σ 2 + G 5 ( u 1 2 + v 1 2 ) G 6 ( u 2 2 + v 2 2 ) G 7 ( u 3 2 + v 3 2 ) ] c 2 2 v 2 f 2 4 , d u 3 d τ = v 3 [ σ 3 + G 8 + G 9 ( u 2 2 + v 2 2 ) ] c 3 2 u 3 , d v 3 d τ = u 3 [ σ 3 G 8 G 9 ( u 2 2 + v 2 2 ) ] c 3 2 v 3 f 3 4 ω 2 ,
where:
G 4 = β 2 ω 1 5 4 ( ω 1 2 ω 2 2 ) 3 , G 5 = 3 ω 1 4 ( ω 2 2 2 ω 1 2 + 1 ) 2 ( ω 1 2 4 ) ( ω 1 2 ω 2 2 ) 2 , G 6 = ( ω 2 4 25 ω 2 2 + 24 ) 4 ω 2 2 ( ω 2 2 4 ) , G 7 = 3 ( ω 2 2 1 ) 8 ( ω 2 2 4 ) , G 8 = β 2 ω 2 5 4 ( ω 1 2 ω 2 2 ) 3 , G 9 = [ 3 + 2 ω 2 ( ω 2 + 2 ) ] 8 ω 2 2 ( ω 2 + 2 ) .
The modified amplitudes can be subsequently rationalized through the time in various parametric regions; the properties of the amplitudes are presented in phase plane curves shown in Figure 15, Figure 16, Figure 17, Figure 18, Figure 19 and Figure 20, taking into account the data from previously used parameters.
Parts of Figure 15, Figure 16, Figure 17, Figure 18, Figure 19 and Figure 20 explore the changing of the new justified phases u j and v j ( j = 1 , 2 , 3 ) , as stated in the system of Equation (54), via τ as shown in parts (a) and (b). The projections of the modulation equation trajectories on the phase planes u j v j are graphed in part (c) of these figures. The indicated curves of Figure 15, Figure 16, Figure 17, Figure 18, Figure 19 and Figure 20 are calculated when ω ¯ 1 = ( 3.02 , 3.43 , 4.08 ) and ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) , respectively, besides the values of the detuning parameters σ 1 = 0.01 , σ 2 = 0.02 , and σ 3 = 0.03 .
The inspection of parts of Figure 15 and Figure 18 shows the time histories of u 1 , v 1 and the projections of the modulation equation trajectory on the phase plane u 1 v 1 . It is clear that u 1 , v 1 and u 1 v 1 increase gradually when time increases, and they behave in the form of a straight line due to the system of Equation (54). Conversely, the decay behavior of the waves describing the variation of u 2 , v 2 , u 3 , and v 3 is shown in Figure 16a,b, Figure 17a,b, Figure 19a,b, and Figure 20a,b, for different values of ω ¯ 1 and ω ¯ 2 . Spiral curves are plotted in the planes u 2 v 2 and u 3 v 3 in Figure 16c. Figure 16c in line with the equations of the system of Equation (54).

6. Conclusions

This work focused on the dynamical motion of a novel three DOF auto-parametric system consisting of the main motion of a non-linear Duffing oscillator and the secondary motion of a spring pendulum. The governing EOM of the system were obtained using Lagrange’s equations and solved using PTMS up to the third approximation. The solvability conditions and the equations of the modulations were obtained in view of the examined primary external resonance case. The curves of time histories and resonance were drawn, taking into consideration the selected values of some data. The non-linear analysis of the solutions in the steady state was presented and interpreted. The stability and instability areas were plotted to clarify their corresponding regions. The properties of the non-linear amplitudes of the modulation system and their stabilities were clarified and examined. The achieved results of this work can be applied to deal with the vibration of swaying buildings and to reduce the vibration of rotor dynamics, as well as in the fields of mechanics and space engineering. These outcomes generalized those of previous works in [5,31] for the absence of an applied force and moment in the unstarched arm of a pendulum.

Author Contributions

Conceptualization, T.S.A. and R.S.; formal analysis, T.S.A., A.A. and A.S.E.; funding acquisition, R.S.; investigation, T.S.A., A.A. and A.S.E.; methodology, T.S.A. and R.S.; project administration, R.S.; resources, T.S.A.; software, T.S.A. and A.S.E.; supervision, T.S.A. and A.A.; validation, T.S.A. and R.S.; Writing—original draft, A.S.E.; Writing—review and editing, T.S.A., R.S. and A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Ministry of Science and Higher Education in Poland, 0612/SBAD/3576.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing does not apply to this article as no datasets were generated or analyzed during the current study.

Conflicts of Interest

The authors declare that they have no conflict of interest.

Appendix A

d 1 = ω 1 { 2 ω 1 4 ( ω 2 2 ) ω 2 + 2 ω 1 3 ( ω 2 1 ) 2 ( ω 2 + 1 ) + ω 2 ( ω 2 2 ) ( 2 ω 2 1 ) × ( ω 2 2 + 1 ) + ω 1 2 [ 5 + ω 2 ( 24 + ω 2 ( 21 2 ( ω 2 4 ) ω 2 ) ) ] 2 ω 1 ( ω 2 1 ) × [ 1 + ω 2 ( 7 + ω 2 ( 6 + ( ω 2 3 ) ω 2 ) ) ] } , d 2 = ( ω 1 2 ) 2 [ 1 ( 1 ω 1 ω 2 ) 2 ] ω 2 ( ω 1 2 ω 2 2 ) , d 3 = ω 1 { 2 ω 1 4 ( ω 2 2 ) ω 2 2 ω 1 3 ( ω 2 1 ) 2 ( ω 2 + 1 ) + ω 2 ( ω 2 2 ) ( 2 ω 2 1 ) × ( ω 2 2 + 1 ) + 2 ω 1 ( ω 2 1 ) [ 1 + ω 2 ( 7 + ω 2 ( 6 + ( ω 2 3 ) ω 2 ) ) ] } , d 4 = ( ω 1 + 2 ) ( ω 2 2 ) ω 2 ( ω 1 2 ω 2 2 ) [ 1 ( 1 + ω 1 ω 2 ) 2 ] 2 , d 5 = ω 1 3 [ 1 6 ω 1 10 ω 1 2 + 4 ω 1 4 ( 1 + 6 ω 1 4 ω 1 2 ) ω 2 2 ] , d 6 = ( ω 1 + 2 ) ( 1 + 2 ω 1 ) 2 ( ω 1 2 ω 2 2 ) 2 [ 1 ( 1 + 2 ω 1 ) 2 ] , d 7 = ω 1 { 2 ω 1 5 ( ω 2 1 ) ( ω 2 + 1 ) 2 + 2 ω 1 6 ω 2 ( ω 1 + 2 ) + ω 2 3 ( ω 2 2 + 1 ) [ 1 + ω 2 ( 2 + β ˜ ω 2 ) ] ω 1 4 { 5 + ω 2 [ 24 + ω 2 ( 21 + 4 ω 2 ( ω 2 + 3 ) ) ] } + 2 ω 1 3 ( ω 2 + 1 ) { 1 + ω 2 [ 7 + ω 2 ( 5 + ω 2 × ( 3 + 2 ω 2 ) ) ] } + ω 1 2 ω 2 { 2 + ω 2 2 [ 20 + ω 2 [ 16 + ω 2 ( β ˜ ( ω 2 + 2 ) + 2 ( 3 + ω 2 ) ) ] ] } 2 ω 2 2 ω 1 [ 1 + ω 2 [ 8 + ω 2 ( 13 + ω 2 ( 9 + β ˜ ( ω 2 + 2 ) + ω 2 ( ω 2 + 4 ) ) ) ] ] } ,
d 8 = ( ω 1 2 ) ( ω 2 + 2 ) ω 2 ( ω 1 2 ω 2 2 ) 2 [ 1 ( 1 ω 1 + ω 2 ) 2 ) ] , d 9 = ω 1 { 2 ω 1 5 ( ω 2 1 ) ( ω 2 + 1 ) 2 + 2 ω 1 6 ω 2 ( ω 1 + 2 ) + ω 2 3 ( ω 2 2 + 1 ) ( 1 + ω 2 ( 2 + β ˜ ω 2 ) ) ω 1 4 [ 5 + ω 2 ( 24 + ω 2 ( 21 + 4 ω 2 ( ω 2 + 3 ) ) ) ] 2 ω 1 3 ( ω 2 + 1 ) [ 1 + ω 2 ( 7 + ω 2 × ( 5 + ω 2 ( 3 + 2 ω 2 ) ) ) ] + ω 1 2 ω 2 [ 2 + ω 2 2 [ 20 + ω 2 ( 16 + ω 2 ( β ˜ ( ω 2 + 2 ) + 2 ( 3 + ω 2 ) ) ) ] ] + 2 ω 2 2 ω 1 [ 1 + ω 2 ( 8 + ω 2 ( 13 + ω 2 ( 9 + β ˜ ( ω 2 + 2 ) + ω 2 ( ω 2 + 4 ) ) ) ) ] } , d 10 = ( ω 2 + 2 ) ω 2 ( ω 2 2 ω 1 2 ) 2 [ 1 ( 1 + ω 1 + ω 2 ) 2 ] , d 11 = ω 1 [ 3 4 ω 1 + ω 1 2 + 4 ω 1 3 + ω 1 4 ( 1 + ω 1 ) ( 3 + ω 1 ) ω 2 2 ] , d 12 = ( ω 1 + 2 ) ( ω 2 2 ω 1 2 ) 2 , d 13 = ω 1 { 1 + 1 2 ω 1 [ 4 ω 1 ( 4 + ( 2 + ω 1 ) ω 1 ) ] + 1 2 [ 2 + ( 2 + ω 1 ) ω 1 ] ω 2 2 } , d 14 = ( ω 1 2 ) ( ω 1 2 ω 2 2 ) [ ω 2 2 ( 2 ω 1 ) 2 ] , d 15 = ω 1 { 1 + ω 1 2 [ 4 + ω 1 ( 4 + ( 2 ω 1 ) ω 1 ) ] + 1 2 [ 2 + ( ω 1 2 ) ω 1 ] ω 2 2 } , d 16 = ( ω 1 + 2 ) ( ω 1 2 ω 2 2 ) [ ω 2 2 ( 2 + ω 1 ) 2 ] , d 17 = [ 1 + ( 2 3 ω 2 ) ω 2 ] , d 18 = ω 2 ( ω 2 2 ) [ ω 2 2 ( 2 ω 2 ) 2 ] , d 19 = 1 + ( 2 + 3 ω 2 ) ω 2 , d 20 = ω 2 ( ω 2 + 2 ) [ ω 2 2 ( 2 + ω 2 ) 2 ] .

Appendix B

Γ 1 = 1 2 [ c 1 + c 2 + c 3 + f 1 sin θ 10 h 10 ω 1 + f 2 sin θ 20 h 20 + f 3 sin θ 30 h 30 ω 2 ] , Γ 2 = 1 h 10 h 20 h 30 ω 1 ω 2 { f 1 [ 1 4 sin θ 10 sin θ 30 h 20 f 3 + h 30 ( 1 4 sin θ 10 sin θ 20 f 2 + h 20 ( 1 4 sin θ 10 c 1 + 1 4 sin θ 10 c 2 + 1 4 sin θ 10 c 3 1 2 cos θ 10 ν 1 1 2 cos θ 10 σ 1 ) ) ω 2 ] + h 10 ω 1 [ f 2 ( 1 4 sin θ 20 sin θ 30 f 3 + h 30 ( 1 4 sin θ 20 c 1 + 1 4 sin θ 20 c 2 + 1 4 sin θ 20 c 3 + 1 2 cos θ 20 ν 2 h 30 2 1 2 cos θ 20 ν 3 h 20 2 1 2 cos θ 20 ν 4 h 10 2 1 2 cos θ 20 σ 2 ) ω 2 ) + h 20 ( f 3 ( 1 4 sin θ 30 c 3 1 2 cos θ 30 ( ν 5 h 20 2 + v 6 + σ 3 ) ) + c 1 ( 1 4 sin θ 30 f 3 + 1 4 ( c 2 + c 3 ) h 30 ω 2 + 1 4 c 2 ( sin θ 30 f 3 + c 3 h 30 ω 2 ) ) ] } ,
Γ 3 = 1 8 h 10 h 20 h 30 ω 1 ω 2 { f 1 ( f 2 ( sin θ 10 sin θ 20 sin θ 30 f 3 + h 30 ( sin θ 10 sin θ 20 c 1 + h 30 ( sin θ 10 sin θ 20 c 2 + h 30 ( sin θ 10 sin θ 20 c 3 2 cos θ 10 sin θ 20 ν 1 + 2 cos θ 10 sin θ 20 h 30 2 ν 2 2 cos θ 20 sin θ 10 h 20 2 ν 3 2 cos θ 20 sin θ 10 h 10 2 ν 4 2 cos θ 10 sin θ 20 σ 1 2 cos θ 20 sin θ 10 σ 2 ) ω 2 ) + h 20 ( sin θ 10 sin θ 30 c 3 f 3 2 cos θ 10 sin θ 30 f 3 ν 1 2 cos θ 30 sin θ 10 f 3 h 20 2 ν 5 2 cos θ 30 sin θ 10 f 3 ν 6 2 cos θ 10 sin θ 30 f 3 σ 1 2 cos θ 30 sin θ 10 f 3 σ 3 2 cos θ 10 c 3 h 30 σ 1 ω 2 2 cos θ 10 c 3 h 30 ν 1 ω 2 + c 3 sin θ 10 ( sin θ 30 f 3 + ( c 2 + c 3 ) h 30 ω 2 ) + c 2 ( sin θ 10 sin θ 30 f 3 + h 30 ( c 3 sin θ 10 2 cos θ 10 ( ν 1 + σ 1 ) ) ω 2 ) ) ) + h 10 ω 1 ( c 2 ( f 3 h 20 ( sin θ 30 c 3 - 2 cos θ 30 ( h 20 2 ν 5 + ν 6 + σ 3 ) ) + sin θ 20 f 2 ( sin θ 30 f 3 + c 3 h 30 ω 2 ) ) + f 2 ( f 3 ( 2 cos θ 20 sin θ 30 h 30 2 ν 2 2 cos θ 20 sin θ 30 h 10 2 ν 4 + h 20 2 ( 2 cos θ 20 sin θ 30 ν 3 2 cos θ 30 × sin θ 20 ν 5 ) 2 cos θ 30 sin θ 20 ν 6 2 cos θ 20 sin θ 30 σ 2 2 cos θ 30 sin θ 20 σ 3 ) + c 3 [ sin θ 20 × sin θ 30 f 3 + 2 cos θ 20 h 30 ( h 30 2 ν 2 h 20 2 ν 3 h 10 2 ν 4 σ 2 ) ω 2 ] ) + c 1 ( f 2 ( sin θ 20 sin θ 30 f 3 + h 30 × ( sin θ 20 c 2 + sin θ 20 c 3 + 2 cos θ 20 ( h 30 2 ν 2 h 20 2 ν 3 h 10 2 ν 4 σ 2 ) ω 2 ) ) + h 20 ( f 3 ( sin θ 30 c 3 + 2 cos θ 30 ( h 20 2 ν 5 ν 6 σ 3 ) ) + 2 c 2 ( sin θ 30 f 3 + c 3 h 30 ω 2 ) ) ) } ,
Γ 4 = 1 h 10 h 20 h 30 ω 1 ω 2 { 1 16 sin θ 10 sin θ 20 sin θ 30 c 3 f 1 f 2 f 3 1 8 [ ν 1 f 1 f 3 cos θ 10 sin θ 30 × ( f 2 sin θ 20 c 3 h 20 ) + f 1 f 2 f 3 cos θ 20 sin θ 10 sin θ 30 ( h 30 2 ν 2 h 20 2 ν 3 ) f 1 f 2 f 3 sin θ 10 × ( cos θ 20 sin θ 30 h 10 2 ν 4 + cos θ 30 sin θ 20 h 20 2 ν 5 ) + f 1 f 3 cos θ 30 ( cos θ 10 h 20 3 ν 1 ν 5 sin θ 10 × sin θ 20 f 2 ν 6 ) + f 1 f 3 cos θ 10 ( h 20 cos θ 30 ( ν 1 ν 6 + ν 5 h 20 2 σ 1 + ν 6 σ 1 ) σ 1 sin θ 30 ( f 2 sin θ 20 + c 3 h 20 ) ) f 1 f 2 f 3 sin θ 10 ( cos θ 20 sin θ 30 σ 2 + cos θ 30 sin θ 20 σ 3 ) ] + 1 4 f 1 f 3 h 20 ν 1 σ 3 × cos θ 10 cos θ 30 ( 1 + σ 1 ) ] + 1 8 f 2 f 3 h 10 ω 1 c 3 cos θ 20 sin θ 30 [ 2 h 30 2 ν 2 h 20 2 ( ν 3 + h 10 2 ν 4 ) ]
+ 3 4 cos θ 20 cos θ 30 f 2 f 3 h 10 h 20 2 ν 2 ν 5 ω 1 + 1 4 f 2 f 3 h 10 ω 1 cos θ 20 cos θ 30 [ h 20 2 ( h 10 2 ν 4 ν 5 + h 20 2 ν 3 ν 5 + ν 3 ν 6 ) h 30 2 ν 2 ν 6 + h 10 2 ν 4 ν 6 ] 1 8 cos θ 20 sin θ 30 c 3 f 2 f 3 h 10 σ 2 ω 1 + 1 4 f 2 f 3 h 10 ω 1 cos θ 20 cos θ 30 [ h 20 2 ( σ 2 ν 5 + ν 3 σ 3 ) + ν 6 σ 2 + σ 3 ( h 10 ν 4 σ 2 h 30 2 ν 2 ) ] 1 8 c 3 f 1 f 2 h 30 ω 2 ( ν 1 cos θ 10 sin θ 20 + h 30 2 v 2 cos θ 20 sin θ 10 ) 1 8 f 1 f 2 ω 2 cos θ 20 [ c 3 sin θ 10 × ( h 10 2 h 30 v 4 + h 20 2 h 30 v 3 ) + 2 v 1 h 30 cos θ 10 ( h 30 2 v 2 h 20 2 v 3 h 10 2 v 4 ) ] 1 8 cos θ 10 sin θ 20 × c 3 f 1 f 2 h 30 σ 1 ω 2 1 4 f 1 f 2 ω 2 cos θ 10 cos θ 20 [ h 30 σ 1 ( h 30 2 v 2 + h 20 2 v 3 + h 10 2 v 4 ) + h 30 v 1 σ 2 ] + 1 32 c 3 f 1 cos θ 20 [ σ 2 ω 2 f 2 h 30 sin θ 10 + 4 f 3 h 10 ω 1 ( 2 h 30 2 ν 2 h 20 2 ν 3 h 10 2 h 20 2 ν 4 ) sin θ 30 ]
+ 1 16 f 2 [ 2 h 10 ω 1 f 3 ( 2 h 20 2 ν 2 ν 5 cos θ 20 cos θ 30 c 3 σ 2 cos θ 20 sin θ 30 ) sin 2 θ 10 c 3 f 1 h 30 ν 1 ω 2 ] + 1 4 f 2 f 3 h 10 ω 1 cos θ 20 cos θ 30 [ h 20 4 ν 5 ( ν 3 + h 10 2 ν 4 ) + ν 6 ( h 20 2 ν 3 h 30 2 ν 2 + h 10 2 ν 4 ) ] + 1 4 f 2 f 3 h 10 ω 1 cos θ 20 cos θ 30 [ h 20 2 ( σ 2 ν 5 + ν 3 σ 3 ) + σ 2 ( ν 6 σ 3 ) + σ 3 ( h 10 ν 4 h 30 2 ν 2 ) ] + 1 8 f 1 f 2 h 30 3 ω 2 v 2 cos θ 20 ( c 3 sin θ 10 2 v 1 cos θ 10 ) 1 8 f 1 f 2 h 30 ω 2 [ c 3 cos θ 20 sin θ 10 × ( h 20 2 v 3 h 10 2 v 4 σ 2 ) cos θ 10 sin θ 20 c 3 σ 1 + cos θ 20 h 20 2 v 1 v 3 ] + 1 4 f 1 f 2 ω 2 h 30 cos θ 10 × cos θ 20 [ h 10 2 v 4 ( v 1 + σ 1 ) + σ 1 ( h 30 2 v 2 + h 20 2 v 3 ) + σ 2 ( v 1 + σ 1 ) ] + c 2 ( sin θ 10 f 2 f 3 h 10 × ( 1 16 sin θ 30 c 3 1 8 ω 1 f 1 cos θ 30 ( h 20 2 ν 5 + ν 6 + σ 3 ) ( f 2 ( 1 8 cos θ 10 sin θ 20 sin θ 30 f 3 + h 30 ( 1 16 sin 2 θ 10 c 3 + sin 2 θ 10 [ 1 16 ( ν 1 + σ 1 ) ] ) ω 2 ) + h 30 ( f 3 [ 1 8 ( cos θ 10 sin θ 30 ν 1 + cos θ 30 sin θ 10 h 20 2 ν 5 ] + c 3 ( 1 16 sin θ 10 sin θ 30 f 3 + cos θ 10 h 30 [ 1 8 ( ν 1 + σ 1 ) ] ω 2 ) ) ) ) + c 1 ( sin θ 10 f 1 ( f 2 ( 1 16 sin θ 10 sin θ 30 f 3 + h 30 [ 1 16 ( sin θ 10 c 2 + sin θ 10 c 3 ) ] 1 8 cos θ 20 ( h 20 2 ν 3 + h 10 2 ν 4 + σ 2 h 30 2 ν 2 ) ) ω 2 ) 1 8 h 20 ( f 3 ( 1 16 sin θ 30 c 3 cos θ 30 ( h 20 2 ν 5 + ν 6 + σ 3 ) ) + 1 16 c 2 × ( sin θ 30 f 3 + c 3 h 30 ω 2 ) ) ) + h 10 ω 1 ( c 2 ( f 3 h 30 ( 1 16 sin θ 30 c 3 1 8 cos θ 30 ( h 20 2 ν 5 + ν 6 + σ 3 ) )
+ 1 16 sin θ 10 f 2 ( sin θ 30 f 3 + h 30 ω 2 c 3 ) ) + f 2 ( f 3 ( 1 8 ( cos θ 20 sin θ 30 h 30 2 ν 2 cos θ 20 sin θ 30 h 10 2 ν 4 h 20 2 ( cos θ 20 sin θ 30 ν 3 + cos θ 30 sin θ 10 ν 5 ) cos θ 30 sin θ 10 ν 6 cos θ 20 sin θ 30 σ 2 cos θ 30 × sin θ 10 σ 3 ) ) + c 3 ( 1 16 sin θ 10 sin θ 30 f 3 + cos θ 20 h 30 ( 1 8 ( h 30 2 ν 2 h 20 2 ν 3 h 10 2 ν 4 σ 2 ) ω 2 ) ) ) ) ) } ,
Γ 5 = 1 h 10 h 20 h 30 ω 1 ω 1 { f 1 ( f 2 ( f 3 ( 1 8 [ cos θ 20 ( cos θ 10 sin θ 30 h 10 2 v 1 v 4 + cos θ 30 sin θ 10 h 20 4 v 3 v 5 ) + cos θ 10 ( cos θ 30 sin θ 20 v 1 v 6 + h 10 2 σ 1 v 4 cos θ 20 sin θ 30 ) + v 6 cos θ 30 ( h 10 2 v 4 cos θ 20 sin θ 10 + σ 1 cos θ 10 sin θ 20 ) + cos θ 10 cos θ 20 sin θ 30 v 1 σ 2 + cos θ 20 cos θ 30 sin θ 10 σ 2 v 6 + cos θ 10 ( σ 1 σ 2 cos θ 20 sin θ 30 + v 1 σ 3 cos θ 30 sin θ 20 ) + σ 3 cos θ 30 ( h 10 2 v 4 cos θ 20 sin θ 10 + σ 1 cos θ 30 sin θ 20 ) + σ 2 σ 3 cos θ 20 cos θ 30 sin θ 10 ] + cos θ 20 h 30 2 v 2 ( 1 8 v 1 cos θ 10 × sin θ 30 + 15 40 cos θ 30 sin θ 10 h 20 2 v 5 1 8 [ cos θ 30 ( sin θ 10 v 6 σ 3 sin θ 10 ) σ 1 cos θ 10 sin θ 30 + h 20 2 ( cos θ 20 cos θ 30 sin θ 10 h 10 2 v 4 v 5 + cos θ 10 v 1 ( cos θ 20 sin θ 30 v 3 + cos θ 30 sin θ 20 v 5 ) + v 3 cos θ 20 ( v 6 cos θ 30 sin θ 10 + σ 1 cos θ 10 sin θ 30 ) + cos θ 10 cos θ 30 sin θ 20 v 5 σ 1 + cos θ 20 cos θ 30 ( v 5 σ 2 sin θ 10 + v 3 σ 3 sin θ 10 ) ) ) ] + c 3 ( sin θ 30 f 3 [ 1 16 ( cos θ 10 sin θ 20 v 1 + cos θ 20 sin θ 10 ( h 30 2 v 2 h 20 2 v 3 h 10 2 v 4 σ 2 ) cos θ 10 sin θ 20 σ 1 + cos θ 10 cos θ 20 h 30 × ( 1 8 [ h 10 v 1 v 4 h 30 2 v 2 ( v 1 + σ 1 ) + h 20 2 v 3 ( v 1 + σ 1 ) + h 10 2 v 4 σ 1 + v 1 σ 2 + σ 1 σ 2 ) ω 2 ) ) + c 2 ( f 3 ( cos θ 10 cos θ 30 h 20 ( h 20 2 v 5 ( v 1 + σ 1 ) + ( v 1 + σ 1 ) ( v 6 + σ 3 ) ) ] + f 2 sin θ 20 × [ 1 16 ( cos θ 10 ( v 1 sin θ 30 + σ 1 sin θ 30 ) sin θ 10 h 20 2 v 5 ( cos θ 30 + v 6 cos θ 30 ) cos θ 30 sin θ 10 σ 3 ) ) c 3 ( cos θ 10 sin θ 30 f 3 θ 20 ( v 1 + σ 1 ) ] + sin θ 20 f 2 × ( 1 32 sin θ 10 sin θ 30 f 3 1 16 cos θ 10 h 30 ( v 1 + σ 1 ) ω 2 ) ) ) ) + c 1 ( c 2 ( sin θ 20 f 2 f 3 h 10 × 1 32 ( sin θ 30 c 3 2 cos θ 30 ( h 20 2 v 5 + v 5 + σ 3 ) ) ω 1 + sin θ 10 f 1 ( 1 32 f 3 h 20 ( sin θ 30 c 3 2 cos θ 30 ( h 20 2 v + v 6 + σ 3 ) ) + 1 32 sin θ 20 f 2 ( sin θ 30 f 3 + c 3 h 30 ω 2 ) ) ) + f 2 ( f 3 ( sin θ 10 f 1 [ 1 16 ( cos θ 20 sin θ 30 ( h 30 2 v 2 h 10 2 v 4 ) h 20 2 ( v 3 cos θ 20 sin θ 30 + v 5 cos θ 30 sin θ 20 ) cos θ 30 sin θ 20 v 6 cos θ 20 sin θ 30 σ 2 ] 1 16 σ 3 h 10 cos θ 30 × sin θ 20 cos θ 20 [ 1 8 ( h 10 2 v 4 v 6 + v 6 σ 2 h 30 2 v 2 ( v 1 + σ 1 ) h 10 2 v 4 σ 3 + σ 2 σ 3
+ h 20 2 ( h 30 2 v 2 v 5 + h 10 2 v 4 v 5 + v 3 v 6 + σ 2 v 5 + v 3 σ 3 ) ] ) ω 1 ) + 1 16 c 3 ( f 3 h 10 ω 1 cos θ 20 sin θ 30 × ( h 30 2 v 2 h 20 2 v 3 h 10 2 v 4 σ 2 ) + ( 1 32 sin θ 10 f 1 ( 1 32 sin θ 20 sin θ 30 f 3 + cos θ 20 h 30 × [ 1 16 ( h 30 2 v 2 h 20 2 v 3 h 10 2 v 4 σ 2 ) ] ω 2 ) ) ) ) } ,
Γ 6 = 1 h 10 h 20 h 30 ω 1 ω 2 f 1 f 2 f 3 { c 2 sin θ 10 sin θ 20 ( 0.01652 c 3 sin θ 30 0.03152 cos θ 30 ( ν 5 h 20 2 + ν 6 + σ 3 ) ) + cos θ 20 ( sin θ 30 c 3 [ 0.03125 ( ν 2 h 30 2 ν 3 h 20 2 ν 4 h 10 2 σ 2 ) ] + 0.0625 cos θ 30 × ( 0.0625 ν 3 ν 5 h 20 2 σ 2 ν 2 h 30 2 ( ν 6 + σ 3 ) + ν 4 h 10 2 σ 3 + σ 2 σ 3 + 0.0625 h 20 2 ν 2 ν 5 h 30 2 + ν 4 ν 5 h 10 2 + ν 3 ν 6 + ν 5 σ 2 + ν 3 σ 3 ) ) ) ) + cos θ 30 ( sin θ 20 c 2 ( sin θ 30 c 3 [ 0.03125 ( ν 1 + σ 1 ) ] + cos θ 30 ( ν 5 h 20 2 [ 0.0625 ( ν 1 + σ 1 ) ] + 0.0625 ( ν 1 + σ 1 ) ( ν 6 + σ 3 ) ) ) + cos θ 20 ( sin θ 30 c 3 × ( 0.0625 ν 1 ν 4 h 10 2 0.0625 ( ν 2 h 30 2 ( ν 1 + σ 1 ) + ν 4 h 10 2 σ 1 + ν 1 σ 2 + σ 1 σ 2 ) ) ) ) } .

References

  1. El-Barki, F.A.; Ismail, A.I.; Shaker, M.O.; Amer, T.S. On the motion of the pendulum on an ellipse. ZAMM 1999, 79, 65–72. [Google Scholar] [CrossRef]
  2. Leven, R.W.; Koch, B.P. Chaotic behavior of a parametrically excited damped pendulum. Phys. Lett. A 1981, 86, 71–74. [Google Scholar] [CrossRef]
  3. Leven, R.W.; Pompe, B.; Wilke, C.; Koch, B.P. Experiments on periodic and chaotic motions of a parametrically forced pendulum. Phys. D Nonlinear Phenom. 1985, 16, 371–384. [Google Scholar] [CrossRef]
  4. Amer, T.S.; Bek, M.A. Chaotic responses of a harmonically excited spring pendulum moving in circular path. Nonlinear Anal. Real World Appl. 2009, 10, 3196–3202. [Google Scholar] [CrossRef]
  5. Hatwal, H.; Mallik, A.; Ghosh, A. Non-linear vibrations of a harmonically excited autoparametric system. J. Sound Vib. 1982, 81, 153–164. [Google Scholar] [CrossRef]
  6. Nayfeh, A.H. Perturbations Methods; WILEY-VCH Verlag GmbH and Co. KGaA: Weinheim, Germany, 2004. Available online: https://onlinelibrary.wiley.com/doi/book/10.1002/9783527617609 (accessed on 19 December 2021).
  7. Miles, J. Resonance and symmetry breaking for the pendulum. Phys. D Nonlinear Phenom. 1988, 16, 252–268. [Google Scholar] [CrossRef]
  8. Bryant, P.J.; Miles, J.W. On a periodically forced, weakly damped pendulum. Part 1: Applied torque. J. Austral. Math. Soc. Ser. B 1990, 32, 1–22. [Google Scholar] [CrossRef] [Green Version]
  9. Bryant, P.J.; Miles, J.W. On a periodically forced, weakly damped pendulum. Part 2: Horizontal forcing. J. Austral. Math. Soc. Ser. B 1990, 32, 23–41. [Google Scholar] [CrossRef] [Green Version]
  10. Bryant, P.J.; Miles, J.W. On a periodically forced, weakly damped pendulum. Part 3: Vertical forcing. J. Austral. Math. Soc. Ser. B 1990, 32, 42–60. [Google Scholar] [CrossRef] [Green Version]
  11. Kapitaniak, M.; Perlikowski, P.; Kapitaniak, T. Synchronous motion of two vertically excited planar elastic pendula. Commun. Nonlinear Sci. Numer. Simulat. 2013, 18, 2088–2096. [Google Scholar] [CrossRef] [Green Version]
  12. Warminski, J.; Kecik, K. Autoparametric vibrations of a nonlinear system with pendulum. Math. Problems Eng. 2006, 2006, 19. [Google Scholar] [CrossRef]
  13. Starosta, R.; Kamińska, G.S.; Awrejcewicz, J. Parametric and external resonances in kinematically and externally excited nonlinear spring pendulum. Int. J. Bifurc. Chaos 2011, 21, 3013–3021. [Google Scholar] [CrossRef]
  14. Starosta, R.; Kamińska, G.S.; Awrejcewicz, J. Asymptotic analysis of kinematically excited dynamical systems near resonances. Nonlinear Dyn. 2012, 68, 459–469. [Google Scholar] [CrossRef]
  15. Amer, T.S.; Bek, M.A.; Hamada, I. On the motion of harmonically excited spring pendulum in elliptic path near resonances. Adv. Math. Phys. 2016, 2016, 15. [Google Scholar] [CrossRef] [Green Version]
  16. Awrejcewicz, J.; Starosta, R.; Kamińska, G.S. Asymptotic analysis of resonances in nonlinear vibrations of the 3-DOF pendulum. Differ. Equ. Dyn. Syst. 2013, 21, 123–140. [Google Scholar] [CrossRef]
  17. Amer, T.S.; Bek, M.A.; Abouhmr, M.K. On the vibrational analysis for the motion of a harmonically damped rigid body pendulum. Nonlinear Dyn. 2018, 91, 2485–2502. [Google Scholar] [CrossRef]
  18. Amer, T.S.; Bek, M.A.; Abohamer, M.K. On the motion of a harmonically excited damped spring pendulum in an elliptic path. Mech. Res. Commun. 2019, 95, 23–34. [Google Scholar] [CrossRef]
  19. El-Sabaa, F.M.; Amer, T.S.; Gad, H.M.; Bek, M.A. On the motion of a damped rigid body near resonances under the influence of harmonically external force and moments. Results Phys. 2020, 19, 103352. [Google Scholar] [CrossRef]
  20. Bek, M.A.; Amer, T.S.; Almahalawy, A.; Elameer, A.S. The asymptotic analysis for the motion of 3DOF dynamical system close to resonances. Alex. Eng. J. 2021, 60, 3539–3551. [Google Scholar] [CrossRef]
  21. Amer, T.S.; Bek, M.A.; Hassan, S.S. The dynamical analysis for the motion of a harmonically two degrees of freedom damped spring pendulum in an elliptic trajectory. Alex. Eng. J. 2022, 2, 1715–1733. [Google Scholar] [CrossRef]
  22. Kamińska, G.S.; Awrejcewicz, J.; Kamiński, H. Resonance study of spring pendulum based on asymptotic solutions with polynomial approximation in quadratic means. Meccanica 2020, 56, 753–767. [Google Scholar]
  23. Bek, M.A.; Amer, T.S.; Sirwah, M.A.; Awrejcewicz, J.; Arab, A.A. The vibrational motion of a spring pendulum in a fluid flow. Results Phys. 2020, 19, 103465. [Google Scholar] [CrossRef]
  24. Amer, W.S.; Bek, M.A.; Abohamer, M.K. On the motion of a pendulum attached with tuned absorber near resonances. Results Phys. 2018, 11, 291–301. [Google Scholar] [CrossRef]
  25. Song, Y.; Sato, H.; Iwata, Y.; Komatsuzaki, T. The response of a dynamic vibration absorber system with a parametrically excited pendulum. J. Sound Vib. 2003, 259, 747–759. [Google Scholar] [CrossRef]
  26. Warminski, J.; Kecik, K. Instabilities in the main parametric resonance area of a mechanical system with a pendulum. J. Sound Vib. 2009, 322, 612–628. [Google Scholar] [CrossRef]
  27. Warminski, J. Regular and chaotic vibrations of a parametrically and self-excited system under internal resonance condition. Meccanica 2005, 40, 181–202. [Google Scholar] [CrossRef]
  28. Kecik, K.; Warminski, J. Dynamics of an autoparametric pendulum-like system with a nonlinear semiactive suspension. Math. Probl. Eng. 2011, 2011, 15. [Google Scholar] [CrossRef] [Green Version]
  29. Kecik, K.; Mitura, A.; Warmiński, J. Efficiency analysis of an autoparametric pendulum vibration absorber. Eksploat. I Niezawodn.-Maint. Reliab. 2013, 15, 221–224. [Google Scholar]
  30. Kecik, K.; Mitura, A. Energy recovery from a pendulum tuned mass damper with two independent harvesting sources. Int. J. Mech. Sci. 2020, 174, 105568. [Google Scholar] [CrossRef]
  31. Brzeski, P.; Perlikowski, P.; Yanchuk, S.; Kapitaniak, T. The dynamics of the pendulum suspended on the forced Duffing oscillator. J. Sound Vib. 2012, 331, 5347–5357. [Google Scholar] [CrossRef] [Green Version]
  32. Gantmacher, F.R. Applications of the Theory of Matrices; John Wiley & Sons: New York, NY, USA, 2005. [Google Scholar]
  33. Amer, T.S.; Bek, M.A.; Hassan, S.S.; Elbendary, S. The stability analysis for the motion of a nonlinear damped vibrating dynamical system with three-degrees-of-freedom. Results Phys. 2021, 28, 104561. [Google Scholar] [CrossRef]
  34. Abady, I.M.; Amer, T.S.; Gad, H.M.; Bek, M.A. The asymptotic analysis and stability of 3DOF non-linear damped rigid body pendulum near resonance. Ain Shams Eng. J. 2022, 13, 101554. [Google Scholar] [CrossRef]
  35. Abohamer, M.K.; Awrejcewicz, J.; Starosta, R.; Amer, T.S.; Bek, M.A. Influence of the motion of a spring pendulum on energy-harvesting devices. Appl. Sci. 2021, 11, 8658. [Google Scholar] [CrossRef]
  36. Amer, W.S.; Amer, T.S.; Hassan, S.S. Modeling and stability analysis for the vibrating motion of three degrees-of-freedom dynamical system near resonance. Appl. Sci. 2021, 11, 11943. [Google Scholar] [CrossRef]
  37. Amer, W.S.; Amer, T.S.; Starosta, R.; Bek, A.M. Resonance in the cart-pendulum system-an asymptotic approach. Appl. Sci. 2021, 11, 11567. [Google Scholar] [CrossRef]
  38. Amer, T.S.; Starosta, R.; Elameer, A.S.; Bek, M.A. Analyzing the stability for the motion of an unstretched double pendulum near resonance. Appl. Sci. 2021, 11, 9520. [Google Scholar] [CrossRef]
  39. Abdelhfeez, S.A.; Amer, T.S.; Elbaz, R.F.; Bek, M.A. Studying the influence of external torques on the dynamical motion and the stability of a 3DOF dynamic system. Alex. Eng. J. 2022, 61, 6695–6724. [Google Scholar] [CrossRef]
Figure 1. Description of the system.
Figure 1. Description of the system.
Applsci 12 01737 g001
Figure 2. The amplitude’s modulation h 1 and modified phase θ 1 : (a,b) at ω ¯ 1 = ( 3.02 , 3.43 , 4.08 ) , (c,d) at ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) , (e,f) at ω ¯ 3 = ( 7.07 , 7.75 , 8.37 ) .
Figure 2. The amplitude’s modulation h 1 and modified phase θ 1 : (a,b) at ω ¯ 1 = ( 3.02 , 3.43 , 4.08 ) , (c,d) at ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) , (e,f) at ω ¯ 3 = ( 7.07 , 7.75 , 8.37 ) .
Applsci 12 01737 g002
Figure 3. The variation of h 2 and θ 2 with τ : (a,b) at ω ¯ 1 = ( 3.02 , 3.43 , 4.08 ) , (c,d) at ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) , (e,f) at ω ¯ 3 = ( 7.07 , 7.75 , 8.37 ) .
Figure 3. The variation of h 2 and θ 2 with τ : (a,b) at ω ¯ 1 = ( 3.02 , 3.43 , 4.08 ) , (c,d) at ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) , (e,f) at ω ¯ 3 = ( 7.07 , 7.75 , 8.37 ) .
Applsci 12 01737 g003
Figure 4. The time histories of the amplitude h 3 and the phase θ 3 : (a,b) at ω ¯ 1 = ( 3.02 , 3.43 , 4.08 ) , (c,d) at ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) , (e,f) at ω ¯ 3 = ( 7.07 , 7.75 , 8.37 ) .
Figure 4. The time histories of the amplitude h 3 and the phase θ 3 : (a,b) at ω ¯ 1 = ( 3.02 , 3.43 , 4.08 ) , (c,d) at ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) , (e,f) at ω ¯ 3 = ( 7.07 , 7.75 , 8.37 ) .
Applsci 12 01737 g004
Figure 5. The change of the attained Y , X , and ϕ with time τ : (ac) when ω ¯ 1 = ( 3.02 , 3.43 , 4.08 ) .
Figure 5. The change of the attained Y , X , and ϕ with time τ : (ac) when ω ¯ 1 = ( 3.02 , 3.43 , 4.08 ) .
Applsci 12 01737 g005
Figure 6. The time histories of the approximate solutions Y , X , and ϕ : (ac) when ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) .
Figure 6. The time histories of the approximate solutions Y , X , and ϕ : (ac) when ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) .
Applsci 12 01737 g006
Figure 7. The variation of Y , X , and ϕ versus τ : (ac) when ω ¯ 3 = ( 7.07 , 7.75 , 8.37 ) .
Figure 7. The variation of Y , X , and ϕ versus τ : (ac) when ω ¯ 3 = ( 7.07 , 7.75 , 8.37 ) .
Applsci 12 01737 g007
Figure 8. The amplitude modulation of h 1 , h 2 , h 3 as a function of θ 1 , θ 2 , θ 3 : (ac) at ω ¯ 1 = 4.08 .
Figure 8. The amplitude modulation of h 1 , h 2 , h 3 as a function of θ 1 , θ 2 , θ 3 : (ac) at ω ¯ 1 = 4.08 .
Applsci 12 01737 g008
Figure 9. The amplitude modulation of h 1 , h 2 , h 3 as a function of θ 1 , θ 2 , θ 3 : (ac) at ω ¯ 2 = 5.72 .
Figure 9. The amplitude modulation of h 1 , h 2 , h 3 as a function of θ 1 , θ 2 , θ 3 : (ac) at ω ¯ 2 = 5.72 .
Applsci 12 01737 g009
Figure 10. The amplitude modulation of h 1 , h 2 , h 3 as a function of θ 1 , θ 2 , θ 3 : (ac) at ω ¯ 3 = 8.37 .
Figure 10. The amplitude modulation of h 1 , h 2 , h 3 as a function of θ 1 , θ 2 , θ 3 : (ac) at ω ¯ 3 = 8.37 .
Applsci 12 01737 g010
Figure 11. The change of the frequency responses of h 1 , h 2 , h 3 with σ 1 at σ 2 = 0 , σ 3 = 0 : (ac) when ω ¯ 1 = ( 3.02 , 3.43 , 4.08 ) .
Figure 11. The change of the frequency responses of h 1 , h 2 , h 3 with σ 1 at σ 2 = 0 , σ 3 = 0 : (ac) when ω ¯ 1 = ( 3.02 , 3.43 , 4.08 ) .
Applsci 12 01737 g011
Figure 12. The variation of the frequency response h 1 with σ 1 at σ 2 = 0 , σ 3 = 0 when ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) .
Figure 12. The variation of the frequency response h 1 with σ 1 at σ 2 = 0 , σ 3 = 0 when ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) .
Applsci 12 01737 g012
Figure 13. The variation of the frequency response h 2 with σ 1 at σ 2 = 0 , σ 3 = 0 when ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) .
Figure 13. The variation of the frequency response h 2 with σ 1 at σ 2 = 0 , σ 3 = 0 when ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) .
Applsci 12 01737 g013
Figure 14. The variation of the frequency responses h 3 with σ 1 at σ 2 = 0 , σ 3 = 0 when ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) .
Figure 14. The variation of the frequency responses h 3 with σ 1 at σ 2 = 0 , σ 3 = 0 when ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) .
Applsci 12 01737 g014
Figure 15. The modified amplitudes via τ and the projections of the modulation equation tracks on the phase plane u 1 v 1 : (ac) when ω ¯ 1 = ( 3.02 , 3.43 , 4.08 ) .
Figure 15. The modified amplitudes via τ and the projections of the modulation equation tracks on the phase plane u 1 v 1 : (ac) when ω ¯ 1 = ( 3.02 , 3.43 , 4.08 ) .
Applsci 12 01737 g015
Figure 16. The modified amplitudes via τ and the projections of the modulation equation tracks on the phase plane u 2 v 2 : (ac) when ω ¯ 1 = ( 3.02 , 3.43 , 4.08 ) .
Figure 16. The modified amplitudes via τ and the projections of the modulation equation tracks on the phase plane u 2 v 2 : (ac) when ω ¯ 1 = ( 3.02 , 3.43 , 4.08 ) .
Applsci 12 01737 g016
Figure 17. The modified amplitudes via τ and the projections of the modulation equation tracks on the phase plane u 3 v 3 : (ac) when ω ¯ 1 = ( 3.02 , 3.43 , 4.08 ) .
Figure 17. The modified amplitudes via τ and the projections of the modulation equation tracks on the phase plane u 3 v 3 : (ac) when ω ¯ 1 = ( 3.02 , 3.43 , 4.08 ) .
Applsci 12 01737 g017
Figure 18. The modified amplitudes via τ and the projections of the modulation equation tracks on the phase plane u 1 v 1 : (ac) when ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) .
Figure 18. The modified amplitudes via τ and the projections of the modulation equation tracks on the phase plane u 1 v 1 : (ac) when ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) .
Applsci 12 01737 g018
Figure 19. The modified amplitudes via τ and the projections of the modulation equation tracks on the phase plane u 2 v 2 : (ac) when ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) .
Figure 19. The modified amplitudes via τ and the projections of the modulation equation tracks on the phase plane u 2 v 2 : (ac) when ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) .
Applsci 12 01737 g019
Figure 20. The modified amplitudes via τ and the projections of the modulation equation tracks on the phase plane u 3 v 3 : (ac) when ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) .
Figure 20. The modified amplitudes via τ and the projections of the modulation equation tracks on the phase plane u 3 v 3 : (ac) when ω ¯ 2 = ( 4.43 , 5.72 , 6.68 ) .
Applsci 12 01737 g020
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Amer, T.S.; Starosta, R.; Almahalawy, A.; Elameer, A.S. The Stability Analysis of a Vibrating Auto-Parametric Dynamical System Near Resonance. Appl. Sci. 2022, 12, 1737. https://doi.org/10.3390/app12031737

AMA Style

Amer TS, Starosta R, Almahalawy A, Elameer AS. The Stability Analysis of a Vibrating Auto-Parametric Dynamical System Near Resonance. Applied Sciences. 2022; 12(3):1737. https://doi.org/10.3390/app12031737

Chicago/Turabian Style

Amer, Tarek S., Roman Starosta, Ashraf Almahalawy, and Abdelkarim S. Elameer. 2022. "The Stability Analysis of a Vibrating Auto-Parametric Dynamical System Near Resonance" Applied Sciences 12, no. 3: 1737. https://doi.org/10.3390/app12031737

APA Style

Amer, T. S., Starosta, R., Almahalawy, A., & Elameer, A. S. (2022). The Stability Analysis of a Vibrating Auto-Parametric Dynamical System Near Resonance. Applied Sciences, 12(3), 1737. https://doi.org/10.3390/app12031737

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop