%~Mouliné par MaN_auto v.0.40.4 (9e9f6e09) 2026-06-09 18:00:51
\documentclass[CRMECA,Unicode,XML,screen,biblatex]{cedram}
\addbibresource{CRMECA_WELLER_20260252.bib}
\COI{The authors do not work for, advise, own shares in, or receive funds from any organization that could benefit from this article, and have declared no affiliations other than their research organizations.}

\usepackage{float}
\usepackage{mathtools}
\usepackage{scalerel}
\usepackage{stackengine}
\usepackage{multirow}
\usepackage[mathscr]{euscript}
\usepackage{tensor}





%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%USAGE ?
\DeclareMathAlphabet{\mathcalboondox}{U}{BOONDOX-calo}{m}{n}
%\DeclareMathAlphabet{\mathbcalboondox}{U}{BOONDOX-calo}{b}{n}
\allowdisplaybreaks
\DeclareFontFamily{U}{BOONDOX-calo}{\skewchar\font=45 }
\DeclareFontShape{U}{BOONDOX-calo}{m}{n}{ <-> s*[1.05] BOONDOX-r-calo}{}
\DeclareFontShape{U}{BOONDOX-calo}{b}{n}{ <-> s*[1.05] BOONDOX-b-calo}{}

\SetMathAlphabet{\mathcalboondox}{bold}{U}{BOONDOX-calo}{b}{n}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


\newcommand{\q}[1]{``#1''}



\stackMath
\newcommand\reallywidehat[1]{
\savestack{\tmpbox}{\stretchto{
\scaleto{
\scalerel*[\widthof{\ensuremath{#1}}]{\kern-.6pt\bigwedge\kern-.6pt} {\rule[-\textheight/2]{1ex}{\textheight}} }{\textheight} }{0.5ex}}
\stackon[1pt]{#1}{\tmpbox} }
\parskip 1ex


\DeclareFontFamily{U}{mathx}{}
\DeclareFontShape{U}{mathx}{m}{n}{<-> mathx10}{}
\DeclareSymbolFont{mathx}{U}{mathx}{m}{n}
\DeclareMathAccent{\widecheck}{0}{mathx}{"71}

\let\oldcheck\check
\renewcommand*{\check}[1]{\mathchoice{\widecheck#1}{\widecheck{#1}}{\widecheck{#1}}{\oldcheck{#1}}}
%\makeatletter
%%\newcommand{\leqnomode}{\tagsleft@true}
%%\newcommand{\reqnomode}{\tagsleft@false}
%\makeatother





\newcommand{\otS}{\otimes_\textup{\!\tiny{s}}}

\newcommand{\wht}{\widehat}
%\newcommand{\rwht}{\reallywidehat}
\newcommand{\ovl}{\bar}

\newcommand{\clbdW}{\mathcalboondox{W}}
\newcommand{\clbdJ}{\mathcalboondox{J}}

\newcommand{\clG}{\mathcal{G}}
\newcommand{\clT}{\mathcal{T}}
\newcommand{\clV}{\mathcal{V}}
\newcommand{\clS}{\mathcal{S}}
\newcommand{\clU}{\mathcal{U}}
\newcommand{\clR}{\mathcal{R}}

\newcommand{\mSe}{{\clS_\eps}}

\newcommand{\scH}{\mathscr{H}}
\newcommand{\eps}{\varepsilon}
\newcommand{\Omg}{\Omega}
\newcommand{\omg}{\omega}
%\newcommand{\gma}{\gamma}


\newcommand{\mbR}{\mathbb{R}}
\newcommand{\mbS}{\mathbb{S}}

\newcommand{\fkq}{\mathfrak{q}}
\newcommand{\fkc}{\mathfrak{c}}
\newcommand{\fkp}{\mathfrak{p}}
\newcommand{\fkr}{\mathfrak{r}}
\newcommand{\fki}{\mathfrak{i}}
\newcommand{\fkS}{\mathfrak{S}}
\newcommand{\fkDeux}{\mathfrak{2}}


\newcommand{\os}{\bar{s}}
\newcommand{\ou}{\bar{u}}

\newcommand{\tD}{\textup{D}}
\newcommand{\tT}{\textup{T}}
\newcommand{\tN}{\textup{N}}
\newcommand{\tS}{\textup{S}}
\newcommand{\tA}{\textup{A}}
\newcommand{\tP}{\textup{P}}
\newcommand{\tQ}{\textup{Q}}
\newcommand{\tH}{\textup{H}}
\newcommand{\tF}{\textup{F}}

\newcommand{\lims}{\lim_{s \to \os}}
\newcommand{\liminfs}{\liminf_{s \to \os}}
\newcommand{\limsups}{\limsup_{s \to \os}}
\DeclareMathOperator{\Lin}{Lin}
\DeclareMathOperator{\Min}{Min}
\DeclareMathOperator{\supp}{supp}
\DeclareMathOperator{\BN}{BN}
\DeclareMathOperator{\beam}{beam}
\DeclareMathOperator{\Argmin}{Argmin}
\DeclareMathOperator{\dist}{dist}

\DeclareMathOperator{\eff}{eff}
\DeclareMathOperator{\sgn}{sgn}
\DeclareMathOperator{\upp}{up}
\DeclareMathOperator{\low}{low}
\DeclareMathOperator{\sym}{sym}
\DeclareMathOperator{\Minn}{Min}
\DeclareMathOperator{\Max}{Max}
\DeclareMathOperator{\plate}{plate}

\DeclareMathOperator{\itd}{\mathit{d}}%\,d

\newcommand{\muk}{^{k}\!\bar{\mu}}
\newcommand{\muo}{^{1}\!\!\bar{\mu}}
\newcommand{\mut}{^{2}\!\bar{\mu}}
\newcommand{\mur}{^{\fkr}\!\bar{\mu}}
\newcommand{\mui}{^{\fki}\!\bar{\mu}}

\newcommand{\Wqrb}{\clbdW_{^{\fkq, \fkr}\tilde{b}}}
\newcommand{\Wqib}{\clbdW_{^{\fkq, \fki}\tilde{b}}}


\newcommand{\ts}[4]{\tensor*[^{#1}]{#2}{^{#3}_{#4}}}

\newcommand{\pus}{\ts{\fkc}{u}{s}{}}
\newcommand{\pPs}{\ts{\fkc}{\tP}{s}{}}

\newcommand{\pvxs}{\ts{\fkc}{v}{}{s}}
\newcommand{\puxs}{\ts{\fkc}{u}{}{s}}
\newcommand{\pPxs}{\ts{\fkc}{\tP}{}{s}}
\newcommand{\pQxs}{\ts{\fkq}{\tQ}{}{s}}
\newcommand{\pTxh}{\ts{\fkc}{\clT}{}{h}}
\newcommand{\pVxh}{\ts{\fkc}{\clV}{}{h}}
\newcommand{\ptauxh}{\ts{\fkp}{\tau}{}{h}}

\newcommand{\pfeh}{\ts{\fkp}{f}{\eps}{h}}
\newcommand{\pgeh}{\ts{\fkc}{g}{\eps}{h}}
\newcommand{\paeh}{\ts{\fkp\!\!}{a}{\eps}{h}}
\newcommand{\pVeh}{\ts{\fkc}{\clV}{\eps}{h}}


\newcommand{\pIGmaehN}{\ts{\fkc}{\Gamma}{\eps}{h\tN}}
\newcommand{\pIGmaxD}{\ts{\fkq}{\Gamma}{}{\tD}}
\newcommand{\pIGmaxN}{\ts{\fkq}{\Gamma}{}{\tN}}
\newcommand{\pGmaehD}{\ts{\fkc}{\Gamma}{\eps}{h\tD}}

\newcommand{\pUpD}{\ts{\fkq}{\clU}{+}{\tD}}
\newcommand{\pUmD}{\ts{\fkq}{\clU}{-}{\tD}}
\newcommand{\pUzpm}{\ts{\fkq}{\clU}{0\pm}{}}
\newcommand{\pUonep}{\ts{\fkq}{\clU}{1+}{}}
\newcommand{\pUonem}{\ts{\fkq}{\clU}{1-}{}}
\newcommand{\pUonepm}{\ts{\fkq}{\clU}{1\pm}{}}
\newcommand{\pUtwop}{\ts{\fkq}{\clU}{2+}{}}
\newcommand{\pUtwom}{\ts{\fkq}{\clU}{2-}{}}
\newcommand{\pUtwopm}{\ts{\fkq}{\clU}{2\pm}{}}
\newcommand{\pU}{\ts{\fkq}{\clU}{}{}}
\newcommand{\pUi}{\ts{\fkq}{\clU}{i}{}}
\newcommand{\pUzero}{\ts{\fkq}{\clU}{0}{}}
\newcommand{\pUone}{\ts{\fkq}{\clU}{1}{}}
\newcommand{\pUtwo}{\ts{\fkq}{\clU}{2}{}}

\newcommand{\pRp}{\ts{\fkq}{\clR}{+}{}}
\newcommand{\pRm}{\ts{\fkq}{\clR}{-}{}}
\newcommand{\pRpm}{\ts{\fkq}{\clR}{\pm}{}}

\newcommand{\pIOxh}{\ts{\fkp}{\Omg}{}{h}}
\newcommand{\pIOeh}{\ts{\fkp}{\Omg}{\eps}{h}}
\newcommand{\pIGmaxh}{\ts{\fkp}{G}{}{h}}
\newcommand{\pIGeh}{\ts{\fkp}{G}{\eps}{h}}
\newcommand{\pIGxh}{\ts{\fkp}{G}{}{h}}

\newcommand{\pqV}{\ts{\fkq,\fkr}{\clV}{}{}}
\newcommand{\pqJ}{\ts{\fkq,\fkr}{\clbdJ}{}{}}
\newcommand{\pqvxs}{\ts{\fkc,\fkr}{v}{}{s}}
\newcommand{\twotauxh}{\ts{2}{\tau}{}{h}}


\stackMath
\newcommand\rwh[1]{
\savestack{\tmpbox}{\stretchto{
\scaleto{
\scalerel*[\widthof{\ensuremath{#1}}]{\kern-.6pt\bigwedge\kern-.6pt} {\rule[-\textheight/2]{1ex}{\textheight}} }{\textheight} }{0.5ex}}
\stackon[1pt]{#1}{\tmpbox} }


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Better widebar

\makeatletter
\let\save@mathaccent\mathaccent
\newcommand*\if@single[3]{%
  \setbox0\hbox{${\mathaccent"0362{#1}}^H$}%
  \setbox2\hbox{${\mathaccent"0362{\kern0pt#1}}^H$}%
  \ifdim\ht0=\ht2 #3\else #2\fi
  }
%The bar will be moved to the right by a half of \macc@kerna, which is computed by amsmath:
\newcommand*\rel@kern[1]{\kern#1\dimexpr\macc@kerna}
%If there's a superscript following the bar, then no negative kern may follow the bar;
%an additional {} makes sure that the superscript is high enough in this case:
\newcommand*\widebar[1]{\@ifnextchar^{{\wide@bar{#1}{0}}}{\wide@bar{#1}{1}}}
%Use a separate algorithm for single symbols:
\newcommand*\wide@bar[2]{\if@single{#1}{\wide@bar@{#1}{#2}{1}}{\wide@bar@{#1}{#2}{2}}}
\newcommand*\wide@bar@[3]{%
  \begingroup
  \def\mathaccent##1##2{%
%Enable nesting of accents:
    \let\mathaccent\save@mathaccent
%If there's more than a single symbol, use the first character instead (see below):
    \if#32 \let\macc@nucleus\first@char \fi
%Determine the italic correction:
    \setbox\z@\hbox{$\macc@style{\macc@nucleus}_{}$}%
    \setbox\tw@\hbox{$\macc@style{\macc@nucleus}{}_{}$}%
    \dimen@\wd\tw@
    \advance\dimen@-\wd\z@
%Now \dimen@ is the italic correction of the symbol.
    \divide\dimen@ 3
    \@tempdima\wd\tw@
    \advance\@tempdima-\scriptspace
%Now \@tempdima is the width of the symbol.
    \divide\@tempdima 10
    \advance\dimen@-\@tempdima
%Now \dimen@ = (italic correction / 3) - (Breite / 10)
    \ifdim\dimen@>\z@ \dimen@0pt\fi
%The bar will be shortened in the case \dimen@<0 !
    \rel@kern{0.6}\kern-\dimen@
    \if#31
      \overline{\rel@kern{-0.6}\kern\dimen@\macc@nucleus\rel@kern{0.4}\kern\dimen@}%
      \advance\dimen@0.4\dimexpr\macc@kerna
%Place the combined final kern (-\dimen@) if it is >0 or if a superscript follows:
      \let\final@kern#2%
      \ifdim\dimen@<\z@ \let\final@kern1\fi
      \if\final@kern1 \kern-\dimen@\fi
    \else
      \overline{\rel@kern{-0.6}\kern\dimen@#1}%
    \fi
  }%
  \macc@depth\@ne
  \let\math@bgroup\@empty \let\math@egroup\macc@set@skewchar
  \mathsurround\z@ \frozen@everymath{\mathgroup\macc@group\relax}%
  \macc@set@skewchar\relax
  \let\mathaccentV\macc@nested@a
%The following initialises \macc@kerna and calls \mathaccent:
  \if#31
    \macc@nested@a\relax111{#1}%
  \else
%If the argument consists of more than one symbol, and if the first token is
%a letter, use that letter for the computations:
    \def\gobble@till@marker##1\endmarker{}%
    \futurelet\first@char\gobble@till@marker#1\endmarker
    \ifcat\noexpand\first@char A\else
      \def\first@char{}%
    \fi
    \macc@nested@a\relax111{\first@char}%
  \fi
  \endgroup
}
\makeatother

\let\oldbar\bar
\renewcommand*{\bar}[1]{{\mathchoice{\widebar{#1}}{\widebar{#1}}{\widebar{#1}}{\oldbar{#1}}}}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\graphicspath{{./figures/}}

\newcommand*{\mk}{\mkern -1mu}
\newcommand*{\Mk}{\mkern -2mu}
\newcommand*{\mK}{\mkern 1mu}
\newcommand*{\MK}{\mkern 2mu}

\hypersetup{urlcolor=purple, linkcolor=blue, citecolor=red}

\newcommand*{\relabel}{\renewcommand{\labelenumi}{(\theenumi)}}
\newcommand*{\relabell}{\renewcommand{\labelenumii}{(\theenumii)}}
\newcommand*{\romanenumi}{\renewcommand*{\theenumi}{\roman{enumi}}\relabel}
\newcommand*{\romanenumii}{\renewcommand*{\theenumii}{\roman{enumii}}\relabell}
\newcommand*{\Romanenumi}{\renewcommand*{\theenumi}{\Roman{enumi}}\relabel}
\newcommand*{\alphenumi}{\renewcommand*{\theenumi}{\alph{enumi}}\relabel}
\newcommand*{\Alphenumi}{\renewcommand*{\theenumi}{\Alph{enumi}}\relabel}
\let\oldtilde\tilde
\renewcommand*{\tilde}[1]{\mathchoice{\widetilde{#1}}{\widetilde{#1}}{\oldtilde{#1}}{\oldtilde{#1}}}
\let\oldhat\hat
\renewcommand*{\hat}[1]{\mathchoice{\widehat{#1}}{\widehat{#1}}{\oldhat{#1}}{\oldhat{#1}}}
\let\oldexists\exists
\renewcommand*{\exists}{\mathrel{\oldexists}}
\let\oldforall\forall
\renewcommand*{\forall}{\mathrel{\oldforall}}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\renewcommand*{\to}{\mathchoice{\longrightarrow}{\rightarrow}{\rightarrow}{\rightarrow}}


\let\oldmapsto\mapsto
\renewcommand*{\mapsto}{\mathchoice{\longmapsto}{\oldmapsto}{\oldmapsto}{\oldmapsto}}

%~ This fixes spacing issues around \left( and \right) and other similar symbols
\let\originalleft\left 
\let\originalright\right
\renewcommand{\left}{\mathopen{}\mathclose\bgroup\originalleft}
\renewcommand{\right}{\aftergroup\egroup\originalright}

\newenvironment{alcases}{\left\{\begingroup\arraycolsep=0pt\renewcommand{\arraystretch}{1.2}\begin{array}{r@{\;}l@{\quad}l}}{\end{array}\endgroup\right.}

 \newenvironment{alrcases}{\left.\begingroup\arraycolsep=0pt\renewcommand{\arraystretch}{1.2}\begin{array}{r@{\;}l@{\;}l}}{\end{array}\endgroup\right\}}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\title{An assembly of two slender anisotropic beams through a very thin adhesive layer}
\alttitle{Assemblage de deux poutres élancées anisotropes via une couche adhésive très mince}

\keywords{
\kwd{assemblies of slender anisotropic beams}
\kwd{asymptotic analysis}
\kwd{variational convergence}
\kwd{mechanical constraints}
}
\altkeywords{
\kwd{assemblages de poutres minces anisotropes}
\kwd{analyse asymptotique}
\kwd{convergence variationnelle}
\kwd{liaisons mécaniques}
}


\author{\firstname{Christian} \lastname{Licht}}
\address{Department of Mathematics, Faculty of Science, Mahidol University, Bangkok,Thailand}
\address{Centre of Excellence in Mathematics, CHE, Bangkok, Thailand}
\address{LMGC, Université de Montpellier, CNRS, Montpellier, France}
\email[C. Licht]{christian.licht@umontpellier.fr}



\author{\firstname{Somsak} \lastname{Orankitjaroen}}
\addressSameAs{1}{Department of Mathematics, Faculty of Science, Mahidol University, Bangkok,Thailand}
\addressSameAs{2}{Centre of Excellence in Mathematics, CHE, Bangkok, Thailand}
\email[S. Orankitjaroen]{somsak.ora@mahidol.ac.th}


\author{\firstname{Thibaut} \lastname{Weller}\IsCorresp}
\addressSameAs{3}{LMGC, Université de Montpellier, CNRS, Montpellier, France}
\email[T. Weller]{thibaut.weller@umontpellier.fr}




\begin{abstract}
We derive various models of assemblies of slender anisotropic linearly elastic beams through an asymptotic analysis taking into account a triplet of small parameters associated with the slenderness of the beams but also the thinness and the stiffness of a very thin third body which connects them. Our models allow the description of the mechanical constraint of the linkage between two beams which strongly depends on the relative orders of magnitude of the previous parameters.
\end{abstract}


\begin{altabstract}
Nous proposons divers modèles d'assemblages de poutres élancées anisotropes linéairement élastiques par l'intermédiaire d'une analyse asymptotique qui prend en compte un triplet de petits paramètres liés à l'élancement des poutres mais aussi à l'épaisseur et à la rigidité d'un troisième corps très mince assurant leur connexion. Nos modèles permettent de décrire le comportement mécanique de la liaison entre les deux poutres, qui dépend fortement des ordres de grandeur relatifs de ces paramètres.
\end{altabstract}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\begin{document}
\maketitle

\section{Setting the problem}



Using the same asymptotic analysis by variational convergence, we continue the studies~\cite{Licht1997,Licht2007} by considering the case of several kinds of assemblies of two slender anisotropic beams through a very thin layer. As usual we do not distinguish between $\mbR^3$ and the Euclidean physical space whose basis is $\{e_1, e_2, e_3 \}$. Throughout the paper, lower Greek (resp. Latin) indices run from $1$ to $2$ (resp. $1$ to $3$) and for every $\xi = (\xi_1, \xi_2, \xi_3)$ of $\mbR^3$, $\wht{\xi}$ denotes $(\xi_1, \xi_2)$.
The space of symmetric $3 \times 3$ matrices is denoted by $\mbS^3$ and we write $\otS$ for the symmetrized tensor product, i.e. given $\xi$ and $\zeta$ in $\mbR^3$, the element $\xi \otS \zeta$ of $\mbS^3$ is defined componentwise by $(\xi \otS \zeta)_{ij} := \frac{1}{2} (\xi_i \zeta_j + \xi_j \zeta_i )$. Finally $\Lin(\mbS^3)$ stands for the space of linear mappings from $\mbS^3$ into itself. Seeking to address several geometric and mechanical situations in a unified manner, we opted for the introduction of upper left superscripts. This makes it possible to consider various cases indexed by $\fkc = (\fkp, \fkq)$ in~$\{1,2\}^2$.

Given two positive real numbers $L^\pm$ and a domain $\omg$ of $\mbR^2$ with a Lipschitz-continuous boundary $\partial \omg$ such that
\[
\int_\omg x_\alpha \itd \wht{x} = \int_\omg x_1x_2 \itd \wht{x} = 0,
\]
for all pairs of numbers $(\eps, h)$ in $(0, +\infty) \times (0, \Minn(L^+, L^-)$), let
\begin{equation}\label{eq1_}
\begin{alcases}
\Omg^{\eps} &:= \eps\omg \times (-L^-, L^+),  &\;\;\:\Omg^{\eps\pm} := \Omg^\eps \cap \left\{x^\eps \in \mbR^3;\, \pm x^\eps_3 > 0\right \} \\
B^\eps_h &:= \eps\omg \times (-h, h) &\\
\pIOeh &:= \ts{\fkp}{\Omg}{\eps+}{h} \cup \ts{\fkp}{\Omg}{\eps-}{h} &\\
\ts{1}{\Omg}{\eps\pm}{h} &:= \Omg^{\eps\pm} \cap \left\{\pm x^\eps_3 > h\right\}, \quad &\ts{2}{\Omg}{\eps\pm}{h} := \Omg^{\eps\pm} \pm he_3 \\
S^\eps &:= \eps\omg \times \{0\},  &\quad\; S^\eps_\eta := \eps \omg \times \{\eta\} \quad \forall \eta \in \mbR\setminus\{0\} \\
\ts{\fkp}{G}{\eps}{h} &:= \eps\omg \times \left\{\pm x_3 \in \left[0, L^\pm + (\fkp-1)h\right)\right\}. &
\end{alcases}
\end{equation}
The domains $\ts{\fkp}{\Omg}{\eps\pm}{h}$ are occupied by two linearly elastic slender beams which are connected by a very thin linearly elastic layer occupying $B^\eps_h$. The case $\fkp = 1$ corresponds to the analysis of the influence of the stiffness and thickness of an inclusion in a slender structure, whereas the case $\fkp = 2$ corresponds to the analysis of the influence of the stiffness and the thickness of an adhesive joint on the bonding of two given slender beams.

We assume that the ``upper part'' $S^\eps_{L^+ + (\fkp-1)h} =: {}^{\fkp}\!S^\eps_{\upp} $ is clamped on a rigid body while the ``lower part'' $S^\eps_{-(L^- + (\fkp-1)h)} =: \,^{\fkp}\!S^\eps_{\low}$ is either clamped on a rigid body $(\fkq=1)$ or subjected to surface forces $(\fkq=2)$! Moreover the beams are subjected to body forces in $\pIOeh$ and surface forces on $\eps \partial \omg \times \{h < \pm x^\eps_3 < L^\pm + (\fkp-1)h\}$. More precisely, if
\begin{equation}
\ts{(\fkp,1)}{\Gamma}{\eps}{h\tD} := {}^{\fkp}\!S^\eps_{\upp} \cup {}^{\fkp}\!S^\eps_{\low}, \quad
\ts{(\fkp,2)}{\Gamma}{\eps}{h\tD} := {}^{\fkp}\!S^\eps_{\upp},
\end{equation}
we assume that the assembly is clamped along $\pGmaehD$, subjected to body forces of density $\pfeh$ in $L^2(\pIOeh, \mbR^3)$ and surface forces of density $\pgeh$ in $L^2(\pIGmaehN, \mbR^3)$ with $\pIGmaehN := \partial \, \pIOeh \setminus (\pGmaehD \cup S^\eps_{\pm h})$ whereas the thin layer is free of loading and perfectly stuck to the beams along $S^\eps_{\pm h}$. The following Figure~\ref{F:1} summarizes the above information:

\begin{figure}[!htbp]
\centering
\includegraphics[scale=.28]{schema_poutres_vectorise}
\caption{Four configurations of abutting slender beams with diameter of order $\varepsilon$. The left panels correspond to slender beams containing a very small inclusion, while the right panels represent a junction through a very small adhesive layer. The traction-free boundary is shown in red; its size is discussed in assumption~\eqref{eq:q1} below.}\label{F:1}
\end{figure}


Lastly we denote the elasticity tensors of the beams and of the layer by $\paeh$, $\mu b$ respectively. They satisfy
\begin{equation}\label{3f}
\begin{cases}
\bullet\; \mu > 0 \\
\bullet\; \left(\paeh, b\right) \in L^\infty
\left(\pIOeh, \Lin\left(\mbS^3\right)\right) \times \Lin\left(\mbS^3\right)  \\
\quad\text{with }\; \paeh, b \text{ symmetric }
\; \exists \, \alpha > 0 \;\text{ s.t.} \\
\qquad
\begin{alrcases}
2 \clbdW_{\paeh}(\xi) &:= \paeh(x^\eps) \xi \cdot \xi \ge \alpha |\xi|^2_{\mbS^3} \quad \text{a.e. }\; x^\eps \in \pIOeh \\
2\clbdW_b(\xi) &:= b\xi \cdot \xi \ge \alpha |\xi|^2_{\mbS^3}
\end{alrcases}
\forall \xi \in \mbS^3.
\end{cases}
\end{equation}

In the following for any domain $\clG$ of $\mbR^d, 1 \le d \le 3$, $H^1_\gamma(\clG, \mbR^{d'}), 1 \le d' \le 3,$ stands for the subspace of $H^1(\clG, \mbR^{d'})$ of the elements with vanishing trace on a subset $\gamma$ of the boundary $\partial \clG$. Hence the problem of the equilibrium of the assembly involves a triplet $s := (\eps, h, \mu)$ of data and can be formulated as:
\begin{multline}\label{eq:CPs}
%\tag{$\pPs$}
\tag{$\,\!^\mathfrak{c}\textup{P}^{s}$}
\Minn
\left\{ \int_{\pIOeh} \clbdW_{\paeh} (e^\eps(v)) \itd x^\eps + \mu \int_{B^\eps_h} \clbdW_b(e^\eps(v)) \itd x^\eps - \int_{\pIOeh} \pfeh \cdot v \itd x^\eps \right. \\
\left.- \int_{\pIGmaehN} \pgeh \cdot v \itd \scH^\eps_2; \quad v \in \ts{\fkc}{\clV}{\eps}{h} := H^1_{\pGmaehD} \left(\pIGeh, \mbR^3\right)
\right\}
\end{multline}
where $e^\eps_{ij}(v) = (\partial v_i/\partial x^\eps_j + \partial v_j/\partial x^\eps_i)/2$ is the strain tensor associated with the displacement $v$ and $\scH^\eps_2$ is the $2$-dimensional Hausdorff-measure. Obviously~\eqref{eq:CPs} has a unique solution $\pus$.

\section{Asymptotic behavior}

As is customary in the analysis of slender beams (see, e.g., \cite{Murat1999, trabu}), we introduce:
\begin{itemize}
\item a change of coordinates $\pi^\eps$:
\begin{equation}
x \in \mbR^3 \mapsto x^\eps = \pi^\eps x := (\eps \wht{x}, x_3) \in \mbR^3
\end{equation}
and drop the index $\eps$ for the image by $(\pi^\eps)^{-1}$ of the previous domains while, in the sequel, the space variables $x^\eps$ in $\ts{\fkp}{G}{\eps}{h}$ and $x$ in $\ts{\fkp}{G}{}{h}$ are systematically connected by $x^\eps = \pi^\eps x$.
\item An operator $\mSe$ which transforms any element $w$ of $\pVeh$ into an element of $\pVxh:=H^1_{ \ts{\fkc}{\Gamma}{}{h\tD} } (\pIGmaxh, \mbR^3)$ such that:
\begin{equation}
\wht{\mSe w}(x) = \wht{w}(x^\eps), \quad (\mSe w)_3(x) = \frac{1}{\eps} w_3(x^\eps)
\end{equation}
and
\begin{equation}
\begin{array}{r@{}l@{}r}
e^\eps(w)(x^\eps) \; =& \; \eps \; e^{\beam}(\eps, \mSe \, w)(x), & \\[2mm] e^{\beam}(\eps, z) \; :=& \;
\begin{bmatrix}
e_{\alpha \beta}(z) / \varepsilon^2 & e_{\alpha 3}(z) / \varepsilon \\
\sym & e_{33}(z)
\end{bmatrix} & \multicolumn{1}{c}{\multirow{2}{*}{ $\raisebox{0pt}[0cm][0cm]{$\left.\rule{0pt}{0.9cm}\right\}$}
\quad \forall z \in H^1 \left({ }^{\fkp}G_h, \mbR^3\right )$. }}
\\[3mm] e_{ij}(z) \; :=& \; \frac{1}{2} \left(\partial_i z_j + \partial_j z_i\right) \;, \quad 1 \le i, j \le 3 &
\end{array}
\end{equation}

\item Some assumptions on the loading and the stiffness of the beams. To this end we introduce:
\[
\ts{1}{\Gamma}{}{\tD} := S_{L^+} \cup S_{L^-}, \quad
\ts{2}{\Gamma}{}{\tD} := S_{L^+}, \quad
\ts{\fkq}{\Gamma}{}{\tN} := \partial \, \Omega \setminus \ts{\fkq}{\Gamma}{}{\tD}, \quad
\ts{\fkc}{\Gamma}{}{h\tN} := \partial \, \ts{\fkp}{G}{}{h} \setminus \ts{\fkc}{\Gamma}{}{h\tD}
\]
and, for any function $\varphi$ defined on a subset of $\ovl{\Omg}$, the mapping $\ts{\fkp}{\tau}{}{h} \, \varphi$ such that:
\begin{equation}
\ts{\fkp}{\tau}{}{h} \, \varphi(\wht{x}, x_3) = \varphi \bigl(\wht{x}, x_3 - \sgn (x_3)(\fkp-1)h \bigr) \quad \forall x \in \pIOxh
\end{equation}
where $\sgn(x_3)$ denotes the sign of $x_3$. We assume that:
%revoir cela.
\begin{equation}\label{eq:q1}
\tag{$\tH_1$}
\begin{dcases}
\text{There exist $h_0 > 0$, $f$ in $L^2(\Omg, \mbR^3)$, $g$ in $L^2(\pIGmaxN, \mbR^3)$ with $(\supp \,(f), \supp \, (g)) \cap B_{h_0} = \emptyset
$ } \\
\text{and $a$ symmetric in $L^\infty(\Omg, \Lin(\mbS^3))$ such that} \\
\exists \, \alpha > 0;\quad 2\clbdW_a(\xi) := a(x) \xi \cdot \xi \ge \alpha |\xi|^2 \quad \text{a.e. } \; x \text{ in }\; \Omg,\quad \forall \xi \in \mbS^3 \\
\bullet\; \wht{\pfeh}(x^\eps) = \eps^2\; \ts{\fkp}{\tau}{}{h} \wht{f}(x),\quad \ts{\fkp}{f}{\eps}{h3}(x^\eps) = \eps\; \ts{\fkp}{\tau}{}{h} f_3(x) \quad \forall x \in \pIOxh \\
\bullet\; \wht{\pgeh}(x^\eps) = \eps^3\; \ts{\fkp}{\tau}{}{h} \wht{g}(x),\quad \ts{\fkc}{g}{\eps}{h3}(x^\eps) = \eps^2\; \ts{\fkp}{\tau}{}{h} g_3(x)
\;\;\: \forall x \in \ts{\fkc}{\Gamma}{}{h\tN} \cap \left\{\pm x_3 \in \left[0, L^{\pm} + (\fkp-1)h\right)\right\} \\
\bullet\; \wht{\pgeh}(x^\eps) = \eps^2\; \ts{\fkp}{\tau}{}{h} \wht{g}(x),\quad \ts{\fkc}{g}{\eps}{h3}(x^\eps) = \eps\; \ts{\fkp}{\tau}{}{h} g_3(x)
\quad \forall x \in \ts{\fkc}{\Gamma}{}{h\tN} \cap \left\{x_3 = -\left(L^{-} + (\fkp-1)h\right) \right\} \\
\bullet\; \paeh(x^\eps) = \ts{\fkp}{\tau}{}{h} a(x) \quad \forall x \in \pIOxh \\
\end{dcases}
\end{equation}
\end{itemize}
Hence $\puxs := \mSe \pus$ is the unique solution of
\begin{multline}\label{eq:cPs}
\tag{$\ts{\fkc}{\tP}{}{s}$}
\Minn
\left\{\int_{\pIOxh} \clbdW_{\ptauxh a} \left(e^{\beam}(\eps, v)\right) \itd x + \mu\int_{B_h} \clbdW_b\left(e^{\beam}(\eps, v)\right) \itd x - \int_{\pIOxh} \ts{\fkp}{\tau}{}{h} \, f \cdot v \itd x \right. \\
\left. - \int_{\ts{\fkc}{\Gamma}{}{h\tN}} \ts{\fkp}{\tau}{}{h} \, g \cdot v \itd \scH_2; \quad v \in \pVxh := H^1_{\ts{\fkc}{\Gamma}{}{h\tD}} \left(\pIGxh, \mbR^3\right)
\;\right\}.
\end{multline}
We will consider $s = (\eps, h, \mu)$ as a triplet of \emph{parameters} taking values in a countable set $\fkS$ of $(0, +\infty) \times (0, h_0) \times (0, +\infty) $ with a unique cluster point $\bar{s}$ in $\{0\} \times \{0\} \times [0, +\infty]$ and distinguish various cases of relative behaviors of the elements of $s$ characterized by $\fkr$:
\begin{equation}\label{eq:q2}
\begin{dcases}
\text{For all $k$ in $\{1, 2\}$ there exists $\muk$ in $[0, +\infty]$ such that $\lim_{s \to \bar{s}} \, \mu/2 h\eps^{2(2-k)} =  \muk,$}\\
\bullet\; \fkr = 1 \quad \text{if }\; \left(\muo, \, \mut\right) \in [0, +\infty) \times \{0\}, \\
\bullet\; \fkr = 2 \quad \text{if }\; \left(\muo, \, \mut\right) \in \{+\infty\} \times [0, +\infty), \\
\bullet\; \fkr = 3 \quad \text{if }\; \left(\muo, \, \mut\right) \in \{+\infty\}^2,
\end{dcases}
\tag{$\tH_2$}
\end{equation}
with the additional assumption:
\begin{equation}\label{eq:q3}
\text{If } \; \fkq = 2, \quad \int_{\Omg^-} f \cdot \rho \itd x + \int_{\ts{\fkq}{\Gamma}{}{\tN} \cap
\{x_3 < 0\}} g \cdot \rho \itd \scH_2 = 0, \quad \forall \rho \in \clR
\tag{$\tH_3$}
\end{equation}
where $\clR$ denotes the set of rigid displacements ($\rho \in \clR \iff e(\rho) = 0$). Of course when $\mu$ remains bounded from below, this assumption may be dropped.

In the sequel, $C$ will denote various constants (independent of $s$!) which may vary from lines to lines. So, classically, $\puxs$ does satisfy:
\begin{align}
\int_{\pIOxh} \left|e^{\beam}\left(\eps, \puxs\right)\right|^2 \itd x &\le C \label{E:estimatione}
\\
\mu \int_{B_h}\left|e^{\beam}\left(\eps, \puxs\right)\right|^2 \itd x &\le C \label{E:11}
\end{align}

In the following, for any measurable function $\varphi$ defined in a subset $\phi$ of $\mbR^3$, we will denote the restriction of $\varphi$ to $\{\pm x_3 > 0\} \cap \phi$ by $\varphi^\pm$. As estimation~\eqref{E:estimatione} involves a domain which varies with $h$, it is convenient to associate with all $v$ in $\pVxh$, an element $\pTxh v$ in $H^1_{\pIGmaxD} (\Omg\setminus S, \mbR^3)$ given by:
\begin{equation}
\begin{dcases}
\bullet\; \text{when } \fkc = (1, 1), \\
\quad \left(\pTxh  v\right)^\pm(x) :=
\begin{cases}
v\left(\wht{x}, x_3 + \sgn(x_3)h\right) &\text{ if } \pm x_3 \in \left(0, L^\pm - h\right), \\
0 &\text{ if } \pm x_3 \in \left(L^\pm - h, L^\pm\right),
\end{cases} \\
\bullet\; \text{when } \fkc = (1, 2), \\
\quad \left(\pTxh v\right)^+(x) \text{ is as $\fkc = (1,1)$}, \\
\quad \left(\pTxh v\right)^-(x) :=
\begin{cases}
v(\wht{x}, x_3-h) & \text{ if } x_3 \in (-L^-+h, 0), \\
v(\wht{x}, -2L^- - x_3+h) & \text{ if } x_3 \in (-L^-, -L^-+h),
\end{cases} \\
\bullet\; \text{when } \fkc = (2,1) \text{ or } (2,2), \\
\quad\left(\pTxh v\right)^\pm(x) := v\left(\wht{x}, x_3 + \sgn(x_3)h\right).
\end{dcases}
\end{equation}
Let
\begin{equation}
\begin{dcases}
x^R := (-x_2, x_1), \quad \forall x \in \mbR^3, \\
\wht{e}_{\alpha\beta} := e_{\alpha\beta} \quad \forall e \in \mbS^3, \\
H^1_m(w) := \left\{w \in H^1(\omg);\; \int_\omg w(\wht{x}) \itd \wht{x} = 0 \right\}, \\
\clV^{\pm}_{\BN} := \left\{w \in H^1\left(\Omg^\pm, \mbR^3\right);\; e_{\alpha\beta}(w) = e_{\alpha 3}(w) = 0 \right\}, \\
\pUpD := \left\{w \in H^1\left(\Omg^+, \mbR^3\right);\; w = 0 \;\text{ on }\;S_{L^+}\right\}, \\
\pUmD := \left\{w \in H^1\left(\Omg^-, \mbR^3\right); w = 0\;\text{ on }\;S_{L^-}\; \text{if }\;\fkq=1\right\}, \\
\pUzpm := V^{\pm}_{\BN} \cap \ts{\fkq}{\clU}{\pm}{\tD}, \\
\pUonep := \biggl\{ w;\; \exists c^+ \in H^1_{L^+}(0, L^+) \;\text{ s.t. }\;\wht{w}(x) = c^+(x_3)x^R,\; w_3 \in L^2\left(0, L^+; H^1_m(\omg)\right)\biggr\}, \\
\pUonem := \biggl\{w;\; \exists c^- \in H^1(-L^-, 0),\; c^-(-L^-) = 0 \\
\mkern182mu \text{if }\;\fkq = 1\;\text{ s.t. }\;
\wht{w}(x) = c^-(x_3)x^R, w_3 \in L^2\left(-L^-, 0; H^1_m(\omg)\right) \biggr\}, \\
\pUtwop := \Biggl\{w;\; \wht{w} \in L^2\left(0, L^+;\; H^1_m\left(\omg, \mbR^2\right)\right) \text{ s.t. } \\
\mkern255mu \int_\omg x^R \cdot \wht{w}(\wht{x}) \itd \wht{x} = 0\quad\text{a.e. }\;x_3 \in (0, L^+),\; w_3 = 0 \Biggr\}, \\
\pUtwom := \Biggl\{w;\; \wht{w} \in L^2\left(-L^-, 0; H^1_m\left(\omg, \mbR^2\right)\right) \;\text{ s.t.} \\
\mkern240mu\int_\omg x^R \cdot \wht{w}(\wht{x}) \itd \wht{x} = 0 \quad\text{a.e. }\; x_3 \in (-L^-, 0),\; w_3 = 0 \Biggr\}, \\
\pUi := \biggl\{u; u^\pm \in \ts{\fkq}{\clU}{i\pm}{} \,\biggr\}, \quad \pU := \pUzero \times \pUone \times \pUtwo,\\
E^\pm_u :=
\begin{bmatrix}
\widehat{e}\left(u^{2\pm}\right) & e_{\alpha 3} \left(u^{1\pm}\right) \\
\sym & e_{33}\left(u^{0\pm}\right)
\end{bmatrix}, \quad \forall u^\pm = \left(u^{0\pm}, u^{1\pm}, u^{2\pm}\right) \in \pUzpm \times \pUonepm \times \pUtwopm.
\end{dcases}
\end{equation}
where there is no doubt we will write indifferently $E_{u^\pm}$ or $E^\pm_u$.

Recalling that $\clR$ stands for the space of the rigid displacements, we let:
\begin{equation}
\begin{dcases}
\pRp = \{0\}, \quad \forall \fkq \in \{1,2\}, \\
\pRm =
\begin{cases}
\{0\} & \text{if $\fkq = 1$}, \\
\left\{ \rho^{-} \text{ s.t. $\; \rho \in \clR$}\right\} & \text{if $\fkq = 2$}.
\end{cases}
\end{dcases}
\end{equation}

To obtain our main results (see Theorems~\ref{T:1}--\ref{T:2} below) we will use the method of variational convergence by first establishing a lower bound for the strain energy of the assembly and next an upper bound for the strain energy of the assembly for sequences of scaled displacements satisfying~\eqref{E:estimatione}--\eqref{E:11}. This will be done under the additional assumption
\begin{equation}\label{Q4}
\sup_{s \in \fkS} h/\eps^2 < +\infty \tag{$\tH_4$}
\end{equation}
which is precisely the premise that underpins the title of our article.

\begin{lemm}\label{L:1}
For all sequences $\{\pvxs\}_{s \in \fkS}$ in $\pVxh$ such that $\int_{\pIOxh}|e^{\beam}(\eps, \pvxs)|^2 \itd x \le C$, there exists $v = (v^0, v^1, v^2)$ in $\pU$ and a not relabeled subsequence such that $(\pTxh \; \pvxs)^\pm$ converges weakly in $H^1(\Omega^\pm, \mbR^3) / \pRpm$ and
\begin{enumerate}\romanenumi
\item\label{lemm1.1} $ \liminfs \int_{\pIOxh} \clbdW_{\ptauxh a} (e^{\beam}(\eps, \pvxs)) \itd x \ge \int_{\Omg \setminus S} \clbdW_a(E_{v}) \itd x, $
\item\label{lemm1.2} the traces $\gamma_{\pm h} (\pvxs)$ on $S_{\pm h}$ converge strongly in $L^2(S, \mbR^3) / \pRpm$ toward the traces $\gamma_{0\pm}(v^{0\pm})$ of $v^{0\pm}$ on $S$.
\end{enumerate}
\end{lemm}

\begin{proof}
Clearly there exists some $v^{0+}$ in $H^1_{S_{L^+}} (\Omg^+, \mbR^3)$ and a not relabeled subsequence such that $(\pTxh \; \pvxs)^+$ weakly converges in $H^1(\Omega^+, \mbR^3)$ toward $v^{0+} $. Note that when $\fkp = 1$, there exists $w$ in $H^1(\Omega^+_\eta, \mbR^3)$ such that, up to a new subsequence, $(\pvxs)^+$ converges weakly in $H^1(\Omega^+_\eta, \mbR^3)$ toward $w$ for all $\eta$ in $(0, L^+)$. As $w = v^{0+}$ in $\Omega_\eta^+$, one has $w = v^{0+}$ which then belongs to $H^1_{S_{L^+}} (\Omg^+, \mbR^3)$. Moreover the trace on $S$ of $(\pTxh \, \pvxs)^{+}$ converges strongly in $L^2(S, \mbR^3)$ toward the trace $\gamma_{0+}(v^{0+})$ of $v^{0+}$ on $S$. Lastly as $\int_{\Omg^+} |e^{\beam} (\eps, (\pTxh \, \pvxs)^+)|^2 \itd x = \int_{\ts{\fkp}{\Omg}{+}{h}} |e^{\beam}(\eps, \pvxs)|^2 \itd x$, by proceeding as in~\cite{Gaudiello2007,Murat1999, Murat2000}, one deduces that there exists $v^+ = (v^{0+}, v^{1+}, v^{2+})$ in $\ts{\fkq}{\clU}{0+}{} \times \pUonep \times \pUtwop$ such that $e^{\beam}(\eps, (\pTxh \, \pvxs)^+)$ converges weakly in $L^2(\Omg^+, \mbS^3)$ toward some $E_{v^+}$ so that
\[
\liminf_{s \to \bar{s}} \quad \int_{\ts{\fkp}{\Omg}{+}{h}} \clbdW_{\ptauxh a}\left(e^{\beam}\left(\eps, \pvxs\right)\right) \itd x \ge \int_{\Omg^+} \clbdW_a(E_{v^+}) \itd x.
\]

We may proceed in the same way in $\ts{\fkp}{\Omg}{-}{h}$ when $\fkq = 1$, while when $\fkq=2$, one has to consider a rigid displacement.
\end{proof}


For any $v$ of $H^1(\Omg\setminus S, \mbR^d)$, $1 \le d \le 3$, let $[v] = \gamma_{0+}(v^+) - \gamma_{0-}(v^-)$ be the jump of $v$ across $S$ where $\gamma_{0\pm}$ is the trace operator from $H^1(\Omg^\pm, \mbR^d)$ into $L^2(S, \mbR^d)$ and
\begin{equation}\label{eq13_}
\pqV := \left\{w \in H^1_{\pIGmaxD} \left(\Omg\setminus S; \mbR^3\right);\quad [\wht{w}] = 0 \;\text{ if }\;\fkr = 2,\quad [w] = 0 \;\text{ if }\; \fkr =3\right\}.
\end{equation}
Let $\pqJ$ and $\Wqrb$ respectively defined on $\pqV$ and $\mbR^3$ by:
\begin{align}\label{14f}
\pqJ(w) &:=
\begin{cases}
\displaystyle\mur \int_S \Wqrb \left([w] \right) \, d\wht{x} & \text{if }\;(\fkq, \fkr) \in \{1, 2 \}^2, \\
0 & \text{if }\;(\fkq, \fkr) \in \{1, 2 \} \times \{3\},
\end{cases} \\
\Wqrb (z) &:= \Minn
\left\{ \clbdW_{b} \left(\zeta \otS e_3\right);\;
\begin{cases} 
\zeta = (\hat{z}, 0) & \text{if } \; (\fkq, \fkr) \; = (1,1)\\
\hat{\zeta} = \hat{z} & \text{if } \; (\fkq, \fkr) \; = (2,1)\\
\zeta_3 = z_3 & \text{if } \; \fkr=2\\
\end{cases}
\right\}.
\end{align}

Under the convention $\infty \times 0 = 0$, one has:
\begin{equation}\label{eq18}
\pqJ(w) = \sum_{\fki=1}^{2} \; \mui \int_S \Wqib ([w]) \, d\hat{x}.
\end{equation}

\begin{lemm}\label{L:2}
For all sequences $\{\pqvxs\}_{s \in \fkS}$ satisfying~\eqref{E:estimatione}--\eqref{E:11}, the field $v^0$ supplied by Lemma~\ref{L:1} belongs to $\pqV$ and $v^{0-}$ is defined up to an element of
\begin{equation}\label{17_}
{}^{\fkq, \fkr}\clR^-_0 := \left\{ \rho^- \in \pRm ; \quad \gamma_{0-} (\widehat{\rho^-}) = 0 \;\text{ if }\; \fkr = 1 \quad \text{and} \quad \muo >0, \; \rho^- = 0 \;d\text{ if }\; \fkr = 2, 3 \right\}.
\end{equation}
We have
\begin{equation}\label{E:18}
\liminfs \mu \int_{B_h} \clbdW_b\left(e^{\beam}\left(\eps, \pqvxs\right)\right) \itd x \ge \pqJ\left(v^0\right).
\end{equation}
\end{lemm}

Note that as $v^{0\pm}$ belongs to $\pUzpm$, $\Wqrb ([v^0])$ reduces either to a constant or to a quadratic function if $\fkr = 1$ or $2$, respectively.
\begin{proof}
Following~\cite{ciarlet} we proceed to the following change of coordinates and fields:
\begin{equation}\label{19b}
\begin{alcases}
 x = (\hat{x}, x_3) \in B_h &\mapsto y = (\hat{x}, x_3 / h) \in B_1 := \omega \times (-1, 1),\\
 v \in H^1(B_h,\mbR^3) &\mapsto V \in H^1\left(B_1,\mbR^3\right) \; \text{ s.t. } \; v(x) = \left(h \hat{V}(y), V_3(y) \right),
\end{alcases}
\end{equation}
so that
\begin{equation}\label{20b}
\begin{dcases}
e^{\beam}(\varepsilon,v)(x) = h / \varepsilon^2 \; e^{\plate}(h / \varepsilon, V)(y),\\
e^{\plate}(t, V) :=
\begin{bmatrix}
\hat{e}(V) & e_{\alpha 3}(V)/t \\
\sym & e_{33}(V)/t^2
\end{bmatrix},
\quad \forall V \in H^1\left(B_1, \mbR^3\right), \\
e_{ij} (V) := \frac{1}{2} \left(\frac{\partial V_i}{\partial y_j} +\frac{\partial V_j}{\partial y_i}\right), \quad 1 \leq i, j \leq 3,\\
\int_{B_h} \clbdW_b \left(e^{\beam}(\varepsilon,v) \right) \, dx = \frac{h}{\varepsilon^4} \int_{B_1} \clbdW_b \left(e^{\plate}(h/\varepsilon, W)\right)\, dy, \\
W := h V.
\end{dcases}
\end{equation}
Let $^{\fkc,\fkr} W_s := h \;\, ^{\fkc,\fkr}V_s$, its traces $\gamma_{\pm 1}(^{\fkc,\fkr} W_s)$ on $S_{\pm 1}$ satisfy:
\begin{equation}\label{21b}
\begin{cases}
\gamma_{\pm 1}\left(\widehat{{}^{\fkc,\fkr} W_s}\right) = h \gamma_{\pm 1}\left(\widehat{^{\fkc,\fkr} V_s}\right) = h \gamma_{\pm h} \left(\widehat{{}^{\fkc,\fkr} v_s} / h\right) \longrightarrow \gamma_{0\pm} \bigl(\widehat{v^{0\pm}}\bigr), &\text{with }\; \partial_\alpha \widehat{v^{0\pm}} = 0 \\
&\text{in }\; L^2\left(S, \mbR^3\right) / \pRpm.\\
\gamma_{\pm 1}\left({}^{\fkc,\fkr} W_{s_3}\right) = h  \gamma_{\pm 1}\left({}^{\fkc,\fkr} V_{s_3}\right) = h\gamma_{\pm h} \left({}^{\fkc,\fkr} v_{s_3}\right) \longrightarrow 0
\end{cases}
\end{equation}

\begin{proof}[{$\bullet$ Case $\fkr=1 \quad \bigl((\muo, \,\mut) \in [0,+\infty) \times \{0\} \bigr)$}]
The result being obvious if $\muo=0$, we assume $\muo >0$. Denoting $^{\fkc, 1} \Psi_s := h / \varepsilon \; ^{\fkc, 1} W_s$, \eqref{E:11}--\eqref{20b} imply:
\begin{equation}\label{22b}
C \frac{2h \varepsilon^2}{\mu} \geq \int_{B_1}
\left|e^{\plate} \left(h / \varepsilon, {}^{\fkc, 1} \Psi_s\right)\right|^2 \, dy.
\end{equation}
Hence, up to a non-relabeled subsequence, $e^{\plate}(h/\varepsilon, {}^{\fkc, 1} \Psi_s)$ weakly converges in $L^2 (B_1, \mbS^3)$ so that by Jensen inequality one has:
\begin{equation}\label{23b}
\liminf_{s \to \bar{s}} \; \mu \int_{B_h} \clbdW_b \left(e^{\beam}\left(\varepsilon, {}^{\fkc, 1} v_s\right) \right) \, dx \ge \muo \int_{S} \clbdW_b \left({}^{\fkc, 1} \ell \right)
\, d\hat{x},
\end{equation}
where $^{\fkc, 1} \ell$ is the weak limit in $L^2(S, \mbS^3)$ of $\int_{-1}^{+1} e^{\plate} (h/\varepsilon, {}^{\fkc, 1} \Psi_s) \, dy_3$. Moreover, because of\linebreak\eqref{20b}--\eqref{21b}, ${}^{\fkc, 1} \Psi_s$ does converge weakly in $H^1 (B_1, \mbR^3)$ toward some Kirchhoff--Love displacement ${}^{\fkc, 1} \Psi$ with vanishing traces on $S_{1}$, so that ${}^{\fkc, 1} \Psi=0$. Hence ${}^{\fkc, 1} \ell_{\alpha \beta}$ vanishes. Next, let $\tau$ in $C^\infty_0 (S, \mbS^3)$ such that $\tau_{\alpha \beta} = \tau_{33} = 0$. The following estimation (see~\cite[eq.~(4.20) p.~695]{lichtetal}):
\begin{equation}
\int_{B_h} \left|{}^{\fkc,1}\!v_{s3}\right|^2 \itd x \le C \left[h\int_{S_h} \left|{}^{\fkc,1}\!v_{s3}\right|^2 \itd \wht{x} + h^2 \int_{B_h} \left|e\left({}^{\fkc,1}\!v_s\right)\right|^2 \itd x + h^2 \int_{\ts{\fkp}{\Omg}{+}{h}} \left|e\left({}^{\fkc,1}\!v_s\right)\right|^2 \itd x
\right]
\end{equation}
and an integration by parts in $\int_S \tau(\hat{x}) \int_{-1}^{+1} e^{\plate} (h/\varepsilon, {}^{\fkc, 1} \Psi_s) \, dy_3 
\itd  \hat{x}$ implies that $\gamma_{+1}
(\widehat{{}^{\fkc, 1} W}_s) - \gamma_{-1} (\widehat{{}^{\fkc, 1} W}_s)$ has a weak limit in $L^2 (S, \mbR^2)$. Therefore, modulo a modification of the field $v^{0-}$ implied in the proof of Lemma~\ref{L:1} (and defined up to a rigid displacement), we get $^{\fkc, 1} \ell_{\alpha 3}= \frac{1}{2}[v^0_\alpha]$ by due account of~\eqref{21b}. Now by taking $\tau$ in $C^\infty_0 (S, \mbS^3)$ such that $\tau_{\alpha \beta} = \tau_{\alpha 3} = 0$, one deduces that $^{\fkc, 1} \ell_{33}=0$ if $q=1$ and the results stated in Lemma~\ref{L:2} stem from the very definition of $\clbdW_{^{\fkq, 1}\tilde{b}}$.
\let\qed\relax
\end{proof}

\begin{proof}[{$\bullet$ Case $\fkr=2 \quad \bigl((\muo, \, \mut) \in \{+\infty\} \times [0,+\infty) \bigr)$}]
We first establish that $\muo = +\infty$ implies $[\widehat{v^0}] = 0$. Inequality~\eqref{22b} implies that $e^{\plate}(h / \varepsilon, \, {}^{\fkc, 1} \Psi_s)$ converges strongly to zero in $L^2(B_1,\mbS^3)$ and ${}^{\fkc, 1} \ell = 0$, thus $[v^0_\alpha]$ which has been seen equal to ${}^{\fkc, 1} \ell_{\alpha 3}$ vanishes. This proves the result when $\mut = 0$.

Next when $\mut > 0$, \eqref{E:11}, \eqref{20b}, \eqref{21b} and assumption~\eqref{Q4} imply that $e^{\plate}(h / \varepsilon, {}^{\fkc, 2} \Psi_s)$ is bounded in $L^2(B_1,\mbS^3)$ if ${}^{\fkc, 2} \Psi_s := h / \varepsilon^2 {}^{\fkc, 2} W_s$ and, up to a non relabeled subsequence, ${}^{\fkc, 2} \Psi_s$ converges in $H^1(B_1,\mbR^3)$ toward some Kirchhoff--Love displacement ${}^{\fkc, 2} \Psi$ whose third component of the traces on $S_{+ 1}$ vanishes, so that the flexural part ${}^{\fkc, 2} \Psi^{\tF}$ of ${}^{\fkc, 2} \Psi$ is equal to zero. Moreover the trace of $\widehat{{}^{\fkc, 2} \Psi}_s$ on $S_{+1}$ converges in $L^2(S, \mbR^2)$ toward a constant by~\eqref{21b} hence $e({}^{\fkc, 2} \Psi)$ vanishes.

Finally~\eqref{20b}--\eqref{21b} imply that $ \int_{-1}^{+1} (\frac{\varepsilon}{h})^2 \frac{\partial}{\partial y_3}\bigl(\frac{h}{\varepsilon^2} {}^{\fkc, 2} W_{s_3}\bigr) \itd y_3$ which is equal to $\gamma_{+h} ({}^{\fkc, 2} v_{s_3}) - \gamma_{-h} ({}^{\fkc, 2} v_{s_3})$ converges strongly in $L^2(S)$ toward $[v^0_3]$ and the result is obtained through the definition of $\clbdW_{{}^{\fkq, 2}\tilde{b}}$.
\let\qed\relax
\end{proof}

\begin{proof}[{$\bullet$ Case $\fkr=3 \quad \bigl(\muo  =  \mut = +\infty  \bigr)$}]
As $e^{\plate} (h / \varepsilon, h / \varepsilon^2 {}^{\fkc, 2} W_s)$ converges strongly to zero in $L^2(B_1,\mbS^3)$, one has $[v_3^0] = 0$ and the proof of Lemma~\ref{L:2} is complete.
\end{proof}\let\qed\relax
\end{proof}



\begin{lemm}\label{L:3}
For all $v = (v^0, v^1, v^2)$ in $\pU$ such that $v^0$ belongs to $\pqV$ there exists a sequence $\{\pqvxs\}_{s \in \fkS}$ in $\pVxh$ such that
\begin{multline}\label{E:21}
\limsups \int_{\pIOxh} \clbdW_{\ptauxh a} \left(e^{\beam}\left(\eps, \pqvxs\right) \right) \itd x + \mu \int_{B_h} \clbdW_b \left(e^{\beam}\left(\eps, \pqvxs\right) \right) \itd x\\
\leq \int_{\Omega \setminus S} \clbdW_a (E_v)
\, dx + \pqJ(v^0).
\end{multline}
\end{lemm}
\begin{proof}
As $\pqJ$ is continuous on $\pqV$, it is enough to assume that $(v^{0 \pm}, \wht{v^{1 \pm}})$ belongs to $C^\infty (\bar{\Omega}^\pm ; \mbR^5)$ and that $(v^1_3, \wht{v^2})$ belongs to $H^1 (-L^-, L^+ ; H^1_m (\omega, \mbR^3) )$ with a compact support in $(-L^-, 0) \cup (0, L^+)$. For any $w$ in $H^1_m(\Omg; \mbR^3)$ let $R_h \, w$ be defined by
\begin{equation}
R_hw =
\begin{dcases}
w^{\tS}(x) + \frac{|x_3|}{h} w^{\tA}(x), & \forall x \in B_h, \\
w(x), & \forall x \in \pIOxh,
\end{dcases}
\end{equation}
where
\begin{equation}
w^{\tS}(x) = \frac{1}{2}\bigl(w(\wht{x}, x_3) + w(\wht{x}, -x_3) \bigr), \quad w^{\tA}(x) = \frac{1}{2} \bigl(w(\wht{x}, x_3) - w(\wht{x}, -x_3) \bigr), \quad \forall x \in B_h.
\end{equation}
Let $\eta$ in $C^\infty_0(-L^-, L^+)$ with $\eta = 1$ if $|x_3| < h_0$, $\eta = 0$ if $|x_3| > (1/2) (h_0+\Minn(L^+, L^-))$ and $\rho$ in $C^\infty_0 (-1, 1)$ such that $ 2 \int_{-1}^{+1} (\rho^{\prime}(y) )^2 \, dy = 1$. Considering first the cases $\fkr = 1, 2$, as there exists some $(\chi, \xi)$ in $L^2 (S, \mbR^2) \times L^2(S)$ such that $\clbdW_{^{\fkq, 1}\tilde{b}} ([v^0]) = \clbdW_b (([\widehat{v^0}], \xi) \otS e_3)$, $\clbdW_{^{\fkq, 2} \tilde{b}} ([v^0]) = \clbdW_b ((\chi, [v^0_3]) \otS e_3 )$, let $^{\fkc, \fkr} \tilde{v}_s$ be such that
\begin{equation}\label{E:24}
^{\fkc, \fkr} \tilde{v}_s :=
\begin{dcases}
R_h \left(v^0 + \varepsilon v^1\right) +
{}^{\fkr} \theta_s&\text{in }\;B_h,\\
{}^{1} \theta_s := \varepsilon^{-1} \; \rho (\cdot / h) (0, \xi), \quad {}^{2} \theta_s := \varepsilon \; \rho (\cdot / h)\left(\chi - [\widehat{v^1}], 0\right),\\
\twotauxh \Bigl(\eta \left(v^0 + \varepsilon v^1 + \varepsilon^2 v^2\right) \Bigr) + \ptauxh \Bigl((1-\eta) \left(v^0 + \varepsilon v^1 + \varepsilon^2 v^2\right) \Bigr) &\text{in }\;\pIOxh.
\end{dcases}
\end{equation}
Clearly $\int_{{}^{\fkp} \Omega_h} \clbdW_{\ptauxh a} (e^{\beam} (\varepsilon, {}^{\fkc, \fkr} \tilde{v}_s))\, dx$ converges to $ \int_{\Omega \setminus S} \clbdW_a (E_v)\, dx$ while a rather tedious computation using \emph{all the previous assumptions} on $v, \eta, \rho$ and the proof of~\cite[Lemma~4.1]{lichtetal} shows that:
\begin{equation}
\lim_{s \to \bar{s}} \quad \mu \int_{B_h} \clbdW_b \left(e^{\beam} \left(\varepsilon, {}^{\fkc, \fkr} \tilde{v}_s\right) \right) \, dx =  {}^{\fkc, \fkr} \clbdJ\left(v^0\right).
\end{equation}

When $\fkr=3$ (and also $\omg \times (-h,h) \cap \supp (v^2) = \emptyset $), taking into account that $v^{0\pm}$ belongs to $\clV^{\pm}_{\BN}$ with $[v^0] = 0$ and that consequently there exists $(v^{0T}, v^{0L})$ in $H^2(-L^-, L^+; \mbR^2) \times H^1(-L^-, L^+)$\linebreak such that
\[
v^0(x) = \left(v^{0T}(x_3),\; v^{0L}(x_3) - x_\alpha \frac{dv^{0T}_\alpha}{d x_3} (x_3) \right),
\]
we introduce $^{\fkc, 3} \tilde{v}_s$ as:
\begin{equation}\label{30_f}
^{\fkc, 3} \tilde{v}_s :=
\begin{dcases}
\left(v^{0T}(0) + x_3 \left(\frac{dv^{0T}}{dx_3}\right)(0), v^{0L}(0) - x_\alpha \frac{dv^{0T}_\alpha}{dx_3}(0) \right), & \text{if }\;x \in B_h, \\
\twotauxh \Bigl(\eta \left(v^0 + \varepsilon v^1 + \varepsilon^2 v^2\right) \Bigr) + \ptauxh \Bigl((1-\eta) \left(v^0 + \varepsilon v^1+ \eps^2v^2\right) \Bigr) \\
\mkern200mu + h \sgn (x_3) \left(0, \frac{dv^{0T}}{dx_3}(0) \eta(x_3) \right) & \text{if }\; x \in \pIOxh.
\end{dcases}
\end{equation}
As $e^{\beam}(\eps, {}^{\fkc, 3} \tilde{v}_s)$ vanishes in $B_h$, one has the expected result.
\end{proof}

We therefore classically deduce:
\begin{theo} \label{T:1}
Problem~\eqref{eq:theo4} defined by
\begin{equation}\label{eq:theo4}
%\tag{$\,^{\fkq,\fkr}\tP_{\overline{\!s\!}}$}
\tag{$\ts{\fkq,\fkr}{\tP}{}{\overline{\!s\!}}$}
\Minn \left\{\int_{\Omg\setminus S} \clbdW_a(E_v) \itd x + \pqJ(v^0) - \int_\Omg f \cdot v^0 \itd x - \int_{\pIGmaxN} g \cdot v^0 \itd \scH_2 ;\; v \in \pU \; \text{ s.t. } \; v^0 \in \pqV\right\}
\end{equation}
has a unique solution $\ts{\fkq,\fkr}{u}{}{\bar{s}}$ with $\ts{\fkq,\fkr}{u}{0-}{\bar{s}}$ defined up to an element of ${}^{\fkq, \fkr} \clR_0^-$ (see~\eqref{17_}) and when $s$ goes to $\bar{s}$ the solution $\puxs$ of ($\pPxs$) converges toward $\ts{\fkq,\fkr}{u}{}{\bar{s}}$ in the sense
\begin{multline}\label{eq27}
\lims \int_{\pIOxh} \clbdW_{\ptauxh a} \left(e^{\beam}\left(\eps, \puxs\right) - \ts{\fkp}{\tau}{}{h} \; E_{\ts{\fkq,\fkr}{u}{}{\bar{s}}} \right) \, dx \\
+ \mu \int_{B_h} \clbdW_b \left(e^{\beam} \left(\eps, \puxs\right) - \left(\frac{\left[\left(\reallywidehat{\ts{\fkq,\fkr}{u}{0}{\bar{s}}} , / \eps, \, \ts{\fkq,\fkr}{u}{0}{\bar{s}_3}c\right) \right] + \left({}^{\fkr} \bar{\chi}, {}^{\fkq,\fkr} \bar{\xi} / \varepsilon \right)}{2h} \right) \otS e_3 \right) \itd x =0,
\end{multline}
where
\begin{align*}
{}^{\fkr} \bar{\chi} = 0 \quad\text{if}\quad\fkr \in \{ 1, 3 \},\quad {}^{\fkDeux} \bar{\chi} &= \Argmin \Bigl\{ \clbdW_b \bigl(\bigl(\chi, \bigl[\ts{\fkq, 2}{u}{0}{\bar{s}_3}\bigr]\bigr) \otS e_3 \bigr), \; \chi \in \mbR^2 \Bigr\},\\
 {}^{1, 1} \bar{\xi} = 0,\quad {}^{2, 1} \bar{\xi}
&= \Argmin \Bigl\{ \clbdW_b \bigl(\bigl(\bigl[\reallywidehat{\ts{2, 1}{u}{0}{\bar{s}}}\bigr], \xi\bigr) \otS e_3 \bigr), \; \xi \in \mbR \Bigr\} \intertext{and} {}^{\fkq, \fkr} \bar{\xi} &= 0 \quad\text{if}\quad (\fkq, \fkr) \in \{ 1, 2 \} \times \{2, 3\}.
\end{align*}
\end{theo}


Even if~\eqref{eq:theo4} involves \emph{abstract} fields defined in ``abstract beams'' occupying $\Omega \setminus S$, we will use the language of Mechanics to comment it. Problem~\eqref{eq:theo4} accounts for the equilibrium of two solids filling the cylindrical domains $\Omega^+$ and $\Omega^-$ made of a ``generalized linearly elastic'' material whose state is described by its field of displacement $v^0$ and two additional variables $v^1$ and $v^2$. Each element of $v := (v^0, v^1, v^2)$ has a specific kinematics, especially $v^0$ which is of Bernoulli--Navier type, while the stress tensor $\Sigma$ does satisfy:
\begin{equation}
\Sigma = a \, E_v, \quad E_v =
\begin{bmatrix}
\hat{e}(v^2) & e_{\alpha 3}(v^1) \\
\sym & e_{33}(v^0)
\end{bmatrix}.
\end{equation}
These solids are subjected to forces $f$, $g$ and clamped on $\ts{\fkq}{\Gamma}{}{\tD} $ with a \emph{constraint} on $S$ which reads as:
\begin{itemize}
\item $\muo = 0$: free to separate.
\item $\muo > 0$: elastic pull-back with $(2\,\Sigma_{\tT}, \Sigma_{\tN}) = \, \muo \; D\clbdW_{^{\fkq, 1} \tilde{b}} ([v^0])$, with:
\begin{itemize}
\item $\Sigma_{\tT} = \widehat{\Sigma \, e_3}$: the tangential stress,
\item $\Sigma_{\tN} = (\Sigma \, e_3)_3$: the normal stress.
\end{itemize}
\item $(\muo, \, \mut) \in \{ +\infty \} \times [0, +\infty)$:
\begin{enumerate}\romanenumii
\item tangential adhesion i.e. $[\widehat{v^0}] = 0$,
\item elastic normal pull-back i.e. $\Sigma_\tN = \mut \; D\clbdW_{^{\fkq, 2} \tilde{b}} ([v^0])$.
\end{enumerate}
\item $(\muo, \mut) \in \{ +\infty \}^2$: perfect adhesion i.e. $[v^0] = 0$.
\end{itemize}
Roughly speaking, the thin layer $B_h$ at the limit shrinks to the surface $S$ and is replaced by a mechanical constraint which does not involve $\fkp$ which only accounts for the length of the cylindrical domains occupied by the scaled beams.


We moreover obviously have:
\begin{theo} \label{T:2}
When $s$ goes to $\os$, the solution $\pus$ of $\ (\pPs)$ satisfies:
\begin{multline*}
\lims \frac{1}{\eps^4}
\left[\int_{ \pIOeh} \clbdW_{\paeh} \left(e^\eps(\pus) - {}^{\fkp}E^\eps_{\ts{\fkq,\fkr}{u}{}{\os}} \right) \itd x^\eps\right. \\
+\left. \mu \int_{B^\eps_h} \clbdW_b \left(e^\eps(\pus) -\frac{\left[\ts{\fkq,\fkr}{\ou}{0\eps}{}\right] + \left({}^{\fkr} \bar{\chi}^\varepsilon, {}^{\fkq, \fkr} \bar{\xi}^\varepsilon\right)}{2h} \otS e_3 \right) \itd x^\eps \right] = 0,
\end{multline*}
where $^{\fkp}E^\eps_{\ts{\fkq,\fkr}{u}{}{\os}} (x^\varepsilon) = \varepsilon \; \ts{\fkp}{\tau}{}{h} \, E_{\ts{\fkq,\fkr}{u}{}{\os}} (x)$ for all $x$ in $^{\fkp} \Omega_h$, $\ts{\fkq,\fkr}{\ou}{0\eps}{} := \clS_{\varepsilon}^{-1} \; \ts{\fkq,\fkr}{u}{0}{\bar{s}}$ and $({}^{\fkr} \bar{\chi}^\varepsilon, {}^{\fkq, \fkr} \bar{\xi}^\varepsilon) (\hat{x}^\varepsilon) = \varepsilon \; ({}^{\fkr} \bar{\chi}, {}^{\fkq, \fkr} \bar{\xi}) (\hat{x})$ for all $\hat{x}$ in $\omega$.
\end{theo}

It should be noted that similarly to the modeling of a single slender beam~\cite{terap}, it is the limit ${}^{\fkq, \fkr} u_{\bar{s}}^0$ of the scaled displacement which is a Bernoulli--Navier displacement and that the true strain tensor is not close to a Bernoulli--Navier strain tensor but should involve ${}^{\fkq, \fkr} u_{\bar{s}}^1$ and ${}^{\fkq, \fkr} u_{\bar{s}}^2$. Here however it is the \emph{sole} jump of the displacement ${}^{\fkq, \fkr} u_{\bar{s}}^0$ which is involved in the energetic relative approximation of the real strain tensor in the very thin layer, with two critical sizes corresponding to the jump of the $~\widehat{\cdot}~$ component and the third component respectively, a ``natural extension'' of the case of the assembly of two three-dimensional linearly elastic bodies linked by a soft layer~\cite{Licht1997}. But contrary to the case of an assembly of two $3$-dimensional bodies linked by a stiff layer~\cite{Licht2016}, $ 1/\eps^4 \int_{B^\eps_h} \clbdW_b (e^\eps(\pus)) \itd x^\eps$ vanishes when $\mu/2h \to +\infty$!



Finally, we might also consider a more complex situation where the elasticity tensor of the joint presents some partial dependencies with respect to a quadruple of parameters $\tilde{s} := (\varepsilon, h, \mu_1, \mu_2)$. We may assume as an example that $\clbdW_b := \mu_1 \clbdW_{b_1} + \mu_2 \clbdW_{b_2} $, with at least one of the two quadratic forms $\clbdW_{b_1}$ and $\clbdW_{b_2}$ being positive definite. Introducing $^k \bar{\mu}_\alpha$ as in~\eqref{eq:q2}, we therefore obtain up to 9 different limit energy densities, i.e. such a case leads to additional asymptotic regimes in which the constraint on $S$ would be of a mixed nature: adhesion in certain directions but elastic pull-back or free separation in others.



\section{A proposal of model}

Similarly to~\cite{terap} we do not use $^{\fkq, \fkr} \bar{u}^{0 \varepsilon}$ as a simple and accurate enough model for the physical situation described by~\eqref{eq:CPs} but we have to take into account two new facts. First the field in question has to be an element of $^{\fkc} V_h^\varepsilon$, the space of admissible displacements for~\eqref{eq:CPs}. Second,\pagebreak{} in practice we do not know what are the limits $(\muo, \, \mut)$, we only know what is $s = (\varepsilon, h, \mu)$! We therefore follow a two-step process:
\begin{enumerate}
\item Taking into account~\eqref{14f}, \eqref{eq18} and~\eqref{E:21}, the problem~\eqref{eq:qQs}
\begin{multline}\label{eq:qQs}
\tag{$\pQxs$}
\Minn \left\{\int_{\Omg\setminus S} \clbdW_a(E_v) \itd x + \sum_{\fkr = 1}^{2}\frac{\mu}{2h\eps^{2(2-\fkr)}} \int_S \clbdW_{^{\fkq, \fkr}\tilde{b}} \left(\left[v^0\right]\right) \itd \wht{x} - \int_\Omg f \cdot v^0 \itd x - \int_{\pIGmaxN} g \cdot v^0 \itd \scH_2 \right.\\
\left. \text{with } \; v = \left(v^0, v^1, v^2\right) \in \pU \; \text{ s.t. } \; v^0 \in H^1_{\pIGmaxD}\left(\Omg\setminus S, \mbR^3\right) \right\}
\end{multline}
has a unique solution (from which it is easy to get a numerical approximation by standard finite element methods) $\ts{\fkq}{\widecheck{u}}{}{s}$ which is obviously close to $\ts{\fkq,\fkr}{u}{}{\os}$.
\item One approximates $\ts{\fkq}{\widecheck{u}}{}{s}$ by a smooth enough field in order to proceed to the construction used in the proof of Lemma~\ref{L:3} to obtain through $\clS_\varepsilon^{-1}$ a field $^{\fkc}\tilde{u}_s^\varepsilon$ in $^{\fkc} \clV_h^\varepsilon$ which is close to $\pus$ in terms of relative energy gap (and similarly for the associated stresses due to~\eqref{3f}\linebreak and~\eqref{eq:q1}). A ``numerical version'' of this step is again easy to perform.
\end{enumerate}



\section{A slight variant: two hollow beams connected by a solid joint}


In the previous analysis we did not require the cross-section $\omega$ to be simply-connected. To handle more realistic cases, we now assume that the joint still occupies the domain $B_h^\varepsilon := \varepsilon \omega \times (-h, h)$, while the two beams occupy the domains $^\fkp\Omega_h^\varepsilon$ defined via~\eqref{eq1_} from $\Omega^\varepsilon$ which now reads as $\varepsilon \, (\omega \setminus \kappa) \times (-L^-, L^+)$. Here, $\kappa$ is a compact subset of $\omega$ whose interior $\mathring{\kappa}$ is a simply-connected domain with a Lipschitz boundary and satisfying $\int_\kappa x_\alpha \, d\hat{x} = \int_\kappa x_1 x_2 \, d\hat{x} = 0$, three conditions henceforth also met by $\omega$.



Concerning the limit problem~\eqref{eq:theo4}, which is key to understanding the asymptotic behavior of the mechanical system, the sole change implied by this new configuration lies in the substitutions of $\Omega$ by $\Omega^{\eff} := (\omega \setminus \kappa) \times (-L^-,L^+)$ and of $S$ by $S^{\eff} := (\omega \setminus \kappa) \times \{ 0 \}$ in the integrals involved in the formulation of~\eqref{eq:theo4} and the definition~\eqref{14f} of $\pqJ$! The proof of Lemma~\ref{L:1} is indeed left unchanged while in the proof of Lemma~\ref{L:2} it suffices to bound the energy from below by zero in $B_h \setminus B_h^{\eff} := \kappa \times (-h,h)$. As to the proof of Lemma~\ref{L:3}, it has to be adapted somewhat by introducing, for $\delta$ positive and small enough, the cut-off function $\varphi_\delta$ defined on $\omega$ by $\varphi_\delta(\hat{x}) = \Max (1 - \dist(\hat{x},\omega \setminus \kappa) / \delta, \, 0 )$. The fields $^{\fkc, \fkr} \bar{v}_s$ defined in $B_h^{\eff}$ are then extended to $B_h \setminus B_h^{\eff}$ through formulas such as~\eqref{E:24} and~\eqref{30_f} so that introducing
\begin{equation}
{}^{\fkc, 1} \bar{\bar{v}}_s := \varphi_\delta {}^{\fkc,1} \bar{v}_s, \quad {}^{\fkc, 2} \bar{\bar{v}}_s := \varphi_\delta \; \left({}^{\fkc,2} \bar{v}_s - \left(\rwh{v^0(\cdot,0)}, 0\right) \right) + \left(\rwh{v^0(\cdot,0)}, 0 \right), \quad {}^{\fkc, 3} \bar{\bar{v}}_s :=  {}^{\fkc, 3} \bar{v}_s
\end{equation}
leads to
\begin{equation}
\lim_{s \to \bar{s}} \; \mu \int_{B_h \setminus B_h^{\eff}} \clbdW_b \left(e^{\beam} \left(\eps, {}^{\fkc, \fkr} \bar{\bar{v}}_s\right) \right) \, dx = 0
\end{equation}
by choosing $\delta$ of order $\varepsilon$.


Thus, although the ``scaled joint'' geometrically shrinks to the entire surface $S$, from a mechanical standpoint the limit surface energy remains concentrated on the effective part $S^{\eff}$. In particular, the limit surface energy density is the same as the one identified in the previous section and corresponds to the \emph{effective limit constraint} between the beams discussed after Theorem~\ref{T:1}.



\printCOI

\printbibliography
\end{document}
