---
title: "同样一片膜可以有几种稳定形状"
description: "在固定面积与体积下比较轴对称变分分支，学习约束寻优、曲率积分、收敛与形状反演的不可辨识性。"
canonical_url: "https://sandboxphysics.com/zh/stories/membrane-mechanics/shape-minima/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/membrane-mechanics/shape-minima/"
last_modified: "2026-10-28"
---

[膜与囊泡力学](https://sandboxphysics.com/zh/stories/membrane-mechanics/)

L3 · 第 3 章

# 同样一片膜可以有几种稳定形状

Sandbox Physics · 2026-10-28

相同面积和体积可以容许多种形状。膜会在允许的形状中降低能量，实验却可能停在某个局部极小值附近。要判断一张长形囊泡照片说明了什么，先要分清几何约束、能量极小值和拍照之前的历史。

## 开始前需要什么

前两章的几何与 Helfrich 能。需要理解一元函数的极小值；这里介绍把曲面限制到少量参数后的数值寻优。完整形状方程还需要变分法和微分方程。

- 说明固定面积与体积如何减少变分参数，并设计约束残差与收敛检查。
- 解释为什么一张平衡形状图不能单独测出弯曲刚度。

## 先选一族能算清的曲面

给囊泡选一根对称轴，用极角 $\theta$ 从北极往南极走。径向距离 $r(\theta)$ 决定曲面。我们用两个低阶 Legendre 多项式描述它：二阶项主要区分长与扁，四阶项允许轮廓进一步调整。指数让半径始终为正。

整体尺度 $s$ 用来满足面积，系数 $a$ 和 $b$ 改变形状。这个有限变分族光滑、轴对称且相对中心呈星形：每条射线只穿过膜一次。反复穿过同一射线的凹陷、极窄颈和不对称形状超出了这个范围；受限的长形和扁形仍能提供可控的能量比较。

$$
r(\theta)=s\exp\!\left[aP_2(\cos\theta)+bP_4(\cos\theta)\right]
$$

这里 $P_2(x)=(3x^2-1)/2$，$P_4(x)=(35x^4-30x^2+3)/8$。参数无量纲，尺度 $s$ 以长度计。

来源：U. Seifert · W. Helfrich

## 把面积和体积真的固定住

先临时令整体尺度为一，沿极角积分求面积和体积。再乘一个尺度，让面积恰好等于目标值。整体缩放不会改变约化体积，因此剩下的工作是找到满足目标约化体积的形状系数。

固定一个 $b$，沿正、负两侧分别寻找 $a$。二分法每次取一个区间的中点，判断约化体积位于目标哪边，再保留夹住解的半个区间。找到两支允许形状后，才比较它们的弯曲能；不能拿体积不同的两个漂亮轮廓直接排序。

最后在四阶系数的有限范围内寻找每一支的最低能量。画出的结果是这族曲面、这个搜索范围内的受限极小值。真实囊泡的全局最低能量还要允许更多形状，并检查自接触与稳定性。

$$
A=2\pi\int_0^\pi r\sin\theta\sqrt{r^2+r_\theta^2}\,\mathrm d\theta,\qquad V=\frac{2\pi}{3}\int_0^\pi r^3\sin\theta\,\mathrm d\theta
$$

面积来自小圆环的周长乘以母线长度；体积来自每个立体角内的径向体积。

$$
s=\sqrt{\frac{A_{\rm target}}{A_{s=1}}},\qquad v(a,b)=v_{\rm target}
$$

第一式消去尺度，第二式消去一个形状自由度。剩下的参数沿受限曲线变化。

## 完整例题：把约束残差换成物理大小

取 $R_A=10\,\mu\mathrm m$、$v=0.80$。目标面积为 $1.26\times10^3\,\mu\mathrm m^2$，目标体积约为 $3.35\times10^3\,\mu\mathrm m^3$。假设某个候选形状的相对体积残差为 $10^{-4}$，对应水量误差约 $0.335\,\mu\mathrm m^3$。

这一步先检查它是不是同一个问题的解。然后把角向积分节点加倍，比较能量、约化体积和高斯曲率积分。数值不断接近同一个答案，才有理由认为差别来自模型，而不是轮廓离散得太粗。

球面还给出独立标尺：零自发曲率时，平均曲率能除以 $\kappa$ 应接近 $8\pi$，高斯曲率积分应接近 $4\pi$。这两个检查分别检验局部曲率与闭合表面的整体几何。

## 形状照片缺少哪个尺度

假设能量只有平均曲率项，面积、体积和自发曲率都已给定。把 $\kappa$ 乘以任意正数，会把所有候选能量同时乘大，但不会移动极小值。即使轮廓测得完美，形状也无法单独确定弯曲刚度。

若再把未知张力、压力或蛋白分布一起拟合，可互相补偿的参数会更多。拟合程序返回一个数字，只表示它在当前假设下找到了一个答案。测热涨落、施加已知力或加入独立压力测量，才能补上形状里缺少的信息。

这也解释了为什么换一套材料参数有时几乎看不出轮廓差异。更可靠的反演方案先找出观测对哪些参数敏感，再选择第二种实验。下一段课程就用热运动提供一个独立能量标尺。

$$
F_b=\kappa\,\mathcal E[\text{shape}],\qquad \arg\min F_b=\arg\min\mathcal E\quad(\kappa>0)
$$

这个结论限于刚度作为整体正因子的模型；加入其他独立能量尺度后，需要重新分析。

## 把推理带进实验

### [比较同一组约束下的两个分支](https://sandboxphysics.com/zh/labs/vesicle-shapes/)

1. 选择中等放瘪的囊泡，先预测长形和扁形是否能同时满足约束。
2. 比较两支的体积、面积和能量，再改变自发曲率，看能量排序何时发生变化。
3. 打开极区曲率偏好，解释轮廓变化。此处极区分布事先指定；思考若蛋白能够移动，还需增加什么方程。

## 练习与解释性答案

1. 两个候选形状的能量相差很小，但体积残差比能量差的相对量级还大。应先改进什么？

   展开解释

   先收紧几何约束并提高积分精度。体积不同会把不同物理问题的解混在一起；在约束和离散误差受控后，才讨论小能量差是否稳定。

2. 一支在有限参数族里达到极小值，能称为稳定囊泡吗？

   展开解释

   可以称为该参数族内的受限极小值。还要允许被排除的扰动，检查二阶能量变化、自接触和模型适用范围，才判断更一般的稳定性。

3. 设计一个补充实验，让只看形状时无法确定的刚度变得可测。

   展开解释

   例如在已知温度下记录膜涨落，再用相机响应校正后的谱拟合刚度和张力；或同时测膜管半径与拉力。两者都引入了独立的能量或力标尺。

[上一章：弯一片膜要付出多少能量](https://sandboxphysics.com/zh/stories/membrane-mechanics/bending-energy/) [下一章：膜的抖动怎样显出弯曲刚度](https://sandboxphysics.com/zh/stories/membrane-mechanics/thermal-modes/)

## 来源与复算范围

1. [W. Helfrich, Elastic Properties of Lipid Bilayers: Theory and Possible Experiments, Z. Naturforsch. C 28, 693–703 (1973)](https://doi.org/10.1515/znc-1973-11-1209)

   曲率弹性原论文，693—703 页。课程采用局域二次曲率能，并另行声明球面为正的曲率约定；球面、膜管和小斜率例题在这里推导，没有复制原图或实测材料参数。

2. [U. Seifert, K. Berndl and R. Lipowsky, Shape transformations of vesicles, Physical Review A 44, 1182–1202 (1991)](https://doi.org/10.1103/PhysRevA.44.1182)

   已发表论文：自发曲率模型、双层耦合模型与轴对称形状分支。配套实验在有限径向变分族中固定面积和体积寻优；完整相图、极窄颈极限和口形囊泡不属于这次复算。