# Asymptotic model of linearly visco-elastic Kelvin-Voigt type plates via Trotter theory

Yotsawat Terapabkajornded, Somsak Orankitjaroen, Christian Licht

► To cite this version:

Yotsawat Terapabkajornded, Somsak Orankitjaroen, Christian Licht. Asymptotic model of linearly visco-elastic Kelvin-Voigt type plates via Trotter theory. Advances in Difference Equations, 2019, 2019, pp.186. 10.1186/s13662-019-2104-6. hal-02318639

HAL Id: hal-02318639

https://hal.science/hal-02318639

Submitted on 17 Oct 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

RESEARCH

Open Access

<figure><img src="image_3.png" /></figure>

# Asymptotic model of linearly visco-elastic Kelvin-Voigt type plates via Trotter theory

Yotsawat Terapabkajornded<fnref n="1" />, Somsak Orankitjaroen<fnref n="1" />,<fnref n="3" />* and Christian Licht<fnref n="1" />,²,<fnref n="3" />

somsak.ora@mahidol.ac.th

Full list of author information is
available at the end of the article

**Abstract**

We confirm the study (Licht in C. R., Méc. 341:697–700, 2013) devoted to the quasi-static response for a visco-elastic Kelvin–Voigt plate whose thickness goes to zero. For each thickness parameter, the quasi-static response is given by a system of partial differential equations with initial and boundary conditions. Reformulating scaled systems into a family of evolution equations in Hilbert spaces of possible states with finite energy, we use Trotter theory of convergence of semi-groups of linear operators to identify the asymptotic behavior of the system. The asymptotic model we obtain and the genuine one have the same structure except an occurrence of a new state variable. Eliminating the new state variable from our asymptotic model leads to the asymptotic model in (Licht in C. R., Méc. 341:697–700, 2013) which involves an integro-differential system.

**MSC:** 74B99

**Keywords:** Asymptotic model; Thin visco-elastic plates; Kelvin–Voigt visco-elasticity; Trotter theory

## 1 Introduction

In a recent study [2], Licht and Weller promoted an old but not so well-known convergence tool, namely Trotter theory of convergence of semi-groups of linear operators acting on variable Hilbert spaces, in determining the asymptotic modeling in physics of continuous media. They provided various asymptotic models through the lens of Trotter theory as a comparison to other classical methods. One of the models mentioned is a reduction of the dimension problem on thin linear visco-elastic Kelvin–Voigt type plates. Licht [1] studied this problem before in 2013 and derived the asymptotic model with Laplace transform technique. He found that the mechanical behavior of the limit model is no longer of Kelvin–Voigt type, because a term of fading memory appears like in the homogenization problem. However, with Trotter theory of convergence, Licht and Weller suggested that the mechanical behaviors of limit and genuine models are the same except for the appearance of a new state variable. It is well known that to have the same structure in both limit and genuine models is useful for numerical computations.

In this study we aim to justify and confirm their suggestion. We reconsider a reduction of the dimension problem of thin linearly visco-elastic Kelvin–Voigt plates in Sect. 2. By defining a small parameter $\varepsilon$, referred to as the thickness of the plate, each problem is expressed as an initial-boundary value problem (2.1). Under suitable assumptions, we

<figure><img src="image_4.png" /></figure>

<footnote>*Correspondence:</footnote>

<footnote n="1">¹Department of Mathematics,
Faculty of Science, Mahidol
University, Bangkok, Thailand</footnote>

<footnote n="3">³Centre of Excellence in
Mathematics, CHE, Bangkok,
Thailand</footnote>
rescale the problems and reformulate them in terms of a family of transient problems. In Sect. 3 we discuss our convergence tool, that is, Trotter theory of the convergence of semi-groups of linear operators. We then follow a Trotter theory approach by letting $\epsilon$ tend to zero to derive the limit model in Sect. 4. As will be seen in the last section, the limit model contains an additional state variable but with the structure like that of the original one. Eliminating the new state variable from the limit model recovers an additional term of fad-ing memory in the limit model derived by the Laplace transform technique, which implies integro-differential equations involving partial derivatives of the field of displacement.

## 2 Setting the problem

Customarily, we assimilate the physical Euclidean space to $\mathbb{R}^3$; the orthonormal basis of which is denoted by {$e_1, e_2, e_3$}, and for all $\xi = (\xi_1, \xi_2, \xi_3)$ in $\mathbb{R}^3$, $\hat{\xi} := (\xi_1, \xi_2)$. We will study the quasi-static response of a thin linearly visco-elastic Kelvin–Voigt plate subjected to a given load. Like the problem setting in [1], a thin linearly visco-elastic Kelvin–Voigt plate occupies a domain $\Omega^\epsilon := \omega \times (-\epsilon, \epsilon)$, where $\omega$ is a bounded domain in $\mathbb{R}^2$ with a Lipschitz-continuous boundary $\partial\omega$ and $\epsilon$ is the small thickness of the plate. The upper, lower surfaces $\omega \times \{\pm\epsilon\}$ and the lateral part of the plate $\omega \times (-\epsilon, \epsilon)$ are referred to as $\Gamma_\pm^\epsilon$ and $\Gamma_\text{lat}^\epsilon$, respectively. The plate is clamped along a portion of the lateral part $\Gamma_D^\epsilon := \gamma_D \times (-\epsilon, \epsilon)$ where $\gamma_D$ is of positive length. Moreover, it is subjected to body forces of density $f^\epsilon$ and surface forces of density $g^\epsilon$ on the upper, lower surfaces together with the rest of its lateral part $\Gamma_N^\epsilon := \partial\Omega^\epsilon \setminus \Gamma_D^\epsilon$ during the time interval $[0, T]$. We denote $\Gamma_{N\pm}^\epsilon = \Gamma_N^\epsilon \cap \Gamma_\pm^\epsilon$ and $\Gamma_{N,\text{lat}}^\epsilon = \Gamma_N^\epsilon \cap \Gamma_\text{lat}^\epsilon$.

The problem of determining the quasi-static evolution of the plate involves the parameter $\epsilon$ of data and the equations satisfied by the fields of displacement $u^\epsilon$ and stress $\sigma^\epsilon$ are:

$$ egin{cases} 
-\operatorname{div} \sigma^{\epsilon} = f^{\epsilon} & \text{in } \Omega^{\epsilon} \times (0, T), \\ 
\sigma^{\epsilon} = a^{\epsilon} e(u^{\epsilon}) + b^{\epsilon} e(\dot{u}^{\epsilon}) & \text{in } \Omega^{\epsilon} \times (0, T), \\ 
e_{ij}(u^{\epsilon}) := (\partial_i u_j^{\epsilon} + \partial_j u_i^{\epsilon})/2, & \text{on } \Gamma_N^{\epsilon} \times (0, T), \quad u^{\epsilon} = 0 \text{ on } \Gamma_D^{\epsilon} \times (0, T), \\ 
u^{\epsilon}(\cdot, 0) = u_0^{\epsilon} & \text{in } \Omega^{\epsilon}, 
\end{cases} 
\quad (2.1) $$

where $e(u^\epsilon)$, $n^\epsilon$, $a^\epsilon$ and $b^\epsilon$ are the linearized strain tensor, the outward unit normal vector, the elasticity and viscosity tensor fields, respectively, while the upper dot represents the time derivative.

To obtain a simpler but precise enough model we study the quasi-static response for the plate as its thickness tends to zero. Following [1, 3], first we rescale the domain $\Omega^\epsilon$ into a fixed domain $\Omega := \omega \times (-1, 1)$ through the mapping $\pi^\epsilon$:

$$ x = (\hat{x}, x_3) \in \overline{\Omega} \quad \mapsto \quad \pi^\epsilon x = (\hat{x}, \epsilon x_3) := x^\epsilon \in \overline{\Omega}^\epsilon. $$

Next we add two hypotheses to the data. Hypothesis (H₁) is on the real loading which has to be connected to fixed quantities defined on $\overline{\Omega}$, they concern their intensity and
horizontality:

$$ (H_1): \left\{ \begin{array}{l} \exists(f,g) \in C^{0,1}([0, T]; L^2(\Omega; \mathbb{R}^3)) \times C^{1,1}([0, T]; L^2(\Gamma_N; \mathbb{R}^3)) \text{ such that} \\ \hat{f}^\epsilon(\pi^\epsilon x, t) = \epsilon \hat{f}(x, t), \quad f_3^\epsilon(\pi^\epsilon x, t) = \epsilon^2 f_3(x, t), \quad \forall (x, t) \in \Omega \times [0, T], \\ \hat{g}^\epsilon(\pi^\epsilon x, t) = \epsilon^2 \hat{g}(x, t), \quad g_3^\epsilon(\pi^\epsilon x, t) = \epsilon^3 g_3(x, t), \quad \forall (x, t) \in \Gamma_{N\pm} \times [0, T], \\ \hat{g}^\epsilon(\pi^\epsilon x, t) = \epsilon \hat{g}(x, t), \quad g_3^\epsilon(\pi^\epsilon x, t) = \epsilon^2 g_3(x, t), \quad \forall (x, t) \in \Gamma_{N,\text{lat}} \times [0, T], \end{array} \right. $$

while hypothesis (H2) is on the elasticity and viscosity tensors:

$$ (H_2): \left\{ \begin{array}{l} a^\epsilon(\pi^\epsilon x) = a(x), \quad b^\epsilon(\pi^\epsilon x) = b(x) \quad \text{with } a, b \in L^\infty(\Omega; \text{Lin}(\mathbb{S}^3)), \quad \text{and} \\ \exists \kappa_A, \kappa_a > 0 : \kappa_a |e|^2 < a(x) e \cdot e < \kappa_A |e|^2, \quad \forall e \in \mathbb{S}^3, \quad \text{a.e. } x \in \Omega, \\ \exists \kappa_B, \kappa_b > 0 : \kappa_b |e|^2 < b(x) e \cdot e < \kappa_B |e|^2, \quad \forall e \in \mathbb{S}^3, \quad \text{a.e. } x \in \Omega, \end{array} \right. $$

where Lin($\mathbb{S}^3$) denotes the space of linear mapping from $\mathbb{S}^3$ into $\mathbb{S}^3$, $\mathbb{S}^3$ being the space of 3 × 3 symmetric matrices. It will be convenient to write $\mathbb{S}^3 = \hat{\mathbb{S}} \oplus \mathbb{S}^\perp$ with $\hat{\mathbb{S}} := \{e \in \mathbb{S}^3; e_{i3} = 0, 0 \le i \le 3\}$, $\mathbb{S}^\perp := \{e \in \mathbb{S}^3; e_{\alpha\beta} = 0, 1 \le \alpha, \beta \le 2\}$, and to denote the projection of $e$ on $\hat{\mathbb{S}}$ and $\mathbb{S}^\perp$ by $\hat{e}$ and $e^\perp$, respectively. Then one associates a scaled displacement $u_\epsilon = S(\epsilon)u^\epsilon$, defined on $\bar{\Omega} \times [0, T]$, with the true physical displacement $u^\epsilon$, defined on $\bar{\Omega}^\epsilon \times [0, T]$, by

$$ \hat{u}^\epsilon(x^\epsilon, t) = \epsilon \hat{u}_\epsilon(x, t), \quad u_3^\epsilon(x^\epsilon, t) = (u_\epsilon)_3(x, t) \quad \forall (x^\epsilon, t) \in \bar{\Omega}^\epsilon \times [0, T]. $$

Using hypotheses (H1) and (H2) and the scaling of displacement, we can formulate the initial-boundary value problem (2.1) in terms of the scaled variational problem

$$ (\mathcal{P}_\epsilon): \left\{ \begin{array}{l} u_\epsilon \in C^{1,1}([0, T]; H_\epsilon), \quad u_\epsilon(\cdot, 0) = u_\epsilon^0 := S(\epsilon)u_\epsilon^0; \\ \int_{\Omega} [a(x)e(\epsilon, u_\epsilon)(x, t) + be(\epsilon, \dot{u}_\epsilon)(x, t)] \cdot e(\epsilon, v)(x) dx \\ = \int_{\Omega} f(x, t) \cdot v(x) dx + \int_{\Gamma_N} g(x, t) \cdot v(x) d\mathcal{H}_2 \quad \forall v \in H_\epsilon, \end{array} \right. $$

where $H_\epsilon$ is the subspace of $H^1(\Omega, \mathbb{R}^3)$ whose elements have a vanishing trace on $\Gamma_D$, $e_{\alpha\beta}(\epsilon, v) = e_{\alpha\beta}(v)$, $e_{\alpha 3}(\epsilon, v) = \epsilon^{-1}e_{\alpha 3}(v)$, $1 \le \alpha\beta \le 2$, $e_{33}(\epsilon, v) = \epsilon^{-2}e_{33}(v)$ and $\mathcal{H}_2$ is the two-dimensional Hausdorff measure. Subspace $H_\epsilon$ is equipped with the inner product:

$$ (u, v)_{H_\epsilon} = \int_{\Omega} ae(\epsilon, u) \cdot e(\epsilon, v) dx \quad \forall u, v \in H_\epsilon. $$

To deal with non-vanishing external loadings, it suffices to split $u_\epsilon$ into $u_\epsilon^e + u_\epsilon^r$, where $u_\epsilon^e$ solves a static problem associated with the evolution problem under consideration and involving the surface loading only. Then $u_\epsilon^r$ does solve an evolution equation with a second member which is a continuous function of $u_\epsilon^e$ and hence a continuous function of the loading. The static problem is

$$ (\mathcal{P}_\epsilon^e): \left\{ \begin{array}{l} u_\epsilon^e(t) \in H_\epsilon; \\ \int_{\Omega} ae(\epsilon, u_\epsilon^e(t)) \cdot e(\epsilon, v) dx = \int_{\Gamma_N} g(x, t) \cdot v(x) d\mathcal{H}_2 \quad \forall v \in H_\epsilon, \end{array} \right. $$
which has a unique solution by the Lax–Milgram lemma. As $g \mapsto u_{\varepsilon}^e$ is linear continuous from $L^2(\Gamma_N; \mathbb{R}^3)$ to $H_{\varepsilon}$, we have $u_{\varepsilon}^e$ in $C^{1,1}([0, T]; H_{\varepsilon})$. The evolution problem is

$$ (\mathcal{P}_{\varepsilon}^r): \left\{ \begin{array}{l} u_{\varepsilon}^r \in C^{1,1}([0, T]; H_{\varepsilon}); \\ \int_{\Omega}[ae(\varepsilon, u_{\varepsilon}^r) + be(\varepsilon, \dot{u}_{\varepsilon}^r)] \cdot e(\varepsilon, v) dx \\ = \int_{\Omega}f \cdot v dx - \int_{\Omega}be(\varepsilon, \dot{u}_{\varepsilon}^e) \cdot e(\varepsilon, v) dx \quad \forall v \in H_{\varepsilon}, \end{array} \right. $$

which can be reformulated in terms of

$$ \left\{ \begin{array}{ll} \frac{du_{\varepsilon}^r}{dt} - A_{\varepsilon}u_{\varepsilon}^r = q_{\varepsilon} & \text{in } H_{\varepsilon}, \\ u_{\varepsilon}^r(0) = u_{\varepsilon}^0 := u_{\varepsilon}^0 - u_{\varepsilon}^e(0), & \end{array} \right. \qquad (2.2) $$

with $D(A_{\varepsilon}) = H_{\varepsilon}$,

$$ A_{\varepsilon}u_{\varepsilon} = w_{\varepsilon} \in H_{\varepsilon}; \quad \int_{\Omega}[ae(\varepsilon, u_{\varepsilon}) + be(\varepsilon, w_{\varepsilon})] \cdot e(\varepsilon, v) dx = 0 \quad \forall v \in H_{\varepsilon}, $$

$$ (q_{\varepsilon}(t), v)_{H_{\varepsilon}} = \int_{\Omega} f(x,t) \cdot v(x) dx - \int_{\Omega} e(\varepsilon, \dot{u}^e(t)) \cdot e(\varepsilon, v) dx \quad \forall v \in H_{\varepsilon}. $$

Clearly $A_{\varepsilon}$ is bounded, selfadjoint and $m$-dissipative so that the evolution equation (2.2) has a unique solution $u_{\varepsilon}^r$ in $C^{1,1}([0, T]; H_{\varepsilon})$.

The crucial point to prove that the asymptotic model has the same structure as the genuine one is to apply the tool of some convergence which is a not so well-known result, involving two fields, in reduction of dimension similar to two-scale convergence in periodic homogenization. The proposed asymptotic models of $(\mathcal{P}_{\varepsilon}^e)$ and $(\mathcal{P}_{\varepsilon}^r)$ will be given by the following problems:

$$ (\mathcal{P}^e): \left\{ \begin{array}{l} U^e := (u^e, u^{e1}) \in H; \\ \int_{\Omega} a(\hat{e}(u^e) + \partial_3 u^{e1} \otimes_s e_3) \cdot (\hat{e}(v) + \partial_3 v^1 \otimes_s e_3) dx = \int_{\Gamma_N} g \cdot v dH_2 \\ \forall V = (v, v^1) \in H, \end{array} \right. $$

set in a Hilbert space $H := V_{KL}(\Omega) \times H^1(-1, 1; L^2(\omega; \mathbb{R}^3))/L^2(\omega; \mathbb{R}^3)$, with $V_{KL} = \{u \in H_{\Gamma_D}^1(\Omega; \mathbb{R}^3); e(u)^\perp = 0\}$, where $H_{\Gamma_D}^1(\Omega; \mathbb{R}^3)$ is the subspace of $H^1(\Omega; \mathbb{R}^3)$ whose elements have a vanishing trace on $\Gamma_D$, and equipped with the inner product:

$$ (U, V)_H = \int_{\Omega} a(\hat{e}(u) + \partial_3 u^1 \otimes_s e_3) \cdot (\hat{e}(v) + \partial_3 v^1 \otimes_s e_3) dx $$

and

$$ (\mathcal{P}^r): \left\{ \begin{array}{l} U^r := (u^r, u^{r1}) \in C^{1,1}([0, T]; H); \\ \int_{\Omega}[a(\hat{e}(u^r) + \partial_3 u^{r1} \otimes_s e_3) + b(\hat{e}(\dot{u}^r) + \partial_3 \dot{u}^{r1} \otimes_s e_3)] \cdot (\hat{e}(v) + \partial_3 v^1 \otimes_s e_3) dx \\ = \int_{\Omega} f \cdot v dx - \int_{\Omega} b(\hat{e}(\dot{u}^e) + \partial_3 \dot{u}^{e1} \otimes_s e_3) \cdot (\hat{e}(v) + \partial_3 v^1 \otimes_s e_3) dx \\ \forall V = (v, v^1) \in H, \end{array} \right. $$
which can be reformulated in terms of

$$ \begin{cases} \frac{dU^r}{dt} - AU^r = Q \quad \text{in } H, \\ U^r(0) = U^{r0} := U^0 - U^e(0), \end{cases} \qquad (2.3) $$

with $D(A) = H$,

$$ \begin{aligned} AU = W &= (w, w^1) \in H; \int_{\Omega} [a(\hat{e}(u) + \partial_3 u^1 \otimes_s e_3) + b(\hat{e}(w) + \partial_3 w^1 \otimes_s e_3)] \\ & \qquad \cdot (\hat{e}(v) + \partial_3 v^1 \otimes_s e_3) dx = 0 \quad \forall V = (v, v^1) \in H, \\ (Q(t), V)_H &= -\int_{\Omega} (\hat{e}(\dot{u}^e(t)) + \partial_3 \dot{u}^{e1}(t) \otimes_s e_3) \cdot (\hat{e}(v) + \partial_3 v^1 \otimes_s e_3) dx \\ & \qquad + \int_{\Omega} f(x, t) \cdot v(x) dx \quad \forall V = (v, v^1) \in H. \end{aligned} $$

Obviously A is bounded, selfadjoint and m-dissipative. Problem (P$^e$) has a unique solution U$^e$ by the Lax-Milgram lemma. and the evolution problem (2.3) has a unique solution U$^r$ in C$^{1,1}$([0, T]; H).

At this point, we have solutions $(u_e^e)$, $(u_e^r)$, $U^e$, and $U^r$ to the problems $(P_e^e)$, $(P_e^r)$, $(P^e)$, and $(P^r)$, respectively. We claim that the sequences $(u_e^e)$, $(u_e^r)$ converge to $U^e$, $U^r$ in the sense of Trotter. Before we prove our claim, let us discuss this type of convergence in detail.

### 3 Trotter theory of convergence of semi-groups of linear operators

In the framework of Trotter theory of convergence of semi-groups of linear operators, we have a family of evolution equations

$$ (P_n): \begin{cases} \frac{du_n}{dt} - A_n u_n = q_n & \text{in } H_n \\ u_n(0) = u_n^0, \end{cases} $$

governed by a sequence of m-dissipative operators $A_n$, with domains $D(A_n)$, defined in a sequence of Hilbert spaces ($H_n$) with norms $\| \cdot \|_{H_n}$, while initial data $u_n^0$ and $q_n$ are in $D(A_n) \times C^{0,1}([0, T]; H_n)$. Each problem $P_n$ thus has a unique solution in $C^{1,1}([0, T]; H_n) \cap C^{0,1}([0, T]; D(A_n))$. On the other hand we have an evolution equation

$$ (P): \begin{cases} \frac{du}{dt} - Au = q & \text{in } H \\ u(0) = u^0, \end{cases} $$

governed by an m-dissipative operator A, with domain $D(A)$, defined in a Hilbert spaces $H$ with a norm $\| \cdot \|_H$, while initial data $u^0$ and $q$ in $D(A) \times C^{0,1}([0, T]; H)$. Problem (P) also has a unique solution in $C^{1,1}([0, T]; H) \cap C^{0,1}([0, T]; D(A))$. The nature of the spaces $H_n$ and $H$ may be very different, but the sequence of spaces $H_n$ has to converge to $H$. This requirement paves the way to introduce a kind of sequence of representative operators $P_n$ between $H$ and $H_n$.
**Definition 1** We say that a sequence of Hilbert spaces $H_n$ approximates a Hilbert space $H$ in the sense of Trotter, if there exists a representative operator $P_n \in \mathcal{L}(H; H_n)$ satisfying two conditions of uniform continuity ($T_1$) and good energetic representation ($T_2$):

($T_1$) $\exists C > 0$ such that $\|P_n u\|_{H_n} \le C\|u\|_H, \forall u \in H, \forall n,$

($T_2$) $\lim_{n \to \infty} \|P_n u\|_{H_n} = \|u\|_H, \forall u \in H,$

where $\mathcal{L}(H; H_n)$ is a space of all continuous linear mappings from $H$ into $H_n$.

**Definition 2** A sequence $(u_n)$ in $H_n$ converges toward an element $u$ in $H$ in the sense of Trotter if and only if $\lim_{n \to \infty} \|P_n u - u_n\|_{H_n} = 0$.

The fundamental result of Trotter theory of convergence of semi-groups of linear operators is the following.

**Theorem 1** Let $u_n$, $u$ be the solutions of $(P_n)$ and $(\mathcal{P})$, respectively, if

(i) $\|P_n u^0 - u_n^0\|_{H_n} \to 0,$

(ii) $\int_0^T \|P_n q(t) - q_n(t)\|_{H_n} dt \to 0,$

(iii) $\forall y \in X, \text{ dense in } H, \|P_n(I-A)^{-1}y - (I-A_n)^{-1}P_n y\|_{H_n} \to 0$

then, uniformly on $[0, T]$, $\|P_n u(t) - u_n(t)\|_{H_n} \to 0$ and $\|u_n(t)\|_{H_n} \to \|u(t)\|_H$.

Now we are ready to prove our assertion.

## 4 A convergence result

First we have to show that the sequence of Hilbert spaces $H_\epsilon$ approximate the Hilbert space $H$ in the sense of Trotter. Let a representative operator $P_\epsilon$ be defined by $u'_\epsilon = P_\epsilon U$ for all $U = (u, u^1)$ that satisfies

$$ \int_{\Omega} ae(\epsilon, u'_{\epsilon}) \cdot e(\epsilon, v) dx = \int_{\Omega} a(\hat{e}(u) + \partial_3 u^1 \otimes_s e_3) \cdot e(\epsilon, v) dx \quad \forall v \in H_{\epsilon}. \quad (4.1) $$

The operator $P_\epsilon$ satisfies ($T_1$) because

$$ \|P_{\epsilon} U\|_{H_{\epsilon}}^{2} = \int_{\Omega} a(\hat{e}(u) + \partial_{3} u^{1} \otimes_{s} e_{3}) \cdot e(\epsilon, u'_{\epsilon}) dx \le C \|U\|_{H} \|P_{\epsilon} U\|_{H_{\epsilon}}. $$

Now choose $v = u'_\epsilon$ in (4.1). The sequences $(u'_\epsilon)$ and $e(\epsilon, u'_\epsilon)$ are bounded in $H_{\Gamma_D}^1(\Omega; \mathbb{R}^3)$ and $L^2(\Omega; \mathbb{S}^3)$, respectively. Then there exist $(u', \kappa')$ in $H_{\Gamma_D}^1(\Omega; \mathbb{R}^3) \times L^2(\Omega; \mathbb{S}^3)$ such that $(u'_\epsilon, e(\epsilon, u'_\epsilon))$ weakly converge toward $(u', \kappa')$ in $H_{\Gamma_D}^1(\Omega; \mathbb{R}^3) \times L^2(\Omega; \mathbb{S}^3)$. It follows that $\hat{e}(\epsilon, u'_\epsilon)$ weakly converge to $\hat{e}(u') = \hat{\kappa}'$ in $L^2(\Omega; \mathbb{S}^3)$. Now we are in the position to introduce in reduction of dimension problem a kind of convergence involving a couple of limit fields similar to the two-scaled convergence in periodic homogenization problem. There exists $u'^1$ in $H^1(-1, 1; L^2(\omega; \mathbb{R}^3))/L^2(\omega; \mathbb{R}^3)$ such that $e(\epsilon, u'^1) \perp$ weakly converge to $\kappa'^1 = \partial_3 u'^1 \otimes_s e_3$ in $L^2(\Omega; \mathbb{S}^3)$. Hence, $\kappa' = \hat{e}(u') + \partial_3 u'^1 \otimes_s e_3$. Next we choose $v = (w + \epsilon(\hat{w}^1, \epsilon w^3))$ where $w$ in $V_{KL}, w^1$ in $H^1(-1, 1; L^2(\omega; \mathbb{R}^3))/L^2(\omega; \mathbb{R}^3)$ in (4.1). As $\epsilon$ tends to zero, (4.1) becomes

$$ (U', W)_H = (U, W)_H \quad \forall W \in H, $$

where $U' = (u', u'^1)$, $U = (u, u^1)$, and $W = (w, w^1)$. It follows that $U' = U$.
Hence,

$$
\begin{align*}
\lim_{\epsilon \to 0} \| P_{\epsilon} U \|_{H_{\epsilon}}^{2} &= \lim_{\epsilon \to 0} \int_{\Omega} a e(\epsilon, u'_{\epsilon}) \cdot e(\epsilon, u'_{\epsilon}) \, dx \\
&= \lim_{\epsilon \to 0} \int_{\Omega} a(\hat{e}(u) + \partial_{3} u^{1} \otimes_{s} e_{3}) \cdot e(\epsilon, u'_{\epsilon}) \, dx \\
&= \int_{\Omega} a(\hat{e}(u) + \partial_{3} u^{1} \otimes_{s} e_{3}) \cdot (\hat{e}(u') + \partial_{3} u'^{1} \otimes_{s} e_{3}) \, dx \\
&= \int_{\Omega} a(\hat{e}(u) + \partial_{3} u^{1} \otimes_{s} e_{3}) \cdot (\hat{e}(u) + \partial_{3} u^{1} \otimes_{s} e_{3}) \, dx = \| U \|_{H}^{2}.
\end{align*}
$$

That is, operator $P_\epsilon$ satisfies $(T_2)$. Thus, the sequence of Hilbert spaces $H_\epsilon$ approximates the Hilbert space $H$ in the sense of Trotter.

To prove our convergence result, the following lemma is important.

**Lemma 1** For all $z_\epsilon$ in $H_\epsilon$, $Z$ in $H$,

$$
\|P_{\epsilon}Z - z_{\epsilon}\|_{H_{\epsilon}} \to 0 \iff e(\epsilon, z_{\epsilon}) \to \hat{e}(z) + \partial_{3}z^{1} \otimes_{s} e_{3} \text{ in } L^{2}(\Omega; \mathbb{S}^{3}).
$$

*Proof* Observe that

$$
\begin{align*}
\| P_{\epsilon} Z - z_{\epsilon} \|_{H_{\epsilon}}^{2} &= \| P_{\epsilon} Z \|_{H_{\epsilon}}^{2} - 2 (P_{\epsilon} Z, z_{\epsilon})_{H_{\epsilon}} + \| z_{\epsilon} \|_{H_{\epsilon}}^{2} \\
&= \| P_{\epsilon} Z \|_{H_{\epsilon}}^{2} - 2 \int_{\Omega} a e(\epsilon, z'_{\epsilon}) \cdot e(\epsilon, z_{\epsilon}) \, dx + \| z_{\epsilon} \|_{H_{\epsilon}}^{2} \\
&= \| P_{\epsilon} Z \|_{H_{\epsilon}}^{2} - 2 \int_{\Omega} a (\hat{e}(z) + \partial_{3} z^{1} \otimes_{s} e_{3}) \cdot e(\epsilon, z_{\epsilon}) \, dx + \| z_{\epsilon} \|_{H_{\epsilon}}^{2} \\
&= \| P_{\epsilon} Z \|_{H_{\epsilon}}^{2} - 2 (\hat{e}(z) + \partial_{3} z^{1} \otimes_{s} e_{3}, e(\epsilon, z_{\epsilon}))_{L^{2}(\Omega; \mathbb{S}^{3})} + \| z_{\epsilon} \|_{H_{\epsilon}}^{2} \\
&= \| P_{\epsilon} Z \|_{H_{\epsilon}}^{2} + \| e(\epsilon, z_{\epsilon}) - (\hat{e}(z) + \partial_{3} z^{1} \otimes_{s} e_{3}) \|_{L^{2}(\Omega; \mathbb{S}^{3})}^{2} - \| Z \|_{H}^{2},
\end{align*}
$$

we have

$$
\| P_{\epsilon} Z - z_{\epsilon} \|_{H_{\epsilon}}^{2} - (\| P_{\epsilon} Z \|_{H_{\epsilon}}^{2} - \| Z \|_{H}^{2}) = \| e(\epsilon, z_{\epsilon}) - (\hat{e}(z) + \partial_{3} z^{1} \otimes_{s} e_{3}) \|_{L^{2}(\Omega; \mathbb{S}^{3})}^{2},
$$

and the statement follows from condition (T_2).

<figure><img src="image_5.png" /></figure>

To prove that $(u_{\epsilon}^{e})$ converges to $U^{e} = (u^{e}, u^{e1})$ in the sense of Trotter, it is enough to show the strong convergence of the scaled strain tensor according to Lemma 1. Like the proof of condition (T_2) of $P_{\epsilon}$, we have $e(\epsilon, u_{\epsilon}^{e}) \to \hat{e}(u^{e}) + \partial_{3}u^{e1} \otimes_{s} e_{3}$ in $L^{2}(\Omega; \mathbb{S}^{3})$. This weak convergence and the fact that $u_{\epsilon}^{e}, U^{e}$ solve $(\mathcal{P}_{\epsilon}^{e}), (\mathcal{P}^{e})$ lead to

$$
\begin{align*}
& \| e(\epsilon, u_{\epsilon}^{e}) - (\hat{e}(u^{e}) + \partial_{3}u^{e1} \otimes_{s} e_{3}) \|_{L^{2}(\Omega; \mathbb{S}^{3})}^{2} \\
&= \int_{\Gamma_{\mathbb{N}}} g \cdot u_{\epsilon}^{e} d\mathcal{H}_{2} - 2 \int_{\Omega} a e(\epsilon, u_{\epsilon}^{e}) \cdot (\hat{e}(u) + \partial_{3}u^{e1} \otimes_{s} e_{3}) \, dx + \int_{\Gamma_{\mathbb{N}}} g \cdot u^{e} d\mathcal{H}_{2} \\
&\to 2 \int_{\Gamma_{\mathbb{N}}} g \cdot u^{e} d\mathcal{H}_{2} - 2 \int_{\Omega} a(\hat{e}(u) + \partial_{3}u^{e1} \otimes_{s} e_{3}) \cdot (\hat{e}(u) + \partial_{3}u^{e1} \otimes_{s} e_{3}) \, dx = 0.
\end{align*}
$$
Therefore, $e(\varepsilon, u_{\varepsilon}^e)$ converge strongly to $\hat{e}(u^e) + \partial_3 u^{e1} \otimes_s e_3$ in $L^2(\Omega; \mathbb{S}^3)$. Note that this convergence is uniform on $[0, T]$.

To prove that $(u_{\varepsilon}^r)$ converges uniformly on $[0, T]$ to $U^r = (u^r, u^{r1})$ in the sense of Trotter, we have to prove the three conditions of Theorem 1.

Concerning initial data $u_{\varepsilon}^{r0}$ and $U^{r0}$, we make an additional assumption on $u_{\varepsilon}^0$. That is,

$$ (H_3): \quad \exists U^0 = (u^0, u^{01}) \text{ in } U^e(0) + D(A) \text{ such that } \|P_{\varepsilon} U^0 - u_{\varepsilon}^0\|_{H_{\varepsilon}} \to 0. $$

Hence,

$$ \begin{aligned} \|P_{\varepsilon} U^{r0} - u_{\varepsilon}^{r0}\|_{H_{\varepsilon}} &= \|P_{\varepsilon} U^0 - P_{\varepsilon} U^e(0) - u_{\varepsilon}^0 + u_{\varepsilon}^e(0)\|_{H_{\varepsilon}} \\ &\le \|P_{\varepsilon} U^0 - u_{\varepsilon}^0\|_{H_{\varepsilon}} + \|P_{\varepsilon} U^e(0) - u_{\varepsilon}^e(0)\|_{H_{\varepsilon}}. \end{aligned} $$

The first term on the right hand side of the inequality goes to zero by the additional assumption on $u_{\varepsilon}^0$, while the second term goes to zero, as previously shown.

In regard to the second members $q_{\varepsilon} = \dot{u}_{\varepsilon}^e$ and $Q = \dot{U}^e$, we again use Lemma 1 by proving $e(\varepsilon, \dot{u}_{\varepsilon}^e) \to \hat{e}(\dot{u}^e) + \partial_3 \dot{u}^{e1} \otimes_s e_3$ in $L^2(\Omega; \mathbb{S}^3)$ uniformly on $[0, T]$ because of the smoothness of $g$ with respect to the time $t$.

In the matter of resolvants $(I - A_{\varepsilon})^{-1}, (I - A)^{-1}$, we have to show that

$$ \forall Y \in X, \text{ dense in } H, \|P_{\varepsilon}(I - A)^{-1}Y - (I - A_{\varepsilon})^{-1}P_{\varepsilon}Y\|_{H_{\varepsilon}} \to 0. $$

We set $z_{\varepsilon} = (I - A_{\varepsilon})^{-1}P_{\varepsilon}Y$ and $Z = (I - A)^{-1}Y$ that satisfy the problems

$$ \left\{ \begin{array}{l} z_{\varepsilon} \in H_{\varepsilon}; \\ \int_{\Omega} (a + b)e(\varepsilon, z_{\varepsilon}) \cdot e(\varepsilon, v) dx = \int_{\Omega} be(\varepsilon, P_{\varepsilon}Y) \cdot e(\varepsilon, v) dx \quad \forall v \in H_{\varepsilon}, \end{array} \right. $$

and

$$ \left\{ \begin{array}{l} Z = (z, z^1) \in H; \\ \int_{\Omega} (a + b)(\hat{e}(z) + \partial_3 z^1 \otimes_s e_3) \cdot (\hat{e}(v) + \partial_3 v^1 \otimes_s e_3) dx \\ = \int_{\Omega} b(\hat{e}(y) + \partial_3 y^1 \otimes_s e_3) \cdot (\hat{e}(v) + \partial_3 v^1 \otimes_s e_3) dx \quad \forall V = (v, v^1) \in H. \end{array} \right. $$

Now we consider $\|P_{\varepsilon}Z - z_{\varepsilon}\|_{H_{\varepsilon}}$. They involve a kind of static problem to which the previous two-fields result on reduction of dimension applies. The proof is then complete and we now have the main theorem.

**Theorem 2** Under assumptions (H$_1$)-(H$_3$), when $\varepsilon$ goes to zero, the sequence of solutions $u_{\varepsilon}$ in $H_{\Gamma_D}^1(\Omega; \mathbb{R}^3)$ to problem $(P_{\varepsilon})$ converges uniformly on $[0, T]$ in the sense of Trotter to the solution $U = (u, u^1)$ in $V_{KL} \times H^1(-1, 1; L^2(\omega; \mathbb{R}^3))/L^2(\omega; \mathbb{R}^3)$ to

$$ (\bar{P}): \left\{ \begin{array}{l} U = (u, u^1) \in C^{1,1}([0, T]; V_{KL} \times H^1(-1, 1; L^2(\omega; \mathbb{R}^3))/L^2(\omega; \mathbb{R}^3)); \\ \int_{\Omega} [a(\hat{e}(u) + \partial_3 u^1 \otimes_s e_3) + b(\hat{e}(\dot{u}) + \partial_3 \dot{u}^1 \otimes_s e_3)] \cdot (\hat{e}(v) + \partial_3 v^1 \otimes_s e_3) dx \\ = \int_{\Omega} f \cdot v dx + \int_{\Gamma_N} g \cdot v dH_2 \\ \forall V = (v, v^1) \in V_{KL} \times H^1(-1, 1; L^2(\omega; \mathbb{R}^3))/L^2(\omega; \mathbb{R}^3). \end{array} \right. $$
Hence *U* solves a problem of visco-elasticity with short memory of Kelvin-Voigt type, but involving a *couple* (u, u¹) of state variables. Clearly the field u¹ can be eliminated as in [1] so that *u* does solve a problem of visco-elasticity with long (but fading) memory.

**Acknowledgements**

The first author was supported by DPST scholarship (Development and Promotion of Science and Technology Talents Project) for doing this research.

**Funding**

The first author was supported by DPST scholarship (Development and Promotion of Science and Technology Talents Project) for doing this research.

**Competing interests**

The authors declare that they have no competing interests.

**Authors' contributions**

CL conceived of the research work and verified the result. SO and YT proved the main result, interpreted the result and wrote a paper. All authors read and approved the final manuscript.

**Authors' information**

Yotsawat Terapabkajornded, email: ongyosawat@gmail.com; Somsak Orankitjaroen, email: somsak.ora@mahidol.ac.th; Christian Licht, email: christian.licht@umontpellier.fr.

**Author details**

¹Department of Mathematics, Faculty of Science, Mahidol University, Bangkok, Thailand. ²LMGC-UMR 5508, Université de Montpellier-CC048, Montpellier, France. ³Centre of Excellence in Mathematics, CHE, Bangkok, Thailand.

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Received: 29 January 2019 Accepted: 14 April 2019 Published online: 14 May 2019

**References**

1. Licht, C.: Thin linearly visco-elastic Kelvin-Voigt plates. C. R., Méc. **341**, 697–700 (2013)

2. Licht, C., Weller, T.: Approximation of semi-groups in the sense of Trotter and asymptotic mathematical modeling in physics of continuous media. Discrete Contin. Dyn. Syst., Ser. S **12**, 1709–1741 (2019)

3. Ciarlet, P.G.: Mathematical elasticity, vol. II. In: Lions, J.L., Papanicolaou, G., Fujita, H., Keller, H.B. (eds.) Study in Mathematics and Its Applications, vol. 27, pp. 14–80. North-Holland, Amsterdam (1997)

**Submit your manuscript to a SpringerOpen**
**journal and benefit from:**

► Convenient online submission

► Rigorous peer review

► Open access: articles freely available online

► High visibility within the field

► Retaining the copyright to your article

Submit your next manuscript at ► springeropen.com