Stretching and bending of (really) thick elastic plates

技术摘要：厚（实）弹性板的拉伸与弯曲

问题陈述
弹性板的机械能传统上被理解为可分解为截然不同的拉伸和弯曲分量。对于渐近薄板（$h = O(\varepsilon)$，其中 $\varepsilon \ll 1$），这种分离是一个严谨的结果，它是基于尺度假设推导出的：中面拉伸极小（$E = O(\varepsilon)$）且弯曲有限（$K = O(1)$）。然而，这种“薄板图景”对于厚板（即使其厚度 $h = O(1)$ 且基尔霍夫假设——即法线保持法向——失效时）也能提供令人惊讶的准确近似。本研究旨在解决的核心难题是：对于“真正的厚板”，在小变形情况下，这种能量在拉伸模式与弯曲模式之间的分离是否在形式上仍然成立？如果是，其有效模量将如何表现？

方法论
作者通过求解一个二维、不可压缩、新胡克尔（neo-Hookean）厚度为 $h = O(1)$ 的板的弹性方程来研究这一问题。他们对板的中面施加了渐近小变形，定义为拉伸应变 $E = O(\varepsilon)$ 和曲率 $K = O(\varepsilon)$。

1. 运动学与本构关系： 该板使用笛卡尔坐标 $(x, z)$ 进行描述，其中 $z \in [-h/2, h/2]$。变形由中面的位移以及通过厚度方向的位移所构成的变形梯度张量 $F$ 表征。材料被建模为不可压缩的新胡克尔材料，其弹性能量密度正比于 $\text{tr}(F^T F) - 2$。
2. 渐近展开： 作者利用小参数 $\varepsilon$ 进行正则摄动展开。他们在 $\varepsilon$ 的幂次项中对位移场（$\tilde{u}, \tilde{w}$）和压力 $p$ 进行展开。
3. 边界条件： 至关重要的是，作者对顶层和底层表面（$z = \pm h/2$）施加了无力边界条件，并在中面（$z=0$）处规定了变形。不同于假设中面连续性的薄板理论，作者明确允许解在中面（$z=0$）处存在不连续性，以满足表面的边界条件，这承认了基尔霍夫假设在厚板中并不成立。
4. 求解策略： 利用傅里叶模态分解（$e^{ikx}$）求解线性化的一阶问题。这使得偏微分方程简化为关于 $z$ 的常微分方程组，并进行了解析求解。结果通过 FEniCSx 进行的有限元模拟得到了验证。

核心贡献与结果

• 能量的形式分离： 计算表明，对于发生小中面变形的厚板（$h = O(1)$），总弹性能量密度仍然可以分离为拉伸项和弯曲项。能量形式如下：
  $$E \propto \sum_k \left( h A(kh) E_k E_{-k} + h^3 B(kh) K_k K_{-k} \right)$$
  其中 $E_k$ 和 $K_k$ 分别是拉伸和弯曲应变的傅里叶振幅。
• 有效模量： 作者推导出了依赖于无量纲参数 $\eta = kh$（即变形波长与板厚度的比值）的有效拉伸（$A$）和弯曲（$B$）模量。
  • 当 $\eta \to 0$（薄板极限）时，$A \to 2$ 且 $B \to 1/6$，恢复了经典的薄板结果。
  • 当 $\eta$ 增大时（板更厚或波长更短），$A$ 和 $B$ 均会减小，表明板在拉伸和弯曲方面都表现出软化特性。
• 恒定的模量比例： 一个引人注目的结果是，通过“拉伸数” $S = A/12B$ 量化的拉伸与弯曲贡献之比在 $\eta \lesssim 1$ 时保持近似常数。具体而言，$S - 1 \approx \eta^2/30$。这意味着，虽然板的绝对刚度随厚度增加而降低，但拉伸相对于弯曲的相对能量成本与薄板几乎完全相同。
• 验证： 渐近解与有限元模拟表现出极佳的一致性，证实了所推导的有效模量的有效性以及中面处存在不连续性的事实。

意义
本文为“薄板理论通常能准确描述厚板力学行为”这一经验观察提供了理论依据。作者认为，这些理论之所以具有惊人的准确性，是因为即使在厚板的小变形情况下，拉伸与弯曲能量的分离现象依然存在。此外，在波长与板厚度相当的尺度范围内，拉伸与弯曲的相对重要性也保持不变。

因此，通过对外部作用力进行重新标度以补偿有效模量的软化，厚弹性板的力学行为可以由薄板理论来捕捉。这一发现对于形态弹性学（morphoelasticity）以及生物系统（例如 Volux 细胞片）的研究尤为重要，因为这些系统的材料是厚实的，且存在生长诱导的变形。这项工作表明，广泛应用于薄壳理论的内在拉伸与弯曲的概念框架，对于厚生物结构仍然适用。