Next Article in Journal
Algebraic Numbers as Product of Powers of Transcendental Numbers
Next Article in Special Issue
Application of Differential Transform to Multi-Term Fractional Differential Equations with Non-Commensurate Orders
Previous Article in Journal
On Some Formulas for Kaprekar Constants
Previous Article in Special Issue
Periodic Solution of the Strongly Nonlinear Asymmetry System with the Dynamic Frequency Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Periodic Solution of a Non-Smooth Double Pendulum with Unilateral Rigid Constrain

1
Department of Mathematics and Information Science, Hebei Normal University, Shijiazhuang 050000, China
2
Department of Mathematics and Physics, Shijiazhuang Tiedao University, Shijiazhuan 050043, China
3
Department of Mechanics Engineering, Shijiazhuang Tiedao University, Shijiazhuang 050043, China
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(7), 886; https://doi.org/10.3390/sym11070886
Submission received: 17 June 2019 / Revised: 1 July 2019 / Accepted: 3 July 2019 / Published: 6 July 2019
(This article belongs to the Special Issue Nonlinear Oscillations and Boundary Value Problems)

Abstract

:
In this paper, a double pendulum model is presented with unilateral rigid constraint under harmonic excitation, which leads to be an asymmetric and non-smooth system. By introducing impact recovery matrix, modal analysis, and matrix theory, the analytical expressions of the periodic solutions for unilateral double-collision will be discussed in high-dimensional non-smooth asymmetric system. Firstly, the impact laws are classified in order to detect the existence of periodic solutions of the system. The impact recovery matrix is introduced to transform the impact laws of high-dimensional system into matrix. Furthermore, by use of modal analysis and matrix theory, an invertible transformation is constructed to obtain the parameter conditions for the existence of the impact periodic solution, which simplifies the calculation and can be easily extended to high-dimensional non-smooth system. Hence, the range of physical parameters and the restitution coefficients is calculated theoretically and non-smooth analytic expression of the periodic solution is given, which provides ideas for the study of approximate analytical solutions of high-dimensional non-smooth system. Finally, numerical simulation is carried out to obtain the impact periodic solution of the system with small angle motion.

1. Introduction

As the most important actuator among all robot mechanisms, the mechanical arm was an important research subject of robot technology [1,2,3]. The researches were focused on the drive system, sensing system, shape design of the mechanical arm, and inertial impact influence on the dynamic response at the link [3,4]. Many results showed that the impact double pendulum could be used to simulate the bionic manipulator [5,6,7]. In fact, with the motion refinement of the manipulator and the installation requirements of the external drive, the gap and damping at the link of the manipulator should be considered, as well as the type and installation mode of the external drive, which can be simplified as a collision double pendulum with external simple harmonic excitation. The collision double pendulum with external excitation can explore the dynamic characteristics of the joint of the external wearable manipulator. Specifically, the non-smooth periodic motion of the system is a hot issue. Through theoretical research, suitable physical and geometric parameters will be provided for the design of the manipulator to improve the application comfort, safety, and service life of the manipulator.
Impact pendulum system is a typical non-smooth dynamic system. In low-dimensional systems most studies focused on the periodic solution, bifurcation, chaos, and so on [8,9,10,11,12]. For non-smooth high-dimensional system, numerical simulation was applied to detect the dynamics [13,14,15,16,17,18]. The oblique impact vibration between the double pendulum and the unilateral rigid constraint surface was studied by numerical method [14]. The influence of excitation parameters and system physical parameters on the steady-state behavior of the system was given. Although the numerical solutions of unilateral constrained systems were discussed, the analytic form and existence conditions of periodic solutions still needed further detection. According to the periodic solution, there are many results of analytical methods. Furthermore, the existence conditions of periodic solutions of single-degree-of-freedom collision systems with unilateral single-impact and bilateral single-impact were obtained in [19,20], which provided a theoretical analytical method for studying periodic solutions of impact systems. However, its expression was so complex that it is difficult to generalize to practical application engineering. In terms of theoretical analysis, it is also difficult to extend the calculation method to multiple degrees of freedom. Hence, in this paper, by introducing impact recovery matrix and coupling the modal analysis and matrix theory, the analytical expressions of the periodic solutions for unilateral double-impact will be discussed in high-dimensional non-smooth system.
This paper is organized as follows. In the next section, a non-smooth double pendulum model with unilateral double-impact is constructed. In Section 3, the types of periodic solutions of two-degree-of-freedom multi-point collisions are studied due to the different boundary conditions. By using the matrix theory [21] and introducing the invertible transformation, the theoretical conditions for the existence of the periodic solutions of collisions and the analytic expressions of the periodic solutions of collisions are discussed in Section 4. Numerical simulations are also carried out to confirm these analytical predictions in Section 5.

2. System Modeling

Based on the double pendulum, we constructed a physical model of the non-smooth double pendulum with unilateral rigid constraint as shown in Figure 1. The model consists of a smooth hinged double pendulum with a rigid wall fixed to the base with a horizontal simple harmonic excitation A cos ( Ω t ) . A , Ω are the amplitude and frequency of the horizontal simple harmonic excitation. Here, only consider the small angular motion of the two pendulums in the vertical plane, which means θ 0 θ * θ 1 θ 0 and θ 0 θ * θ 2 θ 0 ( θ * is a smaller positive number). Assume that the angular velocity is positive in the counter clockwise direction.
Where l 1 , l 2 ; m 1 , m 2 ; c 1 , c 2 are the lengths of the two pendulums, the mass of two pendulums, and the damper of the hinge joint, respectively. Assuming the two pendulums can only swing at a small angle is considered and the elastic collision of the pendulum occurs at the rigid wall θ i = θ 0 , i = 1 , 2 . That is to say, the deformation and time during the pulse excitation are ignored.
When the pendulum collides with the rigid wall, the vector field of the system will have a sudden change in momentum. The phase diagram and time history diagram will be asymmetrical discontinuity, as shown in Figure 2a,b.
In the system, the pendulum angles θ 1 and θ 2 of the double pendulum are taken as generalized coordinates. According to Lagrange function method, the equation of the system can be obtained
{ ( m 1 + m 2 ) l 1 2 θ ¨ 1 + m 2 l 1 l 2 θ ¨ 2 cos ( θ 1 θ 2 ) + m 2 l 1 l 2 θ ˙ 2 2 sin ( θ 1 θ 2 ) + ( m 1 + m 2 ) g l 1 sin θ 1 A Ω 2 ( m 1 + m 2 ) l 1 cos ( Ω t ) cos θ 1 = c 1 θ ˙ 1 c 2 ( θ ˙ 1 θ ˙ 2 ) θ ˙ 1 + = ( 1 γ ) θ ˙ 1 | θ = θ 0 π 2 < θ 1 < θ 0 m 2 l 2 2 θ ¨ 2 + m 2 l 1 l 2 θ ¨ 1 cos ( θ 1 θ 2 ) m 2 l 1 l 2 θ ˙ 1 2 sin ( θ 1 θ 2 ) + m 2 g l 1 sin θ 2 A Ω 2 m 2 l 2 cos ( Ω t ) cos θ 2 = c 2 ( θ ˙ 1 θ ˙ 2 ) θ ˙ 2 + = ( 1 γ ) θ ˙ 2 | θ = θ 0 π 2 < θ 2 < θ 0   ,
where θ 0 [ 0 , π 2 ] , θ ˙ and θ ˙ + represents the angular velocity before and after the impact, respectively. γ ( 0 < γ < 1 ) is a coefficient of restitution and ( 1 γ ) ( 0 , 1 ) .
Let
m = m 2 m 1 , ξ = m 1 + m , ω 1 2 = g l 1 , ω 2 2 = g l 2 , μ = l 1 l 2 , α = ξ μ 1 , a = A l 1 , ω τ = Ω t , p 1 = c 1 ω 1 ( 1 + m ) 1 2 ( m 1 + m 2 ) l 1 2 , p 2 = c 2 ω 1 ( 1 + m ) 1 2 ( m 1 + m 2 ) l 1 2 , p 3 = c 2 ω 1 ( 1 + m ) 1 2 m 2 l 2 2
The non-dimensional Equation (1) is written as
{ θ 1 + α θ 2 cos ( θ 1 θ 2 ) + α θ 2 2 sin ( θ 1 θ 2 ) + 1 1 + m sin θ 1 + p 1 θ 1 + p 2 ( θ 1 θ 2 ) a ω 2 cos ( ω τ ) cos θ 1 = 0 θ 1 + = ( 1 γ ) θ 1 | θ = θ 0 π 2 < θ 1 < θ 0 θ 2 + μ θ 1 cos ( θ 1 θ 2 ) μ θ 1 2 sin ( θ 1 θ 2 ) + μ 1 + m sin θ 2 + p 3 ( θ 2 θ 1 ) μ a ω 2 cos ( ω τ ) cos θ 2 = 0 θ 2 + = ( 1 γ ) θ 2 | θ = θ 0 π 2 < θ 2 < θ 0 ,
where θ i is a time derivation d θ i d τ , i = 1 , 2 .

3. Classification of Periodic Solution

Take t = 0 ( θ 1 = θ 2 = θ 0 ) as the initial time and T 0 = 2 n π ω ( θ 1 = θ 2 = θ 0 ) as the instantaneous of the next collision. This periodic motion, i.e., [ 0 , T 0 ] , is called a one-touch periodic motion.
s 1 : The top pendulum collides with the rigid constraint:
θ 1 ( 0 ) = θ 1 ( T 0 ) = θ 10 , θ 1 ( 0 ) = θ 1 ( T 0 + ) , θ 1 ( T 0 + ) = ( 1 γ ) θ 1 ( T 0 ) .
s 0 : The top pendulum did not collide with the rigid constraint:
θ 1 ( 0 ) = θ 1 ( T 0 ) = θ 10 , θ 1 ( 0 ) = θ 1 ( T 0 ) = θ 10 .
x 1 : The bottom pendulum collides with the rigid constraint:
θ 2 ( 0 ) = θ 2 ( T 0 ) = θ 20 , θ 2 ( 0 ) = θ 2 ( T 0 + ) , θ 2 ( T 0 + ) = ( 1 γ ) θ 2 ( T 0 ) .
x 0 : The bottom pendulum did not collide with the rigid constraint:
θ 2 ( 0 ) = θ 2 ( T 0 ) = θ 20 , θ 2 ( 0 ) = θ 2 ( T 0 ) = θ 20 .
Note that we consider periodic solutions, which implies θ 10 = θ 20 = θ 0 .
We consider the existence conditions of different types of periodic solutions. The classification of periodic solutions is shown as Table 1.
Because of the strong nonlinearity and non-smoothness of the system, it is difficult to find the periodic solution of the original system directly. It is especially difficult to find the periodic analytical solutions of different boundary conditions by selecting reasonable physical parameters and coefficient of restitution of the model. So, considering the small angle oscillation of the manipulator arm (leg) in engineering practice, we study the periodic solution of the small angle motion of the system. Because of non-smoothness, the system is still strongly nonlinear and non-smooth even though it is moving at a small angle. Hence, we will detect the conditions for the existence of parameters of periodic solutions of non-smooth systems and the reference parameters range for the existence of non-smooth periodic solutions can be given.

4. Periodic Solution

In order to facilitate the subsequent analysis, the following approximations
cos ( θ 1 θ 2 ) 1 , sin ( θ 1 θ 2 ) θ 1 θ 2 , sin θ i θ i , cos θ i 1 , i = 1 , 2 , θ 1 2 0 , θ 2 2 0
are given based on the fact that the swings of top and bottom pendulum are enough small during the whole motion process. Substituting (4) into (3), (3) is rewritten as:
{ θ 1 + α θ 2 + 1 1 + m θ 1 + p 1 θ 1 + p 2 ( θ 1 θ 2 ) a ω 2 cos ( ω τ ) = 0 θ 1 + = ( 1 γ ) θ 1 | θ = θ 0 π 2 < θ 1 < θ 0 θ 2 + μ θ 1 + μ 1 + m θ 2 + p 3 ( θ 2 θ 1 ) μ a ω 2 cos ( ω τ ) = 0 θ 2 + = ( 1 γ ) θ 2 | θ = θ 0 π 2 < θ 2 < θ 0 ,
i.e.,
{ θ 1 + θ 1 ξ θ 2 + k 1 θ 1 k 2 θ 2 = a ω 2 cos ( ω τ ) θ 1 + = ( 1 γ ) θ 1 | θ = θ 0 π 2 < θ 1 < θ 0   θ 2 μ θ 1 + μ θ 2 k 3 θ 1 + k 4 θ 2 = 0 θ 2 + = ( 1 γ ) θ 2 | θ = θ 0 π 2 < θ 2 < θ 0 ,
where
k 1 = ( 1 + m ) ( p 1 + p 2 + α p 3 ) , k 2 = ( 1 + m ) ( p 2 + α p 3 ) , k 3 = ( 1 + m ) ( μ p 1 + μ p 2 + p 3 ) , k 4 = ( 1 + m ) ( μ p 2 + p 3 ) .
Equation (6) can be rewritten as
( 1 0 0 1 ) ( θ 1 θ 2 ) + ( k 1 k 2 k 3 k 4 ) ( θ 1 θ 2 ) + ( 1 ξ μ μ ) ( θ 1 θ 2 ) = ( f 10 cos ( ω τ ) 0 ) ,
π 2 < θ i < θ 0 , i = 1 , 2 , θ i + = ( 1 γ ) θ i | θ = θ 0 i = 1 , 2 .
Clearly, Equation (8) represents collision constraint conditions and impact equation.
Consider the initial conditions
( θ 1 ( 0 ) θ 1 ( 0 ) θ 2 ( 0 ) θ 2 ( 0 ) ) = ( θ 0 θ 10 θ 0 θ 20 ) ( θ 0 0 , θ 10 0 , θ 20 0 ) ,
Let ω 1 and ω 2 be the natural frequencies of system (7). Without collision, they can be written as
ω 1 = 1 + μ + ( 1 + μ ) 2 4 μ ( 1 ξ ) 2 , ω 2 = 1 + μ ( 1 + μ ) 2 4 μ ( 1 ξ ) 2 .
Assume that
Φ = ( φ 11 φ 12 φ 21 φ 22 ) = ( 1 1 + ( 1 ω 1 2 ) 2 1 1 + ( 1 ω 2 2 ) 2 1 ω 1 2 1 + ( 1 ω 1 2 ) 2 1 ω 2 2 1 + ( 1 ω 2 2 ) 2 )
is the regular mode matrix of system (7).
Take
X = Φ y ,
where X = ( θ 1 , θ 2 ) T , y = ( y 1 , y 2 ) T . Equation (7) can be decoupled to
I y ¨ + C y ˙ + Λ y = F cos ( ω t ) ,
where
I = ( 1 0 0 1 ) , C = ( 2 ζ ω 1 2 0 0 2 ζ ω 2 2 ) , Λ = ( ω 1 2 0 0 ω 2 2 ) ,
F = ( f 1 f 2 ) T = Φ T P , P = ( f 10 0 ) T .
Using the modal superposition method and considering the transformation (12), the solutions of system (7) are obtained as follows
θ i ( t ) = j = 1 2 φ i j ( e η j t ( a j cos ( ω d j t ) + b j sin ( ω d j t ) ) + A j sin ω t + B j cos ω t )
θ i ( t ) = j = 1 2 φ i j ( e η j t ( b j ω d j η j a j ) cos ( ω d j t ) ( a j ω d j + η j b j ) sin ( ω d j t ) ) + A j ω cos ( ω t ) B j ω sin ( ω t ) )
where ω d j = ω j 2 η j 2 , η j = ζ ω j 2 , a j , b j is the integral constant determined by the initial condition of system (6), i = 1 , 2 .
Let C ( t ) and D are matrices defined as follows
C ( t ) = ( φ 11 e 1 c d 1 φ 12 e 2 c d 2 φ 11 e 1 s d 1 φ 12 e 2 s d 2 φ 11 e 1 g d 1 φ 12 e 2 g d 2 φ 11 e 1 h d 1 φ 12 e 2 h d 2 φ 21 e 1 c d 1 φ 22 e 2 c d 2 φ 21 e 1 s d 1 φ 22 e 2 s d 2 φ 21 e 1 g d 1 φ 22 e 2 g d 2 φ 21 e 1 h d 1 φ 22 e 2 h d 2 ) ,
D = ( φ 11 φ 12 0 0 0 0 φ 11 φ 12 φ 21 φ 22 0 0 0 0 φ 21 φ 22 ) ,
(14) and (15) can be rewritten as
( θ 1 ( t ) θ 1 ( t ) θ 2 ( t ) θ 2 ( t ) ) = C ( t ) ( a 1 a 2 b 1 b 2 ) + D ( A 1 sin ω t + B 1 cos ω t A 2 sin ω t + B 2 cos ω t ω ( B 1 sin ω t A 1 cos ω t ) ω ( B 2 sin ω t A 2 cos ω t ) ) ,
where c d 1 = cos ω d 1 t , c d 2 = cos ω d 2 t , s d 1 = sin ω d 1 t , s d 2 = sin ω d 2 t , e 1 = e η 1 t , e 2 = e η 2 t , g d 1 = η 1 c d 1 + ω d 1 s 1 , g d 2 = η 2 c 2 + ω d 2 s 2 , h d 1 = ω d 1 c 1 η 1 s 1 , h d 2 = ω d 2 c 2 η 2 s 2 ,
A j = ω j 2 ω 2 ( ω j 2 ω 2 ) 2 + 4 ζ 2 ω j 4 ω 2 f j , j = 1 , 2
B j = 2 ζ ω j 2 ω ( ω j 2 ω 2 ) 2 + 4 ζ 2 ω j 4 ω 2 f j , j = 1 , 2
According to Table 1, considering the first kind of collision periodic solutions, the impact periodic solution should satisfy conditions s 1 and x 1 .
From (9), (16) and t = 0 , we have
( θ 10 θ 10 θ 10 θ 20 ) = C ( 0 ) ( a 1 a 2 b 1 b 2 ) + D ( B 1 B 2 ω A 1 ω A 2 ) = ( θ 0 θ 10 θ 0 θ 20 ) ,
where
C ( 0 ) = ( φ 11 φ 12 0 0 η 1 φ 11 η 2 φ 12 φ 11 ω d 1 φ 12 ω d 2 φ 21 φ 22 0 0 η 1 φ 21 η 2 φ 22 φ 21 ω d 2 φ 22 ω d 2 ) , D = ( φ 11 φ 12 0 0 0 0 φ 11 φ 12 φ 21 φ 22 0 0 0 0 φ 21 φ 22 ) .
Moreover, θ 0 0 , θ 10 0 , θ 20 0 should be satisfied.
So, the integral constant a j , b j ( j = 1 , 2 ) should yield
{ φ 11 a 1 + φ 12 a 2 + φ 11 B 1 + φ 12 B 2 0 φ 21 a 1 + φ 22 a 2 + φ 21 B 1 + φ 22 B 2 0 η 1 φ 11 a 1 η 2 φ 12 a 2 + φ 11 ω d 1 b 1 + φ 12 ω d 2 b 2 + ω φ 11 A 1 + ω φ 12 A 2 0 η 1 φ 21 a 1 η 2 φ 22 a 2 + φ 21 ω d 1 b 1 + φ 22 ω d 2 b 2 + ω φ 21 A 1 + ω φ 22 A 2 0 .
Note T 0 and T 0 + are the moments before and after the collision. Since the collision is instantaneous, there are T 0 = T 0 = T 0 + . Collision condition can be expressed as
( θ 1 ( T 0 + ) θ 1 ( T 0 + ) θ 2 ( T 0 + ) θ 2 ( T 0 + ) ) = ( θ 1 ( 0 ) θ 1 ( 0 ) θ 2 ( 0 ) θ 2 ( 0 ) ) = ( θ 0 θ 10 θ 0 θ 20 )
( θ 1 ( T 0 + ) θ 1 ( T 0 + ) θ 2 ( T 0 + ) θ 2 ( T 0 + ) ) = R 1 ( θ 1 ( T 0 ) θ 1 ( T 0 ) θ 2 ( T 0 ) θ 2 ( T 0 ) )
where
R 1 = ( 1 0 0 0 0 ( 1 γ ) 0 0 0 0 1 0 0 0 0 ( 1 γ ) ) , R 1 1 = ( 1 0 0 0 0 1 ( 1 γ ) 0 0 0 0 1 0 0 0 0 1 ( 1 γ ) ) .
So, from (21) and (22), we can obtain
( θ 1 ( T 0 ) θ 1 ( T 0 ) θ 2 ( T 0 ) θ 2 ( T 0 ) ) = R 1 1 ( θ 1 ( T 0 + ) θ 1 ( T 0 + ) θ 2 ( T 0 + ) θ 2 ( T 0 + ) ) = R 1 1 ( θ 0 θ 10 θ 0 θ 20 ) .
Moreover, R 1 is called the corresponding collision recovery matrix for the collision conditions s 1 and x 1 .
According to Equation (16),
( θ 1 ( T 0 ) θ 1 ( T 0 ) θ 2 ( T 0 ) θ 2 ( T 0 ) ) = C ( T 0 ) ( a 1 a 2 b 1 b 2 ) + D ( A 1 sin ( ω T 0 ) + B 1 cos ( ω T 0 ) A 2 sin ( ω T 0 ) + B 2 cos ( ω T 0 ) ω ( B 1 sin ( ω T 0 ) A 1 cos ( ω T 0 ) ) ω ( B 2 sin ( ω T 0 ) A 2 cos ( ω T 0 ) ) )
where
C ( T 0 ) = ( φ 11 e 1 c d 1 φ 12 e 2 c d 2 φ 11 e 1 s d 1 φ 12 e 2 s d 2 φ 11 e 1 g 1 φ 12 e 2 g 2 φ 11 e 1 h 1 φ 12 e 2 h 2 φ 21 e 1 c d 1 φ 22 e 2 c d 2 φ 21 e 1 s d 1 φ 22 e 2 s d 2 φ 21 e 1 g 1 φ 22 e 2 g 2 φ 21 e 1 h 1 φ 22 e 2 h 2 ) , g 1 = η 1 c 1 + ω d 1 s 1 , g 2 = η 2 c 2 + ω d 2 s 2 , h 1 = ω d 1 c 1 η 1 s 1 , h 2 = ω d 2 c 2 η 2 s 2 , c 1 = cos ( 2 n π ω d 1 ω ) ,
s 1 = sin ( 2 n π ω d 1 ω ) , c 2 = cos ( 2 n π ω d 2 ω ) , s 2 = sin ( 2 n π ω d 2 ω ) , D = ( φ 11 φ 12 0 0 0 0 φ 11 φ 12 φ 21 φ 22 0 0 0 0 φ 21 φ 22 )
From (24) and (25),
R 1 1 ( θ 0 θ 10 θ 0 θ 20 ) = C ( T 0 ) ( a 1 a 2 b 1 b 2 ) + D ( B 1 B 2 ω A 1 ω A 2 )
can be known.
From (26),
( θ 0 θ 10 θ 0 θ 20 ) = R 1 C ( T 0 ) ( a 1 a 2 b 1 b 2 ) + R 1 D ( B 1 B 2 ω A 1 ω A 2 )
can be obtained.
From (19) and (27), we can obtain
C ( 0 ) ( a 1 a 2 b 1 b 2 ) + D ( B 1 B 2 ω A 1 ω A 2 ) = R 1 C ( T 0 ) ( a 1 a 2 b 1 b 2 ) + R 1 D ( B 1 B 2 ω A 1 ω A 2 ) ,
i.e.,
( C ( 0 ) R 1 C ( T 0 ) ) ( a 1 a 2 b 1 b 2 ) = ( R 1 E ) D ( B 1 B 2 ω A 1 ω A 2 ) ,
where A 1 , A 2 , B 1 , B 2 are the constants determined by the physical parameters of the system. The existence problem of non-smooth periodic solution of system (6) is transformed into the existence of solutions of (20) and (29). In summary, we can draw a conclusion as follows.
Theorem 1.
If the parameters of the collision system (6) satisfy the condition (20) and (29), system (6) has a one-impact periodic motion.
Let
k 1 = 1 e 1 c 1 , k 2 = 1 e 2 c 2 , m 1 = e 1 s 1 , m 2 = e 2 s 2 , n 1 = η 1 ( 1 γ ) e 1 g 1 , n 2 = η 2 ( 1 γ ) e 2 g 2 , p 1 = ω d 1 ( 1 γ ) e 1 h 1 , p 2 = ω d 2 ( 1 γ ) e 2 h 2 , u 1 = ω ( γ 2 ) ( φ 11 A 1 + φ 12 A 2 ) , u 2 = ω ( γ 2 ) ( φ 21 A 1 + φ 22 A 2 ) ,
Equation (29) can be expressed as follows
( φ 11 k 1 φ 12 k 2 φ 11 m 1 φ 12 m 2 φ 11 n 1 φ 12 n 2 φ 11 p 1 φ 12 p 2 φ 21 k 1 φ 22 k 2 φ 21 m 1 φ 22 m 2 φ 21 n 1 φ 22 n 2 φ 21 p 1 φ 22 p 2 ) ( a 1 a 2 b 1 b 2 ) = ( 0 u 1 0 u 2 ) .
When
k 1 0 , k 2 0 , k 1 p 1 m 1 n 1 0 , k 2 p 2 m 2 n 2 0 ,
by constructing an appropriate invertible transformation matrix
P = ( φ 22 p 1 Ρ 1 φ 22 m 1 Ρ 1 φ 12 p 1 Ρ 1 φ 12 m 1 Ρ 1 φ 21 p 2 Ρ 2 φ 21 m 2 Ρ 2 φ 11 p 2 Ρ 2 φ 11 m 2 Ρ 2 φ 22 n 1 Ρ 1 φ 22 k 1 Ρ 1 φ 12 n 1 Ρ 1 φ 12 k 1 Ρ 1 φ 21 n 1 Ρ 1 φ 21 k 2 Ρ 1 φ 11 n 2 Ρ 1 φ 11 k 2 Ρ 1 ) ,
Equation (31) can be taken as
( 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 ) ( a 1 a 2 b 1 b 2 ) = ( m 1 ( φ 22 u 1 φ 12 u 2 ) P 1 m 2 ( φ 11 u 2 φ 12 u 1 ) P 2 k 1 ( φ 22 u 1 φ 12 u 2 ) P 1 k 1 ( φ 11 u 2 φ 12 u 1 ) P 2 ) ,
where P 1 = ( φ 11 φ 22 φ 12 φ 21 ) ( k 1 p 1 m 1 n 1 ) , P 2 = ( φ 11 φ 22 φ 12 φ 21 ) ( k 2 p 2 m 2 n 2 ) .
If the physical parameters and the collision recovery coefficient of system (6) are given, the unique integral constant a j , b j ( j = 1 , 2 ) can be calculated from (32) as follows
a 1 = m 1 ( φ 22 u 1 φ 12 u 2 ) ( φ 11 φ 22 φ 12 φ 21 ) ( k 1 p 1 m 1 n 1 ) , a 2 = m 2 ( φ 11 u 2 φ 12 u 1 ) ( φ 11 φ 22 φ 12 φ 21 ) ( k 2 p 2 m 2 n 2 )
b 1 = k 1 ( φ 22 u 1 φ 12 u 2 ) ( φ 11 φ 22 φ 12 φ 21 ) ( k 1 p 1 m 1 n 1 ) , b 2 = k 1 ( φ 11 u 2 φ 12 u 1 ) ( φ 11 φ 22 φ 12 φ 21 ) ( k 2 p 2 m 2 n 2 )
Therefore, from (14), (20), (34), and (35), the first kind of collision periodic solution of system (6) is obtained
θ i ( t ) = j = 1 2 φ i j ( e η j t ( a j cos ( ω d j t ) + b j sin ( ω d j t ) ) + A j sin ω t + B j cos ω t ) ( i = 1 , 2 ) ( t = 0 mod 2 n π / ω )
Moreover, when k 1 0 , k 2 0 , k 1 p 1 m 1 n 1 0 , k 2 p 2 m 2 n 2 0 are not held, (29) has an infinite number of solutions or no solutions, which leads to many periodic or no periodic solutions in the original system (6).
For example, if k 1 0 , k 2 0 , k 1 p 1 m 1 n 1 0 , k 2 p 2 m 2 n 2 = 0 , (31) is equivalent to
( φ 11 k 1 φ 12 k 2 φ 11 m 1 φ 12 m 2 0 0 φ 11 ( p 1 m 1 k 1 n 2 ) φ 12 ( p 2 m 2 k 2 n 2 ) 0 1 0 m 2 k 2 0 0 0 p 2 m 2 k 2 n 2 ) ( a 1 a 2 b 1 b 2 ) = ( 0 u 1 0 φ 11 u 2 φ 21 u 1 φ 11 φ 22 φ 12 φ 21 ) .
When φ 11 u 2 φ 21 u 1 0 or u 1 0 is true, Equation (37) has no solution. Hence, system (6) has no collision periodic solution.
Note 1 Conditions (20) and (32) are sufficient and unnecessary conditions for system (6) with the first kind of collision periodic solution.
Note 2 The collision recovery matrices R 2 , R 3 , R 4 of the second type, the third type, and the fourth type of periodic solutions can be shown as
R 2 = ( 1 0 0 0 0 ( 1 γ ) 0 0 0 0 1 0 0 0 0 1 ) , R 3 = ( 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 ( 1 γ ) ) , R 4 = ( 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 ) .
In order to get the second type, the third type and the fourth type of periodic solutions, R 1 can be replaced by R i , i = 2 , 3 , 4 in Equation (28), we have
( C ( 0 ) R i C ( T 0 ) ) ( a 1 a 2 b 1 b 2 ) = ( R i E ) D ( B 1 B 2 ω A 1 ω A 2 ) , i = 2 , 3 , 4
Clearly, for specific matrices, invertible transformation matrices can be constructed. Similar methods can be used to obtain the existence conditions and analytical expressions of periodic solutions.

5. Numerical Simulation

In order to verify the results of theoretical analysis, the suitable physical parameters and coefficient of restitution are selected to solve the periodic solutions according to (20) and (32). From (17), (18), and (20), we know A 2 = 0 , B 2 = 0 and
{ φ 11 a 1 + φ 12 a 2 + φ 11 B 1 0 φ 21 a 1 + φ 22 a 2 + φ 21 B 1 0 η 1 φ 11 a 1 η 2 φ 12 a 2 + φ 11 ω d 1 b 1 + φ 12 ω d 2 b 2 + ω φ 11 A 1 0 η 1 φ 21 a 1 η 2 φ 22 a 2 + φ 21 ω d 1 b 1 + φ 22 ω d 2 b 2 + ω φ 21 A 1 0 .
Let m = 3 , μ = 3.24 , ξ = 0.75 , ω = 1.5 , f 10 = 0.162 , γ = 0.1 , from (20), (30), (31), and (32), the corresponding parameter values are calculated. The phase diagram and Poincare section of the periodic 3 motion of the upper (lower) pendulum of system (6) are shown in Figure 3 and Figure 4. The blue solid line on the right side of the phase diagram is the displacement of the collision point. The red phase point is the projection of the Poincare section on the phase diagram. As can be seen in Figure 3 and Figure 4, the phase diagram and the velocity at the collision point have a momentary jump, indicating that upper (lower) pendulum and the rigid wall collided.
When the parameter value is m = 1 , μ = 3.24 , ξ = 0.5 , ω = 1.5 , f 10 = 0.22 , γ = 0.1 , from (37) and R 2 , a j , b j ( j = 1 , 2 ) can be obtained. Hence, we get the second kind of collision periodic solution. The phase diagram and Poincare section of the periodic 2 motion of the upper pendulum are shown in Figure 5a and Figure 6a. As can be seen in Figure 5a,c, the phase diagram and the velocity at the collision point have a momentary jump, indicating that upper pendulum and the rigid wall collided. But in Figure 6a, there is a small gap between the phase diagram curve and the blue solid line, which indicates that the upper pendulum does not collide with the right wall. As can be seen from Figure 6c, the velocity diagram is continuous, which further confirms the fact that the lower pendulum does not collide.
When the parameter value is m = 1 , μ = 3.2 , ξ = 0.5 , ω = 1.5 , f 10 = 0.1 , γ = 0.1 , from (37) and R 3 , aj, bj(j = 1, 2) can be obtained. Hence, we can detect the third kind of collision periodic solution. The phase diagram and Poincare section of the periodic 2 motion are shown in Figure 7a and Figure 8a. As can be seen in Figure 7a, there is a big gap between the phase diagram curve and the blue solid line, indicating that the upper pendulum does not collide with the right wall and the velocity diagram further confirms the fact in Figure 7c. In Figure 8a,c, the phase diagram and the velocity at the collision point have a momentary jump, showing that lower pendulum and the rigid baffle collided with the right wall.
The simulation results show that the collision periodic solution of the system can be found based on the condition of the existence of the collision analytical periodic solution, which provides a theoretical method for studying the periodic solution of the non-smooth system. The case of collision with left wall can be discussed similarly. Due to the complex type of periodic solution, further analysis will be conducted for the case of bilateral constraint, which could be carried out in a separate paper.

6. Conclusions

A kind of asymmetry system has been constructed by use of the non-smooth unilateral rigid constrained double pendulum, which is subjected to harmonic excitation. Using the mode superposition method and matrix theory, the existence conditions of periodic solutions and the expressions of periodic analytical solutions were deduced for small angle vibration of the system. The numerical simulation showed that this method can predict the existence of periodic solutions of collisions. From the previous calculation, it can be seen that the integral constants and existence conditions of three kinds of periodic solutions of collisions can be calculated by introducing matrix tools. This method provided a matrix computing tool for solving the periodic solutions of high-dimensional non-smooth systems, aiming at multi-point collisions. As long as the appropriate collision recovery matrix was found, the calculation process was similar, which provided a mechanized calculation for simplifying the manipulator into several degrees of freedom problems in the research and design of the manipulator. Thus, a method for the fault research of the fine manipulator was provided to facilitate the engineers to use this method to find the periodic solution of the high-dimensional non-smooth system.

Author Contributions

Conceptualization, X.G. and G.Z.; Methodology, X.G., R.T.; Software, X.G.; Validation, X.G., G.Z., and R.T.; Formal Analysis, X.G.; Investigation, X.G., R.T.; Resources, G.Z., R.T.; Data Curation, R.T.; Writing-Original Draft Preparation, X.G., R.T.; Writing-Review & Editing, X.G.; Visualization, X.G.; Supervision, X.G., R.T.; Project Administration, X.G., R.T.; Funding Acquisition, R.T.

Funding

This research was funded by Natural Science Foundation of China grant number [11872253] and The APC was funded by [11872253].

Acknowledgments

This work was supported by the Natural Science Foundation of China (Nos. 11872253, 11602151), Natural Science Foundation for Outstanding Young Researcher in Hebei Province of China (No. A2017210177), Hundred Excellent Innovative Talents in Hebei Province (SLRC2019037), the Key Project of Education Department in Hebei Province of China (No. ZD2016151) and Basic Research Team Special Support Projects (No. 311008).

Conflicts of Interest

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

  1. Guan, T.M.; Li, N.; Xuan, L.; Wang, M.L. Kinematics Simulation Analysis of Shoulder Joint of New Bionic Robot Arm. J. Dalian Jiaotong Univ. 2018, 39, 65–70. [Google Scholar]
  2. Yang, Q.Z.; Cao, D.F.; Zhao, J.H. Analysis on State of the Art of Upper Limb Rehabilitation Robots. Robot 2013, 35, 631–640. [Google Scholar] [CrossRef]
  3. Hong, Z.R. Design Analysis of Mechanism Structure of Bionic Robotic Arm. J. Lanzhou Petrochem. Polytech. 2017, 17, 1–4. [Google Scholar]
  4. Qian, Z.J.; Zhang, D.G.; Jin, C.Q. A regularized approach for frictional impact dynamics of flexible multi-link manipulator arms considering the dynamic stiffening effect. Multibody Syst. Dyn. 2018, 43, 229–255. [Google Scholar] [CrossRef]
  5. Kumar, R.; Gupta, S.; Ali, S.F. Energy harvesting from chaos in base excited double pendulum. Mech. Syst. Signal Process. 2019, 124, 49–64. [Google Scholar] [CrossRef]
  6. Wojna, M.; Wijata, A.; Wasilewski, G.; Awrejcewicz, J. Numerical and experimental study of a double physical pendulum with magnetic interaction. J. Sound Vib. 2018, 430, 214–230. [Google Scholar] [CrossRef]
  7. Izadgoshasb, I.; Lim, Y.Y.; Tang, L.; Padilla, R.V.; Tang, Z.S.; Sedighi, M.; Sedighi, M. Improving efficiency of piezoelectric based energy harvesting from human motions using double pendulum system. Energy Convers. Manag. 2019, 184, 559–570. [Google Scholar] [CrossRef]
  8. Battelli, F.; Feckan, M. Nonsmooth homoclinic orbits, Melnikov functions and chaos in discontinuous systems. J. Phys. D 2012, 241, 1962–1975. [Google Scholar] [CrossRef]
  9. Tian, R.L.; Zhao, Z.J.; Yang, X.W.; Zhou, Y.F. Subharmonic bifurcation for non-smooth oscillator. Int. J. Bifurc. Chaos 2017, 27, 17501631. [Google Scholar] [CrossRef]
  10. Tian, R.L.; Zhou, Y.F.; Wang, Y.Z.; Feng, W.J.; Yang, X.W. Chaotic threshold for non-smooth system with multiple impulse effect. Nonlinear Dyn. 2016, 85, 1849–1863. [Google Scholar] [CrossRef]
  11. Tian, R.L.; Zhou, Y.F.; Zhang, B.L.; Yang, X.W. Chaotic threshold for a class of impulsive differential system. Nonlinear Dyn. 2016, 83, 2229–2240. [Google Scholar] [CrossRef]
  12. Bi, Q.S.; Chen, Y.S. Bifurcation Analysis of a Double Pendulum With Internal Resonance. Appl. Math. Mech. 2000, 21, 226–234. [Google Scholar]
  13. Han, W.; Hu, H.Y.; Jin, D.P.; Hou, Z.Q. The oblique-impact vibration of a double compound pendulum with the end displacement limit. J. Dyn. Control 2004, 2, 24–30. [Google Scholar]
  14. Han, W.; Hu, H.Y.; Jin, D.P. Oblique impact analysis of two-degree-of-freedom vibration systems. J. Mech. 2003, 35, 723–729. [Google Scholar]
  15. Zhu, X.F.; Cao, X.X. Flutter motion and transition law of two-degree-of-freedom elastic collision system. J. Lanzhou Jiaotong Univ. 2014, 33, 191–195. [Google Scholar]
  16. Hu, H.Y. Numerical analysis of periodic response of high-dimensional non-smooth dynamic systems. Acta Solid Mech. Sin. 1994, 15, 135–145. [Google Scholar]
  17. Li, Q.H.; Lu, Q.S. Analysis to motions of a two-degree-of-freedom vibro-impact system. Acts Mech. Sin. 2001, 33, 776–786. [Google Scholar]
  18. Jin, D.P.; Hu, H.Y. Vibro-impacts and their typical behaviors of mechanical systems. Adv. Mech. 1999, 29, 155–164. [Google Scholar]
  19. Luo, G.W.; Chu, Y.D.; Gou, X.F. Quasi-periodic impact motions and routes to chaos of an impact system with two singer-degree–freedom oscillators. J. Mech. Strength 2005, 27, 445–451. [Google Scholar]
  20. Luo, G.W.; Xie, J.H.; Sun, X.F. Periodic motions and stability of a two-degree-of-freedom vibro-impact system. J. Lanzhou Railw. Inst. 1999, 18, 47–54. [Google Scholar]
  21. Cao, D.Q.; Shu, Z.Z. Periodic motions and robust stability of the multi-degree-of-freedom systems with clearances. Acts Mech. Sin. 1997, 29, 74–83. [Google Scholar]
Figure 1. Collision double pendulum model.
Figure 1. Collision double pendulum model.
Symmetry 11 00886 g001
Figure 2. Phase diagram and time-velocity diagram.
Figure 2. Phase diagram and time-velocity diagram.
Symmetry 11 00886 g002
Figure 3. The first kind of collision periodic solutions: (a) Phase diagram and Poincare section of periodic 3 motion of upper pendulum; (b) Time-displacement diagram of the upper pendulum; (c) Time-velocity diagram of upper pendulum.
Figure 3. The first kind of collision periodic solutions: (a) Phase diagram and Poincare section of periodic 3 motion of upper pendulum; (b) Time-displacement diagram of the upper pendulum; (c) Time-velocity diagram of upper pendulum.
Symmetry 11 00886 g003
Figure 4. The first kind of collision periodic solution: (a) Phase diagram and Poincare section of periodic 3 motion of lower pendulum; (b) The time-displacement diagram of the lower pendulum; (c) The time-velocity diagram of the lower pendulum.
Figure 4. The first kind of collision periodic solution: (a) Phase diagram and Poincare section of periodic 3 motion of lower pendulum; (b) The time-displacement diagram of the lower pendulum; (c) The time-velocity diagram of the lower pendulum.
Symmetry 11 00886 g004
Figure 5. The second kind of collision periodic solutions: (a) Phase diagram and Poincare section of periodic 2 motion of upper pendulum; (b) Time-displacement diagram of the upper pendulum; (c) Time-velocity diagram of upper pendulum.
Figure 5. The second kind of collision periodic solutions: (a) Phase diagram and Poincare section of periodic 2 motion of upper pendulum; (b) Time-displacement diagram of the upper pendulum; (c) Time-velocity diagram of upper pendulum.
Symmetry 11 00886 g005
Figure 6. The second kind of collision periodic solution: (a) phase diagram and Poincare section of periodic 2 motion of lower pendulum; (b) the time-displacement diagram of the lower pendulum; (c) the time-velocity diagram of the lower pendulum.
Figure 6. The second kind of collision periodic solution: (a) phase diagram and Poincare section of periodic 2 motion of lower pendulum; (b) the time-displacement diagram of the lower pendulum; (c) the time-velocity diagram of the lower pendulum.
Symmetry 11 00886 g006
Figure 7. The third kind of collision periodic solutions: (a) Phase diagram and Poincare section of periodic 2 motion of upper pendulum; (b) Time-displacement diagram of the upper pendulum; (c) Time-velocity diagram of upper pendulum.
Figure 7. The third kind of collision periodic solutions: (a) Phase diagram and Poincare section of periodic 2 motion of upper pendulum; (b) Time-displacement diagram of the upper pendulum; (c) Time-velocity diagram of upper pendulum.
Symmetry 11 00886 g007
Figure 8. The third kind of collision periodic solution: (a) Phase diagram and Poincare section of periodic 2 motion of lower pendulum; (b) The time-displacement diagram of the lower pendulum; (c) The time-velocity diagram of the lower pendulum.
Figure 8. The third kind of collision periodic solution: (a) Phase diagram and Poincare section of periodic 2 motion of lower pendulum; (b) The time-displacement diagram of the lower pendulum; (c) The time-velocity diagram of the lower pendulum.
Symmetry 11 00886 g008
Table 1. This classification of periodic solutions.
Table 1. This classification of periodic solutions.
Type of SolutionConstraint Condition
The first kind of periodic solution ( s 1 , x 1 )
The second kind of periodic solution ( s 1 , x 0 )
The third kind of periodic solution ( s 0 , x 1 )
The forth kind of periodic solution ( s 0 , x 0 )

Share and Cite

MDPI and ACS Style

Guo, X.; Zhang, G.; Tian, R. Periodic Solution of a Non-Smooth Double Pendulum with Unilateral Rigid Constrain. Symmetry 2019, 11, 886. https://doi.org/10.3390/sym11070886

AMA Style

Guo X, Zhang G, Tian R. Periodic Solution of a Non-Smooth Double Pendulum with Unilateral Rigid Constrain. Symmetry. 2019; 11(7):886. https://doi.org/10.3390/sym11070886

Chicago/Turabian Style

Guo, Xiuying, Gang Zhang, and Ruilan Tian. 2019. "Periodic Solution of a Non-Smooth Double Pendulum with Unilateral Rigid Constrain" Symmetry 11, no. 7: 886. https://doi.org/10.3390/sym11070886

APA Style

Guo, X., Zhang, G., & Tian, R. (2019). Periodic Solution of a Non-Smooth Double Pendulum with Unilateral Rigid Constrain. Symmetry, 11(7), 886. https://doi.org/10.3390/sym11070886

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