[论文解读] Periodically wrinkled plate of the Föppl-von Kármán type

Igor Velčić|arXiv (Cornell University)|Apr 4, 2011

Structural Analysis and Optimization参考文献 11被引用 3

一句话总结

本文通过三维非线性弹性理论的 $\Gamma$-收敛，从 Föppl-von Kármán 类型推导出周期性褶皱板模型。假设板厚为 $h^2$，中面具有周期性微结构 $h^2\theta(x_1/h, x_2/h)$，结果表明应变能按 $h^8$ 缩放，从而得到一个耦合弯曲（$\nabla^2 v$）与拉伸（$\mathop{\rm sym}\nabla{\bf u} + \frac{1}{2}\nabla v \otimes \nabla v$）能的极限模型，准确捕捉褶皱薄板的力学行为。

ABSTRACT

In this paper we derive, by means of $Γ$-convergence, the periodically wrinkled plate model starting from three dimensional nonlinear elasticity. We assume that the thickness of the plate is $h^2$ and that the mid-surface of the plate is given by $(x_1,x_2) o (x_1,x_2,h^2θ(\per))$, where $θ$ is $[0,1]^2$ periodic function. We also assume that the strain energy of the plate has the order $h^8=(h^2)^4$, which corresponds to the Föppl-von Kármán model in the case of the ordinary plate. The obtained model mixes the bending part of the energy with the stretching part.

研究动机与目标

• 从三维非线性弹性理论出发，严格推导出周期性褶皱薄板的二维板模型。
• 分析厚度为 $h^2$、中面变形为 $h^2\theta(x_1/h, x_2/h)$ 的板的渐近行为，其中 $\theta$ 是 $[0,1]^2$-周期函数。
• 建立一个耦合弯曲与拉伸贡献的极限能量泛函，将经典 Föppl-von Kármán 模型推广至褶皱构型。
• 证明所得模型依赖于预变形 $\theta$ 本身，而不仅其导数，凸显初始几何缺陷在力学响应中的关键作用。

提出的方法

• 对三维非线性弹性能量进行 $\Gamma$-收敛分析，能量按 $h^8$ 缩放，并令 $h \to 0$。
• 利用多尺度收敛方法处理中面中周期性微结构 $\theta(x_1/h, x_2/h)$ 的影响。
• 将参考构型定义为厚度为 $h^2$、中面为 $\hat{\Omega}^h = \Theta^h(\omega \times (-h^2/2, h^2/2))$ 的壳体，其中 $\Theta^h$ 包含周期性形状函数 $\theta$。
• 引入位移假设 $\tilde{\bf y}^h = (\bar{{\bf R}}^h)^T {\bf y}^h - {\bf c}^h$，并将变形分解为宏观与微观分量。
• 推导出极限能量泛函 $J_0^L({\bf u}, v, \bar{{\bf R}})$，其同时包含弯曲项（$\nabla^2 v$）与拉伸项（$\mathop{\rm sym}\nabla{\bf u} + \frac{1}{2}\nabla v \otimes \nabla v$）。
• 在空间 $W^{1,2}(\omega;\mathbb{R}^2) \times W^{2,2}(\omega) \times \mathop{\rm SO}(3)$ 中证明极限泛函的极小化解存在。

实验结果

研究问题

• RQ1 在薄板中存在周期性微屈曲（褶皱）时，其从三维弹性理论推导出的有效二维能量模型如何变化？
• RQ2 为获得一个非平凡的极限模型，该模型能耦合褶皱板中的弯曲与拉伸能，应变能的缩放比例（$h^8$）应如何确定？
• RQ3 为何极限能量泛函依赖于预变形 $\theta$ 本身，而非仅其导数，这与标准板模型形成对比？
• RQ4 如何通过 $\Gamma$-收敛与多尺度收敛方法，为厚度为 $h^2$ 且中面具有周期性变形的板推导出一致的极限模型？
• RQ5 极限能量中旋转 $\bar{{\bf R}}$ 的作用是什么？为何当横向载荷 $f_3 \neq 0$ 时，有 $\bar{{\bf R}}_{33} = \pm 1$？

主要发现

• 极限能量泛函 $J_0^L$ 耦合了弯曲项（$\nabla^2 v$）与拉伸项（$\mathop{\rm sym}\nabla{\bf u} + \frac{1}{2}\nabla v \otimes \nabla v$），反映了褶皱与整体变形之间的力学相互作用。
• 应变能按 $h^8 = (h^2)^4$ 缩放，对应于厚度为 $h^2$ 的板的 Föppl-von Kármán 模型阶次，验证了该缩放选择的合理性。
• 极限泛函的极小化解存在于空间 $W^{1,2}(\omega;\mathbb{R}^2) \times W^{2,2}(\omega) \times \mathop{\rm SO}(3)$ 中，确保了该模型的数学适定性。
• 当横向载荷 $f_3$ 非零时，极限旋转满足 $\bar{{\bf R}}_{33} = 1$ 或 $-1$，表明板的取向在极限下趋于刚体旋转。
• 该模型显式依赖于形状函数 $\theta$ 本身，而不仅其导数，说明在褶皱构型中，初始几何形状具有显著的力学影响。
• 通过 $\Gamma$-收敛与多尺度收敛的推导，严格证明了该模型是三维非线性弹性理论在 $h \to 0$ 极限下的一致极限，确立了其理论一致性。