Global stability of the Euler-Bernoulli beams excited by
multiplicative white noises
Zhenzhen Li
Division of Dynamics and Control, School of Mathematics and Statistics, Shandong University of Technology, ZiBo
255000, China
E-mail: zzlicncn@163.com
Kun Zhao
Beijing Electro-Mechanical Engineering Institute, Beijing 100074, P. R. China
E-mail: zhaokunhit@yeah.net
Hongkui Li
Division of Dynamics and Control, School of Mathematics and Statistics, Shandong University of Technology, ZiBo
255000, China
E-mail: Lhk8068@163.com
Juan L.G. Guirao
Department of Applied Mathematics and Statistics, Technical University of Cartagena, Hospital de Marina, Cartagena
30203, Spain
E-mail: juan.garcia@upct.es
Huatao Chen*
Division of Dynamics and Control, School of Mathematics and Statistics, Shandong University of Technology, ZiBo
255000, China
E-mail: htchencn@sdut.edu.cn
Reception: 05/01/2023 Acceptance: 10/03/2023 Publication: 30/03/2023
Suggested citation:
Zhenzhen Li, Kun Zhao, Hongkui Li, Juan L.G. Guirao and Huatao Chen. (2023). Global stability of the Euler-
Bernoulli beams excited by multiplicative white noises. 3C Tecnología. Glosas de innovación aplicada a la pyme,
12 (01), 386-412. https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
386
Global stability of the Euler-Bernoulli beams excited by
multiplicative white noises
Zhenzhen Li
Division of Dynamics and Control, School of Mathematics and Statistics, Shandong University of Technology, ZiBo
255000, China
E-mail: zzlicncn@163.com
Kun Zhao
Beijing Electro-Mechanical Engineering Institute, Beijing 100074, P. R. China
E-mail: zhaokunhit@yeah.net
Hongkui Li
Division of Dynamics and Control, School of Mathematics and Statistics, Shandong University of Technology, ZiBo
255000, China
E-mail: Lhk8068@163.com
Juan L.G. Guirao
Department of Applied Mathematics and Statistics, Technical University of Cartagena, Hospital de Marina, Cartagena
30203, Spain
E-mail: juan.garcia@upct.es
Huatao Chen*
Division of Dynamics and Control, School of Mathematics and Statistics, Shandong University of Technology, ZiBo
255000, China
E-mail: htchencn@sdut.edu.cn
Reception: 05/01/2023 Acceptance: 10/03/2023 Publication: 30/03/2023
Suggested citation:
Zhenzhen Li, Kun Zhao, Hongkui Li, Juan L.G. Guirao and Huatao Chen. (2023). Global stability of the Euler-
Bernoulli beams excited by multiplicative white noises. 3C Tecnología. Glosas de innovación aplicada a la pyme,
12 (01), 386-412. https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
ABSTRACT
This paper considers the global stability of the Euler-Bernoulli beam excited by multiplicative white
noise. Based on the theory of global random attractors, the Hausdorff dimensions of the global random
attractors for the system is obtained. According to the relationship between Hausdorff dimensions and
global Lyapunov exponents, the global stability of the stochastic beam is derived.
KEYWORDS
Global stochastic stability, Stochastic Euler-Bernoulli beam, Global Lyapunov exponents
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
387
1 INTRODUCTION
The dynamics of the beams are important problem in the elastokinetics [1]. Tajik [2] proposed the
stability analysis of motion equations of unbalanced spinning pre-twisted beam. Dai [3] investigated the
limit point bifurcations and jump of cantilevered microbeams according to Galerkin method and modal
truncation. The bifurcation and chaos for transverse motion of axially accelerating viscoelastic beams
was studied by Chen [4]. Pellicano [5] analyzed the linear subcritical behavior, bifurcation analysis
and stability of a simply supported beam subjected to an axial transport of mass. Based on the direct
method of multiple scales, Mao [6] dealt with problems of stability and saddle-node bifurcations for
supercritically moving beam. Using Melnikov method [7], Zhang [8, 9] investigated the multi-pulse
global bifurcations of a cantilever beam. Zhou [10] studied the chaos and subharmonic bifurcation of
a composite laminated buckled beam with a lumped mass. Applying the phase plane and positive
position feedback approach, Hamed [11] investigated the stability and bifurcation of the cantilever
beam system which carrying an intermediate lumped mass, to name but a few. For more details, one
can see refer to Refs [12–15] and the references therein.
Stochastic stability is one of the most important issue in the research areas of stochastic dynamics [16].
It is well known that, when the problem associated with stability are considered, there exists a very
useful tool named Lyapunov exponents which can be distinguished as local type and global type [17].
Generally speaking, the solution of system is stability when biggest Lyapunov exponent associated is less
than 0. With respect to the investigations on the stability and other dynamics behaviors for the beams
by using local Lyapunov exponents, one can refer to the Refs [18
20] and the references therein. Similarly,
the global Lyapunov exponents are also the powerful tools in studying the global dynamics for stochastic
beams. Unfortunately, calculating the global Lyapunov exponents is not an easy thing, but if we only
consider the global stability of the system, we can use the Hausdorff dimension of global attractors
associated with system to describe the signs of the biggest global Lyapunov exponent. The method
which can be used to get the Hausdorff dimension estimations associated with the global Lyapunov
exponents was due to Debussche [21]. Employing this method, together with the support relationship
between global random attractors and probability invariant measures proposed by Crauel [22,23]. Chen
et al [24] consider the global dynamics of the Euler-Bernoulli beams with additive white noises. With
respect to investigation on the global dynamics of the nonautonomous the Euler-Bernoulli beams by
global attractors theory, see Chen et al [25].
Let
D
= (0
,L
), this paper consider the Euler-Bernoulli beam equation excited multiplicative white
noise in the following form
utt +α(ut)2u+β∥∇u2p(∆)u=σu ˙
W, (1)
with the hinged boundary condition
x=0:u=uxx = 0; x=L:u=uxx =0,(2)
and the initial value
t=τ:u=u0,u
t=u1,(3)
where
u
=
t
(
t, x
)
,xD
is the lateral displacement of the beam,
α
(
ut
)denotes the damping,
β>
0
are constants, the negative and positive of
pR
can show the stretch and compress of the beam.
∥∇u2
denotes the geometry of the beam bending for its elongation.
W
is the one dimensional two-sided
real-valued standard Wiener process, σu ˙
Wrepresents the multiplicative white noise.
Let
u∥≡∥uL2(D)
,
us≡∥uHs
0(D)
,(
u, v
)
(
u, v
)
L2(D)
,(
u, v
)
s
=(
u, v
)
Hs
0(D)
, where
Hs
(
D
),
Hs
0
(
D
),
sR
are the usual Sobolev Spaces,for more detailed, see [26].
A
=∆
2
with boundary
condition
(2)
, then
D
(
A
)=
{u|uH4
(
D
)
H1
0
(
D
)
,
u
=0
}
, and then
A
is self-adjoint, positive,
unbounded linear operators and
A1L
(
L2
(
D
)) is compact. then, their eigenvalues
{λi}iN
satisfy
0
1λ2≤···→∞
and the corresponding eigenvalues
{ei}
i=1
form an orthonormal basis in
L2
(
D
).
Following the the mechanism in [27] p55, the power of (
∆)
s,sR
can also be defined,particularly,
D
(
A1
2
)=
H1
0
(
D
)
H2
(
D
). Moreover, for any
s1,s
2R,s
1>s
2
,
D
(
As1
)can be compact imbedding in
D(As2),and the following holds
us1λ
s1s2
2
1us2,u∈D(As1).(4)
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
388
1 INTRODUCTION
The dynamics of the beams are important problem in the elastokinetics [1]. Tajik [2] proposed the
stability analysis of motion equations of unbalanced spinning pre-twisted beam. Dai [3] investigated the
limit point bifurcations and jump of cantilevered microbeams according to Galerkin method and modal
truncation. The bifurcation and chaos for transverse motion of axially accelerating viscoelastic beams
was studied by Chen [4]. Pellicano [5] analyzed the linear subcritical behavior, bifurcation analysis
and stability of a simply supported beam subjected to an axial transport of mass. Based on the direct
method of multiple scales, Mao [6] dealt with problems of stability and saddle-node bifurcations for
supercritically moving beam. Using Melnikov method [7], Zhang [8, 9] investigated the multi-pulse
global bifurcations of a cantilever beam. Zhou [10] studied the chaos and subharmonic bifurcation of
a composite laminated buckled beam with a lumped mass. Applying the phase plane and positive
position feedback approach, Hamed [11] investigated the stability and bifurcation of the cantilever
beam system which carrying an intermediate lumped mass, to name but a few. For more details, one
can see refer to Refs [12–15] and the references therein.
Stochastic stability is one of the most important issue in the research areas of stochastic dynamics [16].
It is well known that, when the problem associated with stability are considered, there exists a very
useful tool named Lyapunov exponents which can be distinguished as local type and global type [17].
Generally speaking, the solution of system is stability when biggest Lyapunov exponent associated is less
than 0. With respect to the investigations on the stability and other dynamics behaviors for the beams
by using local Lyapunov exponents, one can refer to the Refs [18
20] and the references therein. Similarly,
the global Lyapunov exponents are also the powerful tools in studying the global dynamics for stochastic
beams. Unfortunately, calculating the global Lyapunov exponents is not an easy thing, but if we only
consider the global stability of the system, we can use the Hausdorff dimension of global attractors
associated with system to describe the signs of the biggest global Lyapunov exponent. The method
which can be used to get the Hausdorff dimension estimations associated with the global Lyapunov
exponents was due to Debussche [21]. Employing this method, together with the support relationship
between global random attractors and probability invariant measures proposed by Crauel [22,23]. Chen
et al [24] consider the global dynamics of the Euler-Bernoulli beams with additive white noises. With
respect to investigation on the global dynamics of the nonautonomous the Euler-Bernoulli beams by
global attractors theory, see Chen et al [25].
Let
D
= (0
,L
), this paper consider the Euler-Bernoulli beam equation excited multiplicative white
noise in the following form
utt +α(ut)2u+β∥∇u2p(∆)u=σu ˙
W, (1)
with the hinged boundary condition
x=0:u=uxx = 0; x=L:u=uxx =0,(2)
and the initial value
t=τ:u=u0,u
t=u1,(3)
where
u
=
t
(
t, x
)
,xD
is the lateral displacement of the beam,
α
(
ut
)denotes the damping,
β>
0
are constants, the negative and positive of
pR
can show the stretch and compress of the beam.
∥∇u2
denotes the geometry of the beam bending for its elongation.
W
is the one dimensional two-sided
real-valued standard Wiener process, σu ˙
Wrepresents the multiplicative white noise.
Let
u∥≡uL2(D)
,
us≡∥uHs
0(D)
,(
u, v
)
(
u, v
)
L2(D)
,(
u, v
)
s
=(
u, v
)
Hs
0(D)
, where
Hs
(
D
),
Hs
0
(
D
),
sR
are the usual Sobolev Spaces,for more detailed, see [26].
A
=∆
2
with boundary
condition
(2)
, then
D
(
A
)=
{u|uH4
(
D
)
H1
0
(
D
)
,
u
=0
}
, and then
A
is self-adjoint, positive,
unbounded linear operators and
A1L
(
L2
(
D
)) is compact. then, their eigenvalues
{λi}iN
satisfy
0
1λ2≤···→
and the corresponding eigenvalues
{ei}
i=1
form an orthonormal basis in
L2
(
D
).
Following the the mechanism in [27] p55, the power of (
∆)
s,sR
can also be defined,particularly,
D
(
A1
2
)=
H1
0
(
D
)
H2
(
D
). Moreover, for any
s1,s
2R,s
1>s
2
,
D
(
As1
)can be compact imbedding in
D(As2),and the following holds
us1λ
s1s2
2
1us2,u∈D(As1).(4)
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
Suppose
α
(
ut
)=
αA1
2ut
, here
α>
0is a constant,
A1
2ut
is called the damping with strongly form, then
following abstract form of Euler-Bernoulli beam equation excited multiplicative white noises
utt +αA1
2ut+Au +β∥∇u2p(∆)u=σu ˙
W,
x=0:u=uxx = 0; x=L:u=uxx =0,
u(x, τ)=u0(x),∂
tu(x, τ)=u1(x).
(5)
The rest of paper is organized as follows. Preliminaries and main lemmas are listed in Section 2. Section
3 is devoted to derive main results and proof associated. The proofs for the main lemmas are given in
Section 4. Finally, the conclusions are presented in section (4).
2 PRELIMINARIES AND MAIN LEMMAS
2.1 PRELIMINARIES
Let
E1
=
D
(
A1
2
)
×L2
(
D
)
,E
2
=
D
(
A3
4
)
×H1
0
(
D
)equipped with Graph norms and the induced inner
products, then they are all Hilbert spaces. Let
(X, ∥·∥
X)
be a complete separable metric space with Borel
σ
-algebra
B
(
X
)and (Ω
,F,P
)be a probability space. We consider Ω=
{ω|ω
(
·
)
C
(
R,R
)
(0) = 0
},F
is the
σ
-algebra and
P
is the Wiener measure. Set a family of measure preserving and ergodic
transformations θtω(·)=ω(·+t)ω(·),tR, Consider the following system
dz +µzdt =dW
z(−∞)=0,(6)
and the solution of system (6) is given by
z(θtω) := µ0
−∞
eµτ (θtω)(τ). (7)
z
(
θtω
)is Ornstein-Uhlenbeck process (in Short O-U process). The following results on O-U process
belong to Fan [28].
Lemma 1. The Ornstein-Uhlenbeck process z(θtω)defined in system (7) satisfies
E[|z(θtω)|]= 1
πµ,E|z(θtω)|2=1
2µ,(8)
and there exists a constant t1(ω)>0satisfying
0
t|z(θsω)|ds < 1
πµt, 0
t|z(θeω)|2ds < 1
2µt, tt1,(9)
and the mapping t→ z(θtω)grows sublinearly, i.e.
lim
t→±∞
z(θtω)
t=0.
Moreover If µ2β, β > 0, then
Eeβs+t
s|z(θτω)|2 e
βt
µ,sR,t0,(10)
when µ3r2,r 0, the following holds
Eers+t
s|z(θτω)| e
rt
µ,sR,t0.(11)
The following is Random dynamical system which is due to Aronld [29].
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
389
Definition 1. The flow
{θt:Ω,tR}
is a family of measure preserving transformations in
probability space, such that (
t, ω
)
θtω
is measurable,
θ0
=
id, θt+s
=
θtθs
for all
s, t R
. Then
,F, P, {θt}tRis called a metric dynamical system.
Definition 2. A random dynamical system (RDS) on Polish space (
X, d
)with Borel
σ
algebra
B
(
X
)
on ,F, P, {θt}tRis a measurable mapping
ϕ:R+××XX
(t,x)→ ϕ(t, ω)x
such that Pa.s.
1. ϕ(0)=id on X.
2. ϕ(t+s, ω)=ϕ(t, θsω)ϕ(s, ω)for all s, t R+.
The theory of global random attractors is as following, one can refer to Crauel and Flandoli [22, 30] and
Schmalfuss [31].
Definition 3. A random set
K
(
ω
)is said to absorb the set
BX
for a RDS
ϕ
, if
Pa.s
. there
exists tB(ω)such that
ϕ(t, θtω)BK(ω),ttB(ω).
Definition 4. Let
B⊂
2
X
is a collection of subsets of
X
, then a closed random set
A
(
ω
)is called
random attractor associated with the RDS ϕ, if Pa.s.
1. A(ω)is a random compact set.
2. A(ω)is invariant i.e. ϕ(t, ω),A(ω)=A(θtω)for all t0.
3. For every B∈B,
lim
t→∞ dist (ϕ(t, θtω)B,A(ω)) = 0,
where dist (·,·)denotes the Hausdorff semidistance defined by
dist(A, B) = sup
xA
inf
yBd(x, y), A, B X.
The following random fixed point is important in investigating the global stability.
Definition 5. [29, 32] Let
φ
(
t, ω
)be a RDS,
a
(
ω
)is a random set and consists of one point (
Pa.s
).
a(ω)is called the random fixed point if the following holds
φ(t, ω)a(ω)=a(θtω),tR+.
The coming theorem is very useful to verify the existence of global random attractors in this paper.
Theorem 1. [24] Suppose
Sε
(
t, ω
)is a RDS on Polish space (
X, d
), and suppose that
ϕ
possesses an
absorbing set in
X
and for any nonrandom bounded set
BX, lim
t+Sε(t, θtω)B
is relative compact
P-a.s. Then ϕpossesses uniqueness random attractors defined by the following
A(ω)=
BX
ΛB(ω),
where union is taken over all bounded BX, and ΛB(ω)given by
ΛB(ω)=
s0
ts
ϕ(t, θtω)B.
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
390
Definition 1. The flow
{θt:Ω,tR}
is a family of measure preserving transformations in
probability space, such that (
t, ω
)
θtω
is measurable,
θ0
=
id, θt+s
=
θtθs
for all
s, t R
. Then
,F, P, {θt}tRis called a metric dynamical system.
Definition 2. A random dynamical system (RDS) on Polish space (
X, d
)with Borel
σ
algebra
B
(
X
)
on ,F, P, {θt}tRis a measurable mapping
ϕ:R+××XX
(t,x)→ ϕ(t, ω)x
such that Pa.s.
1. ϕ(0)=id on X.
2. ϕ(t+s, ω)=ϕ(t, θsω)ϕ(s, ω)for all s, t R+.
The theory of global random attractors is as following, one can refer to Crauel and Flandoli [22, 30] and
Schmalfuss [31].
Definition 3. A random set
K
(
ω
)is said to absorb the set
BX
for a RDS
ϕ
, if
Pa.s
. there
exists tB(ω)such that
ϕ(t, θtω)BK(ω),ttB(ω).
Definition 4. Let
B⊂
2
X
is a collection of subsets of
X
, then a closed random set
A
(
ω
)is called
random attractor associated with the RDS ϕ, if Pa.s.
1. A(ω)is a random compact set.
2. A(ω)is invariant i.e. ϕ(t, ω),A(ω)=A(θtω)for all t0.
3. For every B∈B,
lim
t dist (ϕ(t, θtω)B,A(ω)) = 0,
where dist (·,·)denotes the Hausdorff semidistance defined by
dist(A, B) = sup
xA
inf
yBd(x, y), A, B X.
The following random fixed point is important in investigating the global stability.
Definition 5. [29, 32] Let
φ
(
t, ω
)be a RDS,
a
(
ω
)is a random set and consists of one point (
Pa.s
).
a(ω)is called the random fixed point if the following holds
φ(t, ω)a(ω)=a(θtω),tR+.
The coming theorem is very useful to verify the existence of global random attractors in this paper.
Theorem 1. [24] Suppose
Sε
(
t, ω
)is a RDS on Polish space (
X, d
), and suppose that
ϕ
possesses an
absorbing set in
X
and for any nonrandom bounded set
BX, lim
t+Sε(t, θtω)B
is relative compact
P-a.s. Then ϕpossesses uniqueness random attractors defined by the following
A(ω)=
BX
ΛB(ω),
where union is taken over all bounded BX, and ΛB(ω)given by
ΛB(ω)=
s0
ts
ϕ(t, θtω)B.
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
Some powerful transformation on system (5) are derived as following. Set v1=ut, we get
ut=v1,
v1t=αA1
2v1Au β∥∇u2p(∆)u+σu ˙
W,
x=0:u=ux= 0; x=L:u=uxx =0,
u(x, τ)=u0(x),v
1(x, τ)=u1(x),
(12)
let v2=v1+εu, here ε>0, we have
dU
dt =QU +X1(ω, U),Uτ=(u0,u
1+εu0)T,(13)
where
U=U1
U2=u
v2,Q=εI, I
A+εαA1
2ε,αA1
2εI,
X1(ω, U)=X11(ω, U)
X12(ω, U)=0
β∥∇u2p(∆)u+σu ˙
W.
System
(5)
and system
(13)
are equivalent, thus, the dynamical behavior of system
(5)
can be reflected
by system (13).
Set v=ut+εu σuz (θtω), where z(θtω)O-U process formulated by (6), it gives that
dV
dt =QV +X2(θtω)+X3(V),Vτ=(u0,u
1+εu0σu0z(θtω))T,(14)
where
V=V1
V2=u
v,X3(V)=X31(V)
X32(V)=0
β∥∇u2p(∆)u,
X2(θtω)=X21 (θtω)
X22 (θtω)=σuz (θtω)
σµαA1
2+2εσz (θtω)uz (θtω)σvz (θtω).
System (14) is a system with random coefficient, which can be studied ωby ω.
2.2 MAIN LEMMAS
In order to obtain the global stochastic stability of the system
(13)
based on the global random
attractors theory, the first step should be utilized to verify the system (13) can induce a RDS.
Lemma 2. For any
τR
and initial value
VτE1
, system
(14)
possesses a unique local mild solution
V(t;Vτ)C([τ,τ +T],E
1),t[τ,τ +T],T>0.
Let
φ
(
t
)be the solution mapping determined by system
(14)
, which means
V(t;Vτ)
=
φ(t)Vτ, we have φ(t, 0)=φ(0,t, θtω),t0,ω. Defining
S(t, ω) := φ(t, 0),t0,ω,
gives the RDS associated with system
(14)
, which together with the relationship between
ιε
:
(u, v1)T
[
y, v1σuz (θtω)
]
T
=
(u, v)T
implies that system
(13)
can also generate a RDS
Sε
(
t, ω
)with the
following from
Sε(t, ω)=ι1(θtω),S(t, ω)ι(ω):E1E1,
is the RDS induced by system
(13)
. The following are very important in proof for the existence of the
global random attractors for the system (14).
Let
p2
3λ
1
4
1,0= min
1,
α8λ
1
4
112p
λ
1
4
110λ
1
4
115p+4α2λ
1
4
1
,(15)
then
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
391
Lemma 3. For any given U=[U1,U
2]TE1, the following holds
(QU ,U)E1≤−ε
2U2
E1ε
4U223
4A1
4U12.
Based on the Lemma 3, we have
Lemma 4. For any given bounded set
BE1
, there exists a random variable
r1
(
ω
)
>
0and
TB
(
ω
)
0,
for tTB(ω), the following holds,
Sε(t, θtω)BE1r1(ω).
Furthermore,
Er2
1(ω)<.
Lemma 4 shows that
Sε
(
t, ω
)has a global absorbing set. In order to obtain the existence of the global
random attractors for the system
Sε
(
t, ω
), besides the existence of the global absorbing set, we need
verify that the Sε(t, ω)is asymptotically compact.
To begin with, we decompose the solution
V
generated by system
(13)
with the initial value
Uτ
(
ω
)=
(u0,u
1+εu0)T
into two parts
U
=
Ua
+
Ub
=
(ua,u
a
t+εua)T
+
ub,u
b
t+εubT, where Ua
solves
dUa
dt =QU a+0, σua˙
WT
Ua
τ(ω)=0,
(16)
and Ubsolves
dUb
dt =QU b+0,β∥∇u2p(∆)u+σub˙
WT
Ub
τ=(u0,u
1+εu0)T.
(17)
Split the solution
V
of system
(14)
with the initial value
Vτ
(
ω
)=
(u0,u
1+εu0σu0z(θτω))T
into
two parts
V
=
Va
+
Vb
=
(ua,v
a)T
+
ub,v
bT
=
(ua,u
a
t+εua)T
+
ub,u
b
t+εubσuz (θtω)T,where Va
solves
dVa
dt =QV a+σuaz(θtω)µαA 1
2+2εσz (θtω)uaz(θtω)σvaz(θtω)T
Va
τ(ω)=0,
(18)
and Vbsolves
dVb
dt =QV b+σubz(θtω)µαA 1
2+2εσz (θtω)ubz(θtω)σvbz(θtω)T
+0,β∥∇u2puT
Vb
τ=(u0,u
1+εu0σu0z(θτω))T.
(19)
For the solutions of system (18) and (19), we have the following priori estimates respectively.
Lemma 5. The solution Uaof system (16) satisfies
lim
t+
A1
4Ua
E1
=0.
Lemma 6. There exists r2(ω)>,T
2(ω)0, the solutions Ubsatisfies
Ub
E1r2(ω),tT2(ω).
Lemma 7. There exists r3(ω)>0,T
3(ω)0, the following holds
A1
4Ub
E1r3(ω),tT3(ω).
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
392
Lemma 3. For any given U=[U1,U
2]TE1, the following holds
(QU ,U)E1≤−ε
2U2
E1ε
4U223
4A1
4U12.
Based on the Lemma 3, we have
Lemma 4. For any given bounded set
BE1
, there exists a random variable
r1
(
ω
)
>
0and
TB
(
ω
)
0,
for tTB(ω), the following holds,
Sε(t, θtω)BE1r1(ω).
Furthermore,
Er2
1(ω)<.
Lemma 4 shows that
Sε
(
t, ω
)has a global absorbing set. In order to obtain the existence of the global
random attractors for the system
Sε
(
t, ω
), besides the existence of the global absorbing set, we need
verify that the Sε(t, ω)is asymptotically compact.
To begin with, we decompose the solution
V
generated by system
(13)
with the initial value
Uτ
(
ω
)=
(u0,u
1+εu0)T
into two parts
U
=
Ua
+
Ub
=
(ua,u
a
t+εua)T
+
ub,u
b
t+εubT, where Ua
solves
dUa
dt =QU a+0, σua˙
WT
Ua
τ(ω)=0,
(16)
and Ubsolves
dUb
dt =QU b+0,β∥∇u2p(∆)u+σub˙
WT
Ub
τ=(u0,u
1+εu0)T.
(17)
Split the solution
V
of system
(14)
with the initial value
Vτ
(
ω
)=
(u0,u
1+εu0σu0z(θτω))T
into
two parts
V
=
Va
+
Vb
=
(ua,v
a)T
+
ub,v
bT
=
(ua,u
a
t+εua)T
+
ub,u
b
t+εubσuz (θtω)T,where Va
solves
dVa
dt =QV a+σuaz(θtω)µαA 1
2+2εσz (θtω)uaz(θtω)σvaz(θtω)T
Va
τ(ω)=0,
(18)
and Vbsolves
dVb
dt =QV b+σubz(θtω)µαA 1
2+2εσz (θtω)ubz(θtω)σvbz(θtω)T
+0,β∥∇u2puT
Vb
τ=(u0,u
1+εu0σu0z(θτω))T.
(19)
For the solutions of system (18) and (19), we have the following priori estimates respectively.
Lemma 5. The solution Uaof system (16) satisfies
lim
t+
A1
4Ua
E1
=0.
Lemma 6. There exists r2(ω)>,T
2(ω)0, the solutions Ubsatisfies
Ub
E1r2(ω),tT2(ω).
Lemma 7. There exists r3(ω)>0,T
3(ω)0, the following holds
A1
4Ub
E1r3(ω),tT3(ω).
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
The variation equations of system (13)
d
U
dt =Q
U+
X(U)
U,(20)
where
X(U)=
X1(U),
X3(U)
X2(U),
X4(U)
=0,0
β∥∇u2p(∆) ·+2β(u, ∇·)(∆)u+σ˙
W,0.
On the other hand
d
V
dt =Q
V+
X2(V)
V+
X3(V)
V,(21)
here
V={
V1,
V2},
X2((θtω)) =
X21((θtω)),
X23((θtω))
X22((θtω)),
X24((θtω))
=σz (θtω),0
σ(µα+2εσz (θtω)) z(θtω),σz (θtω),
X3(V)=
X31(V),
X33(V)
X32(V),
X34(V)
=0,0
β∥∇u2p(∆) ·+2β(u, ∇·)(∆)u, 0.
Let
U
=
U1,
U2
be the solution of system
(20)
with initial
t
=0:
U0E1
and
V
=
V1,
V2
be the
solution of system (21) with initial value t=0:
V0=IE1.
Let U(1) =U(1)
1,U(1)
2,U(2) =U(2)
1,U(2)
2are two solutions of system (13), then
dU(1) U(2)
dt =QU(1) U(2)+X1U(1)X1U(2),(22)
Analogously, let
V(1)
=
V(1)
1,V(1)
2
=
[u1,v
1],V(2)
=
V(2)
1,V(2)
2
=
[u2,v
2]
are two solutions of
system (14) with initial values V(1)
0,V(2)
0, where V(2)
0=V(1)
0+I,I=[I1,I
2]E1, then
dV(1) V(2)
dt =QV(1) V(2)+X2V(1)X2V(2)
+X3V(1)X3V(2).
(23)
In addition, set Γ= 1,Γ2]T=U(1) U(2)
U, we have the following two Lemmas.
Lemma 8. For u1,u
2∈D(A1
2), there exist constants c1(ω),c
2(ω)in (51), such that
X32 (u1)
X32 (u2)
LD(A1
2),L2(D)c1(ω)
A1
2u1A1
2u2
,
X32 (u1)
LD(A1
2),L2(D)c2(ω).
(24)
Furthemnore,
E(c1(ω)) <,E(c2(ω)) <.
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
393
Lemma 9. Let Sε(ω) := Sε(t, ω), then Sε(ω)is uniformly quasidifferentiable on A(ω).
The following relationship between Hausdorff dimension and global Lyapunov exponents which will
be used to assert the global stochastic stability.
Theorem 2. [21]Suppose that, for
k
=1
,...,d, sup
uA(ω)
ωk
(
DSε
(
ω, u
)) is integrable. Then we can
generalize the notion of global Lyapunov exponents introduced in Ref. [17] by setting
Λk=Eln sup
uA(ω)
ωk(DSε(ω, u))Eln sup
uA(ω)
ωk1(DSε(ω, u)),
for k2and
Λ1=Eln sup
uA(ω)
ω1(DSε(ω, u)).
And it is easy to see that there exists ωdsatisfying (26)(27) if and only if
Λ1+...
d<0.
3 MAIN RESULTS AND PROOFS
By Lemma 4, we get that
Sε
(
t, ω
)has a global absorbing set, which along with Lemma 5, Lemma 6 and
Lemma 7 gives that
Sε
(
t, ω
)is asymptotically compact. Thus, employing the Theorem 1, it is asserted
that there exist global random attractors A(ω) for system 13.
Based on the results and the estimations on global random attractors, the Hausdorff dimension of
the global random attractors
A
(
ω
)can be got by the method proposed by Debussche [21]. The outline
of this method is as follows.
Firstly, it is verified that
Sε
(
ω
)is almost surely uniformly differentiable on
A
(
w
), i.e. for
uA
(
w
),
there exist a linear operator
DSε
(
ω, u
)in
L
(
H
), the space of continuous linear operator from
H
to
H
,
such that if uand u+hare in A(w):
|Sε(ω)(u+h)Sε(ω)uDSε(ω, u)h|≤K(ω)|h|1+α,(25)
where K(w)1,w is a random variable, and α>0is a constant.
Secondly, there exists an integrable random variable ωd, such that
ωd(DSε(ω, u)) ωd(ω),uA(ω),Pa.s, (26)
and
E(ln (ωd)) <0.(27)
Thirdly, there exists a random variable α1such Palmost surely
α11
1(DSε(ω, u)) α1(ω),uA(ω),Pa.s, (28)
and
E(ln α1)<.(29)
Finally,
E(ln K)<.(30)
If the (26),(27),(28)(29) and (30) hold, the Hausdorff dimension of A(w)is less than d.
The main results in this paper is given in the following Theorem.
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
394
Lemma 9. Let Sε(ω) := Sε(t, ω), then Sε(ω)is uniformly quasidifferentiable on A(ω).
The following relationship between Hausdorff dimension and global Lyapunov exponents which will
be used to assert the global stochastic stability.
Theorem 2. [21]Suppose that, for
k
=1
,...,d, sup
uA(ω)
ωk
(
DSε
(
ω, u
)) is integrable. Then we can
generalize the notion of global Lyapunov exponents introduced in Ref. [17] by setting
Λk=Eln sup
uA(ω)
ωk(DSε(ω, u))Eln sup
uA(ω)
ωk1(DSε(ω, u)),
for k2and
Λ1=Eln sup
uA(ω)
ω1(DSε(ω, u)).
And it is easy to see that there exists ωdsatisfying (26)(27) if and only if
Λ1+...
d<0.
3 MAIN RESULTS AND PROOFS
By Lemma 4, we get that
Sε
(
t, ω
)has a global absorbing set, which along with Lemma 5, Lemma 6 and
Lemma 7 gives that
Sε
(
t, ω
)is asymptotically compact. Thus, employing the Theorem 1, it is asserted
that there exist global random attractors A(ω) for system 13.
Based on the results and the estimations on global random attractors, the Hausdorff dimension of
the global random attractors
A
(
ω
)can be got by the method proposed by Debussche [21]. The outline
of this method is as follows.
Firstly, it is verified that
Sε
(
ω
)is almost surely uniformly differentiable on
A
(
w
), i.e. for
uA
(
w
),
there exist a linear operator
DSε
(
ω, u
)in
L
(
H
), the space of continuous linear operator from
H
to
H
,
such that if uand u+hare in A(w):
|Sε(ω)(u+h)Sε(ω)uDSε(ω, u)h|≤K(ω)|h|1+α,(25)
where K(w)1,w is a random variable, and α>0is a constant.
Secondly, there exists an integrable random variable ωd, such that
ωd(DSε(ω, u)) ωd(ω),uA(ω),Pa.s, (26)
and
E(ln (ωd)) <0.(27)
Thirdly, there exists a random variable α1such Palmost surely
α11
1(DSε(ω, u)) α1(ω),uA(ω),Pa.s, (28)
and
E(ln α1)<.(29)
Finally,
E(ln K)<.(30)
If the (26),(27),(28)(29) and (30) hold, the Hausdorff dimension of A(w)is less than d.
The main results in this paper is given in the following Theorem.
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
Theorem 3. Let
d= min
nZ+|1
n
n
i=1
λ1
2
i<
ε
28
εβr2+|p|22M7
πµ 2M8
2µ
16β2r4
ελ
1
4
1
,
when
M7,M
8
are constants described by
(56)
stated in Section 4, then the Hausdorff dimension of
A
(
ω
)
dH(A(ω)) < d.
Obviously, employing and Theorem 2, we have that when
λ1
2
1
ε
28
εβr2+|p|22M7
πµ 2M8
2µ
16β2r4
ελ
1
4
1
,
the Hausdorff dimension of for system
(13)
is 0, which indicates that the global random attractors
consists of only one random fixed point which is global stability.
Proof. The solution mapping of system
(20)
denoted by
DSε
(
t, ω, U
), and
DS
(
t, ω, V
)signifies the
solution mapping of system
(21)
. Moreover, let
DSε
(
ω, U
) :=
DSε
(1
, ω, U
). By Lemma 9, we attain
Sε
(
ω
)is almost surely uniform differentiable, and the conditions
(25)
and
(30)
hold. In the light of
(90)
and
(91)
stated in the proof of Lemma 9 in Section 4, the
(28)
and
(29)
are satisfied. On the other
hand, we have
ωn(DS(ω, V)) = sup
U(i)
0
<1
i=1,...,n
exp 1
tt
0
Tr(Q+
X2(V)+
X3(V)) Qn(s)ds,
where
Qn
(
s
)=
Qns, τ, Vτ;
V0
1,...,
V0
n
is the orthogonal projector from
E1
onto the space spanned
by
V1
(
t
)
,...,
Vn
(
t
), here
V
=
S
(
t, τ
)
Vτ
and
V1,...,
Vn
are the solution of system
(21)
with initial
values
V0
=
V0
1,...,
V0
n
respectively. For any given time
s
,
Vi
(
s
)=
{µi(s)
i(τ)},i
=1
,...,n
is an
orthonormal basis of Qn(s)E1, then
Tr(Q+
X2(V)+
X3(V)) Qn(s)=
n
i=1 (Q+
X2(V)+
X3(V))Vi(s),Vi(s)E1
.
In term of Lemma 3, we get
(QV i,Vi)E1≤−ε
2Vi2
E1ε
4νi2,
and
X3Vi,ViE1
=�β∥∇u2p(∆)µi
i(2β(u, µi)(∆)u, νi)
β∥∇u2p(∆)µi
νi+2β(u, µi)(∆)u∥∥νi
≤|β∥∇u2p|∥µi∥∥νi+2β
λ
1
8
1u2∥∇µi∥∥νi
8
εβr2+|p|2µi2+16β2r4
ελ
1
4
1∥∇µi2+ε
4νi2,
(31)
moreover, we have
X2
V,
VM7|z(θtω)|+M8|z(θtω)|2
V1
2
H2+
V2
2
H2,(32)
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
395
where M7,M
8defined by (57) stated in Section 4. Hence, we can obtain
Tr(Q+
X2(V)+
X3(V)) Qn(s)
≤−ε
2Vi2
E1+8
εβr2+|p|2µi2+16β2r4
ελ
1
4
1∥∇µi2
+M7|z(θtω)|+M8|z(θtω)|2
V1
2
H2+
V2
2
H2
≤−
2+8n
εβr2+|p|2+16β2r4
ελ
1
4
1
n
i=1
λ1
2
i
+2nM7|z(θtω)|+M8|z(θtω)|2.
Thus,
ωn(DSε(ω, V)) = sup
U(i)
0
<1
i=1,...,n
exp 1
tt
0
Tr(Q+
X2(V)+
X3(V)) Qn(s)ds
= exp
1
tt
0
2+8n
εβr2+|p|2+16β2r4
ελ
1
4
1
n
i=1
λ1
2
i+2nM7|z(θtω)|+M8|z(θtω)|2ds
,
which together with (9) gives that T4(ω)>0,tT4,
ωn(DSε(ω, V))
exp
2+8n
εβr2+|p|2+2nM7
πµ +2nM8
2µ+16β2r4
ελ
1
4
1
n
i=1
λ1
2
i
.(33)
Let
T(θtω)=10
σz (θtω)1
.
Clearly, T(θtω)is a linear operator from E1to itself, then we have
DSε(t, ω)=T(θtω)DS(T,ω).
Let
On
be the space spanned by
e1,...,e
n
for any
nN
, then the quadratic form
χOn→
T
(
θtω
)
χ2
E1
is well defined, continuous, and nonnegative on
On
. Let
α1···αn
be the eigenvalues
associated with χ1
2,...,χ
nOnsatisfying
(T(θtω)χi,T(θtω)χj)E1=α2
iδij .
Set χj=(ξj
j)T,j =1,2,...,n, it can be derived
α2
i=(T(θtω)χi,T(θtω)χi)E1
(ξj
j)H2+(σz (θtω)ξj+ηj, σz (θtω)ξj+ηj)
(ξj
j)H2+(ηj
j)+2σ|z(θtω)|∥ξj∥∥ηj+σ2|z(θtω)|2(ξj
j)
≤∥ξj2
H2+ηj2+1
λ
1
2
1σ|z(θtω)|+σ2|z(θtω)|2ξj2
H2+σ|z(θtω)|∥ηj2
1+ σ
λ
1
2
1|z(θtω)|+σ2
λ
1
2
1|z(θtω)|2+λ1
1+λ1
1σ|z(θtω)|.
So
ln αiln
σ
2λ
1
2
1|z(θtω)|+σ
2λ
1
2
1|z(θtω)|2+1
2λ1
+σ
2λ1|z(θtω)|
,
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
396
where M7,M
8defined by (57) stated in Section 4. Hence, we can obtain
Tr(Q+
X2(V)+
X3(V)) Qn(s)
≤−ε
2Vi2
E1+8
εβr2+|p|2µi2+16β2r4
ελ
1
4
1∥∇µi2
+M7|z(θtω)|+M8|z(θtω)|2
V1
2
H2+
V2
2
H2
≤−
2+8n
εβr2+|p|2+16β2r4
ελ
1
4
1
n
i=1
λ1
2
i
+2nM7|z(θtω)|+M8|z(θtω)|2.
Thus,
ωn(DSε(ω, V)) = sup
U(i)
0
<1
i=1,...,n
exp 1
tt
0
Tr(Q+
X2(V)+
X3(V)) Qn(s)ds
= exp
1
tt
0
2+8n
εβr2+|p|2+16β2r4
ελ
1
4
1
n
i=1
λ1
2
i+2nM7|z(θtω)|+M8|z(θtω)|2ds
,
which together with (9) gives that T4(ω)>0,tT4,
ωn(DSε(ω, V))
exp
2+8n
εβr2+|p|2+2nM7
πµ +2nM8
2µ+16β2r4
ελ
1
4
1
n
i=1
λ1
2
i
.(33)
Let
T(θtω)=10
σz (θtω)1
.
Clearly, T(θtω)is a linear operator from E1to itself, then we have
DSε(t, ω)=T(θtω)DS(T,ω).
Let
On
be the space spanned by
e1,...,e
n
for any
nN
, then the quadratic form
χOn
T
(
θtω
)
χ2
E1
is well defined, continuous, and nonnegative on
On
. Let
α1···αn
be the eigenvalues
associated with χ1
2,...,χ
nOnsatisfying
(T(θtω)χi,T(θtω)χj)E1=α2
iδij .
Set χj=(ξj
j)T,j =1,2,...,n, it can be derived
α2
i=(T(θtω)χi,T(θtω)χi)E1
(ξj
j)H2+(σz (θtω)ξj+ηj, σz (θtω)ξj+ηj)
(ξj
j)H2+(ηj
j)+2σ|z(θtω)|∥ξj∥∥ηj+σ2|z(θtω)|2(ξj
j)
≤∥ξj2
H2+ηj2+1
λ
1
2
1σ|z(θtω)|+σ2|z(θtω)|2ξj2
H2+σ|z(θtω)|∥ηj2
1+ σ
λ
1
2
1|z(θtω)|+σ2
λ
1
2
1|z(θtω)|2+λ1
1+λ1
1σ|z(θtω)|.
So
ln αiln
σ
2λ
1
2
1|z(θtω)|+σ
2λ
1
2
1|z(θtω)|2+1
2λ1
+σ
2λ1|z(θtω)|
,
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
for any 1in. Since
ωn(T(θtω)) = α1α2...α
n,
we obtain
ln ωn(T(θtω)) ln
2λ
1
2
1|z(θtω)|+
2λ
1
2
1|z(θtω)|2+n
2λ1
+
2λ1|z(θtω)|
,(34)
which combine with Lemma 1 gives that
E(ωn(T(θtω))) E
2λ
1
2
1|z(θtω)|+
2λ
1
2
1|z(θtω)|2+n
2λ1
+
2λ1|z(θtω)|
M9,
(35)
where M9given by (58) stated in Section 4. Since (33) and (34), we get that
ln ωn(DSε(ω, U)) = ln ωn(T(θtω)) + ln ωn(DS(T,ω))
ln
2λ
1
2
1|z(θtω)|+
2λ
1
2
1|z(θtω)|2+n
2λ1
+
2λ1|z(θtω)|
2+8n
εβr2+|p|2+2nM7
πµ +2nM8
2µ+16β2r4
ελ
1
4
1
n
i=1
λ1
2
i,
then
ωn(DSε(ω, U))
2λ
1
2
1|z(θtω)|+
2λ
1
2
1|z(θtω)|2+n
2λ1
+
2λ1|z(θtω)|
×exp
2+8n
εβr2+|p|2+2nM7
πµ +2nM8
2µ+16β2r4
ελ
1
4
1
n
i=1
λ1
2
i
.
Employing (35), we obtain
E(ωn(DSε(ω, U)))
M9E
exp
2+8n
εβr2+|p|2+2nM7
πµ +2nM8
2µ+16β2r4
ελ
1
4
1
n
i=1
λ1
2
i
.
Let
ωn(ω)=M9exp
2+8n
εβr2+|p|2+2nM7
πµ +2nM8
2µ+16β2r4
ελ
1
4
1
n
i=1
λ1
2
i
,
then
E(ωn(DSε(ω, U))) E(ωn(ω)) ,(36)
and
E(ln ωn(ω)) ln M9
×E
2+8n
εβr2+|p|2+2nM7
πµ +2nM8
2µ+16β2r4
ελ
1
4
1
n
i=1
λ1
2
i
,(37)
together with Lemma 4 we find
ln M9×E
2+8n
εβr2+|p|2+2nM7
πµ +2nM8
2µ+16β2r4
ελ
1
4
1
n
i=1
λ1
2
i
+.
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
397
Therefore, we get dH(A(ω)) < d. Especially, if
λ1
2
1
ε
28
εβr2+|p|22M7
πµ 2M8
2µ
16β2r4
ελ
1
4
1
,
dH
(
A
(
ω
))=0, which merges with the Theorem 2 shows that the largest global Lyapunov exponent of
Sε(ω)is
λ1=E(ln ω1(ω))
E
ε
2+8
εβr2+|p|2+2M7
πµ +2M8
2µ+16β2r4
ελ
1
4
1
n
i=1
λ1
2
i
<0,(38)
thus, we can conclude there exists a random fixed point of system
(14)
which is global stochastic
stability.
4 PROOFS FOR LEMMAS
This section is intended to complete the proofs of Lemmas listed in subsection 2.2.
Proof for Lemma 2: Firstly, we display
X2
(
U
)+
X3
(
U
):
E1E1
satisfies local Lipschitz
condition. Since X3(U)X3(V)+X2(U, θtω)X2(V,θtω)E1
≤∥X3(U)X3(V)E1+X2(U, θtω)X2(V, θtω)E1,
let
c>
0
R
are given constant, for
U,VE1,UE1c, VE1c
, combining Lemmas 4 and
7, we get there exists a positive constant C1(T,c)such that
X3(U)X3(V)E1
=
β∥∇U12p((∆)U1(∆)V1)+β∥∇U12β∥∇V12(∆)V1
β∥∇U12pU1V1+β(∥∇U1+β∥∇V1)V1∥∥U1V1
C1(T,c)UVE1.
On the other hand, by Lemma 1, there exists a positive constant C2(T,c)which satisfies
X2(U, θtω)X2(V,θtω)E1
=|σz (θtω)|∥U1V1H2+|σz (θtω)|∥V2U2
+σ(µ+2εσz (θtω)) z(θtω)(U1V1)+|σαz (θtω)|
A1
2V1A1
2U1
≤|σz (θtω)|
A1
2V1A1
2U1
+|σαz (θtω)|
A1
2V1A1
2U1
+|σ(µ+2εσz (θtω)) z(θtω)|
λ
1
2
1
A1
2V1A1
2U1
+|σz (θtω)|∥V2U2
C2(T,c)UVE1,
The indicated above conclude X2(U)+X3(U):E1E1satisfies local Lipschitz condition.
The semigroup method (Theorem 2.5.4 in [33]) is employed to achieve the existence and uniqueness
of solution for system
(14)
. Based on Lemma 3.5 in Ref. [24] and Lemma 2.2.3 in Ref. [33], we have
Q
is
m
-accretive in
E1
, then it can induce a linear semigroup of contractions formulated by
eQt,tR+,
which
together with the assertion that
X2(U, θtω)
+
X3
(
U
):
E1E1
satisfies local Lipschitz condition can
guarantee the system (14) possesses a unique local mild solution with the form
V(t;Vτ)=eQ(tτ)Vτ+t
τ
eQ(ts)(X2(θsω)+X3(V)(s)) ds,
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
398
Therefore, we get dH(A(ω)) < d. Especially, if
λ1
2
1
ε
28
εβr2+|p|22M7
πµ 2M8
2µ
16β2r4
ελ
1
4
1
,
dH
(
A
(
ω
))=0, which merges with the Theorem 2 shows that the largest global Lyapunov exponent of
Sε(ω)is
λ1=E(ln ω1(ω))
E
ε
2+8
εβr2+|p|2+2M7
πµ +2M8
2µ+16β2r4
ελ
1
4
1
n
i=1
λ1
2
i
<0,(38)
thus, we can conclude there exists a random fixed point of system
(14)
which is global stochastic
stability.
4 PROOFS FOR LEMMAS
This section is intended to complete the proofs of Lemmas listed in subsection 2.2.
Proof for Lemma 2: Firstly, we display
X2
(
U
)+
X3
(
U
):
E1E1
satisfies local Lipschitz
condition. Since X3(U)X3(V)+X2(U, θtω)X2(V,θtω)E1
≤∥X3(U)X3(V)E1+X2(U, θtω)X2(V, θtω)E1,
let
c>
0
R
are given constant, for
U,VE1,UE1c, VE1c
, combining Lemmas 4 and
7, we get there exists a positive constant C1(T,ω,c)such that
X3(U)X3(V)E1
=
β∥∇U12p((∆)U1(∆)V1)+β∥∇U12β∥∇V12(∆)V1
β∥∇U12pU1V1+β(∥∇U1+β∥∇V1)V1∥∥U1V1
C1(T,ω,c)UVE1.
On the other hand, by Lemma 1, there exists a positive constant C2(T,ω,c)which satisfies
X2(U, θtω)X2(V,θtω)E1
=|σz (θtω)|∥U1V1H2+|σz (θtω)|∥V2U2
+σ(µ+2εσz (θtω)) z(θtω)(U1V1)+|σαz (θtω)|
A1
2V1A1
2U1
≤|σz (θtω)|
A1
2V1A1
2U1
+|σαz (θtω)|
A1
2V1A1
2U1
+|σ(µ+2εσz (θtω)) z(θtω)|
λ
1
2
1
A1
2V1A1
2U1
+|σz (θtω)|∥V2U2
C2(T,ω,c)UVE1,
The indicated above conclude X2(U)+X3(U):E1E1satisfies local Lipschitz condition.
The semigroup method (Theorem 2.5.4 in [33]) is employed to achieve the existence and uniqueness
of solution for system
(14)
. Based on Lemma 3.5 in Ref. [24] and Lemma 2.2.3 in Ref. [33], we have
Q
is
m
-accretive in
E1
, then it can induce a linear semigroup of contractions formulated by
eQt,tR+,
which
together with the assertion that
X2(U, θtω)
+
X3
(
U
):
E1E1
satisfies local Lipschitz condition can
guarantee the system (14) possesses a unique local mild solution with the form
V(t;Vτ)=eQ(tτ)Vτ+t
τ
eQ(ts)(X2(θsω)+X3(V)(s)) ds,
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
where, tτ, t, τ R.
Proof for Lemma 3 : Since
(QU ,U)E1
≤−εU12
H2+εU22+εα
A1
2U1
U2∥−ε2U1∥∥U2∥−α
A1
4U2
2
≤−εU12
H2+εU22+ε
2
13p
2λ
1
4
1
A1
2U1
2+εα2λ
1
4
1
2λ
1
4
13pU22α
A1
4U2
2
≤−ε
2U12
H2+εU22+εα2λ
1
4
1
2λ
1
4
13pU22α
A1
4U2
23
4A1
4U12
≤−ε
2U12
H2+
α
λ
1
4
1
+ε+εα2λ
1
4
1
2λ
1
4
13p
U223
4A1
4U12
≤−ε
2U2
E1ε
4U223
4A1
4U12,
where, pand εsatisfies (15). Thus complete the proof.
Before deriving the other proofs for Lemmas, the following estimations and quantities are introduced
2β
A1
4u
2pA1
2u, v≤−1
2β
d
dt βA1
4u2p2
+εp2
4β+3
2A1
4u2ε
4β[βA1
4u2p]2
+4σ2|z(θtω)|2
7εβ βA1
4u2p2,
2σA1
2uz (θtω),A1
2uε
2A1
2u2+2σ2|z(θtω)|2
εA1
2u2,
(39)
2σµ+2εαA1
2σz (θtω)uz (θtω)σvz (θtω),v
2
εσ2|z(θtω)|2v2+ε
2v2+σ2|z(θtω)|24
εu2+ε
4v2
+ε
2v2+4σ2|z(θtω)|2(µ+2ε)2
εu2+α2
εA1
2u2
(40)
2σA1
2uaz(θtω),A1
2uaε
2A1
2ua2+2σ2|z(θtω)|2
εA1
2ua2,(41)
2σ(µ+2εαA1
2σz (θtω))uaz(θtω)σvaz(θtω),A1
2va
2
εσ2|z(θtω)|2va2+ε
2va2+σ2|z(θtω)|24
εua2+ε
4va2
+ε
4va2+8σ2|z(θtω)|2(µ+2ε)2
εua2+α2
εA1
2ua2,
(42)
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
399
(2βA1
4u2pA1
2ua,A1
2vb)
3β2
λ
1
2
1A1
2u4
A1
2ua
2+3p2Va2
E1
+2d
dt
A1
2ub
2+2ε2
A1
2ub
2+2σ2|z(θtω)|2
A1
2ub
2,
2βA1
4u2p(A1
2ub,A1
2vb)
=2βA1
4u2p(A1
2ub,A1
2ut+εubσuz (θtω))
≤−d
dt βA1
4u2pA1
2ub2+β
2A1
2u3ut,
(43)
2σµ+2εαA1
2σz (θtω)ubz(θtω)σvz (θtω),A1
2vb
ε
2A3
4ub2+2σ2|z(θtω)|2
εA1
4ub2+2
εσ2|z(θtω)|2A1
4vb2+ε
2A1
4vb2
+ε
2A1
4vb2+4σ2|z(θtω)|2(µ+2ε)2
εA1
4ub2+α2
εA3
4ub2
+σ2|z(θtω)|24
εA1
4ub2+ε
4A1
4vb2,
(44)
and β
2A1
2u3ut∥≤β
2V3
E1(v+εu+σuz (θtω))
β
2V3
E1(VE1+εVE1+σ|z(θtω)|∥VE1)
β+βε
2+β
2|σz (θtω)|V4
E1,
By Lemma 4, Lemma 5 and Lemma7,
M(t)defined as follow is bounded,
M(t)=3β2
λ
1
2
1A1
2u4
A1
2ua
2+3p2Va2
E1+β
2A1
2u3ut
+2ε2+3
2+2ε+2|p|A1
2ub2<+,
(45)
M0= max
8α2+2
ε+8(µ+2ε)2+4
ελ
1
2
1
,2
ε+ε
4
,(46)
M= max
8(µ+2ε)2+6
ελ
1
2
1
+8α2
ε,ε2+8
4ε,8
7ε
,(47)
M1=4βr2
1(ω)+3|p|λ
1
4
1
2λ
1
4
1
,(48)
M2=3βr4
1(ω)
λ
1
4
1
+2|p|r2
1(ω)+εp2
4β,(49)
M3= max
1,2
ε+ε
4,4(µ+2ε)2+2
ελ
1
2
1
+4α2
ε+2|p|+4
λ
1
4
1
+ε
4λ
1
2
1
,(50)
c1(ω) = max
1,6βr1(ω)
λ
1
4
1
,c
2(ω) = max
1,r2(ω)
2+2βr1(ω)
λ
1
8
1
,(51)
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
400
(2βA1
4u2pA1
2ua,A1
2vb)
3β2
λ
1
2
1A1
2u4
A1
2ua
2+3p2Va2
E1
+2d
dt
A1
2ub
2+2ε2
A1
2ub
2+2σ2|z(θtω)|2
A1
2ub
2,
2βA1
4u2p(A1
2ub,A1
2vb)
=2βA1
4u2p(A1
2ub,A1
2ut+εubσuz (θtω))
≤−d
dt βA1
4u2pA1
2ub2+β
2A1
2u3ut,
(43)
2σµ+2εαA1
2σz (θtω)ubz(θtω)σvz (θtω),A1
2vb
ε
2A3
4ub2+2σ2|z(θtω)|2
εA1
4ub2+2
εσ2|z(θtω)|2A1
4vb2+ε
2A1
4vb2
+ε
2A1
4vb2+4σ2|z(θtω)|2(µ+2ε)2
εA1
4ub2+α2
εA3
4ub2
+σ2|z(θtω)|24
εA1
4ub2+ε
4A1
4vb2,
(44)
and β
2A1
2u3ut∥≤β
2V3
E1(v+εu+σuz (θtω))
β
2V3
E1(VE1+εVE1+σ|z(θtω)|∥VE1)
β+βε
2+β
2|σz (θtω)|V4
E1,
By Lemma 4, Lemma 5 and Lemma7,
M(t)defined as follow is bounded,
M(t)=3β2
λ
1
2
1A1
2u4
A1
2ua
2+3p2Va2
E1+β
2A1
2u3ut
+2ε2+3
2+2ε+2|p|A1
2ub2<+,
(45)
M0= max
8α2+2
ε+8(µ+2ε)2+4
ελ
1
2
1
,2
ε+ε
4
,(46)
M= max
8(µ+2ε)2+6
ελ
1
2
1
+8α2
ε,ε2+8
4ε,8
7ε
,(47)
M1=4βr2
1(ω)+3|p|λ
1
4
1
2λ
1
4
1
,(48)
M2=3βr4
1(ω)
λ
1
4
1
+2|p|r2
1(ω)+εp2
4β,(49)
M3= max
1,2
ε+ε
4,4(µ+2ε)2+2
ελ
1
2
1
+4α2
ε+2|p|+4
λ
1
4
1
+ε
4λ
1
2
1
,(50)
c1(ω) = max
1,6βr1(ω)
λ
1
4
1
,c
2(ω) = max
1,r2(ω)
2+2βr1(ω)
λ
1
8
1
,(51)
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
C(p)=
1,p2,
(1 p+2
λ
1
4
1
),2<p<λ
1
4
12
2,(52)
q=
σµ+2ε+2λ
1
2
1+αλ
1
2
1
2λ
1
2
1
,(53)
M4=ε+c2(ω),(54)
M5=
σµ+2ε+αλ
1
4
1
2λ
1
4
1
+|σ|,(55)
M6=ε16
εβr2+|p|264β2r4
ελ
1
2
1
,(56)
M7=σ+σ
λ
1
2
1
+
σµ+2ε+αλ
1
2
1
2λ
1
2
1
,M
8=σ2
2λ
1
2
1
,(57)
M9=
2λ1πµ +
2λ1πµ +
4µλ
1
2
1
+n
2λ1
.(58)
Proof for Lemma 4 : Taking the inner product of Vby V=[u, v]Tin E1, we get that
d
dtV2
E12(QV ,V)E1+2σA1
2uz (θtω),A1
2u
+2β
A1
2u
2pA1
2u, v
+2σµ+2εαA 1
2σz (θtω)uz (θtω)σvz (θtω),v
.
(59)
Set
H(t)=H(u, v)= 1
2ββ
A1
4u
2p2
+
A1
2u
2+v2.
By (40) and Lemma 3, we find that
d
dtH(u, v)≤−ε
2H(u, v)+εp2
4β+4σ2|z(θtω)|2
7εβ βA1
4u2p2
+
α2+2
ε+(µ+2ε)2+4
ελ
1
2
1
σ2|z(θtω)|2A1
2u2
+2
ε+ε
4σ2|z(θtω)|2v2
≤−ε
2H(u, v)+2|z(θtω)|2H(u, v)+ εp2
4β,
(60)
where Mis formulated by (47). Then, for any t>0, the following holds
H(0,t, ω)eε
2t+0
t2|z(θsω)|2dsH(t, t, ω)+εp2
4β0
t
eε
2s+0
s2|z(θkω)|2dk ds.
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
401
Since v=ut+εu σuz (θtω), we get
U2
E12
A1
2u
2+2
ut+εu σuz (θtω)2+2σuz (θtω)
2
2
1+ σ2
λ
1
4
1|z(θtω)|2
V2
E1.
(61)
Thus, we have
U(0,t, ω)2
E12
1+ σ2
λ
1
4
1|z(θtω)|2
V(0,t, ω)2
E1
2
1+ σ2
λ
1
4
1|z(θtω)|2
eε
2t+0
t2|z(θsω)|2dsH(t, t, ω)
+
1+ σ2
λ
1
4
1|z(θtω)|2
εp2
2β0
t
eε
2s+0
s2|z(θkω)|2dk ds.
(62)
Since the random variable z(θtω)is tempered, along with Lemma 1, we can infer that
2σ2
λ
1
4
1|z(θtω)|2eε
2t+0
t2|z(θsω)|2dsH(t, t, ω)0,t+,
and 2σ2
λ
1
4
1|z(θtω)|2εp2
4β0
t
eε
2s+0
s2|z(θkω)|2dk ds 0,t+,
similarly
2eε
2t+0
t2|z(θsω)|2dsH(t, t, ω)0,t+,
and εp2
2β0
t
eε
2s+0
s2|z(θkω)|2dk ds < .
Let
ρ1(ω)=εp2
2β0
t
eε
2s+0
s2|z(θkω)|2dk ds. (63)
According to (59) - (63), we obtain that there exists T1(ω)>0, such that
U(t, θtω)E1=U(0,t, ω)E1r1(ω),tT1(ω),(64)
where r1(ω)=ρ1(ω).
On the other hand, utilizing (10), we obtain
E(ρ1(ω)) = εp2
2β0
−∞
eε
2sEe20
s|z(θkω)|2dkds
εp2
2β0
−∞
eε
2ssMσ2
µds
=εp2
2β0
T
eε
2ssMσ2
µds εp2
2βT
−∞
eε+2
µsds
=εp2
2β0
T
eε
2ssMσ2
µds εp2eε+2
µT
2βε+2
µ
<+.
(65)
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
402
Since v=ut+εu σuz (θtω), we get
U2
E12
A1
2u
2+2
ut+εu σuz (θtω)2+2σuz (θtω)
2
2
1+ σ2
λ
1
4
1|z(θtω)|2
V2
E1.
(61)
Thus, we have
U(0,t, ω)2
E12
1+ σ2
λ
1
4
1|z(θtω)|2
V(0,t, ω)2
E1
2
1+ σ2
λ
1
4
1|z(θtω)|2
eε
2t+0
t2|z(θsω)|2dsH(t, t, ω)
+
1+ σ2
λ
1
4
1|z(θtω)|2
εp2
2β0
t
eε
2s+0
s2|z(θkω)|2dk ds.
(62)
Since the random variable z(θtω)is tempered, along with Lemma 1, we can infer that
2σ2
λ
1
4
1|z(θtω)|2eε
2t+0
t2|z(θsω)|2dsH(t, t, ω)0,t+,
and 2σ2
λ
1
4
1|z(θtω)|2εp2
4β0
t
eε
2s+0
s2|z(θkω)|2dk ds 0,t+,
similarly
2eε
2t+0
t2|z(θsω)|2dsH(t, t, ω)0,t+,
and εp2
2β0
t
eε
2s+0
s2|z(θkω)|2dk ds < .
Let
ρ1(ω)=εp2
2β0
t
eε
2s+0
s2|z(θkω)|2dk ds. (63)
According to (59) - (63), we obtain that there exists T1(ω)>0, such that
U(t, θtω)E1=U(0,t, ω)E1r1(ω),tT1(ω),(64)
where r1(ω)=ρ1(ω).
On the other hand, utilizing (10), we obtain
E(ρ1(ω)) = εp2
2β0
−∞
eε
2sEe20
s|z(θkω)|2dkds
εp2
2β0
−∞
eε
2ssMσ2
µds
=εp2
2β0
T
eε
2ssMσ2
µds εp2
2βT
−∞
eε+2
µsds
=εp2
2β0
T
eε
2ssMσ2
µds εp2eε+2
µT
2βε+2
µ
<+.
(65)
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
Proof for Lemma 5 : Taking the inner product (·,·)E1of (18) with A1
2Vagives
d
dtA1
4Va2
E1=2(QV a,A1
2Va)E1+ 2(σA1
2uaz(θtω), Aua)
+2 σ(µ+2εαA1
2σz (θtω))uaz(θtω)σvaz(θtω),A1
2va.
(66)
Substituting (41) and (42) into (66)„ we infer the following by combining with Lemma 3
d
dtA1
4Va2
E1≤−ε
2A1
4Va2
E1+4
ε+ε
4σ2|z(θtω)|2A1
4va2
+
8(µ+2ε)2
ελ
1
2
1
+8α2
ε+4
ελ
1
2
1
+2
ελ
1
2
1
σ2|z(θtω)|2A3
4ua2,
hence d
dtA1
4Va2
E1ε
2+σ2|z(θtω)|2M0A1
4Va2
E1,
where M0is given by (46). Thus, we have
A1
4Va2
E1e0
tε
2t+σ2|z(θsω)|2M0dsA1
4Va(t, t, ω)2
E1.
Since the random variable z(θtω)is tempered, along with Lemma 1 states
e0
tε
2t+σ2|z(θsω)|2M0dsA1
4Va(t, t, ω)2
E10,t+,
Thus, it can be obtained that
lim
t+
A1
4Ua
(t, θtω)E1=0.
Proof for Lemma 6 : Taking the inner product (·,·)E1of (19) with Vbgives
d
dtVb2
E1Υ(1)
(2),(67)
here
Υ(1) =2(QV b,Vb)+2σA1
2ubz(θtω),A1
2ub+2β
A1
4ub
2pA1
2ub,v
b
+2 σµ+2εαA1
2σz (θtω)ubz(θtω)σvz (θtω),v
b,
Υ(2) =2 β
A1
4ub
2pA1
2ub,v
b+2β
A1
4u
2pA1
2u, vb,
(68)
in which Υ(1) is bounded. On the other hand
Υ(2) =2β
A1
4uA1
4ua
2pA1
2u, 1
2vb+2β
A1
4uA1
4ua
2pA1
2ua,1
2vb
+2β
A1
4u
2pA1
2u, vb
β
λ
1
4
1
2V2+2Va2+|p|λ
1
4
1
β
(V2+Va2)+
β
λ
1
4
1V2+|p|
V2
+β
2λ
1
4
1
4V2+2Va2+3|p|λ
1
4
1
β
H(ub,v
b).
(69)
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
403
Set
H(t)=H(ub,v
b)= 1
2ββ
A1
4ub
2p2
+
A1
2ub
2+vb2,
taking account into (67), (68) and (69), we can get
d
dtH(ub,v
b)≤−ε
2H(ub,v
b)+
β
λ
1
4
1V2+|p|
V2+εp2
4β
+
2|z(θtω)|2+2β
λ
1
4
1V2+β
λ
1
4
1Va2+3|p|
H(ub,v
b)
+2βV2+2βVa2+|p|λ
1
4
1
λ
1
4
1
(V2+Va2).
(70)
Thus,
H(0,t, ω)eε
2t+0
t2|z(θsω)|2+β
2λ
1
4
1
4V2+2Va2+3|p|λ
1
4
1
β
ds
H(t, t, ω)
+0
t
eε
2s+0
s2|z(θkω)|2+β
2λ
1
4
1
4V2+2Va2+3|p|λ
1
4
1
β
dk
×
2βV2+2βVa2+|p|λ
1
4
1
λ
1
4
1
(V2+Va2)+βV4
λ
1
4
1
+|p|∥V2+εp2
4β
ds.
Applying Lemma 4 and Lemma 5, we obtain
H(0,t, ω)eε
2t+0
t2|z(θsω)|2+M1dsH(t, t, ω)
+M20
t
eε
2s+0
s2|z(θkω)|2+M1dk ds,
where
M1,M
2
are defined by
(48)
and
(49)
correspondingly. Since the random variable
z
(
θtω
)is
tempered, together with Lemma 1, we find
eε
2t+0
t2|z(θsω)|2+M1dsH(t, t, ω)0,t+,
and
M20
t
eε
2s+0
s2|z(θkω)|2+M1dk ds < .
Let
ρ2(ω)=M20
t
eε
2s+0
s2|z(θkω)|2+M1dk ds. (71)
As indicated above, we obtain that there exists T2(ω)>0such that
UE1(t, θtω)=UE1(0,t, ω)r2(ω),tT2(ω),
where r2(ω)=ρ2(ω).
Proof for Lemma 7 : Taking the inner product (·,·)E1of (19) with A1
2Vbgives
d
dtA1
4Vb2
E12(QV b,A1
2Vb)+2σA1
2ubz(θtω), Aub
+2 β
A1
4u
2pA1
2u, A1
2vb
+2 σµ+2εαA1
2σz (θtω)ubz(θtω)σvz (θtω),A1
2vb.
(72)
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
404
Set
H(t)=H(ub,v
b)= 1
2ββ
A1
4ub
2p2
+
A1
2ub
2+vb2,
taking account into (67), (68) and (69), we can get
d
dtH(ub,v
b)≤−ε
2H(ub,v
b)+
β
λ
1
4
1V2+|p|
V2+εp2
4β
+
2|z(θtω)|2+2β
λ
1
4
1V2+β
λ
1
4
1Va2+3|p|
H(ub,v
b)
+2βV2+2βVa2+|p|λ
1
4
1
λ
1
4
1
(V2+Va2).
(70)
Thus,
H(0,t, ω)eε
2t+0
t2|z(θsω)|2+β
2λ
1
4
1
4V2+2Va2+3|p|λ
1
4
1
β
ds
H(t, t, ω)
+0
t
eε
2s+0
s2|z(θkω)|2+β
2λ
1
4
1
4V2+2Va2+3|p|λ
1
4
1
β
dk
×
2βV2+2βVa2+|p|λ
1
4
1
λ
1
4
1
(V2+Va2)+βV4
λ
1
4
1
+|p|∥V2+εp2
4β
ds.
Applying Lemma 4 and Lemma 5, we obtain
H(0,t, ω)eε
2t+0
t2|z(θsω)|2+M1dsH(t, t, ω)
+M20
t
eε
2s+0
s2|z(θkω)|2+M1dk ds,
where
M1,M
2
are defined by
(48)
and
(49)
correspondingly. Since the random variable
z
(
θtω
)is
tempered, together with Lemma 1, we find
eε
2t+0
t2|z(θsω)|2+M1dsH(t, t, ω)0,t+,
and
M20
t
eε
2s+0
s2|z(θkω)|2+M1dk ds < .
Let
ρ2(ω)=M20
t
eε
2s+0
s2|z(θkω)|2+M1dk ds. (71)
As indicated above, we obtain that there exists T2(ω)>0such that
UE1(t, θtω)=UE1(0,t, ω)r2(ω),tT2(ω),
where r2(ω)=ρ2(ω).
Proof for Lemma 7 : Taking the inner product (·,·)E1of (19) with A1
2Vbgives
d
dtA1
4Vb2
E12(QV b,A1
2Vb)+2σA1
2ubz(θtω), Aub
+2 β
A1
4u
2pA1
2u, A1
2vb
+2 σµ+2εαA1
2σz (θtω)ubz(θtω)σvz (θtω),A1
2vb.
(72)
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
Set
H1(t)=H1(u, v)=βA1
4u2p2A1
2ub2+A1
4Vb2
E1,(73)
we find
H1(t)=A1
4Vb2
E1+βA1
2ub2A1
4u2pA1
2ub22A1
2ub2
≥∥A1
4ub2+βA1
2ub2A1
4u2+ (1 p+2
λ
1
4
1
)A3
4ub2
≥∥A1
4ub2+βA1
2ub2A1
4u2+C(p)A3
4ub2
0,
where C(p)is given by (52), which along with (43), (44) and (72) states that
d
dtH1(u, v)≤−ε
2H1(u, v)+H1(u, v)+
M(t)
+σ2|z(θtω)|2βA1
4u2p2A1
2ub2
+
4(µ+2ε)2+2
ελ
1
2
1
+4α2
ε+2|p|+4
λ
1
4
1
+ε
4λ
1
2
1
σ2|z(θtω)|2A3
4ub2
+2
ε+ε
4σ2|z(θtω)|2A1
4vb2
≤−ε
2H1(u, v)+1+M3σ2|z(θtω)|2H1(u, v)+
M(t),
where
M
(
t
)and
M3
are defined by
(45)
and
(50)
respectively. Then by Lemma 4, Lemma 5 and
Lemma7, the following holds
H1(0,t, ω)eε
2t+0
t(1+M3σ2|z(θsω)|2)dsH1(t, t, ω)
+r0
t
eε
2s+0
s(1+M3σ2|z(θkω)|2)dkds.
Since the random variable z(θtω)is tempered, applying Lemma 1, we find
eε
2t+0
t(1+M3σ2|z(θsω)|2)dsH1(t, t, ω)0,t+,
r0
t
eε
2s+0
s(1+M3σ2|z(θkω)|2)dkds < .
then we get
H1(0,t, ω)ρ30(ω),
where ρ30(ω)=r0
teε
2s+0
s(1+M3σ2|z(θkω)|2)dkds.
Exploiting (61) and (73), we can obtain
A1
4Vb2
E1(0,t)+βA1
2ub2(0,t)A1
4u2(0,t)
ρ30(ω)+ p+2
λ
1
4
1A3
4ub2(0 t)
ρ30(ω)+ p+2
λ
1
4
1A1
4Vb2
E1(0,t),
then
(1 p+2
λ
1
4
1
)A1
4Vb2
E1(0,t)ρ30(ω),
utilizing (52), we get
A1
4Vb2
E1(0,t)λ
1
4
1
λ
1
4
1p2
ρ30(ω).(74)
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
405
By (72)-(74), we obtain that there exists T30(ω)>0, such that
A1
4Vb2
E1(0,t)ρ3(ω),tT30(ω),(75)
where, ρ3(ω)= λ
1
4
1
λ
1
4
1p+2
ρ30(ω). Since the relation between A1
4UbE1and A1
4VbE1,
A1
4Ub2
E1(0,t, ω)
2
1+ σ2
λ
1
4
1|z(θtω)|2
A1
4Vb(0,t, ω)2
E1
2
1
4
1
λ
1
4
1p+2
1+ σ2
λ
1
4
1|z(θtω)|2
0
t
eε
2s+0
s(1+M3σ2|z(θkω)|2)dkds.
(76)
Since the random variable z(θtω)is tempered, we get the following holds by Lemma 1,
22
λ
1
4
1p+2|z(θtω)|20
t
eε
2s+0
s(1+M3σ2|z(θkω)|2)dkds 0,t +,
and
2
1
4
1
λ
1
4
1p+20
t
eε
2s+0
s(1+M3σ2|z(θkω)|2)dkds < .
Thus, there exists a constant T3(ω)>0such that
A1
4Ub
E1
(t, θtω)=
A1
4Ub
E1
(0,t, ω)r3(ω),t>T
3(ω),
Proof for Lemma 8 : For κ∈D(A1
2), we find that
X32 (u1)
X32 (u2)κ
β∥∇u12β∥∇u22(∆)κ
+2β(u1,κ)(∆)u12β(u2,κ)(∆)u2
β(∥∇u1 ∥∇u2)(∥∇u1+∥∇u2)(∆)κ
+2β(u1,κ)∥∥(∆)u1(∆)u2+2β|(u1−∇u2,κ)∥∥(∆)u2
β(∥∇u1 ∥∇u2)(∥∇u1+∥∇u2)(∆)κ
+2β∥∇u1∥∇κ∥∥(∆)u1(∆)u2+2β∥∇u1−∇u2∥∇κ∥∥(∆)u2,
Merging Lemma 4 with (51), we find
X32 (u1)
X32 (u2)κ
β(∥∇u1 ∥∇u2)(∥∇u1+∥∇u2)(∆)κ
β
λ
1
4
1
(u1∥−∥u2)2r1(ω)(∆)κ
+2β
λ
1
4
1
r1(ω)κ∥∥(∆)u1+∆u2+2β
λ
1
4
1u1u2∥∥κr1(ω)
c1(ω)u1u2∥∥κ,u1,u
2∈D(A1
2),
along with Lemma 4 demonstrates E(c1(ω)) ≤∞, hence
X32 (u1)
X32 (u2)
LD(A1
2),L2(D)= sup
κ∈D(A1
2)
X32 (u1)
X32 (u2)κ
κD(A1
2)
c1(ω)
A1
2u1A1
2u2
.
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
406
By (72)-(74), we obtain that there exists T30(ω)>0, such that
A1
4Vb2
E1(0,t)ρ3(ω),tT30(ω),(75)
where, ρ3(ω)= λ
1
4
1
λ
1
4
1p+2
ρ30(ω). Since the relation between A1
4UbE1and A1
4VbE1,
A1
4Ub2
E1(0,t, ω)
2
1+ σ2
λ
1
4
1|z(θtω)|2
A1
4Vb(0,t, ω)2
E1
2
1
4
1
λ
1
4
1p+2
1+ σ2
λ
1
4
1|z(θtω)|2
0
t
eε
2s+0
s(1+M3σ2|z(θkω)|2)dkds.
(76)
Since the random variable z(θtω)is tempered, we get the following holds by Lemma 1,
22
λ
1
4
1p+2|z(θtω)|20
t
eε
2s+0
s(1+M3σ2|z(θkω)|2)dkds 0,t +,
and
2
1
4
1
λ
1
4
1p+20
t
eε
2s+0
s(1+M3σ2|z(θkω)|2)dkds < .
Thus, there exists a constant T3(ω)>0such that
A1
4Ub
E1
(t, θtω)=
A1
4Ub
E1
(0,t, ω)r3(ω),t>T
3(ω),
Proof for Lemma 8 : For κ∈D(A1
2), we find that
X32 (u1)
X32 (u2)κ
β∥∇u12β∥∇u22(∆)κ
+2β(u1,κ)(∆)u12β(u2,κ)(∆)u2
β(∥∇u1 ∥∇u2)(∥∇u1+∥∇u2)(∆)κ
+2β(u1,κ)∥∥(∆)u1(∆)u2+2β|(u1−∇u2,κ)∥∥(∆)u2
β(∥∇u1 ∥∇u2)(∥∇u1+∥∇u2)(∆)κ
+2β∥∇u1∥∇κ∥∥(∆)u1(∆)u2+2β∥∇u1−∇u2∥∇κ∥∥(∆)u2,
Merging Lemma 4 with (51), we find
X32 (u1)
X32 (u2)κ
β(∥∇u1 ∥∇u2)(∥∇u1+∥∇u2)(∆)κ
β
λ
1
4
1
(u1∥−u2)2r1(ω)(∆)κ
+2β
λ
1
4
1
r1(ω)κ∥∥(∆)u1+∆u2+2β
λ
1
4
1u1u2∥∥κr1(ω)
c1(ω)u1u2∥∥κ,u1,u
2∈D(A1
2),
along with Lemma 4 demonstrates E(c1(ω)) ≤∞, hence
X32 (u1)
X32 (u2)
LD(A1
2),L2(D)= sup
κ∈D(A1
2)
X32 (u1)
X32 (u2)κ
κD(A1
2)
c1(ω)
A1
2u1A1
2u2
.
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
On the other hand, since
X32 (u)
LD(A1
2),L2(D)= sup
κ∈D(A1
2)
X32 (u)κ
κD(A1
2)
=β∥∇u2p(∆)κ+2β(u, κ)(∆)u
κ
β∥∇u2p(∆)κ+2β(u, κ)(∆)u
κ
[β∥∇u2p]2+2β
λ
1
8
1u2.
In term of Lemma 4 and (51), we find
X32(u)
LD(A1
2),L2(D)r2(ω)
2+2β
λ
1
8
1
r1(ω)=c2(ω),
and
E(c2(ω)) ≤∞.
Proof for Lemma 9 : Taking the inner product of (23) by V(1) V(2) in E1illustrates
1
2
d
V(1) V(2)
2
E1
dt =QV(1) V(2),V(1) V(2)E1
+X2V(1)X2V(2),V(1) V(2)E1
+(X32 (ω, u1)X32 (ω, u2),v
1v2),
(77)
Merging with (24), we have
X32 (ω, u1)X32 (ω, u2)∥≤c2
A1
2(u1u2)
.(78)
Furthermore, we get
(X32 (ω, u1)X32 (ω, u2),v
1v2)≤∥X32 (ω, u1)X32 (ω, u2)∥∥v1v2
c2(ω)
A1
2(u1u2)
v1v2
c2(ω)
2
V(1) V(2)
2
E1
,
(79)
and X2V(1)X2V(2),V(1) V(2)E1
|σz (θtω)(µ+2εσz (θtω)) |+2λ
1
2
1|σz (θtω)|
2λ
1
2
1
(unum2
H2+vnvm2)
+|ασz (θtω)|
2
A1
2(unum)
2+vnvm2
|σz (θtω)(qσz (θtω)) |+2λ
1
2
1|σz (θtω)|
2λ
1
2
1VnVm2
E1
q|z(θtω)|+σ2|z(θtω)|2
2λ
1
2
1
VnVm2
E1,
(80)
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
407
here qis defined by (53), and
QV(1) V(2),V(1) V(2)≤−ε
2V(1) V(2)2
E1ε
4v1v22.(81)
Substituting (79), (80) and (81) into (77), we get
d
V(1) V(2)
2
E1
dt
ε+c2(ω)+2q|z(θtω)|+σ2
λ
1
2
1|z(θtω)|2
V(1) V(2)
2
E1
.
thus, we have
V(1) V(2)
2
E1
(0,t, ω)e0
tM4+2q|z(θsω)|+σ2
λ
1
2
1
|z(θsω)|2ds
I2
E1,(82)
taking into account the between
V
2
E1
and
U
2
E1
,
U(1) U(2)
2
E12
1+ σ2
λ
1
4
1|z(θtω)|2
e0
tM4+2q|z(θsω)|+σ2
λ
1
2
1
|z(θsω)|2ds
I2
E1,
here
M4>
0is formulated by
(54)
. Combine assertion that random variable
z
(
θtω
)is tempered with
Lemma 1, we find
2σ2
λ
1
4
1|z(θtω)|2e0
tM4+2q|z(θsω)|+σ2
λ
1
2
1
|z(θsω)|2ds
0,t+,
and
2e0
tM4+2q|z(θsω)|+σ2
λ
1
2
1
|z(θsω)|2ds
<,
then, t
0,tt
0,
U(1) U(2)
2
E1p
1,set p1= max{1,p
1},we obtain
U(1) U(2)
2
E1
(0,t, ω)=
U(1) U(2)
2
E1
(t, θtω)p1IE1.
and
p1(ω)1,E(p1(ω)) <,E(ln (p1(ω))) <.(83)
On the other hand, taking the inner product of
(21)
by
V
in
E1
, which together with Lemma 3 shows
(Q
V,
V)E1≤−ε
2
V2
E1ε
4
V2
2,(84)
and
X3
V,
VE1
=β∥∇u2p(∆)
V1,
V22βu,
V1(∆)u,
V2
≤|β∥∇u2p|
V1
V2
+2β
λ
1
4
1u2
V1
V2
8
ε|β∥∇u2p|
V1
2+8
ε
2β
λ
1
4
1u2
V1
2
+ε
4
V2
2
8
εβr2+|p|2
V1
2+64β2r4
ελ
1
2
1
V1
2+ε
4
V2
2,
(85)
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
408
here qis defined by (53), and
QV(1) V(2),V(1) V(2)≤−ε
2V(1) V(2)2
E1ε
4v1v22.(81)
Substituting (79), (80) and (81) into (77), we get
d
V(1) V(2)
2
E1
dt
ε+c2(ω)+2q|z(θtω)|+σ2
λ
1
2
1|z(θtω)|2
V(1) V(2)
2
E1
.
thus, we have
V(1) V(2)
2
E1
(0,t, ω)e0
tM4+2q|z(θsω)|+σ2
λ
1
2
1
|z(θsω)|2ds
I2
E1,(82)
taking into account the between
V
2
E1
and
U
2
E1
,
U(1) U(2)
2
E12
1+ σ2
λ
1
4
1|z(θtω)|2
e0
tM4+2q|z(θsω)|+σ2
λ
1
2
1
|z(θsω)|2ds
I2
E1,
here
M4>
0is formulated by
(54)
. Combine assertion that random variable
z
(
θtω
)is tempered with
Lemma 1, we find
2σ2
λ
1
4
1|z(θtω)|2e0
tM4+2q|z(θsω)|+σ2
λ
1
2
1
|z(θsω)|2ds
0,t+,
and
2e0
tM4+2q|z(θsω)|+σ2
λ
1
2
1
|z(θsω)|2ds
<,
then, t
0,tt
0,
U(1) U(2)
2
E1p
1,set p1= max{1,p
1},we obtain
U(1) U(2)
2
E1
(0,t, ω)=
U(1) U(2)
2
E1
(t, θtω)p1IE1.
and
p1(ω)1,E(p1(ω)) <,E(ln (p1(ω))) <.(83)
On the other hand, taking the inner product of
(21)
by
V
in
E1
, which together with Lemma 3 shows
(Q
V,
V)E1≤−ε
2
V2
E1ε
4
V2
2,(84)
and
X3
V,
VE1
=β∥∇u2p(∆)
V1,
V22βu,
V1(∆)u,
V2
≤|β∥∇u2p|
V1
V2
+2β
λ
1
4
1u2
V1
V2
8
ε|β∥∇u2p|
V1
2+8
ε
2β
λ
1
4
1u2
V1
2
+ε
4
V2
2
8
εβr2+|p|2
V1
2+64β2r4
ελ
1
2
1
V1
2+ε
4
V2
2,
(85)
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
with the similar calculation of (80), we have
X2
V,
V
M5|z(θtω)|+σ2
2λ
1
4
1|z(θtω)|2
V
2
E1
,(86)
where M5is given by (55). From (84), (85) and (86), we get
d
V2
E1
dt ≤−ε
V2
E1+2
M5|z(θtω)|+σ2
2λ
1
4
1|z(θtω)|2
V
2
E1
+2
8
εβr2+|p|2
V1
2+64β2r4
ελ
1
2
1
V1
2
M6+2M5|z(θtω)|+σ2
λ
1
4
1|z(θtω)|2
V
2
E1
,
(87)
where M6>0is defined by (56). Thus
VE1e0
tM6+2M5|z(θsω)|+σ2
λ
1
4
1
|z(θsω)|2ds
IE1,
by the relation between
V
2
E1
and
U
2
E1
, we have
U
2
E12
1+ σ2
λ
1
4
1|z(θtω)|2
V2
E1
2
1+ σ2
λ
1
4
1|z(θtω)|2
e0
tM6+2M5|z(θsω)|+σ2
λ
1
4
1
|z(θsω)|2ds
IE1.
The random variable z(θtω)is tempered, which together with Lemma 1 gives that
2σ2
λ
1
4
1|z(θtω)|2e0
tM6+2M5|z(θsω)|+σ2
λ
1
4
1
|z(θsω)|2ds
0,t+,
and
2e0
tM6+2M5|z(θsω)|+σ2
λ
1
4
1
|z(θsω)|2ds
<,
then, t
2,tt
2,
U
2
E1p
2, set p2= max{1,p
2},we obtain
U
2
E1
(0,t, ω)=
U
2
E1
(t, 0
tω)p2IE1,(88)
then we can get
E(p2(ω)) <.(89)
In the rest paper, the value of
U
at
t
=1is still denoted by
U
, and
U
(1) =
DSεω, U(1)I
, where
DSεω, U(1)is the linear solution mapping of system (20).
According to (84) and (89), we get that
DSεω, U(1)
L(E1,E1)= sup
IE1
DSεω, U(1)IE1
IE1p2(ω),(90)
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
409
and
p2(ω)1,E(ln (p2(ω))) <.(91)
On the other hand, Γsatisfies
Γ(ω)E1≤∥U(1) U(2)E1+∥−
UE1
(p1+p2)IE1
c(ω)IE1,
where c(ω) = max{1,p
1+p2}. Obviously, when t=1
Γ(1)E1c(ω)IE1.
Combining (91) with (89), we have
E(ln c(ω)) <,E(c(ω)) <.
Since
Sε
(
ω
) :=
Sε
(1
), merging with
(88)
and
(25)
, we can conclude that
Sε
(
ω
)is almost surly uniform
differentiable on A(ω).
CONCLUSIONS
This paper consider global stochastic stability of the Euler-Bernoulli beam equations excited by
multiplicative white noise. The system can induce a RDS which owns global random attractors,
moreover, Hausdorff dimension of the attractor is finite. Specially, when
λ1
2
1
ε
28
ε(βr2+|p|)22M7
πµ 2M8
2µ
16β2r4
ελ
1
4
1
,
the Hausdorff dimension is 0, which indicates that the stochastic Euler-Bernoulli beam possesses a
random fixed point which is global stochastic stability.
ACKNOWLEDGEMENT
This work is supported by Natural Science Foundation of Shandong Province (No. ZR2020MA054) and
National Natural Science Foundation of China (No. 12072178; No. ).
REFERENCES
[1]
Crespo da Silva M, Glynn C. Nonlinear flexural-flexural-torsional dynamics of inextensional beams.
i. equations of motion. Journal of Structural Mechanics 1978; 6(4):437–448.
[2]
Tajik M, Karami Mohammadi A. Nonlinear vibration, stability, and bifurcation analysis of
unbalanced spinning pre-twisted beam. Mathematics and Mechanics of Solids 2019; 24(11):3514–
3536.
[3]
Dai H, Wang Y, Wang L. Nonlinear dynamics of cantilevered microbeams based on modified couple
stress theory. International Journal of Engineering Science 2015; 94:103–112.
[4]
Chen LQ, Yang XD. Transverse nonlinear dynamics of axially accelerating viscoelastic beams
based on 4-term galerkin truncation. Chaos, Solitons & Fractals 2006; 27(3):748–757.
[5]
Pellicano F, Vestroni F. Nonlinear dynamics and bifurcations of an axially moving beam. Journal
of Vibration and Acoustics 2000; 122(1):21–30.
[6]
Mao XY, Ding H, Chen LQ. Internal resonance of a supercritically axially moving beam subjected
to the pulsating speed. Nonlinear Dynamics 2019; 95(1):631–651.
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
410
and
p2(ω)1,E(ln (p2(ω))) <.(91)
On the other hand, Γsatisfies
Γ(ω)E1≤∥U(1) U(2)E1+∥−
UE1
(p1+p2)IE1
c(ω)IE1,
where c(ω) = max{1,p
1+p2}. Obviously, when t=1
Γ(1)E1c(ω)IE1.
Combining (91) with (89), we have
E(ln c(ω)) <,E(c(ω)) <.
Since
Sε
(
ω
) :=
Sε
(1
), merging with
(88)
and
(25)
, we can conclude that
Sε
(
ω
)is almost surly uniform
differentiable on A(ω).
CONCLUSIONS
This paper consider global stochastic stability of the Euler-Bernoulli beam equations excited by
multiplicative white noise. The system can induce a RDS which owns global random attractors,
moreover, Hausdorff dimension of the attractor is nite. Specially, when
λ1
2
1
ε
28
ε(βr2+|p|)22M7
πµ 2M8
2µ
16β2r4
ελ
1
4
1
,
the Hausdorff dimension is 0, which indicates that the stochastic Euler-Bernoulli beam possesses a
random fixed point which is global stochastic stability.
ACKNOWLEDGEMENT
This work is supported by Natural Science Foundation of Shandong Province (No. ZR2020MA054) and
National Natural Science Foundation of China (No. 12072178; No. ).
REFERENCES
[1]
Crespo da Silva M, Glynn C. Nonlinear flexural-flexural-torsional dynamics of inextensional beams.
i. equations of motion. Journal of Structural Mechanics 1978; 6(4):437–448.
[2]
Tajik M, Karami Mohammadi A. Nonlinear vibration, stability, and bifurcation analysis of
unbalanced spinning pre-twisted beam. Mathematics and Mechanics of Solids 2019; 24(11):3514–
3536.
[3]
Dai H, Wang Y, Wang L. Nonlinear dynamics of cantilevered microbeams based on modified couple
stress theory. International Journal of Engineering Science 2015; 94:103–112.
[4]
Chen LQ, Yang XD. Transverse nonlinear dynamics of axially accelerating viscoelastic beams
based on 4-term galerkin truncation. Chaos, Solitons & Fractals 2006; 27(3):748–757.
[5]
Pellicano F, Vestroni F. Nonlinear dynamics and bifurcations of an axially moving beam. Journal
of Vibration and Acoustics 2000; 122(1):21–30.
[6]
Mao XY, Ding H, Chen LQ. Internal resonance of a supercritically axially moving beam subjected
to the pulsating speed. Nonlinear Dynamics 2019; 95(1):631–651.
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
[7]
Wiggins S. Global bifurcations and chaos: analytical methods, vol. 73. Springer Science & Business
Media: Heidenburg, 2013.
[8]
Zhang W, Yao M, Zhang J. Using the extended melnikov method to study the multi-pulse global
bifurcations and chaos of a cantilever beam. Journal of Sound and Vibration 2009; 319(1-2):541–
569.
[9]
ZHANG W, Yao M. Theories of multi-pulse global bifurcations for high-dimensional systems and
application to cantilever beam. International Journal of Modern Physics B 2008; 22(24):4089–4141.
[10]
Zhou L, Ji P, Chen F. Chaos and subharmonic bifurcation of a composite laminated buckled beam
with a lumped mass. Chaos, Solitons & Fractals 2021; 147:110 933.
[11]
Hamed Y, El Shehry A, Sayed M. Nonlinear modified positive position feedback control of
cantilever beam system carrying an intermediate lumped mass. Alexandria Engineering Journal
2020; 59(5):3847–3862.
[12]
Emam SA, Nayfeh AH. On the nonlinear dynamics of a buckled beam subjected to a primary-
resonance excitation. Nonlinear Dynamics 2004; 35(1):1–17.
[13]
Ghayesh MH. Nonlinear forced dynamics of an axially moving viscoelastic beam with an internal
resonance. International Journal of Mechanical Sciences 2011; 53(11):1022–1037.
[14]
Ghayesh MH. Subharmonic dynamics of an axially accelerating beam. Archive of Applied Mechanics
2012; 82(9):1169–1181.
[15]
Zhang W. Global and chaotic dynamics for a parametrically excited thin plate. Journal of Sound
and vibration 2001; 239(5):1013–1036.
[16] Mao X. Stochastic differential equations and applications. Elsevier, 2007.
[17]
Constantin P, Foiaş C, Temam R. Attractors representing turbulent flows, vol. 314. American
Mathematical Soc., 1985.
[18]
Pezeshki C, Dowell EH. Generation and analysis of lyapunov exponents for the buckled beam.
International journal of non-linear mechanics 1989; 24(2):79–97.
[19]
Wedig WV. Dynamic stability of beams under axial forces–lyapunov exponents for general fluctua-
ting loads. Structural Dynamics. Routledge, 2022; 141–148.
[20]
Janevski G, Kozić P, Pavlović R. Moment lyapunov exponents and stochastic stability of a thin-
walled beam subjected to eccentric axial loads. Journal of Theoretical and Applied Mechanics 2012;
50(1):61–83.
[21]
Debussche A. Hausdorff dimension of a random invariant set. Journal de mathématiques pures et
appliquées 1998; 77(10):967–988.
[22]
Crauel H. Random point attractors versus random set attractors. Journal of the London Mathe-
matical Society 2001; 63(2):413–427.
[23]
Crauel H. Global random attractors are uniquely determined by attracting deterministic compact
sets. Annali di Matematica pura ed applicata 1999; 176(1):57–72.
[24]
Chen H, Guirao JLG, Cao D, Jiang J, Fan X. Stochastic euler–bernoulli beam driven by additive
white noise: global random attractors and global dynamics. Nonlinear Analysis 2019; 185:216–246.
[25]
Chen H, Guirao JL, Jiang J, Cao D, Fan X. Kernel sections and global dynamics of nonautonomous
euler–bernoulli beam equations. The European Physical Journal Plus 2020; 135(1):1–36.
[26]
Chueshov I, Lasiecka I. Von Karman Evolution Equations: Well-posedness and Long Time Dynamics.
Springer Science & Business Media, 2010.
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
411
[27]
Temam R. Infinite-dimensional dynamical systems in mechanics and physics. Springer-Verlag: New
York, 1997.
[28]
Fan X. Attractors for a damped stochastic wave equation of sine–gordon type with sublinear
multiplicative noise. Stochastic Analysis and Applications 2006; 24(4):767–793.
[29] Arnold L. Random dynamical systems. Springer-Verlag, 1998.
[30]
Crauel H, Flandoli F. Attractors for random dynamical systems. Probability Theory and Related
Fields 1994; 100(3):365–393.
[31]
Schmalfuss B. Measure attractors and stochastic attractors, institut for dynamische systeme.
Technical Report, Bermen University 1995.
[32]
Crauel H, Imkeller P, Steinkamp M. Bifurcations of one-dimensional stochastic differential equations.
Stochastic dynamics 1999; :27–47.
[33] Zheng S. Nonlinear evolution equations. CRC Press: Florida, 2004.
https://doi.org/10.17993/3ctecno.2023.v12n1e43.386-412
3C Tecnología. Glosas de innovación aplicadas a la pyme. ISSN: 2254-4143
Ed.43 | Iss.12 | N.1 January - March 2023
412