跳转至

第五章：弹性细杆（Elastic Rods）读书笔记

概述

本章系统介绍了弹性细杆的理论框架，这是研究生物生长形态发生的核心数学工具。弹性细杆理论的核心思想是将三维管状结构内部的应力替换为作用于曲线上的合力和合力矩向量，这一过程导致基尔霍夫方程（Kirchhoff equations）的建立。弹性细杆理论已成功应用于DNA、蛋白质、神经元、弹簧、树枝和水下通信电缆等多种一维结构的建模。

5.1 曲线与细杆的运动学

5.1.1 曲线与Frenet标架

空间曲线的参数化：设 r(S, T) 为定义在区间 [0, L] × [T₀, T₁] 上的 C³类（对S）和 C²类（对T）的空间曲线。若曲线上每一点的两个S导数线性无关，则称该曲线为正则曲线。

弧长定义：弧长 s 作为一种方便的参数化方式定义为：

s(S,t) = ∫₀ˢ |∂r(σ,T)/∂σ| dσ

其中 |r| = √(r·r) 为通常的欧几里得范数。

单位切向量：τ = ∂r/∂s ≡ r′，即沿弧长s的切向量。

Frenet标架 {τ, ν, β} 由三个正交单位向量组成：切向量τ、法向量ν和次法向量β。

  • 曲率κ：衡量切向量沿弧长的转动速率，κ = |∂τ/∂s| = |r″|
  • 挠率τ（torsion）：衡量Frenet标架绕切线τ旋转的度量
  • 次法向量：β = τ × ν

Frenet方程（闭合的一阶微分方程组）：

∂τ/∂s = κ ν
∂ν/∂s = τ β − κ τ
∂β/∂s = −τ ν

空间曲线基本定理：给定可微函数 κ(s) > 0 和 τ(s)，存在一条以s为弧长、κ(s)为曲率、τ(s)为挠率的唯一曲线（在整体旋转、平移和定向意义上）。

5.1.2 细杆与广义标架

细杆的几何定义：细杆由中心线 r(S, T) 和两个额外的单位正交向量场 {d₁(S, T), d₂(S, T)} 组成。这两个单位向量表示材料截面在S处的方向。约定 d₃ = d₁ × d₂，构成右手正交基。

拉伸向量：v = ∂r/∂S，其范数 α = |v| = |v| 表示相对于参考配置的拉伸。

Darboux向量u：描述标架沿弧长的演化： 

∂dᵢ/∂S = u × dᵢ, i = 1, 2, 3

  • u₁, u₂：与弯曲相关
  • u₃：与扭转相关（绕d₃轴的旋转）

自旋向量w： 

∂dᵢ/∂T = w × dᵢ, i = 1, 2, 3

矩阵形式：引入正交矩阵 D = [d₁, d₂, d₃]，Darboux矩阵U和自旋矩阵W为反对称矩阵：

U = [ 0  -u₃  u₂ ]
    [ u₃  0  -u₁ ]
    [-u₂  u₁   0 ]

W = [ 0  -w₃  w₂ ]
    [ w₃  0  -w₁ ]
    [-w₂  w₁   0 ]

兼容性条件（Lax对结构）： 

∂U/∂T − ∂W/∂s = UW − WU

5.1.3 不可伸长、无剪切细杆

当 v₁ = v₂ = 0, v₃ = 1 时，α = 1，s = S，无拉伸且无剪切。此时d₃为切向量，d₁和d₂位于法平面内，通过登记角φ与Frenet标架关联：

d₁ = ν cos φ + β sin φ
d₂ = −ν sin φ + β cos φ

Darboux向量为： 

u = (κ sin φ, κ cos φ, τ + ∂φ/∂S)

过量扭转（∂φ/∂S）：仅是细杆的性质，表示局部基相对于Frenet标架的旋转。

  • Frenet细杆：∂φ/∂S = 0，φ为常数
  • 扭转密度u₃：同时描述空间曲线和细杆的性质

5.2 弹性细杆的力学

5.2.1 线性动量平衡

考虑位于截面S(S)和S(S−ΔS)之间、长度为ΔS的一段细杆。线性动量守恒给出：

∂n/∂S + f = ρA · ∂²r/∂T²

其中： - n(S,T)：截面S处的合内力（由相邻截面作用） - f(S)：作用于截面的体力（单位参考长度） - ρA：线密度（单位参考长度的质量）

5.2.2 角动量平衡

关于原点的总力矩和总角动量守恒。经推导（利用d₁和d₂沿截面主惯性方向的选择），得到：

∂m/∂S + ∂r/∂S × n + l = ρI₂ d₁ × ∂²d₁/∂T² + ρI₁ d₂ × ∂²d₂/∂T²

其中： - m(S,T)：截面S处的合内力矩 - l(S)：外部体力偶 - I₁, I₂：截面的二次面积矩（惯性矩）

5.2.3 细杆局部力学总结

基本平衡方程：

∂n/∂S + f = ρA · ∂²r/∂T²                          (线性动量)
∂m/∂S + ∂r/∂S × n + l = ρI₂ d₁ × ∂²d₁/∂T² + ρI₁ d₂ × ∂²d₂/∂T²  (角动量)

5.3 弹性细杆的本构方程

5.3.1 可伸长、可剪切弹性细杆

假设细杆是超弹性的：存在应变能密度函数 W = W(u−û, v−v̂, S)，使得：

m = ∂W/∂(u−û)
n = ∂W/∂(v−v̂)

其中 v̂, û 是无应力参考配置的应变。

5.3.2 不可伸长、无剪切细杆

取 v = d₃，S = s（弧长），应变能仅是(u−û)的函数：

m = ∂W/∂(u−û)

对于二次应变能 W = yᵀKy：

m = K(u − û),  K = [K₁₁ K₁₂ K₁₃; K₁₂ K₂₂ K₂₃; K₁₃ K₂₃ K₃₃]

简化对角形式（最常用的情况）：

W = K₁(u₁−û₁)² + K₂(u₂−û₂)² + K₃(u₃−û₃)²

其中： - K₁, K₂：主弯曲刚度 - K₃：扭转刚度

线性弹性框架下：

m = E·I₁(u₁−û₁)d₁ + E·I₂(u₂−û₂)d₂ + μJ(u₃−û₃)d₃
  • E：弹性模量（Young's modulus）
  • μ：剪切模量
  • J：取决于截面形状的参数

两个重要无量纲参数：

a = I₂/I₁  （截面不对称性度量，0 < a ≤ 1）
b = μJ/(EI₁) = J/[2I₁(1+σ)]  （标度扭转刚度）
  • σ：泊松比（通常 0 ≤ σ ≤ 1/2）
  • 对于不可压缩材料，b ∈ [2/3, 1]

5.3.3 各向同性、不可伸长但可伸长细杆

对于生长细杆的建模，考虑无剪切但可伸长的情形：

m = E·I₁(u₁−û₁)d₁ + E·I₂(u₂−û₂)d₂ + μJ(u₃−û₃)d₃
n₃ = E·A(α − 1)        （胡克定律）

其中 α = ∂s/∂S 为弹性拉伸。

5.4 标度（Scaling）

为简化系统，采用以下标度方案（对于 I₁ = I₂ = I 的情形）：

[M] = ρAI
[T]² = ρI/AE

t → t·√(ρI/AE)
s → s·√(I/A)
n → n·AE
m → m·√(EI/A)
u → u·√(A/I)
w → w·√(AE/(ρI))

弯曲与扭转刚度比（描述圆截面的单参数）：

ℰ = μJ/(EI₁) = 1/(1+σ) ∈ [2/3, 1]

5.5 弯曲与扭转刚度

截面参数J

在翘曲函数 φ(x₁, x₂) 的框架下，参数J由下式确定：

J = I₁ + I₂ + ∫∫ [x₁∂φ/∂x₂ − x₂∂φ/∂x₁] dx₁dx₂

翘曲函数满足 Neumann 边值问题（拉普拉斯方程）。

应力函数ψ

对于单连通截面，J可表示为：

J = 2∫∫ ψ dx₁dx₂ − I₁ − I₂

其中ψ是φ的共轭调和函数，满足狄利克雷问题。

常见截面形状

椭圆截面（半轴 A < B）：

I₁ = πAB³/4,  I₂ = πA³B/4
J = πA³B³/(A²+B²) = 4I₁I₂/(I₁+I₂)

圆截面（半径 R）：

I₁ = I₂ = πR⁴/4 = J/2
a = 1,  b = ℰ

矩形截面（边长 A < B）：

I₁ = AB³/12,  I₂ = A³B/12
J ≈ A³B/3 − 64A⁴/π⁵ tanh(πB/(2A))  （误差 ±0.502%）

5.6 基尔霍夫弹性细杆模型总结

最常用的模型：无剪切 + 不可伸长 + 对角本构关系

三大类方程：

1. 运动学方程

∂r/∂s = d₃
∂D/∂s = DU
∂D/∂t = DW

2. 力学平衡方程

∂n/∂s + f = ρA·∂²r/∂t²
∂m/∂s + ∂r/∂s × n + l = ρI₂ d₁ × ∂²d₁/∂t² + ρI₁ d₂ × ∂²d₂/∂t²

3. 本构方程

m = E·I₁(u₁−û₁)d₁ + E·I₂(u₂−û₂)d₂ + μJ(u₃−û₃)d₃

模型特点： 1. 材料由三个刚度（EI₁, EI₂, μJ）、密度ρ和截面几何（A, L, û）描述 2. 系统有15个变量，方程组完备可解 3. 旋转惯性项（方程右侧）通常在典型时间尺度分析中被忽略 4. û = 0 表示自然直细杆

基尔霍夫动力学类比（Kirchhoff Kinetic Analogy）

重要发现：平衡方程在形式上等价于定点刚体（陀螺）的欧拉方程！

  • 陀螺轴 ↔ 细杆切向量
  • 时间 ↔ 弧长
  • 当 û₂ = 0，方程对应拉格朗日情况（可积系统）

例子：拉格朗日陀螺轴的运动等价于横向各向同性不可伸长细杆的平衡构型。

5.7 螺旋细杆示例

5.7.1 螺旋线几何

螺旋线参数方程（绕 e₁ 轴，柱半径 R）：

r = Pδse₁ + Rcos(δs)e₂ + Rsin(δs)e₃

其中 δ = ±1/√(P²+R²)，+对应右手螺旋。

  • 螺距（每圈高度）：h = 2π|P|
  • 每圈弧长：l = 2π/|δ|

Frenet标架：

τ(s) = Pδe₁ − Rδsin(δs)e₂ + Rδcos(δs)e₃
ν(s) = −cos(δs)e₂ − sin(δs)e₃
β(s) = −Rδe₁ + Pδsin(δs)e₂ − Pδcos(δs)e₃

曲率和挠率（恒定）：

κ = Rδ² = R/(P²+R²)
τ = Pδ² = P/(P²+R²)

扭曲螺旋细杆（Frenet标架旋转角 γ）：

d₁ = νcos(γs) + βsin(γs)
d₂ = −νsin(γs) + βcos(γs)
d₃ = τ

Darboux向量：

u = (κsin(γs), κcos(γs), τ+γ)

5.7.2 螺旋平衡解

平衡方程（无体力 f = l = 0）：

n′ = 0
m′ + d₃ × n = 0

三类特殊螺旋：

  1. 自然直细杆（κ̂ = τ̂ = 0）：存在任意 ω 的螺旋解
  2. 直细杆（u = 0）：中心线为直线
  3. 零力螺旋（n = 0）：m 平行于 u

一般螺旋平衡（κ ≠ 0, ω = 0）：

n = [K₃κ(τ−τ̂) − K₁τ(κ−κ̂)] (0, 1, τ/κ)
m = (0, K₁(κ−κ̂), K₃(τ−τ̂))

力-矩螺旋（沿螺旋轴的组合）：

M = n·e₁ = δ⁻¹[K₁κ(κ−κ̂) + K₃τ(τ−τ̂)]
N = m·e₁ = δκ⁻¹[K₃κ(τ−τ̂) − K₁τ(κ−κ̂)]

这些即经典的螺旋弹簧公式（Thomson & Tait, 1867）。

5.7.3 过卷或欠卷螺旋

应用：DNA单分子实验（光镊/磁镊拉伸）。

问题：沿轴向拉伸螺旋弹簧时，它是展开还是进一步卷绕？

结论（基于螺旋弹簧公式）：

  • 对于金属弹簧（圆形截面，ℰ = μJ/(EI₁) ∈ [2/3, 1]）
  • 随着轴向张力 N > 0 增加，先过卷（overwinding），后展开（unwinding）
  • 初始过卷条件：ℰ < 1

线性近似（小拉伸）：

θ = N(K₁−K₃)κ̂²τ̂ / [K₁K₃√(κ̂²+τ̂²)] + O(N²)

当初始状态 κ̂²τ̂ ≠ 0 时，若 K₃/K₁ < 1，则 θ > 0（过卷）。

几何解释：M = 0 的平衡在曲率-挠率平面上位于椭圆上：

κ(κ−κ̂) + ℰ·τ(τ−τ̂) = 0

5.8 平面弹性细杆：伯努利-欧拉方程

问题设定

  • 平面、无剪切、不可伸长
  • 圆形截面
  • 自然直（无内禀曲率）
  • 无体力、无体力偶
  • 无扭转、无扭转

此时 Darboux 向量：u = (0, κ, 0)

角度参数化

设细杆位于 e₁-e₂ 平面，切向量与 e₁ 轴夹角为 θ：

τ = d₃ = cosθ·e₁ + sinθ·e₂
κ = θ′

运动方程

F′ + f = ρA·ẍ
G′ + g = ρA·ÿ
EI·θ″ + Gcosθ − Fsinθ = ρI·θ̈

静力学解

静态、无体力：F 和 G 常数，令 G = 0，得：

θ″ + αsinθ = 0

其中 α = −F/(EI)，这就是摆的方程！

基尔霍夫类比：平面情况等价于摆的运动，这是拉格朗日陀螺类比的二维版本（Bernoulli和Euler已知）。

本章重要公式汇总

公式 内容
(5.8)-(5.10) Frenet方程
(5.12)-(5.13) 广义标架演化方程
(5.41)-(5.42) 线性/角动量平衡
(5.81) 对角本构方程
(5.92)-(5.93) 细杆平衡方程
(5.115) 平面弹性细杆的摆方程

物理意义与应用

  1. DNA建模：螺旋结构和扭曲响应
  2. 蛋白质折叠：α-螺旋的几何
  3. 植物生长：细胞壁各向异性与手性
  4. 细菌鞭毛：螺旋形运动
  5. 神经轴突：力学稳定性
  6. 工程应用：弹簧设计、圆柱结构

关键概念中英对照

中文 English
曲率 Curvature (κ)
挠率 Torsion (τ)
Frenet标架 Frenet frame
Darboux向量 Darboux vector
弯矩 Bending moment
扭矩 Twisting moment
本构方程 Constitutive equation
超弹性 Hyperelastic
弧长 Arc length
螺旋弹簧 Helical spring
欧拉方程 Euler equations
摆方程 Pendulum equation

读书笔记完成日期：2026年5月9日