Next Article in Journal
Analytical Model of a Multi-Step Straightening Process for Linear Guideways Considering Neutral Axis Deviation
Next Article in Special Issue
Periodic Orbits of Third Kind in the Zonal J2 + J3 Problem
Previous Article in Journal
Systematic Review of Decision Making Algorithms in Extended Neutrosophic Sets
Previous Article in Special Issue
Evaluating Water Use Efficiency in China’s Western Provinces Based on a Slacks-Based Measure (SBM)-Undesirable Window Model and a Malmquist Productivity Index
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Long Time Behavior and Global Dynamics of Simplified Von Karman Plate Without Rotational Inertia Driven by White Noise

1
Division of Dynamics and Control, School of Astronautics, Harbin Institute of Technology, Harbin 150001, China
2
Division of Dynamics and Control, School of Mathematics and Statistics, Shandong University of Technology, Zibo 255000, China
3
School of Mathematics, Southwest Jiaotong University, Chengdu 610031, China
*
Authors to whom correspondence should be addressed.
Symmetry 2018, 10(8), 315; https://doi.org/10.3390/sym10080315
Submission received: 5 June 2018 / Revised: 20 July 2018 / Accepted: 22 July 2018 / Published: 1 August 2018
(This article belongs to the Special Issue New Trends in Dynamics)

Abstract

:
Without the assumption that the coefficient of weak damping is large enough, the existence of the global random attractors for simplified Von Karman plate without rotational inertia driven by either additive white noise or multiplicative white noise are proved. Instead of the classical splitting method, the techniques to verify the asymptotic compactness rely on stabilization estimation of the system. Furthermore, a clear relationship between in-plane components of the external force that act on the edge of the plate and the expectation of radius of the global random attractors can be obtained from the theoretical results. Based on the relationship between global random attractor and random probability invariant measure, the global dynamics of the plates are analyzed numerically. With increasing the in-plane components of the external force that act on the edge of the plate, global D -bifurcation, secondary global D -bifurcation and complex local dynamical behavior occur in motion of the system. Moreover, increasing the intensity of white noise leads to the dynamical behavior becoming simple. The results on global dynamics reveal that random snap-through which seems to be a complex dynamics intuitively is essentially a simple dynamical behavior.

1. Introduction

1.1. Background and Literatures Review

There exists an essential difference between full Von Karman plate without rotational inertia and simplified Von Karman plate without rotational inertia. From the view of physics, the former takes account into the acceleration in-plane and the latter neglects it [1,2]. In the mathematical standpoint, the governing equations of full Von Karman plates without rotational inertia comprise coupled plate equations and wave equations, while the coupled plate equation and elliptic equation compose the governing equations of Von Karman plate [3].
The definition of global random attractors for random dynamical system (RDS) established by Arnold [4] were proposed by Crauel and Flandoli [5] and Schmalfuss [6]. The former developed the theory of global random attractors in phase space, while the random attractor is seen as a subset in the space of probability measures by Schmalfuss. Afterwards, Crauel et al. [7] introduced a notion of global random attractors which is accessible to the researcher who are not familiar to the probabilistic language. Furthermore, the assertion that global random attractors are uniquely determined by attracting deterministic compact sets in phase space was attained by Crauel in [8]. Invoking these theories, the existence of global random attractors for RDS related to a plenty of mathematical physics problems have been studied by many researchers, (e.g., [9,10,11,12,13,14] and the references therein).
Von Karman plate equation is a well-known model that arise in nonlinear elastodynamics, which can be found in many engineering applications, for instance, wing skin in airplane, vertical fins of High-speed aircraft, etc. For more details, see the Introduction in Monograph [1]. There is an abundant achievements on Von Karman plate and the brief list given below is by no means exhaustive. Invoking the adjoint method, Pappalardo and Guida [15] considered the optimal control problem associated with the vibration of Von Karman plates. Utilizing the theory of plate theory, monohull ship was modelled by Fortuna and Muscato [16]; moreover, the problem of identification and adaptive control were also investigated. To compute the modal parameters of the plates, the system identification algorithm was proposed by Pappalardo and Guida [17]. A survey on computational methods for motion of multibody systems (include plates) was made by Pappalardo and Guida [18]. For the sake of computing the motion of large deformation of the plates, Pappalardo et al. [19,20,21] developed different kinds of plate/shell finite elements.
From the mathematical view, to address the long time behavior of the mathematical physics problems, it must be verified that they can generate a dynamical system, which can be accomplished by achieving the existence and uniqueness of the solution for the systems. Lasiecka [22] was concerned with weak, classical and intermediate solutions to full von Karman equations. With respect to the nonautonomous case, one can refer to Leiva and Sivoli [23] and Abels et al. [24]. According to the proof of a “sharp regularity” estimates of the Von Karman bracket, the consequence of global existence, uniqueness and regularity of solutions for simplified Von Karman plate with nonlinear boundary dissipation can be founded in Favini et al. [25]. For more detail, one can refer to the Monograph [26]. As for the long time behavior of the Von Karman plate equations, the global attractors as well as inertial manifolds for the system in autonomous situation were studied by Chueshov and Lasiecka [27] and Chueshov and Lasiecka [28], respectively. Lasiecka [29] studied the uniform decay rates for thermoelastic full von karman system. Dynamics of a thermoelastic von Karman plate in a subsonic gas flow was addressed by Ryzhkova [30]. For a von Karman plate equation with a boundary memory condition which can even be a fractional damping, Park and Sun [31] tackled the uniform decay of the solution. For more detail of the research status of long time behavior of the Von Karman plate before the year 2010, one can refer to Monograph [26]. The study of long-time dynamics of a von Karman equation with time delay is due to Park [32]. Chueshov [33] investigated questions related to global attractors for delayed, nonrotational von Karman plates without any damping in the status of flow–structure interactions. Eliminating flutter for clamped von Karman plates in subsonic flows was concerned by Lasiecka and Webster [34]. Without assuming large values for the coefficient of damping, Khanmamedov [35] proved the existence of the global random attractors for Von Karman plate equation. With a very strict assumption on the coefficient of the weakly damping, the existence of global random attractors for simplified Von Karman plate without rotational inertia driven by multiplicative white noise was studied by Chen et al. [36].
The study on dynamics of Von Karman plates can be divided into two parts: the investigation of local dynamics and investigation of global dynamics. There exists abundant studies on local dynamics for Von Karman plates. For instance, invoking the Bubnov–Galerkin approach, Awrejcewicz and Krysko [37] analyzed the complex parametric vibrations of plates and shells. Nonlinear vibration and dynamic response as well as Thermal post-buckling of functionally graded thermoelastic Von Karman plates was considered by Huang and Shen [38] and Park and Kim [39], respectively. Employing the Homotopy perturbation technique, Rashidi et al. [40] studied the nonlinear vibration of Von Karman rectangular plate. Ghayesh et al. [41] tackled the nonlinear dynamics of axially moving Von Karmam plates. Ghayesh and Farokhi [42] were devoted to handling the nonlinear dynamics of Von Karman plate in MEMS. The nonlinear vibrations of viscoelastic Von Karman plate was analyzed by Amabili [43]. The global dynamics of nonlinear systems which can reveal more dynamic information than local dynamics are important in engineering applications. Compared with the literature on local dynamics, the investigation related to global dynamics is insufficient. Global dynamics of four-dimensional perturbed Hamiltonian systems and parametrically forced mechanical systems were addressed by Wiggins [44] and Feng and Wiggins [45], respectively. The technique employed in those works is Melnikov method which was invoked by Zhang to tackle a parametrically Von Karman plate in [46]. Due to the lack of analytical tools, numerical method is the main approach to study the global dynamics of nonlinear systems. According to Cell to Cell mapping method proposed by Hsu [47], Xu et al. [48] addressed global stochastic bifurcation in Duffing system.
There exist two standpoints in study on dynamics of random dynamical system, which are equivalent in the deterministic case, the “static” standpoint and the “dynamical” standpoint. However, two views are very different in the stochastic status (see [4,49,50]). The investigation on dynamics of the RDS associated with vibration of Von Karman plates in this paper means study the of “dynamical” dynamics of the systems; alternatively, the global dynamics in this paper are understood as the change in the pattern of existing probability invariant measures of the RDS. There exist some results on global dynamics on RDS. Crauel and Flandoli [49] asserted that additive noise destroys pitchfork bifurcation in one dimensional system. The statement that parametric noise (even a multiplicative white noise) destroys Hopf bifurcation was duo to Arnold et al. [51]. Wang [52] focused on the bifurcation for stochastic parabolic equations. The investigation on stochastic bifurcation in Duffing system by the theory of random attractors was due to Schenk-Hoppé [53]. According to some invariant manifolds to derive the lower bounds on the dimension of global random attractors, Caraballo et al. [54] studied the stochastic pitchfork bifurcation of the reaction diffusion equation with multiplicative white noise.

1.2. Formulation and Contribution of This Investigation

In some circumstances, the Von Karman plate equation that epitomizes certain distinct features and mathematical difficulties which lead to the “splitting method” [55], a traditional approach in the study on the existence of global random attractors, for extensive mathematical problems becomes invalid, such as SAVKP and SMVKP (introduced in Section 2.1). The existence of global attractors for the system in deterministic case (such as Chueshov and Lasiecka [27]) and stochastic case (e.g., [36]) relies on large enough value of damping coefficient. To our best knowledge, there hardly exists results on global random attractors for SAVKP and SMVKP with arbitrary small coefficient of the weakly damping.
Recently, the study on dynamics of Von Karman plate equation mainly focused on the local dynamics, inspired by Crauel [49] and Schenk-Hoppé [53]. Based on the existence of global random attractors and the relationship between invariant measure and global random attractor summarized in Proposition 2 in Section 3.1, the dynamics of Von Karman plate can be accomplished by employing the stochastic subdivision algorithm method proposed by Keller and Ochs [56] to achieve the global random attractors numerically.
As far as we know, the consequence of investigation in this aspect also do not be published in any composition. The purpose of this paper is to investigate the existence of global random attractors for SAVKP and SMVKP and to derive the global dynamics by achieving the structure of their global random attractors.

1.3. Organization of the Paper

The rest of this paper is organized as follows. In Section 2, the mathematical description of model and main results main results are given. Section 3 is intended to provided preliminary results employed in accomplishing the main proof which are given in Section 4. Finally, based on the main results listed in Section 3, summary and conclusions is made in Section 5.
Finally, to express the results and their respective proofs succinctly, the following conventions are made. Unless otherwise stated, in the sequel, the letter C , C ˜ , C ^ , C i , C ˜ i , C ^ i ( i = 1 , 2 , 3 , , ) are positive constants; in addition, C ( a , , b ) , C ˜ ( a , , b ) , C ^ ( a , , b ) and C i ( a , , b ) , C ˜ i ( a , , b ) , C ^ i ( a , , b ) ( i = 1 , 2 , 3 , , ) are positive constants depended on a , , b .

2. The Mathematical Model and Main Results

Section 2.1 is used to make the mathematical description of the model considered in this paper. The main results of this paper are listed in Section 2.2.

2.1. Mathematical Description of the Model

Let D R 2 be a bounded domain with boundary Γ ; without loss of generality, assume the origin 0 belongs to Γ . Suppose ( x , y ) Γ is an arbitrarily given point, while Γ ( x , y ) denotes the arc, oriented in the usual manner, joining the origin 0 to the point ( x , y ) along the boundary of D. For more details, one can refer to Ciarlet [3]. The governing equation of simplified Von Karman plate without rotational inertia is:
ρ h U t t + c h G ( U t ) + D ¯ Δ 2 U = P + V ( U ) , U ,
with the clamped boundary
x D : U = U n = 0 ,
where U is transversal displacement of the plate. [ · , · ] is Von Karman bracket [26] (also known as Monge-Ampère form [3]) with the form of
[ v , u ] = 2 u x 1 2 2 v x 2 2 + 2 v x 1 2 2 u x 2 2 2 2 u x 1 x 2 2 v x 1 x 2 .
V ( U ) is the Airy function satisfies
Δ 2 V = E h 2 [ U , U ] , x Γ : V = ϕ 0 , V n = ϕ 1 ,
in which the physical parameter E , ρ , h , I , D can refer to [1]. ϕ 0 , ϕ 2 is defined as
ϕ 0 ( x , y ) = x Γ ( x , y ) f 2 d s + y Γ ( x , y ) f 1 d s + Γ ( x , y ) x ¯ f 2 y ¯ f 1 d s ,
and
ϕ 1 ( x , y ) = n 1 ( x , y ) Γ ( x , y ) f 2 d s + n 2 ( x , y ) Γ ( x , y ) f 1 d s ,
where f 1 , f 2 are components of the in-plane force on boundary along the direction x , y , which comply with
n 1 N 11 + n 2 N 12 = f 1 , n 1 N 12 + n 2 N 22 = f 2 ,
where N 11 , N 12 , N 22 are membrane forces in the plate; for more details, one can refer to [3]. Moreover, let V = V 0 + V 1 , V 0 be the solution of the following system
Δ 2 V 0 = E h 2 [ w , w ] ,
V 0 = V 0 n = 0 .
V 1 satisfies
Δ 2 V 1 = 0 ,
V 1 = ϕ 0 , V 1 n = ϕ 1 .
Thus, Equation (1) can be rewritten as:
ρ h U t t ρ I Δ U t t + c h G ( U t ) + D ¯ Δ 2 U = F + P + V 1 + V 0 , U .
In some cases, only V 0 is named Airy function, while V 1 is called the in-plane force, (see Chueshov and Lasiecka [26]). This convention is employed in this paper.
To formulate the system tackled in this paper, some spaces are introduced in the following. Let u u L 2 ( D ) , u s u H 0 s ( D ) , ( u , v ) ( u , v ) L 2 ( D ) , ( u , v ) s = ( u , v ) H 0 s ( D ) , where H s ( D ) , H 0 s ( D ) , s R are the usual Sobolev Spaces. Let A = Δ 2 : H 0 2 H 4 ( D ) L 2 ( D ) , then A is self-adjoint, positive, unbounded linear operators and A 1 L ( L 2 ( D ) ) is compact. Therefore, their eigenvalues { Λ i } i N satisfy 0 < Λ 1 Λ 2 and the corresponding eigenvalues { w i } i = 1 forms an orthonormal basis in L 2 ( D ) . Then, we can interpret the power of A s , s R by the method developed by Temam [55]. Specifically, A 1 4 u = u , A 1 2 u = Δ u , u D ( A ) , however, it is mentioned here that A 1 2 Δ , A 1 4 . Nevertheless, it is emphasized here that A 1 2 Δ with the boundary in Equation (2). In fact, the operator A 1 : = Δ with the boundary condition in Equation (2) is a self-adjoint, positive, unbounded linear operators from L 2 ( D ) to L 2 ( D ) and A 1 L ( L 2 ( D ) , L 2 ( D ) ) is compact. Thus, the power of ( Δ ) s , s R can be defined; furthermore, A 1 s u A s 2 u , s [ 0 , 1 ] , u D ( A ) .
Suppose P is a stochastic pressure signified by white noise, then the dynamics equation of abstract dimensionless clamped simplified Von Karman plate without rotational inertia driven by white noise are as follow
u t t + α u t + A u k [ F ( u ) + F 0 , u ] = σ q ( x ) W ˙ ,
u t t + α u t + A u k [ F ( u ) + F 0 , u ] = σ u W ˙ ,
with the clamped boundary condition
x D : u = u n = 0 ,
where u : R × D R is the dimensionless transversal displacement of the plate. σ 0 , k > 0 are a given constant,
q ( x ) ( H 0 2 H 4 ) ( D ) .
W is the one dimensional two-sided real-valued standard Wiener process, and W ˙ is called white noise. α > 0 is the coefficient of damping.
Equation (10a) describes abstract dimensionless clamped simplified Von Karman plate without rotational inertia driven by additive/multiplicative white noise, respectively. Furthermore, F satisfies
Δ 2 F = [ u , u ] ,
x Γ : F = F n = 0
and F 0 is in agreement with
Δ 2 F 0 = 0 ,
x Γ : F 0 = g 1 , F 0 n = g 2
where g 1 H 3 2 ( Γ ) , g 2 H 1 2 ( Γ ) are derived from Equations (3) and (4). By the monograph of Lions and Magenes [57], F 0 H 2 ( D ) , define
N ( u ) = k [ F ( u ) + F 0 , u ] .
The system described by Equations (10a), (10c), and (10e)–(10h) is denoted by SAVKP. The system interpreted by Equations (10b), (10c), and (10e)–(10h) is represented by SMVKP.
Invoking the compactness of A, the Hilbert space D ( A s ) , s R with norm · s and ( · , · ) s can be defined as the mechanism in Temam [55], especially, D ( A 1 2 ) = H 0 2 ( D ) . Moreover, for all s 1 , s 2 R , s 1 s 2 , D ( A s 1 ) can be compact imbedding in D ( A s 2 ) and the following holds
u s 1 Λ 1 s 1 s 2 2 u s 2 , u D ( A s 1 ) .
Let E 1 + 2 s = D ( A 1 2 + s ) × D ( A s ) , 0 s 1 2 equipped with Graph norm and the induced inner product, then they are all Hilbert spaces.
Let ( X , · X ) be a separable space with Borel σ algebra B ( X ) and ( Ω , F , P ) be a probability space. { θ t : Ω Ω , t R } is a family of measure preserving transformations such that ( t , ω ) θ t ω is measurable, θ 0 = id , θ t + s = θ t θ s for all s , t R . Then, the flow θ t together with the probability space ( Ω , F ¯ , P , { θ t } t R ) is called a metric dynamical system. For the particular applications in this paper, the metric dynamical systems generated by a one dimensional two-sided standard Wiener process defined on a Probability space ( Ω , F ¯ , P ) is introduced there. Let Ω = ω | ω ( · ) C ( R , R ) , ω ( 0 ) = 0 , F ¯ is the σ algebra induced by the compact open topology for this set and P is the Wiener measure on F . Set
θ t ω ( · ) = ω ( · + t ) ω ( · ) , t R ,
according to Arnold [4], we have P is ergodic with respect to the flow { θ t } t R . Thus, ( Ω , F ¯ , P , { θ } t R ) is the metric dynamical systems employed in this paper. Moreover, the Ornstein–Uhlenbeck process, which should be used in transforming a stochastic system to a random system, is introduced as follows
d z + μ z d t = d W , z ( ) = 0 .
in which μ > 0 . The general form for the solution of Equation (13) is
z ( t ) = t e μ ( t s ) d W .
Let
z ( θ t ω ) : = μ 0 e μ τ ( θ t ω ) ( τ ) d τ .
where { θ t } t R is defined by Equation (12) in Section 2. Merging with integration by parts, z ( θ t ω ) is the solution for the system in Equation (13).
Although by no means always, it will be convenient to reduce Equations (10a) and (10b) to an evolution equation of the first order in time in the following manner. Let v ̲ = u t + ε u , U = [ U 1 , U 2 ] T = [ u , v ̲ ] T , then Equation (10a) can be transformed to the ensuing form
d U d t = A 1 U + F 1 ( U )
where
A 1 = ε I , I A + ε ( α ε ) , ( α ε ) I , F 1 ( U ) = F 11 ( U ) F 12 ( U )
and
F 11 ( U ) = 0 , F 12 = N ( U 1 ) + σ q ( x ) W ˙ .
The system described by Equations (10c), (15) and (10e)–(10h) denoted by SAVKPT1. Obviously, SAVKP is equivalent to SAVKPT1.
To accomplish the stabilization estimation for the solution of SAVKPT1, the following systems is needed. Suppose U ( 1 ) = [ U 1 ( 1 ) , U 2 ( 1 ) ] = [ u 1 , v ̲ 1 ] , U ( 2 ) = [ U 1 ( 2 ) , U 2 ( 2 ) ] = [ u 2 , v ̲ 2 ] are two solution of SAVKPT1, then
d U ( 1 ) U ( 2 ) d t = A 1 U ( 1 ) U ( 2 ) + F ˜ 1 ( U ( 1 ) U ( 2 ) ) ,
where F ˜ 1 ( U ( 1 ) U ( 2 ) ) = F 1 ( U ( 1 ) ) F 1 ( U ( 2 ) ) .
Furthermore, let v = v ̲ σ q ( x ) z ( θ t ω ) , Φ = [ u , v ] T and
[ u , v ] T = [ u , v ̲ σ q ( x ) z ( θ t ω ) ] T T A ( θ t ω ) [ u , v ̲ ] T , [ u , v ̲ ] E 1 ,
thus
d Φ d t = A 1 Φ d t + F 2 ( z ( θ t ω ) ) + F 3 ( Φ )
where Φ = [ u , v ] T ,
F 2 ( z ( θ t ω ) ) = F 21 ( z ( θ t ω ) ) F 22 ( z ( θ t ω ) ) = σ q ( x ) z ( θ t ω ) σ ( μ α + ε ) q ( x ) z ( θ t ω ) , F 3 ( Φ ) = F 31 ( Φ ) F 32 ( Φ ) ,
here
F 31 ( Φ ) = 0 , F 32 ( Φ ) = N ( U 1 ) .
Equation (18) is a partial differential equations with random coefficient which can be studied ω by ω . Let SAVKPT2 signify the system described by Equations (18), (10c) and (10e)–(10h). It is emphasized that SAVKPT2 is not equivalent to SAVKPT1.
Analogously, let v ̲ = u t + ε u , U = [ U 1 , U 2 ] T = [ u , v ̲ ] T and
[ u , v ] T = [ u , v ̲ σ u z ( θ t ω ) ] T T M ( θ t ω ) [ u , v ̲ ] T , [ u , v ] E 1 , ω Ω ,
define Φ = [ Φ 1 , Φ 2 ] T = [ u , v ] T , the following system associated with SMVKP can be attained.
d U d t = A 1 U + F 1 ( U )
in which
A 1 = ε I , I A + ε ( α ε ) , ( α ε ) I , F 1 = F 11 ( U ) F 12 ( U )
and
F 11 ( U ) = 0 , F 12 ( U ) = N ( U 1 ) + σ U 1 W ˙ .
d Φ d t = A 1 Φ + F 2 ( Φ , θ t ω ) + F 3 ( Φ , θ t ω )
where
F 2 ( Φ , θ t ω ) = F 21 ( Φ , θ t ω ) F 22 ( Φ , θ t ω ) = σ z ( θ t ω ) Φ 1 σ z ( θ t ω ) Φ 2 , F 3 ( Φ , θ t ω ) = F 31 ( Φ , θ t ω ) F 32 ( Φ , θ t ω )
F 31 ( Φ , θ t ω ) = 0 , F 32 ( Φ , θ t ω ) = σ ( μ α + 2 ε ) z ( θ t ω ) Φ 1 σ 2 z 2 ( θ t ω ) Φ 1 + N ( Φ 1 ) ,
SMVKPT1 represents the systems defined by Equations (20), (10c) and (10e)–(10f), then SMVKP and SMVKPT1 are equivalent. The system described by Equations (21), (10c) and (10e)–(10f) is denoted by SMVKPT2.
Furthermore, assume Φ ( 1 ) = [ Φ 1 ( 1 ) , Φ 2 ( 1 ) ] T = [ u 1 , v 1 ] T , Φ ( 2 ) = [ Φ 1 ( 2 ) , Φ 2 ( 2 ) ] T = [ u 2 , v 2 ] T are two solutions of SAVKPT2, thus
d Φ ( 1 ) Φ ( 2 ) d t = A 1 Φ ( 1 ) Φ ( 2 ) + DF 1 ( Φ ( 1 ) Φ ( 2 ) , θ t ω ) + DF 2 ( Φ ( 1 ) Φ ( 2 ) , θ t ω ) ,
in which
DF 1 ( Φ ( 1 ) Φ ( 2 ) , θ t ω ) = D F 11 ( Φ ( 1 ) Φ ( 2 ) , θ t ω ) D F 12 ( Φ ( 1 ) Φ ( 2 ) , θ t ω ) = F 2 ( Φ , θ t ω ) F 2 ( Φ , θ t ω ) , DF 2 ( Φ ( 1 ) Φ ( 2 ) , θ t ω ) = D F 21 ( Φ ( 1 ) Φ ( 2 ) , θ t ω ) D F 22 ( Φ ( 1 ) Φ ( 2 ) , θ t ω ) = F 3 ( Φ , θ t ω ) F 3 ( Φ , θ t ω ) .
Equation (22) is used to obtain the stabilization estimation for SMVKPT1.
Remark 1.
For the sake of brevity, when no ambiguity is possible, the symbols used in SAVKPT and SMVKPT, SAVKPT1 and SMVKPT1, and SAVKPT2 and SMVKPT2 are the same or similar. Since each of the symbols has a clearl explanation, it is not confusing to express the main results in this paper. However, it must be kept in mind that they are not the same.

2.2. Main Results

Approved by the equivalent between SAVKPT and SAVKPT1, and SMVKPT and SMVKPT1, it is enough to only address the dynamical behavior of SAVKPT1 and SMVKPT1. This subsection is used to present the main results of this paper. Let
ε 0 = Λ 1 α 2 Λ 1 + α 2 , 0 < ε ε 0 ,

2.2.1. Random Attractors in Additive White Noise Case

This part is devoted to providing the main results for SAVKPT1. The following Theorem considers the existence and uniqueness of solution for SAVKPT2.
Theorem 1.
For any given initial value t = τ : Φ = Φ τ E 1 , there exists a uniqueness (mild) solution for SAVKPT2 Φ ( t , τ , ω ) C ( [ τ , τ + T ] , E 1 ) , t [ τ , τ + T ] , T > 0 .
Furthermore, let Φ A ( t , τ , ω ; Φ τ ) = φ A ( t , τ , ω ) Φ τ , which means that φ A ( t , τ , ω ) is the solution mapping of SAVKPT2. Then,
φ A ( t τ , 0 , θ τ ω ) = φ A ( t , τ , ω ) , φ A ( τ , τ , ω ) = i d .
Let
S ε A ( t , ω ) = φ A ( t , 0 , ω ) = φ A ( 0 , t , θ t ω ) , t 0 , ω Ω ,
by means of Theorem 1, S ε A ( t , ω ) is the RDS induced by SAVKPT2. Correspondingly, SAVKPT1 can also generate a RDS S A ( t , ω ) which is defined as
S A ( t , ω ) = T A ( θ t ω ) 1 S ε A ( t , ω ) T A ( ω ) : E 1 E 1 .
here T A ( θ t ω ) is defined by Equation (17). Furthermore, the solution mapping determined by SAVKPT1 is denoted by ϕ A ( t , τ , ω ) , then
S A ( t , ω ) = ϕ A ( t , 0 , ω ) = ϕ A ( 0 , t , θ t ω ) , t 0 , ω Ω .
The following turns to the existence of global random attractors for SAVKPT1.
Let δ 0 be any given positive constant, k ¯ = max { 1 , k } and
μ > μ 0 = 64 | σ | 2 q 2 π ε 2 ,
K = K 0 + | z ( θ t ω ) K 1 + | z ( θ t ω ) | 2 K 2 ,
ρ ( ω ) = σ 2 q ( x ) 2 | z ( ω ) | 2 + 2 0 e ϵ t + t 0 2 k ¯ | σ z ( θ s ω ) | q d k K d s + δ 0 ,
r ( ω ) = ρ ( ω ) .
in which
K 0 = 3 k ε 2 C ( 7 6 , C 1 , C 2 ) F 0 2 2 , K 1 = 2 k ¯ k | σ | q C ( 1 4 , C 1 , C 2 ) F 0 2 2 , K 2 = 2 σ 2 ( μ α + ε ) 2 k 0 q ( x ) 2 + 4 k 2 | σ | 2 q 1 2 F 0 2 2 ε + 4 σ 2 ε A 1 2 q ( x ) 2 ,
C ( η , C 1 , C 2 ) , η = 1 4 , 7 6 are constants that satisfy Equation (45) in Section 3.2.1.
The next theorem considers the existence and expectation of radius of global random attractors for S A ( t , ω ) .
Theorem 2.
S A ( t , ω ) possesses the global random attractors A ( ω ) , ω Ω in E 1 which satisfies P -a.s.
A ( ω ) B ( 0 , r ( ω ) ) , E ( ρ ( ω ) ) M < .
where B ( 0 , r ( ω ) ) denotes the open ball centered at the origin with radius is r ( ω ) , while M is given by Equation (A24) in Section 3.2.1.

2.2.2. Random Attractors in Multiplicative White Noise Case

The main results for SMVKPT1 is given in this part.
The form of following theorem is very similar to the Theorem 1.
Theorem 3.
For any given initial value t = τ : Φ = Φ τ E 1 , SMVKPT2 possesses uniqueness (mild) solution Φ M ( t , τ , ω ; Φ τ ) C ( [ τ , τ + T ] , E 1 ) , t [ τ , τ + T ] , T > 0 .
Furthermore, set Φ M ( t , τ , ω ; Φ τ ) = φ M ( t , τ , ω ) Φ τ , which indicates that φ B ( t , τ , ω ) is the solution mapping of SMVKPT2. Thus,
φ M ( t τ , 0 , θ τ ω ) = φ M ( t , τ , ω ) , φ M ( τ , τ , ω ) = i d .
Therefore, set
S ε M ( t , ω ) = φ M ( t , 0 , ω ) = φ M ( 0 , t , θ t ω ) , t 0 , ω Ω ,
then, S ε M ( t , ω ) the RDS induced by SMVKPT2 and RDS S M ( t , ω ) defined as follow
S M ( t , ω ) = T M ( θ t ω ) 1 S ε M ( t , ω ) T M ( ω ) : E 1 E 1 .
is the RDS generate by SMVKPT1, where T M ( θ t ω ) is defined by Equation (19). In additive, let
S M ( t , ω ) = ϕ M ( t , 0 , ω ) = ϕ M ( 0 , t , θ t ω ) , t 0 , ω Ω .
Thus, ϕ ( t , τ , ω ) is the solution mapping of SMVKPT1.
The following turns to the existence of global random attractors for SMVKPT1.
Suppose δ 0 is any given positive constant, μ = α 2 ε and
3 σ 2 μ Λ 1 + 12 σ π μ < ε .
Set
C 0 ( θ t ω ) = 3 2 | σ z ( θ t ω ) | + σ 2 | z ( θ t ω ) | 2 Λ 1 , ρ 0 ( ω ) = t e ε t + t 0 C 0 ( θ s ω ) d s K d s + δ 0 ,
ρ ( ω ) = ρ 0 ( ω ) + k C ( 1 2 , C 1 , C 2 ) F 0 2 2 , K 0 = ε k C ( 3 2 , C 1 , C 2 ) F 0 2 2 , K 1 = 8 k C ( 1 8 , C 1 , C 2 ) F 0 2 2 , K 2 = 4 k σ 2 Λ 1 C ( 1 8 , C 1 , C 2 ) F 0 2 2 ,
K = K 0 + | z ( θ t ω ) | K 1 + | z ( θ t ω ) | 2 K 2 ,
r ( ω ) = | σ z ( ω ) | 2 Λ 1 + 1 ρ ( ω ) ,
where C ( η , C 1 , C 2 ) , η = 1 2 , 3 2 , 1 8 are constants satisfying Equation (45) in Section 3.2.2. Then, the existence and finite expectation of radius of global random attractors for S M ( t , ω ) can be asserted by the next theorem.
Theorem 4.
There exists global random attractors A ( ω ) , ω Ω for S M ( t , ω ) in E 1 ; moreover,
A ( ω ) B ( 0 , r ( ω ) ) , P - a . s . E ( ρ ( ω ) ) M < ,
where B ( 0 , r ( ω ) ) denotes the open ball centered at the origin with radius is r ( ω ) , while M is given by Equation (A43) in Section 3.2.2.
Comparing Theorems 1 and 3, as well as Theorems 2 and 4, their forms are the same or similar, while there exists essential difference between them, which is expounded in Section 5.

2.2.3. Global Dynamics in Both Additive and Multiplicative White Noise Cases

Based on the theoretical results and the relationship between invariant measures and global attractors introduced in Proposition 2 in Section 3.1, the rest of this subsection is dedicated to studying the global dynamics of the stochastic Von Kaman plates, which is accomplished by deriving the components of global random attractor numerically, the main components are referred as global random point attractor and global random basic attractor. The modal equations associated with the stochastic Von Karman plates which are not display here (see Equation (A1) in Appendix A) can be obtained by employing inertial manifold with delay [58] and nonlinear gakerlin method [59].
Let s 1 = 2 , s 2 = 2 , s 0 = 4 , D = [ 0 , 8 ] × [ 0 , 6 ] R 2 in Appendix A, Figure 1 shows the model for vibration of Von Karman plate. The eigenvalues { Λ i } i = 1 4 and eigenvectors { w i } i = 1 4 of operator A and integration with respect to space variable in Equation (A2) listed in Appendix A can be performed by COMSOL with Matlab [60], and then the solution of model equations can be obtained by stochastic Runge–Kutta method [61].
The situation of additive white noise. In this case, the P in Figure 1 is equal to σ q ( x ) W ˙ , let α = 0.5 , k = 2 , q ( x ) = Λ 1 ( 0 . 2 w 1 + 0 . 2 w 2 + 0 . 2 w 3 + 0 . 2 w 4 ) and ( f 1 , f 2 ) are component of the in-plane force on boundary along the direction x , y , thus the g 1 , g 2 in Equation (10h) can be derived by Equations (3) and (4). Since D is rectangle, g 2 = 0 . Furthermore, let f 1 = f 2 , the dynamics of simplified Von Karman plate without rotational inertia driven by additive white noise is signified by the motion of the position x = 3 , y = 2 of the plate, which are studied in the following cases.
Case I. Let f 1 = 1.2 , σ = 0.2 ; the global random basic attractor, global random point attractor and global random attractors for SAVKP are the same. Figure 2 shows that global random basic attractor is a random fixed point which supports a invariant Markov measure δ x 0 ( ω ) . Furthermore, δ x 0 ( ω ) is almost surely global stability.
Case II. Set f 1 = 1.3 , σ = 0.2 , in this situation, global random basic attractor is equivalent to global random point attractor. Global random basic attractor (see Figure 3c) and its section (see Figure 3a) indicate that the system possesses two invariant Markov measures δ x 1 ( ω ) , δ x 2 ( ω ) which are supported by two fixed random points. Figure 3b illustrates the section of the global random attractor of SAVKP.
Case III. let f 1 = 3 , σ = 0.2 , Figure 4 describes the invariant measures and random attractors for SAVKP. Global random basic attractor is equivalent to global random point attractor in this status. Section of global random basic attractor demonstrated by Figure 4a reveals that the steady states of the system comprise four parts, which means that are least four stable invariant Markov measures for SAVKP exist. In addition, in Figure 4a, it can be obtained that the local dynamics of the system may be complex. The sketch of global random basic attractor is shown by Figure 4b,
Case IV. when f 1 = 1.3 , σ = 2 , global random basic attractor is also equivalent to global random point attractor in this case. The numerical results on the invariant measures and random attractors for SAVKP (see Figure 5) expose the system has a almost surely global stability invariant Markov measure supported by a random fixed points.
The situation of multiplicative white noise. In this case, the P in Figure 1 is equal to σ u W ˙ , while the remaining parameters are chosen to be the same as in the situation of additive white noise. The the dynamics of SMVKP are studied in the following cases.
Case I. Let f 1 = 1.1 , σ = 0.02 Λ 1 , similar to the Case I in additive noise, global random basic attractor, global random point attractor and global random attractors for SAVKP are the same in this circumstance. The assertion that there exists an almost surely global stability invariant measure supported by a random fixed points for SMVKP can be attained by Figure 6.
Case II. When f 1 = 1.2 , σ = 0.04 Λ 1 , the numerical results on invariant measures and random attractors were described by Figure 7. Invoking the section of global random basic attractor (see Figure 7a), it can be obtained that the system possesses two local stable invariant Markov measures, which together with the numerical results of section global random point attractor described by Figure 7b give that another invariant measure exists, which could even be a unstable invariant Markov measure.
Case III. Set f 1 = 2.9 , σ = 0.04 Λ 1 , Figure 8 expresses the numerical results on global random attractors for SMVKP. The section of global random basic attractor shown by Figure 8a indicates that there exist four local stable invariant Markov measures. Furthermore, SMVKP has another three invariant measures which are interpreted by section of global random point attractor (see Figure 8b). The sketch of global random point attractor is illustrated by Figure 8c.
Case IV. Let f 1 = 1.2 , σ = 0.4 Λ 1 , the similar results in Case IV can be got, the Figure to describe the invariant measure and random attractor is not displayed here.
Some affirmations can be approved by the aforementioned numerical results. For the clamped irrotational inertia Von Karman driven by additive white noise, fixed σ = 0.2 , let f 1 vary from 1.2 to 1.3 , the global D -bifurcation occurs in the motion of the system. Change the value of f 1 to be a big one, such f 1 = 3 , the dynamical behavior becomes much more interesting. From the view of global dynamics, there exits secondary D -bifurcation. The local dynamics of the system is complex. On the other hand, let f 1 = 1.3 and change the σ from 0.2 to 2, the phenomenon of global D -bifurcation disappears. As for clamped irrotational inertia Von Karman driven by multiplicative white noise, fixed σ = 0.04 Λ 1 , the similar global dynamics of the system can be obtained with varying the f 1 from 1.1 to 2.9 . In addition, once the coefficient of the multiplicative white noise becomes big, global D -bifurcation vanishes. Nevertheless, there exist differences between the two cases above. The multiplicative white noise is more likely to result in the appearance of global D -bifurcation and secondary global D -bifurcation in the motion of clamped Von Karman without rotational inertia than additive noises. when the secondary D -bifurcation occurs, the local dynamics of the system driven by additive white noise is more complex than the situation of multiplicative white noise.

3. Preliminary Results

This section pays attention to give preliminaries and derive certain estimates for solution of SAVKPT1 and SMVKPT1, which are very important to prove main results provided in Section 2.2.

3.1. Basic Theory Related to Global Random Attractors

This subsection is devoted to introduce basic theory related to the theory of random attractors used in this paper.
The next Proposition pertinent to z ( θ t ω ) can be found in References [13,62].
Proposition 1.
z ( θ t ω ) defined by Equation (14) satisfies
E [ z ( θ t ω ) ] = 0 , D [ z ( θ t ω ) ] = 1 2 μ ,
and z ( θ t ω ) is a stationary Gauss Process and Markov Process, its probability-distribution function induce a Markov semigroup. Furthermore,
lim t ± 1 t 0 t | z ( θ τ ω ) | 2 d τ = E [ | z ( ω ) | 2 ] = 1 2 μ ,
lim t ± 1 t 0 t | z ( θ τ ω ) | d τ = E [ | z ( ω ) | ] = 1 π μ ,
and
lim t ± z ( θ t ω ) t = 0 .
The sets { z ( ω ) } is tempered with respect to { θ t } t R . Moreover, if μ 2 β , β > 0 , then
E e β s s + t | z ( θ τ ω ) | 2 d τ e β t μ , s R , t 0 .
When μ 3 r 2 , r 0 , the following holds
E e r s s + t | z ( θ τ ω ) | d τ e r t μ , s R , t 0 .
Moreover,
E | z ( θ s ω ) | r = Γ ( r + 1 2 ) π μ r , r R + ,
where Γ ( · ) is usual Γ function.
To give the notion of global random attractors, the definition of RDS inaugurated by Arnold [4] is needed to be given firstly.
Definition 1.
([4]) A RDS on Polish space ( X , d ) with Borel σ-algebra B ( X ) over a metric dynamical system ( Ω , F , P , { θ t } t R ) is a ( B ( R + ) × F × B ( X ) , B ( X ) ) measurable mapping
ϕ : R + × Ω × X X , ( t , ω , x ) Φ ( t , ω , x )
such that, for P a . s . , ω Ω ,
(i).
ϕ ( 0 , ω ) = id on X.
(ii).
ϕ ( t + s , ω ) = ϕ ( t , θ s ω ) ϕ ( s , ω ) for all s , t R + .
A RDS is continuous or differential if ϕ ( t , ω ) : X X is continuous or differential. Furthermore, ϕ ( t , θ t ω ) can be understood as the solution start from t to 0.
The coming definitions related to random attractors for RDS was established by Crauel et al. [5,63].
Definition 2.
A random set K ( ω ) is said to absorb the set B X for a RDS ϕ , if P - a . s . there exists t B ( ω ) such that
ϕ ( t , θ t ω ) B K ( ω ) for all t t B ( ω )
Definition 3.
Let B 2 X be a collection of subsets of X, then a closed random set A ( ω ) is called B -random attractor associated with the RDS ϕ, if P - a . s .
(i).
A ( ω ) is a random compact set.
(ii).
A ( ω ) is invariant i.e., ϕ ( t , ω ) A ( ω ) = A ( θ t ω ) for all t 0
(iii).
For every B B ,
lim t d i s t ( ϕ ( t , θ t ω ) B , A ( ω ) ) = 0
where dist denotes the Hausdorff semidistance:
d i s ( A , B ) = sup x A inf y B d ( x , y ) , A , B X .
When B is composed of all bounded set of X, then A ( ω ) is the global random attractors for ϕ. If B = { { x } : x X } , A ( ω ) is said to be global random point attractor.
The next theorem, dedicated to verifying the existence of random attractors for SEBT1 directly, can be derived from Theorem 3.11 in Crauel et al. [5].
Theorem 5.
Suppose ϕ is an RDS on connected Polish space ( X , d ) , and suppose that ϕ possesses an absorbing set in X and for any nonrandom bounded set B X , lim t + ϕ ( t , θ t ω ) B is relatively compact P -a.s. Then, ϕ possesses uniqueness global random attractors defined by the following
A ( ω ) = B X Λ B ( ω ) ¯
where union is taken over all bounded B X , and Λ B ( ω ) given by
Λ B ( ω ) = s 0 t s ϕ ( t , θ t ω ) B ¯
Furthermore, A ( ω ) is measurable with respect to F and connected. P - a . s . .
The following assertion provides the relationship between random attractors and invariant measures which is important to expound the global dynamics for RDS was stated in [4,8,63,64].
Proposition 2.
When the RDS ϕ possesses global random attractor comply with Definition 3, by the Corollary 4.4 in Crauel [8], this attractors supports every invariant measures. The random point attractor of ϕ given by Definition 3 always supports at lest one invariant measure which even is a invariant Markov measure (Crauel [63], P423; Arnold [4], Theorem 1.6.13 and Theorem 1.7.5). When ϕ is a white noise RDS or SDS, together with the Theorem 3.6 in Crauel [63] give that every invariant Markov measure is supported by the global point attractor. On the other hand, if the global random attractors for ϕ exists, then ϕ also has the global point attractor. For any fixed ω Ω , taking advantage of pullback mechanism [4], follow the proof of Theorem 5.2 in Birnir [64], the global random attractor can be decomposed into two ingredients, one is random basic attractor which supports all stable invariant Markov measures of white noise RDS, the other is random remainder. For the definition of basic and remainder, we refer to Birnir [64].
With the assertion that the RDS possesses a global random attractors, in light of Proposition 2, the investigation on global dynamics of RDS can be accomplished by exploiting the numerical results on the structure of global random attractor.
The next proposition should be used in checking that the RDS ϕ ( t , ω ) satisfies condition “for any nonrandom bounded set B X , lim t + ϕ ( t , θ t ω ) B is relative compact P -a.s.” in Theorem 5.
Proposition 3.
([26]) Suppose B X is any bounded set, let
α ( B ) = inf d : B has a finite cover of diameter < d
then α ( B ) is known as the Kuratowski’s α-measure of non-compactness of B, in short, α-measure of B, which has the following properties.
(i) 
α ( B ) = 0 if and only if B is pre-compact.
(ii) 
α ( A B ) max α ( A ) , α ( B ) .
(iii) 
α ( A + B ) α ( A ) + α ( B ) .
(iv) 
α ( c o B ) = α ( B ) , where α ( c o B ) is the closed convex hull of B.
(v) 
If B 1 B 2 B 3 are nonempty closed sets in X such that α ( B n ) 0 as n , then n 1 B n is nonempty and compact.

3.2. Main Estimates

This subsection presents the main estimates for the solutions of the systems and some lemmas that are momentous to derive the proof for main results.
Firstly, the following properties on Von Karman bracket were given by Proposition 1.4.2 in Chueshov and Lasiecka [26].
If either at least on of u , v , w belongs to H 0 2 ( D ) or all of them belong to ( H 0 1 H 2 ) ( D ) , then
D [ u , v ] w d x = D [ u , w ] v d x = D [ v , w ] u d x .
If v , w ( H 0 1 H 2 ) ( D ) , then
D [ u , v ] w d x = D [ u , w ] v d x , u H 2 ( D ) .
The next Lemma consider the sharp regularity Von Karman bracket.
Lemma 1.
([25]) Suppose D is a bounded domain with regular boundary or a rectangle in R 2 , the following estimates holds
[ u , w ] j θ C u 2 θ + β w 3 j β , u H 2 θ + β ( D ) , w H 3 j β ( D ) ,
where j = 0 , 1 , 0 < β θ < 1 .
[ u , w ] j C u 2 β w 3 j + β , u H 2 β ( D ) , w H 3 j + β ( D ) ,
in which j = 1 , 2 , 0 β < 1 . Let Δ 2 be the inverse operator of Δ 2 with clamped boundary condition, invoking Equations (41) and (42), we have ( u , w ) G ( u , w ) Δ 2 [ u , w ] is bounded mapping from ( H 2 × H 2 ) ( D ) to ( H 3 W 2 , ) ( D ) which satisfies
[ u , w ] 2 C u 1 w 2 , G ( u , v ) 2 , C u 2 v 2 , [ w , G ( u , v ) ] 2 C u 2 v 2 w 2 .
The assertion that the condition “D is a bounded domain with regular boundary or a rectangle in R n ” can be relaxed to “D is a bounded domain satisfies cone property in R n ”, which can be derived using the property of continuation in Sobolev space [65].
The coming Lemma reveals a relationship between Airy function F and in-plane force F 0 .
Lemma 2.
Let F is the solution for system described by Equations (10e) and (10f); F 0 satisfies Equations (10g) and (10h); and then
D Δ F Δ F 0 d x = 0 .
Proof. 
See the Appendix B.1. ☐
According to Lemma 2, the following estimates which play a crucial role in obtaining the existence of global absorbing set for vibration of Von Karman with a arbitrarily small coefficient of weakly damping.
Let Φ = [ u , v ] T be the solution of SAVKPT2 or SMVKPT2, ε satisfies (23), define
Υ ( F , F 0 ) = D [ u , u ] F 0 d x .
Combining Lemma 2 with Equations (10e)–(10h), we have
Υ ( F , F 0 ) = D Δ 2 F ( u ) F 0 d x = Γ Δ F n F 0 d s Γ F 0 n Δ F d s ,
then
Υ ( F , F 0 ) Δ F n 3 2 , Γ F 0 3 2 , 0 , Γ + F 0 n 1 2 , Γ Δ F 1 2 , 0 , C 1 Δ F n 3 2 , Γ F 0 3 2 , 0 , Γ + C 2 F 0 n 1 2 , Γ Δ F 1 2 , 0 , Γ , η Δ F 2 + C ( η , C 1 , C 2 ) F 0 2 2 .
where η is any given positive, C ( η , C 1 , C 2 ) is a positive, which is inversely proportional to η . In addition, invoking the trace theorem [57], we attain
F 0 2 2 C ˜ g 1 3 2 , Γ 2 + C ^ g 2 1 2 , Γ 2 .
In contrast to Inequality (9.1.17) provided by Lemma 9.1.7 in Chueshov and Lasiecka [26], the inequality in Equation (45) give a clear relationship between Υ ( F , F 0 ) and the in-plane force F 0 , which along with Equation (46) indicates that Υ ( F , F 0 ) is determined by the component of the in-plane force on boundary along the direction x , y .
Since the form of damping is weak and the coefficient of it can be arbitrarily small, it is accomplished by verifying the condition in Theorem 5 that “for any nonrandom bounded set B E 1 , lim t + ϕ ( t , θ t ω ) B is relative compact P -a.s.” relies on the stabilization estimation of the considered systems. The results listed in the next two Lemmas are very important to accomplish the stabilization estimation of SAVKPT1 and SMVKPT1.
Lemma 3.
Suppose { u n } n = 1 + L ( ( τ , t ) , ( H 0 2 H 4 ) ( D ) , { u t n } n = 1 + L ( ( τ , t ) , L 2 ( D ) are two given bounded sequences, which are weakly star convergence to u , u t , respectively. Then, for n , m N , the following holds
lim n lim m τ t F ( u n ) + F 0 , u n F ( u m ) + F 0 , u m , u t n u t m d s = 0 .
Lemma 4.
Suppose { u n } n = 1 + L ( ( τ , t ) , ( H 0 2 H 4 ) ( D ) , { u t n } n = 1 + L ( ( τ , t ) , L 2 ( D ) are two given bounded sequences, which are weakly star convergence to u , u t , respectively. Then, for n , m N , the following holds
lim n lim m τ t F ( u n ) + F 0 , u n F ( u m ) + F 0 , u m , u n u m d s = 0 .
Lemma 3 was proposed by Khanmamedov [35]. With a similar treatment, the results in Lemma 4 can be attained. Hence, it is omitted it here.
It is noticed that the A 1 in Equations (15) and (20) are the same, then the estimates given in the next Lemma, which can be obtained by simple computation, should be employed in both cases of additive white noise and multiplicative white noise.
Lemma 5.
for any U = [ U 1 , U 2 ] T E 1 , the ensuring holds
( A 1 U , U ) E 1 ε 2 U E 1 2 k 0 2 U 2 ,
where 0 < ε < ε 0 , 0 < k 0 α ε and ε , ε 0 satisfies Equation (23).
The following Lemmas are used to prove the existence and uniqueness of solutions for SAVKP and SMVKP by semigroup theory.
Lemma 6.
A 1 L ( E 2 , E 1 ) , A 1 1 L ( E 1 , E 1 ) .
Proof. 
See Appendix B.2. ☐
Base on the Lemma 6, the ensuring results focus on eigenvalues of A 1 can be verified.
Lemma 7.
The eigenvalues of A 1 are as follows.
λ k ( + ) = α + α 2 4 Λ k 2 , λ k ( ) = α α 2 4 Λ k 2 , k = 1 , , n , .
Proof. 
See Appendix B.3. ☐

3.2.1. Main Estimates Only Be Valid in Situation of Additive White Noise

This subsubsection is to give main estimates that only be valid in additive white noise case.
For any given T > 0 , τ R , for U , V E 1 , U E 1 M , V E 1 M , by Reference [26], we have there exists constant M = C ( M ) , such that
[ F ( U 1 ) + F 0 , U 1 ] [ F ( V 1 ) + F 0 , V 1 ] C ( M ) U V E 1
along with Equation (50), we have
F 3 ( U ) F 3 ( V ) E 1 C 1 ( T , τ , ω , M ) U V E 1
where C 1 ( T , τ , ω , M ) is a constant.
Let the value of η in Equation (45) be 7 6 , 1 2 , 1 4 , respectively, we find
3 k ε 2 Υ ( F , F 0 ) Υ ( F , F 0 ) 7 k ε 4 Δ F 2 + 3 k ε 2 C ( 7 6 , C 1 , C 2 ) F 0 2 2 ,
k Υ ( F , F 0 ) Υ ( F , F 0 ) k 2 Δ F 2 + k 2 C ( 1 2 , C 1 , C 2 ) F 0 2 2 ,
k ¯ | σ z ( θ t ω ) | q k Υ ( F , F 0 ) k ¯ | σ z ( θ t ω ) | q k 4 Δ F 2 + k ¯ k | σ z ( θ t ω ) | q C ( 1 4 , C 1 , C 2 ) F 0 2 2 .
Applying Equations (39) and (40) and Lemma 1, we get
D [ F ( u ) , u ] v d x = D [ u t , u ] F ( u ) d x + ε D [ u , u ] F ( u ) d x σ z ( θ t ω ) D [ q , u ] F ( u ) d x 1 2 D ( Δ 2 F ( u ) ) t F ( u ) d x ε D Δ 2 F ( u ) F ( u ) d x + | σ z ( θ t ω ) | D [ q , u ] F ( u ) d x 1 4 d d t Δ F ( u ) 2 ε Δ F ( u ) 2 + | σ z ( θ t ω ) | q 1 k A 1 2 u 2 + k 4 Δ F ( u ) 2
and
D [ F 0 , u ] v d x = D [ u t , u ] F 0 d x + ε D [ u , u ] F 0 d x σ z ( θ t ω ) D [ q , u ] F 0 d x 1 2 D ( Δ 2 F ( u ) ) t F 0 d x ε D Δ 2 F ( u ) F 0 d x + | σ z ( θ t ω ) | D [ q , u ] F 0 d x = 1 2 d d t Υ ( F , F 0 ) ε Υ ( F , F 0 ) + η 2 A u 2 + 1 2 η | σ z ( θ t ω ) | q 1 F 0 2 2 .
where η is positive constant. Set η = ε 4 k , then
D [ F 0 , u ] v d x 1 2 Υ ( F , F 0 ) ε Υ ( F , F 0 ) + ε 8 k A u 2 + 2 k ε | σ z ( θ t ω ) | q 1 F 0 2 2 .
Employing Young’s inequality, it can be obtained that
2 σ A 1 2 q ( x ) z ( θ t ω ) , A 1 2 u ε 4 A 1 2 u 2 + 4 σ 2 | z ( θ t ω ) | 2 ε A 1 2 q ( x ) 2 ,
2 σ ( μ α + ε ) q ( x ) z ( θ t ω ) , v k 0 2 v 2 + σ ( μ α + ε ) z ( θ t ω ) 2 k 0 q ( x ) 2 .
The next lemma shows S A ( t , ω ) possesses global random absorbing set in E 1 . In addition, the expectation of radius of this random set is bounded.
Lemma 8.
For any given non random bounded set B E 1 , there exists T B ( ω ) 0 , such that, for t T B ( ω ) , the following holds P -a.s.
S A ( t , θ t ω ) B E 1 r ( ω )
and E ( ρ ( ω ) ) M < . In which ρ ( ω ) and r ( ω ) are formulated by Equations (26) and (27), respectively.
Proof. 
See Appendix B.4. ☐
Let K ( ω ) = B ( 0 , r ( ω ) ) ; based on Lemma 8, it can be obtained that, for any non-random bounded set B E 1 , there exists T B ( ω ) > 0 , such that, for T 1 T B ( ω ) , the following holds
S A ( T 1 , θ T 1 ω ) B K ( ω ) .
The coming lemma plays a key role in verifying the condition in Theorem 5 that “for any nonrandom bounded set B E 1 , lim t + ϕ ( t , θ t ω ) B is relative compact P -a.s.” which is significant to prove the existence of the global random attractors for SAVKPT1.
Lemma 9.
Suppose { U n ( 0 ) } n = 1 + is any given bounded sequence in any given non-random bounded set B E 1 , T 0 [ 0 , ) is a any given constant, n , m N . Then, for ε ¯ > 0 , there exist T B ( ω , ε ¯ ) T B ( ω ) and T = max T B ( ω , ε ¯ ) , T 0 such that
S A ( T , θ T ω ) U m ( 0 ) S A ( T , θ T ω ) U n ( 0 ) E 1 ε ¯ + K ε ¯ , B , T , ω A ( U m ( 0 ) , U n ( 0 ) ) ,
and
lim inf n + lim inf m + K ε ¯ , B , T , ω A ( U m ( 0 ) , U n ( 0 ) ) = 0 .
Proof. 
See the Appendix B.5. ☐

3.2.2. Main Estimates Only Be Valid in Situation of Multiplicative White Noise

For any given T > 0 , τ R , U , V E 1 , U E 1 M , V E 1 M , M is constant, take account into z ( θ t ω ) is continuous in [ τ , τ + T ] , merging with Equation (50), we have that there exists a constant C ( T , τ , ω , M ) such that
F 2 ( U , θ t ω ) + F 3 ( U ) F 2 ( V , θ t ω , ) F 3 ( V ) E 1 C ( T , τ , ω , M ) U V E 1 ,
which indicates that F 2 ( U , θ t ω ) + F 3 ( U ) : E 1 E 1 satisfies Lipschtiz condition.
Let Φ = [ Φ 1 , Φ 2 ] T = [ u , v ] T E 2 ; the following estimates should be frequently used.
σ z ( θ t ω ) u , u 2 | σ z ( θ t ω ) | u 2 2 ,
( σ z ( θ t ω ) v , v ) | σ z ( θ t ω ) | v 2 2 ,
σ 2 z 2 ( θ t ω ) u , v σ 2 | z ( θ t ω ) | 2 2 λ 1 A u 2 + v 2 .
D [ F ( u ) , u ] v d x = 1 2 D ( Δ 2 F ( u ) ) t F ( u ) d x ε D Δ 2 F ( u ) F ( u ) d x + σ z ( θ t ω ) D [ u , u ] F ( u ) d x = 1 4 d d t Δ F ( u ) 2 ε Δ F ( u ) 2 + σ z ( θ t ω ) Δ F ( u ) 2
D [ F 0 , u ] v d x = D [ u t , u ] F 0 d x + ε D [ u , u ] F 0 d x σ z ( θ t ω ) D [ u , u ] F 0 d x = 1 2 d d t Υ ( F , F 0 ) ε Υ ( F , F 0 ) + σ z ( θ t ω ) Υ ( F , F 0 ) d x
where Υ ( F , F 0 ) is defined by Equation (44).
Let the value of η in Equation (45) be 1 8 and 3 2 , respectively, and the following can be satisfied
4 2 | σ z ( θ t ω ) | + σ 2 | z ( θ t ω ) | 2 Λ 1 k Υ ( F , F 0 ) 1 2 2 | σ z ( θ t ω ) | + σ 2 | z ( θ t ω ) | 2 Λ 1 k Δ F 2
+ 4 2 | σ z ( θ t ω ) | + σ 2 | z ( θ t ω ) | 2 Λ 1 k C ( 1 8 , C 1 , C 2 ) F 0 2 2 ,
ε k Υ ( F , F 0 ) 3 2 ε k Δ F 2 + ε k C ( 3 2 , C 1 , C 2 ) F 0 2 2 ,
The coming Lemma asserts that S M ( t , ω ) possesses global random absorbing set in E 1 .
Lemma 10.
Let B E 1 be any given non-random bounded set, then there exists T B ( ω ) 0 such that
S M ( t , θ t ω ) B E 1 r ( ω ) , t T B ( ω ) , P - a . s .
and
E ( ρ ( ω ) ) M < .
where ρ ( ω ) , r ( ω ) is formulated by Equations (30) and (32), respectively. M is given in Theorem 4.
Proof. 
See Appendix B.6. ☐
Let K ( ω ) = B ( 0 , r ( ω ) ) , invoke Lemma 10, we have that for any given bounded set B E 1 , there exists T B ( ω ) > 0 , such that for T 1 T B ( ω ) , the coming holds
S M ( T 1 , θ T 1 ω ) B K ( ω ) .
Similar to the results presented in Lemma 9, the stabilization estimation of SMVKPT1 provided in the next Lemma can be derived.
Lemma 11.
Suppose { U n ( 0 ) } n = 1 + is any given bounded sequence in any given non-random bounded set B E 1 , T 0 [ 0 , ) is a any given constant, n , m N . Then, for ε ¯ > 0 , there exist T B ( ω , ε ¯ ) T B ( ω ) and T = max T B ( ω , ε ¯ ) , T 0 such that
S M ( T , θ T ω ) U m ( 0 ) S M ( T , θ T ω ) U n ( 0 ) E 1 ε ¯ + K ε ¯ , B , T , ω M ( U m ( 0 ) , U n ( 0 ) ) ,
and
lim inf n + lim inf m + K ε ¯ , B , T , ω M ( U m ( 0 ) , U n ( 0 ) ) = 0 .
Proof. 
See Appendix B.7. ☐

4. Proofs for Main Results

4.1. Proof for Theorem 1

Along with Lemmas 5–7, we have A 1 can induce a linear semigroup of contractions denoted by e A 1 t , t R + . Invoking Equation (51), we get F 2 ( θ t ω ) + F 3 ( U ) : E 1 E 1 satisfies the Lipschtiz condition. z ( θ t ω ) is continuous in [ τ , t ] . Hence, according to Theorem 2.5.1 in Reference [66], we conclude that SAVKP possesses uniqueness (mild) solution with the form
Φ A ( t , τ , ω ; Φ τ ) = e A ( t τ ) Φ τ + τ t e A ( t s ) F 2 ( θ s ω ) + F 3 ( Φ ( s ) ) d s .
where, t τ , t , τ R . Thus,
Φ A ( t τ , 0 , θ τ ω ; Φ τ ) = e A ( t τ ) Φ τ + 0 t τ e A ( t τ s ) F 2 ( θ s + τ ω ) + F 3 ( Φ ( s ) ) d s .
Let s + τ = T , we have
e A ( t τ ) Φ τ + 0 t τ e A ( t τ s ) F 2 ( θ s + τ ω ) + F 3 ( Φ ( s ) ) d s = e A ( t τ ) Φ τ + τ t e A ( t T ) F 2 ( θ T ω ) + F 3 ( Φ ( T τ ) ) d T
Therefore,
Φ A ( t τ , 0 , θ τ ω ; Φ τ ) = e A ( t τ ) Φ τ + τ t e A ( t T ) F 2 ( θ T ω ) + F 3 ( Φ ) ( T τ ) d T .
Hence, Φ A ( t τ , 0 , θ τ ω ; Φ τ ) is the value of a solution for Equation (10a) with initial value t = τ , Φ = Φ τ at time t. Then
Φ A ( t τ , 0 , θ τ ω ; Φ τ ) = Φ A ( t , τ , ω ; Φ τ ) .
Furthermore, the solution mapping φ A ( t , τ , ω ) satisfies
φ A ( t , τ , ω ) = φ A ( t τ , 0 , θ τ ω ) .
On the other hand, Φ A ( τ , τ , ω ; Φ τ ) = Φ τ which means φ A ( τ , τ , ω ) = 0 . Thus, the proof is completedf.

4.2. Proof for Theorem 2

It follows from Lemma 8 that the RDS S ( t , ω ) possesses global absorbing set in E 1 . The estimation of expectation of radius of the global absorbing set is obtained in the proof for Lemma 8 if the condition that “for any nonrandom bounded set B E 1 , lim t + S A ( t , θ t ω ) B is relative compact P -a.s.” can be satisfied, which is achieved in the following. Then, according the Theorem 5, we can complete the proof of this theorem.
Suppose B E 1 is any bounded set, let T ˜ 0 T P and
A n B ( ω ) = s n T ˜ 0 S A ( s , θ s ω ) B ¯ , n N ,
From Reference [67], we have
A n + 1 B ( ω ) A n B ( ω ) , n N ,
and A n B ( ω ) , n N is nonempty closed set.
Let A B ( ω ) = lim n A n B ( ω ) , A B ( ω ) = lim T + S A ( T , θ T ω ) B , then
A ¯ B ( ω ) = A ¯ B ( ω )
In the following, we verify A B ( ω ) is relative compact. It is enough to prove
lim T + α K ( S A ( T , θ T ω ) B ) = 0 ,
where α K ( · ) denote the α measure.
Taking advantage of the contradiction method, suppose Equation (70) is not true, then there exists ε ^ 0 > 0 , for T ¯ 1 0 and T ¯ 0 > T ¯ 1 such that
α K ( S A ( T ¯ 0 , θ T ¯ 0 ω ) B ) 3 ε ^ 0 .
Let ε ¯ 0 = ε ^ 0 , T ¯ 1 = T B ( ε ¯ , ω ) , T 0 = T ¯ 0 , T = max T B ( ε ¯ , ω ) , T 0 , obviously, T = T ¯ 0 . Hence, α K ( S ( T ¯ 0 , θ T ¯ 0 ω ) B ) 3 ε ¯ 0 which reveals that there exists sequence { U ( 0 ) } n = 1 B such that
sup τ R S A ( T ¯ 0 , θ T ¯ 0 ω ) U m ( 0 ) S A ( T ¯ 0 , θ T ¯ 0 ω ) U n ( 0 ) E 1 2 ε ¯ 0 , m n , n , m = 1 , 2 ,
which contradicts Equation (59). Thus, the supposed assumption is invalid which means Equation (70) holds. Then, α K A B ( ω ) = 0 which means lim n α K ( A n B ( ω ) ) = 0 , together with (69), by the Proposition 3, we have n 1 A n B ( ω ) is compact. On the other hand, since
n 1 A n B ( ω ) = N 0 n N S A ( n T ˜ 0 , θ n T ˜ 0 ω ) B ¯ ,
Let
Λ B ( ω ) = N 0 n N S A ( n T ˜ 0 , θ n T ˜ 0 ω ) B ¯ ,
then Λ B ( ω ) is compact.
According to Theorem 5, it can be found that the RDS S A ( t , ω ) possesses random attractors A ( ω ) E 1 , ω Ω defined by
A ( ω ) = B E 1 Λ B ( ω ) ¯ .

4.3. Proof for Theorem 3

By Equation (60), we have F 2 ( U , θ t ω ) + F 3 ( U ) : E 1 E 1 satisfies Lipschtiz condition, and the rest of proof is very similar to the Proof for Theorem 1 provided in Section 4.4. Thus, it is omitted here.

4.4. Proof for Theorem 4

This process of proof can be completed with replacing Lemma 8 by Lemma 10, Equation (59) by Equation (68) and the RDS S A ( t , ω ) by RDS S M ( t , ω ) in the proof for Theorem 4. Thus, it is omitted here.

5. Summary and Conclusions

The affirmation that expectation of radius of global random attractors for nonlinear stochastic plates considered in the paper is in direct proportion to the intensity of white noise and F 0 can be derived by Theorem 2 and Theorem 4. On the other hand, it follows from Theorem 4 and the assumption in Equation (28) that, for the clamped Von Karman plate without rotational inertia driven by multiplicative white noise, values for the coefficient of the noise that are too large might result in the non-existence of global random attractor for the system, the phenomenon of which cannot be obtained in the status of additive white noise. The estimate of Equation (45) is derived in Section 3.2 to achieve the existence of global absorbing set for the system, which can also accomplished by Inequality (9.1.17) provided by Lemma 9.1.7 in Chueshov and Lasiecka [26]. However, the inequality in Equation (45) gives a clear relationship between Υ ( F , F 0 ) defined by Equation (44) in Section 3.2 and the in-plane force F 0 , which along with Equation (46) indicates that Υ ( F , F 0 ) is determined by the component of the in-plane force on boundary along the direction x , y . This assertion can expounds the bucking phenomenon from the theoretical results on global random attractors.
Compared with the results on global attractors provided by Chen et al. [36] in situation of multiplicative noise, the investigation on existence of global random attractors for the systems, which are carried out in this paper, do not need the assumption that coefficient of the damping is big enough. This statement can be achieved by the process of proving the existence of global random absorbing set for clamped Von Karman plate without rotational inertia driven by white noise (see Lemma 8 in Section 3.2.1 and Lemma 10 in Section 3.2.2) and attaining the stabilization estimation of the systems(see Lemma 9 in Section 3.2.1 and Lemma 11 in Section 3.2.2). It is noticed that the condition for the existence of global random attractors presented by Chen et al. [36] is so conservative that, once there does not exist a steady ingredient of the loading, which means p = 0 in the governing equations of Von Karman plate considered in [36], then the global random attractors only comprise the random fixed points (even the trivial solution). It is led by the too big coefficient of damping. The results obtained in this paper indicate that there exists global random attractors for the systems with the small coefficient of weak damping. Moreover, the statement that global random attractors are composed of more than a fixed random point is illustrated by the numerical results on global dynamics. On the other hand, the conclusion that multiplicative white noise is more likely to result in the appearance of global D -bifurcation and secondary global D -bifurcation in the motion of clamped Von Karman without rotational inertia than additive noises can also be validated from the numerical results on global dynamics of the systems.
In engineering applications, ( f 1 , f 2 ) are the in-plane components of the external force that act on the edge of the plate. In light of Equations (3) and (4) given in Section 2.1, with Equation (46) attained in Section 3.2.1, the boundary value of F 0 and F 0 n are determined by ( f 1 , f 2 ) . Thus, the aforementioned qualitative results on long time behavior for the two kinds of nonlinear stochastic Von Karman plates are demonstrated by the numerical results on global dynamics. Big enough value of external force that acts on the edge of the plate leads to the appearance of global D -bifurcation in the motion of the simplified Von Karman plate without rotational inertia, which coincides with the phenomena named bucking [1]. Furthermore, the fact mentioned above together with assertion that global D -bifurcation would disappear with increasing the value of intensity of the white noise can give a reasonable explanation of how the random snap-through [68] occurs in motion of the stochastic Von Karman plate and how to eliminate it in the following manner. When the global random attractor only comprise a fixed random point, then random snap-through occurs. Once the global D -bifurcation appears in the motion of the Von Karman plate, the random snap-through will disappears. Alternatively, the large value of the intensity of the white noise results in the occurring of the random snap-through which can be eliminated by increasing the value of external force that act on the edge of the plate. As indicated above, it can be concluded that the random snap-through which seems to be a complex dynamics intuitively is essentially a simple dynamical behavior.

Author Contributions

Conceptualization, Investigation, Writing—Original Draft Preparation, and Writing—Review and Editing, H.C.; Investigation, Supervision, Project Administration, Funding Acquisition, and Writing—Review and Editing, D.C.; Investigation, Writing—Review and Editing, J.J.; Resources; and Writing—Review and Editing, X.F.

Funding

This work was supported by the Key Project of National Natural Science Foundation of China (No. 11732005) and the National Natural Science Foundation of China (No. 91216106, 11472089).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

RDSRandom Dynamical System
SAVKPThe system described by Equations (10a), (10c), (10e)–(10h)
SAVKPT1The system described by Equations (10c), (15) and (10e)–(10h)
SAVAPT2The system described by Equations (18), (10c) and (10e)–(10h)
SMVKPThe system interpreted by Equations (10b), (10c), (10e)–(10h)
SMVKPT1The system interpreted by Equations (20), (10c) and (10e)–(10f)
SMVAPT2The system interpreted by Equations (21), (10c) and (10e)–(10f)

Appendix A. The Model Equation

Suppose u = i = 1 s l i ( t ) w i , F = i = 1 s 0 a i ( t ) w i , s = s 1 + s 2 , s 0 , s 1 , s 2 N , where { w i } i = 1 4 are eigenvectors related to eigenvalues { Λ i } i = 1 k of operator A. Taking the inner product of Equation (10e) by w i , i = 1 , , k in L 2 , we have
a i ( t ) = 1 Λ i [ i = 1 s l i ( t ) w i , i = 1 s l i ( t ) w i ] , w i , i = 1 , , s 0 ,
then
F = i = 1 s 0 1 Λ i [ i = 1 s l i ( t ) w i , i = 1 s l i ( t ) w i ] , w i w i .
Thus, the model equations associated with SAVKP have the following form
d m ̲ = M ̲ 1 m ̲ + M ̲ 2 ( m ̲ , m ¯ ) d t + σ M ̲ 3 d W , d m ¯ = M ¯ 2 m ¯ + M ¯ 2 ( m ̲ , m ¯ ) d t + σ M ¯ 3 d W , t ( τ h h , τ h ] ; m ¯ = m ¯ τ h h , t ( τ h h , τ h ]
where l ̲ = ( l ̲ 1 , , l ̲ k 1 ) T = ( l 1 , , l k 1 ) T represents the low-frequency modal and the other is high-frequency modal l ¯ = ( l ¯ 1 , , l ¯ k 2 ) T = ( l k 1 , , l k 1 + k 2 ) T , l i is the value of l at time i. g ̲ = l ̲ t , g ¯ = l ¯ t , m ̲ = [ l ̲ 1 , g ̲ 1 , , l ̲ k , g ̲ k 1 ] T , m ¯ = [ l ¯ 1 , g ¯ 1 , , l ¯ k , g ¯ k 2 ] T . h is step size of numerical integration. τ N is an undetermined constant.
M ̲ 1 = 0 1 0 0 0 Λ 1 α 0 0 0 0 0 0 0 0 0 0 Λ s 1 α , M ̲ 2 ( m ̲ , m ¯ ) = 0 F 1 ( l ) 0 F s 1 ( l ) , M ̲ 3 = 0 q ( x ) , w 1 0 q ( x ) , w s 1
and
M ¯ 1 = 0 1 0 0 0 Λ k 1 + 1 α 0 0 0 0 0 0 0 0 0 0 Λ k α , M ¯ 2 ( m ̲ , m ¯ ) = 0 F s 1 + 1 ( l ) 0 F s ( l ) , M ¯ 3 ( m ¯ ) = 0 q ( x ) , w s 1 + 1 0 q ( x ) , w s
in which
F ¯ = F ¯ 1 F ¯ s = σ q ( x ) W ˙ , w 1 σ q ( x ) W ˙ , w s F ( l ) = F 1 ( l ) F 2 ( l ) = k [ F , i = 1 s l i w i ] , w 1 [ F 0 , i = 1 s l i w i ] , w 1 [ F , i = 1 s l i w i ] , w s [ F 0 , i = 1 s l i w i ] , w s ,
Furthermore, let
M ̲ ^ 3 = 0 l ̲ 1 0 l ̲ s 1 , M ¯ ^ 3 = 0 l ¯ 1 0 l ¯ s 2 ,
then, respectively, replacing M ̲ 3 , M ¯ 3 in Equation (A1) by M ̲ ^ 3 , M ¯ ^ 3 leads to the modal equations for SMVKP.

Appendix B. Proof of the Lemmas in Section 3

Appendix B.1. Proof of Lemma 2

Since
D F Δ 2 F 0 d x = Γ Δ F 0 n F D ( Δ F 0 ) · F d x , = Γ Δ F 0 n F d s Γ F n Δ F 0 d s + D Δ F Δ F 0 d x ,
together with Equation (10f) gives
D F Δ 2 F 0 d x = D Δ F Δ F 0 d x .
Taking into account Equation (10g), we have
D Δ F Δ F 0 d x = 0 .

Appendix B.2. Proof of Lemma 6

Let E ¯ = ( H 4 × L 2 ) ( D ) , define
a ( U , V ) = < A 1 U , V > ,
where < , > = D · d x . By Cauchy inequality, there exists M 0 such that
< A 1 U , V > = ( ε U 1 + U 2 , U 1 ) + A U 1 + ε ( α ε ) U 1 , U 2 ( ( α ε ) U 2 , U 2 ) , M U E ¯ V E ¯ . U , V E ¯
From the Proofs of Lax–Milgram given by Temam [55], we have
A 1 L ( E ¯ , E ¯ * ) ,
where E ¯ * is the conjugate space of E ¯ . Moreover, Equation (49) yields
( A 1 U , U ) E 1 ε 2 U E 1 2 , U E 2 .
Let
D 0 = U U E ¯ , A 1 U E 1 ,
then D 0 = E 2 . Since E 1 E ¯ * , then (A5) yields A 1 L ( E 2 , E 1 ) .
Furthermore, Invoking Equation (A6), this gives A 1 : E 2 E 1 is injection. On the other hand, since E 1 E ¯ * , by applying (A5), we have A 1 : E 2 E 1 is surjection. By open mapping theorem, we obtain that the inverse of A 1 exists and A 1 1 L ( E 1 , E 1 ) .

Appendix B.3. Proof of Lemma 7

Without loss of generality, suppose [ w k , e k ] T , k N is eigenvector with respect to eigenvalues λ k , then
A 1 [ w k , e k ] T = λ k [ w k , e k ] T ,
Thus,
ε w k + e k = λ k w k ,
A w k + ε ( α ε ) w k ( α ε ) e k = λ k e k .
Equation (A8) yields e k = ( λ k + ε ) w k ; substituting it into Equation (A9), we have
A w k + ε ( α ε ) w k ( α ε ) ( l k + ε ) w k = λ k ( λ k + ε ) w k ,
then
A w k α λ k w k = λ k 2 w k .
Since A w k = Λ k w k , we obtain
Λ k w k α λ k w k = λ k 2 w k .
Hence,
λ k ( + ) = α + α 2 4 Λ k 2 , λ k ( ) = α α 2 4 Λ k 2 .

Appendix B.4. Proof of Lemma 8

Taking the inner product of Equation (18) by Φ in E 1 gives
d d t Φ E 1 2 = 2 ( A 1 Φ , Φ ) E 1 + 2 σ A 1 2 q ( x ) z ( θ t ω ) , A 1 2 Φ 1 2 σ ( μ α + ε ) q ( x ) z ( θ t ω ) , Φ 2 + 2 k [ F ( Φ 1 ) + F 0 , Φ 1 ] , Φ 2 .
According to Lemma 5, we find
2 ( A 1 Φ , Φ ) E 1 ε Φ E 1 2 k 0 Φ 2 ,
It follows from Equations (57) and (58) that the third and fourth term on the right of the inequality sign in Equation (A11) can be controlled by
ε 4 A 1 2 Φ 1 2 + k 0 2 Φ 2 2 + 4 σ 2 | z ( θ t ω ) | 2 ε A 1 2 q ( x ) 2 + σ ( μ α + ε ) q ( x ) z ( θ t ω ) 2 k 0 q ( x ) 2 ,
Since Φ = [ Φ 1 , Φ 2 ] T = [ u , v ] T , utilizing Equations (55) and (56), we have
2 k [ F ( Φ 1 ) + F 0 , Φ 1 ] , Φ 2 k 2 d d t Δ F ( Φ 1 ) 2 2 k ε Δ F ( Φ 1 ) 2 + 2 k ¯ | σ z ( θ t ω ) | q A 1 2 Φ 1 2 + k 4 Δ F ( Φ 1 ) 2 k Υ ( F , F 0 ) 2 k ε Υ ( F , F 0 ) + ε 4 A Φ 1 2 + 4 k 2 ε | σ z ( θ t ω ) | q 1 F 0 2 2 .
as indicated above, we get
d d t Φ E 1 2 3 ε 4 Φ E 1 2 + 4 σ 2 | z ( θ t ω ) | 2 ε A 1 2 q ( x ) 2 + σ ( μ α + ε ) q ( x ) z ( θ t ω ) 2 k 0 q ( x ) 2 k 2 d d t Δ F ( Φ 1 ) 2 2 k ε Δ F ( Φ 1 ) 2 + 2 k ¯ | σ z ( θ t ω ) | q A 1 2 Φ 1 2 + k 4 Δ F ( Φ 1 ) 2 k d d t Υ ( F , F 0 ) 2 k ε Υ ( F , F 0 ) + ε 4 A Φ 1 2 + 4 k 2 ε | σ z ( θ t ω ) | q 1 F 0 2 2 ε 2 Φ E 1 2 k 2 d d t Δ F ( Φ 1 ) 2 k d d t Υ ( F , F 0 ) k ε 4 Δ F ( Φ 1 ) 2 k ε 2 Υ ( F , F 0 ) + 2 k ¯ | σ z ( θ t ω ) | q A 1 2 Φ 1 2 + k 4 Δ F ( Φ 1 ) 2 + 4 k 2 ε | σ z ( θ t ω ) | q 1 F 0 2 2 + 4 σ 2 | z ( θ t ω ) | 2 ε A 1 2 q ( x ) 2 + σ ( μ α + ε ) q ( x ) z ( θ t ω ) 2 k 0 q ( x ) 2 7 k ε 4 ε Δ F ( Φ 1 ) 2 3 k ε 2 | Υ ( F , F 0 ) | .
By means of Equation (52), we catch
d d t Φ E 1 2 k 2 d d t Δ F ( Φ 1 ) 2 k d d t Υ ( F , F 0 ) ε 2 Φ E 1 2 k ε 4 ε Δ F ( Φ 1 ) 2 k ε 2 Υ ( F , F 0 ) + 2 k ¯ | σ z ( θ t ω ) | q A 1 2 Φ 1 2 + k 4 Δ F ( Φ 1 ) 2 + 3 k ε 2 C ( 7 6 , C 1 , C 2 ) F 0 2 2 + σ ( μ α + ε ) q ( x ) z ( θ t ω ) 2 k 0 q ( x ) 2 + 4 σ 2 | z ( θ t ω ) | 2 ε A 1 2 q ( x ) 2 + 4 k 2 ε | σ z ( θ t ω ) | q 1 F 0 2 2 .
It can be derived from Equation (54) that
d d t Φ E 1 2 k 2 d d t Δ F ( Φ 1 ) 2 k d d t Υ ( F , F 0 ) ε 2 Φ E 1 2 k ε 4 ε Δ F ( Φ 1 ) 2 k ε 2 Υ ( F , F 0 ) + 2 k ¯ | σ z ( θ t ω ) | q A 1 2 Φ 1 2 + k 4 Δ F ( Φ 1 ) 2 + k Υ ( F , F 0 ) + σ ( μ α + ε ) q ( x ) z ( θ t ω ) 2 k 0 q ( x ) 2 + 2 k ¯ k | σ z ( θ t ω ) | q C ( 1 4 , C 1 , C 2 ) F 0 2 2 + 4 σ 2 | z ( θ t ω ) | 2 ε A 1 2 q ( x ) 2 + 4 k 2 ε | σ z ( θ t ω ) | q 1 F 0 2 2 + 3 k ε 2 C ( 7 6 , C 1 , C 2 ) F 0 2 2 .
Let
E 1 ( t , τ , ω ) = E ¯ 1 ( Φ 1 , Φ 2 , ω ) = Φ 2 2 + A 1 2 Φ 1 2 + k 2 Δ F ( Φ 1 ) 2 + k Υ ( F , F 0 ) ,
E ( t , τ , ω ) = E ¯ ( Φ 1 , Φ 2 , ω ) = Φ 2 2 + A 1 2 Φ 1 2 + k 2 Δ F ( Φ 1 ) 2 ,
Invoking Equation (53), we have
E ( Φ 1 , Φ 2 , ω ) E 1 ( Φ 1 , Φ 2 , ω ) + k C ( 1 2 , C 1 , C 2 ) F 0 2 2 ,
which along with Equations (25) and (A11) gives
d d t E ¯ 1 ( Φ 1 , Φ 2 , ω ) E 1 2 ε 2 E ¯ 1 ( Φ 1 , Φ 2 , ω ) + 2 k ¯ | σ z ( θ t ω ) | q E ¯ 1 ( Φ 1 , Φ 2 , ω ) + K ¯ .
Setting ϵ = ε 2 and then employing Equation (A13), it can be asserted that, for any give τ R , t τ , the next holds
E 1 ( t , τ , ω ) e ϵ ( t τ ) + τ t 2 k ¯ | σ z ( θ s ω ) | q d s E 1 ( τ , τ , ω ) + τ t e ϵ ( t s ) + s t 2 k ¯ | σ z ( θ s ω ) | q d k K d s .
which means that for any t < 0 ,
E 1 ( 0 , t , ω ) e ϵ t + t 0 2 k ¯ | σ z ( θ s ω ) | q d s E 1 ( t , t , ω ) + t 0 e ϵ t + t 0 2 k ¯ | σ z ( θ s ω ) | q d k K d s ,
Incorporating Equation (A12), we obtain
E ( 0 , t , ω ) e ϵ t + t 0 2 k ¯ | σ z ( θ s ω ) | q d s E 1 ( t , t , ω ) + t 0 e ϵ t + t 0 2 k ¯ | σ z ( θ s ω ) | q d k K d s + k C ( 1 2 , C 1 , C 2 ) F 0 2 2 .
Since Φ = [ Φ 1 , Φ 2 ] T = [ u , v ] T , v = u t + ε u σ q ( x ) z ( θ t ω ) , then
A 1 2 u 2 + u t + ε u 2 ( 0 , t , ω ) e ϵ t + t 0 2 k ¯ | σ z ( θ s ω ) | q d s E ( t , t , ω ) + t 0 e ϵ t + t 0 2 k ¯ | σ z ( θ s ω ) | q d k K d s + e ϵ t + t 0 2 k ¯ | σ z ( θ s ω ) | q d s σ 2 q ( x ) 2 | z ( θ t ω ) | 2 + σ 2 q ( x ) 2 | z ( ω ) | 2
By Equation (34), we have that there exists T 1 < 0 , such that for t < T 1 , the following holds ( P -a.s.)
1 t 0 t 2 k ¯ | σ z ( θ s ω ) | q d s 4 k ¯ | σ | q π μ ,
furthermore, according to Equation (24), we find ϵ 0 ϵ 4 k ¯ | σ | q π μ > 0 , then
ϵ 1 t 0 t 2 k ¯ | σ z ( θ s ω ) | q d s t ϵ 0 t , t < T 1 .
Since z ( ω ) is tempered, we have
e ϵ t + t 0 2 k ¯ | σ z ( θ s ω ) | q d s E ( t , t , ω ) 0 , t ,
and
e ϵ t + t 0 2 k ¯ | σ z ( θ s ω ) | q d s σ 2 q ( x ) 2 | z ( θ t ω ) | 2 0 , t .
Merging with Equation (25), we attain
0 e ϵ t + t 0 2 k ¯ | σ z ( θ s ω ) | q d k K d s T 1 0 e ϵ s + s 0 2 k ¯ | σ z ( θ s ω ) | q d k K 0 d s + T 1 e ϵ 0 s K 0 d s + T 1 0 e ϵ s + s 0 2 k ¯ | σ z ( θ s ω ) | q d k | z ( θ s ω ) | K 1 d s + T 1 e ϵ 0 s | z ( θ s ω ) | K 1 d s . + T 1 0 e ϵ s + s 0 2 k ¯ | σ z ( θ s ω ) | q d k | z ( θ s ω ) | 2 K 2 d s + T 1 e ϵ 0 s | z ( θ s ω ) | 2 K 2 d s .
then ρ ( ω ) defined by Equation (26) is bounded ( P -a.s.).
It follows from Equations (A15) to (A18) that there exists T B ¯ ( ω ) < 0 , such that, for t ¯ T B ¯ ( ω ) , the following holds
A 1 2 u 2 + u t + ε u 2 ( 0 , t ¯ , ω ) ρ ( ω ) .
On the other hand, since
A 1 2 u 2 + u t + ε u 2 ( 0 , t ¯ , ω ) = A 1 2 u 2 + u t + ε u 2 ( t ¯ , θ t ¯ ω ) ,
let t = t ¯ , T B ( ω ) = T B ¯ ( ω ) , then, for any non-random bounded set B E 1 , there exists T B ( ω ) > 0 , such that, for t > T B ( ω ) , the following holds ( P -a.s.)
S A ( t , θ t ω ) B E 1 r ( ω ) .
The rest is intended to estimate the expectation of ρ ( ω ) . Invoking Equation (36), we have
E e ϵ s + s 0 2 k ¯ | σ z ( θ k ω ) | q d k e ϵ s 2 k ¯ | σ | q μ s ,
on the other hand,
E | z ( θ s ω ) | 4 1 2 = Γ ( 5 2 ) μ 2 π 1 2 , E | z ( θ s ω ) | 2 1 2 = Γ ( 3 2 ) μ 2 π 1 2 ,
where Γ ( · ) is Γ function, which together with Equation (A21) and Cauchy equality gives
E e ϵ s + s 0 2 k | σ z ( θ k ω ) | q d k | z ( θ s ω ) | 2 e ϵ s 2 k σ | q μ s Γ ( 5 2 ) μ 2 π 1 2 ,
E e ϵ s + s 0 2 k | σ z ( θ k ω ) | q d k | z ( θ s ω ) | e ϵ s 2 k σ | q μ s Γ ( 3 2 ) μ 2 π 1 2 .
Define
M = M 0 + M 1 + M 2 + M 3 + δ 0 ,
in which
M 0 = T 1 0 e ϵ s 2 k σ | q μ s K 0 d s + e ϵ 0 T 1 K 0 ϵ 0 , M 1 = K 1 Γ ( 3 2 ) μ 2 π 1 2 T 1 0 e ϵ s 2 k σ | q μ s d s + e ϵ 0 T 1 K 2 ϵ 0 π μ , M 2 = K 2 Γ ( 5 2 ) μ 2 π 1 2 T 1 0 e ϵ s 2 k σ | q μ s d s + e ϵ 0 T 1 K 3 2 ϵ 0 μ , M 3 = σ 2 q ( x ) 2 2 μ .
Obviously, M < . In addition, merging with Equations (25) and (A21)–(A23), we have
E ( ρ ( ω ) ) M < .

Appendix B.5. Proof of Lemma 9

For t 0 , let
{ U n } n = 1 + = { ϕ ( τ , t , ω ) U n ( 0 ) } n = 1 + = { [ u n , v ̲ n ] T } n = 1 + , τ [ t , 0 ] ,
By Lemma 8, we find { U n } n = 1 + has a weakly star subsequence in L ( [ t , 0 ] , E 1 ) , still denoted by { U n } n = 1 + .
Taking the inner product of Equation (16) by U n U m in E 1 , which, merging with Equation (5) and v ̲ = u t + ε u , gives
d d t U n U m E 1 2 ε U n U m E 1 2 + 2 K ¯ ( U m , U n ) ,
where
K ¯ ( U m , U n ) = F ( u n ) + F 0 , u n F ( u m ) + F 0 , u m , u t n u t m + F ( u n ) + F 0 , u n F ( u m ) + F 0 , u m , ε ( u n u m ) .
Since ϕ ( 0 , t , ω ) = S ( t , θ t ω ) , t R + , ω Ω , then
U n U m E 1 2 ( t , θ t ω ) e ε t U n ( 0 ) U m ( 0 ) E 1 2 + 2 t 0 e ε s K ¯ ( U m , U n ) ( s ) d s .
Notice that U n ( 0 ) , U m ( 0 ) B , then e ε t U n ( 0 ) U m ( 0 ) E 1 2 0 , ( t ) . Let
T B ( ε ¯ ) = max 1 ε log ε ¯ 2 B E 1 2 , T B ( ω ) .
thus, for ε ¯ > 0 , there exists
T = max T B ( ε ¯ ) , T 0 ,
such that
U n U m E 1 2 ( T , θ T ω ) ε ¯ + 2 T 0 e ε s K ¯ ( U m , U n ) ( s ) d s .
Set
K ε ¯ , B , T , ω A ( U m ( 0 ) , U n ( 0 ) ) = 2 T 0 e ε s K ¯ ( U m , U n ) ( s ) d s ,
Since e ε s L ( [ T , 0 ] ) , along with Equation (A26), and Lemmas 3 and 4, we obtain
lim inf n + lim inf m + K ε ¯ , B , T , ω A ( U m ( 0 ) , U n ( 0 ) ) = 0 ,

Appendix B.6. Proof of Lemma 10

Taking the inner product of Equation (21) by Φ = [ Φ 1 , Φ 2 ] T = [ u , v ] T in E 1 , we have
d d t Φ E 1 2 = 2 ( A 1 Φ , Φ ) E 1 + 2 ( F 2 Φ , θ t ω ) + F 3 ( Φ , θ t ω ) , Φ E 1 .
Invoking Lemma 5, we find
2 ( A 1 Φ , Φ ) E 1 ε Φ E 1 2 k 0 Φ 2 .
Since
2 ( F 2 ( Φ , θ t ω ) + F 3 ( Φ , θ t ω ) , Φ ) E 1 = 2 σ z ( θ t ω ) Φ 1 , Φ 1 2 + 2 ( σ z ( θ t ω ) Φ 2 , Φ 2 ) + 2 σ 2 z 2 ( θ t ω ) Φ 1 , Φ 2 + 2 ( N ( Φ 1 ) , Φ 2 )
by Equations (61)–(65), it can be obtained that
( F 2 ( Φ , θ t ω ) + F 3 ( Φ , θ t ω ) , Φ ) E 1 = 2 | σ z ( θ t ω ) | A 1 2 Φ 1 2 + 2 | σ z ( θ t ω ) | Φ 2 2 + σ 2 | z ( θ t ω ) | 2 Λ 1 A 1 2 Φ 1 2 + Φ 2 2 k 2 d d t Δ F ( Φ 1 ) 2 2 k ε Δ F ( Φ 1 ) 2 + 2 k σ z ( θ t ω ) Δ F ( Φ 1 ) 2 k d d t Υ ( F ( Φ 1 ) , F 0 ) 2 k ε Υ ( F ( Φ 1 ) , F 0 ) + 2 k σ z ( θ t ω ) Υ ( F ( Φ 1 ) , F 0 ) .
Substituting Equations (A28)–(A29) into Equation (A27), we obtain
d d t Φ 1 E 1 2 + k 2 Δ F ( Φ 1 ) 2 + k Υ ( F ( Φ 1 ) , F 0 ) ε Φ E 1 2 ε k 2 Δ F ( Φ 1 ) 2 3 ε k 2 Δ F ( Φ 1 ) 2 ε k Υ ( F ( Φ 1 ) , F 0 ) ε k Υ ( F ( Φ 1 ) , F 0 ) 2 | σ z ( θ t ω ) | + σ 2 | z ( θ t ω ) | 2 Λ 1 A 1 2 Φ 1 2 + Φ 2 2 + k Δ F ( Φ 1 ) 2 + 3 k Υ ( F ( Φ 1 ) , F 0 ) + 4 2 | σ z ( θ t ω ) | + σ 2 | z ( θ t ω ) | 2 Λ 1 k | Υ ( F ( Φ 1 ) , F 0 ) | .
along with Equations (66) and (67), we have
d d t Φ 1 E 1 2 + k 2 Δ F ( Φ 1 ) 2 + k Υ ( F ( Φ 1 ) , F 0 ) ε Φ 1 E 1 2 + k 2 Δ F ( Φ 1 ) 2 + k Υ ( F ( Φ 1 ) , F 0 ) + 3 2 | σ z ( θ t ω ) | + σ 2 | z ( θ t ω ) | 2 Λ 1 A 1 2 Φ 1 2 + Φ 2 2 + k 2 Δ F ( Φ 1 ) 2 + k Υ ( F ( Φ 1 ) , F 0 ) + ε k C ( 3 2 , C 1 , C 2 ) F 0 2 2 + 4 k 2 | σ z ( θ t ω ) | + σ 2 | z ( θ t ω ) | 2 Λ 1 C ( 1 8 , C 1 , C 2 ) F 0 2 2
Let
E 1 ( t , τ , ω ) = E ¯ 1 ( Φ 1 , Φ 2 , ω ) = Φ 2 2 + A 1 2 Φ 1 2 + k 2 Δ F ( Φ 1 ) 2 + k Υ ( F , F 0 ) ,
E ( t , τ , ω ) = E ¯ ( Φ 1 , Φ 2 , ω ) = Φ 2 2 + A 1 2 Φ 1 2 + k 2 Δ F ( Φ 1 ) 2 ,
applying Equation (53), we attain
E ¯ ( Φ 1 , Φ 2 , ω ) E ¯ 1 ( Φ 1 , Φ 2 , ω ) + k C ( 1 2 , C 1 , C 2 ) F 0 2 2 ,
according to Equation (A30), we have
d d t E ¯ 1 ( Φ 1 , Φ 2 , ω ) ε E ¯ 1 ( Φ 1 , Φ 2 , ω ) + C 0 ( θ t ω ) E ¯ 1 ( Φ 1 , Φ 2 , ω ) + K
where C 0 ( θ t ω ) is given by Equation (29) and K is denoted by Equation (31).
It follows form Equation (A32) that
E ¯ 1 ( Φ 1 , Φ 2 , ω ) e ε ( t τ ) + τ t C 0 ( θ s ω ) d s E ¯ 1 ( Φ 1 , Φ 2 , ω ) + τ t e ε ( t s ) + s t C 0 ( θ k ω ) d k K d s ,
which means that for any given t τ , t , τ R , ω Ω , the following holds
E 1 ( t , τ , ω ) e ε ( t τ ) + τ t C 0 ( θ s ω ) d s E 1 ( τ , τ , ω ) + τ t e ε ( t s ) + s t C 0 ( θ k ω ) d k K d s .
Since
lim t 1 t 0 t C 0 ( θ s ω ) d s = 3 σ 2 2 μ λ 1 6 σ π μ ,
then, there exits T 1 < 0 such that for any t < T 1 , the ensuing can be satisfied
lim t 1 t 0 t C 0 ( θ s ω ) d s 3 σ 2 μ λ 1 12 σ π μ ,
which together with Equation (28) gives that ϵ = ε 3 σ 2 μ λ 1 12 σ π μ > 0 ; hence, when t < T 1 ,
ε t + t 0 C 0 ( θ s ω ) d s ϵ t < 0 ,
therefore
e ε t + t 0 C 0 ( θ s ω ) d s E ( t , t , ω ) 0 , t . 0 e ε t + t 0 C 0 ( θ s ω ) d s K 0 d t
= T 1 0 e ε t + t 0 C 0 ( θ s ω ) d s K 0 d s + e ϵ T 1 K 0 ϵ , 0 e ε t + t 0 C 0 ( θ s ω ) d s | z ( θ t ω ) | K 1 d t
= T 1 0 e ε t + t 0 C 0 ( θ s ω ) d s | z ( θ t ω ) | K 1 d s + T 1 e ϵ s | z ( θ s ω ) | K 1 d s . t e ε t + t 0 C 0 ( θ s ω ) d s | z ( θ s ω ) | 2 K 2 d t
= T 1 t e ε t + t 0 C 0 ( θ s ω ) d s | z ( θ t ω ) | 2 K 2 d s + T 1 e ϵ s | z ( θ s ω ) | 2 K 2 d s .
Combining Equation (31) with the fact that z ( θ t ω ) is tempered, we have that ρ 0 ( ω ) defined by Equation (29) is a bounded random variable ( P -a.s.). Furthermore, merging with Equation (A31), we find the random variable ρ 0 ( ω ) in Equation (30) is also bounded.
Incorporating Equation (A40), it can be derived that. for any initial value Φ B , there exists T ¯ B ( ω ) < 0 such that
E ( 0 , t , ω ) ρ ( ω ) , t T ¯ B ( ω ) .
Based on the relationship between v and u, we have that. for any initial value [ u , v ] T E 1 , the following holds
A 1 2 u 2 + u t + ε u 2 E 1 ( 0 , t , ω ) | σ z ( ω ) | 2 Λ 1 + 1 ρ ( ω ) , t T ¯ B ( ω ) .
Let t = t ¯ , T B ( ω ) = T ¯ B ( ω ) ; along with Equation (32), it can be asserted that, for any given non-random bounded set B E 1 , there exists T B ( ω ) > 0 such that for t > T B ( ω ) , the following holds P -a.s.
S ( t , θ t ω ) B E 1 r ( ω ) .
Thus, S ( t , ω ) possesses global random absorbing set in E 1 .
The remainder of this proof is intended to estimate the expectation of ρ ( ω ) . It follows from Equations (35), (36) and (29) that
E e ϵ s + s 0 C 0 ( θ k ω ) d k e ϵ s 6 | σ | μ s 3 σ 2 μ Λ 1 s ,
Merging with Equation (A21) and Cauchy inequality, we have
E e ε s + s 0 C 0 ( θ k ω ) d k | z ( θ s ω ) | 2 e ϵ s 6 | σ | μ s 3 σ 2 μ Λ 1 s Γ ( 5 2 ) μ 2 π 1 2 ,
E e ε s + s 0 C 0 ( θ k ω ) d k | z ( θ s ω ) | e ϵ s 6 | σ | μ s 3 σ 2 μ Λ 1 s Γ ( 3 2 ) μ 2 π 1 2 .
Let
M = M 0 + M 1 + M 2 + k C ( 1 2 , C 1 , C 2 ) F 0 2 2 + δ 0 ,
in which
M 0 = T 1 0 e ϵ s 6 | σ | μ s 3 σ 2 μ Λ 1 s K 0 d s + e ϵ 0 T 1 K 0 ϵ , M 1 = K 1 Γ ( 3 2 ) μ 2 π 1 2 T 1 0 e ϵ s 6 | σ | μ s 3 σ 2 μ Λ 1 s d s + e ϵ 0 T 1 K 2 ϵ π μ , M 2 = K 2 Γ ( 5 2 ) μ 2 π 1 2 T 1 0 e ϵ s 6 | σ | μ s 3 σ 2 μ Λ 1 s d s + e ϵ 0 T 1 K 3 2 ϵ μ .
Together with Equations (30), (A41), (A42) and (A37)–(A39), we have
E ( ρ ( ω ) ) M < .

Appendix B.7. Proof of Lemma 11

For t 0 , τ [ t , 0 ] , let { Φ n ( 0 ) } n = 1 + = { T ( θ t ω ) U n ( 0 ) } n = 1 + ,
{ U n } n = 1 + = { ϕ ( τ , t , ω ) U n ( 0 ) } n = 1 + = { [ u n , v ̲ n ] T } n = 1 + , { Φ n } n = 1 + = { φ ( τ , t , ω ) Φ n ( 0 ) } n = 1 + = { [ u n , v n ] T } n = 1 + .
By Lemma 10, we find that { U n } n = 1 + , { Φ n } n = 1 + have weakly star convergence in L ( [ t , 0 ] , E 1 ) , still denoted by { U n } n = 1 + , { Φ n } n = 1 + . Obviously Φ n Φ m satisfies Equation (22). Taking the inner product of Equation (22) by Φ n Φ m in E 1 , we have
d d t Φ n Φ m E 1 2 = 2 ( A 1 ( Φ n Φ m ) , Φ n Φ m ) E 1 + 2 F 2 ( Φ n , θ t ω ) F 2 ( Φ m , θ t ω ) , Φ n Φ m E 1 + 2 F 3 ( Φ n , θ t ω ) F 3 ( Φ m , θ t ω ) , Φ n Φ m E 1
In light of Lemma 5, we have
2 A 1 ( Φ n Φ m ) , Φ n Φ m E 1 ε Φ n Φ m E 1 2 k 0 Φ 2 , n Φ 2 , m .
In term of Equations (61)–(63), we find
2 ( F 2 ( Φ n , θ t ω ) F 2 ( Φ m , θ t ω ) + F 3 ( Φ n , θ t ω ) F 3 ( Φ m , θ t ω ) , Φ ) E 1 2 | σ z ( θ t ω ) | A 1 2 Φ 1 2 + 2 | σ z ( θ t ω ) | Φ 2 2 + σ 2 | z ( θ t ω ) | 2 Λ 1 A 1 2 Φ 1 2 + Φ 2 2 + ( ( N ( u n ) N ( u m ) , v n v m ) ,
Merging with Equation (29), we obtain
2 ( F 2 ( Φ n , θ t ω ) F 2 ( Φ m , θ t ω ) + F 3 ( Φ n , θ t ω ) F 3 ( Φ m , θ t ω ) , Φ ) E 1 1 3 C 0 ( θ t ω ) Φ n Φ m E 1 2 + ( ( N ( u n ) N ( u m ) , v n v m ) .
Since v = u t + ε u t σ u z ( θ t ω ) , then
2 ( F 2 ( Φ n , θ t ω ) F 2 ( Φ m , θ t ω ) + F 3 ( Φ n , θ t ω ) F 3 ( Φ m , θ t ω ) , Φ ) E 1 1 3 C 0 ( θ t ω ) Φ n Φ m E 1 2 + K ¯ ( U m , U n ) + σ z ( θ t ω ) K ^ ( U m , U n ) .
where K ¯ ( U m , U n ) is denoted by Equation (A41)
K ^ ( U m , U n ) = F ( u n ) + F 0 , u n F ( u m ) + F 0 , u m , u n u m .
Substituting Equations (A45) and (A46) into Equation (A44), it can be obtained
d d t Φ n Φ m E 1 2 ε Φ n Φ m E 1 2 + 2 3 C 0 ( θ t ω ) Φ n Φ m E 1 2 + 2 K ¯ ( U m , U n ) + 2 σ z ( θ t ω ) K ^ ( U m , U n ) ,
hence
Φ n Φ m E 1 2 ( 0 , t , ω ) e ε t + t 0 2 3 C 0 ( θ k ω ) d k Φ n Φ m E 1 2 ( t , t , ω ) + 2 t 0 e ε s + s 0 2 3 C 0 ( θ k ω ) d k K ¯ ( U m , U n ) + σ z ( θ t ω ) K ^ ( U m , U n ) d s .
By Equation (A35), we have that for any given T 0 T ¯ , t > T 0 , the following holds
Φ n Φ m E 1 2 ( 0 , t , ω ) e ϵ t Φ n Φ m E 1 2 ( t , t , ω ) + 2 t 0 e ε s + s 0 2 3 C 0 ( θ k ω ) d k K ¯ ( U m , U n ) + σ z ( θ s ω ) K ^ ( U m , U n ) d s ,
and
e ϵ t Φ n Φ m E 1 2 ( t , t , ω ) 0 , t + .
Hence, for any ε ¯ > 0 , let
T B ( ε ¯ ) = max 1 ϵ log ε ¯ 4 | σ z ( ω ) | 2 Λ 1 + 1 B E 1 2 , T B ( ω ) .
then, it follows from Equations (A50) and (A51) that there exists
T = max T B ( ε ¯ ) , T 0 ,
such that
Φ n Φ m E 1 2 ( T , θ T ω ) ε ¯ | σ z ( ω ) | 2 Λ 1 + 1 + 2 T 0 e ε s + s 0 2 3 C 0 ( θ k ω ) d k K ( s , Φ m , Φ n ) d s .
where T B ( ω ) is defined in Lemma 10,
K ( T , Φ m , Φ n ) = T 0 e ε s + s 0 2 3 C 0 ( θ k ω ) d k K ¯ ( U m , U n ) + σ z ( θ s ω ) K ^ ( U m , U n ) d s ,
employ the relationship between { U n } n = 1 + and { Φ n } n = 1 + , we derive
U n U m E 1 2 ( T , θ T ω ) | σ z ( ω ) | 2 Λ 1 + 1 ε ¯ | σ z ( ω ) | 2 Λ 1 + 1 + K ( T , Φ m , Φ n ) , = ε ¯ + 2 | σ z ( ω ) | 2 Λ 1 + 1 K ( T , Φ m , Φ n ) .
Let
K ε ¯ , B , T , ω M ( U m ( 0 ) , U n ( 0 ) ) = 2 | σ z ( ω ) | 2 Λ 1 + 1 K ( T , Φ m , Φ n ) ,
then
U n U m E 1 2 ( T , θ T ω ) ε ¯ + K ε ¯ , B , T , ω M ( U m ( 0 ) , U n ( 0 ) ) .
along with Equation (1), we have e ε s + s 0 2 3 C 0 ( θ k ω ) d k L ( [ T , 0 ] ) , z ( θ t ω ) C ( [ T , 0 ] ) , which together with Lemmas 3 and 4 gives
lim inf n + lim inf m + K ε ¯ , B , T , ω M ( U m ( 0 ) , U n ( 0 ) ) = 0 ,

References

  1. Amabili, M. Nonlinear Vibrations and Stability of Shells and Plates; Cambridge University Press: Cambridge, UK, 2008. [Google Scholar]
  2. Lagnese, J.E. Boundary Stabilization of Thin Plates; SIAM: Philadelphia, PA, USA, 1989. [Google Scholar]
  3. Ciarlet, P.G. Mathematical Elasticity: Theory of Plates; Elsevier: Amsterdam, The Netherlands, 1997. [Google Scholar]
  4. Arnold, L. Random Dynamical Systems; Springer: Berlin, Germany, 1998. [Google Scholar]
  5. Crauel, H.; Flandoli, F. Attractors for random dynamical systems. Probab. Theory Relat. Fields 1994, 100, 365–393. [Google Scholar] [CrossRef]
  6. Schmalfuss, B. Measure Attractors and Stochastic Attractors, Institut for Dynamische Systeme; Technical Report; Bermen University: Bermen, Germany, 1995. [Google Scholar]
  7. Crauel, H.; Debussche, A.; Flandoli, F. Random attractors. J. Dyn. Differ. Equ. 1997, 9, 307–341. [Google Scholar] [CrossRef]
  8. Crauel, H. Global random attractors are uniquely determined by attracting deterministic compact sets. Ann. Mat. Pura Appl. 1999, 176, 57–72. [Google Scholar] [CrossRef]
  9. Bates, P.W.; Lu, K.; Wang, B. Random attractors for stochastic reaction–diffusion equations on unbounded domains. J. Differ. Equ. 2009, 246, 845–869. [Google Scholar] [CrossRef]
  10. Caraballo, T.; Langa, J.A.; Robinson, J.C. Stability and random attractors for a reaction-diffusion equation with multiplicative noise. Discret. Contin. Dyn. Syst. 2000, 6, 875–892. [Google Scholar]
  11. Zhou, S.; Yin, F.; Ouyang, Z. Random attractor for damped nonlinear wave equations with white noise. SIAM J. Appl. Dyn. Syst. 2005, 4, 883–903. [Google Scholar] [CrossRef]
  12. Fan, X.; Chen, H. Attractors for the stochastic reaction–diffusion equation driven by linear multiplicative noise with a variable coefficient. J. Math. Anal. Appl. 2013, 398, 715–728. [Google Scholar] [CrossRef]
  13. Fan, X. Attractors for a damped stochastic wave equation of Sine–Gordon type with sublinear multiplicative noise. Stochastic Anal. Appl. 2006, 24, 767–793. [Google Scholar] [CrossRef]
  14. You, Y. Global Attractor for Nonlinear Wave Equations with Critical Exponent on Unbounded Domain. Appl. Math. Nonlinear Sci. 2016, 2, 581–602. [Google Scholar] [CrossRef]
  15. Pappalardo, C.M.; Guida, D. Use of the Adjoint Method for Controlling the Mechanical Vibrations of Nonlinear Systems. Machines 2018, 6, 19. [Google Scholar] [CrossRef]
  16. Fortuna, L.; Muscato, G. A roll stabilization system for a monohull ship: Modeling, identification, and adaptive control. IEEE Trans. Control Syst. Technol. 1996, 4, 18–28. [Google Scholar] [CrossRef]
  17. Pappalardo, C.M.; Guida, D. System Identification Algorithm for Computing the Modal Parameters of Linear Mechanical Systems. Machines 2018, 6, 12. [Google Scholar] [CrossRef]
  18. Pappalardo, C.M.; Guida, D. On the Computational Methods for Solving the Differential-Algebraic Equations of Motion of Multibody Systems. Machines 2018, 6, 20. [Google Scholar] [CrossRef]
  19. Pappalardo, C.M.; Zhang, Z.; Shabana, A.A. Use of independent volume parameters in the development of new large displacement ANCF triangular plate/shell elements. Nonlinear Dyn. 2018, 91, 2171–2202. [Google Scholar] [CrossRef]
  20. Pappalardo, C.M.; Wallin, M.; Shabana, A.A. A new ANCF/CRBF fully parameterized plate finite element. J. Comput. Nonlinear Dyn. 2017, 12, 031008. [Google Scholar] [CrossRef]
  21. Pappalardo, C.M.; Yu, Z.; Zhang, X.; Shabana, A.A. Rational ANCF thin plate finite element. J. Comput. Nonlinear Dyn. 2016, 11, 051009. [Google Scholar] [CrossRef]
  22. Lasiecka, I. Weak, classical and intermediate solutions to full von Kármán system of dynamic nonlinear elasticity. Appl. Anal. 1998, 68, 121–145. [Google Scholar]
  23. Leiva, H.; Sivoli, Z. Existence, stability and smoothness of a bounded solution for nonlinear time-varying thermoelastic plate equations. J. Math. Anal. Appl. 2003, 285, 191–211. [Google Scholar] [CrossRef]
  24. Abels, H.; Mora, M.G.; Müller, S. The time-dependent von Kármán plate equation as a limit of 3d nonlinear elasticity. Calc. Var. Partial Differ. Equ. 2011, 41, 241–259. [Google Scholar] [CrossRef]
  25. Favini, A.; Horn, M.A.; Lasiecka, I.; Tataru, D. Global existence, uniqueness and regularity of solutions to a von Karman system with nonlinear boundary dissipation. Differ. Integral Equ. 1996, 9, 267–294. [Google Scholar]
  26. Chueshov, I.; Lasiecka, I. Von Karman Evolution Equations: Well-Posedness and Long Time Dynamics; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  27. Chueshov, I.; Lasiecka, I. Attractors for second-order evolution equations with a nonlinear damping. J. Dyn. Differ. Equ. 2004, 16, 469–512. [Google Scholar] [CrossRef]
  28. Chueshov, I.; Lasiecka, I. Inertial manifolds for von Kármán plate equations. Appl. Math. Optim. 2002, 46, 179–206. [Google Scholar]
  29. Lasiecka, I. Uniform decay rates for full von karman system of dynamic theromelasticity with free boundary conditions and partial boundary dissipation. Commun. Partial Differ. Equ. 1999, 24, 1801–1847. [Google Scholar] [CrossRef]
  30. Ryzhkova, I. Dynamics of a thermoelastic von Karman plate in a subsonic gas flow. Zeitschrift für angewandte Mathematik und Physik 2007, 58, 246–261. [Google Scholar] [CrossRef]
  31. Yeoul Park, J.; Hye Park, S. Uniform decay for a von Karman plate equation with a boundary memory condition. Math. Methods Appl. Sci. 2005, 28, 2225–2240. [Google Scholar] [CrossRef]
  32. Park, S.H. Long-time dynamics of a von Karman equation with time delay. Appl. Math. Lett. 2018, 75, 128–134. [Google Scholar] [CrossRef]
  33. Chueshov, I.; Lasiecka, I.; Webster, J.T. Attractors for delayed, nonrotational von Karman plates with applications to flow-structure interactions without any damping. Commun. Partial Differ. Equ. 2014, 39, 1965–1997. [Google Scholar] [CrossRef]
  34. Lasiecka, I.; Webster, J.T. Eliminating flutter for clamped von Karman plates immersed in subsonic flows. Commun. Pure Appl. Anal. 2014, 13, 1935–1969. [Google Scholar] [CrossRef] [Green Version]
  35. Khanmamedov, A.K. Global attractors for von Karman equations with nonlinear interior dissipation. J. Math. Anal. Appl. 2006, 318, 92–101. [Google Scholar] [CrossRef]
  36. Chen, H.; Cao, D.; Jiang, J. Random Attractors for Von Karman Plates Subjected to Multiplicative White Noise Loadings. In Dynamical Systems: Theoretical and Experimental Analysis; Springer: New York, NY, USA, 2016; pp. 59–70. [Google Scholar]
  37. Awrejcewicz, J.; Krysko, A. Analysis of complex parametric vibrations of plates and shells using Bubnov-Galerkin approach. Arch. Appl. Mech. 2003, 73, 495–504. [Google Scholar] [CrossRef]
  38. Alijani, F.; Bakhtiari-Nejad, F.; Amabili, M. Nonlinear vibrations of FGM rectangular plates in thermal environments. Nonlinear Dyn. 2011, 66, 251. [Google Scholar] [CrossRef]
  39. Park, J.S.; Kim, J.H. Thermal postbuckling and vibration analyses of functionally graded plates. J. Sound Vib. 2006, 289, 77–93. [Google Scholar] [CrossRef]
  40. Rashidi, M.; Shooshtari, A.; Bég, O.A. Homotopy perturbation study of nonlinear vibration of Von Karman rectangular plates. Comput. Struct. 2012, 106, 46–55. [Google Scholar] [CrossRef]
  41. Ghayesh, M.H.; Amabili, M.; Païdoussis, M.P. Nonlinear dynamics of axially moving plates. J. Sound Vib. 2013, 332, 391–406. [Google Scholar] [CrossRef]
  42. Ghayesh, M.H.; Farokhi, H. Nonlinear dynamics of microplates. Int. J. Eng. Sci. 2015, 86, 60–73. [Google Scholar] [CrossRef]
  43. Amabili, M. Nonlinear vibrations of viscoelastic rectangular plates. J. Sound Vib. 2016, 362, 142–156. [Google Scholar] [CrossRef]
  44. Wiggins, S. Global Bifurcations and Chaos: Analytical Methods; Springer Science & Business Media: Heidenburg, Germany, 2013; Volume 73. [Google Scholar]
  45. Feng, Z.; Wiggins, S. On the existence of chaos in a class of two-degree-of-freedom, damped, strongly parametrically forced mechanical systems with brokenO (2) symmetry. Zeitschrift für angewandte Mathematik und Physik ZAMP 1993, 44, 201–248. [Google Scholar] [CrossRef]
  46. Zhang, W. Global and chaotic dynamics for a parametrically excited thin plate. J. Sound Vib. 2001, 239, 1013–1036. [Google Scholar] [CrossRef]
  47. Hsu, C. A theory of cell-to-cell mapping dynamical systems. J. Appl. Mech. 1980, 47, 931–939. [Google Scholar] [CrossRef]
  48. Xu, W.; He, Q.; Fang, T.; Rong, H. Stochastic bifurcation in Duffing system subject to harmonic excitation and in presence of random noise. Int. J. Non-Linear Mech. 2004, 39, 1473–1479. [Google Scholar] [CrossRef]
  49. Crauel, H.; Flandoli, F. Additive noise destroys a pitchfork bifurcation. J. Dyn. Differ. Equ. 1998, 10, 259–274. [Google Scholar] [CrossRef]
  50. Esteban, M.; Núñez, E.P.; Torres, F. Bifurcation Analysis of Hysteretic Systems with Saddle Dynamics. Appl. Math. Nonlinear Sci. 2016, 2, 449–464. [Google Scholar]
  51. Arnold, L.; Bleckert, G.; Schenk-Hoppé, K.R. The Stochastic Brusselator: Parametric Noise Destroys Hoft Bifurcation. In Stochastic Dynamics; Springer: New York, NY, USA, 1999; pp. 71–92. [Google Scholar]
  52. Wang, B. Existence, stability and bifurcation of random complete and periodic solutions of stochastic parabolic equations. Nonlinear Anal. Theory Methods Appl. 2014, 103, 9–25. [Google Scholar] [CrossRef] [Green Version]
  53. Schenk-Hoppé, K.R. Random attractors—General properties, existence and applications to stochastic bifurcation theory. Discret. Contin. Dyn. Syst.-A 1998, 4, 99–130. [Google Scholar] [CrossRef]
  54. Caraballo, T.; Langa, J.A.; Robinson, J.C. A stochastic pitchfork bifurcation in a reaction-diffusion equation. In Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences; The Royal Society: London, UK, 2001; Volume 457, pp. 2041–2061. [Google Scholar]
  55. Temam, R. Infinite-Dimensional Dynamical Systems in Mechanics And Physics; Springer: New York, NY, USA, 1997. [Google Scholar]
  56. Keller, H.; Ochs, G. Numerical approximation of random attractors. In Stochastic Dynamics; Springer: New York, NY, USA, 1999; pp. 93–115. [Google Scholar]
  57. Lions, J.L.; Magenes, E. Non-Homogeneous Boundary Value Problems and Applications; Springer Science & Business Media: Paris, France, 2012; pp. 7–8, 112–131, 275–276. [Google Scholar]
  58. Debussche, A.; Temam, R. Some new generalizations of inertial manifolds. Discret. Contin. Dyn. Syst.-A 1996, 2, 543–558. [Google Scholar]
  59. Marion, M.; Temam, R. Nonlinear galerkin methods. SIAM J. Numerical Anal. 1989, 26, 1139–1157. [Google Scholar] [CrossRef]
  60. Multiphysics, A. COMSOL Multiphysics 3.5 a Reference Manual, PDE Mode Equation Based Modeling; Multiphysics Ltd.: Stohkholm, Sweden, 2008. [Google Scholar]
  61. Kloeden, P.; Eckhard, P. Numerical Solution of Stochastic Differential Equations; Springer: Berlin, Germany, 1992. [Google Scholar]
  62. Schmalfuss, B. Measure attractors and random attractors for stochastic partial differential equations. Stochastic Anal. Appl. 1999, 17, 1075–1101. [Google Scholar] [CrossRef]
  63. Crauel, H. Random point attractors versus random set attractors. J. Lond. Math. Soc. 2001, 63, 413–427. [Google Scholar] [CrossRef]
  64. Birnir, B. Basic Attractors and Control; Springer: New York, NY, USA, 2015. [Google Scholar]
  65. Adams, R.A.; Fournier, J.J. Sobolev Spaces; Academic Press: Cambridge, MA, USA, 2003. [Google Scholar]
  66. Zheng, S. Nonlinear Evolution Equations; CRC Press: Boca Raton, FL, USA, 2004. [Google Scholar]
  67. Chepyzhov, V.; Vishik, M. A Hausdorff dimension estimate for kernel sections of non-autonomous evolution equations. Indiana Univ. Math. J. 1993, 42, 1057–1076. [Google Scholar] [CrossRef]
  68. Crisfield, M. A fast incremental/iterative solution procedure that handles “snap-through”. In Computational Methods in Nonlinear Structural and Solid Mechanics; Elsevier: New York, NY, USA, 1981; pp. 55–62. [Google Scholar]
Figure 1. The model for vibration of Von Kaman plate under the random loading.
Figure 1. The model for vibration of Von Kaman plate under the random loading.
Symmetry 10 00315 g001
Figure 2. Invariant measures and random attractors for SAVKP in Case I.: (a) invariant measure and section of global random basic attractor; and (b) global random basic attractor.
Figure 2. Invariant measures and random attractors for SAVKP in Case I.: (a) invariant measure and section of global random basic attractor; and (b) global random basic attractor.
Symmetry 10 00315 g002
Figure 3. Invariant measures and random attractors for SAVKP in Case II: (a) invariant measures and global random basic attractor; (b) section of global random attractor; and (c) global random basic attractor.
Figure 3. Invariant measures and random attractors for SAVKP in Case II: (a) invariant measures and global random basic attractor; (b) section of global random attractor; and (c) global random basic attractor.
Symmetry 10 00315 g003
Figure 4. Invariant measures and random attractors for SAVKP in Case III: (a) invariant measure and section of global random basic attractor; and (b) global random basic attractor.
Figure 4. Invariant measures and random attractors for SAVKP in Case III: (a) invariant measure and section of global random basic attractor; and (b) global random basic attractor.
Symmetry 10 00315 g004
Figure 5. Invariant measures and random attractors for SAVKP in Case IV: (a) invariant measure and section of global random basic attractor; and (b) global random basic attractor.
Figure 5. Invariant measures and random attractors for SAVKP in Case IV: (a) invariant measure and section of global random basic attractor; and (b) global random basic attractor.
Symmetry 10 00315 g005
Figure 6. Invariant measures and random attractors for SMVKP in Case I: (a) invariant measure and section of global random basic attractor; and (b) global random basic attractor.
Figure 6. Invariant measures and random attractors for SMVKP in Case I: (a) invariant measure and section of global random basic attractor; and (b) global random basic attractor.
Symmetry 10 00315 g006
Figure 7. Invariant measures and random attractors for SMVKP in Case II: (a) invariant measure and section of global random basic attractor; and (b) global random basic attractor.
Figure 7. Invariant measures and random attractors for SMVKP in Case II: (a) invariant measure and section of global random basic attractor; and (b) global random basic attractor.
Symmetry 10 00315 g007
Figure 8. Invariant measures and random attractors for SMVKP in Case III: (a) invariant measures and global random basic attractor; (b) section of global random point attractor; and (c) global random point attractor.
Figure 8. Invariant measures and random attractors for SMVKP in Case III: (a) invariant measures and global random basic attractor; (b) section of global random point attractor; and (c) global random point attractor.
Symmetry 10 00315 g008

Share and Cite

MDPI and ACS Style

Chen, H.; Cao, D.; Jiang, J.; Fan, X. Long Time Behavior and Global Dynamics of Simplified Von Karman Plate Without Rotational Inertia Driven by White Noise. Symmetry 2018, 10, 315. https://doi.org/10.3390/sym10080315

AMA Style

Chen H, Cao D, Jiang J, Fan X. Long Time Behavior and Global Dynamics of Simplified Von Karman Plate Without Rotational Inertia Driven by White Noise. Symmetry. 2018; 10(8):315. https://doi.org/10.3390/sym10080315

Chicago/Turabian Style

Chen, Huatao, Dengqing Cao, Jingfei Jiang, and Xiaoming Fan. 2018. "Long Time Behavior and Global Dynamics of Simplified Von Karman Plate Without Rotational Inertia Driven by White Noise" Symmetry 10, no. 8: 315. https://doi.org/10.3390/sym10080315

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop