Stretching and bending of (really) thick elastic plates

Technische Samenvatting: Rekken en Buigen van (Echt) Dikke Elastische Platen

Probleemstelling
De mechanische energie van een elastische plaat wordt traditioneel begrepen als een splitsing in afzonderlijke rek- en buigcomponenten. Deze scheiding is een rigoureus resultaat voor asymptotisch dunne platen ($h = O(\varepsilon)$ met $\varepsilon \ll 1$), afgeleid van schaalingsaannames waarbij de rek van de middenlijn klein is ($E = O(\varepsilon)$) en het buigen eindig is ($K = O(1)$). Deze "dunne-plaat-visie" biedt echter vaak verrassend nauwkeurige benaderingen voor dikke platen, zelfs wanneer de dikte $h = O(1)$ (orde grootheid) is en de Kirchhoff-hypothese (normalen blijven normaal) niet langer standhoudt. De centrale kwestie die in dit werk wordt behandeld, is of deze scheiding van energie in rek- en buigmodi formeel verschijnt voor "echt dikke" platen onder kleine vervormingen, en zo ja, hoe de effectieve moduli zich gedragen.

Methodologie
De auteurs onderzoeken dit probleem door de vergelijkingen van de elasticiteit op te lossen voor een tweedimensionale, incompressibele, neo-Hookeaanse plaat met een dikte $h = O(1)$. Zij leggen asymptotisch kleine vervormingen op aan de middenlijn van de plaat, gedefinieerd door een rekspanning $E = O(\varepsilon)$ en een kromming $K = O(\varepsilon)$.

1. Kinematica en Constitutieve Relaties: De plaat wordt beschreven met behulp van Cartesiaanse coördinaten $(x, z)$ waarbij $z \in [-h/2, h/2]$. De vervorming wordt gekenmerkt door de vervormingsgradiënttensor $F$, afgeleid van de verplaatsing van de middenlijn en de verplaatsing van punten door de dikte heen. Het materiaal wordt gemodelleerd als incompressibel neo-Hookeaans, waarbij de elastische energiedichtheid proportioneel is aan $\text{tr}(F^T F) - 2$.
2. Asymptotische Expansie: De auteurs voeren een reguliere perturbatie-expansie uit in de kleine parameter $\varepsilon$. Zij expanderen de verplaatsingsvelden ($\tilde{u}, \tilde{w}$) en de druk $p$ in machten van $\varepsilon$.
3. Randvoorwaarden: Cruciaal is dat de auteurs krachtvrije randvoorwaarden opleggen aan de boven- en onderzijde ($z = \pm h/2$) en een voorgeschreven vervorming op de middenlijn ($z=0$). In tegenwijl bij dunne-plaattheorieën de continuïteit van de middenlijn wordt verondersteld, staan de auteurs expliciet discontinuïteiten in de oplossing over de middenlijn toe ($z=0$) om aan de randvoorwaarden aan de oppervlakken te voldoen, waarbij zij erkennen dat de Kirchhoff-hypothese niet geldt voor dikke platen.
4. Oplossingsstrategie: Het lineaire leidende probleem wordt opgelost met behulp van een Fourier-modusdecompositie ($e^{ikx}$). Dit reduceert de partiële differentiaalvergelijkingen tot een systeem van gewone differentiaalvergelijkingen in $z$, die analytisch worden opgelost. De resultaten worden gevalideerd tegen eindige-elementen-simulaties uitgevoerd in FEniCSx.

Belangrijkste Bijdragen en Resultaten

• Formele Scheiding van Energie: De berekeningen tonen aan dat voor een dikke plaat ($h=O(1)$ die kleine vervormingen van de middenlijn ondergaat), de totale elastische energiedichtheid nog steeds splitst in een rekterm en een buigterm. De energie neemt de vorm aan:
  $$E \propto \sum_k \left( h A(kh) E_k E_{-k} + h^3 B(kh) K_k K_{-k} \right)$$
  waarbij $E_k$ en $K_k$ de Fourier-amplitudes zijn van de rek- en buigspanningen, respectievelijk.
• Effectieve Moduli: De auteurs leiden effectieve rek- ($A$) en buigmoduli ($B$) af die afhankelijk zijn van de dimensieloze parameter $\eta = kh$ (de verhouding tussen de plaatdikte en de golflengte van de vervorming).
  • Wanneer $\eta \to 0$ (de dunne-plaatlimiet), geldt $A \to 2$ en $B \to 1/6$, waarmee de klassieke dunne-plaatresultaten worden teruggevonden.
  • Wanneer $\eta$ toeneemt (dikker plaat of kortere golflengtes), nemen zowel $A$ als $B$ af, wat aangeeft dat de plaat verzacht in zowel rek als buiging.
• Constante Verhouding van Moduli: Een opvallend resultaat is dat de verhouding tussen de bijdragen van rek en buiging, gekwantificeerd door het "rekgetal" $S = A/12B$, ongeveer constant blijft voor $\eta \lesssim 1$. Specifiek: $S - 1 \approx \eta^2/30$. Dit impliceert dat hoewel de absolute stijfheid van de plaat afneemt met de dikte, de relatieve energetische kosten van rek versus buiging bijna identiek blijven aan die van een dunne plaat.
• Validatie: De asymptotische oplossing vertoont een uitstekende overeenkomst met de eindige-elementen-simulaties, wat de geldigheid van de afgeleide effectieve moduli en de aanwezigheid van discontinuïteiten over de middenlijn bevestigt.

Betekenis
Het artikel biedt een theoretische rechtvaardiging voor de empirische observatie dat dunne-plaattheorieën vaak de mechanica van dikke platen accuraat beschrijven, zoals biologische celbladen. De auteurs stellen dat de verrassende nauwkeurigheid van deze theorieën voortvloeit uit het feit dat de scheiding van rek- en buigenergie persisteert, zelfs in de limiet van kleine vervormingen van een dikke plaat. Bovendien blijft de relatieve belangrijkheid van rek versus buiging behouden tot op lengteschalen die vergelijkbaar zijn met de dikte van de plaat.

Concluderend kan het mechanische gedrag van dikke elastische platen worden gevangen door dunne-plaattheorieën, mits de externe krachten worden geschaald om rekening te houden met het verzachten van de effectieve moduli. Deze bevinding is bijzonder relevant voor de morfofelasticiteit en de studie van biologische systemen (bijv. Volvox celbladen) waar het materiaal dik is en door groei geïnduceerde vervormingen optreden. Het werk suggereert dat het conceptuele kader van intrinsieke rek en buiging, breed gebruikt in dunne-schiltheorieën, ook toepasbaar blijft op dikke biologische structuren.