Stretching and bending of (really) thick elastic plates

Résumé technique : Étirement et flexion de plaques élastiques (vraiment) épaisses

Énoncé du problème
L'énergie mécanique d'une plaque élastique est traditionnellement comprise comme se séparant en deux composantes distinctes : l'étirement et la flexion. Cette séparation est un résultat rigoureux pour les plaques asymptotiquement minces ($h = O(\varepsilon)$ avec $\varepsilon \ll 1$), dérivé d'hypothèses d'échelle où l'étirement de la ligne médiane est faible ($E = O(\varepsilon)$) et la flexion est finie ($K = O(1)$). Cependant, cette « image de la plaque mince » fournit souvent des approximations étonnamment précises pour les plaques épaisses, même lorsque l'épaisseur $h$ est de l'ordre de l'unité ($h = O(1)$) et que l'hypothèse de Kirchhoff (les normales restent normales) ne tient plus. Le dilemme central abordé par ce travail est de savoir si cette séparation de l'énergie en modes d'étirement et de flexion émerge formellement pour les plaques « vraiment épaisses » sous de petites déformations, et si oui, comment les modules effectifs se comportent.

Méthodologie
Les auteurs étudient ce problème en résolvant les équations de l'élasticité pour une plaque néo-Hookéenne incompressible en deux dimensions d'épaisseur $h = O(1)$. Ils imposent des déformations asymptotiquement petites sur la ligne médiane de la plaque, définies par une déformation d'étirement $E = O(\varepsilon)$ et une courbure $K = O(\varepsilon)$.

1. Cinématique et relations de comportement : La plaque est décrite à l'aide de coordonnées cartésiennes $(x, z)$ où $z \in [-h/2, h/2]$. La déformation est caractérisée par un tenseur de gradient de déformation $F$ dérivé du déplacement de la ligne médiane et du déplacement des points à travers l'épaisseur. Le matériau est modélisé comme un néo-Hookéen incompressible, avec une densité d'énergie élastique proportionnelle à $\text{tr}(F^T F) - 2$.
2. Développement asymptotique : Les auteurs effectuent un développement en perturbations régulières dans le paramètre de petite taille $\varepsilon$. Ils développent les champs de déplacement ($\tilde{u}, \tilde{w}$) et la pression $p$ en puissances de $\varepsilon$.
3. Conditions aux limites : Crucialement, les auteurs imposent des conditions aux limites de force nulle sur les surfaces supérieure et inférieure ($z = \pm h/2$) et une déformation prescrite sur la ligne médiane ($z=0$). Contrairement aux théories des plaques minces où la continuité de la ligne médiane est supposée, les auteurs permettent explicitement des discontinuités dans la solution à travers la ligne médiane ($z=0$) pour satisfaire les conditions aux limites aux surfaces, reconnaissant que l'hypothèse de Kirchhoff ne s'applique pas aux plaques épaisses.
4. Stratégie de résolution : Le problème linéaire de premier ordre est résolu en utilisant une décomposition en modes de Fourier ($e^{ikx}$). Cela réduit les équations aux dérivées partielles à un système d'équations différentielles ordinaires en $z$, qui sont résolues analytiquement. Les résultats sont validés par des simulations par éléments finis réalisées sous FEniCSx.

Contributions clés et résultats

• Séparation formelle de l'énergie : Les calculs démontrent que pour une plaque épaisse ($h = O(1)$ subissant de petites déformations de la ligne médiane), la densité d'énergie élastique totale se sépare toujours en un terme d'étirement et un terme de flexion. L'énergie prend la forme :
  $$E \propto \sum_k \left( h A(kh) E_k E_{-k} + h^3 B(kh) K_k K_{-k} \right)$$
  où $E_k$ et $K_k$ sont les amplitudes de Fourier des déformations d'étirement et de flexion, respectivement.
• Modules effectifs : Les auteurs dérivent des modules d'étirement ($A$) et de flexion ($B$) effectifs qui dépendent du paramètre sans dimension $\eta = kh$ (le rapport entre l'épaisseur de la plaque et la longueur d'onde de la déformation).
  • Lorsque $\eta \to 0$ (la limite de la plaque mince), $A \to 2$ et $B \to 1/6$, retrouvant les résultats classiques des plaques minces.
  • Lorsque $\eta$ augmente (plaques plus épaisses ou longueurs d'onde plus courtes), $A$ et $B$ diminuent tous deux, indiquant que la plaque s'adoucit tant en étirement qu'en flexion.
• Rapport constant des modules : Un résultat frappant est que le rapport entre les contributions d'étirement et de flexion, quantifié par le « nombre d'étirement » $S = A/12B$, reste approximativement constant pour $\eta \lesssim 1$. Plus précisément, $S - 1 \approx \eta^2/30$. Cela implique que bien que la rigidité absolue de la plaque diminue avec l'épaisseur, le coût énergétique relatif de l'étirement par rapport à la flexion reste presque identique à celui d'une plaque mince.
• Validation : La solution asymptotique montre un excellent accord avec les simulations par éléments finis, confirmant la validité des modules effectifs dérivés et la présence de discontinuités à travers la ligne médiane.

Signification
Cet article fournit une justification théorique à l'observation empirique selon laquelle les théories des plaques minces décrivent souvent avec précision la mécanique des plaques épaisses, telles que les feuillets cellulaires biologiques. Les auteurs soutiennent que la précision surprenante de ces théories provient du fait que la séparation des énergies d'étirement et de flexion persiste même dans la limite de petites déformations d'une plaque épaisse. De plus, l'importance relative de l'étirement par rapport à la flexion est préservée jusqu'à des échelles de longueur comparables à l'épaisseur de la plaque.

Par conséquent, le comportement mécanique des plaques élastiques épaisses peut être capturé par des théories de plaques minces, à condition de rescaler les forces externes pour tenir compte de l'adoucissement des modules effectifs. Cette conclusion est particulièrement pertinente pour la morphoélasticité et l'étude des systèmes biologiques (par exemple, les feuillets cellulaires de Volvox) où le matériau est épais et où les déformations induites par la croissance interviennent. Le travail suggère que le cadre conceptuel de l'étirement et de la flexion intrinsèques, largement utilisé dans les théories des coques minces, reste applicable aux structures biologiques épaisses.