The Long-Term Dynamic Behavior of Solutions to a Class of Generalized Higher-Order Kirchhoff-Type Coupled Wave Equations

Guoguang Lin, Min Shao

Department of Mathematics, Yunnan University, Kunming, China.

DOI: 10.4236/jamp.2022.107150                

Abstract

In this paper, we study the long-term dynamic behavior of a class of generalized high-order Kirchhoff-type coupled wave equations. Firstly, the existence of uniqueness global solution of this kind of equations in E_k space is proved by prior estimation and Galerkin method; Then, through using Rellich-Kondrachov compact embedding theorem, it is proved that the solution semigroup S(t) has the family of the global attractors A_k in space E_k; Finally, through linearization method, proves that the operator semigroup S(t) Frechet differentiable and the attenuation of linearization problem volume element. Furthermore, we can obtain the finite Hausdorff dimension and Fractal dimension of the family of the global attractors A_k.

Keywords

Kirchhoff Equation, Existence and Uniqueness of Solutions, Global Attractor Family, Dimension Estimation

1. Introduction

In this paper, we study the long-term dynamic behavior of a class of generalized high-order Kirchhoff-type coupled wave equations:

$u_{tt}+M\left({\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{m}v\Vert }^2\right){\left(-\Delta \right)}^{m}u+\beta {\left(-\Delta \right)}^m{u}_t+g\left(u_t,v\right)=f_1\left(x\right),$ (1)

$v_{tt}+M\left({\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{m}v\Vert }^2\right){\left(-\Delta \right)}^{2m}v+\beta {\left(-\Delta \right)}^{2m}{v}_t+g\left(u,v_t\right)=f_2\left(x\right),$ (2)

the boundary conditions:

$\frac{{\partial }^{i}u}{\partial n^i}=0,i=0,1,2,\cdots ,m-1,x\in \partial \Omega ,t>0,$ (3)

$\frac{{\partial }^{j}v}{\partial v^j}=0,j=0,1,2,\cdots ,2m-1,x\in \partial \Omega ,t>0,$ (4)

the initial conditions:

$u\left(x,0\right)=u_0\left(x\right),u_t\left(x,0\right)=u_1\left(x\right),v\left(x,0\right)=v_0\left(x\right),v_t\left(x,0\right)=v_1\left(x\right),x\in \Omega ,$ (5)

where $\Omega$ is a bounded domain in $R^n$ with smooth boundary $\partial \Omega$, $u_0\left(x\right),u_1\left(x\right)$ is a known function, $g\left(u,v\right),f_i\left(x\right),i=1,2$ are nonlinear source term and the external force interference terms, $m>1,\beta$ is real number.

Recently, the global attractor and its dimension estimation for Kirchhoff type equations have been favored by many scholars. Many scholars have done a lot of research on this kind of problems and obtained good results [1] [2] [3].

Lin Guoguang, Gao Yunlong [1] studied the longtime behavior of solution to initial boundary value problem for a class of strongly damped higher-order Kirchhoff type equation:

$u_{tt}+{\left(-\Delta \right)}^m{u}_t+{\left(\alpha +\beta {\Vert {\nabla }^{m}u\Vert }^2\right)}^q{\left(-\Delta \right)}^{m}u+g\left(u\right)=f\left(x\right),\left(x,t\right)\in \Omega \times \left[0,+\infty \right),$

they got the existence and uniqueness of the solution by the Galerkin method and obtained the existence of the global attractor in $H_0^m\left(\Omega \right)\times L^2\left(\Omega \right)$ according to the attractor theorem, besides, the estimation of the upper bound of Hausdorff dimension for the attractor was established.

Guoguang Lin, Ming Zhang [2] studied the initial boundary value problem for a class of Kirchhoff-type coupled equations:

$u_{tt}-M\left({\Vert \nabla u\Vert }^2+{\Vert \nabla v\Vert }^2\right)\Delta u-\beta \Delta u_t+g_1\left(u,v\right)=f_1\left(x\right),$

$v_{tt}-M\left({\Vert \nabla u\Vert }^2+{\Vert \nabla v\Vert }^2\right)\Delta v-\beta \Delta v_t+g_2\left(u,v\right)=f_2\left(x\right),$

they obtained the existence of the global attractor and a precise estimate of upper bound of Hausdorff dimension.

Lin Guoguang, Yang Lujiao [3] studied the long-time properties of solutions of generalized Kirchhoff-type equation with strongly damped terms:

$u_{tt}+M\left({\Vert {\nabla }^{m}u\Vert }_p^p\right){\left(-\Delta \right)}^{2m}u+\beta {\left(-\Delta \right)}^{2m}{u}_t+g\left(u\right)=f\left(x\right),$

by assuming the nonlinear source terms $g\left(u\right)$ and Kirchhoff stress term $M\left(s\right)$, the author verified the appropriateness of the solution and proved the existence of the global attractor, obtained the upper boundary estimation of the Hausdorff dimension and Fractal dimension of a family of the global attractor.

For more significant research results about the global attractor and its dimension estimation of Kirchhoff equation, please refer to the literature [4] - [18].

2. Existence and Uniqueness of Solutions

The following symbols and assumptions are introduced for the convenience of statement:

$\begin{array}{l}{V}_m=H^m\left(\Omega \right)\cap H_0^1\left(\Omega \right),V_{2m}=H^m\left(\Omega \right)\cap H_0^1\left(\Omega \right),V_{4m}=H^{4m}\left(\Omega \right)\cap H_0^1\left(\Omega \right),\\ \Vert \cdot \Vert ={\Vert \cdot \Vert }_{L^2\left(\Omega \right)},E_0=V_m\times V_0\times V_{2m}\times V_0,E_1=V_{2m}\times V_0\times V_{4m}\times V_0,\\ E_k=V_{m+k}\times V_k\times V_{2m+2k}\times V_{2k},V_0=L^2\left(\Omega \right).\end{array}$

In order to obtain our results, we consider system (1)-(5) under some assumptions on $M\left(s\right)$ and $g\left(u,v\right)$. Precisely, we state the general assumptions:

(A1) $M\left(s\right)\in C^2\left(\left[0,+\infty \right),R\right)$ is not decreasing function and for positive constants ${\delta }_0,{\delta }_1$,

(1) $\epsilon +1\le {\delta }_0\le M\left(s\right)\le {\delta }_1$,

(2) $M\left(s\right)$ is a non-negative Lipschitz function, L is associated with the Lipschitz constant $M\left(s\right)$, $M\left(s\right)=M\left({\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{m}v\Vert }^2\right)$.

(A2) For any $u,v,p,q\in V$, $g\left(u_t,v\right),g\left(u,v_t\right)\in C^2\left(ℝ\right)$, there exist $\alpha \ge 0$, $\epsilon \ge 0$, $C\left(\alpha ,\epsilon \right)\ge 0$, such that

$\begin{array}{l}\left(g\left(u_t,v\right),p\right)+\left(g\left(u,v_t\right),q\right)\\ \ge \alpha \left({\Vert p\Vert }^2+{\Vert q\Vert }^2\right)-\epsilon \left({\Vert u\Vert }^2+{\Vert v\Vert }^2\right)-C\left(\alpha ,\epsilon \right)\left({\Vert u_t\Vert }^2+{\Vert v_t\Vert }^2\right).\end{array}$

Lemma 1 Assuming (A1)-(A2) are true, letting $\left(u_0,p_0,v_0,q_0\right)\in E_0$, $f_1\left(x\right),f_2\left(x\right)\in L^2\left(\Omega \right)$, then there is a solution $\left(u,p,v,q\right)$ for problem (1)-(5) which has the following properties:

(i) $\left(u,p,v,q\right)\in L^{\infty }\left(\left(0,+\infty \right);E_0\right)$ ;

(ii)

$y\left(t\right)\le y\left(0\right){\text{e}}^{-\frac{2{\epsilon }_0}{k_1}t}+\frac{k_1}{2\epsilon {\epsilon }_0}\left({\Vert f_1\left(x\right)\Vert }^2+{\Vert f_2\left(x\right)\Vert }^2\right),$ (6)

where $y\left(t\right)={\Vert {\nabla }^{m}u\Vert }^2+{\Vert p\Vert }^2+{\Vert {\nabla }^{2m}v\Vert }^2+{\Vert q\Vert }^2$.

(iii) There are normal numbers $C\left(R_0\right)$ and $t_0=t_0\left(\Omega \right)>0$, such that

${\Vert \left(u,p,v,q\right)\Vert }_{E_0}^2={\Vert {\nabla }^{m}u\Vert }^2+{\Vert p\Vert }^2+{\Vert {\nabla }^{2m}v\Vert }^2+{\Vert q\Vert }^2\le C\left(R_0\right).$ (7)

Proof: Let $p=u_t+\epsilon u$ inner product with Equation (1),

$\left(u_{tt}+M\left({\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{m}v\Vert }^2\right){\left(-\Delta \right)}^{m}u+\beta {\left(-\Delta \right)}^m{u}_t+g\left(u_t,v\right),p\right)=\left(f_1\left(x\right),p\right),$ (8)

according to the hypothesis (A1), using the Young inequality, Holder inequation, Poincare inequality, etc., there are

$\left(u_{tt},p\right)=\frac{1}{2}\frac{\text{d}}{\text{d}t}{\Vert p\Vert }^2-\epsilon {\Vert p\Vert }^2+{\epsilon }^2\left(u,p\right),$ (9)

$\left(M\left({\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{m}v\Vert }^2\right){\left(-\Delta \right)}^{m}u,p\right)\ge \frac{\delta }{2}\frac{\text{d}}{\text{d}t}{\Vert {\nabla }^{m}u\Vert }^2+\epsilon {\delta }_0{\Vert {\nabla }^{m}u\Vert }^2,$ (10)

$\left(\beta {\left(-\Delta \right)}^m{u}_t,p\right)=\beta {\Vert {\nabla }^m{u}_t\Vert }^2+\frac{\beta \epsilon }{2}\frac{\text{d}}{\text{d}t}{\Vert {\nabla }^{m}u\Vert }^2,$ (11)

$\left(f_1\left(x\right),p\right)\le \frac{\epsilon }{2}{\Vert p\Vert }^2+\frac{1}{2\epsilon }{\Vert f_1\left(x\right)\Vert }^2,$ (12)

where $\delta ={\delta }_0$ or ${\delta }_1$.

Similarly, letting $q=v_t+\epsilon v$ inner product with Equation (2), the treatment of each item is similar to (9)-(11), and the above results are sorted out,

$\begin{array}{l}\frac{1}{2}\frac{\text{d}}{\text{d}t}\left(\left(\delta +\beta \epsilon \right)\left({\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{2m}v\Vert }^2\right)+{\Vert p\Vert }^2+{\Vert q\Vert }^2\right)-\frac{3\epsilon }{2}\left({\Vert p\Vert }^2+{\Vert q\Vert }^2\right)\\ +{\epsilon }^2\left[\left(u,p\right)+\left(v,q\right)\right]+\epsilon {\delta }_0\left({\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{2m}v\Vert }^2\right)+\beta \left({\Vert {\nabla }^m{u}_t\Vert }^2+{\Vert {\nabla }^{2m}{v}_t\Vert }^2\right)\\ +\left(g\left(u_t,v\right),p\right)+\left(g\left(u,v_t\right),q\right)=\frac{1}{2\epsilon }\left({\Vert f_1\left(x\right)\Vert }^2+{\Vert f_2\left(x\right)\Vert }^2\right),\end{array}$ (13)

using the Young inequality, the Poincare inequality, and the assumptions (A2), The individual items in Equation (13) are treated as follows:

${\epsilon }^2\left[\left(u,p\right)+\left(v,q\right)\right]\ge -\epsilon \left({\Vert u\Vert }^2+{\Vert v\Vert }^2\right)-\frac{{\epsilon }^3}{4}\left({\Vert p\Vert }^2+{\Vert q\Vert }^2\right),$ (14)

by the Poincare inequality has

$-\epsilon \left({\Vert u\Vert }^2+{\Vert v\Vert }^2\right)\ge -\epsilon \left({\lambda }_1^{-m}+{\lambda }_1^{-2m}\right)\left({\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{2m}v\Vert }^2\right),$ (15)

where ${\lambda }_1$ is the first eigenvalue with homogeneous Dirichlet boundary conditions of $-\Delta$, in the same way

$-C\left(\alpha ,\epsilon \right)\left({\Vert u_t\Vert }^2+{\Vert v_t\Vert }^2\right)\ge -C\left(\alpha ,\epsilon \right)\left({\lambda }_1^{-m}+{\lambda }_1^{-2m}\right)\left({\Vert {\nabla }^m{u}_t\Vert }^2+{\Vert {\nabla }^{2m}{v}_t\Vert }^2\right),$ (16)

$\begin{array}{l}\frac{1}{2}\frac{\text{d}}{\text{d}t}\left(\left(\delta +\beta \epsilon \right)\left({\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{2m}v\Vert }^2\right)+{\Vert p\Vert }^2+{\Vert q\Vert }^2\right)\\ +\left(\alpha -\frac{3\epsilon }{2}-\frac{{\epsilon }^3}{4}\right)\left({\Vert p\Vert }^2+{\Vert q\Vert }^2\right)+\epsilon \left({\delta }_0-2\left({\lambda }_1^{-m}+{\lambda }_1^{-2m}\right)\right)\left({\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{2m}v\Vert }^2\right)\\ +\left(\beta -C\left(\alpha ,\epsilon \right)\left({\lambda }_1^{-m}+{\lambda }_1^{-2m}\right)\right)\left({\Vert {\nabla }^m{u}_t\Vert }^2+{\Vert {\nabla }^{2m}{v}_t\Vert }^2\right),\\ \le \frac{1}{2\epsilon }\left({\Vert f_1\left(x\right)\Vert }^2+{\Vert f_2\left(x\right)\Vert }^2\right),\end{array}$ (17)

where $\alpha >\frac{3\epsilon }{2}+\frac{{\epsilon }^3}{4}$, ${\delta }_0>2\left({\lambda }_1^{-m}+{\lambda }_1^{-2m}\right)$, $\beta >C\left(\alpha ,\epsilon \right)\left({\lambda }_1^{-m}+{\lambda }_1^{-2m}\right)$.

Let $\bar{y}\left(t\right)=\left(\delta +\beta \epsilon \right)\left({\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{2m}v\Vert }^2\right)+{\Vert p\Vert }^2+{\Vert q\Vert }^2$,

$y\left(t\right)={\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{2m}v\Vert }^2+{\Vert p\Vert }^2+{\Vert q\Vert }^2$, there are normal numbers $k_1=\min \left\{1,\delta +\beta \epsilon \right\}$, such that

$\bar{y}\left(t\right)\ge k_{1}y\left(t\right)\ge 0,$ (18)

let ${\epsilon }_0=\min \left\{\alpha -\frac{3\epsilon }{2}-\frac{{\epsilon }^3}{4},\epsilon \left({\delta }_0-2\left({\lambda }_1^{-m}+{\lambda }_1^{-2m}\right)\right)\right\}$, there are

$\frac{\text{d}}{\text{d}t}y\left(t\right)+\frac{2{\epsilon }_0}{k_1}y\left(t\right)\le \frac{1}{\epsilon }\left({\Vert f_1\left(x\right)\Vert }^2+{\Vert f_2\left(x\right)\Vert }^2\right),$ (19)

by Gronwall inequality,

$y\left(t\right)\le y\left(0\right){\text{e}}^{-\frac{2{\epsilon }_0}{k_1}t}+\frac{k_1}{2\epsilon {\epsilon }_0}\left({\Vert f_1\left(x\right)\Vert }^2+{\Vert f_2\left(x\right)\Vert }^2\right),$ (20)

so there are normal numbers $C\left(R_0\right)$ and $t_0=t_0\left(\Omega \right)>0$, such that

${\Vert \left(u,p,v,q\right)\Vert }_{E_0}^2={\Vert {\nabla }^{m}u\Vert }^2+{\Vert p\Vert }^2+{\Vert {\nabla }^{2m}v\Vert }^2+{\Vert q\Vert }^2\le C\left(R_0\right).$ (21)

Lemma 1 is proved.

Lemma 2 Assuming (A1)-(A2) are true, $\left(u_0,p_0,v_0,q_0\right)\in E_k$, $f_1\left(x\right)\in H^m\left(\Omega \right)$, $f_2\left(x\right)\in H^{2m}\left(\Omega \right)$, then there is a solution $\left(u,p,v,q\right)$ for problem (1)-(5), which has the following properties:

(i) $\left(u,p,v,q\right)\in L^{\infty }\left(\left(0,+\infty \right);E_k\right)$ ;

(ii)

$y_1\left(t\right)\le y_1\left(0\right){\text{e}}^{-\frac{2{\epsilon }_1}{k_2}t}+\frac{k_2}{2\epsilon {\epsilon }_1}\left({\Vert {\nabla }^k{f}_1\left(x\right)\Vert }^2+{\Vert {\nabla }^{2k}{f}_2\left(x\right)\Vert }^2\right),$ (22)

where $y\left(t\right)={\Vert {\nabla }^{m}u\Vert }^2+{\Vert p\Vert }^2+{\Vert {\nabla }^{2m}v\Vert }^2+{\Vert q\Vert }^2$.

(iii) There are normal numbers $C\left(R_k\right)$ and $t_{0k}$, such that

${\Vert \left(u,p,v,q\right)\Vert }_{E_k}^2={\Vert {\nabla }^{m+k}u\Vert }^2+{\Vert {\nabla }^{2m+2k}v\Vert }^2+{\Vert {\nabla }^{k}p\Vert }^2+{\Vert {\nabla }^{2k}q\Vert }^2\le C\left(R_k\right).$ (23)

Proof: Let ${\left(-\Delta \right)}^{k}p$ inner product with Equation (1),

$\begin{array}{l}\left(u_{tt}+M\left({\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{m}v\Vert }^2\right){\left(-\Delta \right)}^{m}u+\beta {\left(-\Delta \right)}^m{u}_t+g\left(u_t,v\right),{\left(-\Delta \right)}^{k}p\right)\\ =\left(f_1\left(x\right),{\left(-\Delta \right)}^{k}p\right),\end{array}$ (24)

according to the hypothesis (A1), using the Young inequality, Holder inequation, Poincare inequality, etc., there are

$\left(u_{tt},{\left(-\Delta \right)}^{k}p\right)\ge \frac{1}{2}\frac{\text{d}}{\text{d}t}{\Vert {\nabla }^{k}p\Vert }^2-\left(\epsilon +\frac{{\epsilon }^2}{2}\right){\Vert {\nabla }^{k}p\Vert }^2-\frac{{\epsilon }^2{\lambda }_1^{-m}}{2}{\Vert {\nabla }^{m+k}u\Vert }^2,$ (25)

$\left(M\left({\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{m}v\Vert }^2\right){\left(-\Delta \right)}^{m}u,{\left(-\Delta \right)}^{k}p\right)\ge \frac{\delta }{2}\frac{\text{d}}{\text{d}t}{\Vert {\nabla }^{m+k}u\Vert }^2+\epsilon {\delta }_0{\Vert {\nabla }^{m+k}u\Vert }^2,$ (26)

$\begin{array}{l}\left(\beta {\left(-\Delta \right)}^m{u}_t,{\left(-\Delta \right)}^{k}p\right)\\ =\beta {\Vert {\nabla }^{m+k}{u}_t\Vert }^2+\frac{\beta \epsilon }{2}\frac{\text{d}}{\text{d}t}{\Vert {\nabla }^{m+k}u\Vert }^2,\end{array}$ (27)

where $\delta ={\delta }_0$ or ${\delta }_1$.

Similarly, letting ${\left(-\Delta \right)}^{2m}q$ inner product with Equation (2), the treatment of each item is similar to (25)-(27), using Young inequality, Poincare inequality, and assumption (A2), the individual terms are treated as follows:

$\begin{array}{l}\left(g\left(u_t,v\right),{\left(-\Delta \right)}^{k}p\right)+\left(g\left(u,v_t\right),{\left(-\Delta \right)}^{2k}q\right)\\ \le \frac{C_1}{2\epsilon }\left({\lambda }_1^{-m}+{\lambda }_1^{-2m}\right)\left({\Vert {\nabla }^{m+k}{u}_t\Vert }^2+{\Vert {\nabla }^{2m+2k}{v}_t\Vert }^2\right)+\frac{\epsilon C_1}{2}\left({\Vert {\nabla }^{k}p\Vert }^2+{\Vert {\nabla }^{2k}q\Vert }^2\right),\end{array}$ (28)

$\begin{array}{l}\left(f_1\left(x\right),{\left(-\Delta \right)}^{k}p\right)+\left(f_2\left(x\right),{\left(-\Delta \right)}^{2k}2q\right)\\ \le \frac{\epsilon }{2}\left({\Vert {\nabla }^{k}p\Vert }^2+{\Vert {\nabla }^{2k}q\Vert }^2\right)+\frac{1}{2\epsilon }\left({\Vert {\nabla }^k{f}_1\left(x\right)\Vert }^2+{\Vert {\nabla }^{2k}{f}_2\left(x\right)\Vert }^2\right).\end{array}$ (29)

Then sort out the result of appeal and other inner product items, get

$\begin{array}{l}\frac{1}{2}\frac{\text{d}}{\text{d}t}\left(\left(\delta +\beta \epsilon \right)\left({\Vert {\nabla }^{m+k}u\Vert }^2+{\Vert {\nabla }^{2m+2k}v\Vert }^2\right)+{\Vert {\nabla }^{k}p\Vert }^2+{\Vert {\nabla }^{2k}q\Vert }^2\right)\\ +\frac{{\epsilon }^2-\epsilon C_1-3\epsilon }{2}\left({\Vert {\nabla }^{k}p\Vert }^2+{\Vert {\nabla }^{2k}q\Vert }^2\right)\\ +\left(\epsilon {\delta }_0-\frac{{\epsilon }^2}{2}\left({\lambda }_1^{-m}+{\lambda }_1^{-2m}\right)\right)\left({\Vert {\nabla }^{m+k}u\Vert }^2+{\Vert {\nabla }^{2m+2k}v\Vert }^2\right)\\ +\left(\beta -\frac{C_1}{2\epsilon }\left({\lambda }_1^{-m}+{\lambda }_1^{-2m}\right)\right)\left({\Vert {\nabla }^{m+k}{u}_t\Vert }^2+{\Vert {\nabla }^{2m+2k}{v}_t\Vert }^2\right)\\ \le \frac{1}{2\epsilon }\left({\Vert {\nabla }^k{f}_1\left(x\right)\Vert }^2+{\Vert {\nabla }^{2k}{f}_2\left(x\right)\Vert }^2\right),\end{array}$ (30)

using the Young inequality, Poincare inequality, and assumption (A2), the individual terms in Equation (30) are treated as follows:

$\begin{array}{l}\left(g\left(u_t,v\right),{\left(-\Delta \right)}^{k}p\right)+\left(g\left(u,v_t\right),{\left(-\Delta \right)}^{2k}q\right)\\ \le \frac{C_1}{2\epsilon }\left({\lambda }_1^{-m}+{\lambda }_1^{-2m}\right)\left({\Vert {\nabla }^{m+k}{u}_t\Vert }^2+{\Vert {\nabla }^{2m+2k}{v}_t\Vert }^2\right)+\frac{\epsilon C_1}{2}\left({\Vert {\nabla }^{k}p\Vert }^2+{\Vert {\nabla }^{2k}q\Vert }^2\right),\end{array}$ (31)

$\begin{array}{l}\left(f_1\left(x\right),{\left(-\Delta \right)}^{k}p\right)+\left(f_2\left(x\right),{\left(-\Delta \right)}^{2k}2q\right)\\ \le \frac{\epsilon }{2}\left({\Vert {\nabla }^{k}p\Vert }^2+{\Vert {\nabla }^{2k}q\Vert }^2\right)+\frac{1}{2\epsilon }\left({\Vert {\nabla }^k{f}_1\left(x\right)\Vert }^2+{\Vert {\nabla }^{2k}{f}_2\left(x\right)\Vert }^2\right),\end{array}$ (32)

where $\epsilon -C_1-3>0$, ${\delta }_0>\frac{\epsilon }{2}\left({\lambda }_1^{-m}+{\lambda }_1^{-2m}\right)$, $\beta >\frac{C_1}{2\epsilon }\left({\lambda }_1^{-m}+{\lambda }_1^{-2m}\right)$.

Let ${\bar{y}}_1\left(t\right)=\left(\delta +\beta \epsilon \right)\left({\Vert {\nabla }^{m+k}u\Vert }^2+{\Vert {\nabla }^{2m+2k}v\Vert }^2\right)+{\Vert {\nabla }^{k}p\Vert }^2+{\Vert {\nabla }^{2k}q\Vert }^2$,

$y_1\left(t\right)={\Vert {\nabla }^{m+k}u\Vert }^2+{\Vert {\nabla }^{2m+2k}v\Vert }^2+{\Vert {\nabla }^{k}p\Vert }^2+{\Vert {\nabla }^{2k}q\Vert }^2$, there are normal numbers $k_2=\min \left\{1,\delta +\beta \epsilon \right\}$, such that

${\bar{y}}_1\left(t\right)\ge k_2{y}_1\left(t\right)\ge 0,$ (33)

let ${\epsilon }_1=\min \left\{\frac{{\epsilon }^2-\epsilon C_1-3\epsilon }{2}-\frac{{\epsilon }^3}{4},\epsilon {\delta }_0-\frac{{\epsilon }^2}{2}\left({\lambda }_1^{-m}+{\lambda }_1^{-2m}\right)\right\}$, get

$\frac{\text{d}}{\text{d}t}{y}_1\left(t\right)+\frac{2{\epsilon }_1}{k_2}{y}_1\left(t\right)\le \frac{1}{\epsilon }\left({\Vert {\nabla }^k{f}_1\left(x\right)\Vert }^2+{\Vert {\nabla }^{2k}{f}_2\left(x\right)\Vert }^2\right),$ (34)

by Gronwall inequality,

$y_1\left(t\right)\le y_1\left(0\right){\text{e}}^{-\frac{2{\epsilon }_1}{k_2}t}+\frac{k_2}{2\epsilon {\epsilon }_1}\left({\Vert {\nabla }^k{f}_1\left(x\right)\Vert }^2+{\Vert {\nabla }^{2k}{f}_2\left(x\right)\Vert }^2\right),$ (35)

so there are normal numbers $C\left(R_k\right)$ and $t_{0k}>0$, such that

${\Vert \left(u,p,v,q\right)\Vert }_{E_k}^2={\Vert {\nabla }^{m+k}u\Vert }^2+{\Vert {\nabla }^{2m+2k}v\Vert }^2+{\Vert {\nabla }^{k}p\Vert }^2+{\Vert {\nabla }^{2k}q\Vert }^2\le C\left(R_k\right)$ (36)

Lemma 1 is proved.

Theorem 1 Assuming (A1)-(A2) is true, $\left(u_0,p_0,v_0,q_0\right)\in E_k$, $f_1\left(x\right)\in V_k$, $f_2\left(x\right)\in V_{2k}$, then the initial boundary value problem (1)-(5) has a unique solution

$\left(u\left(x,t\right),p\left(x,t\right),v\left(x,t\right),q\left(x,t\right)\right)\in L^{\infty }\left(\left(0,+\infty \right);E_k\right).$ (37)

Proof: According to literature [9] and Galerkin method, combining with lemma 1 and lemma 2, we can easily obtain the existence of solutions.

Next, prove the uniqueness of the solution:

Assuming $\left(u_1,p_1,v_1,q_1\right),\left(u_2,p_2,v_2,q_2\right)\in E_k$ are the two solutions of the problem (1)-(5), letting $\bar{u}=u_1-u_2,\bar{v}=v_1-v_2$, obtain that

$\{\begin{array}{l}{\bar{u}}_{tt}+M\left(s_1\right){\left(-\Delta \right)}^m{u}_1-M\left(s_2\right){\left(-\Delta \right)}^m{u}_2+\beta {\left(-\Delta \right)}^m{\bar{u}}_t+g\left(u_{1t},v_1\right)-g\left(u_{2t},v_2\right)=0,\\ {\bar{v}}_{tt}+M\left(s_1\right){\left(-\Delta \right)}^{2m}{v}_1-M\left(s_2\right){\left(-\Delta \right)}^{2m}{v}_2+\beta {\left(-\Delta \right)}^{2m}{\bar{v}}_t+g\left(u_1,v_{1t}\right)-g\left(u_2,v_{2t}\right)=0,\\ \bar{u}\left(x,0\right)=0,{\bar{u}}_t\left(x,0\right)=0,\bar{v}\left(x,0\right)=0,{\bar{v}}_t\left(x,0\right)=0,x\in \Omega ,\\ \frac{{\partial }^i\bar{u}}{\partial n^i}=0,i=0,1,2,\cdots ,m-1,x\in \partial \Omega ,t>0,\\ \frac{{\partial }^j\bar{v}}{\partial v^j}=0,j=0,1,2,\cdots ,2m-1,x\in \partial \Omega ,t>0,\end{array}$(38)

where $s_1={\Vert {\nabla }^m{u}_1\Vert }^2-{\Vert {\nabla }^m{v}_1\Vert }^2,s_2={\Vert {\nabla }^m{u}_2\Vert }^2-{\Vert {\nabla }^m{v}_2\Vert }^2$.

Let ${\bar{u}}_t,{\bar{v}}_t$ inner product with the first two equations in (38) and obtain,

$\{\begin{array}{l}\frac{1}{2}\frac{\text{d}}{\text{d}t}{\Vert {\bar{u}}_t\Vert }^2+\left(M\left(s_1\right){\left(-\Delta \right)}^m{u}_1-M\left(s_2\right){\left(-\Delta \right)}^m{u}_2,{\bar{u}}_t\right)+\beta {\Vert {\nabla }^m{\bar{u}}_t\Vert }^2=\left(g\left(u_{2t},v_2\right)-g\left(u_{1t},v_1\right),{\bar{u}}_t\right),\hfill \\ \frac{1}{2}\frac{\text{d}}{\text{d}t}{\Vert {\bar{v}}_t\Vert }^2+\left(M\left(s_1\right){\left(-\Delta \right)}^{2m}{v}_1-M\left(s_2\right){\left(-\Delta \right)}^{2m}{v}_2,{\bar{v}}_t\right)+\beta {\Vert {\nabla }^{2m}{\bar{v}}_t\Vert }^2=\left(g\left(u_2,v_{2t}\right)-g\left(u_1,v_{1t}\right),{\bar{v}}_t\right),\hfill \end{array}$(39)

using the Young inequality, Poincare inequality, as well as the assumptions (A2), for processing on the type of individual items as follows:

$\begin{array}{l}\left(M\left(s_1\right){\left(-\Delta \right)}^m\bar{u}+M\left(s_1\right){\left(-\Delta \right)}^m{u}_2-M\left(s_2\right){\left(-\Delta \right)}^m{u}_2,{\bar{u}}_t\right)\\ \ge \frac{M\left(s\right)}{2}\frac{\text{d}}{\text{d}t}{\Vert {\nabla }^{m}u\Vert }^2-\left(M\left(s_1\right){\left(-\Delta \right)}^m{u}_2-M\left(s_2\right){\left(-\Delta \right)}^m{u}_2,{\bar{u}}_t\right)\\ \ge \frac{M\left(s\right)}{2}\frac{\text{d}}{\text{d}t}{\Vert {\nabla }^{m}u\Vert }^2-L(\left(\Vert {\nabla }^m{u}_1\Vert +\Vert {\nabla }^m{u}_2\Vert \right)\Vert {\nabla }^m\bar{u}\Vert \\ +L\left(\Vert {\nabla }^m{v}_1\Vert +\Vert {\nabla }^m{v}_2\Vert \right)\Vert {\nabla }^m\bar{v}\Vert )\Vert {\left(-\Delta \right)}^m{\bar{u}}_2\Vert \Vert {\bar{u}}_t\Vert \\ \ge C_2\left(\Vert {\nabla }^m\bar{u}\Vert +\Vert {\nabla }^{2m}\bar{v}\Vert \right)\Vert {\bar{u}}_t\Vert \\ \ge C_2\left({\Vert {\nabla }^m\bar{u}\Vert }^2+{\Vert {\nabla }^{2m}\bar{v}\Vert }^2+\frac{{\Vert {\bar{u}}_t\Vert }^2}{2}\right),\end{array}$ (40)

similarly, we obtain

$\left(M\left(s_1\right){\left(-\Delta \right)}^{2m}{v}_1-M\left(s_2\right){\left(-\Delta \right)}^{2m}{v}_2,{\bar{v}}_t\right)\ge C_3\left({\Vert {\nabla }^m\bar{u}\Vert }^2+{\Vert {\nabla }^{2m}\bar{v}\Vert }^2+\frac{{\Vert {\bar{v}}_t\Vert }^2}{2}\right),$ (41)

$\begin{array}{l}\left|\left(g\left(u_{1t},v_1\right)-g\left(u_{2t},v_1\right)+g\left(u_{2t},v_1\right)-g\left(u_{2t},v_2\right),{\bar{u}}_t\right)\right|\\ \le \left|\left(g^{\prime }\left({\xi }_t,v_1\right)\left|{\bar{u}}_t\right|+g^{\prime }\left(u_{2t},\eta \right)\left|\bar{v}\right|,{\bar{u}}_t\right)\right|\\ \le {\Vert g^{\prime }\left({\xi }_t,v_1\right)\Vert }_{\infty }{\Vert {\bar{u}}_t\Vert }^2+\Vert g^{\prime }{\left(u_{2t},\eta \right)}_{\infty }\Vert \Vert \left|\bar{v}\right|\Vert \Vert {\bar{u}}_t\Vert \le C_4\left(\frac{{\Vert \bar{v}\Vert }^2}{2}+\frac{3{\Vert {\bar{u}}_t\Vert }^2}{2}\right),\end{array}$ (42)

$\left|\left(g\left(u_1,v_{1t}\right)-g\left(u_2,v_{2t}\right),{\bar{v}}_t\right)\right|\le C_5\left(\frac{3{\Vert {\bar{v}}_t\Vert }^2}{2}+\frac{{\Vert \bar{u}\Vert }^2}{2}\right),$ (43)

Through the (40)-(43), finally will become

$\begin{array}{l}\frac{1}{2}\frac{\text{d}}{\text{d}t}\left({\Vert {\bar{u}}_t\Vert }^2+{\Vert {\bar{v}}_t\Vert }^2+M\left(s_1\right)\left({\Vert {\nabla }^m\bar{u}\Vert }^2+{\Vert {\nabla }^{2m}\bar{v}\Vert }^2\right)\right)+\beta \left({\Vert {\nabla }^m{\bar{u}}_t\Vert }^2+{\Vert {\nabla }^{2m}{\bar{v}}_t\Vert }^2\right)\\ \le \left(C_2+C_3\right)\left({\Vert {\nabla }^m\bar{u}\Vert }^2+{\Vert {\nabla }^{2m}\bar{v}\Vert }^2\right)+\frac{C_2+3C_4}{2}{\Vert {\bar{u}}_t\Vert }^2\\ +\frac{C_3+3C_5}{2}{\Vert {\bar{v}}_t\Vert }^2+\frac{C_4}{2}{\Vert v\Vert }^2+\frac{C_5}{2}{\Vert u\Vert }^2\\ \le C_6\left({\Vert {\nabla }^m\bar{u}\Vert }^2+{\Vert {\nabla }^{2m}\bar{v}\Vert }^2\right)+C_7\left({\Vert {\bar{u}}_t\Vert }^2+{\Vert {\bar{v}}_t\Vert }^2\right).\end{array}$ (44)

Let $y_2\left(t\right)={\Vert {\bar{u}}_t\Vert }^2+{\Vert {\bar{v}}_t\Vert }^2+M\left(s_1\right)\left({\Vert {\nabla }^m\bar{u}\Vert }^2+{\Vert {\nabla }^{2m}\bar{v}\Vert }^2\right)$, there are normal numbers $k_3=\min \left\{C_7,C_8\right\}$, where $C_6\le C_{8}M\left(s_1\right)$, such that

$\frac{\text{d}}{\text{d}t}{y}_2\left(t\right)\le k_3{y}_2\left(t\right),$ (45)

using the Gronwall inequality, we have

$y_2\left(t\right)\le y_2\left(0\right){\text{e}}^{k_{4}t}=0,$ (46)

$y_2\left(t\right)={\Vert {\bar{u}}_t\Vert }^2+{\Vert {\bar{v}}_t\Vert }^2+M\left(s_1\right)\left({\Vert {\nabla }^m\bar{u}\Vert }^2+{\Vert {\nabla }^{2m}\bar{v}\Vert }^2\right)=0,$ (47)

so

$\bar{u}=\bar{v}=0.$ (48)

Theorem 1 is proved.

3. The Family of Global Attractors and Dimension Estimation

Theorem 2 [9] Assume $E_0$ is a Banach space, ${\left\{S\left(t\right)\right\}}_{t\ge 0}$ is the operator semigroup, $S\left(t\right):E\to E$, $S\left(t+\tau \right)=S\left(t\right)\cdot S\left(\tau \right)\left(\forall t,\tau \ge 0\right)$, $S\left(0\right)=I$, where I is the identity operator. If $S\left(t\right)$ satisfies

1) Semigroup $S\left(t\right)$ is uniformly bounded in E, i.e. $\forall R>0$, exists a constant $C\left(R\right)$ such that when ${\Vert u\Vert }_E\le R$, there is ${\Vert S\left(t\right)u\Vert }_E\le C\left(R\right)\left(\forall t\in \left[0,+\infty \right)\right)$ ;

2) There exists a bounded absorbing set B in E, that is, for any bounded set $B\subset E$, there exists a constant $t_0>0$, such that

$S\left(t\right)B\subset B\left(t>t_0\right)$, (49)

3) ${\left\{S\left(t\right)\right\}}_{t\ge 0}$ is completely continuous operator.

Then operator semigroup $S\left(t\right)$ has compact global attractor A.

Theorem 3 Let $S\left(t\right)$ is a solution semigroup generated by the initial boundary value problems (1)-(5) under the hypothesis of lemma 1 and lemma 2, then the initial boundary value problems (1)-(5) have the family of global attractors. There are compact sets satisfying:

$A_k\subset E_k\subset E_0$ and $A_k=\omega \left(B_{0k}\right)=\underset{\tau \ge 0}{\cap }\overline{\underset{t\ge \tau }{\cup }S\left(t\right)B_{0k}},k=1,2,\cdots ,m.$

where $B_{0k}=\left\{\left(u,p,v,q\right)\in E_k:{\Vert {\nabla }^{m+k}u\Vert }^2+{\Vert {\nabla }^{k}p\Vert }^2+{\Vert {\nabla }^{2m+2k}v\Vert }^2+{\Vert {\nabla }^{2k}q\Vert }^2\le R_{0k}\right\}$,

1) $S\left(t\right)A_k=A_k,t>0$,

2) $A_k$ attracts all bounded sets of $E_k$, that is, any bounded set $B_{0k}\subset E_k$, having

$\underset{t\to \infty }{\lim }dist\left(S\left(t\right)B_k,A_k\right)=0$, where $dist\left(S\left(t\right)B_k,A_k\right)=\underset{x\in B_{0k}}{\sup }\underset{y\in A_k}{\inf }{\Vert S\left(t\right)-y\Vert }_{E_k}$.

then compact set $A_k$ are called the family of global attractors of semigroup $S\left(t\right)$.

Proof: From lemma 1, lemma 2, for any bounded set $B_{0k}\subset E_k$ and $B_{0k}\subset \left\{{\Vert \left(u,p,v,q\right)\Vert }_{E_k}\le R_k\right\}$, the equation has solution semigroups $S\left(t\right):E_k\to E_k$, and

${\Vert S\left(t\right)\left(u_0,p_0,v_0,q_0\right)\Vert }_{E_k}^2={\Vert u\Vert }_{V_{m+k}}^2+{\Vert p\Vert }_{V_k}^2+{\Vert v\Vert }_{V_{2m+2k}}^2+{\Vert q\Vert }_{V_{2k}}^2\le C\left(R_k\right),$ (50)

where $t_0\ge 0$, $\left(u_0,v_0\right)\in B_{0k}$, shows that $S{\left(t\right)}_{t\ge 0}$ is uniformly bounded in $E_k$ ;

Further,

$B_{0k}=\left\{\left(u,p,v,q\right)\in E_k:{\Vert {\nabla }^{m+k}u\Vert }^2+{\Vert {\nabla }^{k}p\Vert }^2+{\Vert {\nabla }^{2m+2k}v\Vert }^2+{\Vert {\nabla }^{2k}q\Vert }^2\le R_{0k}\right\}$ is a

bounded absorption set of semigroup $S\left(t\right)$ ; $E_k$ is compactly embedded in $E_0$, i.e., the bounded set in $E_k$ is a compact set in $E_0$, so the operator semigroup $S\left(t\right)$ is completely continuous operator. Then there exists a global attractor family of equations

$A_k=\omega \left(B_{0k}\right)=\underset{\tau \ge 0}{\cap }\overline{\underset{t\ge \tau }{\cup }S\left(t\right)B_{0k}},k=1,2,\cdots ,m.$

Theorem 3 is proved.

After the family of global attractors are obtained, in order to estimate the Hausdroff dimension and Fractal dimension of the family of global attractors, the initial boundary value problem (1)-(5) is linearized and obtain that

$\{\begin{array}{l}{U}_{tt}+M^{\prime }\left({\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{m}v\Vert }^2\right)\left[\left({\nabla }^{m}u,{\nabla }^{m}U\right)+\left({\nabla }^{m}v,{\nabla }^{m}V\right)\right]{\left(-\Delta \right)}^{m}u\\ +M\left({\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{m}v\Vert }^2\right){\left(-\Delta \right)}^{m}U+\beta {\left(-\Delta \right)}^m{U}_t+\frac{\partial g\left(u_t,v\right)}{\partial u_t}{U}_t+\frac{\partial g\left(u_t,v\right)}{\partial v}V=0,\\ U_{tt}+M^{\prime }\left({\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{m}v\Vert }^2\right)\left[\left({\nabla }^{m}u,{\nabla }^{m}U\right)+\left({\nabla }^{m}v,{\nabla }^{m}V\right)\right]{\left(-\Delta \right)}^{m}u\\ +M\left({\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{m}v\Vert }^2\right){\left(-\Delta \right)}^{m}U+\beta {\left(-\Delta \right)}^m{U}_t+\frac{\partial g\left(u_t,v\right)}{\partial u_t}{U}_t+\frac{\partial g\left(u_t,v\right)}{\partial v}V=0,\\ U\left(x,0\right)={\xi }_1,U_t\left(x,0\right)={\xi }_2,V\left(x,0\right)={\eta }_1,V_t\left(x,0\right)={\eta }_2,\\ {U\left(x,0\right)|}_{x\in \Omega }={V\left(x,0\right)|}_{x\in \Omega }=0,t>0.\end{array}$(51)

where $\left({\xi }_1,{\xi }_2,{\eta }_1,{\eta }_2\right)\in E_0$, $\left(u,p,v,q\right)=S\left(t\right)\left(u_0,p_0,v_0,q_0\right)$ is the solution of the initial boundary value problem (51).

Given $\left(u_0,p_0,v_0,q_0\right)\in A_k$, $S\left(t\right):E_k\to E_k$, for any $\left({\xi }_1,{\xi }_2,{\eta }_1,{\eta }_2\right)\in E_k$, there exists a unique solution $\left(U\left(t\right),P\left(t\right),V\left(t\right),Q\left(t\right)\right)\in L^{\infty }\left(0,T;E_k\right)$ to the linear initial boundary value problem (51).

Lemma 3 For any $t>0$, $R>0$, the mapping $S\left(t\right):E_k\to E_k$ is Frechet differentiable. The derivative on ${\rho }_0={\left(u_0,p_0,v_0,q_0\right)}^{\text{T}}$ is a linear operator on $E_k$,

$DS\left(t\right){\rho }_0:\left(\xi ,\zeta ,\eta ,\sigma \right)\to \left(U,P,V,Q\right),$

where $\left(U\left(t\right),P\left(t\right),V\left(t\right),Q\left(t\right)\right)$ is the solution of the problem (51).

Proof: suppose ${\rho }_0=\left(u_0,p_0,v_0,q_0\right)\in E_k$, ${\tilde{\rho }}_0=\left(u_0+\xi ,p_0+\zeta ,v_0+\eta ,q_0+\sigma \right)\in E_k$ with ${\Vert {\rho }_0\Vert }_{E_k}\le R_k$, ${\Vert {\tilde{\rho }}_0\Vert }_{E_k}\le R_k$, we denote

$S\left(t\right){\rho }_0=\rho =\left(u_1,p_1,v_1,q_1\right),S\left(t\right){\tilde{\rho }}_0=\left(u_2,p_2,v_2,q_2\right).$

First, we can prove a Lipschitz property of $S\left(t\right)$ on the bounded sets on $E_k$, that is

${\Vert S\left(t\right){\rho }_0-S\left(t\right){\tilde{\rho }}_0\Vert }_{E_k}^2\le {\text{e}}^{ct}{\Vert \left(\xi ,\zeta ,\eta ,\sigma \right)\Vert }_{E_k}^2.$ (52)

We now consider the difference $\theta =u_2-u_1-U,\omega =v_2-v_1-V$ is the solution to problem (53),

$\{\begin{array}{l}{\theta }_{tt}+M\left({\Vert {\nabla }^m{u}_1\Vert }^2+{\Vert {\nabla }^m{v}_1\Vert }^2\right){\left(-\Delta \right)}^m\theta +\beta {\left(-\Delta \right)}^m{\theta }_t=h_1,\\ {\omega }_{tt}+M\left({\Vert {\nabla }^m{u}_1\Vert }^2+{\Vert {\nabla }^m{v}_1\Vert }^2\right){\left(-\Delta \right)}^{2m}\omega +\beta {\left(-\Delta \right)}^{2m}{\omega }_t=h_2,\\ \omega \left(0\right)={\omega }_t\left(0\right)=0,\theta \left(0\right)={\theta }_t\left(0\right)=0.\end{array}$ (53)

where

$\begin{array}{c}{h}_1=\left[M\left(s\right)-M\left(\tilde{s}\right)\right]{\left(-\Delta \right)}^m{u}_2\\ +2M^{\prime }\left(s\right)\left[\left({\nabla }^m{u}_1,{\nabla }^{m}U\right)+\left({\nabla }^m{v}_1,{\nabla }^{m}V\right)\right]{\left(-\Delta \right)}^m{u}_1\\ +\frac{\partial g\left(u_{1t},v_1\right)}{\partial u_{1t}}{U}_t+\frac{\partial g\left(u_{1t},v_1\right)}{\partial v_1}V+g\left(u_{1t},v_1\right)-g\left(u_{2t},v_2\right),\end{array}$ (54)

$\begin{array}{c}{h}_2=\left[M\left(s\right)-M\left(\tilde{s}\right)\right]{\left(-\Delta \right)}^{2m}{v}_2\\ +2M^{\prime }\left(s\right)\left[\left({\nabla }^m{u}_1,{\nabla }^{m}U\right)+\left({\nabla }^m{v}_1,{\nabla }^{m}V\right)\right]{\left(-\Delta \right)}^{2m}{v}_1\\ +\frac{\partial g\left(u_1,v_{1t}\right)}{\partial u_1}U+\frac{\partial g\left(u_1,v_{1t}\right)}{\partial v_{1t}}{V}_t+g\left(u_1,v_{1t}\right)-g\left(u_2,v_{2t}\right),\end{array}$ (55)

where $s={\Vert {\nabla }^m{u}_1\Vert }^2+{\Vert {\nabla }^m{v}_1\Vert }^2,\tilde{s}={\Vert {\nabla }^m{u}_2\Vert }^2+{\Vert {\nabla }^m{v}_2\Vert }^2$.

Some items are treated as follows

$\begin{array}{l}\left[M\left(s\right)-M\left(\tilde{s}\right)\right]{\left(-\Delta \right)}^m{u}_2+2M^{\prime }\left(s\right)\left[\left({\nabla }^{m}U,{\nabla }^m{u}_1\right)+\left({\nabla }^{m}V,{\nabla }^m{v}_1\right)\right]{\left(-\Delta \right)}^m{u}_1\\ =\left(M\left(s\right)-M\left(\tilde{s}\right)\right){\left(-\Delta \right)}^m{u}_2+2M^{\prime }\left(s\right)\left[\left({\nabla }^m{u}_1,{\nabla }^{m}U\right)+\left({\nabla }^m{v}_1,{\nabla }^{m}V\right)\right]{\left(-\Delta \right)}^m{u}_1\\ =-2M^{\prime }\left(s\right)\left[\left(-{\nabla }^{m}U,{\nabla }^m{u}_1\right)+\left(-{\nabla }^{m}V,{\nabla }^m{v}_1\right)\right]{\left(-\Delta \right)}^m{u}_1+M^{\prime }\left(\alpha s+\left(1-\alpha \right)\tilde{s}\right)\\ \times \left[\left({\nabla }^m\left(u_1-u_2\right),{\nabla }^m\left(u_2+u_1\right)\right)+\left({\nabla }^m\left(v_1-v_2\right),{\nabla }^m\left(v_2+v_1\right)\right)\right]{\left(-\Delta \right)}^m{u}_2\end{array}$

$\begin{array}{l}\le M^{″}\left(\xi \right){\left[\left({\nabla }^m\left(u_1-u_2\right),{\nabla }^m\left(u_2+u_1\right)\right)+\left({\nabla }^m\left(v_1-v_2\right),{\nabla }^m\left(v_2+v_1\right)\right)\right]}^2{\left(-\Delta \right)}^m{u}_2\\ +M^{\prime }\left(s\right)\left[\left({\nabla }^m\left(u_1-u_2\right),{\nabla }^m\left(u_2-u_1\right)\right)+\left({\nabla }^m\left(v_1-v_2\right),{\nabla }^m\left(v_2-v_1\right)\right)\right]{\left(-\Delta \right)}^m{u}_2\\ +2M^{\prime }\left(s\right)\left[\left({\nabla }^m\left(u_1-u_2\right),{\nabla }^m{u}_1\right)+\left({\nabla }^m\left(v_1-v_2\right),{\nabla }^m{v}_1\right)\right]{\left(-\Delta \right)}^m\left(u_2-u_1\right)\\ -2M^{\prime }\left(s\right)\left[\left({\nabla }^m\theta ,{\nabla }^m{u}_1\right)+\left({\nabla }^m\omega ,{\nabla }^m{v}_1\right)\right]{\left(-\Delta \right)}^m{u}_1\\ =I_1+I_2+I_3+I_4.\end{array}$(56)

Let ${\left(-\Delta \right)}^k{\theta }_t$ inner product with the first equation in (53), ${\left(-\Delta \right)}^{2k}{\omega }_t$ and inner product with the second equation in (53), obtain

$\begin{array}{c}\left(I_1,{\left(-\Delta \right)}^k{\theta }_t\right)\le \left(2M^{″}\left(\xi \right)\left[{\left({\nabla }^m\left(u_1-u_2\right),{\nabla }^m\left(u_2+u_1\right)\right)}^2+{\left({\nabla }^m\left(v_1-v_2\right),{\nabla }^m\left(v_2+v_1\right)\right)}^2\right]{\nabla }^{m+k}{u}_2,{\nabla }^{m+k}{\theta }_t\right)\\ \le 2c_1\left(\frac{4}{{\epsilon }^2}{\Vert {\nabla }^m\left(u_2-u_1\right)\Vert }^4+\frac{4}{{\epsilon }^2}{\Vert {\nabla }^m\left(v_2-v_1\right)\Vert }^4+\frac{{\epsilon }^2{\Vert {\nabla }^{m+k}{\theta }_t\Vert }^2}{8}\right),\\ \le 2c_1\left(\frac{4}{{\epsilon }^2}\left({\lambda }_1^{-2k}+{\lambda }_1^{-2m-4k}\right)\left({\Vert {\nabla }^{m+k}\left(u_2-u_1\right)\Vert }^4+{\Vert {\nabla }^{2m+2k}\left(v_2-v_1\right)\Vert }^4\right)+\frac{{\epsilon }^2{\Vert {\nabla }^{m+k}{\theta }_t\Vert }^2}{8}\right),\end{array}$

$\begin{array}{c}\left(I_2,{\left(-\Delta \right)}^k{\theta }_t\right)=\left(M^{\prime }\left(s\right)\left[\left({\nabla }^m\left(u_1-u_2\right),{\nabla }^m\left(u_2-u_1\right)\right)+\left({\nabla }^m\left(v_1-v_2\right),{\nabla }^m\left(v_2-v_1\right)\right)\right]{\left(-\Delta \right)}^{m+k}{u}_2,{\nabla }^{m+k}{\theta }_t\right)\\ \le c_2\left(\frac{4}{{\epsilon }^2}{\Vert {\nabla }^m\left(u_2-u_1\right)\Vert }^4+\frac{4}{{\epsilon }^2}{\Vert {\nabla }^m\left(v_2-v_1\right)\Vert }^4+\frac{{\epsilon }^2{\Vert {\nabla }^{m+k}{\theta }_t\Vert }^2}{8}\right)\\ \le c_2\left(\frac{4}{{\epsilon }^2}\left({\lambda }_1^{-2k}+{\lambda }_1^{-2m-4k}\right)\left({\Vert {\nabla }^{m+k}\left(u_2-u_1\right)\Vert }^4+{\Vert {\nabla }^{2m+2k}\left(v_2-v_1\right)\Vert }^4+\right)\frac{{\epsilon }^2{\Vert {\nabla }^{m+k}{\theta }_t\Vert }^2}{8}\right),\end{array}$

$\begin{array}{c}\left(I_3,{\left(-\Delta \right)}^k{\theta }_t\right)\le \left(2M^{\prime }\left(s\right)\left[\left({\nabla }^m\left(u_1-u_2\right),{\nabla }^m{u}_1\right)+\left({\nabla }^m\left(v_1-v_2\right),{\nabla }^m{v}_1\right)\right]{\nabla }^{m+k}\left(u_2-u_1\right),{\nabla }^{m+k}{\theta }_t\right)\\ \le 2c_3\left(\frac{4}{{\epsilon }^2}{\Vert {\nabla }^m\left(u_2-u_1\right)\Vert }^4+\frac{4}{{\epsilon }^2}{\Vert {\nabla }^m\left(v_2-v_1\right)\Vert }^4+\frac{4}{{\epsilon }^2}{\Vert {\nabla }^{m+k}\left(u_2-u_1\right)\Vert }^4+\frac{3{\epsilon }^2{\Vert {\nabla }^{m+k}{\theta }_t\Vert }^2}{16}\right)\\ \le 2c_3\left(\frac{4}{{\epsilon }^2}\left({\lambda }_1^{-2k}+{\lambda }_1^{-2m-4k}+1\right)\left({\Vert {\nabla }^{m+k}\left(u_2-u_1\right)\Vert }^4+{\Vert {\nabla }^{2m+2k}\left(v_2-v_1\right)\Vert }^4\right)+\frac{3{\epsilon }^2{\Vert {\nabla }^{m+k}{\theta }_t\Vert }^2}{16}\right),\end{array}$

$\begin{array}{c}\left(I_4,{\left(-\Delta \right)}^k{\theta }_t\right)=\left(-2M^{\prime }\left(s\right)\left[\left({\nabla }^m\theta ,{\nabla }^m{u}_1\right)+\left({\nabla }^m\omega ,{\nabla }^m{v}_1\right)\right]{\nabla }^{m+k}{u}_1,{\nabla }^{m+k}{\theta }_t\right)\\ \le 2c_4\left(\frac{4{\Vert {\nabla }^{m+k}\theta \Vert }^2}{{\epsilon }^2}+\frac{4{\Vert {\nabla }^{m+k}\omega \Vert }^2}{{\epsilon }^2}+\frac{{\epsilon }^2{\Vert {\nabla }^{m+k}{\theta }_t\Vert }^2}{8}\right)\\ \le 2c_4\left(\frac{4}{{\epsilon }^2}\left({\lambda }_1^{-2k}+{\lambda }_1^{-2m-4k}\right)\left({\Vert {\nabla }^{m+k}\theta \Vert }^2+{\Vert {\nabla }^{2m+2k}\omega \Vert }^2\right)+\frac{{\epsilon }^2{\Vert {\nabla }^{m+k}{\theta }_t\Vert }^2}{8}\right),\end{array}$

which implies that

$\begin{array}{l}\left|\left(\left[M\left(s\right)-M\left(\tilde{s}\right)\right]{\left(-\Delta \right)}^m{u}_2+2M^{\prime }\left(s\right)\left[\left({\nabla }^{m}U,{\nabla }^m{u}_1\right)+\left({\nabla }^{m}V,{\nabla }^m{v}_1\right)\right]{\left(-\Delta \right)}^m{u}_1,{\left(-\Delta \right)}^k{\theta }_t\right)\right|\\ \le \frac{c_5}{{\epsilon }^2}\left({\lambda }_1^{-2k}+{\lambda }_1^{-2m-4k}+1\right)({\Vert {\nabla }^{m+k}\left(u_2-u_1\right)\Vert }^4+{\Vert {\nabla }^{2m+2k}\left(v_2-v_1\right)\Vert }^4\\ +{\Vert {\nabla }^{m+k}\theta \Vert }^2+{\Vert {\nabla }^{2m+2k}\omega \Vert }^2)+c_6{\epsilon }^2{\Vert {\nabla }^{m+k}{\theta }_t\Vert }^2.\end{array}$(57)

Analogously,

$\begin{array}{l}\left|\left(\left[M\left(s\right)-M\left(\tilde{s}\right)\right]{\left(-\Delta \right)}^{2m}{v}_2+2M^{\prime }\left(s\right)\left[\left({\nabla }^m{u}_1,{\nabla }^{m}U\right)+\left({\nabla }^m{v}_1,{\nabla }^{m}V\right)\right]{\left(-\Delta \right)}^{2m}{v}_1,{\left(-\Delta \right)}^{2k}{\omega }_t\right)\right|\\ \le \frac{c_7}{{\epsilon }^2}\left({\lambda }_1^{-2k}+{\lambda }_1^{-2m-4k}+1\right)({\Vert {\nabla }^{m+k}\left(u_2-u_1\right)\Vert }^4+{\Vert {\nabla }^{2m+2k}\left(v_2-v_1\right)\Vert }^4\\ +{\Vert {\nabla }^{m+k}\theta \Vert }^2+{\Vert {\nabla }^{2m+2k}\omega \Vert }^2)+c_8{\epsilon }^2{\Vert {\nabla }^{2m+2k}{\omega }_t\Vert }^2.\end{array}$(58)

Further,

$\begin{array}{l}\left(\frac{\partial g\left(u_{1t},v_1\right)}{\partial u_{1t}}{U}_t+\frac{\partial g\left(u_{1t},v_1\right)}{\partial v_1}V+g\left(u_{1t},v_1\right)-g\left(u_{2t},v_2\right),{\left(-\Delta \right)}^k{\theta }_t\right)\\ =\left(\frac{\partial g\left({\xi }_1,v_1\right)}{\partial {\xi }_1}\left(u_{1t}-u_{2t}\right)+\frac{\partial g\left(u_{2t},{\eta }_1\right)}{\partial {\eta }_1}\left(v_1-v_2\right)+\frac{\partial g\left(u_{1t},v_1\right)}{\partial u_{1t}}{U}_t+\frac{\partial g\left(u_{1t},v_1\right)}{\partial v_1}V,{\left(-\Delta \right)}^k{\theta }_t\right)\\ \le (\frac{{\partial }^{2}g\left({\xi }_2,v_1\right)}{\partial {\xi }_2}\left(u_{1t}-u_{2t}\right)\left({\xi }_1-u_{1t}\right)+\frac{\partial g\left(u_{1t},v_1\right)}{\partial u_{1t}}\left(-{\theta }_t\right)+\frac{{\partial }^{2}g\left(u_{2t},{\eta }_2\right)}{\partial ,{\eta }_2^2}\left(v_1-v_2\right)\left({\eta }_1-v_1\right)\\ +\frac{{\partial }^{2}g\left({\xi }_3,v_1\right)}{\partial {\xi }_3\partial v}\left(v_1-v_2\right)\left(u_{2t}-u_{1t}\right)+\frac{\partial g\left(u_{1t},v_1\right)}{\partial v_1}\left(-\omega \right),{\left(-\Delta \right)}^k{\theta }_t)\end{array}$

$\begin{array}{l}\le c_9{\lambda }_1^{-\frac{k}{2}-m}{\Vert {\nabla }^{m+k}\left(u_{1t}-u_{2t}\right)\Vert }^2\Vert {\nabla }^k{\theta }_t\Vert +c_{10}{\Vert {\nabla }^k{\theta }_t\Vert }^2+c_{11}{\lambda }_1^{-\frac{3k}{2}-2m}{\Vert {\nabla }^{2m+2k}\left(v_1-v_2\right)\Vert }^2\Vert {\nabla }^k{\theta }_t\Vert \\ +c_{12}{\lambda }_1^{-k-\frac{3m}{2}}\Vert {\nabla }^{m+k}\left(u_{1t}-u_{2t}\right)\Vert \Vert {\nabla }^{2k+2m}\left(v_1-v_2\right)\Vert \Vert {\nabla }^k{\theta }_t\Vert +c_{13}{\lambda }_1^{-\frac{k}{2}-m}\Vert {\nabla }^{2m+2k}\omega \Vert \Vert {\nabla }^k{\theta }_t\Vert \\ \le \frac{c_{14}}{2}\left(\frac{3}{2{\epsilon }^2}\left({\Vert {\nabla }^{m+k}\left(u_{1t}-u_{2t}\right)\Vert }^4+{\Vert {\nabla }^{2m+2k}\left(v_1-v_2\right)\Vert }^4\right)+{\epsilon }^2\left({\Vert {\nabla }^k{\theta }_t\Vert }^2+5{\Vert {\nabla }^{2m+2k}\omega \Vert }^2\right)\right),\end{array}$(59)

where,

$\begin{array}{l}{\xi }_1=r_1{u}_{1t}-\left(1-r_1\right)u_{2t},{\eta }_1=r_2{v}_1-\left(1-r_2\right)v_2,{\xi }_2=r_3{\xi }_1-\left(1-r_3\right)u_{1t},\\ {\eta }_2=r_4{\eta }_1-\left(1-r_4\right)v_1,{\xi }_3=r_5{u}_{2t}-\left(1-r_5\right)u_{1t},\\ c_{14}=\max \left\{c_9{\lambda }_1^{-\frac{k}{2}-m},\frac{2c_{10}}{{\epsilon }^2},c_{11}{\lambda }_1^{-\frac{3k}{2}-2m},c_{12}{\lambda }_1^{-k-\frac{3m}{2}},c_{13}{\lambda }_1^{-\frac{k}{2}-m}\right\}.\end{array}$

Similarly

$\begin{array}{l}\left(\frac{\partial g\left(u,v_t\right)}{\partial u}U+\frac{\partial g\left(u,v_t\right)}{\partial v_t}{V}_t+g\left(u,v_t\right)-g\left(\tilde{u},{\tilde{v}}_t\right),{\omega }_t\right)\\ \le \frac{c_{15}}{2}(\frac{3}{2{\epsilon }^2}\left({\Vert {\nabla }^{m+k}\left(u_1-u_2\right)\Vert }^4+{\Vert {\nabla }^{2m+2k}\left(v_{1t}-v_{2t}\right)\Vert }^4\right)\\ \stackrel{}{}+{\epsilon }^2\left({\Vert {\nabla }^{m+k}\theta \Vert }^2+5{\Vert {\nabla }^{2k}{\omega }_t\Vert }^2\right)).\end{array}$ (60)

Based on the above Equations (57)-(60), it is sorted out that

$\begin{array}{l}\frac{1}{2}\frac{\text{d}}{\text{d}t}\left\{{\Vert {\nabla }^k{\theta }_t\Vert }^2+{\Vert {\nabla }^{2k}{\omega }_t\Vert }^2+\mu \left({\Vert {\nabla }^{m+k}\theta \Vert }^2+{\Vert {\nabla }^{2m+2k}\omega \Vert }^2\right)\right\}\\ +\left(\beta -\left(c_6+c_8\right){\epsilon }^2\right)\left({\Vert {\nabla }^{m+k}{\theta }_t\Vert }^2+{\Vert {\nabla }^{2m+2k}{\omega }_t\Vert }^2\right)\\ \le c_{16}({\Vert {\nabla }^{m+k}\left(u_1-u_2\right)\Vert }^4+{\Vert {\nabla }^{m+k}\left(u_{1t}-u_{2t}\right)\Vert }^4+{\Vert {\nabla }^{2m+2k}\left(v_1-v_2\right)\Vert }^4\\ +{\Vert {\nabla }^{2m+2k}\left(v_{1t}-v_{2t}\right)\Vert }^4)+c_{17}\left({\Vert {\nabla }^k{\theta }_t\Vert }^2+{\Vert {\nabla }^{2k}{\omega }_t\Vert }^2+{\Vert {\nabla }^{m+k}\theta \Vert }^2+{\Vert {\nabla }^{2m+2k}\omega \Vert }^2\right),\end{array}$ (61)

where $\beta \ge \left(c_6+c_8\right){\epsilon }^2$,

by Gronwall inequality

$\begin{array}{l}{\Vert {\theta }_t\Vert }^2+{\Vert {\omega }_t\Vert }^2+\mu \left({\Vert {\nabla }^m\theta \Vert }^2+{\Vert {\nabla }^{2m}\omega \Vert }^2\right)\\ \le c_{18}{\text{e}}^{C_{19}t}{\displaystyle {\int }_0^t({\Vert {\nabla }^{m+k}\left(u_1-u_2\right)\Vert }^4}+{\Vert {\nabla }^{m+k}\left(u_{1t}-u_{2t}\right)\Vert }^4\\ +{\Vert {\nabla }^{2m+2k}\left(v_1-v_2\right)\Vert }^4+{\Vert {\nabla }^{2m+2k}\left(v_{1t}-v_{2t}\right)\Vert }^4)\text{d}\tau \\ \le c_{20}{\text{e}}^{c_{21}t}{\Vert \left(\xi ,\zeta ,\eta ,\sigma \right)\Vert }_{E_k}^4,\end{array}$ (62)

so as ${\Vert \left(\xi ,\zeta ,\eta ,\sigma \right)\Vert }_{E_k}^2\to 0$ in $E_k$, there are

$\frac{{\Vert S\left(t\right){\tilde{\rho }}_0-S\left(t\right){\rho }_0-\left(DS\left(t\right){\rho }_0\right)\left({\xi }_1,{\xi }_2,{\eta }_1,{\eta }_2\right)\Vert }_{E_k}^2}{{\Vert \left(\xi ,\zeta ,\eta ,\sigma \right)\Vert }_{E_k}^2}\le c_{20}{\text{e}}^{c_{21}t}{\Vert \left(\xi ,\zeta ,\eta ,\sigma \right)\Vert }_{E_k}^2\to 0.$

The differentiability of S(t) is proved.

The next step will be used in demonstrating the process of dimension estimation. it seems obvious that the Equations (1)-(5) also can be written as

${\phi }_t+H\left(\phi \right)=F\left(\phi \right),$ (63)

where $\phi ={\left(u,p,v,q\right)}^{\text{T}}$, $p=u_t+\epsilon u$, $q=v_t+\epsilon v$,

$H\left(\phi \right)=\left(\begin{array}{c}\epsilon u-p\\ {\left(-\Delta \right)}^{m}u-\epsilon \left(\beta {\left(-\Delta \right)}^m-\epsilon \right)u+\left(\beta {\left(-\Delta \right)}^m-\epsilon \right)p\\ \epsilon v-q\\ {\left(-\Delta \right)}^{2m}v-\epsilon \left(\beta {\left(-\Delta \right)}^{2m}-\epsilon \right)v+\left(\beta {\left(-\Delta \right)}^{2m}-\epsilon \right)q\end{array}\right),$ (64)

$F\left(\phi \right)=\left(\begin{array}{c}0\\ f_1\left(x\right)-g\left(u_t,v\right)+\left(1-M\left({\Vert {\nabla }^{m}u\Vert }^2+{\Vert {\nabla }^{m}v\Vert }^2\right)\right){\left(-\Delta \right)}^{m}u\\ 0\\ f_2\left(x\right)-g\left(u,v_t\right)+\left(1-M\left({\Vert {\nabla }^{m}u\Vert }^2+\Vert {\nabla }^{m}v\Vert \right)\right){\left(-\Delta \right)}^{2m}v\end{array}\right).$ (65)

Consider the first variation equation of (63)

${\Psi }^{\prime }+P\left(\phi \right)\Psi ={\Gamma }_1\left(\phi \right)\Psi +{\Gamma }_2\left(\phi \right)\Psi ,$ (66)

where $\Psi ={\left(U,P,V,Q\right)}^{\text{T}}$, $P=U_t+\epsilon U$, $Q=V_t+\epsilon V$, and $\phi ={\left(u,p,v,q\right)}^{\text{T}}$ is the solution to problem (63), and

$P\left(\phi \right)=\left(\begin{array}{cccc}\epsilon I& -I& 0& 0\\ \left(1-\beta \epsilon \right){\left(-\Delta \right)}^m+{\epsilon }^{2}I& \beta {\left(-\Delta \right)}^m-\epsilon I& 0& 0\\ 0& 0& \epsilon I& -I\\ 0& 0& \left(1-\beta \epsilon \right){\left(-\Delta \right)}^{2m}+{\epsilon }^{2}I& \beta {\left(-\Delta \right)}^{2m}-\epsilon I\end{array}\right),$(67)

${\Gamma }_1\left(\phi \right)=\left(\begin{array}{cccc}0& 0& 0& 0\\ \epsilon \frac{\partial g\left(u_t,v\right)}{\partial u_t}& -\frac{\partial g\left(u_t,v\right)}{\partial u_t}& -\frac{\partial g\left(u_t,v\right)}{\partial v}& 0\\ 0& 0& 0& 0\\ -\frac{\partial g\left(u,v_t\right)}{\partial u}& 0& \epsilon \frac{\partial g\left(u,v_t\right)}{\partial v_t}& -\frac{\partial g\left(u,v_t\right)}{\partial v_t}\end{array}\right),$ (68)

${\Gamma }_2\left(\phi \right)\Psi =\left(\begin{array}{c}0\\ \left(1-M\left(s\right)\right){\left(-\Delta \right)}^{m}U-2M^{\prime }\left(s\right)\left[\left({\nabla }^{m}u,{\nabla }^{m}U\right)+\left({\nabla }^{m}v,{\nabla }^{m}V\right)\right]{\left(-\Delta \right)}^{m}u\\ 0\\ \left(1-M\left(s\right)\right){\left(-\Delta \right)}^{2m}V-2M^{\prime }\left(s\right)\left[\left({\nabla }^{m}u,{\nabla }^{m}U\right)+\left({\nabla }^{m}v,{\nabla }^{m}V\right)\right]{\left(-\Delta \right)}^{2m}v\end{array}\right).$(69)

Theorem 4 Under the condition of Theorem 3, the global attractors of initial boundary value problems (1)-(5) have finite dimensional Hausdroff dimension

and fractal dimension, and then $d_H\left(A_k\right)<\frac{2}{3}{n}_0$, $d_F\left(A_k\right)<\frac{4}{3}{n}_0$.

Proof: For any fixed $\left(u_0,p_0,v_0,q_0\right)\in E_k$, assume that ${\chi }_1,{\chi }_2,\cdots ,{\chi }_n$ are $n_0$ elements in $E_k$, and ${\psi }_1\left(t\right),{\psi }_2\left(t\right),\cdots ,{\psi }_{n_0}\left(t\right)$ are $n_0$ solutions of the linearized Equation (66) with an initial value ${\psi }_1\left(0\right)={\chi }_1,{\psi }_2\left(0\right)={\chi }_2,\cdots ,{\psi }_{n_0}\left(0\right)={\chi }_{n_0}$, where $n_0$ is a natural number. It can be obtained by calculation

$\begin{array}{l}{\Vert {\psi }_1\left(t\right)\wedge \cdots \wedge {\psi }_{n_0}\left(t\right)\Vert }_{{\wedge }^{n_0}{E}_k}\\ \le {\Vert {\chi }_1\wedge \cdots \wedge {\chi }_n\Vert }_{{\wedge }^{n_0}{E}_k}\exp {\displaystyle {\int }_0^{t}tr{F}^{\prime }\left(\phi \left(\tau \right)\right)\circ Q_{n{\text{}}_0}\left(\tau \right)\text{d}\tau }.\end{array}$ (70)

where $\wedge$ represents the outer product, tr represents the trace of the operator, $Q_{n_0}\left(\tau \right)=Q_{n_0}\left(\tau ,{\phi }_0;{\chi }_1,{\chi }_2,\cdots ,{\chi }_{n_0}\right)$ represents the orthogonal projection from $E_k$ to $span\left\{{\psi }_1\left(t\right),{\psi }_2\left(t\right),\cdots ,{\psi }_{n_0}\left(t\right)\right\}$.

At a given time $\tau$, let $h_j\left(\tau \right)={\left({\xi }_j\left(\tau \right),{\zeta }_j\left(\tau \right),{\eta }_j\left(\tau \right),{\sigma }_j\left(\tau \right)\right)}^{\text{T}},j=1,2,\cdots ,n_0$

are the standard orthogonal basis of space $span\left\{{\psi }_1\left(t\right),{\psi }_2\left(t\right),\cdots ,{\psi }_{n_0}\left(t\right)\right\}$, then define the inner product on $E_k$

$\begin{array}{c}\left(h_j,\overline{h_j}\right)=\left({\nabla }^{m+k}{\xi }_j,{\nabla }^{m+k}\overline{{\xi }_j}\right)+\left({\nabla }^k{\zeta }_j,{\nabla }^k\overline{{\zeta }_j}\right)\\ +\left({\nabla }^{2m+2k}{\eta }_j,{\nabla }^{2m+2k}\overline{{\eta }_j}\right)+\left({\nabla }^{2k}{\sigma }_j,{\nabla }^{2k}\overline{{\sigma }_j}\right),\end{array}$

${\Vert h_j\Vert }_{E_k}^2={\left(h_j,h_j\right)}_{E_k}={\Vert {\nabla }^{m+k}{\xi }_j\Vert }^2+{\Vert {\nabla }^k{\zeta }_j\Vert }^2+{\Vert {\nabla }^{2m+2k}{\eta }_j\Vert }^2+{\Vert {\nabla }^{2k}{\sigma }_j\Vert }^2=1.$

Through the above conditions, can get

$\begin{array}{c}tr{F}^{\prime }\left(\phi \left(\tau \right)\right)\circ Q_{n_0}\left(\tau \right)={\displaystyle \underset{j=1}{\overset{+\infty }{\sum }}{\left(F^{\prime }\left(\phi \left(\tau \right)\right)\circ Q_{n_0}\left(\tau \right)h_j\left(\tau \right),h_j\left(\tau \right)\right)}_{E_k}}\\ ={\displaystyle \underset{j=1}{\overset{n_0}{\sum }}{\left(F^{\prime }\left(\phi \left(\tau \right)\right)h_j\left(\tau \right),h_j\left(\tau \right)\right)}_{E_k}},\end{array}$

by the Holder inequality, Young and Poincare inequality

$\begin{array}{l}\left(P\left(\phi \right)h_j,h_j\right)\\ =\epsilon \left({\nabla }^{m+k}{\xi }_j,{\nabla }^{m+k}{\xi }_j\right)-\left({\nabla }^{m+k}{\zeta }_j,{\nabla }^{m+k}{\xi }_j\right)+\left(1-\beta \epsilon \right)\left({\nabla }^{m+k}{\xi }_j,{\nabla }^{m+k}{\zeta }_j\right)\\ +{\epsilon }^2\left({\nabla }^k{\xi }_j,{\nabla }^k{\zeta }_j\right)+\beta \left({\nabla }^{m+k}{\zeta }_j,{\nabla }^{m+k}{\zeta }_j\right)-\epsilon \left({\nabla }^k{\zeta }_j,{\nabla }^k{\zeta }_j\right)\\ +\epsilon \left({\nabla }^{2m+2k}{\eta }_j,{\nabla }^{2m+2k}{\eta }_j\right)-\left({\nabla }^{2m+2k}{\sigma }_j,{\nabla }^{2m+2k}{\eta }_j\right)+{\epsilon }^2\left({\nabla }^{2k}{\eta }_j,{\nabla }^{2k}{\sigma }_j\right)\\ +\beta \left({\nabla }^{2m+2k}{\sigma }_j,{\nabla }^{2m+2k}{\sigma }_j\right)-\epsilon \left({\nabla }^{2k}{\sigma }_j,{\nabla }^{2k}{\sigma }_j\right)\end{array}$

$\begin{array}{l}+\left(1-\beta \epsilon \right)\left({\nabla }^{2m+2k}{\eta }_j,{\nabla }^{2m+2k}{\sigma }_j\right)\\ \ge \epsilon {\Vert {\nabla }^{m+k}{\xi }_j\Vert }^2-\left(2-\beta \epsilon \right)\Vert {\nabla }^{m+k}{\zeta }_j\Vert \Vert {\nabla }^{m+k}{\xi }_j\Vert -{\epsilon }^2\Vert {\nabla }^k{\xi }_j\Vert \Vert {\nabla }^k{\zeta }_j\Vert \\ +\beta {\Vert {\nabla }^{m+k}{\zeta }_j\Vert }^2-\epsilon {\Vert {\nabla }^k{\zeta }_j\Vert }^2-\left(2-\beta \epsilon \right)\Vert {\nabla }^{2m+2k}{\sigma }_j\Vert \Vert {\nabla }^{2m+2k}{\eta }_j\Vert \\ -{\epsilon }^2\Vert {\nabla }^{2k}{\eta }_j\Vert \Vert {\nabla }^{2k}{\sigma }_j\Vert +\beta {\Vert {\nabla }^{2m+2k}{\sigma }_j\Vert }^2-\epsilon {\Vert {\nabla }^{2k}{\sigma }_j\Vert }^2+\epsilon {\Vert {\nabla }^{2m+2k}{\eta }_j\Vert }^2\end{array}$

$\begin{array}{l}\ge \left(\beta {\lambda }_1^m-\frac{{\epsilon }^2+2\epsilon }{2}-\frac{c_{21}\left(2-\beta \epsilon \right){\lambda }_1^m}{2}\right){\Vert {\nabla }^k{\zeta }_j\Vert }^2-\frac{{\epsilon }^2}{2}{\Vert {\nabla }^k{\xi }_j\Vert }^2-\frac{{\epsilon }^2}{2}{\Vert {\nabla }^{2k}{\eta }_j\Vert }^2\\ +\left(\epsilon -\frac{c_{21}\left(2-\beta \epsilon \right)}{2}\right){\Vert {\nabla }^{m+k}{\xi }_j\Vert }^2+\left(\epsilon -\frac{c_{22}\left(2-\beta \epsilon \right)}{2}\right){\Vert {\nabla }^{2m+2k}{\eta }_j\Vert }^2\\ +\left(\beta {\lambda }_1^{2m}-\frac{{\epsilon }^2+2\epsilon }{2}-\frac{c_{22}\left(2-\beta \epsilon \right){\lambda }_1^{2m}}{2}\right){\Vert {\nabla }^{2k}{\sigma }_j\Vert }^2,\end{array}$(71)

$\begin{array}{l}\left({\Gamma }_1\left(\phi \right)h_j,h_j\right)\\ =\epsilon \left(\frac{\partial g\left(u_t,v\right)}{\partial u_t}{\nabla }^k{\xi }_j,{\nabla }^k{\zeta }_j\right)-\left(\frac{\partial g\left(u_t,v\right)}{\partial u_t}{\nabla }^k{\zeta }_j,{\nabla }^k{\zeta }_j\right)\\ -\left(\frac{\partial g\left(u_t,v\right)}{\partial v}{\nabla }^k{\eta }_j,{\nabla }^k{\zeta }_j\right)-\left(\frac{\partial g\left(u,v_t\right)}{\partial u}{\nabla }^{2k}{\xi }_j,{\nabla }^{2k}{\sigma }_j\right)\\ +\epsilon \left(\frac{\partial g\left(u,v_t\right)}{\partial v_t}{\nabla }^{2k}{\eta }_j,{\nabla }^{2k}{\sigma }_j\right)-\left(\frac{\partial g\left(u,v_t\right)}{\partial v_t}{\nabla }^{2k}{\sigma }_j,{\nabla }^{2k}{\sigma }_j\right)\end{array}$

$\begin{array}{l}\le \epsilon {\Vert \frac{\partial g\left(u_t,v\right)}{\partial u_t}\Vert }_{\infty }\Vert {\nabla }^k{\xi }_j\Vert \Vert {\nabla }^k{\zeta }_j\Vert -{\Vert \frac{\partial g\left(u_t,v\right)}{\partial u_t}\Vert }_{\infty }{\Vert {\nabla }^k{\zeta }_j\Vert }^2\\ -{\Vert \frac{\partial g\left(u_t,v\right)}{\partial v}\Vert }_{\infty }\Vert {\nabla }^k{\eta }_j\Vert \Vert {\nabla }^k{\zeta }_j\Vert -{\Vert \frac{\partial g\left(u,v_t\right)}{\partial u}\Vert }_{\infty }\Vert {\nabla }^{2k}{\xi }_j\Vert \Vert {\nabla }^{2k}{\sigma }_j\Vert \\ +\epsilon {\Vert \frac{\partial g\left(u,v_t\right)}{\partial v_t}\Vert }_{\infty }\Vert {\nabla }^{2k}{\eta }_j\Vert \Vert {\nabla }^{2k}{\sigma }_j\Vert -{\Vert \frac{\partial g\left(u,v_t\right)}{\partial v_t}\Vert }_{\infty }{\Vert {\nabla }^{2k}{\sigma }_j\Vert }^2\\ \le \frac{{\epsilon }^2{c}_{23}}{2}{\Vert {\nabla }^k{\xi }_j\Vert }^2+\frac{c_{23}-c_{24}}{2}{\Vert {\nabla }^k{\zeta }_j\Vert }^2+\frac{{\epsilon }^2{c}_{25}}{2}{\Vert {\nabla }^{2k}{\eta }_j\Vert }^2+\frac{c_{23}-c_{26}}{2}{\Vert {\nabla }^{2k}{\sigma }_j\Vert }^2,\end{array}$(72)

$\begin{array}{l}\left({\Gamma }_2\left(\phi \right)h_j,h_j\right)\\ =\left(1-M\left(s\right)\right)\left({\nabla }^{m+k}{\xi }_j,{\nabla }^{m+k}{\zeta }_j\right)-2M^{\prime }\left(s\right)\left({\nabla }^{m}u,{\nabla }^m{\xi }_j\right)\left({\nabla }^{2m+k}u,{\nabla }^k{\zeta }_j\right)\\ -2M^{\prime }\left(s\right)\left({\nabla }^{m}v,{\nabla }^m{\eta }_j\right)\left({\nabla }^{2m+k}u,{\nabla }^k{\zeta }_j\right)+\left(1-M\left(s\right)\right)\left({\nabla }^{2m+2k}{\eta }_j,{\nabla }^{2m+2k}{\sigma }_j\right)\\ -2M^{\prime }\left(s\right)\left({\nabla }^{m}u,{\nabla }^m{\xi }_j\right)\left({\nabla }^{4m+2k}v,{\nabla }^{2k}{\sigma }_j\right)\\ -2M^{\prime }\left(s\right)\left({\nabla }^{m}v,{\nabla }^m{\eta }_j\right)\left({\nabla }^{4m+2k}v,{\nabla }^{2k}{\sigma }_j\right)\end{array}$

$\begin{array}{l}\le \left(1-{\delta }_0\right){\lambda }_1^{\frac{m}{2}}\left(\Vert {\nabla }^{m+k}{\xi }_j\Vert \Vert {\nabla }^k{\zeta }_j\Vert +\Vert {\nabla }^{2m+2k}{\eta }_j\Vert \Vert {\nabla }^{2k}{\sigma }_j\Vert \right)\\ +2c_{26}\left({\lambda }_1^{-\frac{k}{2}}+{\lambda }_1^{-\frac{m}{2}-k}\right)\left(\Vert {\nabla }^{m+k}{\xi }_j\Vert \Vert {\nabla }^k{\zeta }_j\Vert +\Vert {\nabla }^{2m+2k}{\eta }_j\Vert \Vert {\nabla }^k{\zeta }_j\Vert \right)\\ +2c_{27}\left({\lambda }_1^{-\frac{k}{2}}+{\lambda }_1^{-\frac{m}{2}-k}\right)\left(\Vert {\nabla }^{m+k}{\xi }_j\Vert \Vert \Vert {\nabla }^{2k}{\sigma }_j\Vert \Vert +\Vert {\nabla }^{2m+2k}{\eta }_j\Vert \Vert {\nabla }^{2k}{\sigma }_j\Vert \right)\end{array}$

$\begin{array}{l}\le \frac{1-{\delta }_0}{2}{\lambda }_1^{\frac{m}{2}}\left({\Vert {\nabla }^{m+k}{\xi }_j\Vert }^2+{\Vert {\nabla }^k{\zeta }_j\Vert }^2+{\Vert {\nabla }^{2m+2k}{\eta }_j\Vert }^2+{\Vert {\nabla }^{2k}{\sigma }_j\Vert }^2\right)\\ +2c_{28}\left({\lambda }_1^{-\frac{k}{2}}+{\lambda }_1^{-\frac{m}{2}-k}\right)\left({\Vert {\nabla }^{m+k}{\xi }_j\Vert }^2+{\Vert {\nabla }^k{\zeta }_j\Vert }^2+{\Vert {\nabla }^{2m+2k}{\eta }_j\Vert }^2+{\Vert {\nabla }^{2k}{\sigma }_j\Vert }^2\right)\\ \le \frac{1-{\delta }_0}{2}{\lambda }_1^{\frac{m}{2}}+2c_{28}\left({\lambda }_1^{-\frac{k}{2}}+{\lambda }_1^{-\frac{m}{2}-k}\right).\end{array}$(73)

Based on the above Equations (71)-(73), it is sorted out that

$\begin{array}{l}{\left(F^{\prime }\left(\phi \left(\tau \right)\right)h_j\left(\tau \right),h_j\left(\tau \right)\right)}_{E_k}=\left(\left(-P\left(\phi \right)+{\Gamma }_1\left(\phi \right)+{\Gamma }_2\left(\phi \right)\right)h_j,h_j\right)\\ \le -\left(\beta {\lambda }_1^m-\frac{{\epsilon }^2+2\epsilon }{2}-\frac{c_{21}\left(2-\beta \epsilon \right){\lambda }_1^m}{2}-\frac{c_{23}-c_{24}}{2}\right){\Vert {\nabla }^k{\zeta }_j\Vert }^2\\ -\left(\beta {\lambda }_1^{2m}-\frac{{\epsilon }^2+2\epsilon }{2}-\frac{c_{22}\left(2-\beta \epsilon \right){\lambda }_1^{2m}}{2}-\frac{c_{25}-c_{26}}{2}\right){\Vert {\nabla }^{2k}{\sigma }_j\Vert }^2\\ -\left(\epsilon -\frac{c_{21}\left(2-\beta \epsilon \right)}{2}\right){\Vert {\nabla }^{m+k}{\xi }_j\Vert }^2-\left(\epsilon -\frac{c_{22}\left(2-\beta \epsilon \right)}{2}\right){\Vert {\nabla }^{2m+2k}{\eta }_j\Vert }^2\\ +\frac{{\epsilon }^2\left(c_{23}-1\right)}{2}{\Vert {\nabla }^k{\xi }_j\Vert }^2+\frac{{\epsilon }^2\left(c_{25}-1\right)}{2}{\Vert {\nabla }^{2k}{\eta }_j\Vert }^2+\frac{1-{\delta }_0}{2}{\lambda }_1^{\frac{m}{2}}\\ +2c_{28}\left({\lambda }_1^{-\frac{k}{2}}+{\lambda }_1^{-\frac{m}{2}-k}\right),\end{array}$ (74)

let

$\begin{array}{c}b=\min \{\beta {\lambda }_1^m-\frac{{\epsilon }^2+2\epsilon }{2}-\frac{c_{21}\left(2-\beta \epsilon \right){\lambda }_1^m}{2}-\frac{c_{23}-c_{24}}{2},\epsilon -\frac{c_{21}\left(2-\beta \epsilon \right)}{2},\\ \epsilon -\frac{c_{22}\left(2-\beta \epsilon \right)}{2},\beta {\lambda }_1^{2m}-\frac{{\epsilon }^2+2\epsilon }{2}-\frac{c_{22}\left(2-\beta \epsilon \right){\lambda }_1^{2m}}{2}-\frac{c_{25}-c_{26}}{2},\\ -\frac{1-{\delta }_0}{2}{\lambda }_1^{\frac{m}{2}}-2c_{28}\left({\lambda }_1^{-\frac{k}{2}}+{\lambda }_1^{-\frac{m}{2}-k}\right)\},\end{array}$

$a=\max \left\{\frac{{\epsilon }^2\left(c_{23}-1\right)}{2},\frac{{\epsilon }^2\left(c_{25}-1\right)}{2}\right\},$

we can obtain

${\displaystyle \underset{j=1}{\overset{n_0}{\sum }}{\left(F^{\prime }\left(\phi \left(\tau \right)\right)h_j\left(\tau \right),h_j\left(\tau \right)\right)}_{E_k}}\le -n_{0}b+a{\displaystyle \underset{j=1}{\overset{n_0}{\sum }}\left({\Vert {\nabla }^k{\xi }_j\Vert }^2+{\Vert {\nabla }^{2k}{\eta }_j\Vert }^2\right)},$ (75)

for almost all times t, there is ${\displaystyle \underset{i=1}{\overset{n_0}{\sum }}{\Vert {\nabla }^k{\xi }_i\Vert }^2}\le {\displaystyle \underset{i=1}{\overset{n_0}{\sum }}{\lambda }_i^{s^{\prime }-1}}$, ${\displaystyle \underset{j=1}{\overset{n_0}{\sum }}{\Vert {\nabla }^{2k}{\eta }_j\Vert }^2}\le {\displaystyle \underset{j=1}{\overset{n_0}{\sum }}{\lambda }_j^{s^{″}-1}}$, so

$tr{F}^{\prime }\left(\phi \left(\tau \right)\right)\circ Q_{n_0}\left(\tau \right)\le -n_{0}b+a\left({\displaystyle \underset{i=1}{\overset{n_0}{\sum }}{\lambda }_i^{s^{\prime }-1}}+{\displaystyle \underset{j=1}{\overset{n_0}{\sum }}{\lambda }_j^{s^{″}-1}}\right),$ (76)

because of

$q_{n_0}\left(t\right)=\underset{{\phi }_0\in B_{0k}}{\sup }\underset{{\Psi }_j(0)\in E_k}{\sup }\left\{\frac{1}{t}{\displaystyle {\int }_0^{t}tr{F}^{\prime }\left(\phi \left(\tau \right)\right)\circ Q_{n_0}\left(\tau \right)\text{d}\tau }\right\},q_{n_0}=\underset{t\to \infty }{\lim }{q}_{n_0}\left(t\right),$ (77)

$q_{n_0}\left(t\right)\le -n_{0}b+a\left({\displaystyle \underset{i=1}{\overset{n_0}{\sum }}{\lambda }_i^{s^{\prime }-1}}+{\displaystyle \underset{j=1}{\overset{n_0}{\sum }}{\lambda }_j^{s^{″}-1}}\right),q_{n_0}\le -n_{0}b+a\left({\displaystyle \underset{i=1}{\overset{n_0}{\sum }}{\lambda }_i^{s^{\prime }-1}}+{\displaystyle \underset{j=1}{\overset{n_0}{\sum }}{\lambda }_j^{s^{″}-1}}\right),$ (78)

Therefore, the Lyapunov exponent $K_1,K_2,\cdots ,K_{n_0}$ on set $B_{0k}$ is uniformly bounded, and

$K_1+K_2+\cdots +K_{n_0}\le -n_{0}b+a\left({\displaystyle \underset{i=1}{\overset{n_0}{\sum }}{\lambda }_i^{s^{\prime }-1}}+{\displaystyle \underset{j=1}{\overset{n_0}{\sum }}{\lambda }_j^{s^{″}-1}}\right),$ (79)

so

${\left(q_{ij}\right)}_{+}\le -n_{0}b+a\left({\displaystyle \underset{i=1}{\overset{n_0}{\sum }}{\lambda }_i^{s^{\prime }-1}}+{\displaystyle \underset{j=1}{\overset{n_0}{\sum }}{\lambda }_j^{s^{″}-1}}\right)\le a\left({\displaystyle \underset{i=1}{\overset{n_0}{\sum }}{\lambda }_i^{s^{\prime }-1}}+{\displaystyle \underset{j=1}{\overset{n_0}{\sum }}{\lambda }_j^{s^{″}-1}}\right)\le \frac{2}{5}{n}_{0}b,$ (80)

$q_{n_0}\le -n_{0}b\left(1-\frac{a}{n_{0}b}\left({\displaystyle \underset{i=1}{\overset{n_0}{\sum }}{\lambda }_i^{s^{\prime }-1}}+{\displaystyle \underset{j=1}{\overset{n_0}{\sum }}{\lambda }_j^{s^{″}-1}}\right)\right)\le -\frac{3}{5}{n}_{0}b,$ (81)

further

$\underset{1\le i,j\le n_0}{\max }\frac{{\left(q_{ij}\right)}_{+}}{\left|q_{n_0}\right|}\le \frac{2}{3}.$ (82)

Thus, can obtain $d_H\left(A_k\right)<\frac{2}{3}{n}_0$, $d_F\left(A_k\right)<\frac{4}{3}{n}_0$.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Lin, G.G. and Gao, Y.L. (2017) The Global and Exponential Attractors for the Higher-Order Kirchhoff-Type Equation with Strong Linear Damping. Journal of Mathematics Research, 9, 145-167. https://doi.org/10.5539/jmr.v9n4p145
[2]	Lin, G.G. and Zhang, M. (2018) The Estimates of the Upper Bounds of Hausdorff Dimensions for the Global Attractor for a Class of Nonlinear Coupled Kirchhoff-Type Equations. Advances in Pure Mathematics, 8, 1-10. https://doi.org/10.4236/apm.2018.81001
[3]	Lin, G.G. and Yang, L.J. (2021) Global Attractors and Their Dimension Estimates for a Class of Generalized Kirchhoff Equations. Advances in Pure Mathematics, 11, 317-333. https://doi.org/10.4236/apm.2021.114020
[4]	Lin, G.G. and Chen, Y.H. (2020) A Family of Global Attractors for the Higher-Order Kirchhoff-Type Equations and Its Dimension Estimation. American Journal of Applied Mathematics, 8, 300-310. https://doi.org/10.11648/j.ajam.20200806.12
[5]	Lin, G.G. and Zhou, C.M. (2021) The Dynamic Behavior of a Class of Kirchhoff Equations with High Order Strong Damping. Journal of Applied Mathematics and Physics, 9, 1041-1055. https://doi.org/10.4236/jamp.2021.95071
[6]	Lin, G.G. and Hu, L.J. (2018) Estimate on the Dimension of Global Attractor for Nonlinear Higher-Order Coupled Kirchhoff Type Equations. Advances in Pure Mathematics, 8, 11-24.
[7]	Cavalcanti, M.M., Domingos Cavalcanti, V.N., Prates Filho, J.S. and Soriano, J.A. (1998) Existence and Exponential Decay for a Kirchhoff-Carrier Model with Viscosity. Journal of Mathematical Analysis and Applications, 226, 40-60. https://doi.org/10.1006/jmaa.1998.6057
[8]	Ono, K. (1997) Global Existence, Decay, and Blowup of Solutions for Some Mildly Degenerate Nonlinear Kirchhoff Strings. Journal of Differential Equations, 137, 273-301. https://doi.org/10.1006/jdeq.1997.3263
[9]	Lin, G.G. (2011) Nonlinear Evolution Equations. Yunnan University Press, Kunming.
[10]	Nakao, M. (2008) An Attractor for a Nonlinear Dissipative Wave Equation of Kirchhoff Type. Journal of Mathematical Analysis and Applications, 353, 652-659. https://doi.org/10.1016/j.jmaa.2008.09.010
[11]	Lin, G.G. (2019) Dynamic Properties of Several Kinds of the Kirchhoff Equations. Chongqing University Press, Chongqing.
[12]	Lin, G.G. and Guan, L.P. (2019) Random Attractor Family for a Class of Stochastic Higher-Order Kirchhoff Equations. Journal of Mathematics Research, 11, 23-33. https://doi.org/10.5539/jmr.v11n3p23
[13]	Wang, M.X., Tian, C.C. and Lin, G.G. (2014) Global Attractor and Dimension Estimation for a 2D Generalized Anisotropy Kuramoto-Sivashinsky Equation. International Journal of Modern Nonlinear Theory and Application, 3, 163-172. https://doi.org/10.4236/ijmnta.2014.34018
[14]	Tchaho, C.T.D., Omanda, H.M., Mbourou, G.N., Bogning, J.R. and Kofané, T.C. (2021) Higher Order Solitary Wave Solutions of the Standard KdV Equations. Open Journal of Applied Sciences, 11, 103-125. https://doi.org/10.4236/ojapps.2021.111008
[15]	Qin, C.L., Du, J.J. and Lin. G.G. (2017) Random Attractors for the Kirchhoff-Type Suspension Bridge Equations with Strong Damping and White Noises. International Journal of Modern Nonlinear Theory and Application, 6, 134-147. https://doi.org/10.4236/ijmnta.2017.64012
[16]	Liao, Y.H., Lin, G.G. and Liu, J. (2022) A Family of Global Attractors for a Class of Generalized Kirchhoff-Beam Equations. Journal of Applied Mathematics and Physics, 10, 930-951. https://doi.org/10.4236/jamp.2022.103064
[17]	Kirchhoff, G. (1883) Vorlesungen Fiber Mechanic. Tenbner, Stuttgart.
[18]	Hosoya, M. and Yamada, Y. (1991) On Some Nonlinear Wave Equations II: Global Existence and Energy Decay of Solutions. Journal of the Faculty of Science, the University of Tokyo. Section 1A, Mathematics, 38, 239-250.