# Study on Chaotic Peculiarities of Magnetic-Mechanical Coupled System of Giant Magnetostrictive Actuator

Authors: Hongbo Yan; Enzuo Liu; Pengbo Zhao; Pei Liu; Rui Cao
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

We studied the chaotic peculiarities of magnetic-mechanical coupled system of GMA. Based on the working principle of GMA and according to Newton’s second law of motion, first piezomagnetic equation, disk spring design theory, and structural dynamics principle of GMA, the present study established a GMA magnetic-mechanical coupled system model. By carrying out data modeling of this coupled system model, the bifurcation chart of the system with the variation of damping factor, excitation force, and exciting frequency parameters as well as the homologous offset oscillogram, phase plane trace chart, and Poincaré diagram was obtained, and the chaotic peculiarities of the system were analyzed. The influence of parametric errors on the coupled system was studied. The analytical results showed that the oscillation equation of the GMA magnetic-mechanical coupled system had nonlinearity and the movement morphology was complicated and diversified. By adjusting the damping factor, exciting frequency, and excitation force parameters of the system, the system could work under the stable interval, which provided theoretical support for the stability design of GMA.


---

## Body

## 1. Introduction

Due to its advantages of large displacement, fast response speed, large magnetic-mechanical coupling coefficient, and high reliability [1, 2], the giant magnetostrictive material (GMM) has a broad application prospect in aerospace, national defense and military industry, machinery, precision machining, medical equipment, and other fields [3–8]. The giant magnetostrictive actuator (GMA) is becoming an outstanding representative of new driving mechanism because its core component GMM has superior performance, which can better realize the mutual transformation of magnetic-mechanical energy.In recent years, many scholars have studied the magnetic-mechanical coupled system. Yoo et al. [9] used the three-dimensional finite element method to carry out analysis of the GMA magnetic-mechanical coupled system and verified the reliability of this method through experiments. Niu et al. [10] proposed a GMA dynamic model based on the J-A model, which considered the magnetic-mechanical coupling effect and could improve the forecast accuracy of GMA output displacement. Ebrahimi et al. [11] established a magnetic-mechanical coupled system model by conducting change calculation of continuous energy function for the minimized system. This model studied the relationship between stress and magnetic conductivity. Yang et al. [12] established a magnetic-mechanical strong coupled system model for the force sensor made by GMM and compared it with the measuring result to examine the effectiveness of the model. Zhao [13] carried out experimental research on the dynamic and static peculiarities of GMM under the impact of a magnetic-mechanical coupled system. Li [14] proposed a GMA hysteresis nonlinear model based on the magnetic-mechanical coupling principle and carried out experimental analysis on the dynamic and static peculiarities of GMA force and output displacement.However, as a result of the hysteresis nonlinearity in the input and output relation of GMA, it not only makes GMA easy to produce positioning errors in the open-loop system but also causes unstable closed-loop system in complex tracking problems, which seriously affected the adhibition of GMA [15]. Therefore, in order to improve the performance and enhance the actual control and application of GMA, enhancing the stability of the GMA nonlinear system has become a hot topic in recent years.If the GMA magnetic-mechanical coupled system is under the working condition of chaotic movement, the actual control and stability of the GMA nonlinear system will be seriously damaged. Therefore, the study of chaotic peculiarities of the magnetic-mechanical coupled system has occupied a very important position in the stability of the GMA nonlinear system.Zhu et al. [16] found that under the action of in-plane harmonic wave and random excitation, the giant magnetostrictive thin film shape memory alloy composite cantilever slab system had random bifurcation and chaos phenomena with the change of parameters. Xu et al. [17] discussed that the giant magnetostrictive piezoelectric composite sensor system had multifrequency response in the harmonic magnetic scene, and the system movement changed from periodical orbit to chaos when increasing the external excitation. Bowong and Kakmeni [18] proposed a robust control scheme with simple structure to suppress the chaotic movement of the system. Sarmah and Ananthakrishna [19] established a magnetic-strain coupling nonlinear oscillator model, which effectively predicted the chaotic quasiperiodic path of the system. Sun [20] established a mechanical model of a GMA system based on the piezomagnetic theory and analyzed the chaotic phenomenon that may exist in the model. Yan et al. [21] established the mathematical model of a GMA hysteresis oscillatory system and studied the bifurcation and chaos phenomena of the GMA system. At present, there are fewer research achievements of parameters on the chaotic peculiarities of the GMA magnetic-mechanical coupled system and the systematic theoretical basis is insufficient. However, the extensive research has been made on the chaotic peculiarities of other nonlinear systems [22–26].First of all, based on the GMA working principle, this paper set up the GMA magnetic-mechanical coupled system model and obtained the coupled equation of this system. Then, it used the Melnikov method to analyze the threshold value of chaos phenomenon generated by the system. In addition, on the basis of bifurcation chart, offset oscillogram, phase plane trace chart, and Poincaré diagram, it explained the chaotic peculiarities of this coupled system. Finally, the influence of parametric errors on the coupled system was studied. The research result is beneficial to optimize the stability design and improvement plan of the GMA magnetic-mechanical coupled system.


## 2. Working Principle of the GMA

The GMA working principle is as follows: when the alternating current is implemented in the driving coil, the alternating magnetic scene will be generated. Under the effect of the driving magnetic scene, the GMM rod has extension changes along the direction of the rod, thus generating the output displacement. The bias current is applied in the bias coil, which can generate the bias magnetic field so as to avoid the occurrence of GMA frequency-doubled effect and improve the sensitivity and linear peculiarities of GMA. The pretensioner is composed of an outer wall sleeve, disk spring, upper end cover, and output rod, which can offer some prestressing forces so that the GMM rod can work in an optimal state. The closed magnetic circuit is mainly composed of upper and lower magnetic conductive loops, flux sleeve, and flux sheet, which can lower magnetic leakage and effectively alleviate the jam of magnetic scene to external equipment. Figure1 shows the GMA structure chart.Figure 1
GMA structure. 1, base; 2, outer wall socket; 3, flux sleeve; 4, drive coil; 5, upper end cover; 6, lower magnetic conductive ring; 7, coil former; 8, flux sheet; 9,GMM rod; 10, bias coil; 11, upper magnetic conductive ring; 12, disc ring; 13, output rod.


## 3. GMA Magnetic-Mechanical Coupled System Model

We can see from the GMA working principle that under the impact of steady-state magnetic scene, the displacement of upper end face of GMM rod is the same as that of the output rod. In the movement process, GMM output forceT, disk spring prepressure Px, damping force of output rod Qx˙, and external excitation Fcosωt of output rod have a combined action on the GMA output rod. Its equivalent mechanical model is shown in Figure 2, where x is the displacement of the output rod.Figure 2
Equivalent mechanical model of GMA output rod.In light of Newton’s second law of motion, the dynamic differential equation of the GMA magnetic-mechanical coupled system is as follows:(1)T−P−Q−Fcosωt=mx¨,where m is the quality of the output rod.According to the first piezomagnetic equation, when the excitation magnetic scene is uniform to the axial direction of the rod, the simplified piezomagnetic equation along the direction of the rod can be indicated as follows [27, 28] by ignoring the influence of shear stress, transverse stress, and transverse excitation:(2)ζ=S33Hσ+d33H3,(3)B=d33σ+μ33TH3,where ζ, σ, and H3 are the strain intensity, stress intensity, and magnetic field intensity of GMM; d33 is the piezomagnetic coefficient; S33H is the compliant mechanism; B is the magnetic scene intensity; and μ33T is the magnetic conductivity.The GMA system has been analyzed and the output forceT of GMM rod is as follows:(4)T=EAζ,where E and A are the resilience modulus and cross-sectional area of GMM rod.Under the condition of steady-state magnetic scene and ignoring magnetic leakage flux, the magnetic scene strengthH3 of GMM is as follows according to the magnetomotive force of magnetic circuit:(5)H3=NIμ33TPA−d33σμ33T,where N and I are the count of turns and current of coil and P is the total magnetic resistance of magnetic circuit.It is obtained from the simultaneous formulas (2)–(5) that the output force T of GMM rod is(6)T=EAS33H−d332μ33Tσ+d33NIμ33TPA.Under the impact of steady-state magnetic scene, the number of turns and current of coil and the cross-sectional area of the GMM rod are all constants, and when the magnetic circuit is the closed circuit, the GMM magnetic resistance is only related to its relative magnetic conductivity and cross-sectional area, so it can be regarded as the constant. Without considering the influence of system magnetic resistance and magnetic scene changes on various parameters [29], the GMM rod’s output force T is(7)T=EAl1−Ed332μ33Tx,where l is the GMM rod’s length.We can see from the disk spring design theory [30] that the prepressure of the disk spring P is(8)P=4E1−μ2d4K1D2K42fzdK42hd−fzdhd−fz2d+1,where μ is Poisson’s ratio and d and D are the thickness and outer diameter of disk spring. K1 and K4 are determined by disk spring design handbook. fz is the deformation amount of combination disk spring. If fz=2f, f is the deformation amount of a single piece disk spring. fz in this paper is taken as x, and h is the calculated value of normal deformation amount of disk spring.It can be obtained through simplifying formula (8):(9)P=kx+u1x2+u2x3,where(10)k=4EK44d3K1D21−μ21K42+hd2,u1=−6EK44hdK1D21−μ2,u2=2EK44dK1D21−μ2.Through the analysis of the GMA system, it can be seen that the damping forceQ of output rod is as follows:(11)Q=cx˙,where c is the damping factor of the GMA system.Formulas (7), (9), and (11) are substituted into formula (1) to obtain the equation of the GMA magnetic-mechanical coupled system:(12)mx¨+cx˙+kex+u1x2+u2x3−Fcosωt=0,where(13)ke=k+EAlEd332μ33T−1.Dividing both sides of formula (12) by m, the formula can be sorted out as follows:(14)x¨+Cx˙+αx+ω0x2+βx3−F0cosωt=0,where(15)C=cm,α=kem,ω0=u1m,β=u2m,F0=Fm.


## 4. Chaotic Peculiarity Analysis of GMA Magnetic-Mechanical Coupled System

When carrying out chaotic response analysis of the GMA magnetic-mechanical coupled system, formula (14) needs to be simplified into the standard form. If x˙=y, C=εr, F0=εf, and 0<ε≪1, the standard form of formula (14) shall be(16)x˙=y,y˙=−αx−ω0x2−βx3+ε−rx˙+fcosωt.The function is defined as follows:(17)lx=αx+ω0x2+βx3,mx˙,x¨,ωt=−rx˙+fcosωt.Whenε=0, formula (16) is(18)x˙=y,y˙=−αx−ω0x2−βx3.Whenε=0, system (18) is the Hamilton system without disturbance, and its Hamilton action shall be(19)h=12y2+α2x2+ω03x3+β4x4.It can be obtained from formula (18) that its balance points are O(0, 0), A−ω0+ω02−4αω0/2β,0, and B−ω0−ω02−4αω0/2β,0.The featured equation of the system shall be(20)λ2=−α−2ω0x2−3βx2.The balance point (0, 0) is substituted into formula (20) to obtain λ=±iα, which are the conjugate complex roots. Therefore, the balance point O (0, 0) is the center. Similarly, the balance point B−ω0−ω02−4αω0/2β,0 is also the center while A−ω0+ω02−4αω0/2β,0 is the saddle point. As a result, the system has the homoclinic orbit I that passes through the saddle point A−ω0+ω02−4αω0/2β,0, as shown in Figure 3.Figure 3
Homoclinic orbit.Therefore, the parameter equation of homoclinic orbit is(21)xt=±22ω0e±ω0t/3β3βe±2ω0t/3β+1+ω03β,yt=±22ω021−e±2ω0t/3β9β5/2e±2ω0t/3β+12.The Melnikov function of homoclinic orbit is(22)Mt0=∫−∞+∞mxt,x˙t,ωt+t0ytdt=∫−∞+∞−rx˙t+fcosωt+t0ytdt.Formula (21) is substituted into formula (22) to obtain the integral:(23)Mt0=4ω03r81β5/2+22πωfe−2πωβ1/2/ω0sinωt0β1/21−cosh3πωβ1/2/ω0.WhenMt0=0,(24)sinωt0=−2ω03r1−cosh3πωβ1/2/ω081πωfβ2e−2πωβ1/2/ω0.Assinωt0≤1,(25)−2ω03r1−cosh3πωβ1/2/ω081πωfβ2e−2πωβ1/2/ω0≤1.r and f are substituted into formula (25). According to the Melnikov method, as we can see from Mt0=0, the threshold condition that may make the system generate chaos in the sense of Smale horseshoe transformation is as follows:(26)2ω03C1−cosh3πωβ1/2/ω081πωF0β2e−2πωβ1/2/ω0≤1.Whenε is small enough, if Mt0 has a straightforward null point within the range of t0∈0,N, the system may have the bifurcation of homoclinic orbit, which means that the chaotic behavior may occur.


## 5. Numerical Modeling

For equation (14), each parameter value is shown in Table 1. The fourth-order Runge–Kutta method was used to solve the equation of the GMA magnetic-mechanical coupled system. Based on the numerical solution, the bifurcation chart of the system under different parameters was obtained, and the chaotic peculiarities of the system under the changes of damping factor, excitation force, and exciting frequency were studied combined with the offset oscillogram, phase plane trace chart, and Poincaré diagram.Table 1
Each parameter value of the GMA system.
ParameterValueEquivalent damping factor of GMA systemC0.1Equivalent stiffness coefficient of GMA systemα−0.2Modulus of secondary rigidity term of disk springω00.1Modulus of third rigidity term of disk springβ1Excitation force of GMA systemF0.22Exciting frequency of GMA systemω1.2


### 5.1. The Influence of Damping Factor on the Chaotic Peculiarities of the System

Figure4 shows the bifurcation chart of the GMA magnetic-mechanical coupled system changing with the damping factor of the system. The initial condition is [0, 0, 0.06]. When the damping factor of the system changes, the path to chaos of the system is as follows:(27)Chaotic movement⟶intermittent chaos⟶period2⟶periodic movementFigure 4
The bifurcation chart of the GMA magnetic-mechanical coupled system changing with the damping factor of the system.We can see from Figure4 that when the damping factor of the system is ∈(0.06, 0.14), the system carried out chaotic movement, as shown in Figure 5 (c = 0.1). When the damping factor of the system is (0.14, 0.32), the system degrades into period 2 movement with inverse bifurcation, as shown in Figure 6 (c = 0.2). When the damping factor of the system is 0.2, the offset oscillogram is regularly distributed and has a constant period. The phase trajectory chart is a closed curve. Poincaré diagram shows two attractive points, so the system carries out a typical doubly periodic movement and has intermittent chaos when the damping factor of the system is around 0.16. When the damping factor of the system is ∈(0.32, 1.00), the system carries out the steady-state movement with the period of 1, as shown in Figure 7 (c = 0.5), that is, when the damping factor of the system is 0.5, the offset oscillogram is regularly distributed and has a constant period. The phase trajectory chart is a closed curve. Poincaré diagram shows one attractive point, so the system carries out steady-state movement with the single period. Table 2 shows the form of movement of damping factor in different areas.Figure 5
(a) The offset oscillogram, (b) phase trajectory chart, and (c) Poincaré diagram of the system whenC = 0.1.
(a)(b)(c)Figure 6
(a) The offset oscillogram, (b) phase trajectory chart, and (c) Poincaré diagram of the system whenC = 0.2.
(a)(b)(c)Figure 7
(a) The offset oscillogram, (b) phase trajectory chart, and (c) Poincaré diagram of the system whenC = 0.5.
(a)(b)(c)Table 2
The form of movement of damping factor in different areas.
Ranges ofCMovement form(0.06, 0.14)Chaotic movement (Figure5)(0.14, 0.32)Period 2 (Figure6)(0.32, 1.00)Periodic movement (Figure7)Figure5 shows the offset oscillogram, phase trajectory chart, and Poincaré diagram of the system when the damping factor of the system is 0.1. It can be seen that the offset oscillogram in Figure 5(a) has an irregular distribution with an unstable period. The phase trajectory chart in Figure 5(b) cannot be closed and has filled the phase space area for a long time. The phase trajectory movement is the reciprocating movement with an infinite period. In Figure 5(c), the Poincaré diagram is neither a limited point set nor a closed curve. Therefore, the system can be determined to carry out chaotic movement.


### 5.2. The Influence of Excitation Force on the Chaotic Peculiarities of the System

Figure8 shows the bifurcation chart of the GMA magnetic-mechanical coupled system changing with the excitation force. The initial condition is [0 0 0.02]. When the excitation force alters, the path to chaos of the system is as follows:(28)Periodic movement⟶Period2⟶chaotic movement⟶multiple intermittent chaost⟶period2⟶periodic movement⟶chaotic movement⟶multiple intermittent chaos⟶quasi−periodic movement⟶periodic movementFigure 8
The bifurcation chart of the GMA magnetic-mechanical coupled system changing with the excitation force.We can see from Figure8 that when the excitation force is ∈(0.02, 0.045), the system carries out the steady-state movement with the period of 1. When the excitation force is ∈(0.045, 0.14), the system carries out the movement with the period of 2 accompanied with the tangent bifurcation, as shown in Figure 9 (F = 0.1). When the excitation force is 0.1, the offset oscillogram is regular with a constant period, while the phase trajectory chart is a closed curve. Poincaré diagram shows two attractive points, so the system carries out a typical doubly periodic movement. When the excitation force is ∈(0.14, 0.315), the system carries out the chaotic movement as shown in Figure 10 (F = 0.3) and has multiple intermittent chaos around 0.155 and 0.315. When the excitation force is ∈(0.315, 0.36), the system degrades into period 2 movement with inverse bifurcation. When the excitation force is ∈(0.36, 0.39), the system carries out the movement with the period of 1. When the excitation force is ∈(0.39, 0.59), the system carries out the chaotic movement as shown in Figure 11 (F = 0.4) and has multiple intermittent chaos around 0.52 and 0.545. When the excitation force is ∈(0.59, 0.74), the system carries out quasiperiodic movement and has multiple intermittent chaos around 0.605 and 0.68. When the excitation force is ∈(0.74, 1.00), the system carries out the steady-state periodic movement. Table 3 shows the form of movement of excitation force in different areas.Figure 9
(a) The offset oscillogram, (b) phase trajectory chart, and (c) Poincaré diagram of the system whenF = 0.1.
(a)(b)(c)Figure 10
(a) The offset oscillogram, (b) phase trajectory chart, and (c) Poincaré diagram of the system whenF = 0.3.
(a)(b)(c)Figure 11
(a) The offset oscillogram, (b) phase trajectory chart, and (c) Poincaré diagram of the system whenF = 0.4.
(a)(b)(c)Table 3
The form of movement of excitation force in different areas.
Ranges ofFMovement form(0.02, 0.045)Periodic movement(0.045, 0.14)Period 2 (Figure9)(0.14, 0.315)Chaotic movement (Figure10)(0.315, 0.36)Period 2(0.36, 0.39)Periodic movement(0.39, 0.59)Chaotic movement (Figure11)(0.59, 0.74)Quasiperiodic movement(0.74, 1.00)Periodic movementFigures10 and 11 show the offset oscillogram, phase plane trace chart, and Poincaré diagram of the system when the excitation forces are 0.3 and 0.4. At this time, it can be obtained that the offset oscillograms in Figures 10(a) and 11(a) have an irregular distribution with an unstable period. The phase trajectory charts in Figures 10(b) and 11(b) fill the phase space area and are neither overlapping nor entangled. The Poincaré diagrams in Figures 10(c) and 11(c) are neither a limited point set nor a closed curve, which indicates that the system is in the state of chaotic movement.


### 5.3. The Influence of Exciting Frequency on the Chaotic Peculiarities of the System

Figure12 shows the bifurcation chart of the GMA magnetic-mechanical coupled system changing with the exciting frequency. The initial condition is 000.12. When the exciting frequency changes, the path to chaos of the system is as follows:(29)Chaotic movement⟶intermittent chaos⟶periodic movement⟶chaotic movement⟶intermittent chaos⟶periodic movement⟶skip movement⟶quasi−periodic movement⟶chaotic movement⟶multiple intermittent chaos⟶periodic movementFigure 12
The bifurcation chart of the GMA magnetic-mechanical coupled system changing with the exciting frequency.We can see from Figure12 that when the exciting frequency is ∈(0.12, 0.31), the system carries out the chaotic movement. When the exciting frequency is ∈(0.31, 0.36), the system carries out the movement with the period of 1. When the exciting frequency is ∈(0.36, 0.42), the system carries out the chaotic movement, as shown in Figure 13 (ω = 0.39), and has the intermittent chaos at the exciting frequency of 0.395. When the exciting frequency is ∈(0.42, 0.48), the system carries out the steady-state movement with the single period. When the exciting frequency is greater than 0.48, the system carries out the skip movement. When the exciting frequency is ∈(0.48, 0.525), the system skips from the single period movement to quasiperiodic movement. When the exciting frequency is ∈(0.525, 0.655), the system carries out the chaotic movement, as shown in Figure 14 (ω = 0.6). When the exciting frequency is ∈(0.655, 1.00), the system carries out the steady-state movement with the single period, as shown in Figure 15 (ω = 0.8), and has the intermittent chaos at the excitation frequencies around 0.665, 0.69, and 0.77. Table 4 shows the form of movement of exciting frequency in different areas.Figure 13
(a) The offset oscillogram, (b) phase trajectory chart, and (c) Poincaré diagram of the system whenω = 0.39.
(a)(b)(c)Figure 14
(a) The offset oscillogram, (b) phase trajectory chart, and (c) Poincaré diagram of the system whenω = 0.6.
(a)(b)(c)Figure 15
(a) The offset oscillogram, (b) phase trajectory chart, and (c) Poincaré diagram of the system whenω = 0.8.
(a)(b)(c)Table 4
The form of system movement of exciting frequency in different areas.
Ranges ofωMovement form(0.12, 0.31)Chaotic movement(0.31, 0.36)Periodic movement(0.36, 0.42)Chaotic movement (Figure13)(0.42, 0.48)Periodic movement(0.48, 0.525)Quasiperiodic movement(0.525, 0.655)Chaotic movement (Figure14)(0.655, 1.00)Periodic movement (Figure15)Figures13 and 14 show the offset oscillogram, phase plane trace chart, and Poincaré diagram of the system when the excitation frequencies are 0.39 and 0.6. It can be obtained that the offset oscillograms in Figures 13(a) and 14(a) have an irregular distribution with an unstable period. The phase trajectory charts in Figures 13(b) and 14(b) fill the phase space region and are neither overlapping nor entangled. The Poincaré diagrams in Figures 13(c) and 14(c) are neither a finite point set nor a closed curve, which proves that the system is in the state of chaotic movement.We can see from Figure15 that when the exciting frequency is 0.8, the offset oscillogram has a regular distribution with a constant period. The phase trajectory chart is a closed curve, and Poincaré diagram shows an attractive point, which indicates that the system carries out a typical steady-state movement with the single period.


## 5.1. The Influence of Damping Factor on the Chaotic Peculiarities of the System

Figure4 shows the bifurcation chart of the GMA magnetic-mechanical coupled system changing with the damping factor of the system. The initial condition is [0, 0, 0.06]. When the damping factor of the system changes, the path to chaos of the system is as follows:(27)Chaotic movement⟶intermittent chaos⟶period2⟶periodic movementFigure 4
The bifurcation chart of the GMA magnetic-mechanical coupled system changing with the damping factor of the system.We can see from Figure4 that when the damping factor of the system is ∈(0.06, 0.14), the system carried out chaotic movement, as shown in Figure 5 (c = 0.1). When the damping factor of the system is (0.14, 0.32), the system degrades into period 2 movement with inverse bifurcation, as shown in Figure 6 (c = 0.2). When the damping factor of the system is 0.2, the offset oscillogram is regularly distributed and has a constant period. The phase trajectory chart is a closed curve. Poincaré diagram shows two attractive points, so the system carries out a typical doubly periodic movement and has intermittent chaos when the damping factor of the system is around 0.16. When the damping factor of the system is ∈(0.32, 1.00), the system carries out the steady-state movement with the period of 1, as shown in Figure 7 (c = 0.5), that is, when the damping factor of the system is 0.5, the offset oscillogram is regularly distributed and has a constant period. The phase trajectory chart is a closed curve. Poincaré diagram shows one attractive point, so the system carries out steady-state movement with the single period. Table 2 shows the form of movement of damping factor in different areas.Figure 5
(a) The offset oscillogram, (b) phase trajectory chart, and (c) Poincaré diagram of the system whenC = 0.1.
(a)(b)(c)Figure 6
(a) The offset oscillogram, (b) phase trajectory chart, and (c) Poincaré diagram of the system whenC = 0.2.
(a)(b)(c)Figure 7
(a) The offset oscillogram, (b) phase trajectory chart, and (c) Poincaré diagram of the system whenC = 0.5.
(a)(b)(c)Table 2
The form of movement of damping factor in different areas.
Ranges ofCMovement form(0.06, 0.14)Chaotic movement (Figure5)(0.14, 0.32)Period 2 (Figure6)(0.32, 1.00)Periodic movement (Figure7)Figure5 shows the offset oscillogram, phase trajectory chart, and Poincaré diagram of the system when the damping factor of the system is 0.1. It can be seen that the offset oscillogram in Figure 5(a) has an irregular distribution with an unstable period. The phase trajectory chart in Figure 5(b) cannot be closed and has filled the phase space area for a long time. The phase trajectory movement is the reciprocating movement with an infinite period. In Figure 5(c), the Poincaré diagram is neither a limited point set nor a closed curve. Therefore, the system can be determined to carry out chaotic movement.


## 5.2. The Influence of Excitation Force on the Chaotic Peculiarities of the System

Figure8 shows the bifurcation chart of the GMA magnetic-mechanical coupled system changing with the excitation force. The initial condition is [0 0 0.02]. When the excitation force alters, the path to chaos of the system is as follows:(28)Periodic movement⟶Period2⟶chaotic movement⟶multiple intermittent chaost⟶period2⟶periodic movement⟶chaotic movement⟶multiple intermittent chaos⟶quasi−periodic movement⟶periodic movementFigure 8
The bifurcation chart of the GMA magnetic-mechanical coupled system changing with the excitation force.We can see from Figure8 that when the excitation force is ∈(0.02, 0.045), the system carries out the steady-state movement with the period of 1. When the excitation force is ∈(0.045, 0.14), the system carries out the movement with the period of 2 accompanied with the tangent bifurcation, as shown in Figure 9 (F = 0.1). When the excitation force is 0.1, the offset oscillogram is regular with a constant period, while the phase trajectory chart is a closed curve. Poincaré diagram shows two attractive points, so the system carries out a typical doubly periodic movement. When the excitation force is ∈(0.14, 0.315), the system carries out the chaotic movement as shown in Figure 10 (F = 0.3) and has multiple intermittent chaos around 0.155 and 0.315. When the excitation force is ∈(0.315, 0.36), the system degrades into period 2 movement with inverse bifurcation. When the excitation force is ∈(0.36, 0.39), the system carries out the movement with the period of 1. When the excitation force is ∈(0.39, 0.59), the system carries out the chaotic movement as shown in Figure 11 (F = 0.4) and has multiple intermittent chaos around 0.52 and 0.545. When the excitation force is ∈(0.59, 0.74), the system carries out quasiperiodic movement and has multiple intermittent chaos around 0.605 and 0.68. When the excitation force is ∈(0.74, 1.00), the system carries out the steady-state periodic movement. Table 3 shows the form of movement of excitation force in different areas.Figure 9
(a) The offset oscillogram, (b) phase trajectory chart, and (c) Poincaré diagram of the system whenF = 0.1.
(a)(b)(c)Figure 10
(a) The offset oscillogram, (b) phase trajectory chart, and (c) Poincaré diagram of the system whenF = 0.3.
(a)(b)(c)Figure 11
(a) The offset oscillogram, (b) phase trajectory chart, and (c) Poincaré diagram of the system whenF = 0.4.
(a)(b)(c)Table 3
The form of movement of excitation force in different areas.
Ranges ofFMovement form(0.02, 0.045)Periodic movement(0.045, 0.14)Period 2 (Figure9)(0.14, 0.315)Chaotic movement (Figure10)(0.315, 0.36)Period 2(0.36, 0.39)Periodic movement(0.39, 0.59)Chaotic movement (Figure11)(0.59, 0.74)Quasiperiodic movement(0.74, 1.00)Periodic movementFigures10 and 11 show the offset oscillogram, phase plane trace chart, and Poincaré diagram of the system when the excitation forces are 0.3 and 0.4. At this time, it can be obtained that the offset oscillograms in Figures 10(a) and 11(a) have an irregular distribution with an unstable period. The phase trajectory charts in Figures 10(b) and 11(b) fill the phase space area and are neither overlapping nor entangled. The Poincaré diagrams in Figures 10(c) and 11(c) are neither a limited point set nor a closed curve, which indicates that the system is in the state of chaotic movement.


## 5.3. The Influence of Exciting Frequency on the Chaotic Peculiarities of the System

Figure12 shows the bifurcation chart of the GMA magnetic-mechanical coupled system changing with the exciting frequency. The initial condition is 000.12. When the exciting frequency changes, the path to chaos of the system is as follows:(29)Chaotic movement⟶intermittent chaos⟶periodic movement⟶chaotic movement⟶intermittent chaos⟶periodic movement⟶skip movement⟶quasi−periodic movement⟶chaotic movement⟶multiple intermittent chaos⟶periodic movementFigure 12
The bifurcation chart of the GMA magnetic-mechanical coupled system changing with the exciting frequency.We can see from Figure12 that when the exciting frequency is ∈(0.12, 0.31), the system carries out the chaotic movement. When the exciting frequency is ∈(0.31, 0.36), the system carries out the movement with the period of 1. When the exciting frequency is ∈(0.36, 0.42), the system carries out the chaotic movement, as shown in Figure 13 (ω = 0.39), and has the intermittent chaos at the exciting frequency of 0.395. When the exciting frequency is ∈(0.42, 0.48), the system carries out the steady-state movement with the single period. When the exciting frequency is greater than 0.48, the system carries out the skip movement. When the exciting frequency is ∈(0.48, 0.525), the system skips from the single period movement to quasiperiodic movement. When the exciting frequency is ∈(0.525, 0.655), the system carries out the chaotic movement, as shown in Figure 14 (ω = 0.6). When the exciting frequency is ∈(0.655, 1.00), the system carries out the steady-state movement with the single period, as shown in Figure 15 (ω = 0.8), and has the intermittent chaos at the excitation frequencies around 0.665, 0.69, and 0.77. Table 4 shows the form of movement of exciting frequency in different areas.Figure 13
(a) The offset oscillogram, (b) phase trajectory chart, and (c) Poincaré diagram of the system whenω = 0.39.
(a)(b)(c)Figure 14
(a) The offset oscillogram, (b) phase trajectory chart, and (c) Poincaré diagram of the system whenω = 0.6.
(a)(b)(c)Figure 15
(a) The offset oscillogram, (b) phase trajectory chart, and (c) Poincaré diagram of the system whenω = 0.8.
(a)(b)(c)Table 4
The form of system movement of exciting frequency in different areas.
Ranges ofωMovement form(0.12, 0.31)Chaotic movement(0.31, 0.36)Periodic movement(0.36, 0.42)Chaotic movement (Figure13)(0.42, 0.48)Periodic movement(0.48, 0.525)Quasiperiodic movement(0.525, 0.655)Chaotic movement (Figure14)(0.655, 1.00)Periodic movement (Figure15)Figures13 and 14 show the offset oscillogram, phase plane trace chart, and Poincaré diagram of the system when the excitation frequencies are 0.39 and 0.6. It can be obtained that the offset oscillograms in Figures 13(a) and 14(a) have an irregular distribution with an unstable period. The phase trajectory charts in Figures 13(b) and 14(b) fill the phase space region and are neither overlapping nor entangled. The Poincaré diagrams in Figures 13(c) and 14(c) are neither a finite point set nor a closed curve, which proves that the system is in the state of chaotic movement.We can see from Figure15 that when the exciting frequency is 0.8, the offset oscillogram has a regular distribution with a constant period. The phase trajectory chart is a closed curve, and Poincaré diagram shows an attractive point, which indicates that the system carries out a typical steady-state movement with the single period.


## 6. Coupled System with Parametric Errors

In order to study the influence of parametric uncertainties on the coupled system, the parameters can be compared with the random errors within the range of ±20%, which is similar to the method adopted in [31–33]. Take Figure 15(b) as an example; when ω is the random error within the range of ±20%, as shown in Figure 16, the amplitude of displacement decreases by 0.18%, and a similar periodic orbit also exists. The results show that the coupled system can stay on a periodic orbit with more sensitive parameters. It indicates that the time interval of ω leads the system to periodicity motion. This is also similar to the conclusion from Tusset et al. [31]. Jia et al. [34] designed one self-adaptive algorithm based on the Preisach model and processed dynamic adjustment at the parameters; this algorithm validly compensated nonlinear hysteresis of output displacement of the giant magnetostrictive actuator, reduced the parameter errors of the system, and also improved the control precision of the actuator.Figure 16
Coupled system without and with parametric uncertainties.To make giant magnetostrictive magnetic-mechanical coupled system enter into periodicity motion, we can lead a passive rotational magnetorheological (MR) damper into the system. The MR damper restrains the chaos action of the system through adding different currents. The research in [35] found that the simple introduction of a passive rotational MR damper without electric current has a smaller effect on the chaos action of the system, so we could use the rotational MR damper with energizing discontinuity to make the coupled system keep periodicity vibration. Tusset et al. [31] processed the state-dependent Riccati equation (SDRE) control and the nonlinear saturation control (NSC) which are related to the status of the system; it could also make the coupled system enter into periodicity status.


## 7. Conclusion

(1)
Chaotic peculiarities of the GMA magnetic-mechanical coupled system were studied, and the value range of each parameter in chaotic movement of the system was obtained through numerical modeling: chaotic movement occurs when the damping factor of the system is within the interval of (0.06, 0.14) or when the excitation force is within the intervals of (0.14, 0.315) and (0.39, 0.59) or when the exciting frequency is within the intervals of (0.12, 0.31), (0.36, 0.42) and (0.525, 0.655).(2)
By controlling the damping factor, excitation force, and exciting frequency parameters of the system, the GMA magnetic-mechanical coupled system can avoid chaotic movement purposefully, thus improving the stability of GMA.(3)
The chaotic peculiarities of the GMA magnetic-mechanical coupled system were subjected to theoretical analysis, which could effectively avoid the occurrence of chaotic movement, improve the stability of the GMA system, and contribute to the wide application of GMA.



---

*Source: 6864795-2020-11-25.xml*