On the mixtures of MGT viscoelastic solids

Noelia Bazarra; José R. Fernández; Ramón Quintanilla; Noelia Bazarra; José R. Fernández; Ramón Quintanilla

doi:10.3934/era.2022219

Electronic Research Archive

2022, Volume 30, Issue 12: 4318-4340. doi: 10.3934/era.2022219

Previous Article Next Article

Research article Special Issues

On the mixtures of MGT viscoelastic solids

1.
Departamento de Matemática Aplicada I, Universidade de Vigo, Escola de Enxeñería de Telecomunicación, Campus As Lagoas Marcosende s/n, Vigo 36310, Spain
2.
Departamento de Matemáticas, E.S.E.I.A.A.T.-U.P.C., Colom 11, Terrassa 08222, Barcelona, Spain

Received: 08 July 2022 Revised: 14 September 2022 Accepted: 21 September 2022 Published: 28 September 2022

In this paper, we study, from both analytical and numerical points of view, a problem involving a mixture of two viscoelastic solids. An existence and uniqueness result is proved using the theory of linear semigroups. Exponential decay is shown for the one-dimensional case. Then, fully discrete approximations are introduced using the finite element method and the implicit Euler scheme. Some a priori error estimates are obtained and the linear convergence is derived under suitable regularity conditions. Finally, one- and two-dimensional numerical simulations are presented to demonstrate the convergence, the discrete energy decay and the behavior of the solution.

Keywords:

Citation: Noelia Bazarra, José R. Fernández, Ramón Quintanilla. On the mixtures of MGT viscoelastic solids[J]. Electronic Research Archive, 2022, 30(12): 4318-4340. doi: 10.3934/era.2022219

Related Papers:

[1]	Qiaoping Li, Sanyang Liu . Predefined-time vector-polynomial-based synchronization among a group of chaotic systems and its application in secure information transmission. AIMS Mathematics, 2021, 6(10): 11005-11028. doi: 10.3934/math.2021639
[2]	Minghung Lin, Yiyou Hou, Maryam A. Al-Towailb, Hassan Saberi-Nik . The global attractive sets and synchronization of a fractional-order complex dynamical system. AIMS Mathematics, 2023, 8(2): 3523-3541. doi: 10.3934/math.2023179
[3]	Sukono, Siti Hadiaty Yuningsih, Endang Rusyaman, Sundarapandian Vaidyanathan, Aceng Sambas . Investigation of chaos behavior and integral sliding mode control on financial risk model. AIMS Mathematics, 2022, 7(10): 18377-18392. doi: 10.3934/math.20221012
[4]	Abdulaziz Khalid Alsharidi, Saima Rashid, S. K. Elagan . Short-memory discrete fractional difference equation wind turbine model and its inferential control of a chaotic permanent magnet synchronous transformer in time-scale analysis. AIMS Mathematics, 2023, 8(8): 19097-19120. doi: 10.3934/math.2023975
[5]	Omar Kahouli, Imane Zouak, Ma'mon Abu Hammad, Adel Ouannas . Chaos, control and synchronization in discrete time computer virus system with fractional orders. AIMS Mathematics, 2025, 10(6): 13594-13621. doi: 10.3934/math.2025612
[6]	Pratap Anbalagan, Evren Hincal, Raja Ramachandran, Dumitru Baleanu, Jinde Cao, Michal Niezabitowski . A Razumikhin approach to stability and synchronization criteria for fractional order time delayed gene regulatory networks. AIMS Mathematics, 2021, 6(5): 4526-4555. doi: 10.3934/math.2021268
[7]	Pratap Anbalagan, Evren Hincal, Raja Ramachandran, Dumitru Baleanu, Jinde Cao, Chuangxia Huang, Michal Niezabitowski . Delay-coupled fractional order complex Cohen-Grossberg neural networks under parameter uncertainty: Synchronization stability criteria. AIMS Mathematics, 2021, 6(3): 2844-2873. doi: 10.3934/math.2021172
[8]	Xinna Mao, Hongwei Feng, Maryam A. Al-Towailb, Hassan Saberi-Nik . Dynamical analysis and boundedness for a generalized chaotic Lorenz model. AIMS Mathematics, 2023, 8(8): 19719-19742. doi: 10.3934/math.20231005
[9]	Canhong Long, Zuozhi Liu, Can Ma . Synchronization dynamics in fractional-order FitzHugh–Nagumo neural networks with time-delayed coupling. AIMS Mathematics, 2025, 10(4): 8673-8687. doi: 10.3934/math.2025397
[10]	Honglei Yin, Bo Meng, Zhen Wang . Disturbance observer-based adaptive sliding mode synchronization control for uncertain chaotic systems. AIMS Mathematics, 2023, 8(10): 23655-23673. doi: 10.3934/math.20231203

Abstract

1. Introduction

Since pioneering works of Pecora and Carroll's ^[1], chaos synchronization and control have turned a hot topic and received much attention in various research areas. A number of literatures shows that chaos synchronization can be widely used in physics, medicine, biology, quantum neuron and engineering science, particularly in secure communication and telecommunications ^[1,2,3]. In order to realize synchronization, experts have proposed lots of methods, including complete synchronization and $Q$ - $S$ synchronization ^[4,5], adaptive synchronization ^[6], lag synchronization^[7,8], phase synchronization ^[9], observer-based synchronization ^[10], impulsive synchronization ^[11], generalized synchronization ^[12,13], lag projective synchronization ^[14,15], cascade synchronization et al ^{[16,17,18,19,20]}. Among them, the cascade synchronization method is a very effective algorithm, which is characterized by reproduction of signals in the original chaotic system to monitor the synchronized motions.

It is know that, because of the complexity of fractional differential equations, synchronization of fractional-order chaotic systems is more difficult but interesting than that of integer-order systems. Experts find that the key space can be enlarged by the regulating parameters in fractional-order chaotic systems, which enables the fractional-order chaotic system to be more suitable for the use of the encryption and control processing. Therefore, synchronization of fractional-order chaotic systems has gained increasing interests in recent decades ^{[21,22,23,24,25,26,27,28,29,30,31]}. It is noticed that most synchronization methods mentioned in ^{[4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20]} work for integer-order chaotic systems. Here, we shall extend to cascade synchronization for integer-order chaotic systems to a kind of general form, namely function cascade synchronization (FCS), which means that one chaotic system may be synchronized with another by sending a signal from one to the other wherein a scaling function is involved. The FCS is effective both for the fractional order and integer order chaotic systems. It constitutes a general method, which can be considered as a continuation and extension of earlier works of ^[13,16,19]. The nice feature of our method is that we introduce a scaling function for achieving synchronization of fractional-order chaotic systems, which can be chosen as a constant, trigonometric function, power function, logarithmic and exponential function, hyperbolic function and even combinations of them. Hence, our method is more general than some existing methods, such as the complete synchronization approach and anti-phase synchronization approach et al.

To sum up, in this paper, we would like to use the FCS approach proposed to study the synchronization of fractional-order chaotic systems. We begin our theoretical work with the Caputo fractional derivative. Then, we give the FCS of the fractional-order chaotic systems in theory. Subsequently, we take the fractional-order unified chaotic system as a concrete example to test the effectiveness of our method. Finally, we make a short conclusion.

2. The method of FCS for fractional-order chaotic systems

As for the fractional derivative, there exists a lot of mathematical definitions ^[32,33]. Here, we shall only adopt the Caputo fractional calculus, which allows the traditional initial and boundary condition assumptions. The Caputo fractional calculus is described by

$\begin{eqnarray} \frac{d^q f(t)}{dt^q} = \frac{1}{{\Gamma}(q-n)}\int^t_0 \frac{f^{(n)}(\xi)}{(t-\xi)^{q-n+1}}d\xi, \quad\quad n-1 \lt q \lt n . \end{eqnarray}$

(2.1)

Here, we give the function cascade synchronization method to fractional-order chaotic systems. Take a fractional-order dynamical system:

$\begin{eqnarray} \frac{d^q x}{d t^q} = f(x) = L x+N( x) \end{eqnarray}$

(2.2)

as a drive system. In the above $x = (x_1, x_2, x_3)^T$ is the state vector, $f: R^3 \rightarrow R^3$ is a continuous function, $L x$ and $N(x)$ represent the linear and nonlinear part of $f(x)$ , respectively.

Firstly, on copying any two equations of $(2.2)$ , such as the first two, one will obtain a sub-response system:

$\begin{eqnarray} \frac{d^q y}{d t^q} = L_1 y+N_1(y, x_3)+\tilde{U} \end{eqnarray}$

(2.3)

with $y = (X_1, Z)^T$ . In the above, $x_3$ is a signal provided by (2.2), while $\tilde{U} = (u_1, u_2)^T$ is a controller to be devised.

For the purpose of realizing the synchronization, we now define the error vector function via

$\begin{eqnarray} \tilde{e} = y-\tilde{Q}(\tilde{x})\tilde{x} \end{eqnarray}$

(2.4)

where $\tilde{e} = (e_1, e_2)^T$ , $\tilde{x} = (x_1, x_2)^T$ and $\tilde{Q}(\tilde{x}) = {\rm diag}({Q}_1(x_1), {Q}_2(x_2))$ .

Definition 1. For the drive system (2.2) and response system (2.3), one can say that the synchronization is achieved with a scaling function matrix $\tilde{Q}(\tilde{x})$ if there exists a suitable controller $\tilde{U}$ such that

$\begin{eqnarray} \lim\limits_{t \to \infty }{||\tilde{e}||} = \lim\limits_{t \to \infty }{||y-\tilde{Q}(\tilde{x})\tilde{x}||} = 0. \end{eqnarray}$

(2.5)

Remark 1. We would like to point out that one can have various different choices on the scaling function $\tilde{Q}(\tilde{x})$ , such as constant, power function, trigonometric function, hyperbola function, logarithmic and exponential function, as well as limited quantities of combinations and composite of the above functions. Particularly, when $\tilde{Q}(\tilde{x}) = I$ and $-I$ ( $I$ being a unit matrix), the problem is reducible to the complete synchronization and anti-phase synchronization of fractional-order chaotic systems, respectively. When $\tilde{Q}(\tilde{x}) = \alpha I$ , it becomes to the project synchronization. And when $\tilde{Q}(\tilde{x})$ = ${\rm diag}(\alpha_1, \alpha_2)$ , it turns to the modified projective synchronization. Hence, our method is more general than the existing methods in ^[4,13].

It is noticed from (2.5) that the system (2.3) will synchronize with (2.2) if and only if the error dynamical system (2.5) is stable at zero. For this purpose, an appropriate controller $\tilde{U}$ such that (2.5) is asymptotical convergent to zero is designed, which is described in the following theorem.

Theorem 1. For a scaling function matrix $\tilde{Q}(\tilde{x})$ , the FCS will happen between (2.2) and (2.3) if the conditions:

(i) the controller $\tilde{U}$ is devised by

$\begin{eqnarray} \tilde{U} = \tilde{K}\tilde{e}-N_1(y, x_3)+\tilde{Q}(\tilde{x})N_1(\tilde{x})+\tilde{P}(\tilde{x})\tilde{x} \end{eqnarray}$

(2.6)

(ii) the matrix $\tilde{K}$ is a $2\times 2$ matrix such that

$\begin{eqnarray} L_1+\tilde{K} = -\tilde{C}, \end{eqnarray}$

(2.7)

are satisfied simultaneously. In the above, $\tilde{P}(\tilde{x}) = {\rm {diag}}(\dot{Q}_1(x_1)\frac{ d^q x_1}{d t^q}, \dot{Q}_2(x_2)\frac{ d^q x_2}{d t^q})$ , $\tilde{K}$ is a $2\times 2$ function matrix to be designed. While $\tilde{C} = (\tilde{C}_{ij})$ is a $2\times 2$ function matrix wherein

$\begin{eqnarray} \tilde{C}_{ii} \gt 0 \quad \quad{\rm and} \quad \quad \tilde{C}_{ij} = -\tilde{C}_{ji}, \quad i\neq j. \end{eqnarray}$

(2.8)

Remark 2. It needs to point out that the construction of the suitable controller $\tilde{U}$ plays an important role in realizing the synchronization between (2.2) and (2.3). Theorem 2 provides an effective way to design the controller. It is seen from the theorem that the controller $\tilde{U}$ is closely related to the matrix $\tilde{C}$ . Once the condition $(2.8$ ) is satisfied, one will has many choices on the controller $\tilde{U}$ .

Remark 3. Based on the fact that the fractional orders themselves are varying parameters and can be applied as secret keys when the synchronization algorithm is adopted in secure communications, it is believed that our method will be more suitable for some engineering applications, such as chaos-based encryption and secure communication.

Proof: Let's turn back to the error function given in (2.4). Differentiating this equation with respect to $t$ and on use of the first two equations of (2.2) and (2.3), one will obtain the following dynamical system

$\begin{eqnarray} \frac{d^q \tilde{e}}{d t^q}& = &\frac{d^q y}{d t^q}-\tilde{Q}(\tilde{x})\frac{d^q x}{d t^q}-\tilde{P}(\tilde{x})\tilde{x}\\ & = &L_1 y+N_1(y, x_3)+\tilde{U}-\tilde{Q}(\tilde{x})[L_1 \tilde{x}+N_1(\tilde{x})]-\tilde{P}(\tilde{x})\tilde{x}\\ & = &L_1\tilde{e}+N_1(y, x_3)-\tilde{Q}(\tilde{x})N_1(\tilde{x})-\tilde{P}(\tilde{x})\tilde{x}+\tilde{K}\tilde{e}-N_1(y, x_3)+\tilde{Q}(\tilde{x})N_1(\tilde{x})+\tilde{P}(\tilde{x})\tilde{x}\\ & = &(L_1+\tilde{K})\tilde{e}. \end{eqnarray}$

(2.9)

Assuming that $\lambda$ is an arbitrary eigenvalue of matrix $L_1+\tilde{K}$ and its eigenvector is recorded as $\eta$ , i.e.

$\begin{eqnarray} ( L_1+\tilde{K} )\eta = \lambda \eta, \quad \eta\neq 0. \end{eqnarray}$

(2.10)

On multiplying (2.10) by $\eta^{\mathbf {H}}$ on the left, we obtain that

$\begin{eqnarray} \eta^{\mathbf {H}}( L_1+\tilde{K} )\eta = \lambda \eta^{\mathbf {H}}\eta \end{eqnarray}$

(2.11)

where ${\mathbf {H}}$ denotes conjugate transpose. Since $\bar{\lambda}$ is also an eigenvalue of $L_1+\tilde{K}$ , we have that

$\begin{eqnarray} \eta^{\mathbf {H}}( L_1+\tilde{K} )^{\mathbf {H}} = \bar{\lambda} \eta^{\mathbf {H}}. \end{eqnarray}$

(2.12)

On multiplying (2.12) by $\eta$ on the right, we derive that

$\begin{eqnarray} \eta^{\mathbf {H}}( L_1+\tilde{K} )^{\mathbf {H}}\eta = \bar{\lambda} \eta^{\mathbf {H}}\eta \end{eqnarray}$

(2.13)

From (2.11) and (2.13), one can easily get that

$\begin{eqnarray} \lambda+\bar{\lambda}& = &\eta^{\mathbf {H}}[( L_1+\tilde{K} )^{\mathbf {H}}+( L_1+\tilde{K} )]\eta/\eta^{\mathbf {H}}\eta\\[0.06in] & = &-\eta^{\mathbf {H}}(\tilde{C}+\tilde{C}^{\mathbf {H}})\eta/\eta^{\mathbf {H}}\eta\\[0.06in] & = &-\eta^{\mathbf {H}}\Lambda\eta/\eta^{\mathbf {H}}\eta \end{eqnarray}$

(2.14)

with $\Lambda = \tilde{C}+\tilde{C}^{\mathbf {H}}$ . Since $\tilde{C}$ satisfy the condition (2.8), one can know that $\Lambda$ denotes a real positive diagonal matrix. Thus we have $\eta^{\mathbf {H}}\Lambda\eta > 0$ . Accordingly, we can get

$\begin{eqnarray} \lambda+\bar{\lambda} = 2\mathbf {Re}(\lambda) = -\eta^{\mathbf {H}}\Lambda\eta/\eta^{\mathbf {H}}\eta \lt 0, \end{eqnarray}$

(2.15)

which shows

$\begin{eqnarray} |\mathbf{arg}\lambda| \gt \frac{\pi}{2} \gt \frac{q\pi}{2}. \end{eqnarray}$

(2.16)

According to the stability theorem in Ref. ^[34], the error dynamical system (2.9) is asymptotically stable, i.e.

$\begin{eqnarray} \lim\limits_{t \to \infty }{||\tilde{e}||} = \lim\limits_{t \to \infty }{||y-\tilde{Q}(\tilde{x})\tilde{x}||} = 0, \end{eqnarray}$

(2.17)

which implies that synchronization can be achieved between (2.2) and (2.3). The proof is completed.

Next, on copying the last two equations of (2.2), one will get another sub-response system:

$\begin{eqnarray} \frac{d^q z}{d t^q} = L_2 z+N_2(z, X_1)+\bar{U} \end{eqnarray}$

(2.18)

where $X_1$ is a synchronized variable in (2.3), $z = (X_2, X_3)^T$ and $\bar{U} = (u_3, u_4)^T$ is the controller being designed.

Here, we make analysis analogous to the above. Now we define the error $\bar{e}$ via

$\begin{eqnarray} \tilde{e} = z-\bar{Q}(\bar{x})\bar{x} \end{eqnarray}$

(2.19)

where $\bar{e} = (e_3, e_4)^T$ , $\bar{x} = (x_2, x_3)^T$ and $\bar{Q}(\bar{x}) = {\rm diag}({Q}_3(x_2), {Q}_4(x_3))$ . If devising the the controller $\bar{U}$ as

$\begin{eqnarray} \bar{U} = \bar{K}\bar{e}-N_2(z, X_1)+\bar{Q}(\bar{x})N_2(\bar{x})+\bar{P}(\bar{x})\bar{x} \end{eqnarray}$

(2.20)

and $L_2+\bar{K}$ satisfying

$\begin{eqnarray} L_2+\bar{K} = -\bar{C} \end{eqnarray}$

(2.21)

where $\bar{P}(\bar{x}) = {\rm {diag}}(\dot{Q}_3(x_2)\frac{ d^q x_2}{d t^q}, \dot{Q}_4(x_3)\frac{ d^q x_3}{d t^q})$ , $\bar{C} = (\bar{C}_{ij})$ denotes a $2\times 2$ function matrix satisfying

$\begin{eqnarray} \bar{C}_{ii} \gt 0 \quad \quad{\rm and} \quad \quad \bar{C}_{ij} = -\bar{C}_{ji}, \quad i\neq j, \end{eqnarray}$

(2.22)

then the error dynamical system (2.19) satisfies

$\begin{eqnarray} \lim\limits_{t \to \infty }{||\bar{e}||} = \lim\limits_{t \to \infty }{||z-\bar{Q}(\bar{x})\bar{x}||} = 0. \end{eqnarray}$

(2.23)

Therefore, one achieve the synchronization between the system (2.2) and (2.18). Accordingly, from (2.5) and (2.23), one can obtain that

$\begin{eqnarray} \left\{\begin{array}{l} \lim\limits_{t \to \infty }{||X_1-Q_1(x_1)x_1||} = 0, \\[0.08in] \lim\limits_{t \to \infty }{||X_2-Q_3(x_2)x_2||} = 0, \\[0.08in] \lim\limits_{t \to \infty }{||X_3-Q_4(x_3)x_3||} = 0. \end{array}\right. \end{eqnarray}$

(2.24)

which indicates the FCS is achieved for the fractional order chaotic systems.

3. Applications of FCS method to fractional-order unified chaotic system

In the sequel, we shall extend the applications of FCS approach to the fractional-order unified chaotic system to test the effectiveness.

The fractional-order unified chaotic system is described by:

$\begin{eqnarray} \left\{\begin{array}{l} \\ \frac{d^{q} x_{1}}{d t^q} = (25a+10)(x_2-x_1), \\[0.1in] \\ \frac{d^q x_{2}}{d t^q} = (28-35a)x_1-x_1x_3+(29a-1)x_2, \\[0.1in] \\ \frac{d^q x_{3}}{d t^q } = x_1x_2-\frac{a+8}{3}x_3, \end{array}\right. \end{eqnarray}$

(3.1)

where $x_i, (i = 1, 2, 3)$ are the state parameters and $a\in[0, 1]$ is the control parameter. It is know that when $0\leq a < 0.8$ , the system (3.1) corresponds to the fractional-order Lorenz system ^[35]; when $a = 0.8$ , it is the Lü system ^[36]; while when $0.8 < a < 1$ , it turns to the Chen system ^[37].

According to the FCS method in section 2, we take (3.1) as the drive system. On copying the first two equation, we get a sub-response system of (3.1):

$\begin{eqnarray} \left\{\begin{array}{l} \\ \frac{d^q X_{1}}{d t^q} = (25a+10)(Z-X_1)+u_1, \\[0.12in] \\ \frac{d^q Z}{d t^q} = (28-35a)X_1-Z x_3+(29a-1)Z+u_2, \end{array}\right. \end{eqnarray}$

(3.2)

where $\tilde{U} = (u_1, u_2)^T$ is a controller to be determined. In the following, we need to devise the desired controller $\tilde{U}$ such that (3.1) can be synchronized with (3.2). For this purpose, we set the error function $\tilde{e} = (e_1, e_2)$ via :

$\begin{eqnarray} \tilde{e} = (e_1, e_2) = (X_1-x_1(x_1^2+\alpha_1), Z-x_2\tanh{x_2}). \end{eqnarray}$

(3.3)

On devising the controller $\tilde{U}$ as (2.6), one can get that the error dynamical system is

$\begin{eqnarray} \frac{d^q \tilde{e}}{d t^q} = (L_1+\tilde{K})\tilde{e}, \end{eqnarray}$

(3.4)

where

$\begin{equation} \begin{split} & L_1 = \left( \begin{array}{cc} -10-25a & -10-25a\\ 28-35a & 29a-1 \\ \end{array} \right), \end{split} \quad\quad\quad \begin{split} & N_1(y, x_3) = \left( \begin{array}{c} 0 \\ -X_1x_3 \\ \end{array} \right){.} \end{split} \end{equation}$

(3.5)

If choosing, for example, the matrix $\tilde{K}$ as

$\begin{equation} \begin{split} & \tilde{K} = \left( \begin{array}{cc} -\lambda_1+25a+10 & x_1+x_1x_2-25a\\ -x_1-x_1x_2+35a-38 & -\lambda_2-29a+1 \\ \end{array} \right), \end{split} \end{equation}$

(3.6)

where $\lambda_1 > 0$ and $\lambda_2 > 0$ , then one can obtain that

$\begin{equation} \begin{split} & \tilde{C} = \left( \begin{array}{cc} -\lambda_1 & x_1+x_1x_2+10\\ -x_1-x_1x_2-10 & -\lambda_2 \\ \end{array} \right). \end{split} \end{equation}$

(3.7)

Therefore the dynamical system (3.4) becomes

$\begin{eqnarray} \frac{d^q \tilde{e}}{d t^q} = \begin{split} &\left( \begin{array}{cc} -\lambda_1 & x_1+x_1x_2\\ -x_1-x_1x_2 & -\lambda_2 \\ \end{array} \right)\tilde{e}. \end{split} \end{eqnarray}$

(3.8)

According to Theorem 2, the synchronization is realized in the system (3.1) and (3.2).

Subsequently, on copying the last two equations of (3.1), we get another sub-response system:

$\begin{eqnarray} \left\{\begin{array}{l} \\ \frac{\partial^q X_{2}}{\partial t^q} = (28-35a)X_1-X_1 X_3+(29a-1)X_2+u_3, \\[0.15in] \\ \frac{\partial^q X_{3}}{\partial t^q} = X_1X_2-\frac{a+8}{3}X_3+u_4, \end{array}\right. \end{eqnarray}$

(3.9)

where $\bar{U} = ({u_3, \ u_4})^T$ is the controller needed. When choosing the error function $\bar{e} = (e_3, e_4)$ as:

$\begin{eqnarray} \bar{e} = (e_3, e_4) = (X_2-\alpha_2 x_2, X_3-x_3(\alpha_3+e^{-x_3})), \end{eqnarray}$

(3.10)

and the controller $\bar{U}$ as (2.20), where

$\begin{equation} \begin{split} & L_2 = \left( \begin{array}{cc} 29a-1 & 0\\ 0 & -\frac{a+8}{3} \\ \end{array} \right), \end{split} \quad\quad\quad \begin{split} & N_2(z, X_1) = \left( \begin{array}{c} -X_1X_3 \\ X_1X_2 \\ \end{array} \right), \end{split} \end{equation}$

(3.11)

and the matrix $\bar{K}$ is chosen by

$\begin{equation} \begin{split} & \bar{K} = \left( \begin{array}{cc} -\lambda_3-29a+1 & 1+x_2x_3+e^{-x_3}\\ -1-x_2x_3-e^{-x_3} & -\lambda_4-\frac{a+8}{3} \\ \end{array} \right), \end{split} \end{equation}$

(3.12)

where $\lambda_3 > 0$ and $\lambda_4 > 0$ . Calculations show that the error dynamical system (2.19) becomes

$\begin{eqnarray} \frac{d^q \bar{e}}{d t^q} = \begin{split} &\left( \begin{array}{cc} -\lambda_3 & 1+x_2x_3+e^{-x_3}\\ -1-x_2x_3-e^{-x_3} & -\lambda_4 \\ \end{array} \right)\bar{e}. \end{split} \end{eqnarray}$

(3.13)

which, according to the stability theorem, indicates that $\bar{e}$ will approach to zero with time evolutions. Therefore, the FCS is realized for the fractional-order unified chaotic system.

In the above, we have revealed that the FCS is achieved for the fractional-order unified chaotic system in theory. In the sequel, we shall show that the FCS is also effective in the numerical algorithm.

For illustration, we set the fractional order $q = 0.98$ and the parameters $\lambda_i (i = 1, \cdots, 4)$ as $(\lambda_1, \lambda_2, \lambda_3, \lambda_4) = (2, 3, 0.5, 0.3)$ . It is noticed that when the value of $a\in[0, 1]$ is given, the system (3.1) will be reduced to a concrete system. For example, when $a = 0.2$ , it corresponds to the fractional-order Lorenz system. The chaotic attractors are depicted in Figure 1. Time responses of states variables and synchronization errors of the Lorenz system are showed in and , respectively. When $a = 0.8$ , it is the fractional-order Lü system. The chaotic attractors, time responses of state variables and synchronization errors are exhibited in –, respectively. When $a = 0.95$ , it turns to the fractional-order Chen system. Numerical simulation results are depicted in Figures 7–9. From the chaotic attractors pictures marked by Figures 1, 4 and 5, one can easily see that the trajectories of the response system (colored red) display certain consistency to that of the drive system (colored black) because of the special scaling functions chosen. Meanwhile, one can also see the synchronization is realized from Figures 3, 6 and 9. Therefore, we conclude that the FCS is a very effective algorithm for achieving the synchronization of the fractional-order unified chaotic system.

Figure 1. FCS of the fractional-order Lorenz system. Here we choose

$(\alpha_1, \alpha_2, \alpha_3) = (0.2, 2, -1.5)$ , initial values

$(x_1, x_2, x_3) = (-1, -0.5, -0.2)$ and

$(X_1, X_2, X_3) = (0.2, 0.3, 0.1)$ .

DownLoad: Full-Size Img PowerPoint

Figure 2. Time responses of state variables

$x_i$ and

$X_i (i = 1, 2, 3)$ for the fractional-order Lorenz system.

DownLoad: Full-Size Img PowerPoint

Figure 3. Synchronization errors of the Lorenz system.

DownLoad: Full-Size Img PowerPoint

Figure 4. FCS of the fractional-order Lü system with

$a = 0.8$ . Here we choose

$(\alpha_1, \alpha_2, \alpha_3) = (-0.5, 2.5)$ , initial values

$(x_1, x_2, x_3) = (0.5, -0.5, 0.2)$ and

$(X_1, X_2, X_3) = (0.15, -0.1, -0.1)$ .

DownLoad: Full-Size Img PowerPoint

Figure 5. Time responses of state variables

$x_i$ and

$X_i (i = 1, 2, 3)$ for the Lü system.

DownLoad: Full-Size Img PowerPoint

Figure 6. Synchronization errors of the Lü system.

DownLoad: Full-Size Img PowerPoint

Figure 7. FCS of the fractional-order Chen system with

$a = 0.95$ . Here we choose

$(\alpha_1, \alpha_2, \alpha_3) = (0.5, 1.5, 3)$ , initial values

$(x_1, x_2, x_3) = (1.5, 0.02, -0.01)$ and

$(X_1, X_2, X_3) = (-2, 0.01, -0.05)$ .

DownLoad: Full-Size Img PowerPoint

Figure 8. Time evolutions of state variables

$x_i$ and

$X_i$

$(i = 1, 2, 3)$ for the Chen system.

DownLoad: Full-Size Img PowerPoint

Figure 9. Synchronization errors of the Chen system.

DownLoad: Full-Size Img PowerPoint

4. Conclusions

Chaos synchronization, because of the potential applications in telecommunications, control theory, secure communication et al, has attracted great attentions from various research fields. In the present work, via the stability theorem, we successfully extend the cascade synchronization of integer-order chaotic systems to a kind of general function cascade synchronization algorithm for fractional-order chaotic systems. Meanwhile, we apply the method to the fractional-order unified chaotic system for an illustrative test. Corresponding numerical simulations fully reveal that our method is not only accuracy, but also effective.

It is worthy of pointing out that the scaling function introduced makes the method more general than the complete synchronization, anti-phase synchronization, modified projective synchronization et al. Therefore, in this sense, our method is applicable and representative. However, the present work just study the fractional-order chaotic system without time-delay. It is known that in many cases the time delay is inevitably in the real engineering applications. Lag synchronization seems to be more practical and reasonable. Hence, it will be of importance and interest to study whether the FCS method can be used to realize the synchronization of fractional-order chaotic systems with time-delay. We shall considered it in our future work.

Acknowledgments

The authors would like to express their sincere thanks to the referees for their kind comments and valuable suggestions. This work is supported by the National Natural Science Foundation of China under grant No.11775116 and No.11301269.

Conflict of interest

We declare that we have no conflict of interests.

References

[1]	H. F. Tiersten, M. Jahanmir, A theory of composites modeled as inerpenetreting solid continua, Arch. Ration. Mech. Anal., 65 (1977), 153–192. https://doi.org/10.1007/BF00276554 doi: 10.1007/BF00276554
[2]	R. J. Atkin, R. E. Craine, Continuum theories of mixtures: basic theory and historical development, Q. J. Mech. Appl. Math., 29 (1976), 209–244. https://doi.org/10.1093/qjmam/29.2.209 doi: 10.1093/qjmam/29.2.209
[3]	R. M. Bowen, Theory of mixtures, in Continuum Physics III, Academic Press, New York, (1976), 689–722.
[4]	A. Bedford, D. S. Drumheller, Theories of immiscible and structured materials, Int. J. Eng. Sci., 21 (1983), 863–960. https://doi.org/10.1016/0020-7225(83)90071-X doi: 10.1016/0020-7225(83)90071-X
[5]	A. Bedford, M. Stern, A multi-continuum theory of composite elastic materials, Acta Mech., 14 (1972), 85–102. https://doi.org/10.1007/BF01184851 doi: 10.1007/BF01184851
[6]	A. Bedford, M. Stern, Towards a diffusing continuum theory of composite elastic materials, J. Appl. Mech., 38 (1971), 8–14. https://doi.org/10.1115/1.3408772 doi: 10.1115/1.3408772
[7]	R. M. Bowen, J. C. Wiese, Diffusion in mixtures of elastic materials, Int. J. Eng. Sci., 7 (1969), 689–722. https://doi.org/10.1016/0020-7225(69)90048-2 doi: 10.1016/0020-7225(69)90048-2
[8]	A. C. Eringen, D. J. Ingram, A continuum theory of chemically reacting media, Int. J. Eng. Sci., 3 (1965), 197–212. https://doi.org/10.1016/0020-7225(65)90044-3 doi: 10.1016/0020-7225(65)90044-3
[9]	A. E. Green, P. M. Naghdi, A dynamical theory of interacting continua, Int. J. Eng. Sci., 3 (1965), 231–241. https://doi.org/10.1016/0020-7225(65)90046-7 doi: 10.1016/0020-7225(65)90046-7
[10]	A. E. Green, P. M. Naghdi, A note on mixtures, Int. J. Eng. Sci., 6 (1968), 631–635. https://doi.org/10.1016/0020-7225(68)90064-5 doi: 10.1016/0020-7225(68)90064-5
[11]	J. D. Ingram, A. C. Eringen, A continuum theory of chemically reacting media Ⅱ, Int. J. Eng. Sci., 5 (1967), 289–322. https://doi.org/10.1016/0020-7225(67)90040-7 doi: 10.1016/0020-7225(67)90040-7
[12]	D. Ieșan, R. Quintanilla, On the theory of interacting continua with memory, J. Therm. Stresses, 25 (2002), 1161–1177. https://doi.org/10.1080/01495730290074586 doi: 10.1080/01495730290074586
[13]	P. D. Kelly, A reacting continuum, Int. J. Eng. Sci., 2 (1964), 129–153. https://doi.org/10.1016/0020-7225(64)90001-1 doi: 10.1016/0020-7225(64)90001-1
[14]	K. R. Rajagopal, L. Tao, Mechanics of mixtures, Ser. Adv. Math. Appl. Sci., 35 (1995). https://doi.org/10.1142/2197 doi: 10.1142/2197
[15]	X. Zhang, E. Zhai, Y. Wu, D. Sun, Theoretical and numerical analyses on Hydro–Thermal–Salt–Mechanical interaction of unsaturated salinized soil subjected to typical unidirectional freezing process, Int. J. Geomech., 21 (2021), 04021104. https://doi.org/10.1061/(ASCE)GM.1943-5622.0002036 doi: 10.1061/(ASCE)GM.1943-5622.0002036
[16]	X. Zhang, Y. Wu, E. Zhai, P. Ye, Coupling analysis of the heat-water dynamics and frozen depth in a seasonally frozen zone, J. Hydrol., 593 (2021), 125603. https://doi.org/10.1016/j.jhydrol.2020.125603 doi: 10.1016/j.jhydrol.2020.125603
[17]	N. Bazarra, J. R. Fernández, R. Quintanilla, Analysis of a Moore-Gibson-Thompson thermoelasticity problem, J. Comput. Appl. Math., 382 (2021), 113058. https://doi.org/10.1016/j.cam.2020.113058 doi: 10.1016/j.cam.2020.113058
[18]	N. Bazarra, J. R. Fernández, R. Quintanilla, On the decay of the energy for radial solutions in Moore-Gibson-Thompson thermoelasticity, Math. Mech. Solids, 26 (2021), 1507–1514. https://doi.org/10.1177/1081286521994258 doi: 10.1177/1081286521994258
[19]	M. Conti, V. Pata, M. Pellicer, R. Quintanilla, A new approach to MGT-thermoviscoelasticity, Discrete Contin. Dyn. Syst., 41 (2021), 4645–4666. https://doi.org/10.3934/dcds.2021052 doi: 10.3934/dcds.2021052
[20]	M. Conti, V. Pata, R. Quintanilla, Thermoelasticity of Moore-Gibson-Thompson type with history dependence in the temperature, Asymptotic Anal., 120 (2020), 1–21. https://doi.org/10.3233/ASY-191576 doi: 10.3233/ASY-191576
[21]	J. R. Fernández, R. Quintanilla, Moore-Gibson-Thompson theory for thermoelastic dielectrics, Appl. Math. Mech., 42 (2021), 309–316. https://doi.org/10.1007/s10483-021-2703-9 doi: 10.1007/s10483-021-2703-9
[22]	K. Jangid, S. Mukhopadhyay, A domain of influence theorem for a natural stress-heat-flux problem in the Moore-Gibson-Thompson thermoelasticity theory, Acta Mech., 232 (2021), 177–187. https://doi.org/10.1007/s00707-020-02833-1 doi: 10.1007/s00707-020-02833-1
[23]	K. Jangid, S. Mukhopadhyay, A domain of influence theorem under MGT thermoelasticity theory, Math. Mech. Solids, 26 (2020), 285–295. https://doi.org/10.1177/1081286520946820 doi: 10.1177/1081286520946820
[24]	R. Quintanilla, Moore-Gibson-Thompson thermoelasticity, Math. Mech. Solids, 24 (2019), 4020–4031. https://doi.org/10.1177/1081286519862007 doi: 10.1177/1081286519862007
[25]	J. R. Fernández, R. Quintanilla, On a mixture of an MGT viscous material and an elastic solid, Acta Mech., 233 (2022), 291–297. https://doi.org/10.1007/s00707-021-03124-z doi: 10.1007/s00707-021-03124-z
[26]	Z. Liu, S. Zheng, Semigroups Associated with Dissipative Systems, Chapman and Hall/CRC, Boca Raton, 1999.
[27]	P. G. Ciarlet, Basic error estimates for elliptic problems, Handb. Numer. Anal., 2 (1993), 17–351. https://doi.org/10.1016/S1570-8659(05)80039-0 doi: 10.1016/S1570-8659(05)80039-0
[28]	M. Campo, J. R. Fernández, K. L. Kuttler, M. Shillor, J. M. Viaño, Numerical analysis and simulations of a dynamic frictionless contact problem with damage, Comput. Methods Appl. Mech. Eng., 196 (2006), 476–488. https://doi.org/10.1016/j.cma.2006.05.006 doi: 10.1016/j.cma.2006.05.006

This article has been cited by:

1.	Shengliang Zhang, A meshless multi-symplectic local radial basis function collocation scheme for the “good” Boussinesq equation, 2022, 431, 00963003, 127297, 10.1016/j.amc.2022.127297
2.	Minghung Lin, Yiyou Hou, Maryam A. Al-Towailb, Hassan Saberi-Nik, The global attractive sets and synchronization of a fractional-order complex dynamical system, 2023, 8, 2473-6988, 3523, 10.3934/math.2023179
3.	Yanyun Xie, Wenliang Cai, Jing Wang, Jesus M. Munoz-Pacheco, Stability and Synchronization of a Fractional‐Order Unified System with Complex Variables, 2024, 2024, 1026-0226, 10.1155/2024/2728661
4.	Shaohui Yan, Hanbing Zhang, Defeng Jiang, Jiawei Jiang, Yu Cui, Yuyan Zhang, Finite-time synchronization of fractional-order chaotic system based on hidden attractors, 2023, 98, 0031-8949, 105226, 10.1088/1402-4896/acf308
5.	Haifeng Huang, Investigation of a high-performance control algorithm for a unified chaotic system synchronization control based on parameter adaptive method, 2024, 18724981, 1, 10.3233/IDT-240178

Reader Comments

Your name:*

Email:*
© 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

Electronic Research Archive

1.1 1.7

Metrics

Article views(1793) PDF downloads(67) Cited by(2)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Figures and Tables

Figures(4) / Tables(1)

Electronic Research Archive

On the mixtures of MGT viscoelastic solids

Related Papers:

Abstract

1. Introduction

2. The method of FCS for fractional-order chaotic systems

3. Applications of FCS method to fractional-order unified chaotic system

4. Conclusions

Acknowledgments

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Catalog

Electronic Research Archive

On the mixtures of MGT viscoelastic solids

Related Papers:

Abstract

1. Introduction

2. The method of FCS for fractional-order chaotic systems

3. Applications of FCS method to fractional-order unified chaotic system

4. Conclusions

Acknowledgments

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog