第八章：薄膜与壳体 阅读笔记
概述
本章从硬材料构成的二维结构转向软生物结构的力学建模，将生物结构视为Morphoelastic薄膜、薄板或壳体。主要动机是描述丝状细胞（如真菌、放线菌、根毛、花粉管）的顶端生长模式。本章假设表面始终保持轴对称，从而将可能的变形限制在满足特定对称性的范围内，并基于简单基本原理推导完整的力学理论。
8.1 弹性薄膜（Elastic Membranes）
8.1.1 运动学（Kinematics）
几何参数化： 轴对称薄膜表面S由平面曲线C绕z轴旋转而成。曲线C由材料坐标S参数化，从表面与z轴的交点O开始计数。在任意时刻t，当前构型由曲线C描述：
r(S, t) = r(S, t)er + z(S, t)ez , S ∈ [0, L]
几何关系： 引入角θ(S, t)——曲线C在S处法向与z轴的夹角。弧长s = s(S, t)从O点测量。几何变量满足：
∂r/∂s = cos θ
∂z/∂s = − sin θ
单位切向量和法向量：
τ = cos θer − sin θez
ν = sin θer + cos θez
主曲率： 两个重要几何量是主曲率
κs = ∂θ/∂s
κϕ = sin θ/r
伸长比（Stretches）： - 径向伸长比：αϕ(S, t) = r(S, t)/R(S)，新半径与原始半径之比 - 切向伸长比：αs(S, t) = ∂s(S, t)/∂S，相对于参考弧长的伸长量 - 法向伸长比：αn = h(S, t)/H，薄膜厚度的变化
不可压缩条件： 对于不可压缩材料，αs·αϕ·αn = 1。对于脂质双层，αn = 1且αs·αϕ = 1（面积不可压缩）。
8.1.2 力学（Mechanics）
表面张力定义： - ts：沿切向τ方向的线张力（单位：力/长度） - tϕ：沿方位角ϕ方向的线张力（环向）
力平衡： 考虑面积rdsdϕ的表面微元，受外部体力f = τcτ + Pν，其中τc为周围介质对薄膜的曳引力，P为跨膜压差。
平衡方程推导： 通过力平衡得到三个标量方程：
∂tϕ/∂ϕ = 0
θ′ rts + tϕ sin θ − rP = 0
∂(rts)/∂s − tϕ cos θ + rτc = 0
广义Laplace方程（Young-Laplace方程）：
P = κs ts + κϕ tϕ
第一积分（First Integral）： 在恒压条件下，存在沿曲线C保持不变的积分：
C = r²(2tsκϕ − P)
8.1.3 本构定律（Constitutive Laws）
线性定律： 对于小变形，各向同性不可压缩薄膜：
ts = A(ess + 2eϕϕ)
tϕ = A(eϕϕ + 2ess)
用伸长比表示：
ts = A(2αs² + αϕ² − 3)
tϕ = A(αs² + 2αϕ² − 3)
8.1.4 完整方程组
机械平衡方程（无剪切应力）简化为关于薄膜变量的方程组：
∂s/∂S = αs
∂z/∂S = −αs sin θ
∂r/∂S = αs cos θ
∂θ/∂S = [αs sin θ·r² + R²(S)(2αs² − 3) − pR²(S)] / [r·(2r² − 3R(S)² + αs²R²(S))]
第一积分闭合：
C = 2Ar sin θ·[αs² + 2(r²/R²(S)) − 3] − r²P
αs² = 3 − 2(r²/R²(S)) + p·(r/(2 sin θ))
8.1.4.1 简单例子
球体膨胀： 半径R0的球在压力p下膨胀到半径R1：
R1 = R0√(pR0/12 + √(p²R0²/144 + 1))
无限圆柱： θ = π/2, αs = 1, R = R0, r = R1：
R1 = (1/2)√(p²R0⁴ + 4R0² + pR0²)
顶端生长的启示： 顶端形状可近似为连接在长圆柱支撑上的球冠。球和圆柱的两个简单解表明，为使薄膜连接圆柱和球，参数p需要急剧变化——从球到圆柱需减小约六倍。因此需要弹性参数A相应增大。这说明为形成弹性薄膜支撑的顶端形状，顶端附近必须存在显著的刚度变化。
8.1.4.2 边界条件
- 压力恒定且表面穿过z轴时，使用上述方程组
- S = 0在轴交点处：s(0) = 0
- z(L) = z0设定轴向平移
- 边界处设定r(L) = RL（固定半径）或θ(L) = θL（固定角度）
- r(0) = 0的条件，但由于R(0) = 0，方程(8.37)存在可去奇点，需用小参数ε展开处理
- 若固定体积而非压力，则压力成为体积约束的拉格朗日乘子
8.1.5 应用：海胆的形状
D'Arcy Thompson的洞见： Thompson最早意识到海胆形状与桌上液滴的相似性。"海胆外壳的一系列形态不过是水滴滴状形态的体现。"
力学类比： 海胆形状是体腔压力（coelomic pressure）、管足（podia）和自身重力的力学平衡结果。
海胆骨骼（Test）： 由大量方解石板组成，通过胶原纤维松散连接。生长过程中，新板在顶端插入，各板通过边缘的增生、收缩和 fitting过程扩展。
模型假设： 假设线张力ts和tϕ在所有点相等以最大化强度。此时ts在表面所有点为常数。假设海胆平放，压力随重力线性增加P = ρgz。
毛细长度： 参数a = ρg/ts是毛细长度（表面张力力与重力相同量级的特征长度）的倒数。
数值求解： 从原点附近的一阶展开开始：
θ(ε) = aε, r(ε) = ε, z(ε) = 2, ε ≪ 1
形态空间（Morphospace）： 大a值时形状趋于球形，a减小则越来越扁平。这产生一维形态空间，可与实际物种形状惊人地匹配。
张力相等的困难： 若张力相等则常数为常数，αϕ也常数，意味着r(S) = αs·R(S)。所得形状必须从具有相同形状的无应力薄膜实现。若滴状解确实是海胆骨骼形状的有效表示，则必须同时存在：1）参考形状向当前形状的松弛过程；2）平衡线张力以沉积新材料的生长过程。
8.2 非线性弹性壳体（Nonlinearly Elastic Shells）
现在考虑轴对称壳体——可由表面表示、能抵抗弯曲、具有非平面无应力形状的弹性物体。运动学描述与薄膜相同，直接进入力学部分。
8.2.1 力学
平衡方程： 通过在表面微元上平衡线性和角动量，得到三个方程：
d(rqs)/ds = rqn − rκsts + κϕtϕ (法向应力平衡)
d(rts)/ds = tϕ cos θ + rκsqs − rτs (切向应力平衡)
d(rms)/ds = mϕ cos θ − rqs (角动量平衡)
符号说明： - ts, tϕ：分别为经向和方位角向应力 - ms, mϕ：弯曲力矩 - qs：表面法向剪切应力 - qn：壳体所受总法向应力（压力驱动时qn = P；压力+主动力时qn = P + τn） - τs：外部切向剪切应力（壳体变形时与环境的摩擦）
8.2.1.1 本构关系
应变不变量： 假设材料各向同性，应变能密度函数W（单位体积能量）只依赖于伸长比(αs, αϕ, αn)的三个应变不变量：
I1 = αs² + αϕ² + αn²
I2 = αs²αϕ² + αs²αn² + αϕ²αn²
I3 = αs²αϕ²αn²
不可压缩条件： αn = 1/(αs·αϕ)，I3 = 1
线张力公式：
ts = 2Hαn(αs² − αn²)[∂W/∂I1 + αϕ²·∂W/∂I2]
tϕ = 2Hαn(αϕ² − αn²)[∂W/∂I1 + αs²·∂W/∂I2]
Mooney-Rivlin模型：
Wmr = C1(I1 − 3) + C2(I2 − 3)
Fung模型（软组织）：
Wfu = (C1/γ)[e^(γ(I1−3)) − 1]
弯曲力矩： 各向同性并与表面平均曲率变化成比例：
mϕ = ms = B(κs + κϕ − 2K0)
8.2.1.2 完整壳体方程与薄膜极限
六个常微分方程描述完整系统：
dz/dS = −αs sin(θ)
dr/dS = αs cos(θ)
dθ/dS = αsκs
dκs/dS = αs[(cos θ sin θ)/r − κs] + αs·qs/B
dts/dS = αs·(A/r)cos θ·(fϕ − fs) + αs·κs − qs/A − τs/A
dqs/dS = αs·(qn/A − κs·fs − (sin θ/r)·fϕ − (qs cos θ)/(A·r))
薄膜极限： 忽略弯曲力矩时，qs = 0，此时：
qn = A(κs·fs + κϕ·fϕ) = κs·ts + κϕ·tϕ
8.2.2 标度分析
无量纲数：
ξ = Peff·w / 3A
ξ = Peff·w / HE
物理意义： ξ衡量变形的大小程度： - 高ξ：软材料或高压 → 大弹性变形 - 低ξ：刚性材料或低压 → 小弹性变形
8.2.3 应用：稻瘟病菌（Rice Blast Fungus）
8.2.3.1 背景
** Magnaporthe grisea（稻瘟病菌）：** 攻击水稻叶片的病原菌。形成圆顶状结构——附着胞（appressorium），产生巨大膨压穿透水稻细胞壁感染植物。
侵染过程： 1. 小分生孢子落在叶片表面 2. 萌发管从分生孢子长出，钩入叶片表面 3. 开始形成附着胞 4. 萌发管最终与附着胞分离，附着胞成为独立单元
附着胞结构： 通过基部一圈强力粘合剂附着于植物表面。壁厚约0.1 μm。黑色素层约100 nm厚，允许水分子渗透，与内部约3.2 M甘油浓度结合，可产生约0.3–0.5 MPa（3-5个大气压）的内部压力。黑色素化后，膨压增加至少三倍，可达8 MPa（80个大气压）。
穿透机制： 膨压积累后，穿透钉（penetration peg）穿透叶片表面。力之大甚至可穿透薄聚酯薄膜。
8.2.3.2 模型
几何建模： 附着胞简化为轴对称弹性壳体，特征长度w ≈ 3 μm（附着胞半径）。壳体通过线性基础在靠近表面的小区域粘合到叶片表面：
g = −µz·ez
边界条件：
(τs, τn) = (gs, gn) 当 0 ≤ σ ≤ ηL
= 0 其他
参数估计： - 典型半径：R ≈ 3 μm - 压力：P ≈ 1-8 MPa - 壁厚：h ≈ 0.1 μm - 杨氏模量：E ≈ 10-100 MPa - 弯曲模量（薄板理论）：B = EH³/[12(1−ν²)] ≈ 10-100 × 10⁻¹³ N·m（假设ν = 1/2不可压缩材料）
非线性应力刚化： 若压力增加五倍而半径基本不变，则表观弹性模量也增加五倍。这种刚化行为与Fung模型（γ参数较大）兼容。
8.2.3.3 分析结果
Neo-Hookean模型： 对于黑色素化后的正常压力水平，穿透前，简单的neo-Hookean模型能正确捕捉附着胞形状。
形状稳定性问题： Neo-Hookean模型预测大压力增加会导致大变形，这与实验观察不符（压力大幅增加但尺寸无明显变化）。
Fung模型的优势： 使用Fung应力刚化模型能很好地捕捉这种"高压下形状稳定"的行为。这清楚表明了非线性和大变形理论的重要性。
核心概念总结
薄膜 vs 壳体
- 薄膜：仅能抵抗张力，不能抵抗压缩或弯曲
- 薄板：初始平面结构，能抵抗弯曲、张力和压缩
- 壳体：初始曲面结构，能抵抗弯曲、张力和压缩
Kirchhoff-Love壳体理论
基于Cosserat曲面或带约束的Kirchhoff-Love壳理论，可从标准理性力学推导流体膜（脂质双层）模型。
顶端生长（Tip Growth）
丝状细胞（真菌、放线菌、根毛、花粉管）的共同生长模式： - 生长和新细胞壁材料的加入集中在细胞顶端 - 生长模式反映了细胞壁各向异性、机械应力和生长过程之间的相互作用 - 轴对称假设下可建立完整的力学理论
生长引起的形状变化
- 海胆案例：体腔压力、管足和重力竞争下，柔软板片的松弛和生长使形状呈滴状
- 稻瘟病菌案例：附着胞在巨大内部压力下保持形状稳定，需要应力刚化材料模型（如Fung模型）而非简单neo-Hookean模型
关键方程
- Laplace/Young方程：P = κs·ts + κϕ·tϕ
- 薄膜平衡：六个ODE方程描述几何和力学耦合
- 壳体完整方程：十个方程（六个ODE + 四个代数关系）
生物学启示
本章建立的弹性薄膜和壳体理论框架对理解生物形态发生具有重要意义：
- 刚度梯度：顶端生长需要局部刚度变化（参数p减小/A增大）
- 应力刚化：高压下形状稳定表明细胞壁具有非线性本构响应
- 多尺度结构：从分子（黑色素层）到细胞（附着胞）到组织（叶片穿透）的力学耦合
- 生长与力平衡：生物形状是生长过程与力学平衡共同作用的结果