第七章：增生生长（Accretive Growth）

阅读笔记

---

一、章节概述

本章探讨了可建模为二维物体的生物结构生长问题，即横向长度尺度远小于其他两个典型长度尺度的三维物体。在建立有效的表面生长模型时，区分硬表面和软表面至关重要。贝壳、角、鹿角和牙齿可建模为硬表面，这些结构无法轻易变形；而叶片、花粉管和花朵则可建模为软表面，存在一定程度的弹性变形。本章以贝壳为典型范例，系统阐述了增生生长的几何学与力学机制。

---

二、硬表面与软表面的基本区别

2.1 硬表面增生生长

硬表面结构（如贝壳、角、鹿角、牙齿）的生长通过在结构边界局部添加物质来实现，因为物质本体无法通过变形来添加质量。从连续介质力学的角度来看，表面生长带来了一个特别有趣的挑战：由于新材料点不断沉积在物体表面上，所有物理量定义所在的域也随时间演化，因此不存在固定的参考配置。这使得将增生生长纳入经典连续介质力学需要特别小心，相关理论研究和连续介质平衡定律已有多位学者提出。

2.2 软表面面积生长

软表面结构（如叶片、花粉管、花朵）的生长通常分布在整个区域，形状是局部物质添加与表面弯曲和拉伸特性之间竞争的结果。这种生长类型被称为"面积生长"（area growth），区别于硬表面的局部增生机制。

---

三、内禀增生生长理论（Intrinsic Accretive Growth）

3.1 基本数学框架

本章采用局部构造观点来研究硬体表面的生长问题，遵循Moulton及其合作者建立的理论。考虑从生成曲线开始的表面生长：

• 在时间t=0时，假设曲线为二次连续可微曲线R(S) = r(S, 0) ∈ ℝ³，由弧长S ∈ [0, L]定义
• 在曲线上的每一点S定义一个生长向量q(S, 0)，表示生长过程中材料点的局部演化
• 生长向量场是材料点在S位置的初始生长速度
• 在任意时间t，指定生长场q(S, t)，定义材料点S在时间t的生长速度

生长速度向量描述的是增生生长过程，该曲线可被视为细胞轨迹——即细胞随时间沉积新材料的痕迹。由于生长速度的方向和幅度由局部过程决定，因此假设速度向量是局部坐标系以及物理、化学或生物场的函数是合理的，例如增生率、形态发生素梯度、温度、pH值等。

3.2 活动标架与Darboux矩阵

沿曲线r(S, t)建立一个通用标架D = (d₁, d₂, d₃)，其中d₃是r(S, t)的单位切向量。标架满足正交归一条件：

DᵀD = DDᵀ = I

这意味着D是一个正交矩阵。沿曲线的时间变化和空间变化分别由自旋矩阵W和Darboux矩阵U描述：

• Darboux矩阵U：描述标架在给定时间t沿曲线的旋转
• 自旋矩阵W：描述标架在点S处随时间的旋转

这两个矩阵都是反对称矩阵，其元素与曲率和挠率相关。

3.3 生长向量的相容性条件

假设曲线r至少在t和S上二次可微，二阶混合偏导数的对称性提供了两组相容性条件：

∂²r/∂S∂t = ∂²r/∂t∂S
∂²D/∂S∂t = ∂²D/∂t∂S

利用正交归一条件DᵀD = I，经过简化后得到：

q' - v̇ = Wv - Uq
Ẇ - U̇ = UW - WU

这些方程对于给定的生长向量q构成了七个因变量{u, w, λ}的六个非线性一阶偏微分方程组。额外的自由度来自于通用标架选择的额外自由度。

3.4 Frenet标架的特殊情况

如果假设曲线至少三次可微，则可使用Frenet标架，{d₁, d₂, d₃} = {ν, β, τ}分别表示法向量、副法向量和切向量。此时Darboux矩阵取简单形式：

U = [0  -τ  κ
     τ   0   0
    -κ   0   0]

其中κ ≥ 0和τ分别表示几何曲率和挠率。这使得系统简化为六个方程六个未知数，可以针对初始轮廓和生长场求解。

---

四、简单示例

4.1 圆柱、圆锥和环面

4.1.1 圆柱生长

取生成曲线C₀为半径为1的圆，即r(S, 0) = (cos(S), sin(S), 0)。选择沿副法线方向的恒定生长速度场q = (0, 1, 0)。此时相容性方程和标架演化都是平凡的，曲面由积分得到：

r(S, t) = ∫₀ᵗ q dt + r(S, 0) = (cos(S), sin(S), t)

这产生了圆柱面的参数化。

4.1.2 扭曲圆柱

如果添加沿切向和法向的恒定速度分量，取生长速度场q = (0, c₂, c₃)，其中c₂c₃ ≠ 0，得到：

r(S, t) = (cos(S + c₃t), sin(S + c₃t), c₂t)

这是扭曲圆柱的参数化。这种差异可通过细胞轨迹运动来理解：在(a)中，细胞直接沿圆环的副法线方向沉积新材料；而在(b)中，细胞有额外的切向分量。

4.1.3 扭曲圆锥

进一步添加法向速度分量，取q = (c₁, c₂, c₃)且c₁c₂c₃ ≠ 0。此时由于λ̇ = -c₁，曲面参数化为：

r(S, t) = ((1 - c₁t)cos(S + c₃t), (1 - c₁t)sin(S + c₃t), c₂t)

所得曲面是具有扭曲细胞轨迹的圆锥。

4.2 环形生长

如果局部生长向量不是常数，可获得更有趣的形状。再次取单位圆作为初始曲线，通过取q₂为生长轴平面中的线性函数来实现卷曲：

q₂ = b₁ + b₂cos(S)

其中b₁和b₂为常数。选择q₁ = q₃ = 0并积分控制方程，可得到环面参数化公式。

---

五、形状平面性与形状不变性

5.1 形状平面性条件

考虑任意平面曲线r(S, t)。假设曲线在演化过程中保持平面的条件是其挠率必须恒等于零，即u₃ = u̇₃ = 0。使用这个条件可得到q₂的方程，其通解为：

q₂ = b₁(t) + b₂(t)x(S, t) + b₃(t)y(S, t)

这个形式意味着一个重要的普适特征：如果生成曲线保持平面，则存在一个沿着向量b₂X̂ + b₃Ŷ的区分轴，使得生成曲线上的副法线生长在投影到这个轴上时是一个线性函数。这个轴被称为生长轴。

5.2 形状不变性条件

进一步限制生长速度场，使得平面生成曲线的形状保持相似（即除了整体膨胀外与原始曲线相同）。这意味着曲线只能进行平移、旋转和膨胀。形状不变性要求：

• 在任意时间t，|r'| = λ = λ(t)
• 曲率κ(S, t)是按膨胀因子λ缩放的初始曲率：κ(S, t) = λ⁻¹κ(S, 0)

在全局坐标系中，解的形式为：

q = c₁(t)r + c₂(t)r⊥ + a(t)

其中r⊥是垂直于位置向量r且具有相同大小的向量。第一项描述沿径向的生长速度，第二项给出沿圆周方向的速度，第三项描述刚性平移。

---

六、自相似性

6.1 自相似性的定义

许多贝壳的一个显著特征是尺度不变性，也称为自相似性。自相似性是物体的全局性质，可定义为：空间中连续界标点之间的距离比例恒定。例如，对数螺旋贝壳中连续轮之间的距离比例是常数。

6.2 局部自相似的条件

定义局部自相似性的问题是：哪些局部生长速度场能确保全局自相似性？尺度不变性的关键特性是：对于每个材料点，生长向量的局部方向不随时间变化。在q₃ = 0的情况下，这表达为：

∂t arctan(q₂/q₁) = 0

为了在时间上强制不变性，还要求生长向量的大小与生成曲线大小（即圆截面半径）的比值在时间上不变。这促使了生长定律的提出：

q₁ = -λ(t)c
q₂ = λ(t)(b₁ + b₂cos(S))
q₃ = 0

这个生长定律产生的结构类似于图7.4A，但区别在于其生长率的大小也是尺度不变的。

---

七、贝壳生长的应用

7.1 背景知识

贝壳因其美丽和简单而几个世纪以来一直吸引着人类的想象力。贝壳曾被用作货币、珠宝、装饰品、工具，也是众多收藏家的最爱。大多数贝壳是软体动物门腹足纲的硬化矿化外壳。贝壳是由矿物和生物聚合物组成的微层状复合生物陶瓷，排列成棱柱状方解石晶体或砖墙状文石晶体层。这些层的规则和复杂排列占动物总质量的98%，赋予了贝壳超过矿物质本身的强度。

大多数现代软体动物的贝壳是外套膜产生的重复分泌过程的产物。外套膜是软体动物的一层薄软组织，从贝壳中延伸出来包裹其边缘。在生长过程中，外套膜在分泌角质层的同时向前移动到矿化贝壳边缘之外。矿化发生时，角质层固定在外壳表面。贝壳的形状，即围绕较软外套膜的较硬角质层，在生长过程中逐渐固定在矿化贝壳中。

令人惊讶的是，尽管腹足纲软体动物是一个多样化的后生动物群体，已知有超过100,000种现生和化石物种，但它们的外壳似乎遵循着基于扩展螺旋的简单位设计，与对数螺旋密切相关。Moseley于1838年首次将对称螺旋识别为贝壳的关键特征。

7.2 几何描述

贝壳的理想形式可以在第五章曲线和杆的理论框架内描述。贝壳的中心轴是一条对数螺旋，半径随弧长增加。表示贝壳曲面的一个方便方法是考虑中心线曲线p = p(t)，然后在曲线上每一点附加三个正交向量d₁(s)、d₂(s)和d₃(s)。

7.2.1 对数螺旋

中心曲线取为三维中的对数螺旋：

p(t) ≡ (x, y, z) = Ae^(αt)(sinβcost, -sinβsint, cosβ)

其中A、α和β是实参数。这条曲线在水平x-y平面中是一条等角螺旋（对数螺旋）。等角螺旋的定义性质是：从中心到曲线上任一点的射线与该点切线之间的夹角是常数。如果选择参数t为x-y平面中的角度θ，则径向线在该点与切向量的夹角为γ = arccot(α)。

对数螺旋是自相似的：给定n个点{pᵢ = p(t₀ + 2iπ), i = 1, ..., n}，三个连续点之间的距离比例：

d(pᵢ₋₁, pᵢ) / d(pᵢ, pᵢ₊₁)

与i的选择无关。这个比例的计算也可以访问α的值。

7.2.2 生成曲线

通过考虑围绕中心线的生成曲线来构建实际的贝壳曲面，并使用相同的指数函数缩放其大小。为简单起见，取生成曲线C为一个圆，其半径随时t指数扩展：

ρ(S, t) = Be^(αt)(cos(S)d₁(t) + sin(S)d₂(t))

这里有两种不同的贝壳参数化方法： 1. 选择平面Π使其包含z轴（Moseley 1838年的原始描述） 2. 定义Π为每一点处曲线的法平面（Okamoto的提议）

第二种选择在数学上更自然和优雅，因为它不依赖于对数螺旋的特殊选择。

---

八、贝壳的增生生长

8.1 生长定律

前面的贝壳表示是纯粹几何的，从外部视角描述形状。现在通过将生长描述为局部增生过程来比较这种方法。

首先施加形状不变性、形状平面性和空间自相似性这三个条件。这三个条件导致线性或平面生长的两种可能生长定律：

q₁ = -c, q₂ = b₁ + b₂cos(S), q₃ = 0         （定律1）
q₁ = -λ(t)c, q₂ = λ(t)(b₁ + b₂cos(S)), q₃ = 0  （定律2）

这两种定律产生相同的结构，只有时间演化可以区分不同的生长过程。

对于过程(7.63)，完成一圈的时间是固定的，与大小或时间无关，因此材料沉积速率必须指数增加。这个生长过程似乎与观察到的贝壳生长率随年龄增加而降低的现象不一致。

8.2 常数生长率假设

假设材料率ṁ（表面积增加率）在时间上恒定，得到自相似形式：

c          b₁ + b₂cos(S)
q₁ = -———, q₂ = ——————————, q₃ = 0
 λ                λ

这个形式具有对数时间的卷曲率，而半径按t^(1/2)增加。因此，假设恒定的生长率，圆形开口的自相似增生结构的存在意味着生长速度必须遵循上述规律。

---

九、其他增生结构

一旦理解了局部增生模型的关键成分，就相当容易理解不同生长向量成分之间的平衡，以重现自然界中观察到的不同形状。例如，几种角和贝壳及其相关的生长向量。

另一个具有非圆形生成曲线的例子是双壳贝壳。这种形状通过使用简单的初始曲线和沿给定方向线性的副法线生长以及随时间变化的卷曲率获得。

---

十、力学在形态模式中的作用

10.1 概述

增生生长的运动学分析揭示了贝壳的对数螺旋结构由简单的局部定律控制。它还揭示了沿开口存在生长梯度。如何建立和维持这个生长梯度，以及贝壳如何识别生长方向，仍然是一个谜。

然而，我们可以将这个基本生长结构作为贝壳形状的底层演化，并将注意力集中在叠加在基本螺旋形状上的形态模式上。有两种基本的形态模式叠加在自相似结构上：共缘装饰（commarginal）和反缘装饰（antimarginal）。

10.2 共缘装饰（Commarginal Ornamentation）

共缘装饰或横向肋骨是许多软体动物贝壳的共同特征，其特征是开口膨胀的振荡分量。生成曲线在整体膨胀因子方面保持相似。

力学模型：将外套膜建模为具有给定生长率的圆形生长弹性环。当外套膜粘附在当前贝壳边缘时，它受到压缩/拉伸力和粘附力。假设作用在外套膜上的压缩力足够小，使其保持圆形，平衡直径是这两种力之间的平衡。

反馈机制：在平衡时，外套膜不是处于无应力状态，假设新层的沉积由应力水平控制，通过调整沉积角度。如果外套膜处于压缩状态，它会增加沉积角度以增加新开口的半径，从而减轻应力；反之，如果外套膜处于拉伸状态，角度减小，开口气孔半径增加率也减小。这个简单的定律自然地产生了由于增生和软组织生长之间的不匹配而导致半径轮廓随时间的振荡。

10.3 反缘装饰（Antimarginal Ornamentation）

反缘装饰是垂直于开口的模式，因此它们是生成曲线形状的修饰，不同于简单的膨胀。这包括反缘肋骨以及更复杂的模式，如刺。

力学模型：通过假设外套膜是一个围绕开口生长的不可伸长弹性环来模拟反缘肋骨的形成。进一步假设下一个增生步骤中的开口取决于外套膜的形状。由于生长外套膜将表现出相对于前一生长增量的长度过剩，它在粘附到开口时将受到屈曲。新沉积的贝壳材料层呈现由下一个基础在矿化时设定的变形形状。因此，外套膜的力学变形决定了开口形状的演化。

分叉分析：遵循生长环理论，可以直接耦合到演化基础。外套膜屈曲时第一分叉发生的临界长度可以计算。结果表明，在特定参数范围内，不同模式（模式2和模式3）并不存在于所有k值，模式2有一个最小值。

10.4 刺的形成（Spine Formation）

刺也是反缘装饰，因为它们改变了生成曲线的形状。然而，与前一个例子的区别在于，它们在外套膜长度大量过剩的阶段出现。刺的形成可以根据外套膜的生长和弹性响应来理解。

基本机制：刺形成的核心思想是预测由具有相对于先前形成和固定配置过剩长度的薄弹性外套膜边缘逐步产生的形式。外套膜边缘被建模为不可伸长的弹性杆，假设以固定速率生长。在每个生长增量中，外套膜通过弹性角质层粘附到当前贝壳边缘。由于外套膜自上次分泌以来已经生长，它比贝壳边缘长，粘附时会发生变形。变形外套膜的形状是作为外套膜弯曲矩和附件力（作为机械刚性基础的服务模板）之间平衡而获得的弹性平衡形状。

形态变化：贝壳的发育过程包括：首先逐步增加外套膜边缘的长度；其次通过力学平衡确定外套膜形状；最后将贝壳边缘更新为变形外套膜的形状。出现形式受外套膜边缘的长度尺寸和弹性特性、粘附强度和边缘生长率控制。

---

十一、本章重要公式汇总

生长向量与曲线演化

• 曲线演化方程：ṙ(S, t) = q(S, t)
• 切向量关系：r' = v = v₃d₃ = λd₃

相容性条件

• q' - v̇ = Wv - Uq
• Ẇ - U̇ = UW - WU

Frenet标架下的Darboux矩阵

U = [0  -τ  κ
     τ   0   0
    -κ   0   0]

形状平面性的生长速度形式

q₂ = b₁(t) + b₂(t)x(S, t) + b₃(t)y(S, t)

形状不变性的生长速度

q = c₁(t)r + c₂(t)r⊥ + a(t)

对数螺旋参数化

p(t) ≡ (x, y, z) = Ae^(αt)(sinβcost, -sinβsint, cosβ)

自相似生长速度（常数材料率假设）

c          b₁ + b₂cos(S)
q₁ = -———, q₂ = ——————————, q₃ = 0
 λ                λ

---

十二、关键概念解释

增生生长（Accretive Growth）

指通过在边界局部添加物质而实现的生长，材料本体无法变形以添加新物质。这是硬组织（如贝壳、角、牙齿）的主要生长方式。

Darboux矩阵

描述曲线标架沿空间方向旋转的反对称矩阵，其元素与曲率和挠率相关。

自相似性（Self-Similarity）

指物体在不同尺度上具有相同形状的特性。对数螺旋贝壳中连续轮之间的距离比例恒定。

细胞轨迹（Cell Track）

在增生生长过程中，细胞沉积新材料的痕迹，可视为细胞在表面上留下的路径。

生长轴（Growth Axis）

控制生成曲线上生长速度线性变化方向的特定轴。

共缘装饰（Commarginal Ornamentation）

沿贝壳生长方向延伸的横向肋骨或装饰，是开口膨胀振荡的结果。

反缘装饰（Antimarginal Ornamentation）

垂直于开口方向的装饰，包括肋骨和刺，是生成曲线形状的修饰。

生长率（Growth Rate）

ṁ = ∫_C |q| ds，表示表面积的增加率。

---

十三、章节总结与思考

本章系统地建立了硬表面增生生长的数学理论框架，并以贝壳为典型范例展示了该理论的应用价值。主要贡献包括：

1. 建立了内禀增生生长的完整数学描述：通过活动标架、Darboux矩阵和自旋矩阵，将生长过程用偏微分方程组描述。

2. 揭示了生长过程的普适特征：形状平面性条件表明存在一个区分性的生长轴，使副法线生长投影到该轴上时呈线性函数。

3. 建立了贝壳几何与生长的联系：将对数螺旋的几何描述与局部增生生长理论相结合，解释了贝壳形态的力学基础。

4. 阐明了形态模式的力学机制：通过外套膜-贝壳相互作用的力学模型，解释了共缘装饰、反缘装饰和刺形成的原理。

值得进一步思考的问题： - 生长方向（即生长轴）是如何被贝壳在生物学上识别和维持的？ - 不同物种的贝壳形态多样性背后的遗传控制机制是什么？ - 增生生长理论能否推广到其他具有类似生长模式的生物结构？

---

注：本笔记基于《The Mathematics and Mechanics of Biological Growth》第7章内容整理，章节原名"Accretive Growth"，作者为Alain Goriely。本笔记仅供学习研究使用。