Random Attractors for the Kirchhoff-Type Suspension Bridge Equations with Strong Damping and White Noises

Chuangliang Qin¹, Jinji Du¹, Guoguang Lin^2*

¹School of Mathematics and Information, Xinyang College, Xinyang, China.
²School of Mathematics and Statistics, Yunnan University, Kunming, China.

DOI: 10.4236/ijmnta.2017.64012                

Abstract

In this paper, we investigate the existence of random attractor for the random dynamical system generated by the Kirchhoff-type suspension bridge equations with strong damping and white noises. We first prove the existence and uniqueness of solutions to the initial boundary value conditions, and then we study the existence of the global attractors of the equation.

Keywords

Kirchhoff-Type Suspension Bridge Equations, Random Attractors, Random Dynamical System

1. Introduction

In this paper, we consider the following stochastic Kirchhoff-type suspension bridge equations

$\{\begin{array}{ll}{u}_{tt}+{\Delta }^{2}u+{\Delta }^2{u}_t+\left(p-{\left|\nabla u\right|}^2\right)\Delta u+b{u}^{+}+f\left(u\right)=q\left(x\right)\stackrel{\dot{}}{W},\hfill & \Omega \times \left[\tau ,+\infty \right),\tau \in ℝ,\hfill \\ u\left(x,t\right)=\nabla u\left(x,t\right)=0,\hfill & x\in \partial \Omega ,t\ge \tau ,\hfill \\ u\left(x,\tau \right)=u_0\left(x\right),u_t\left(x,\tau \right)=u_1\left(x\right)\hfill & x\in \Omega ,\hfill \end{array}$ (1.1)

where $u\left(x\mathrm{,}t\right)$ is an unknown function,which represents the downward deflection of the road bed in the vertical plane, $u^{+}=u$ for $u\ge 0$ and $u^{+}=0$ for $u<0$ . $b>0$ denotes the spring constant of the ties, the real constant p represents the axial force acting at the end of the road bed of the bridge in the reference configuration. Namely, p is negative when the bridge is stretched, positive when compressed. $\Omega$ is an open bounded subset of ${ℝ}^2$ with sufficiently smooth boundary $\partial \Omega$ . $q\left(x\right)\in H^3\left(\Omega \right)$ is not identically equal to zero, $f$ is a nonlinear function satisfying certain conditions. $\stackrel{\dot{}}{W}$ is the derivative of a one-dimensional two-valued Wiener process $W\left(t\right)$ and $q\left(x\right)\stackrel{\dot{}}{W}$ formally describes white noise.

We assume that the nonlinear function $f\in C^2\left(ℝ\mathrm{,}ℝ\right)$ with $f\left(0\right)=0$ , which satisfies the following assumptions:

(a) Growth conditions:

$\left|f\left(s\right)\right|\le C_0\left(1+{\left|s\right|}^k\right)\mathrm{,}k\ge \mathrm{1,}\forall s\in ℝ\mathrm{,}$ (1.2)

where $C_0$ is a positive constant. For example, obviously, $f\left(s\right)={\left|s\right|}^{k-1}s$ satisfies (1.2).

(b) Dissipation conditions:

$F\left(s\right)\mathrm{<mo>\; :\; </mo>}={\displaystyle {\int }_0^s}f\left(r\right)\text{d}r\ge C_1\left({\left|s\right|}^{k+1}-1\right)\mathrm{,}\forall s\in ℝ\mathrm{,}$ (1.3)

and

$sf\left(s\right)\ge C_2\left(F\left(s\right)-1\right)\mathrm{,}\forall s\in ℝ\mathrm{,}$ (1.4)

where $C_1\mathrm{,}{C}_2$ are positive constants.

When $f\left(u\right)=0$ and $q\left(x\right)=0$ , Equation (1.1) is regarded as a model of naval structures,which is originally in [1] introduced by Lazer and McKenna. To the best of our knowledge, Qin [2] [3] proved random attractor for stochastic Kirchhoff equation with white noise, Ma [4] investigated the asymptotic behavior of the solution for the floating beam, that is, the “noise” is absent in (1.1). No one else has studied the long-time behavior of the solutions about these problems, it is just our interest in this paper. As far as the other related problems are concerned, we refer the reader to [2] - [7] and the references therein.

It is well known that Crauel and Flandoli originally introduced the random attractor for the infinite-dimensional RDS [8] [9] . A random attractor of RDS is a measurable and compact invariant random set attracting all orbits. It is the appropriate generalization of the now classical attractor exists, it is the smallest attracting compact set and the largest invariant set [10] . Zhou et al. [11] studied random attractor for damped nonlinear wave equation with white noise. Fan [12] proved random attractor for a damped stochastic wave equation with multiplicative noise. These abstract results have been successfully applied to many stochastic dissipative partial differential equations. The existence of a random attractors for the wave equations has been investigated by several authors [8] [9] [10] .

The outline of this paper is as follows: In Section 2, we recall many basic concepts related to a random attractor for genneral random dynamical system. In Section 3, We prove the existence and uniqueness of the solution corresponding to system (1.1) which determines RDS. In Section 4, we prove the existence of random attractor of the random dynamical system.

2. Random Dynamical System

In this section, we recall some basic concepts related to RDS and a random attractor for RDS in [8] [9] [10] , which are important for getting our main results.

Let $\left(X\mathrm{,}{\Vert \cdot \Vert }_X\right)$ be a separable Hilbert space with Borel s-algebra $\mathcal{B}\left(X\right)$ , and let $\left(\Omega \mathrm{,}\mathcal{F}\mathrm{,}\mathcal{P}\right)$ be a probability space. ${\theta }_t\mathrm{<mo>\; :\; </mo>}\Omega \to \Omega \mathrm{,}t\in ℝ$ is a family of measure preserving transformations such that $\left(t\mathrm{,}\omega \right)\mapsto {\theta }_t\omega$ is measurable, ${\theta }_0=id$ and ${\theta }_{t+s}={\theta }_t{\theta }_s$ for all $t\mathrm{,}s\in ℝ$ . The flow ${\theta }_t$ together with the probability space $\left(\Omega \mathrm{,}\mathcal{F}\mathrm{,}\mathcal{P}\mathrm{,}{\left({\theta }_t\right)}_{t\in ℝ}\right)$ is called a metric dynamical system.

Definition 2.1. Let $\left(\Omega \mathrm{,}\mathcal{F}\mathrm{,}\mathcal{P}\mathrm{,}{\left({\theta }_t\right)}_{t\in ℝ}\right)$ be a metric dynamical system. Suppose that the mapping $\varphi \mathrm{<mo>\; :\; </mo>}{ℝ}^{+}\times \Omega \times X\to X$ is $\left(\mathcal{B}\left({ℝ}^{+}\right)\times \mathcal{F}\times \mathcal{B}\left(X\right)\right)$ -measurable and satisfies the following properties:

1) $\varphi \left(0,\omega \right)x=x,x\in X$ and $\omega \in \Omega$ ;

2) $\varphi \left(t+s,\omega \right)=\varphi \left(t,{\theta }_s\omega \right)\circ \varphi \left(s,\omega \right)$ , for all $t\mathrm{,}s\in {ℝ}^{+}\mathrm{,}x\in X$ and $\omega \in \Omega$ .

Then $\varphi$ is called a random dynamical system (RDS). Moreover, $\varphi$ is called a continuous RDS if $\varphi$ is continuous with respect to x for $t\ge 0$ and $\omega \in \Omega$ .

Definition 2.2. A set-valued map $D\mathrm{<mo>\; :\; </mo>}\Omega \to 2^X$ is said to be a closed (compact)random set if $D\left(\omega \right)$ is closed (compact) for $P-a\mathrm{.}s\mathrm{.}\omega \in \Omega$ , and $\omega \mapsto d\left(x\mathrm{,}D\left(\omega \right)\right)$ is $p-a\mathrm{.}s\mathrm{.}$ measurable for all $x\in X$ .

Definition 2.3. If K and B are random sets such that for $P-a\mathrm{.}s\mathrm{.}\omega$ there exists a time $t_B\left(\omega \right)$ such that for all $t\ge t_B\left(\omega \right)$ ,

$\varphi \left(t\mathrm{,}{\theta }_{-t}\omega \right)B\left({\theta }_{-t}\omega \right)\subset K\left(\omega \right)\mathrm{,}$

then K is said to absorb B, and $t_B\left(\omega \right)$ is called the absorption time.

Definition 2.4. A random set $\mathcal{A}=A{\left(\omega \right)}_{\omega \in \Omega }\subset X$ is called a random attractor associated to the RDS $\varphi$ if $P-a\mathrm{.}s\mathrm{.}$ :

1) $\mathcal{A}$ is a random compact set, i.e., $A\left(\omega \right)$ is compact for $P-a\mathrm{.}s\mathrm{.}\omega \in \Omega$ , and the map $\omega \mapsto d\left(x\mathrm{,}A\left(\omega \right)\right)$ is measurable for every $x\in X$ ;

2) $\mathcal{A}$ is $\varphi$ -invariant, i.e., $\varphi \left(t\mathrm{,}\omega \right)A\left(\omega \right)=A\left({\theta }_t\omega \right)$ for all $t\ge 0$ and $P-a\mathrm{.}s\mathrm{.}\omega \in \Omega$ ;

3) $\mathcal{A}$ attracts every set B in X, i.e., for all bounded (and non-random) $B\subset X$ ,

$\underset{t\to \infty }{lim}d\left(\varphi \left(t\mathrm{,}{\theta }_{-t}\omega \right)B\left({\theta }_{-t}\omega \right)\mathrm{,}A\left(\omega \right)\right)=\mathrm{0,}$

where $d\left(\cdot \mathrm{,}\cdot \right)$ denotes the Hausdorff semi-distance:

$d\left(A,B\right)=\underset{x\in A}{\sup }\underset{y\in B}{\inf }d\left(x,y\right),A,B\in X.$

Note that $\varphi \left(t\mathrm{,}{\theta }_{-t}\omega \right)x$ can be interpreted as the position of the trajectory which was in x at time $-t$ . Thus, the attraction property holds from $t=-\infty$ .

Theorem 2.1. [8] (Existence of a random attractor) Let $\varphi$ be a continuous random dynamical system on X over $\left(\Omega \mathrm{,}\mathcal{F}\mathrm{,}\mathcal{P}\mathrm{,}{\left({\theta }_t\right)}_{t\in ℝ}\right)$ . Suppose that there exists a random compact set $K\left(\omega \right)$ absorbing every bounded non-random set $B\subset X$ . Then the set

$\mathcal{A}=A{\left(\omega \right)}_{\omega \in \Omega }=\overline{{\displaystyle \underset{B\subset X}{\cup }}{\wedge }_B\left(\omega \right)}$

is a global random attractor for $\varphi$ , where the union is taken over all bounded $B\subset X$ , and ${\wedge }_B\left(\omega \right)$ is the w-limitsset of B given by

${\wedge }_B\left(\omega \right)={\displaystyle \underset{s\ge 0}{\cap }}\overline{{\displaystyle \underset{t\ge s}{\cup }}\varphi \left(t,{\theta }_{-t}\omega \right)B\left({\theta }_{-t}\omega \right)}$

3. Existence and Uniqueness of Solutions

With the usual notation, we denote

$\begin{array}{l}H=L^2\left(\Omega \right),V=H^2\left(\Omega \right)\\ D\left(A\right)=H^2\left(\Omega \right)\cap H_0^1\left(\Omega \right),D\left(A^2\right)=\left\{u\in H^4\left(\Omega \right):A^{2}u\in L^2\left(\Omega \right)\right\},\end{array}$

where $A=-\Delta ,A^2={\Delta }^2$ . We denote $H\mathrm{,}V$ with the following inner products and norms,respectively:

$\begin{array}{l}\left(u,v\right)={\displaystyle {\int }_{\Omega }}uvdx,{\Vert u\Vert }^2=\left(u,u\right),\forall u,v\in H,\\ \left(\left(u,v\right)\right)={\displaystyle {\int }_{\Omega }}\Delta u\Delta vdx,{\Vert u\Vert }_2^2=\left(\left(u,u\right)\right),\forall u,v\in V.\end{array}$

And we introduce the space $E=D\left(A\right)\times H$ , which is used throughout the paper and endow the space E with the following usual scalar product and norm:

$\begin{array}{l}{\left(y_1,y_2\right)}_E=\left(\left(u_1,u_2\right)\right)+\left(v_1,v_2\right),\forall y_i={\left(u_i,v_i\right)}^{\text{T}}\in E,i=1,2,\\ {\Vert y\Vert }_E^2={\left(y,y\right)}_E,\forall y={\left(u,v\right)}^{\text{T}}\in E\end{array}$

where T denotes the transposition.

More generally, define ${\mathcal{H}}^r=D\left(A^{\frac{r}{2}}\right)$ for $r\in ℝ$ , which turns out to be a Hilbert space with the inner product ${\left(u,v\right)}_r=\left(A^{\frac{r}{2}}u,A^{\frac{r}{2}}v\right)$ , we denote by ${\Vert \cdot \Vert }_r=\Vert A^{\frac{r}{2}}\Vert$ the norm on ${\mathcal{H}}^r$ induced by the above inner product. Let $\lambda >0$

be the first eigenvalue of $A^{2}u=\lambda u,u\left(x,t\right)=\nabla u\left(x,t\right)=0,x\in \partial \Omega$ , by the compact embeddings ${\mathcal{H}}^{r+1}\to {\mathcal{H}}^r$ along with the generalized Poincaré inequality, we have

${\Vert u\Vert }_{r+1}^2\ge \lambda {\Vert u\Vert }_r^2\mathrm{.}$ (3.1)

It is convenient to reduce (1.1) to an evolution of the first order in time

$\{\begin{array}{l}{u}_t=v,\hfill \\ v_t=-A^{2}u-A^{2}v+\left(p-{\left|\nabla u\right|}^2\right)Au-b{u}^{+}-f\left(u\right)+q\left(x\right)\stackrel{\dot{}}{W},\hfill \\ u\left(x,\tau \right)=u_0\left(x\right),u_t\left(x,\tau \right)=u_1\left(x\right),x\in \Omega ,\hfill \end{array}$ (3.2)

whose equivalent Itó equation is

$\{\begin{array}{l}\text{d}u=v\text{d}t,\hfill \\ \text{d}v=-A^{2}u\text{d}t-A^{2}v\text{d}t+\left(p-{\left|\nabla u\right|}^2\right)Au\text{d}t-b{u}^{+}\text{d}t-f\left(u\right)\text{d}t+q\left(x\right)\text{d}W,\hfill \\ u\left(x,\tau \right)=u_0\left(x\right),u_t\left(x,\tau \right)=u_1\left(x\right),x\in \Omega ,\hfill \end{array}$ (3.3)

where $W\left(t\right)$ is a one-dimensional two-sided real-valued Wiener process on $\left(\Omega \mathrm{,}\mathcal{F}\mathrm{,}\mathcal{P}\mathrm{,}{\left({\theta }_t\right)}_{t\in ℝ}\right)$ . Without loss of generality,we can assume that

$\Omega =\left\{\omega \left(t\right)=W\left(t\right)\in C\left(ℝ,ℝ\right):\omega \left(0\right)=0\right\},$

where P is a Wiener measure. We can define a family of measure preserving and ergodic transformations ${\left({\theta }_t\right)}_{t\in ℝ}$ by

${\theta }_t\omega (\cdot )=\omega \left(\cdot +t\right)\omega \left(t\right)\mathrm{,}t\in ℝ\mathrm{,}\omega \in \Omega \mathrm{.}$

Let $z=v-q\left(x\right)W$ , we consider the random partial differential equation equivalent to (3.3)

$\{\begin{array}{l}\frac{\text{d}u}{\text{d}t}=z+q\left(x\right)W,\\ \frac{\text{d}z}{\text{d}t}=-A^{2}u-A^{2}z+\left(p-{\left|\nabla u\right|}^2\right)Au-b{u}^{+}-f\left(u\right)-A^{2}q\left(x\right)W,\\ u\left(x,\tau \right)=u_0\left(x\right),z\left(\tau ,\omega \right)=z\left(x,\tau ,\omega \right)=u_1\left(x\right)-q\left(x\right)W\left(\tau \right),x\in \Omega .\end{array}$ (3.4)

Apparently, there is no stochastic differential in (3.4) by comparing with stochastic differential Equation (3.3). Let

$\phi =\left(\begin{array}{c}u\\ z\end{array}\right),L=\left(\begin{array}{cc}0& I\\ -A^2& -A^2\end{array}\right),F\left(\phi ,\omega \right)=\left(\begin{array}{c}q\left(x\right)W\\ \left(p-{\left|\nabla u\right|}^2\right)Au-b{u}^{+}-f\left(u\right)-A^{2}q\left(x\right)W\end{array}\right),$

then (3.4) can be written as

$\stackrel{\dot{}}{\phi }=L\phi +F\left(\phi ,\omega \right),\phi \left(\tau ,\omega \right)={\left(u_0,z\left(\tau ,\omega \right)\right)}^T.$ (3.5)

From [13] we know that L is the infinitesimal generators of $C_0$ -semigroup ${\text{e}}^{Lt}$ on E. It is not difficult to check that the functions $F\left(\phi \mathrm{,}\omega \right)\mathrm{<mo>\; :\; </mo>}E\to E$ is locally Lipschitz continuous with respect to $\phi$ and bounded for every $\omega \in \Omega$ . By the classical semigroup theory of existence and uniqueness of solutions of evolution differential equations [13] , so we have the following theorem:

Theorem 3.1. Consider (3.5). For each $\omega \in \Omega$ and initial value $\phi \left(\tau ,\omega \right)={\left(u_0,z\left(\tau ,\omega \right)\right)}^{\text{T}}={\left(u_0,u_1\left(x\right)-q\left(x\right)W\left(\tau \right)\right)}^{\text{T}}\in E$ , there exists a unique function $\phi \left(t\mathrm{,}\omega \right)$ such that satisfies the integral equation

$\phi \left(t,\omega \right)={\text{e}}^{L\left(t-\tau \right)}\phi \left(\tau ,\omega \right)+{\displaystyle {\int }_{\tau }^t}F\left(\phi \left(s\right),\omega \right)\text{d}s,\forall \phi \left(\tau ,\omega \right)\in E.$

By theorem 3.1, we can prove that for $P-a\mathrm{.}s\mathrm{.}$ every $\omega \in \Omega$ the following statements hold for all $T>0$ :

1) If $\phi \left(\tau ,\omega \right)\in E$ , then $\phi \left(t\mathrm{,}\omega \right)\in C\left(\left[\tau \mathrm{,}\tau +T\right)\mathrm{<mo>\; ;\; </mo>}D\left(A\right)\right)\times C\left(\left[\tau \mathrm{,}\tau +T\right)\mathrm{<mo>\; ;\; </mo>}H\right)$ .

2) $\phi \left(t\mathrm{,}\phi \left(\tau ,\omega \right)\right)$ is continuous in t and $\phi \left(\tau ,\omega \right)$ .

3) The solution mapping of (3.5) satisfies the properties of RDS.

Equation (3.5) has a unique solution for every $\omega \in \Omega$ . Hence the solution mapping

$\tilde{S}\left(t\mathrm{,}\omega \right)\mathrm{<mo>\; :\; </mo>}\phi \left(\tau \mathrm{,}\omega \right)\mapsto \phi \left(t\mathrm{,}\omega \right)$ (3.6)

generates a random dynamical system, so the transformation

$S\left(t,\omega \right):\phi \left(\tau ,\omega \right)+{\left(0,q\left(x\right)W\left(\tau \right)\right)}^{\text{T}}\mapsto \phi \left(t,\omega \right)+{\left(0,q\left(x\right)W\left(t\right)\right)}^{\text{T}}$ (3.7)

also determines a random dynamical system corresponding to Equation (3.2).

4. Existence of a Random Attractor

In this section,we prove the existence of a random attractor for RDS (3.7) in E. Let $\tilde{z}=z+\epsilon u,\psi ={\left(u,\tilde{z}\right)}^{\text{T}}$ , where

$\epsilon =\frac{{\lambda }^2}{4{\lambda }^2+3\lambda +4}.$ (4.1)

So Equation (3.4) can be written as

$\stackrel{\dot{}}{\psi }+Q\psi =\tilde{F}\left(\psi ,\omega \right),\psi \left(\tau ,\omega \right)={\left(u_0,z\left(\tau ,\omega \right)+\epsilon u_0\right)}^{\text{T}},t\ge \tau .$ (4.2)

where

$Q=\left(\begin{array}{cc}\epsilon I& -I\\ \left(1-\epsilon \right)A^2+{\epsilon }^{2}I& A^2-\epsilon I\end{array}\right),$

$\tilde{F}\left(\psi ,\omega \right)=\left(\begin{array}{c}q\left(x\right)W\\ \left(p-{\left|\nabla u\right|}^2\right)Au-b{u}^{+}-f\left(u\right)+\left(\epsilon -A^2\right)q\left(x\right)W\end{array}\right).$

The mapping

${\tilde{S}}_{\epsilon }\left(t\mathrm{,}\omega \right)\mathrm{<mo>\; :\; </mo>}{\left(u_0\mathrm{,}z\left(\tau \mathrm{,}\omega \right)+\epsilon u_0\right)}^{\text{T}}\mapsto {\left(u\left(t\right)\mathrm{,}z\left(t\right)+\epsilon u\left(t\right)\right)}^{\text{T}}\mathrm{,}E\to E\mathrm{,}t\ge \tau$

is defined by (4.2).

To show the conjugation of the solution of the stochastic partial differential Equation (1) and the random partial differential Equation (4.2), we introduce the homeomorphism

$R_{\epsilon }\mathrm{<mo>\; :\; </mo>}{\left(u\mathrm{,}z\right)}^{\text{T}}\mapsto {\left(u\mathrm{,}z+\epsilon u\right)}^{\text{T}}$

with the inverse homeomorphism $R_{-\epsilon }$ . Then the transformation

${\tilde{S}}_{\epsilon }\left(t\mathrm{,}\omega \right)=R_{\epsilon }S\left(t\mathrm{,}\omega \right)R_{-\epsilon }$ (4.3)

also determines RDS corresponding to Equation (1). Therefore, for RDS (7) we only need consider the equivalent random dynamical system $S_{\epsilon }\left(t\mathrm{,}\omega \right)=R_{\epsilon }S\left(t\mathrm{,}\omega \right)R_{-\epsilon }$ , where $S_{\epsilon }\left(t\mathrm{,}\omega \right)$ is decided by

$\stackrel{\dot{}}{\xi }+Q\xi =G\left(\xi ,\omega \right),\xi \left(\tau ,\omega \right)={\left(u_0,u_1+\epsilon u_0\right)}^{\text{T}},t\ge \tau ,$ (4.4)

where

$\xi =\left(\begin{array}{c}u\left(t\right)\\ u_t\left(t\right)+\epsilon u\left(t\right)\end{array}\right),G\left(\xi ,\omega \right)=\left(\begin{array}{c}0\\ \left(p-{\left|\nabla u\right|}^2\right)Au-b{u}^{+}-f\left(u\right)+q\left(x\right)\stackrel{\dot{}}{W}\end{array}\right).$

Next, we prove a positivity property of the operator Q in E that plays a vital role throughout the paper.

Lemma 4.1. For any $\varphi ={\left(u\mathrm{,}z\right)}^{\text{T}}\in E$ , there holds

${\left(Q\varphi \mathrm{,}\varphi \right)}_E\ge \frac{\epsilon }{2}{\Vert \varphi \Vert }_E^2+\frac{\epsilon }{4}{\Vert u\Vert }_2^2+\frac{\lambda }{2}{\Vert z\Vert }^2\mathrm{.}$

Proof. Since $Q\varphi ={\left(\epsilon u-z\mathrm{,}\left(1-\epsilon \right)A^{2}u+{\epsilon }^{2}u+A^{2}z-\epsilon z\right)}^{\text{T}}$ , by using the Poincaré inequality and the Young inequality, we conclude that

$\begin{array}{c}{\left(Q\varphi ,\varphi \right)}_E=\epsilon {\Vert u\Vert }_2^2-\epsilon \left(Au,Az\right)+{\epsilon }^2\left(u,z\right)+{\Vert Az\Vert }^2-\epsilon {\Vert z\Vert }^2\\ \ge \epsilon {\Vert u\Vert }_2^2-\frac{\epsilon }{8}{\Vert u\Vert }_2^2-2\epsilon {\Vert Az\Vert }^2-\frac{\epsilon }{8}{\Vert u\Vert }_2^2-\frac{2{\epsilon }^3}{\lambda }{\Vert z\Vert }^2+{\Vert Az\Vert }^2-\epsilon {\Vert z\Vert }^2\\ \ge \frac{\epsilon }{2}{\Vert \varphi \Vert }_E^2+\frac{\epsilon }{4}{\Vert u\Vert }_2^2+\left(1-2\epsilon \right)\lambda {\Vert z\Vert }^2-\left(\frac{2\epsilon }{\lambda }+\frac{3\epsilon }{2}\right){\Vert z\Vert }^2\\ =\frac{\epsilon }{2}{\Vert \varphi \Vert }_E^2+\frac{\epsilon }{4}{\Vert u\Vert }_2^2+\frac{\lambda }{2}{\Vert z\Vert }^2.\end{array}$

where $\epsilon =\frac{{\lambda }^2}{4{\lambda }^2+3\lambda +4}$ . ,

Lemma 4.2. Let (1.2)-(1.4) hold, there exist a random variable $r_1\left(\omega \right)>0$ , and a bounded ball $B_0$ of $E$ centered at 0 with random radius $r_0\left(\omega \right)>0$ such that for any bounded non-random set $B$ of $E$ , there exists a deterministic $T\left(B\right)\le -1$ such that the solution $\psi \left(t\mathrm{,}\omega \mathrm{,}\psi \left(\tau \mathrm{,}\omega \right)\right)={\left(u\left(t\mathrm{,}\omega \right)\mathrm{,}\tilde{z}\left(t\mathrm{,}\omega \right)\right)}^{\text{T}}$ of (4.2) with initial value ${\left(u_0\mathrm{,}{u}_1+\epsilon u_0\right)}^{\text{T}}\in B$ satisfies for $P-a\mathrm{.}s\mathrm{.}\omega \in \Omega$ ,

${\Vert \psi \left(-\mathrm{1,}\omega \mathrm{,}\psi \left(\tau \mathrm{,}\omega \right)\right)\Vert }_E\le r_0\left(\omega \right)\mathrm{,}\tau \le T\left(B\right)\mathrm{,}$

and for all $\tau \le t\le 0$

${\Vert \psi \left(t\mathrm{,}\omega \mathrm{,}\psi \left(\tau \mathrm{,}\omega \right)\right)\Vert }_E^2\le R\left(\tau \mathrm{,}\omega \right)\mathrm{,}$ (4.5)

where $\tilde{z}\left(t\omega \right)=u_t\left(t\right)+\epsilon u\left(t\right)-q\left(x\right)W\left(t\right)$ , and $R\left(\tau \mathrm{,}\omega \right)$ is given by

$\begin{array}{l}R\left(\tau ,\omega \right)=2{\text{e}}^{-{\epsilon }_1\left(t-\tau \right)}\left({\Vert u_0\Vert }_2^2+{\Vert u_1+\epsilon u_0\Vert }^2+{\Vert q\Vert }^2{\left|W\right|}^2+{\displaystyle {\int }_{\Omega }}F\left(u_0\right)\text{d}x+{\left|\nabla u_0\right|}^2-p{)}^2\right)\\ +r_1^2\left(\omega \right).\end{array}$

Besides it is easy to deduce a similar absorption result for

$\varphi \left(-1\right)={\left(u\left(-1\right)\mathrm{,}{u}_t\left(-1\right)+\epsilon u\left(-1\right)\right)}^{\text{T}}$

instead of $\psi \left(-1\right)$ .

Proof. We take the inner product in E of (4.2) with $\psi ={\left(u\mathrm{,}\tilde{z}\right)}^{\text{T}}$ , where $\tilde{z}=u_t+\epsilon u-q\left(x\right)W$ , we get

$\frac{1}{2}\frac{\text{d}}{\text{d}t}{\Vert \psi \Vert }_E^2+{\left(Q\psi \mathrm{,}\psi \right)}_E={\left(\tilde{F}\left(\psi \mathrm{,}\omega \right)\mathrm{,}\psi \right)}_E\mathrm{,}\forall t\ge \tau \mathrm{,}$ (4.6)

where

$\begin{array}{c}{\left(\tilde{F}\left(\psi ,\omega \right),\psi \right)}_E=\left(\left(u,q\left(x\right)W\right)\right)-b\left(u^{+},\tilde{z}\right)-\left(f\left(u\right),\tilde{z}\right)\\ -\left(A^{2}q\left(x\right)W,\tilde{z}\right)+\epsilon \left(q\left(x\right)W,\tilde{z}\right)+\left(\left(p-{\left|\nabla u\right|}^2\right)Au,\tilde{z}\right).\end{array}$ (4.7)

We deal with the terms in (4.7) one by one as follows:

$\left(\left(u\mathrm{,}q\left(x\right)W\right)\right)\le \frac{\epsilon }{4}{\Vert u\Vert }_2^2+\frac{1}{\epsilon }{\Vert q\Vert }_2^2{\left|W\left(t\right)\right|}^2\mathrm{<mo>\; ;\; </mo>}$ (4.8)

$\begin{array}{c}-b\left(u^{+},\tilde{z}\right)=-b\left(u^{+},u_t+\epsilon u-q\left(x\right)W\right)\\ =-\frac{1}{2}\frac{\text{d}}{\text{d}t}b{\Vert u^{+}\Vert }^2-\epsilon b{\Vert u^{+}\Vert }^2+b\left(u^{+},q\left(x\right)W\right)\\ \le -\frac{1}{2}\frac{\text{d}}{\text{d}t}b{\Vert u^{+}\Vert }^2-\frac{\epsilon b}{2}{\Vert u^{+}\Vert }^2+\frac{b{\Vert q\Vert }^2}{2\epsilon }{\left|W\right|}^2;\end{array}$ (4.9)

$\epsilon \left(q\left(x\right)W\mathrm{,}\tilde{z}\right)\le \frac{\lambda }{4}{\Vert \tilde{z}\Vert }^2+\frac{{\epsilon }^2}{\lambda }{\Vert q\Vert }^2{\left|W\right|}^2\mathrm{<mo>\; ;\; </mo>}$ (4.10)

$\left|-\left(A^{2}q\left(x\right)W\mathrm{,}\tilde{z}\right)\right|\le \frac{{\Vert q\Vert }_4^2}{\lambda }{\left|W\right|}^2+\frac{\lambda }{4}{\Vert \tilde{z}\Vert }^2.$ (4.11)

By using (1.2)-(1.3) and the Hölder inequality, we get

$\begin{array}{c}\left(f\left(u\right),q\left(x\right)W\right)\le C_0{\displaystyle {\int }_{\Omega }}\left(1+{\left|u\right|}^k\right)q\left(x\right)W\left(t\right)\text{d}x\\ \le C_0\Vert q\Vert \left|W\right|+C_0{\left({\displaystyle {\int }_{\Omega }}{\left|u\right|}^{k+1}\text{d}x\right)}^{\frac{k}{k+1}}{\Vert q\Vert }_{k+1}\left|W\right|\\ \le C_0\Vert q\Vert \left|W\right|+C_0{C}_1^{-\frac{k}{k+1}}{\left({\displaystyle {\int }_{\Omega }}\left(F\left(u\right)+C_1\right)\text{d}x\right)}^{\frac{k}{k+1}}{\Vert q\Vert }_{k+1}\left|W\right|\\ \le C_0\Vert q\Vert \left|W\right|+\frac{\epsilon C_0{C}_1^{-1}}{2}{\displaystyle {\int }_{\Omega }}F\left(u\right)\text{d}x+\frac{C_0}{2\epsilon }{\Vert q\Vert }_{k+1}^{k+1}{\left|W\right|}^{k+1}+\frac{\epsilon C_0\left|\Omega \right|}{2}.\end{array}$ (4.12)

Inequality (4.12) together with (1.4) yields

$\begin{array}{c}-\left(f\left(u\right),\tilde{z}\right)=-\left(f\left(u\right),u_t+\epsilon u-q\left(x\right)W\right)\\ \le -\frac{\text{d}}{\text{d}t}{\displaystyle {\int }_{\Omega }}F\left(u\right)\text{d}x-\epsilon C_2{\displaystyle {\int }_{\Omega }}F\left(u\right)\text{d}x+\epsilon C_2\left|\Omega \right|+\left(f\left(u\right),q\left(x\right)W\right)\\ \le -\frac{\text{d}}{\text{d}t}{\displaystyle {\int }_{\Omega }}F\left(u\right)\text{d}x-\frac{\epsilon \left(2C_2-C_0{C}_1^{-1}\right)}{2}{\displaystyle {\int }_{\Omega }}F\left(u\right)\text{d}x+C_0\Vert q\Vert \left|W\right|\\ +\frac{C_0}{2\epsilon }{\Vert q\Vert }_{k+1}^{k+1}{\left|W\right|}^{k+1}+\frac{\epsilon \left(C_0+2C_2\right)}{2}\left|\Omega \right|.\end{array}$ (4.13)

$\begin{array}{c}\left(\left(p-{\left|\nabla u\right|}^2\right)Au,\tilde{z}\right)=\left(p-{\left|\nabla u\right|}^2\right)\left(Au,u_t+\epsilon u-q\left(x\right)W\right)\\ =\left(p-{\left|\nabla u\right|}^2\right)\left[\frac{1}{2}\frac{\text{d}{\left|\nabla u\right|}^2}{\text{d}t}+\epsilon {\left|\nabla u\right|}^2-\left(Au,q\left(x\right)W\right)\right]\\ =-\frac{1}{4}\frac{\text{d}}{\text{d}t}{\left({\left|\nabla u\right|}^2-p\right)}^2-\epsilon {\left({\left|\nabla u\right|}^2-p\right)}^2-\epsilon p\left({\left|\nabla u\right|}^2-p\right)\\ +\left({\left|\nabla u\right|}^2-p\right)\left(Au,q\left(x\right)W\right).\end{array}$ (4.14)

$\left({\left|\nabla u\right|}^2-p\right)\left(Au\mathrm{,}q\left(x\right)W\right)\le \frac{\epsilon }{4}{\left({\left|\nabla u\right|}^2-p\right)}^2+\frac{\epsilon }{2}{\left|\nabla u\right|}^4+\frac{1}{2{\epsilon }^3}{\left|\nabla q\right|}^4{\left|W\right|}^4\mathrm{.}$ (4.15)

Collecting with (4.6)-(4.15) and Lemma 4.1, we get that

$\begin{array}{l}\frac{\text{d}}{\text{d}t}\left[{\Vert \psi \Vert }_E^2+b{\Vert u^{+}\Vert }^2+2{\displaystyle {\int }_{\Omega }}F\left(u\right)\text{d}x+\frac{1}{2}{\left({\left|\nabla u\right|}^2-p\right)}^2+2C_1\left|\Omega \right|\right]\\ +{\epsilon }_1\left[{\Vert \psi \Vert }_E^2+b{\Vert u^{+}\Vert }^2+2{\displaystyle {\int }_{\Omega }}F\left(u\right)\text{d}x+\frac{1}{2}{\left({\left|\nabla u\right|}^2-p\right)}^2+2C_1\left|\Omega \right|\right]\\ \le M\left(1+\left|W\left(t\right)\right|+{\left|W\left(t\right)\right|}^2+{\left|W\left(t\right)\right|}^4+{\left|W\left(t\right)\right|}^{k+1}\right),\end{array}$

where ${\epsilon }_1=\min \left\{\epsilon ,\frac{2C_2-C_0{C}_1^{-1}}{2}\epsilon \right\}$ , $C_2>\frac{C_0{C}_1^{-1}}{2}$ , and

$\begin{array}{l}M=\max \{\epsilon p^2+\epsilon \left(C_0+2C_2+2C_1\right)\left|\Omega \right|,C_0\Vert q\Vert ,\begin{array}{c}\\ \end{array}\\ \frac{{\Vert q\Vert }_2^2}{\epsilon }+\left(\frac{{\epsilon }^2}{\lambda }+\frac{b}{2\epsilon }\right){\Vert q\Vert }^2+\frac{{\Vert q\Vert }_4^2}{\lambda },\frac{C_0{\Vert q\Vert }_{k+1}^{k+1}}{2\epsilon },\frac{{\Vert \nabla q\Vert }^4}{2{\epsilon }^3}\}.\end{array}$

By the Gronwall lemma, we conclude that

$\begin{array}{l}{\Vert \psi \left(t\mathrm{,}\omega \mathrm{,}\psi \left(\tau \mathrm{,}\omega \right)\right)\Vert }_E^2\\ \le {\text{e}}^{-{\epsilon }_1\left(t-\tau \right)}\left({\Vert \psi \left(\tau \mathrm{,}\omega \right)\Vert }_E^2+b{\Vert u_0^{+}\Vert }^2+2{\displaystyle {\int }_{\Omega }}F\left(u_0\right)\text{d}x+\frac{1}{2}{\left({\left|\nabla u_0\right|}^2-p\right)}^2+2C_1\left|\Omega \right|\right)\\ +M{\displaystyle {\int }_{\tau }^t}{\text{e}}^{-{\epsilon }_1\left(t-s\right)}\left(1+\left|W\left(s\right)\right|+{\left|W\left(s\right)\right|}^2+{\left|W\left(s\right)\right|}^4+{\left|W\left(s\right)\right|}^{k+1}\right)\text{d}s\end{array}$

$\begin{array}{l}\le 2{\text{e}}^{-{\epsilon }_1\left(t-\tau \right)}[{\Vert u_0\Vert }_2^2+{\Vert u_1+\epsilon u_0\Vert }^2+{\Vert q\Vert }^2{\left|W\left(\tau \right)\right|}^2+b{\Vert u_0^{+}\Vert }^2\begin{array}{c}\\ \end{array}\\ +{\displaystyle {\int }_{\Omega }}F\left(u_0\right)\text{d}x+\frac{1}{4}{\left({\left|\nabla u_0\right|}^2-p\right)}^2]\\ +M{\displaystyle {\int }_{\tau }^t}{\text{e}}^{-{\epsilon }_1\left(t-s\right)}\left(1+\left|W\left(s\right)\right|+{\left|W\left(s\right)\right|}^2+{\left|W\left(s\right)\right|}^4+{\left|W\left(s\right)\right|}^{k+1}\right)\text{d}s.\end{array}$ (4.16)

Let

$\begin{array}{c}{r}_0^2\left(\omega \right)=2\left(1+\underset{\tau \le -1}{sup}{\text{e}}^{{\epsilon }_1\tau }{\Vert q\Vert }^2{\left|W\left(\tau \right)\right|}^2\right)+\frac{M}{{\epsilon }_1}\\ +M{\displaystyle {\int }_{-\infty }^{-1}}{\text{e}}^{-{\epsilon }_1\left(-1-s\right)}\left(\left|W\left(s\right)\right|+{\left|W\left(s\right)\right|}^2+{\left|W\left(s\right)\right|}^4+{\left|W\left(s\right)\right|}^{k+1}\right)\text{d}s\mathrm{,}\end{array}$

$r_1^2\left(\omega \right)=\frac{M}{{\epsilon }_1}+M{\displaystyle {\int }_{-\infty }^0}{\text{e}}^{{\epsilon }_{1}s}\left(\left|W\left(s\right)\right|+{\left|W\left(s\right)\right|}^2+{\left|W\left(s\right)\right|}^4+{\left|W\left(s\right)\right|}^{k+1}\right)\text{d}s,$

where $r_0^2\left(\omega \right)$ and $r_1^2\left(\omega \right)$ are finite $P-a\mathrm{.}s\mathrm{.}$ , we get a bounded set B of E, we choose $T\left(B\right)\le -1$ such that

${\text{e}}^{-{\epsilon }_1\left(-1-\tau \right)}\left({\Vert u_0\Vert }_2^2+{\Vert u_1+\epsilon u_0\Vert }^2+b{\Vert u_0^{+}\Vert }^2+{\displaystyle {\int }_{\Omega }}F\left(u_0\right)\text{d}x+C_1\left|\Omega \right|+\frac{1}{4}{\left({\left|\nabla u_0\right|}^2-p\right)}^2\right)\le 1$ (4.17)

for all $\forall {\left(u_0\mathrm{,}{u}_1+\epsilon u_0\right)}^{\text{T}}\in B$ , and

${\text{e}}^{{\epsilon }_1\tau }\left({\Vert u_0\Vert }_2^2+{\Vert u_1+\epsilon u_0\Vert }^2+b{\Vert u_0^{+}\Vert }^2+{\displaystyle {\int }_{\Omega }}F\left(u_0\right)\text{d}x+C_1\left|\Omega \right|+\frac{1}{4}{\left({\left|\nabla u_0\right|}^2-p\right)}^2\right)\le 1$ (4.18)

for all $\forall {\left(u_0\mathrm{,}{u}_1+\epsilon u_0\right)}^{\text{T}}\in B$ , and for all $\tau \le T\left(B\right)$ . ,

Let $u\left(t\right)$ be a solution of problem (1.1) with initial value ${\left(u_0\mathrm{,}{u}_t+\epsilon u_0\right)}^{\text{T}}\in B$ . we make the decomposition $u\left(t\right)=y_1\left(t\right)+y_2\left(t\right)$ , where $y_1\left(t\right)$ and $y_2\left(t\right)$ satisfy

$\{\begin{array}{ll}{y}_{1tt}+{\Delta }^2{y}_1+{\Delta }^2{y}_{1t}+\left(p-{\left|\nabla u\right|}^2\right)\Delta y_1=0,\hfill & \Omega \times \left[\tau ,+\infty \right),\tau \in ℝ,\hfill \\ y_1\left(x,t\right)=\nabla y_1\left(x,t\right)=0,\hfill & x\in \partial \Omega ,t\ge \tau ,\hfill \\ y_1\left(x,\tau \right)=u_0\left(x\right),y_{1t}\left(x,\tau \right)=u_1\left(x\right),\hfill & x\in \Omega .\hfill \end{array}$ (4.19)

and

$\{\begin{array}{ll}{y}_{2tt}+{\Delta }^2{y}_2+{\Delta }^2{y}_t+\left(p-{\left|\nabla u\right|}^2\right)\Delta y_2+b{u}^{+}+f\left(u\right)=q\left(x\right)\stackrel{\dot{}}{W},\hfill & \Omega \times \left[\tau ,+\infty \right),\tau \in ℝ,\hfill \\ y_2\left(x,t\right)=\nabla y_2\left(x,t\right)=0,\hfill & x\in \partial \Omega ,t\ge \tau ,\hfill \\ y_2\left(x,\tau \right)=0,y_{2t}\left(x,\tau \right)=0,\hfill & x\in \Omega .\hfill \end{array}$ (4.20)

Lemma 4.3. Let $p<\frac{\sqrt{\lambda }}{4}$ , $B$ be a bounded non-random subset of $E$ , $\forall {\left(u_0\mathrm{,}{u}_1+\epsilon u_0\right)}^{\text{T}}\in B$

$\begin{array}{c}{\Vert Y_1\left(0\right)\Vert }_E^2={\Vert y_1\left(0\right)\Vert }_2^2+{\Vert y_{1t}\left(0\right)+\epsilon y_1\left(0\right)\Vert }^2\\ \le \frac{{\text{e}}^{\epsilon \tau }}{\left(1-\epsilon \right)C\left(p\right)}\left({\Vert u_0\Vert }_2^2-\left(p-{\left|\nabla u_0\right|}^2\right){\left|\nabla u_0\right|}^2+{\Vert u_1+\epsilon u_0\Vert }^2\right)\mathrm{,}\end{array}$ (4.21)

where $Y_1={\left(y_1\mathrm{,}{y}_{1t}+\epsilon y_1\right)}^{\text{T}}$ satisfies (4.19), $C\left(p\right)=\{\begin{array}{ll}0,\hfill & p\le 0,\hfill \\ 1-\frac{p}{\sqrt{\lambda }},\hfill & p<\frac{\sqrt{\lambda }}{4}.\hfill \end{array}$

Proof. Let $v=y_{1t}+\epsilon y_1$ , taking the scalar product in H of (4.19), we get that

$\begin{array}{l}\frac{1}{2}\frac{\text{d}}{\text{d}t}\left[{\Vert v\Vert }^2+\left(1-\epsilon \right){\Vert y_1\Vert }_2^2+\left({\left|\nabla u_0\right|}^2-p\right){\left|\nabla y_1\right|}^2\right]\\ +{\Vert v\Vert }_2^2-\epsilon {\Vert v\Vert }^2+\epsilon \left(1-\epsilon \right){\Vert y_1\Vert }_2^2+\epsilon \left({\left|\nabla u_0\right|}^2-p\right){\left|\nabla y_1\right|}^2+{\epsilon }^2\left(y_1\mathrm{,}v\right)=\mathrm{0,}\end{array}$ (4.22)

using the Hölder inequality and the Young inequality, we get that

$\begin{array}{l}{\Vert v\Vert }_2^2-\epsilon {\Vert v\Vert }^2+\epsilon \left(1-\epsilon \right){\Vert y_1\Vert }_2^2+{\epsilon }^2\left(y_1\mathrm{,}v\right)\\ \ge \epsilon \left(1-\epsilon \right){\Vert y_1\Vert }_2^2+{\Vert v\Vert }_2^2-\epsilon {\Vert v\Vert }^2-\frac{\epsilon \left(1-\epsilon \right)}{2}{\Vert y_1\Vert }_2^2-\frac{{\epsilon }^3}{2\left(1-\epsilon \right)\lambda }\\ \ge \frac{\epsilon \left(1-\epsilon \right)}{2}{\Vert y_1\Vert }_2^2+\left(\lambda -\epsilon -\frac{\epsilon }{2\lambda }\right){\Vert v\Vert }^2\ge \frac{\epsilon \left(1-\epsilon \right)}{2}{\Vert y_1\Vert }_2^2+\frac{\epsilon }{2}{\Vert v\Vert }^2\mathrm{,}\end{array}$ (4.23)

we have that

$\begin{array}{l}\frac{\text{d}}{\text{d}t}\left[{\Vert v\Vert }^2+\left(1-\epsilon \right){\Vert y_1\Vert }_2^2+\left({\left|\nabla u_0\right|}^2-p\right){\left|\nabla y_1\right|}^2\right]\\ +\epsilon \left[{\Vert v\Vert }^2+\left(1-\epsilon \right){\Vert y_1\Vert }_2^2+2\left({\left|\nabla u_0\right|}^2-p\right){\left|\nabla y_1\right|}^2\right]\le 0.\end{array}$

Let $\vartheta \left(t\right)={\Vert v\Vert }^2+\left(1-\epsilon \right){\Vert y_1\Vert }_2^2+\left({\left|\nabla u_0\right|}^2-p\right){\left|\nabla y_1\right|}^2$ , we can get that

$\vartheta \left(t\right)\ge \left(1-\epsilon \right)C\left(p\right){\Vert y_1\Vert }_2^2+{\Vert v\Vert }^2+{\left|\nabla u\right|}^2{\left|\nabla y_1\right|}^2>\left(1-\epsilon \right)C\left(p\right){\Vert Y_1\Vert }_E^2>0$

and for $p<\frac{\sqrt{\lambda }}{4}$ we get

$\epsilon \left[{\Vert v\Vert }^2+\left(1-\epsilon \right){\Vert y_1\Vert }_2^2+2\left({\left|\nabla u_0\right|}^2-p\right){\left|\nabla y_1\right|}^2\right]\ge \frac{\epsilon }{2}\vartheta \left(t\right)\mathrm{,}$

so

$\frac{\text{d}}{\text{d}t}\vartheta \left(t\right)+\frac{\epsilon }{2}\vartheta \left(t\right)\le \mathrm{0,}$

By the Gronwall lemma, getting that (4.21). ,

Lemma 4.4. Assume that (1.2) holds, there exists a random radius $r_2\left(\omega \right)$ , such that for $P-a\mathrm{.}s\mathrm{.}\omega \in \Omega$ ,

${\Vert A^{\frac{1}{2}}{Y}_2\left(\mathrm{0,}\omega \mathrm{<mo>\; ;\; </mo>}{Y}_2\left(\tau \mathrm{,}\omega \right)\right)\Vert }_E^2\le r_2^2\left(\omega \right)\mathrm{,}$ (4.24)

where $Y_2={\left(y_2\mathrm{,}{y}_{2t}+\epsilon y_2-q\left(x\right)W\right)}^{\text{T}}$ satisfies (4.20).

Proof. Let $Y_2={\left(y_2\mathrm{,}{y}_{2t}+\epsilon y_2-q\left(x\right)W\right)}^{\text{T}}$ , Equation (4.20) can be written to

$Y_{2t}+Q{Y}_2=H\left(Y_2,\omega \right),Y_2\left(\tau \right)={\left(0,-q\left(x\right)W\left(\tau \right)\right)}^{\text{T}},$ (4.25)

where

$H\left(Y_2\mathrm{,}\omega \right)=\left(\begin{array}{c}q\left(x\right)W\\ \left(p-{\left|\nabla u\right|}^2\right)A{y}_2-b{u}^{+}-f\left(u\right)+\left(\epsilon -A^2\right)q\left(x\right)W\end{array}\right)\mathrm{.}$

Taking the scalar product in $E$ of (4.25) with $A{Y}_2$ , we get that

$\frac{1}{2}\frac{\text{d}}{\text{d}t}{\Vert A^{\frac{1}{2}}{Y}_2\Vert }_E^2+{\left(Q{Y}_2\mathrm{,}A{Y}_2\right)}_E={\left(H\left(Y_2\mathrm{,}\omega \right)\mathrm{,}A{Y}_2\right)}_E\mathrm{,}$ (4.26)

where

$\begin{array}{l}{\left(H\left(Y_2,\omega \right),A{Y}_2\right)}_E\\ =\left(\left(A{y}_2,q\left(x\right)W\right)\right)+\left(\left(p-{\left|\nabla u\right|}^2\right)A{y}_2,A\left(y_{2t}+\epsilon y_2-q\left(x\right)W\right)\right)\\ +\left(-b{u}^{+}-f\left(u\right)+\left(\epsilon -A^2\right)q\left(x\right)W,A\left(y_{2t}+\epsilon y_2-q\left(x\right)W\right)\right).\end{array}$ (4.27)

Due to Lemma 4.1, we get that

${\left(Q{Y}_2\mathrm{,}A{Y}_2\right)}_E\ge \frac{\epsilon }{2}{\Vert A^{\frac{1}{2}}{Y}_2\Vert }_E^2+\frac{\epsilon }{4}{\Vert A^{\frac{1}{2}}{y}_2\Vert }_2^2+\frac{\lambda }{2}{\Vert A^{\frac{1}{2}}\left(y_{2t}+\epsilon y_2-q\left(x\right)W\right)\Vert }^2\mathrm{.}$ (4.28)

Using the Young inequality, we have that

$\left(\left(A{y}_2\mathrm{,}q\left(x\right)W\right)\right)\le \frac{\epsilon }{4}{\Vert A^{\frac{1}{2}}{y}_2\Vert }_2^2+\frac{1}{\epsilon }{\Vert A^{\frac{1}{2}}q\Vert }_2^2{\left|W\right|}^2\mathrm{<mo>\; ;\; </mo>}$ (4.29)

$\begin{array}{l}\left|-\left(b{u}^{+}\mathrm{,}A\left(y_{2t}+\epsilon y_2-q\left(x\right)W\right)\right)\right|\\ \le \frac{2b^2}{\lambda }{\Vert A^{\frac{1}{2}}{u}^{+}\Vert }_2^2+\frac{\lambda }{8}{\Vert A^{\frac{1}{2}}\left(y_{2t}+\epsilon y_2-q\left(x\right)W\right)\Vert }^2\mathrm{<mo>\; ;\; </mo>}\end{array}$ (4.30)

$\begin{array}{l}\left|-\left(A^{2}q\left(x\right)W\mathrm{,}A\left(y_{2t}+\epsilon y_2-q\left(x\right)W\right)\right)\right|\\ \le \frac{2}{\lambda }{\Vert A^{\frac{5}{2}}q\Vert }^2{\left|W\right|}^2+\frac{\lambda }{8}{\Vert A^{\frac{1}{2}}\left(y_{2t}+\epsilon y_2-q\left(x\right)W\right)\Vert }^2\mathrm{<mo>\; ;\; </mo>}\end{array}$ (4.31)

$\begin{array}{l}\left(\epsilon q\left(x\right)W\mathrm{,}A\left(y_{2t}+\epsilon y_2-q\left(x\right)W\right)\right)\\ \le \frac{2{\epsilon }^2}{\lambda }{\Vert A^{\frac{1}{2}}q\Vert }^2{\left|W\right|}^2+\frac{\lambda }{8}{\Vert A^{\frac{1}{2}}\left(y_{2t}+\epsilon y_2-q\left(x\right)W\right)\Vert }^2.\end{array}$ (4.32)

By (1.2), (4.5) and Sobolev embedding theorem, we obtain that $f\left(u\right)$ is uniformly bounded in $L^{\infty }$ , that is, there exists a constant $M>0$ such that

${\left|f^{\prime }\left(s\right)\right|}_{L^{\infty }}\le M.$ (4.33)

Combining with (4.33), the Sobolev embedding theorem and the Young inequality, we have that

$\begin{array}{l}\left|-\left(f\left(u\right)\mathrm{,}A\left(y_{2t}+\epsilon y_2-q\left(x\right)W\right)\right)\right|\\ \le \Vert A^{\frac{1}{2}}f\left(u\right)\Vert \Vert A^{\frac{1}{2}}\left(y_{2t}+\epsilon y_2-q\left(x\right)W\right)\Vert \\ \le \frac{2\mu M^2}{\lambda }{\Vert u\Vert }_2^2+\frac{\lambda }{8}{\Vert A^{\frac{1}{2}}\left(y_{2t}+\epsilon y_2-q\left(x\right)W\right)\Vert }^2\mathrm{,}\end{array}$ (4.34)

where $\mu$ is a positive constant.

$\begin{array}{l}-\left(\left(p-{\left|\nabla u\right|}^2\right)A{y}_2,A\left(y_{2t}+\epsilon y_2-q\left(x\right)W\right)\right)\\ =\frac{1}{2}\frac{\text{d}}{\text{d}t}\left({\left|\nabla u\right|}^2-p\right){\Vert y_2\Vert }_2^2+\epsilon \left({\left|\nabla u\right|}^2-p\right){\Vert y_2\Vert }_2^2+\left(p-{\left|\nabla u\right|}^2\right)\left(A{y}_2,q\left(x\right)W\right)\\ \ge \frac{1}{2}\frac{\text{d}}{\text{d}t}\left({\left|\nabla u\right|}^2-p\right){\Vert y_2\Vert }_2^2+\epsilon \left({\left|\nabla u\right|}^2-p\right){\Vert y_2\Vert }_2^2-\left(\left|p\right|+{\left|\nabla u\right|}^2\right){\Vert y_2\Vert }_2\Vert Aq\Vert \left|W\right|\end{array}$

$\begin{array}{l}\ge \frac{1}{2}\frac{\text{d}}{\text{d}t}\left({\left|\nabla u\right|}^2-p\right){\Vert y_2\Vert }_2^2+\frac{\epsilon }{2}\left({\left|\nabla u\right|}^2-p\right){\Vert y_2\Vert }_2^2-\frac{\epsilon }{2}\left|p\right|{\Vert y_2\Vert }_2^2\\ -\frac{1}{2}{\Vert Aq\Vert }^2{\left|W\right|}^2-\frac{1}{2}{\left(\left|p\right|+{\left|\nabla u\right|}^2\right)}^2{\Vert y_2\Vert }_2^2\\ =\frac{1}{2}\frac{\text{d}}{\text{d}t}\left({\left|\nabla u\right|}^2-p\right){\Vert y_2\Vert }_2^2+\frac{\epsilon }{2}\left({\left|\nabla u\right|}^2-p\right){\Vert y_2\Vert }_2^2\\ -\frac{1}{2}\left(\epsilon \left|p\right|+{\left(\left|p\right|+{\left|\nabla u\right|}^2\right)}^2\right){\Vert y_2\Vert }_2^2-\frac{1}{2}{\Vert Aq\Vert }^2{\left|W\right|}^2\end{array}$ (4.35)

Let $E\left(t\right)={\Vert A^{\frac{1}{2}}{Y}_2\Vert }_E^2+\left({\left|\nabla u\right|}^2-p\right){\Vert y_2\Vert }_2^2$ , by the Poincaré and $C\left(p\right)$ , we get

that, $E\left(t\right)\ge C\left(p\right){\Vert A^{\frac{1}{2}}{Y}_2\Vert }_E^2>0$ . Using (4.27)-(4.35) and (4.5), for $\tau \le T\left(\omega \right)$ , from (4.27) we get that

$\begin{array}{l}\frac{\text{d}}{\text{d}t}E\left(t\right)+\epsilon E\left(t\right)\\ \le C\left(\epsilon \left|p\right|+{\left(\left|p\right|+R\left(s,\omega \right)\right)}^2+\frac{4b^2+4\mu M^2}{\lambda }\right)R\left(s,\omega \right)\\ +\left(\frac{2}{\epsilon }{\Vert A^{\frac{1}{2}}q\Vert }_2^2+\frac{4}{\lambda }{\Vert A^{\frac{5}{2}}q\Vert }^2+\frac{4{\epsilon }^2}{\lambda }{\Vert A^{\frac{1}{2}}q\Vert }^2+\frac{1}{2}{\Vert Aq\Vert }^2\right){\left|W\left(t\right)\right|}^2,\tau \le t\le 0.\end{array}$

Using the Gronwall lemma, we get that

$\begin{array}{c}E\left(0\right)\le {\text{e}}^{\epsilon \tau }{\Vert A^{\frac{1}{2}}q\Vert }^2+C{\displaystyle {\int }_{\tau }^0}{\text{e}}^{\epsilon s}\left(\epsilon \left|p\right|+{\left(\left|p\right|+R\left(s\mathrm{,}\omega \right)\right)}^2+\frac{4b^2+4\mu M^2}{\lambda }\right)R\left(s\mathrm{,}\omega \right)\text{d}s\\ +\left(\frac{2}{\epsilon }{\Vert A^{\frac{1}{2}}q\Vert }_2^2+\frac{4}{\lambda }{\Vert A^{\frac{5}{2}}q\Vert }^2+\frac{4{\epsilon }^2}{\lambda }{\Vert A^{\frac{1}{2}}q\Vert }^2+\frac{1}{2}{\Vert Aq\Vert }^2\right){\displaystyle {\int }_{\tau }^0}{\text{e}}^{\epsilon s}{\left|W\left(s\right)\right|}^2\text{d}s.\end{array}$ (4.36)

Set

$\begin{array}{c}{r}_2^2\left(\omega \right)=\left(\frac{2}{\epsilon }{\Vert A^{\frac{1}{2}}q\Vert }_2^2+\frac{4}{\lambda }{\Vert A^{\frac{5}{2}}q\Vert }^2+\frac{4{\epsilon }^2}{\lambda }{\Vert A^{\frac{1}{2}}q\Vert }^2+\frac{1}{2}{\Vert Aq\Vert }^2\right){\displaystyle {\int }_{-\infty }^0}{\text{e}}^{\epsilon s}{\left|W\left(s\right)\right|}^2\text{d}s\\ +{\Vert A^{\frac{1}{2}}q\Vert }^2\underset{\tau \le 0}{sup}{\text{e}}^{\epsilon \tau }{\left|W\left(\tau \right)\right|}^2\\ +C{\displaystyle {\int }_{-\infty }^0}{\text{e}}^{\epsilon s}\left(\epsilon \left|p\right|+{\left(\left|p\right|+R\left(s\mathrm{,}\omega \right)\right)}^2+\frac{4b^2+4\mu M^2}{\lambda }\right)R\left(s\mathrm{,}\omega \right)\text{d}s\mathrm{.}\end{array}$

Since $li{m}_{t\to \infty }\frac{W\left(t\right)}{t}=0$ , $r_2^2\left(\omega \right)$ is finite P-a.s., together with 4.18 and 4.36, we get that

${\Vert A^{\frac{1}{2}}{Y}_2\left(\mathrm{0,}\omega \mathrm{<mo>\; ;\; </mo>}{Y}_2\left(\tau \mathrm{,}\omega \right)\right)\Vert }_E^2\le \frac{r_2^2\left(\omega \right)}{C\left(p\right)}\mathrm{.}$

This completes the proof of Lemma 4.4. ,

Theorem 4.5. Let $p<\frac{\sqrt{\lambda }}{4}$ , (1.2)-(1.4) hold, $q\left(x\right)\in H^3\left(\Omega \right)\cap H_0^2\left(\Omega \right)$ , then the

random dynamical system $S_{\epsilon }\left(t\mathrm{,}\omega \right)$ possesses a nonempty compact random attractor $\mathcal{A}$ .

Proof. Let $B_1\left(\omega \right)$ be the ball of $H^3\left(\Omega \right)\times \left(H^2\left(\Omega \right)\cap H_0^1\left(\Omega \right)\right)$ of radius

$\frac{r_2\left(\omega \right)}{\sqrt{C\left(p\right)}}$ , by the compact embedding $H^3\left(\Omega \right)$ , it follows that $B_1\left(\omega \right)$ is

compact in E. for every bounded non-random set B of E and any $\psi \left(0\right)\in {\tilde{S}}_{\epsilon }\left(t\mathrm{,}{\theta }_{-t}\omega \right)B$ , by Lemma 4.4, we know that $Y_2\left(0\right)\mathrm{<mo>\; =\; </mo>}\psi \left(0\right)-Y_1\left(0\right)\in B_1\left(\omega \right)$ . Therefore, for $\tau \le 0$ ,

$\begin{array}{l}\underset{l\left(0\right)\in B_1\left(\omega \right)}{\inf }{\Vert \psi \left(0\right)-l\left(0\right)\Vert }_E^2\le {\Vert Y_1\left(0\right)\Vert }_E^2\\ \le \frac{{\text{e}}^{\epsilon \tau }}{\left(1-\epsilon \right)C\left(p\right)}\left({\Vert u_0\Vert }_2^2-\left(p-{\left|\nabla u_0\right|}^2\right){\left|\nabla u_0\right|}^2+{\Vert u_1+\epsilon u_0\Vert }^2\right).\end{array}$

So, for all $t\ge 0$ ,

$\begin{array}{l}d\left({\tilde{S}}_{\epsilon }\left(t\mathrm{,}{\theta }_{-t}\omega \right)B\mathrm{,}{B}_1\left(\omega \right)\right)\\ \le \frac{{\text{e}}^{-\epsilon \tau }}{\left(1-\epsilon \right)C\left(p\right)}\left({\Vert u_0\Vert }_2^2-\left(p-{\left|\nabla u_0\right|}^2\right){\left|\nabla u_0\right|}^2+{\Vert u_1+\epsilon u_0\Vert }^2\right)\mathrm{.}\end{array}$

From relation (4.3) between $S_{\epsilon }\left(t\mathrm{,}\omega \right)$ and ${\tilde{S}}_{\epsilon }\left(t\mathrm{,}\omega \right)$ , w can obtain that for any non-random bounded $B\subset EP-a\mathrm{.}s\mathrm{.}$ ,

$d\left({\tilde{S}}_{\epsilon }\left(t\mathrm{,}{\theta }_{-t}\omega \right)B\mathrm{,}{B}_1\left(\omega \right)\right)\to \mathrm{0,}\text{as}t\to +\infty \mathrm{.}$

Hence, the RDS $S_{\epsilon }\left(t\mathrm{,}\omega \right)$ associated with (3.7) possesses a uniformly attracting compact set $B_1\left(\omega \right)\subset E$ . Using Theorem 2.1, we complete the proof.

Acknowledgements

This work is supported by the National Natural Sciences Foundation of People’s Republic of China under Grant 11561076.

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	Lazer, A.C. and McKenna, P.J. (1990) Large-Amplitude Periodic Oscillations in Suspension Bridges: Some New Connection with Nonlinear Analysis. SIAM Review, 32, 537-578. https://doi.org/10.1137/1032120
[2]	Qin, C.L. and Lin, G.G. (2010) Random Attractor for Weakly Damped Stochastic Kirchhoff Equation. Journal of Yunnan University (Natural Sciences Edition), 32, 101-108.
[3]	Qin, C.L. and Du, J.J. (2016) Random Attractor for Strongly Damped Stochastic Kirchhoff Equation. Journal of Science of Teachers College and University, 36, 6-11.
[4]	Ma, Q.Z. (2005) Asymptotic Behavior of the Solution for the Floating Beam. Journal of Northwest Normal University (Natural Science), 41, 4-7. (in Chinese)
[5]	Zhong, C.K., Ma, Q.Z. and Sun, C.Y. (2007) Existence of Strong Solutions and Global Attractors for the Suspension Bridge Equations. Nonlinear Analysis, 67, 442-454.
[6]	Park, J.Y. and Kang, J.R. (2011) Uniform Attractor for Non-Autonomous Suspension Bridge Equations with Localized Damping. Mathematical Methods in the Applied Sciences, 34, 487-496.
[7]	Ma, Q.Z., Wang, S.P. and Chen, X.B. (2011) Uniform Attractors for the Coupled Suspension Bridge Equations. Applied Mathematics and Computation, 217, 6604-6615.
[8]	Crauel, H. and Flandoli, F. (1994) Attractors for Random Dynamical Systems. Probability Theory and Related Fields, 100, 365-393. https://doi.org/10.1007/BF01193705
[9]	Crauel, H., Debussche, A. and Flandoli, F. (1997) Random Attractors. Journal of Dynamics and Differential Equations, 9, 307-314. https://doi.org/10.1007/BF02219225
[10]	Arnold, L. (1998) Random Dynamical Systems. Springer, New York. https://doi.org/10.1007/978-3-662-12878-7
[11]	Zhou, S.F., Yin, F.Q. and Ouyang, Z.G. (2008) Random Attractor for Damped Nonlinear Wave Equation with White Noise. SIAM Journal on Applied Dynamical Systems, 4, 883-903. ttps://doi.org/10.1137/050623097
[12]	Fan, X.M. (2008) Random Attractor for a Damped Stochastic Wave Equation with Multiplicative Noise. International Journal of Mathematics, 19, 421-437. https://doi.org/10.1142/S0129167X08004741
[13]	Pazy, A. (1983) Semigroup of Linear Operators and Applications to Partial Differential Equations. Springer, New York. https://doi.org/10.1007/978-1-4612-5561-1