第16章 生长中的圆柱体：读书笔记
概述
本章探讨了生物体中圆柱状结构的生长力学问题，涵盖多种典型的管状结构如血管、植物茎干、气道、食道等。残余应力在管状结构中起着至关重要的作用：它调节血管中的周向应力梯度、提高生长中植物茎干的刚性、机械稳定动脉、并协助气道和食道的正常功能。本章从空腔形成、生长结构的失稳现象、以及某些管状真菌的旋转等几个方面进行了系统分析。
16.1 生长圆柱体的运动学
基本模型
本章考虑一个由N个不同弹性、各向同性、不可压缩圆柱壳层组成的生长圆柱结构，每个壳层由对角生长张量定义，具有不同但均匀的材料特性。初始构型B0中，结构具有内半径A0和外半径An≡B0，高度H0。该圆柱体进一步细分为n个圆柱壳层，第i个壳层的内外半径分别为Ai-1和Ai。
变形几何
假设圆柱体在所有变形过程中保持其完整性和圆柱几何形状。从初始构型到当前构型的总变形在圆柱坐标{R0,Θ,Z}到{r,θ,z}中给出：
r = r(R0), θ = Θ, z = ζZ
其中ζ是圆柱体的常数轴向拉伸。引入乘法分解F = AG，其中弹性变形张量A由弹性拉伸αi(R0) = (αr, αθ, αz)定义，生长张量由生长向量γ^i = (γi,r, γi,θ, γi,z)定义。
不可压缩条件
由于不可压缩约束，det F = det(AG) = det G = g，即：
dr/dR0 · r/R0 = g(R0)/ζ
应力分析
由于变形在圆柱坐标中是对角的，凯莱应力张量也是对角的：[T] = diag(tr, tθ, tz)。径向应力通过积分给出，在各层界面处连续，但周向和轴向应力分量通常不连续。边界条件包括：内压力P和轴向力N，它们共同决定当前内半径a和轴向拉伸ζ。
16.2 植物中的空化作用
背景
许多植物如蒲公英、细香葱和竹子具有空心茎。许多根部和水生植物的组织中含有大型通气组织（aerenchyma）。这些组织中的气体囊具有多种功能：机械经济性、空心茎的刚性、水生植物的浮力、气体运输和液流增加、叶绿素生产中光散射的变化以及种子传播。
空化的两种机制
十九世纪，Sachs将通气组织的形成分为两类： 1. 裂解性（lysigenous）：通过均匀生长和细胞死亡形成 2. 裂生性（schizogenous）：通过差分生长在组织中产生张力，导致相邻细胞的机械撕裂
Takano的实验
Takano等人分析了赤霉素和机械应力对豆类植物茎空化的影响： - 在天然实心茎的矮生豆类植物中施加赤霉素会增加茎的长度、减小直径并诱导茎空化 - 在天然空心茎的pole豆类植物中，抑制赤霉素或施加机械应力会增加茎的粗度、减少轴向长度并阻止茎空化
模型分析
假设茎是一个不可压缩的neo-Hookean圆柱体，初始外半径为B0，生长遵循各向同性但非均匀的线性规律：
γr = γθ = ν1 + ν2R0
γz = μ1 + μ2(B0 - R0)
计算表明，在柱芯处产生的拉伸径向和周向应力可能足以导致材料破裂。通过与已知植物髓组织的断裂应力比较，证实了机械撕裂导致通气组织形成的假设。
16.3 生长圆柱体的分岔
欧拉屈曲与桶形失稳
在轴向压缩载荷下，弹性杆会发生屈曲。欧拉在1744年首次推导了屈曲临界载荷，拉格朗日在1770年进一步完善。经典欧拉屈曲公式：
NEuler = π²EI/L²
对于短粗结构，桶形失稳（barreling）成为主导模式。周向数n=0对应轴对称桶形模式，n=1对应欧拉屈曲模式。
分岔分析
对生长圆柱体应用分岔方案，包括以下步骤： 1. 运动学：展开变形χ = χ^(0) + χ^(1) 2. 线性化方程：将凯莱应力张量展开 3. 边界值问题：通过傅里叶展开简化依赖关系
有效刚度
生长过程中，由于残余应力场的存在，结构的有效属性发生变化。通过"思想实验"测量生长圆柱体在零载荷下的几何参数，然后施加法向载荷直至屈曲，得到有效杨氏模量：
Eeff = 4Pcrit/π²σ²
其中σ = (b0² + a0²)^(1/2)/l0是短粗比。
轴向生长的影响
对于径向依赖生长的各向同性neo-Hookean圆柱体，轴向生长会降低结构的稳定性。差分生长在圆柱体内部产生大的压缩区域，使圆柱体因生长而预压缩，因此在更小的载荷下就会屈曲。
16.4 植物中的组织张力
背景现象
当切开花茎（如蒲公英）时，它会向后弯曲。这种形状变化是由于原始圆柱构型中存在残余应力。植物生理学家称之为"组织张力"。在生长过程中，外层组织（表皮）的静止长度与内层较软组织（髓）的静止长度不同。
机械意义
分离的单个大黄组织柔软且无法在自身重量下保持稳定。然而，组合构型是刚性的，允许植物向上生长。实验表明，外层组织在整体刚性中起着重要作用，其横截面积不到10%却贡献了高达70%的刚性。
双层模型
将茎建模为具有不同生长和弹性特性的双层圆柱壳。内部组织（髓）建模为neo-Hookean材料，外部组织（表皮加皮质薄壁组织）建模为Fung材料（具有强应变刚化效应）。
关键发现：差分轴向生长结合非线性应变刚化是生物系统中控制结构特性的简单通用机制。外部层处于张应力状态，通常非常薄，可以近似为表面条件。
16.5 动脉的屈曲
动脉的特殊力学性质
动脉系统是一个复杂的机械系统，具有以下独特性质： - 弹性响应是非线性的（Roy, 1881） - 在生理条件下，动脉被束缚在周围组织上，工作在约1.3-1.6的大拉伸状态 - 尽管压力有很大变化，主要动脉的端到端距离基本保持不变（Patel和Fry） - 即使在切除后、无应力和无压力的情况下，动脉仍然有轴向应力 - 动脉也存在径向残余应力，通过"开角"测量
动脉的双层结构
动脉壁由三层组成：内膜（中空）、中膜和外膜。由于内膜的机械贡献可以忽略不计，动脉被简化为双层结构（media和adventitia）。采用Holzapfel-Gasser-Ogden模型，每层由胶原纤维增强，各层具有不同的胶原纤维方向。
生长张量
假设生长张量在每层中是常数和对角的：
[G] = diag(γr, 1, γz) (内层)
[G] = diag(γr, 1, 1) (外层)
利用开角数据确定生长张量的参数值。
分岔分析的主要结论
-
材料非线性的作用：应变刚化行为对于捕捉心动周期中动脉的正确行为非常重要。这些非线性性对于解释动脉在轴向张力下的"反转线"至关重要。
-
径向生长的作用：周向残余应力大大减少了跨壁应力梯度，这是差分生长对动脉影响的著名特征。
-
轴向生长的作用：轴向残余应力在防止屈曲中起重要作用。有残余应力的动脉在轴向拉伸1.5的情况下，能够承受更大的内压才发生屈曲。
-
临床意义：动脉屈曲是动脉扭曲的标志，可能发生在异常生长刺激或遗传缺陷后，可能导致严重的并发症包括中风、眩晕、晕厥和持续性耳鸣。
16.6 周向屈曲与粘膜折叠
环向失稳机制
考虑一个双层圆柱壳，每层具有恒定的径向和周向生长但无轴向生长。取决于每层的几何形状和刚度，会产生各种有趣的失稳模式： - 如果外层是刚性的，内层各向同性生长，内层将发生周向屈曲 - 如果内层是刚性的，外层各向同性生长，将产生花状图案
这些模式被称为"周向屈曲"或"粘膜折叠"，在食道、血管和胃肠道的粘膜中观察到。
单层圆柱体的分析
考虑在生长和外压下、无端位移的单层圆柱壳。通过稳定性分析，得到临界外压作为径向生长γ的函数。径向吸收（γ < 1）是失稳的，而初始的径向生长（γ > 1）有稳定作用，但在γ较大时变为失稳。
双层圆柱体
对于双层圆柱体，临界压力取决于内层与外层生长率的比值γ1/γ2： - 内层生长较快 → 失稳（产生压缩应力） - 外层生长较快 → 稳定（产生张应力）
16.7 哮喘与气道重塑
哮喘的力学特征
哮喘是一种以气道狭窄和肺功能下降为特征的复杂疾病。气道重塑涉及多个时间尺度和空间尺度，涉及许多不同的化学、生化和物理刺激。气道壁厚度增加是哮喘重塑的关键特征。
气道折叠的力学
在某些刺激下，围绕气道壁的平滑肌收缩，腔壁折叠或屈曲。在正常生理条件下，这只会导致适度的气道狭窄。然而，在哮喘患者中，这种狭窄往往发生在较低的刺激下，并导致过度狭窄。
双层模型
将气道建模为双层圆柱形态结构： - 内层（粘膜）：薄而硬，包括基底膜、固有层和上皮 - 外层（粘膜下层）：软而厚，包括疏松结缔组织
平滑肌围绕外层，收缩时提供作用在气道壁-平滑肌界面的法向力。
关键发现
-
各向同性差分生长的作用：如果硬粘膜比软粘膜下层生长更快，气道实际上变得不太稳定。这似乎违反直觉，但差分生长创造了残余应力，导致内层处于压缩状态。
-
各向异性生长的作用：内层的各向异性生长（仅径向生长）可增强气道，但显著降低屈曲模式数。
-
临床意义：这一分析表明，差分生长可能是理解经历重塑后气道力学行为的重要因素。
16.8 纤维收缩产生的残余应力
纤维增强材料
除了通过不均匀或各向异性生长产生残余应力外，另一种可能的机制是各向异性材料中纤维的收缩或伸长。考虑具有单族预压缩纤维的薄壳圆柱体，螺旋纤维缠绕在圆柱体周围。
旋转方向的反转
一个显著发现：对于适当选择的参数，旋转方向可以完全反转。具体来说： - 对于小预压缩，螺旋纤维产生的旋转方向与直觉预期相反 - 随着预压缩增加到临界值，旋转方向发生反转 - 这种反转取决于纤维角度和预压缩程度
这表明：即使螺旋纤维是右旋的，也可能产生左旋或右旋旋转，取决于角度和预压缩水平。
16.9 Phycomyces的螺旋生长
Phycomyces blakesleeanus
这种真菌产生一个显著的单个细胞的气生菌丝（孢子囊柄），长度可达10厘米。其细胞壁是各向异性的，主要由嵌入弹性基质中的几丁质微纤维组成。
生长阶段
该真菌经历四个生长阶段： - Stage I：初始气生生长阶段 - Stages II和III：孢子囊发育 - Stage IV：轴向生长恢复，具有惊人的螺旋生长反转现象
在Stage IVa（约1小时）观察到右旋生长（RH），随后是Stage IVb的左旋生长（LH），持续长达48小时。
纤维沉积模型
基本假设： 1. 新纤维沿现有纤维方向以无应力状态沉积 2. 生长区线性扩展到恒定大小
关键机制：当前构型在给定时间步成为下一时间步的参考构型。因此，材料点中的有效应力是具有不同预压缩的纤维的组合。
旋转反转的物理解释
- Stage IVa：只有生长区顶部附近的纤维起作用，净旋转是逆时针的
- Stage IVb：观察到的旋转由整个区的积分决定，导致左旋螺旋
- 在任何给定时间，生长区中同时存在逆时针旋转（顶部）和顺时针旋转（底部）的单个材料点
本章总结
本章系统分析了生长圆柱体的力学问题，主要结论：
-
残余应力的普遍性：在各种生物圆柱结构中，残余应力起着关键的机械作用，从植物茎干的气道组织张力到血管的应力调节。
-
生长与失稳的关系：差分生长可以在没有外部载荷的情况下引起失稳（自发屈曲），也可以增加或减少结构对外部载荷的稳定性。
-
材料非线性的重要性：应变刚化行为（如Fung模型）对于捕捉生物材料的真实响应至关重要，特别是在动脉和气道中观察到"反转线"现象。
-
理论与临床的关联：动脉屈曲导致扭曲可能是buckling instability的结果，气道重塑中的矛盾现象（壁变厚同时稳定性降低）可以用差分生长解释。
-
螺旋生长的纤维机制：Phycomyces的旋转反转是由纤维沉积和预压缩的演化导致的，无需假设纤维手性本身决定宏观旋转方向。
这些分析为理解生物体中管状结构的生长、形态发生和稳定性提供了坚实的理论基础。