---
title: "弯一片膜要付出多少能量"
description: "推导球面和圆柱的 Helfrich 弯曲能，分清弯曲刚度、张力、面积压缩模量及固定拓扑的高斯项。"
canonical_url: "https://sandboxphysics.com/zh/stories/membrane-mechanics/bending-energy/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/membrane-mechanics/bending-energy/"
last_modified: "2026-10-28"
---

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

L2 · 第 2 章

# 弯一片膜要付出多少能量

Sandbox Physics · 2026-10-28

膜抵抗的是偏离偏爱曲率的弯曲。把每一小块偏离的平方加起来，就得到 Helfrich 弯曲能。它解释了一个初看奇怪的结果：没有自发曲率时，大球和小球的弯曲能相同。

## 开始前需要什么

第一章的主曲率、平均曲率、面积半径和约化体积。会计算平方；本章把积分解释为把小块能量相加，把导数解释为微小改变量的比值。

- 从局域能量密度推导球面的总弯曲能并检查量纲。
- 说明张力与面积压缩模量为什么有相同单位却代表不同的量。

## 从每一小块开始相加

先把膜视为厚度可忽略的光滑曲面。局部偏爱的总曲率记为 $c_0$，单位是长度倒数。膜偏离它时，能量密度按 $(2H-c_0)^2$ 增长。前面的弯曲刚度 $\kappa$ 以能量计，越大就越难弯。

积分符号表示把整个表面划成越来越小的面积块 $\mathrm dA$，计算每块的能量后相加。高斯项还有自己的系数 $\bar\kappa$。两种系数都属于曲率弹性；直接拉长膜的代价还需要另一项弹性能。

$$
F_b=\int_A\left[\frac{\kappa}{2}(2H-c_0)^2+\bar\kappa K\right]\,\mathrm dA
$$

局域曲率平方乘以面积后没有量纲，故两个刚度系数都具有能量单位。这里假定连续、各向同性、流体双层，忽略膜厚尺度上的细节。

来源：W. Helfrich

## 球越大，为什么能量未必越大

球上每一点的 $2H=2/R$ 相同，可以直接从积分中提出。面积 $4\pi R^2$ 随半径平方增加，曲率平方却随半径平方减小，两者恰好抵消。先取 $c_0=0$ 并只计算平均曲率项，就得到 $8\pi\kappa$。

允许自发曲率后，球面结果改为 $2\pi\kappa(2-c_0R)^2$。此时半径重新进入能量，膜存在偏爱的弯曲尺度。翻转曲率符号约定时要同时翻转 $c_0$，物理预测才保持相同。

对于无边界且拓扑固定的闭合曲面，Gauss–Bonnet 定理给出 $\int K\,\mathrm dA=4\pi(1-g)$，$g$ 是孔洞数。球拓扑的高斯项是常数；比较同一拓扑下的形状时它不影响排序。融合或穿孔会改变这个前提。

$$
F_{b,\rm sphere}^{(H)}=2\pi\kappa(2-c_0R)^2,\qquad F_{b,\rm sphere}^{(H)}\big|_{c_0=0}=8\pi\kappa
$$

上标只标识平均曲率部分。球拓扑的高斯贡献还需另加一个与形状无关的常数。

来源：U. Seifert

## 张力和拉伸刚度扮演不同角色

张力是增加少量面积需要付出的能量除以面积，也可理解为沿切口每单位长度的拉力。若把总面积严格固定，张力可以作为保证面积不变的拉格朗日乘子：一个在求极值时使用的辅助量，其数值由约束决定。

面积压缩模量 $K_A$ 描述张力随相对面积伸长变化的快慢。设轻松状态的面积为 $A_0$，小应变为 $\alpha=(A-A_0)/A_0$，二次拉伸能给出下面两式。张力与 $K_A$ 都以能量除以面积计，但前者是当前的应力，后者是应力对应变的斜率。

囊泡常先用展开热皱褶的方式增加可见面积，随后才明显拉开脂质。实验中用哪一种面积定义应变，会影响解释；解释脂质分子之间的距离变化时，需要先分开皱褶展开和真实面积拉伸。

$$
F_s=\frac{K_AA_0}{2}\alpha^2,\qquad \sigma=\frac{\partial F_s}{\partial A}=K_A\alpha
$$

求导就是在同一材料状态附近比较一小段面积变化和能量变化。此式限于小应变的弹性拉伸模型。

## 完整例题：把能量换成热能

取教学参数 $\kappa=20k_BT$，没有自发曲率。球面平均曲率能为 $160\pi k_BT\simeq503k_BT$。在 $T=298\,\mathrm K$ 时，$k_BT\simeq4.11\,\mathrm{pN\,nm}$，总能约为 $2.07\times10^3\,\mathrm{pN\,nm}$。

把球半径从五微米增到十微米，上面这个结果保持不变。若膜从来是平的，要封成球仍须付出这份弯曲代价。球在已闭合、面积与体积固定的其他形状中是否最低，则是另一个比较问题。

再设 $K_A=0.20\,\mathrm{N\,m^{-1}}$、$\alpha=0.01$，小应变模型给出 $\sigma=2.0\,\mathrm{mN\,m^{-1}}$。这些数值只用来演算量纲，并未标定某种脂质。

## 把推理带进实验

### [把球面作为能量标尺](https://sandboxphysics.com/zh/labs/vesicle-shapes/)

1. 把约化体积调到一、自发曲率与极区曲率设为零，预测能量除以弯曲刚度后的数值。
2. 比较数值积分与球面解析值，再降低约化体积，查看两个允许形状的能量。
3. 保持几何约束，只改变自发曲率，解释局部偏好怎样改变形状排序。

## 练习与解释性答案

1. 若球满足 $c_0R=2$，它的平均曲率弯曲能是多少？总能一定为零吗？

   展开解释

   平均曲率项为零，因为整张球面恰好达到偏爱曲率。总能仍可能含高斯项、拉伸项和外力功；必须先说明采用的能量与参考零点。

2. 没有自发曲率时，把同一形状整体缩小十倍，平均曲率弯曲能怎样变？

   展开解释

   保持不变。每处曲率增大十倍、平方增大一百倍，面积减到百分之一。这个缩放结论要求刚度不随位置和尺度改变，且连续膜模型仍适用。

3. 为什么比较固定面积的两种形状时，不能靠张力乘以面积来选胜者？

   展开解释

   两者面积相同，这一项相同。要比较的是弯曲能和其他确实随形状改变的项。求解时的面积乘子仍决定驻点条件与力学响应。

[上一章：一层柔软的膜为什么能撑出形状](https://sandboxphysics.com/zh/stories/membrane-mechanics/curvature-and-constraints/) [下一章：同样一片膜可以有几种稳定形状](https://sandboxphysics.com/zh/stories/membrane-mechanics/shape-minima/)

## 来源与复算范围

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)

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