Next Article in Journal
State-Space Approach to the Time-Fractional Maxwell’s Equations under Caputo Fractional Derivative of an Electromagnetic Half-Space under Four Different Thermoelastic Theorems
Previous Article in Journal
Asymptotic Synchronization for Caputo Fractional-Order Time-Delayed Cellar Neural Networks with Multiple Fuzzy Operators and Partial Uncertainties via Mixed Impulsive Feedback Control
Previous Article in Special Issue
A Low-Computational Burden Closed-Form Approximated Expression for MSE Applicable for PTP with gfGn Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Analysis of Viscoelastic Friction System Characteristics of Robotic Arm Brake Based on Fractional Differential Theory

1
College of Electronic Science and Engineering, Jilin University, Changchun 130012, China
2
State Key Laboratory of Integrated Optoelectronics, JLU Region, Changchun 130012, China
3
State Grid Jibei Electric Power Co., Ltd., Material Branch, Beijing 100075, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(10), 565; https://doi.org/10.3390/fractalfract8100565
Submission received: 5 August 2024 / Revised: 22 September 2024 / Accepted: 26 September 2024 / Published: 28 September 2024
(This article belongs to the Special Issue Fractional Processes and Systems in Computer Science and Engineering)

Abstract

:
In engineering practice, the nonlinear vibration effect can easily lead to chaos in the system, which will not only reduce the performance of the system but also lead to premature fatigue of components, control failure, and increased safety risks. In view of the core position of the robotic arm in modern industry, this study relies on the robotic arm brake system to explore the theoretical basis of integrated viscoelastic materials as a vibration isolation layer. By analyzing the dynamic characteristics of the friction braking system with fractional differential terms, it aims to provide a new perspective for understanding and controlling the chaotic phenomena of a class of nonlinear friction systems. Firstly, we construct a model of a friction system and analyze its dynamic characteristics in detail. The self-excited vibration of the system under disturbance is studied. The relationship between amplitude and frequency is calculated by a nonlinear approximate analytical algorithm, and the accuracy of this relationship is verified by a numerical algorithm. Then, we compare the differences between non-fractional systems and fractional systems. It is found that with the increase in the fractional order term, the vibration amplitude of the system decreases significantly, which helps to reduce the nonlinear characteristics generated by the friction system and narrow the range of unstable solutions. Secondly, we also study the influence of parameter coefficients on the amplitude–frequency characteristics and analyze the local static bifurcation characteristics through singularity theory. Finally, we study the dynamic bifurcation behavior under different parameter perturbations and find that the change in system parameters will lead to the alternation of periodic motion and chaotic motion.

1. Introduction

In today’s society, collaborative robots [1] have become the fastest-growing and most widely marketed type in the global robot industry. The cooperation between collaborative robots and humans in industrial production can give full play to their respective advantages and improve work efficiency, thus forming a healthy and efficient production mode. Thus far, collaborative robots have been applied to various industries, such as manufacturing [2], automobiles [3], electronics [4], biomedicine [5], logistics [6], agriculture [7], aerospace [8], and many other industries. In the process of human–machine cooperation, we must first ensure the safety of both humans and machines. As an important component of robots, brakes can provide good protection for robots, especially when people manipulate machines, which can protect both the robot and the controller. Therefore, it is necessary to study brakes.
A brake is a device that slows down, stops, or maintains the stopped state of moving parts. It is widely used in high-speed rail [9], robotic arms [10], new energy electric vehicles [11], and other fields, attracting many scholars to study it. Yuan [12] used finite element software to compare the stability of brake systems with and without noise reduction pads, studying how the use of noise reduction pads improves brake friction squeal. Additionally, the relationship between the muffler’s structure and brake friction squeal was investigated. Liu et al. [13] conducted a study on the force amplification coefficient that affects the force transmission performance of the dry disk brake cartridge pressurization mechanism during the process and return stroke, aiming at the mobility and braking efficiency of tracked vehicles. Suo et al. [14] obtained the thermal elastoplastic constitutive relation of the brake disk through experiments and simulated the stress–strain response relationship of the lower brake disk using numerical calculations based on the sequential coupling method. In many countries, brakes are used to ensure the safety of joints and equipment, such as the robotic arm of the National Aeronautics and Space Administration [15], the European ERA [16], and the space robotic arm designed by Chinese universities [17].
Fractional calculus is an important branch of mathematics. Fractional calculus not only occupies an inevitable position in the development of mathematics, but its unique memory characteristics and precise description of real-world models, especially its simplicity in describing nonlinear models, have attracted widespread attention in both theory and practice. In response to the needs of applied disciplines, the theory and computational techniques of fractional calculus continue to advance and have become a core tool for solving complex scientific and engineering problems [18,19,20]. In some systems that require high accuracy, researchers have increasingly high requirements for describing systems. Applying fractional calculus to differential equations simplifies the resulting system differential equations, making the calculated results more accurate. Fractional calculus models are applied in various engineering fields [21,22,23,24,25,26,27,28]. To address the issue of low harmonic current at the DC bus in DC microgrid systems, Lin et al. [29] proposed a fractional capacitance-based method to force low harmonic currents. This method can theoretically achieve complete suppression of any or combination of low harmonic currents. Zhu et al. [30] used the backstepping method to design a boundary controller for fractional order reaction diffusion systems with spatially dependent coupling coefficients and proved the fitness of the observation gain and control gain kernel matrix equations. For error systems and closed-loop systems with output feedback, the Mittag–Leffler stability of the system was analyzed using the fractional order Lyapunov method. Moreover, using Wirtinger’s inequality, the stability conditions of the coupled system are improved, and a numerical solution method for the kernel function partial differential equation is given. Numerical simulation verifies the theoretical results. Zhu et al. [31] proposed fractional order power flow and voltage analysis for systems with fractional order inductors and capacitors. The results show that power system analysis based on fractional order component models can not only reduce reactive power losses in branches but also help improve voltage stability in power systems.
Viscoelastic materials [32] are widely used in various industries [33,34,35]. He et al. [36] used laminated piezoelectric actuators as active isolation elements and designed passive isolation elements based on viscoelastic materials. Li et al. [37] proposed a novel hybrid isolator for simulating rigid satellites. Tang et al. [38] studied rubber isolators for auxiliary power units. A nonlinear friction derivative dynamic model of rubber isolators has been established. Chang et al. [39] proposed a viscoelastic constitutive model for metal rubber that includes fractional order differentiation. Experimental results have shown that the proposed nonlinear dynamic system model of metal rubber with fractional order subterms has a continuous mathematical expression and can accurately reflect the complete dynamic performance of the nonlinear system of metal rubber.
In various projects, there are few linear problems but more nonlinear problems. Therefore, there have been many studies on nonlinear problems. Hamed [40] studied the approximate nonlinear dynamic behavior of AFM systems using perturbation methods and frequency response analysis. Kandil et al. [41] studied the oscillation behavior of the bearing system and derived a nonlinear dynamic equation for the control system. In order to eliminate the vibration of nonlinear dynamic beams, Ms et al. [42] added active control to their research and analyzed its stability, obtaining analytical and numerical solutions. Wang [43,44] established a nonlinear dynamic equation for the relative rotation system of two end-face rotating shafts, qualitatively analyzed the dynamic equation of equal torque, and studied the stability and other properties of the equation. Finally, an approximate analytical algorithm was developed to calculate the approximate solution of the equation under specific conditions. In addition, the uniqueness and exact periodic solutions of a class of relative rotational nonlinear dynamic systems were also studied.
In summary, the research on robots and brakes mainly considers the design of robot brakes and the dynamic performance of disk brakes. However, the dynamic performance analysis of nonlinear systems in robot electromagnetic brakes that incorporate fractional differentiation has not yet been fully investigated. In this paper, viscoelastic materials are added to the electromagnetic brake of a robot, which is used as a vibration isolation layer. The stress of the newly added material is represented by a fractional order term, thereby establishing a dynamic model of the electromagnetic brake of a robotic arm with fractional differentiation. We calculated the relationship between the frequency and amplitude under perturbation external conditions. By altering the external conditions, this study examined the characteristics of the friction system, analyzed the conditions for bifurcation behavior, and analyzed the local static and global steady-state characteristics.

2. Calculating an Approximate Solution

The new brake of the robotic arm is composed of a friction plate, brake caliper, electromagnet, spring, brake disk, shaft, and vibration isolation layer. Figure 1 is the functional structure diagram of the new electromagnetic brake of the robotic arm with vibration isolation layer. Figure 2 shows two states of the electromagnetic brake: (a) When the electrical current is switched off, the friction plate contacts the brake disk under the action of spring force, and the torque generated by the system is quickly stopped by the drive shaft. (b) When the current passes through the magnetic coil of the electromagnetic brake, the friction plate is pulled apart under the action of the magnetic force, releasing the brake disk. At this time, the transmission shaft runs normally.
By analyzing the braking process of the electromagnetic brake, obtain the simplified model (Figure 3):
Kinetic equation:
m x ¨ ( t ) + k 1 x ( t ) + c 1 x ˙ ( t ) + k 2 x ( t ) 3 + c 2 x ˙ ( t ) 3 + K D p 1 [ x ( t ) ] F f = F cos ( ω t ) .
In Equation (1), k1 and k2 are linear and nonlinear stiffness coefficients; c1 and c2 are linear and nonlinear damping coefficients, respectively; F cos ( ω t ) is defined as an external disturbance function; D p 1 x t is the fractional derivative; p 1 ( 0 p 1 1 ) is the order of the fractional derivative; and K(K > 0) is the coefficient of the fractional derivative. F f is the friction between the mass and the conveyor belt, which belongs to dry friction.
There are many forms of definitions of fractional differential. In this paper, the expression of fractional order terms is defined using the Caputo [45] definition. Let the domain of the function x(t) be (a, b), and assume that x(t) is continuously differentiable up to order n on (a, b), where n − 1 < p < n and p > 0, then the Caputo fractional derivative is defined as follows:
D t p t 0 C x t = 1 Γ n p t 0 t t u n p 1 x ( n ) u d u , n 1 < p < n ( n N + ) x ( n ) t , p = n ( n N ) .
In Equation (2), Γ t is a Gamma function satisfying Γ z + 1 = z Γ z when z > 0.
In order to study the vibration generated by the brake during braking and its impact on frictional vibration, a static friction model is applied to the motion process of the system as a friction force. For high-performance motion control, in addition to the robot’s own coupling and nonlinear time-varying effects, the effect of nonlinear friction, such as the Stribeck effect, is particularly pronounced at low speeds. The dry friction model with the Stribeck effect exhibits a negative slope characteristic at low relative velocities. When the relative velocity reaches a certain threshold, the dry friction coefficient remains unchanged. In this paper, a dry friction model with the Stribeck effect is selected.
F f = F n μ ( x ˙ ( t ) v c ) ,
where F n denotes external force and μ ( x ˙ ( t ) v c ) expresses the negative slope characteristics of the friction coefficient using a polynomial function model:
μ ( x ˙ ( t ) v c ) = μ s sign ( x ˙ ( t ) v c ) + k 3 ( x ˙ ( t ) v c ) + k 4 ( x ˙ ( t ) v c ) 3 .
where μ represents the coefficient of friction; k3 and k4 are constants; μ s is the maximum static friction coefficient; and vc is the speed of the conveyor belt.
In actual robotic arm operation, the brake usually works in a state close to equilibrium, and the nonlinear effect of its dynamic response is relatively weak. Introducing a small parameter ε can simulate this weak nonlinear effect, making the model closer to actual operating conditions. Therefore, Equation (1) is simplified into the following dimensionless form:
x ¨ ( t ) + ω 0 2 x ( t ) + ε β 1 x ˙ ( t ) + ε α 2 x ( t ) 3 + ε β 2 x ˙ ( t ) 3 + ε γ D p 1 [ x ( t ) ] ε f = ε f 0 cos ( ω t ) ,
ε f [ x ˙ ( t ) ] = ε [ υ sign ( x ˙ ( t ) v c ) + α n 3 x ˙ ( t ) + α n 4 ( x ˙ ( t ) 3 3 x ˙ ( t ) 2 v c + 3 x ˙ ( t ) v c 2 ) ] ,
where ω 0 = k 1 / m , ε α 1 = k 1 / m , ε β 1 = c 1 / m , ε α 2 = k 2 / m , ε β 2 = c 2 / m , ε γ = K / m , ε f 0 = F / m , ε f = F f / m , ε υ = F n μ s / m , ε α n 3 = F n k 3 / m , ε α n 4 = F n k 4 / m , and ε ς = F n μ m / m .
This article uses the average method to obtain the analytical solution of the system. The averaging method is a technique for handling weakly nonlinear system dynamics, which simplifies problem solving by separating the fast oscillations and slow varying parts of the system response and eliminating the fast oscillations. This method transforms complex nonlinear problems into easier-to-handle linear or weakly nonlinear problems, preserving key dynamic characteristics of the system and helping to predict long-term behavior. In this article, the average method can effectively analyze the dynamic behavior of the robotic arm brake and provide theoretical support for design optimization. Let ω 2 = ω 0 2 1 + ε σ 1 and σ 1 be tuning factors, then Equation (5) can be written as follows:
x ¨ t + ω 2 x t = ε ω 0 2 σ 1 x t + ε f 0 cos ( ω t ) ε β 1 x ˙ ( t ) + α 2 x ( t ) 3 + β 2 x ˙ ( t ) 3 + γ D p 1 [ x ( t ) ] f ( x ˙ ( t ) ) .
When ε = 0 , the solution of the derived system and its derivative are as follows:
x t = a cos ω t + θ ,
x ˙ t = a ω sin ω t + θ .
The changes of a and θ are considered here. Let ϕ = ω t + θ a , by differentiating Equation (8) to time t and eliminating Equation (9), we obtain the following:
a ˙ cos ϕ a θ ˙ sin ϕ = 0 ,
a ˙ sin ϕ a θ ˙ cos ϕ = ε ω ω 0 2 σ 1 x t + f 0 cos ( ω t ) β 1 x ˙ ( t ) α 2 x ( t ) 3 β 2 x ˙ ( t ) 3 γ D p 1 [ x ( t ) ] + f ( x ˙ ( t ) ) .
Further obtaining:
a ˙ = 1 ω P 1 a , θ + P 2 a , θ + P 3 a , θ sin ϕ ,
θ ˙ = 1 a ω P 1 a , θ + P 2 a , θ + P 3 a , θ cos ϕ ,
where:
P 1 a , θ = ε γ D p 1 a cos ϕ P 2 a , θ = ε ω 0 2 σ 1 a cos ϕ + f 0 cos ϕ + θ + β 1 a ω sin ϕ α 2 ( a cos ϕ ) 3 + β 2 ( a ω sin ϕ ) 3 P 3 a , θ = ε f ( a ω sin ϕ ) .
Integrate Equations (12) and (13) as follows:
a ˙ = 1 ω T 0 T P 1 a , θ sin ϕ d ϕ 1 ω T 0 T P 2 a , θ sin ϕ d ϕ 1 ω T 0 T P 3 a , θ sin ϕ d ϕ = a ˙ 1 + a ˙ 2 + a ˙ 3 ,
θ ˙ = 1 a ω T 0 T P 1 a , θ cos ϕ d ϕ 1 a ω T 0 T P 2 a , θ cos ϕ d ϕ 1 a ω T 0 T P 3 a , θ cos ϕ d ϕ = θ ˙ 1 + θ ˙ 2 + θ ˙ 3 .
Before calculating the fractional integral, it is important to introduce two basic formulas:
B 1 = lim T 0 T sin ω t t p d t t = s 1 1 p , d t = 1 1 p s p 1 p d s B 1 = 1 1 p lim T 0 T 1 p sin ω s 1 1 p d s ,
B 2 = lim T 0 T cos ω t t p d t t = s 1 1 p , d t = 1 1 p s p 1 p d s B 2 = 1 1 p lim T 0 T 1 p cos ω s 1 1 p d s .
Using the residue theorem and the contour integral:
B 1 = ω p 1 Γ 2 p 1 p cos p π 2 = ω p 1 Γ 1 p cos p π 2 ,
B 2 = ω p 1 Γ 2 p 1 p sin p π 2 = ω p 1 Γ 1 p sin p π 2 .
Calculate the integral of the first part of Equation (15):
a ˙ 1 = lim T 1 T ω 0 T P 1 a , θ sin ϕ d ϕ = lim T ε γ T ω 0 T D p 1 a cos ϕ sin ϕ d ϕ = ε a γ Γ 1 p 1 lim T 1 T 0 T 0 t cos ω s s p 1 d s sin ω t + θ sin ω t + θ d t + ε a γ Γ 1 p 1 lim T 1 T 0 T 0 t sin ω s s p 1 d s cos ω t + θ sin ω t + θ d t = A 1 + A 2 .
A 1 = ε a γ ω p 1 1 2 sin p 1 π 2 .
Similarly, it can be determined that when T , A 2 in Equation (21) approaches 0, so:
a ˙ 1 = A 1 = ε a γ ω p 1 1 2 sin p 1 π 2 .
Using the same method, we can obtain the integral of the first part of Equation (16).
θ ˙ 1 = ε γ ω p 1 1 2 cos p 1 π 2 .
Next:
a ˙ 2 = 1 ω T 0 T P 2 a , θ sin ϕ d ϕ = ε 2 ω f 0 sin θ + β 1 a ω + 3 β 2 a 3 ω 3 4 ,
θ ˙ 2 = 1 a ω T 0 T P 2 a , θ cos ϕ d ϕ = ε 2 ω a ω 0 2 σ 1 a + f 0 cos θ 3 α 2 a 3 4 .
Finally:
a ˙ 3 = 1 ω T 0 T P 3 a , θ sin ϕ d ϕ = ε 2 α n 3 a + α n 4 ( 3 4 a 3 ω 2 + 3 a v c 2 ) ,
θ ˙ 3 = 1 T a ω 0 T P 3 a , θ cos ϕ d ϕ = 0 .
Therefore, it is obtained that:
a ˙ = ε a γ ω p 1 1 2 sin p 1 π 2 ε 2 ω f 0 sin θ + β 1 a ω + 3 β 2 a 3 ω 3 4 + ε 2 α n 3 a + α n 4 ( 3 4 a 3 ω 2 + 3 a v c 2 ) ,
θ ˙ = ε γ ω p 1 1 2 cos p 1 π 2 ε 2 ω a ω 0 2 σ 1 a + f 0 cos θ 3 α 2 a 3 4 .
The singular point a s , θ s in the moving phase plane is the solution of Equations (29) and (30), and eliminating θ s leads to the amplitude frequency curve equation of the nonlinear system:
a γ ω p 1 sin p 1 π 2 β 1 a ω 3 β 2 a 3 ω 3 4 + ω α n 3 a + α n 4 ( 3 4 a 3 ω 2 + 3 a v c 2 ) 2 + a γ ω p 1 cos p 1 π 2 ω 0 2 σ 1 a 3 α 2 a 3 4 2 f 0 2 = W a s , ω .

3. Stability Analysis

Analyze the stability of the system. Let a ˙ = 0 and θ ˙ = 0 to obtain:
0 = ε a γ ω p 1 1 2 sin p 1 π 2 ε 2 ω f 0 sin θ + β 1 a ω + 3 β 2 a 3 ω 3 4 + ε 2 α n 3 a + α n 4 ( 3 4 a 3 ω 2 + 3 a v c 2 ) ,
0 = ε γ ω p 1 1 2 cos p 1 π 2 ε 2 ω a ω 0 2 σ 1 a + f 0 cos θ 3 α 2 a 3 4 ,
Let a = a 0 + Δ a , θ = θ 0 + Δ θ , and substitute Equations (32) and (33):
d Δ a d t = ε γ ω p 1 1 2 sin p 1 π 2 + ε D e Δ a + ε 2 ω f 0 cos θ 0 Δ θ ,
d Δ θ d t = ε γ ω p 1 1 2 cos p 1 π 2 ε 2 ω ω 0 2 σ 1 9 α 2 a 0 2 4 Δ a + ε 2 ω f 0 sin θ 0 Δ θ ,
where D e = 1 2 ω β 1 ω + 3 β 2 3 a 0 2 ω 3 4 + 1 2 α n 3 + α n 4 ( 3 4 3 a 0 2 ω 2 + 3 v c 2 ) .
The characteristic determinant can be obtained from Equations (34) and (35) as follows:
γ ω p 1 1 2 sin p 1 π 2 + D e λ γ ω p 1 1 2 cos p 1 π 2 1 2 ω ω 0 2 σ 1 9 α 2 a 0 2 4 a D e + a γ ω p 1 1 2 sin p 1 π 2 γ ω p 1 1 2 sin p 1 π 2 + D e λ = λ 2 + N 1 λ + N 2 = 0 ,
where λ is the root of the characteristic equation,
N 1 = γ ω p 1 1 2 sin ( p 1 π 2 ) + D e 2 a D e + a γ ω p 1 1 2 sin ( p 1 π 2 ) γ ω p 1 1 2 cos ( p 1 π 2 ) 1 2 ω ω 0 2 σ 1 9 α 2 a 0 2 4 , N 2 = 2 γ ω p 1 1 2 sin ( p 1 π 2 ) + D e .
According to the Lyapunov stability theory and Routh criterion, if N1 > 0 and N2 > 0, the trajectory of the singularity is asymptotically stable. If N2 < 0, the trajectory corresponding to the singularity is unstable. Therefore, N2 = 0 is the critical condition for determining whether the trajectory is stable.

4. Curve of Amplitude Versus Frequency and Verification of Numerical Solution

According to Equation (31), take a set of operating conditions as follows: F0 = 1, Fn = 0.2, m = 1, k1 = 1, k2 = 1, k3 = 1, k4 = 0.1, c1 = 0.1, c2 = 0.01, K = 0.8, p = 0.4, vc = 0.1, and us = 0.1.
Analyze the relationship between the amplitude and frequency of a system with fractional subdivisions (FOS, K = 0.8, p = 0.4) and a system without fractional subdivisions (NFOS, K = 0, p = 0).
Obtained from Figure 4, the amplitude–frequency relationship curve lags behind, and there are unstable regions in both FOS and NFOS. In FOS, it can be seen that due to the existence of fractional subdivisions, the maximum value of the curve is reduced, and the nonlinearity is also reduced. The figure shows that with A as the demarcation point, when the abscissa ω > A , the amplitude of FOS is greater than that of NFOS, and when the abscissa ω > A , the amplitude of FOS is less than that of NFOS.
Using the Runge–Kutta numerical method to verify the accuracy of the approximate analytical solution for the main resonance of the system calculated by the averaging method. Calculate the numerical iteration format of the system:
Z 1 ( t l ) = Z 2 ( t l 1 ) h j = 1 l C j 1 Z 1 ( t l j ) ,
Z 2 ( t l ) = [ k 1 Z 1 ( t l 1 ) c 1 Z 2 ( t l 1 ) k 2 Z 1 ( t l 1 ) 3 c 1 Z 2 ( t l 1 ) 3 + F n ( μ s s i g n ( Z 1 ( t l 1 ) v c ) + k 3 ( Z 1 ( t l 1 ) v c ) + k 4 ( Z 1 ( t l 1 ) v c ) 3 ) + F cos ω h i K 1 Z 3 ( t l 1 ) ] h j = 1 l C j 1 Z 2 ( t l j ) ,
Z 3 ( t l ) = Z 2 ( t l 1 ) h 1 p 1 j = 1 l C j 1 p 1 Z 3 ( t l j ) ,
where Z 1 is displacement; Z 2 is velocity; Z 3 is fractional order sub-item; and h is step size.
The numerical solution of the system was obtained by forward and backward frequency scanning, as shown in Figure 5. Comparing the approximate analytical solution calculated by the averaging method with the solution calculated by the numerical method, it was found that the solutions obtained by the two methods were relatively close, and the correctness of the approximate analytical solution could be obtained.

5. Analysis of Amplitude Frequency Response Characteristics

Analyze the amplitude–frequency response characteristics under different parameter values. The obtained results are represented in Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14 and Figure 15.
From Figure 6, when K increases, the maximum value of the FOS curve decreases, the resonance interval gradually shifts to a small frequency range, and the nonlinearity weakens. It shows that the increase in K not only increases the damping of the system but also makes the system skeleton curve move to the left. When K increases, the maximum value of the relation curve between NFOS amplitude and frequency increases, and the resonance interval also changes to a small frequency range, but the nonlinearity tends to increase. When K takes the same value, the resonance interval of the FOS curve is located at the bottom left of the resonance interval of the NFOS amplitude–frequency curve. In conclusion, the vibration amplitude of FOS is smaller than that of NFOS in the resonance range.
When p is taken as a larger value, there is no unstable solution interval in the amplitude and frequency relationship curve, thus reducing the nonlinear characteristics. As the value of K increases, the peak value gradually decreases. When p is small, when K gradually changes to a larger value, the maximum amplitude gradually decreases. The resonance interval gradually moves to the low-frequency interval.
From Figure 8, we can see that when the linear stiffness gradually increases, the resonance interval gradually moves towards a higher frequency range, and the peak value gradually decreases. From Figure 9, when the value of nonlinear stiffness k2 changes to a larger value, the bending degree of the curve of the system increases, resulting in the resonance interval biased towards a higher frequency range. The results indicate that when the k2 value changes significantly, the nonlinear characteristics of the system become stronger. As shown in Figure 10 and Figure 11, when the values of k3 and k4 change significantly, the relationship curve between amplitude and frequency decreases, and the bending degree decreases gradually. It shows that the increase of k4 weakens the nonlinearity of the system, and the influence on the curve is larger than that of Figure 10 and Figure 11.
As shown in Figure 12, when the value of linear stiffness c1 changes to a larger value. The nonlinear characteristics of the system become stronger. Resonance interval tends to shift to a higher frequency range. Analyzing Figure 13, when Fn = 0, there is an unstable solution region. As Fn increases to 0.6, the system amplitude gradually decreases. When Fn is equal to 0.6, when the value of ω changes to a larger value, the curve of amplitude and frequency decreases. From Figure 14 and Figure 15, as ω increases, the system amplitude increases with the increase of c2. When c2 is greater than 0.025, an independent loop appears, and the nonlinearity is significantly enhanced.

6. Study on the Local Static Bifurcation Response Characteristics of the System

Construct the unfolding function by Equation (31):
G = 9 16 α 2 2 + 3 4 a n 4 ω 3 3 4 b 2 ω 3 2 a 6   + 3 α 2 2 σ 1 + γ ω p 1 c o s ( π p 1 2 ) 3 2 a n 4 ω 3 3 2 β 2 ω 3 β 1 ω ω 3 a n 4 v c 2 + a n 3 + γ ω p 1 s i n ( π p 1 2 ) a 4   + β 1 ω ω 3 a n 4 v c 2 + a n 3 + γ ω p 1 s i n ( π p 1 2 ) 2 + γ ω p 1 c o s ( π p 1 2 ) a 2 f 0 2 .
Bifurcation point set:
B = χ R k | ( a , ξ ) , G = G a = G ω = 0 .
Lag point set:
H = χ R k | ( a , ξ ) , G = G a a = G a = 0 .
Hyperbolic limit point set:
D = χ R k | ( a , ξ ) ( j = 1 , 2 , a 1 a 2 ) , G = G a = 0 .
After calculation:
G a = 2 3 4 a 3 β 2 ω 3 ω a a n 3 + a n 4 3 4 a 3 ω 2 + 3 a v c 2 + a β 1 ω + γ ω p 1 s i n ( π p 1 2 ) β 1 ω ω a n 3 + a n 4 9 4 a 2 ω 2 + 3 v c 2 + 9 4 a 2 β 2 ω 3 + γ ω p 1 s i n ( π p 1 2 ) + 2 a σ 1 3 4 a 3 α 2 + a γ ω p 1 c o s ( π p 1 2 ) σ 1 9 4 a 2 α 2 + γ ω p 1 c o s ( π p 1 2 ) ,
G ω = 2 a p 1 γ ω p 1 1 c o s ( π p 1 2 ) a σ 1 3 4 a 3 α 2 + a γ ω p 1 c o s ( π p 1 2 ) 2 3 4 a 3 β 2 ω 3 ω a a n 3 + a n 4 3 4 a 3 ω 2 + 3 a v c 2 + a β 1 ω + a γ ω p 1 sin ( π p 1 2 ) a a n 3 a β 1 + a n 4 3 4 a 3 ω 2 + 3 a v c 2 + 3 2 a 3 a n 4 ω 2 9 4 a 3 β 2 ω 2 a p 1 γ ω p 1 1 sin ( π p 1 2 ) ,
G a a = 2 σ 1 9 4 a 2 α 2 + γ ω p 1 cos ( π p 1 2 ) 2 2 9 2 9 a a n 4 ω 3 9 2 a β 2 ω 3 3 4 a 3 β 2 ω 3 ω a a n 3 + a n 4 3 4 a 3 ω 2 + 3 a v c 2 + a β 1 ω + a γ ω p 1 sin ( π p 1 2 ) + 2 β 1 ω ω a n 3 + a n 4 9 4 a 2 ω 2 + 3 v c 2 + 9 4 a 2 β 2 ω 3 + γ ω p 1 sin ( π p 1 2 ) 2 9 a α 2 a σ 1 3 4 a 3 β 2 + a γ ω p 1 cos ( π p 1 2 ) .
The system parameters are as follows: F0 = 1, Fn = 0.2, m = 1, k1 = 1, k2 = 1, k3 = 1, k4 = 0.1, c1 = 0.1, c2 = 0.01, K = 4, p = 0.1, vc = 0.1, and us = 0.1. The bifurcation point set and the lag point set are calculated as shown in Figure 16 and Figure 17.
As can be seen from Figure 16, the intersection point of the unfolding surface is the lag point set of the system when the current parameter value is: H = (1.706, 4.066, 3.548). According to the analysis of Figure 17, when the former parameter values are taken, the intersection point of the expanded surface is the bifurcation point set of the system: B = (2.387, 2.201, 0.3124).
By sorting Equation (37), let:
ξ = 9 16 α 2 2 + 3 4 a n 4 ω 3 3 4 b 2 ω 3 2 ξ 1 = 3 α 2 2 σ 1 + γ ω p 1 c o s ( π p 1 2 ) 3 2 a n 4 ω 3 3 2 β 2 ω 3 β 1 ω ω 3 a n 4 v c 2 + a n 3 + γ ω p 1 s i n ( π p 1 2 ) ξ 2 = β 1 ω ω 3 a n 4 v c 2 + a n 3 + γ ω p 1 s i n ( π p 1 2 ) 2 + γ ω p 1 c o s ( π p 1 2 ) ζ = f 0 2 .
Further obtaining:
0 = a 6 + ξ 1 ξ a 4 + ξ 2 ξ a 2 ζ ξ .
Equation (47) is a general expansion of equation 0 = a 6 + ζ ξ ; the codimension is 2, and the system transition set is then calculated. The migration set can be represented as = B H D . Figure 8 is the system transition set and bifurcation diagram.
The parameters F0 = 1, Fn = 0.2, m = 1, k1 = 1, k2 = 1, k3 = 1, k4 = 0.1, c1 = 0.1, c2 = 0.01, K = 4, p = 0.1, vc = 0.1, and us = 0.1 are selected. The analysis of Figure 18 shows that the corresponding bifurcation topology can be obtained by taking the values of ξ 1 / ξ and ξ 2 / ξ in each region. By analyzing Figure 18b, we conclude that when the unfolding parameter passes through the transition set H from region (1) into region (2), the jump phenomenon of the system can be seen from the bifurcation topology diagram. Parts (2)–(4) in Figure 18b show different bifurcation states.

7. Bifurcation Response Characteristics

Take a set of operating parameters of the system and conduct bifurcation research with the following parameter values: F0 = 1, m = 1, k1 = 1, k2 = 1, k3 = 1, k4 = 0.1, c1 = 0.1, p = 0.1, vc = 0.1, and us = 0.1.
From Figure 19, it can be seen that when Fn changes, the motion is different. When the Fn value increases, the system state initially undergoes generalized three-period motion and then becomes a large-scale chaotic state. It was not until Fn became 4.2 that it transitioned to a small degree of chaotic motion. When in a chaotic state, the corresponding LLE graph value is also greater than 0. The Fn values of different motion states are taken below; the x ˙ x diagrams and Poincaré diagrams are drawn to verify different motion states.
The data in Figure 19 is summarized in Table 1.
From Figure 20, it can be seen that when K changes, it shows a different motion state from Figure 19. When the value of Fn changes to a larger value, the state is a single cycle → chaotic state → two period motion → chaotic state → three period motion → chaos → single period motion. Compared with K = 0.5, the system state is more abundant at K = 1.5. In a single cycle state, the LLE value is less than 0, while the value of LLE in other motion states is greater than 0. The Fn values of different motion states are 2, 2.55, 2.6, 2.71, 2.88, 3.5, 3.8, and 4.35, respectively. The x ˙ x diagrams and Poincaré sections are drawn to verify different motion states.
The data in Figure 20 is summarized in Table 2.
It can be seen from Figure 21 that when K is equal to 2, when Fn increases, the system changes in the order of single period → chaotic state → two-period motion. In a single cycle state, the value of LLE is less than 0, while the value of LLE in other motion states is greater than 0. Next, Fn values of different motion states are taken as 2, 3.5, and 4, respectively, and the x ˙ x picture and Poincaré sections are calculated to verify different motion states of the system.
The data in Figure 21 is summarized in Table 3.
Finally, to visually demonstrate the impact of fractional order on bifurcation, we selected K as the bifurcation parameter in our study.
Analyzing Figure 22, the system exhibits a variety of rich motion phenomena, including various periodic motion states. Our ultimate goal is to maintain the motion within a single-cycle state. A bifurcation diagram within this range shows that the system moves in one period, and the maximum Liapunov curve shows that the LLE value at this time is less than 0.

8. Conclusions

In this study, we conducted in-depth research on a type of friction system. By applying fractional calculus to viscoelastic systems, we have successfully described their unique dynamic characteristics. We derived the equations for system amplitude and frequency using a nonlinear approximate analytical algorithm (averaging method) and verified the effectiveness of the analytical solution through power series expansion. By comparing the amplitude frequency curves of FOS and NFOS, we found that adding fractional order terms can effectively reduce the vibration amplitude, reduce the nonlinear characteristics of the system, and narrow down the range of unstable solutions of the system.
We analyzed in detail the effects of damping parameters, stiffness parameters, fractional order parameters, and disturbance parameters on the amplitude frequency curve. Research has found that changes in system parameter values can lead to significant changes in amplitude frequency curves and stability, providing a theoretical basis for material selection in system production. In addition, we investigated the local static bifurcation response and obtained the law of bifurcation characteristics changing with the opening parameters. Furthermore, we investigated the dynamic bifurcation behavior under different parameter perturbations, and the results showed that changing the system parameters would cause the system to alternate between single-period motion, multi-period motion, and chaotic motion. Therefore, by controlling the system parameters, we can maintain the system within the periodic motion range, thereby improving the stability of the system.
In order to better describe the response of the braking system and improve its performance, we propose the following suggestions:
  • Adjust fractional order parameters to optimize the vibration characteristics of the system, thereby reducing unnecessary energy loss.
  • Based on the amplitude frequency curve, select damping and stiffness parameters reasonably to ensure that the system remains stable under predetermined operating conditions.
By controlling the disturbance parameters, the system can avoid entering an unstable or chaotic state and improve the reliability and performance of the braking system.
The advantages and disadvantages of using fractional derivatives for viscoelastic materials in dynamic models are summarized as follows:
  • Advantages: (1) Fractional calculus can more accurately describe the memory characteristics and time dependence of viscoelastic materials. (2) Fractional calculus provides a more flexible tool to simulate the complex dynamic behavior of different materials.
  • Disadvantages: (1) Compared with integer order calculus, fractional order calculus has more complex calculations, especially in numerical calculations. This may lead to a decrease in computational efficiency, especially when processing large-scale data or conducting real-time simulations. (2) The theoretical foundation of fractional calculus is still relatively weak. The approximate processing of fractional calculus may lead to errors in the results, and further verification and improvement are needed in future research.

Author Contributions

Conceptualization, data curation, visualization, writing—original draft, methodology, and formal analysis: W.M. Project administration, investigation, and supervision: Q.D. Software, writing—review and editing, and validation: W.L. Funding acquisition and resources: Z.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Jilin Province Science and Technology Development Plan Project (20240302038GX) and the Beijing Enterprise Horizontal Project (3R2205502419).

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Conflicts of Interest

Author Zhenqi Yang was employed by the company State Grid Jibei Electric Power Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Mahdi, H.; Akgun, S.A.; Saleh, S.; Dautenhahn, K. A survey on the design and evolution of social robots—Past, present and future. Robot. Auton. Syst. 2022, 156, 104193. [Google Scholar] [CrossRef]
  2. Shah, A. Emerging trends in robotic aided additive manufacturing. In Proceedings of the International Conference on Additive Manufacturing and Advanced Materials (AM2), Gandhinagar, India, 6 January 2022. [Google Scholar] [CrossRef]
  3. Zbiss, K.; Kacem, A.; Santillo, M.; Mohammadi, A. Automatic Collision-Free Trajectory Generation for Collaborative Robotic Car-Painting. IEEE Access 2022, 10, 9950–9959. [Google Scholar] [CrossRef]
  4. Ju, H.; Park, H.; Kim, N.; Lim, J.; Jung, D.; Lee, J. A Locally Actuatable Soft Robotic Film for Actively Reconfiguring Shapes of Flexible Electronics. Soft Robot. 2022, 9, 767–775. [Google Scholar] [CrossRef] [PubMed]
  5. Volpe, G.; Cohen, S.; Capps, R.C.; Giacomelli, B.; McManus, R.; Scheckelhoff, K.; Choudhary, K.; Dabestani, A.T.; Hermann, S.; Kuiper, S.; et al. Robotics in acute care hospitals. Am. J. Health Syst. Pharm. 2012, 69, 1601–1603. [Google Scholar] [CrossRef]
  6. Oran, I.B.; Cezayirlioglu, H.R. AI—Robotic Applications in Logistics Industry and Savings Calculation. J. Organ. Behav. Res. 2021, 6, 148–165. [Google Scholar] [CrossRef]
  7. Ghiani, L.; Sassu, A.; Palumbo, F.; Mercenaro, L.; Gambella, F. In-Field Automatic Detection of Grape Bunches under a Totally Uncontrolled Environment. Sensors 2021, 21, 3908. [Google Scholar] [CrossRef] [PubMed]
  8. Liu, X.-F.; Zhang, X.-Y.; Cai, G.-P.; Chen, W.-J. Capturing a Space Target Using a Flexible Space Robot. Appl. Sci. 2022, 12, 984. [Google Scholar] [CrossRef]
  9. Bocii, L.S.; Sanjuan, J.J.V. Determination of the Braking Characteristics in Case of the Variation Depending on the Temperature of the Material Properties of the Friction Coupling of the Disc Brake of Railway Vehicles. Acta Tech. Napoc. Ser. Appl. Math. Mech. Eng. 2022, 65, 125–134. [Google Scholar]
  10. Zhang, W.; Zhao, C.; Chen, P.; Chen, E.; Lei, T. Numerical simulation and experimental research on mechanical behaviour of hydraulic disc brakes based on multi-body dynamics. Sci. Rep. 2022, 12, 18594. [Google Scholar] [CrossRef]
  11. Zhou, S.; Wang, Q.; Liu, J. Control Strategy and Simulation of the Regenerative Braking of an Electric Vehicle Based on an Electromechanical Brake. Trans. Famena 2022, 46, 23–40. [Google Scholar] [CrossRef]
  12. Yuan, Q. Study on the Influence of Damping Shim on Friction Squeal Characteristics of Automobile Disc Brakes. Noise Vib. Control 2022, 42, 201–207. [Google Scholar]
  13. Liu, S.H.; Wang, Z.T.; Ren, Y.R.; Ning, K.Y.; Yang, L.L. Research on full-stroke transfer coefficient of brake’s ball-plate forcing mechanism. J. Mach. Des. 2022, 39, 7–12. [Google Scholar]
  14. Suo, R.; Shi, X. Temperature Field and Stress Field Distribution of Forged Steel Brake Disc for High speed Train. Jordan J. Mech. Ind. Eng. 2022, 16, 113–121. [Google Scholar]
  15. Christian, S. Canadian space robotic activities. Acta Astronaut. 2004, 41, 239–246. [Google Scholar] [CrossRef]
  16. Verzijden, P. ERA performance measurements test results. In Proceedings of the 7th ESA Workshop on Advanced Space Technologies for Robotics and Automation ‘ASTRA 2002′ ESTEC, Noordwijk, The Netherlands, 19–21 November 2002. [Google Scholar]
  17. Shao, Z.Y.; Sun, H.X.; Jia, Q.X. Development of a general 2-dof space module. In Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems, Beijing, China, 9–15 October 2006; pp. 1002–1007. [Google Scholar] [CrossRef]
  18. Podlubny, I. Chapter 3—Existence and Uniqueness Theorems. In Mathematics in Science and Engineering; Elsevier: Amsterdam, The Netherlands, 1999; Volume 198, pp. 121–136. [Google Scholar] [CrossRef]
  19. Podlubny, I. Chapter 10—Survey of Applications of the Fractional Calculus. In Mathematics in Science and Engineering; Elsevier: Amsterdam, The Netherlands, 1999; Volume 198, pp. 261–307. [Google Scholar] [CrossRef]
  20. Podlubny, I. Chapter 9—Fractional-order Systems and Controllers. In Mathematics in Science and Engineering; Elsevier: Amsterdam, The Netherlands, 1999; Volume 198, pp. 243–260. [Google Scholar] [CrossRef]
  21. Song, Y.; Wang, H.; Chang, Y.; Li, Y. Nonlinear creep model and parameter identification of mudstone based on a modified fractional viscous body. Environ. Earth Sci. 2019, 78, 607. [Google Scholar] [CrossRef]
  22. Di Paola, M.; Alotta, G.; Burlon, A.; Failla, G. A novel approach to nonlinear variable-order fractional viscoelasticity. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2020, 378, 20190296. [Google Scholar] [CrossRef] [PubMed]
  23. Mezhoud, R.; Saoudi, K.; Zaraï, A.; Abdelmalek, S. Conditions for the local and global asymptotic stability of the time–fractional Degn–Harrison system. Int. J. Nonlinear Sci. Numer. Simul. 2020, 21, 749–759. [Google Scholar] [CrossRef]
  24. Birs, I.; Muresan, C.; Nascu, I.; Ionescu, C. A survey of recent advances in fractional order control for time delay systems. IEEE Access 2019, 7, 30951–30965. [Google Scholar] [CrossRef]
  25. Yavari, M.; Nazemi, A. On fractional infinite-horizon optimal control problems with a combination of conformable and Caputo–Fabrizio fractional derivatives. ISA Trans. 2020, 101, 78–90. [Google Scholar] [CrossRef]
  26. Pourhashemi, A.; Ramezani, A.; Siahi, M. Dynamic Fractional-Order Sliding Mode Strategy to Control and Stabilize Fractional-Order Nonlinear Biological Systems. IETE J. Research 2020, 68, 2560–2570. [Google Scholar] [CrossRef]
  27. Tran, N.; Van Au, V.; Zhou, Y.; Tuan, N.H. On a final value problem for fractional reaction-diffusion equation with Riemann-Liouville fractional derivative. Math. Methods Appl. Sci. 2020, 43, 3086–3098. [Google Scholar] [CrossRef]
  28. Marco, A.; Prabakaran, B.; Ivan, B. Anisotropic fractional viscoelastic constitutive models for human descending thoracic aortas. J. Mech. Behav. Biomed. Mater. 2019, 99, 186–197. [Google Scholar] [CrossRef]
  29. Lin, Z.L.; Zhang, J.Y.; He, L.Z. Method of Multiple Low-order Harmonic Currents Suppression Based on Fractional-order Capacitor. In Proceedings of the Chinese Society For Electrical Engineering, Wuhan, China, 9 February 2022; Volume 42, pp. 8921–8932. [Google Scholar] [CrossRef]
  30. Zhuang, B.; Cui, B.T.; Lou, X.Y.; Chen, J. Backstepping-based Output Feedback Boundary Control for Coupled Fractional Reaction-diffusion Systems. Acta Autom. Sin. 2022, 48, 2729–2743. [Google Scholar] [CrossRef]
  31. Zhu, H.H.; Zhu, J.Z.; Li, S.L.; Fan, J.W.; Done, H.J.; Wu, W.L. Power flow calculation and voltage analysis of fractional order power system. Electr. Mach. Control 2022, 26, 38–46. [Google Scholar] [CrossRef]
  32. Lu, C.H.; Bai, H.B. Study on Constitutive Model of Viscoelastic Material. Polym. Mater. Sci. Eng. 2007, 23, 28–31. [Google Scholar]
  33. Dai, J.; Han, M.; Ang, K.K. Moving element analysis of partially filled freight trains subject to abrupt braking. Int. J. Mech. Sci. 2019, 151, 85–94. [Google Scholar] [CrossRef]
  34. Wu, Q.; Luo, S.; Cole, C. Longitudinal dynamics and energy analysis for heavy haul trains. J. Mod. Transp. 2014, 22, 127–136. [Google Scholar] [CrossRef]
  35. Mastroddi, F.; Martarelli, F.; Eugeni, M.; Riso, C. Time- and frequency-domain linear viscoelastic modeling of highly damped aerospace structures. Mech. Syst. Signal Process. 2019, 122, 42–55. [Google Scholar] [CrossRef]
  36. He, Y.; Li, H.B.; Du, J. Dynamic Modulus Fitting Method of Viscoelastic Materials Based on Complex Neural Network. Chin. Q. Mech. 2022, 43, 406–415. [Google Scholar] [CrossRef]
  37. Li, M.M.; Fang, B.; Zhen, Y.X.; Zhao, J.X. A hybrid vibration isolator based on piezoelectric and viscoelastic materials. J. Vib. Shock 2017, 36, 134–140. [Google Scholar] [CrossRef]
  38. Tang, Z.H.; Luo, G.H.; Chen, W.; Yang, G.; Fang, J. Parallel dynamic model of rubber isolator about five-parameter fractional derivatives. J. Aerodyn. 2013, 28, 275–282. [Google Scholar] [CrossRef]
  39. Chang, Y.J.; Tian, W.W.; Chen, E.L.; Shen, Y.J.; Xing, W.C. Dynamic model for the nonlinear hysteresis of metal rubber based on the fractional-order derivative. Vib. Shock 2020, 39, 233–241. [Google Scholar] [CrossRef]
  40. Hamed, Y.S.; Albogamy, K.M.; Sayed, M. Nonlinear vibrations control of a contact-mode AFM model via a time-delayed positive position feedback. AEJ Alex. Eng. J. 2021, 60, 963–977. [Google Scholar] [CrossRef]
  41. Kandil, A.; Sayed, M.; Saeed, N.A. On the nonlinear dynamics of constant stiffness coefficients 16-pole rotor active magnetic bearings system. Eur. J. Mech. A/Solids 2020, 84, 104051. [Google Scholar] [CrossRef]
  42. Ms, A.; Aamb, C.; Imd, E. Stability and bifurcation analysis of a buckled beam via active control. Appl. Math. Model. 2020, 82, 649–665. [Google Scholar] [CrossRef]
  43. Wang, K. Stability and approximate solution of nonlinear dynamic system of a cylinder with two end faces in relative rotation. Acta Phys. Sin. 2005, 54, 5530–5533. [Google Scholar] [CrossRef]
  44. Wang, K.; Guan, X.P.; Qiao, J.M. Precise periodic solutions and uniqueness of periodic solutions of some relative rotation nonlinear dynamic system. Acta Phys. Sin. 2010, 59, 3648–3653. [Google Scholar] [CrossRef]
  45. Podlubny, I. Chapter 2—Fractional Derivatives and Integrals. In Mathematics in Science and Engineering; Elsevier: Amsterdam, The Netherlands, 1999; Volume 198, pp. 41–119. [Google Scholar] [CrossRef]
Figure 1. Functional structure diagram of new electromagnetic brake for robotic arm with vibration isolation layer.
Figure 1. Functional structure diagram of new electromagnetic brake for robotic arm with vibration isolation layer.
Fractalfract 08 00565 g001
Figure 2. Different states of the new electromagnetic brake of the robotic arm. (a) Electromagnet power-off braking state. (b) Electromagnet power-on starting state.
Figure 2. Different states of the new electromagnetic brake of the robotic arm. (a) Electromagnet power-off braking state. (b) Electromagnet power-on starting state.
Fractalfract 08 00565 g002
Figure 3. Nonlinear dynamic model of the new electromagnetic brake of a robotic arm.
Figure 3. Nonlinear dynamic model of the new electromagnetic brake of a robotic arm.
Fractalfract 08 00565 g003
Figure 4. System amplitude–frequency response curve.
Figure 4. System amplitude–frequency response curve.
Fractalfract 08 00565 g004
Figure 5. Comparison diagram of analytical and numerical solutions.
Figure 5. Comparison diagram of analytical and numerical solutions.
Fractalfract 08 00565 g005
Figure 6. The effect of K pairs on the amplitude and frequency relation curves of FOS (p = 0.33) and NFOS (p = 0).
Figure 6. The effect of K pairs on the amplitude and frequency relation curves of FOS (p = 0.33) and NFOS (p = 0).
Fractalfract 08 00565 g006
Figure 7. The effect of K on the amplitude–frequency curves at p = 0.9 and p = 0.1.
Figure 7. The effect of K on the amplitude–frequency curves at p = 0.9 and p = 0.1.
Fractalfract 08 00565 g007
Figure 8. Amplitude–frequency curve when k1 changes (K = 1, p = 0.33).
Figure 8. Amplitude–frequency curve when k1 changes (K = 1, p = 0.33).
Fractalfract 08 00565 g008
Figure 9. Amplitude–frequency curve when k2 changes (K = 1, p = 0.33).
Figure 9. Amplitude–frequency curve when k2 changes (K = 1, p = 0.33).
Fractalfract 08 00565 g009
Figure 10. Amplitude–frequency curve when k3 changes (K = 1, p = 0.33).
Figure 10. Amplitude–frequency curve when k3 changes (K = 1, p = 0.33).
Fractalfract 08 00565 g010
Figure 11. Amplitude–frequency curve when k4 changes (K = 1, p = 0.33).
Figure 11. Amplitude–frequency curve when k4 changes (K = 1, p = 0.33).
Fractalfract 08 00565 g011
Figure 12. Amplitude–frequency curve when c1 changes (K = 1, p = 0.33).
Figure 12. Amplitude–frequency curve when c1 changes (K = 1, p = 0.33).
Fractalfract 08 00565 g012
Figure 13. Amplitude–frequency curve when Fn changes (K = 1, p = 0.33).
Figure 13. Amplitude–frequency curve when Fn changes (K = 1, p = 0.33).
Fractalfract 08 00565 g013
Figure 14. Amplitude frequency curve of c2 variation range 0.005-0.02 (K = 1, p = 0.33).
Figure 14. Amplitude frequency curve of c2 variation range 0.005-0.02 (K = 1, p = 0.33).
Fractalfract 08 00565 g014
Figure 15. Amplitude frequency curve of c2 variation range 0.02-0.035 (K = 1, p = 0.33).
Figure 15. Amplitude frequency curve of c2 variation range 0.02-0.035 (K = 1, p = 0.33).
Fractalfract 08 00565 g015
Figure 16. Set of lag points, red represents G, green represents the first derivative of G with respect to a, and blue represents the second derivative of G with respect to a.
Figure 16. Set of lag points, red represents G, green represents the first derivative of G with respect to a, and blue represents the second derivative of G with respect to a.
Fractalfract 08 00565 g016
Figure 17. Set of bifurcation points, red represents G, green represents the first derivative of G with respect to a, and yellow represents the first derivative of G with respect to ω .
Figure 17. Set of bifurcation points, red represents G, green represents the first derivative of G with respect to a, and yellow represents the first derivative of G with respect to ω .
Fractalfract 08 00565 g017
Figure 18. (a) Transition set. (b) Bifurcation diagram.
Figure 18. (a) Transition set. (b) Bifurcation diagram.
Fractalfract 08 00565 g018
Figure 19. When c2 = 0.01, ω = 1.8, and K = 0.5. (a) Bifurcation diagram; (b) Largest Lyapunov Exponent (LLE); (c) phase diagram at Fn = 2; (d) Poincaré sections at Fn = 2; (e) phase diagram at Fn = 2.5; (f) Poincaré sections at Fn = 2.5; (g) phase diagram at Fn = 4; (h) Poincaré sections at Fn = 4; (i) phase diagram at Fn = 5; and (j) Poincaré sections at Fn = 5.
Figure 19. When c2 = 0.01, ω = 1.8, and K = 0.5. (a) Bifurcation diagram; (b) Largest Lyapunov Exponent (LLE); (c) phase diagram at Fn = 2; (d) Poincaré sections at Fn = 2; (e) phase diagram at Fn = 2.5; (f) Poincaré sections at Fn = 2.5; (g) phase diagram at Fn = 4; (h) Poincaré sections at Fn = 4; (i) phase diagram at Fn = 5; and (j) Poincaré sections at Fn = 5.
Fractalfract 08 00565 g019aFractalfract 08 00565 g019b
Figure 20. When c2 = 0.01, ω = 1.8, and K = 1.5. (a) Bifurcation diagram; (b) LLE; (c) phase diagram at Fn = 2; (d) Poincaré sections at Fn = 2; (e) phase diagram at Fn = 2.55; (f) Poincaré sections at Fn = 2.55; (g) phase diagram at Fn = 2.6; (h) Poincaré sections at Fn = 2.6; (i) phase diagram at Fn = 2.71; (j) Poincaré sections at Fn = 2.71; (k) phase diagram at Fn = 2.88; (l) Poincaré sections at Fn = 2.88; (m) phase diagram at Fn = 3.5; (n) Poincaré sections at Fn = 3.5; (o) phase diagram at Fn = 3.8; (p) Poincaré sections at Fn = 3.8; (q) phase diagram at Fn = 4.35; (r) Poincaré sections at Fn = 4.35.
Figure 20. When c2 = 0.01, ω = 1.8, and K = 1.5. (a) Bifurcation diagram; (b) LLE; (c) phase diagram at Fn = 2; (d) Poincaré sections at Fn = 2; (e) phase diagram at Fn = 2.55; (f) Poincaré sections at Fn = 2.55; (g) phase diagram at Fn = 2.6; (h) Poincaré sections at Fn = 2.6; (i) phase diagram at Fn = 2.71; (j) Poincaré sections at Fn = 2.71; (k) phase diagram at Fn = 2.88; (l) Poincaré sections at Fn = 2.88; (m) phase diagram at Fn = 3.5; (n) Poincaré sections at Fn = 3.5; (o) phase diagram at Fn = 3.8; (p) Poincaré sections at Fn = 3.8; (q) phase diagram at Fn = 4.35; (r) Poincaré sections at Fn = 4.35.
Fractalfract 08 00565 g020aFractalfract 08 00565 g020b
Figure 21. When c2 = 0.01, ω = 1.8, and K = 2. (a) Bifurcation diagram; (b) LLE; (c) phase diagram at Fn = 2; (d) Poincaré sections at Fn = 2; (e) phase diagram at Fn = 3.5; (f) Poincaré sections at Fn = 3.5; (g) phase diagram at Fn = 4; (h) Poincaré sections at Fn = 4.
Figure 21. When c2 = 0.01, ω = 1.8, and K = 2. (a) Bifurcation diagram; (b) LLE; (c) phase diagram at Fn = 2; (d) Poincaré sections at Fn = 2; (e) phase diagram at Fn = 3.5; (f) Poincaré sections at Fn = 3.5; (g) phase diagram at Fn = 4; (h) Poincaré sections at Fn = 4.
Fractalfract 08 00565 g021aFractalfract 08 00565 g021b
Figure 22. When Fn = 0.5 and p = 0.33. (a) Bifurcation diagram; (b) LLE.
Figure 22. When Fn = 0.5 and p = 0.33. (a) Bifurcation diagram; (b) LLE.
Fractalfract 08 00565 g022
Table 1. The motion state when Fn takes different values.
Table 1. The motion state when Fn takes different values.
FnKPoincareMotion State
20.5Three mapping pointsThree cycles
2.50.5Mapping point setChaos
40.5Mapping point setChaos
50.5Mapping point setChaos
Table 2. Each typical point and motion state when Fn takes different values.
Table 2. Each typical point and motion state when Fn takes different values.
FnKPoincareMotion State
21.5A mapping pointSingle cycle
2.551.5Mapping point setChaos
2.61.5Three-period mappingThree cycles
2.711.5Mapping point setChaos
2.881.5Biperiodic mappingTwo cycles
3.51.5Mapping point setChaos
3.81.5Three-period mappingThree cycles
4.351.5Mapping point setChaos
Table 3. Each typical point and motion state when Fn takes different values.
Table 3. Each typical point and motion state when Fn takes different values.
FnKPoincareMotion State
22A mapping pointSingle cycle
3.52Mapping point setChaos
42Biperiodic mappingTwo cycles
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ma, W.; Du, Q.; Li, W.; Yang, Z. Theoretical Analysis of Viscoelastic Friction System Characteristics of Robotic Arm Brake Based on Fractional Differential Theory. Fractal Fract. 2024, 8, 565. https://doi.org/10.3390/fractalfract8100565

AMA Style

Ma W, Du Q, Li W, Yang Z. Theoretical Analysis of Viscoelastic Friction System Characteristics of Robotic Arm Brake Based on Fractional Differential Theory. Fractal and Fractional. 2024; 8(10):565. https://doi.org/10.3390/fractalfract8100565

Chicago/Turabian Style

Ma, Wenli, Qiaoling Du, Wenhao Li, and Zhenqi Yang. 2024. "Theoretical Analysis of Viscoelastic Friction System Characteristics of Robotic Arm Brake Based on Fractional Differential Theory" Fractal and Fractional 8, no. 10: 565. https://doi.org/10.3390/fractalfract8100565

APA Style

Ma, W., Du, Q., Li, W., & Yang, Z. (2024). Theoretical Analysis of Viscoelastic Friction System Characteristics of Robotic Arm Brake Based on Fractional Differential Theory. Fractal and Fractional, 8(10), 565. https://doi.org/10.3390/fractalfract8100565

Article Metrics

Back to TopTop