第四章阅读笔记：线上的生长（Gowing on a Line）
概述
本章探讨了一维丝状结构（filamentary structures）的生长问题，从微观分子链到宏观植物茎干、神经轴突等都属于丝状结构的范畴。这些结构虽然在尺寸和微观组成上差异巨大，但其生长、运动和形状变化遵循着普遍的物理规律。本章建立了形态弹性杆理论（morphoelastic rod theory）的基本框架，系统分析了生长与弹性变形之间的耦合关系。
1. 丝状结构的三种生长方式
丝状结构的生长可以归纳为三种基本模式：
1.1 轴向生长（Axial Growth）
最常见的生长方式，通过增加丝状结构中心线的长度来实现。这种生长可以发生在顶端（如根尖、芽尖），也可以沿整个长度分布。轴向生长是植物根系、茎干伸长的主要机制。
1.2 径向生长（Thickness/Density Growth）
通过增加丝状结构的厚度或密度实现，此时长度保持不变。在树木中形成年轮的过程称为次生生长（secondary growth），通常发生在初生伸长之后。
1.3 差异生长（Differential Growth）
当丝状结构横截面上不同位置的轴向生长速率不同时，会产生曲率和扭转。植物的向重性（gravitropism）反应就是差异生长的典型例子——将植物水平放置时，下侧生长快于上侧，直至恢复垂直方向。
差异生长在三维情况下还会引发扭转和螺旋生长，常见于缠绕藤本植物、某些细菌和真菌中。
2. 纯弹性变形（Purely Elastic Deformations）
在一维约束下，杆件只能发生拉伸或压缩变形。
2.1 拉伸 stretch 的定义
对于初始弧长为 \(S_0\) 的杆件，当前构型中的位置为 \(s = s(S_0, t)\)，局部拉伸（stretch）定义为： $\(\lambda(S_0, t) = \frac{\partial s}{\partial S_0}\)$
当 \(\lambda > 1\) 时表示伸长，\(\lambda < 1\) 表示压缩。
2.2 本构关系
胡克定律（Hookean law）： $\(\sigma = E(\alpha - 1) = E\epsilon\)$
其中 \(E\) 为杨氏模量，\(\epsilon = \alpha - 1\) 为应变。
胡克定律仅在小变形下有效。当 \(\sigma = -E\) 时，预测任何线元会被压缩至一点，这显然是非物理的。
neo-Hookean 非线性定律（适用于大变形）： $\(\sigma = \mu(\alpha^2 - \alpha^{-1})\)$
其中 \(\mu\) 为剪切模量。在 \(\alpha \approx 1\) 附近，可恢复胡克行为，且 \(E = 3\mu\)。
3. 无弹性变形的生长（Growth Without Elastic Deformations）
3.1 生长拉伸（Gouth Stretch）
生长可以用生长拉伸 \(\gamma\) 来描述，类似于弹性拉伸。若生长过程在时间上演化，满足演化方程： $\(\frac{\partial \gamma}{\partial t} = G(\gamma, s, S_0)\)$
均匀线性生长：\(G = 1\)，则 \(s(S_0, t) = t S_0\)
指数生长（细胞以恒定速率分裂）： $\(\frac{\partial \gamma}{\partial t} = k\gamma \implies \gamma(t) = e^{kt}\)$ $\(s(S_0, t) = S_0 e^{kt}, \quad l(t) = L_0 e^{kt}\)$
其中 \(k\) 为相对生长速率。
3.2 顶端生长（Tip Growth）
顶端生长是许多生物丝状结构的通用机制，物质仅在末端添加。考虑一个长度为 \(a\) 的生长区：
通过数值积分（前进欧拉法）可以模拟生长过程。重要结论： - 尽管顶端指数增长，整体长度呈线性增长：\(l(t) = L_0 + akt\) - 这是因为扩展区长度保持恒定
3.3 应用：球形肿瘤生长（Spheroid Tumor Growth）
作为局部化生长的应用示例，考虑球形肿瘤的生长。假设生长是各向同性的，体积元满足： $\(dv = \eta(r, R_0) dV_0\)$
营养物质浓度 \(u(r, t)\) 满足扩散方程。当肿瘤半径超过临界值 \(b_{cr}\) 时，中心区域营养耗尽形成坏死核（necrotic core）。
临界半径： $\(b_{cr} = \frac{6D u_b}{Q}\)$
关键发现： - \(t < t_{cr}\)：肿瘤呈指数生长 - \(t > t_{cr}\)：坏死核形成，仅外围壳层增殖，生长转变为线性
这与顶端生长系统的行为相似。
4. 带弹性变形的生长（Growth With Elastic Deformation）
4.1 形态弹性理论的基本假设
核心假设：全轴向拉伸 \(\lambda\) 是生长拉伸 \(\gamma\) 和弹性拉伸 \(\alpha\) 的乘积： $\(\lambda = \alpha \gamma\)$
这一看似简单的乘法分解实际上非常精妙，在非线性弹性理论的框架下才能完全理解和证明。
4.2 两板间生长的杆件
考虑一个在两块不可移动板之间生长的杆件。初始长度 \(L_0\) 等于板间距，因此总长度保持不变 \(\lambda = 1\)。
若杆件均匀线性生长 \(\gamma = 1 + t\)，则： $\(\alpha = \frac{1}{\gamma} = \frac{1}{t+1}\)$
胡克定律下的应力： $\(\sigma = -\frac{Et}{t+1} \xrightarrow{t \to \infty} \sigma_\infty = -E\)$
应力最终饱和至 \(-E\)。
neo-Hookean 定律下： $\(\sigma = \mu\left(\frac{1}{(1+t)^2} - t - 1\right)\)$
应力发散至 \(-\infty\)。
4.3 虚拟配置（Virtual Configuration）
生长的杆件在两块板之间的约束可以分解为两个过程： 1. 生长变形：杆件在无约束状态下生长至长度 \(L(t) = L_0(1+t)\)，此为虚拟配置（应力自由） 2. 弹性变形：将生长后的杆件压缩回长度 \(l = L_0\) 以满足边界条件
关键概念：应变必须相对于参考配置来定义。对于生长材料，参考配置是不断演化的。虚拟配置是一个数学上的有用构造，通常无法物理实现。
4.4 动态平衡生长（Homeostatic Growth）
假设丝状结构有一个动态平衡应力 \(\sigma^*\)——即通过主动生长过程趋向恢复的应力。
演化定律： $\(\frac{\partial \gamma}{\partial t} = \gamma(\sigma - \sigma^*)\)$
当 \(\sigma = \sigma^*\) 时，生长停止。
考虑一个初始处于平衡状态的丝状物，突然被拉伸或压缩至长度 \(l\)，其长度将指数趋近于新的平衡长度。
4.5 应用：神经元的生长（Growth of Neurons）
神经元是应力介导生长的典型例子。轴突的延长依赖于机械力： - 生长锥（growth cone）通过肌动蛋白-马达-离合系统产生张力 - 外部施加的张力可以控制轴突的伸长速率
实验表明，在临界张力以上，生长速率与施加的张力呈线性关系。
简化模型： - 轴突被视为一维弹性结构 - 存在粘附力 \(f = K(s(S) - S)\) - 生长速率取决于张力是否超过临界值
5. 生长与塑性（Is This Just Plasticity?）
5.1 相似性与差异
对塑性理论熟悉的读者可能会注意到生长变形分解与塑性变形的相似性。然而，生物生长与金属塑性有本质区别：
| 方面 | 塑性理论 | 生长理论 |
|---|---|---|
| 微观机制 | 滑移线、剪切应力 | 细胞分裂、壁松弛 |
| 进化规律 | 由微观结构决定 | 由生化调控 |
| 应力依赖 | 临界屈服应力 | 生理条件 |
5.2 不同的材料模型
理想流体：\(\dot{\epsilon}_{an} = \sigma / \eta\)，仅有不可逆流动
Maxwell 流体：弹性 + 粘性响应 $\(\dot{\epsilon} = \frac{\sigma}{\eta} + \frac{\dot{\sigma}}{E}\)$
Kelvin 固体：粘弹性，响应发生在有限时间内 $\(\sigma = E\epsilon + \eta\dot{\epsilon}\)$
Bingham 材料：存在临界屈服应力 $\(\dot{\epsilon}_{an} = \frac{1}{\nu}(\sigma - \sigma^*)H(\sigma - \sigma^*)\)$
Maxwell-Bingham 材料：弹性 + 塑性屈服
Lockhart 植物生长模型本质上等同于 Maxwell-Bingham 系统。
6. 植物细胞生长的应用（The Growth of Plant Cells）
6.1 Lockhart 模型
植物器官（如茎、胚芽鞘）的伸长主要由细胞扩张驱动。植物细胞被细胞壁包围，细胞壁是含有纤维素微纤维的凝胶复合结构。
渗透压与膨压： $\(\Delta P = P_{osm} - P_{tur}\)$
当细胞充分膨胀时，\(\Delta P = 0\)。
Lockhart 方程（可逆水通透）： $\(\frac{1}{\lambda}\frac{d\lambda}{dt} = K \Delta P\)$
6.2 不可逆细胞生长
细胞体积的不可逆增大（生长）涉及细胞壁松弛。当膨压超过临界值 \(P^*\) 时，发生不可逆延伸：
其中 \(\phi\) 为延伸性系数。
完整的 Lockhart 方程： $\(\frac{1}{\lambda}\frac{d\lambda}{dt} = \frac{\phi K}{\phi + K}(P_{osm} - P^*)H(P_{osm} - P^*)\)$
6.3 模型应用
通过控制膨压的不同技术（渗透溶液、压力探针等），可以测量生长速率与膨压的关系，从斜率获得 \(\phi\)，从截距获得 \(P^*\)。
6.4 扩展 Lockhart 模型
考虑弹性和生长的耦合： $\(\frac{1}{\lambda}\frac{d\lambda}{dt} = \phi(P_{tur} - P^*)H(P_{tur} - P^*) + \frac{1}{E}\frac{dP}{dt}\)$
该方程在形式上等价于 Maxwell-Bingham 材料。
7. 关键概念总结
- 丝状结构：一个维度远大于其他两个维度的三维物体，普遍存在于自然界
- 乘积分解：\(\lambda = \alpha \gamma\) 是形态弹性理论的核心
- 三种生长模式：轴向生长、径向生长、差异生长
- 虚拟配置：生长产生的应力自由构型，是理解应变的关键
- 动态平衡：生物系统趋向于维持特定应力状态
- Lockhart 模型：植物细胞生长的经典定量模型
- 生长与塑性的区别：生长是由生化调控的主动过程，而非简单的物理塑性流动
8. 重要公式一览
| 公式 | 描述 |
|---|---|
| \(\lambda = \partial s / \partial S_0\) | 拉伸定义 |
| \(\sigma = E(\alpha - 1)\) | 胡克定律 |
| \(\sigma = \mu(\alpha^2 - \alpha^{-1})\) | neo-Hookean 定律 |
| \(\lambda = \alpha \gamma\) | 形态弹性分解 |
| \(\partial\gamma/\partial t = \gamma(\sigma - \sigma^*)\) | 动态平衡生长 |
| \(\frac{1}{\lambda}\frac{d\lambda}{dt} = \frac{\phi K}{\phi + K}(P_{osm} - P^*)\) | Lockhart 方程 |
9. 读后感
本章建立了研究一维生物丝状结构生长的基本数学框架。通过将总变形分解为弹性分量和生长分量，我们能够分别处理物理变形和生物过程。这一分解看似简单，但蕴含着深刻的物理洞察——参考构型随着生长不断演化，虚拟构型的概念对于理解应力状态至关重要。
植物细胞生长的 Lockhart 模型是一个特别成功的应用范例，它将力学与生理学完美结合，揭示了膨压、渗透压和细胞壁性质在细胞伸长中的相互作用。神经元生长的例子则展示了力学信号如何调控细胞发育。
与塑性理论的对比提醒我们，生物生长是一个高度调控的主动过程，不能简单套用金属塑性的唯象理论。