第十章：形态弹性板（Morphoelastic Plates）

阅读笔记

1. 概述

本章由 Alain Goriely 撰写，收录于《The Mathematics and Mechanics of Biological Growth》（Interdisciplinary Applied Mathematics 45）一书，是该著作关于生物生长数学力学研究的第十章。本章系统性地将前面章节中发展的轴对称壳与膜理论推广到更为一般的板壳结构，建立了生长与重塑的完整板壳理论框架。

本章的核心目标是建立一套适用于生长弹性板的数学理论。在生物系统中，薄板结构广泛存在——从植物叶片到细胞外基质，从昆虫角质层到软组织表层，因此发展一套统一的薄板生长理论具有重要的生物学意义。本章首先回顾了经典弹性板理论的基本框架，重点介绍了 Föppl-von Kármán 方程，随后将这些方程推广以纳入生长的效应，并探讨了生长导致的板结构失稳与屈曲现象。

从历史脉络来看，弹性壳理论的研究始于 Chladni、Germain 和 Kirchhoff 关于振动板的实验研究（参考文献 [203, 437, 698]），后经 Love（参考文献 [805]）发展建立了小应变线性弹性壳的一致性一般理论。Cosserat 兄弟（参考文献 [235]）首先推导了以应力和偶应力张量表示的弹性壳控制方程，随后 Synge 和 Chien（参考文献 [1204]）提出了无需坐标的理论框架。此后的研究者们，包括 [501, 719, 932, 966] 等，提出了众多替代性的推导和公式化方法。

在生长板的理论方面，已有若干研究尝试将板理论推广以包含生长效应，包括 [93, 323, 705, 851, 852] 等工作。本章主要采用 Ben Amar 及合作者（参考文献 [275, 277, 923]）提出的 Föppl-von Kármán 方程生长推广形式。

2. 弹性板理论基础

2.1 壳与板的基本概念

在进入板理论细节之前，有必要明确壳与板的定义与区分。壳（shell）是二维弹性体，其参考形状不必是平的；而板（plate）是壳的一种特例，其初始参考形状是平的。更为一般地，形态弹性壳（morphoelastic shell）是能够随时间重塑和生长的弹性壳。这些理想化对象是许多物理、工程和生物系统的合适模型。

从几何角度而言，三维弹性体在一定条件下可以简化为弹性板或弹性壳。这种简化基于"薄"的假设——当物体的厚度远小于其横向特征尺度时，三维弹性理论可以在厚度方向上降维，得到等价的三维曲面（板的中面或壳的参考面）上的二维方程。

2.2 Kirchhoff-Love 假设

经典板理论建立在 Kirchhoff-Love 假设之上，这些假设大大简化了三维弹性问题：

小厚度假设：板的厚度 H 远小于其横向尺寸 L，即 H/L << 1。这意味着可以将板视为二维曲面而非三维实体。
小位移假设：变形足够小，可以完全由中面的垂直位移 ζ(X, Y ) 表征。
小曲率假设：板的曲率保持小量。
无剪切假设：厚度方向无剪应变，即变形后法线保持与中面垂直（只是略微偏离 Z 方向）。
无正应力假设：Cauchy 应力张量的 Z 分量趋于零（T_zz ≈ 0）。

这些假设使得分析大幅简化，同时保留了弹性板行为的物理本质。

2.3 Airy 势能

在 Kirchhoff-Love 假设下二维应力状态可以用单个标量函数——Airy 势能（Air function）χ(X, Y ) ——来表示：

T_xx = ∂²χ/∂Y²
T_yy = ∂²χ/∂X²
T_xy = -∂²χ/∂X∂Y

这种表示自动满足二维应力平衡方程。Airy 势能是一种数学工具，它将两个应力分量打包成一个标量函数，从而简化了方程系统。

3. 曲面几何基础：平均曲率与高斯曲率

3.1 参数化曲面

本章从几何角度建立曲面的数学描述。考虑一个可定向的参数化曲面 Σ，由位置向量方程定义：

x = x(ξ¹, ξ²) ∈ ℝ³, (ξ¹, ξ²) ∈ D ⊆ ℝ²

假设 x 至少是 C² 类的，且切向量 r_i = ∂x/∂ξⁱ 在 D 上线性独立。由于 Σ 可定向，可以定义单位法向量：

n = (r₁ × r₂)/|r₁ × r₂|

{r₁, r₂, n} 形成一个基，但该基不一定是正交的。

3.2 第一基本形式（度规张量）

曲面的面积元素涉及度规张量 M，其分量为 M_ij = r_i · r_j。曲面面积由以下积分给出：

A = ∫_Σ dS = ∬_D det(M) dξ¹ dξ²

曲面上一条曲线 ρ 的弧长为：

L = ∫_ρ ds = ∫_I √(M_ij ξ̇ⁱ ξ̇ʲ) dt

第一基本形式为 ds² = M_ij dξⁱ dξʲ，它描述了曲面的内禀几何性质。

3.3 第二基本形式与曲率

法向曲率是曲面的内禀性质。对于曲面上一点 P 处的曲线，其曲率向量分解为法向曲率向量 k_n 和测地曲率向量 k_g：

k = k_n + k_g

其中 k_n 沿法向量方向，k_g 沿曲面切向。测地曲率对于测地线恒为零。

法向曲率通过第二基本形式给出：

k_n = K_ij ξ̇ⁱ ξ̇ʲ

其中 K_ij = -n · (∂r_j/∂ξⁱ) 是外在曲率张量（第二基本形式的系数）。

3.4 主曲率、平均曲率与高斯曲率

主曲率（principal curvatures）k₁ 和 k₂ 是主曲率矩阵 L = M⁻¹ K 的特征值，表示法向曲率的极值。

平均曲率（mean curvature）为：

K_M = (1/2) tr(L) = (k₁ + k₂)/2

高斯曲率（Gaussian curvature）为：

K_G = det(L) = k₁ × k₂

3.5 Gauss 的 Theodora Egregium 定理

高斯曲率是曲面的内禀量——它只依赖于度规（第一基本形式），而不依赖于法向量。这是 Gauss 的 Theodora Egregium（奇妙定理）的核心内容（参考文献 [169, 496]）。这一重要结果表明，曲面的弯曲程度完全由其内禀几何决定，与外部嵌入方式无关。

根据高斯曲率的符号，曲面可分为： - 椭圆型（elliptic）：K_G > 0，如球面 - 双曲型（hyperbolic）：K_G < 0，如马鞍面 - 抛物型（parabolic）：K_G = 0，如平面

极小曲面（minimal surface）是平均曲率 K_M 在所有点上恒为零的曲面，这类曲面在数学和物理学中具有重要意义（参考文献 [1001]）。

4. Föppl-von Kármán 方程

4.1 方程的推导背景

Föppl-von Kármán 方程是经典板理论的核心方程。它是通过对三维弹性连续体的控制方程进行系统的量纲展开而得到的，展开参数为小量 H/L（厚度与特征长度之比）。

在 Kirchhoff-Love 假设下，变形完全由中面垂直位移 ζ(X, Y ) 表征，结合 Airy 势能 χ(X, Y )，可以得到两个耦合的偏微分方程——Föppl-von Kármán 方程。

4.2 控制方程

生长弹性板的 Föppl-von Kármán 方程为：

方程一（弯曲方程）：

D∇⁴ζ - ∇²_K M - 2H[χ, ζ] = P

方程二（Air 函数方程）：

∇²χ + E[ζ, ζ] - K_G = 0

其中： - E 是杨氏模量（Young's modulus） - D = EH³/12 是板的弯曲刚度（bending rigidity） - P 是作用于板法向的压强 - ∇⁴ 是双调和算子 - ∇² 是拉普拉斯算子 - H 是厚度（注意：与平均曲率 K_M 区分使用） - [f, g] 是非线性双线性括号，定义为：

[f, g] = (1/2)(∂²f/∂X²)(∂²g/∂Y²) + (1/2)(∂²f/∂Y²)(∂²g/∂X²) - (∂²f/∂X∂Y)(∂²g/∂X∂Y)

4.3 方程的物理意义

第一个方程描述了板的弯曲行为：左边第一项 D∇⁴ζ 代表弯曲弹性恢复力，第二项 -∇²_K M 代表平均曲率变化的贡献，第三项 -2H[χ, ζ] 代表面内应力与弯曲变形的耦合（几何非线性），右边 P 是外加载荷。

第二个方程是相容方程：左边第一项 ∇²χ 代表面内应力的弹性贡献，第二项 E[ζ, ζ] 来自几何非线性（由弯曲引起的膜应变），第三项 -K_G 来自高斯曲率变化。在无生长情况下 K_G = 0，方程简化为经典的 Föppl-von Kármán 相容方程。

4.4 经典 Föppl-von Kármán 方程

在没有生长的情况下（即 K_M = K_G = 0），上述方程退化为经典的 Föppl-von Kármán 方程，这些方程最初由 Föppl 和 von Kármán 正式推导（参考文献 [372, 1286]）。此后，许多研究提出了严格的推导方法，包括 [74, 209, 210, 394, 788, 934] 等。

经典 Föppl-von Kármán 方程成功描述了薄板的弹性弯曲与面内应力的耦合，是结构力学中的基石方程。

5. 生长弹性板的推广

5.1 生长的几何解释

本章采用 Ben Amar 及合作者提出的推广框架，将生长效应纳入 Föppl-von Kármán 方程。关键创新在于将生长解释为内禀曲率的来源。

具体而言，生长通过两个几何参数进入方程： - 平均曲率 K_M：描述中面的平均弯曲程度 - 高斯曲率 K_G：描述高斯曲率（内禀曲率）

这两个曲率可以从生长应变直接计算得到。这意味着生长改变了板的参考构型——板"记得"它曾经生长的方式，并据此产生内禀应力。

5.2 生长的物理机制

当板发生生长时，新生的材料会"记住"其形成时的几何构型。如果生长在板内是非均匀的（即不同点的生长量不同），就会在板内产生残余应力。这种残余应力来源于几何约束——生长区域想要扩张或收缩，但被周围的材料所约束。

生长导致的内禀曲率可以理解为：生长过程使板的中面发生了永久的几何改变。当生长发生时，板的参考状态（无应力状态）与当前的中面构型不再一致——板"想要"处于一个弯曲的构型，但弹性恢复力试图将其拉回平的构型。这种竞争产生了内禀应力分布。

5.3 生长的分类

根据生长方式的不同，可以分为：

各向同性生长：板在所有方向均匀生长，若生长后板的面积增大但厚度不变，则板保持平直（K_G = 0）。
各向异性生长：板在不同方向上的生长率不同。例如：
径向生长慢、周向生长快 → 正高斯曲率（碗状）
边缘生长快、中心生长慢 → 负高斯曲率（波浪边缘）

6. 生长与屈曲

6.1 生长诱导的失稳

当生长产生的内禀应力超过板的弹性承载能力时，板会发生失稳（屈曲）。这是一个典型的力学生物学问题：生物体利用生长诱导的应力来实现特定的形态发生。

本章给出了几个典型的生物示例：

示例一：藻类 Acetabularia acetabulum

Acetabularia 是一种单细胞绿藻，具有伞状的"帽"结构。在形态发生过程中，其形状发生变化（图 10.3A）。研究表明，通过考虑盘状结构在不同方向上的差异生长，可以重现这种平衡形状（图 10.3B，参考文献 [275]）。

示例二：植物叶片

叶片发育中的波纹或平整控制与高斯曲率直接相关。差分生长导致叶片产生内禀曲率： - 叶片边缘与中央的生长差异产生高斯曲率 - 正高斯曲率 → 碗状叶片 - 负高斯曲率 → 波纹边缘

这一机制解释了植物叶片如何通过调控局部生长来实现其特定形状（参考文献 [830, 937, 1097, 1141, 1142]）。

示例三：草叶的褶皱

平行夹持的生长条带的失稳可以解释草叶的褶皱现象（图 10.5）。生长在条带中产生压缩应力，当压缩应力超过临界值时，条带发生屈曲形成Rippling（波纹）模式。这一模型很好地解释了植物叶片的皱褶形态（参考文献 [277]）。

6.2 负高斯曲率的产生

当盘的边缘生长快于中央部分时，会产生负高斯曲率结构。负高斯曲率的特点是主曲率符号相反（如马鞍面）。为释放残余应力，板会形成波浪边缘的形状（图 10.4，参考文献 [1032, 1033]）。

这一机制在自然界中广泛存在：从植物叶片边缘的波浪到花朵的皱褶，都可能是差分生长导致负高斯曲率的结果。

6.3 屈曲的数学描述

屈曲本质上是一个分岔问题。当外加载荷或内禀生长应变达到某个临界值时，板的平直构型失去稳定性，出现新的非平凡解（屈曲构型）。

在 Föppl-von Kármán 方程框架下，屈曲表现为方程的解从唯一的平直解 ζ = 0 分岔到非零的弯曲解。临界条件由线性稳定性分析确定：当作用在板上的有效压缩应力达到屈曲临界值时，系统发生分岔。

7. Cosserat 板与 Kirchhoff 板理论

7.1 理论扩展

Föppl-von Kármán 方程是弹性板的简化模型。放松一些模型假设会得到更复杂的模型。

Cosserat 板理论和 Kirchhoff 板理论约束三维体的可能变形以适应二维几何（参考文献 [33]）。这种方法可以探索材料的完整非线性行为，不限于薄体情况。这些理论通过直接使用后续章节将讨论的三维体生长理论，可以推广到生长体的情况（参考文献 [851, 852]）。

7.2 壳的一般理论

当参考构型是壳而非板时，可以使用基于微分几何的一般理论来推导生长壳的基本方程（参考文献 [1110]）。这一一般理论特别有意义，因为它明确展示了曲面内禀几何与生长内禀几何之间的耦合（将在第 12.5 节讨论）。

8. 关键公式总结

8.1 几何量

度规张量：M_ij = r_i · r_j
外在曲率张量：K_ij = -n · (∂r_j/∂ξⁱ)
主曲率矩阵：L = M⁻¹ K
平均曲率：K_M = (1/2) tr(L) = (k₁ + k₂)/2
高斯曲率：K_G = det(L) = k₁ × k₂

8.2 Föppl-von Kármán 方程

弯曲方程：

D∇⁴ζ - ∇²_K M - 2H[χ, ζ] = P

相容方程：

∇²χ + E[ζ, ζ] - K_G = 0

非线性括号：

[f, g] = (1/2)(∂²f/∂X²)(∂²g/∂Y²) + (1/2)(∂²f/∂Y²)(∂²g/∂X²) - (∂²f/∂X∂Y)(∂²g/∂X∂Y)

弯曲刚度：

D = EH³/12

8.3 Airy 势能与应力分量

T_xx = ∂²χ/∂Y²
T_yy = ∂²χ/∂X²
T_xy = -∂²χ/∂X∂Y

9. 生物学应用与意义

9.1 形态发生的力学控制

本章建立的形态弹性板理论揭示了生物形态发生中力学的核心作用。通过生长实现形态变化是生物界的基本策略：从细胞分裂到器官发育，从叶片展开到花朵绽放，都涉及生长介导的力学过程。

关键认识是：生长本身可以产生力。当生长受到几何约束时，生长过程会产生内禀应力。这些应力进而导致屈曲、褶皱、波纹等丰富的形态发生模式。

9.2 形态-功能关系

生物结构的功能往往与其形态密切相关。植物叶片需要平整以最大化光合作用面积，同时又需要一定的机械强度和水分保持能力。波纹边缘可以增加叶片的机械稳定性，同时减少风阻。通过差分生长调控高斯曲率，植物可以实现这些多目标优化。

9.3 生长与应力的耦合

形态弹性板理论的核心是生长与弹性的耦合。这种耦合是双向的： - 生长 → 内禀应力 → 弹性变形 - 弹性应力 → 生长速率的调控

后者的机制被称为机械敏感生长（mechanosensitive growth），是力学生物学的重要研究方向。

10. 与前序章节的联系

本章是前面关于轴对称壳和膜理论的发展与深化。在第 5-9 章中发展的概念——生长应变、应力-生长耦合、形态稳定性——在板理论框架下得到了更为一般的表述。

本章也起到了从二维板壳理论向三维体生长理论过渡的桥梁作用。本章末尾明确指出，后续章节将直接使用三维体生长理论来处理生长板问题，这构成了本书第四部分"体生长：三维理论"的基础。

11. 参考文献要点

本章引用的关键文献包括：

Ben Amar 及合作者 ([275, 277, 923])：Föppl-von Kármán 方程的生长推广
Föppl, von Kármán ([372, 1286])：经典 Föppl-von Kármán 方程
Dumais 及合作者 ([1032, 1033])：负高斯曲率与差分生长
Goriely 及合作者 ([830, 937, 1097, 1141, 1142])：植物叶片的形态发生
Conti 及 Ciarletta ([851, 852])：生长的 Cosserat 板理论
Efrati 及 Klein ([1110])：生长壳的一般微分几何理论

12. 阅读思考

高斯曲率的内禀性：Gauss 的 Theorema Egregium 表明高斯曲率是内禀量，仅由度规决定。这意味着无论板如何嵌入三维空间，其高斯曲率变化都完全由生长应变决定。这一事实对理解生长的几何约束有何启示？
屈曲的适应性意义：生物体为何经常利用屈曲来实现形态？屈曲作为一种形态发生机制有什么优势？
多尺度耦合：从细胞水平的生长到器官水平的形态，涉及多个尺度的耦合。如何在形态弹性板理论框架下纳入这些多尺度效应？
生长的时间依赖性：本章主要关注准静态平衡，没有深入讨论生长的动力学过程。如何将生长的速率方程纳入 Föppl-von Kármán 框架？

13. 结论

第十章系统地建立了生长弹性板的数学力学框架。通过将经典 Föppl-von Kármán 方程推广以纳入生长效应，本章揭示了生长作为内禀曲率来源的重要概念。高斯曲率和平均曲率作为生长的几何表征，将生长过程与弹性响应统一在同一个数学框架中。

生长诱导的屈曲失稳是形态发生的核心机制之一。通过差分生长控制局部曲率，生物体可以实现复杂而精确的形态构建。从藻类的伞状结构到植物叶片的波纹边缘，形态弹性板理论为理解这些现象提供了统一的力学基础。

本章建立的框架不仅适用于板结构，也为后续章节讨论一般三维生长理论奠定了基础，是连接二维壳膜理论与三维体生长理论的关键桥梁。

阅读日期：2026年5月9日 章节字数：约 4500 字