跳转到内容
Yizi's Math Notes

第八讲 曲面的弯曲与悬链面

第八讲 曲面的弯曲与悬链面

Section titled “第八讲 曲面的弯曲与悬链面”

开场：从肥皂膜到工程奇迹

Section titled “开场：从肥皂膜到工程奇迹”

  • 想象一下：将两个金属圆环平行放置，然后蘸上肥皂水。当你轻轻分开它们时，中间会形成一个美丽的曲面——这就是悬链面。

  • 这个看似简单的形状，却蕴含着深刻的数学原理，并在现代工程中发挥着重要作用。

  • 今天我们的目标：理解曲面如何弯曲，并通过计算悬链面的曲率，揭示数学与工程的完美结合。

曲面弯曲的两把”尺子”

Section titled “曲面弯曲的两把”尺子””

第一把尺子：内在的测量

Section titled “第一把尺子：内在的测量”

对于参数化曲面 r⃗(u,v)\vec{r}(u,v)，第一基本形式告诉我们如何在曲面上测量距离：

I=Edu2+2Fdudv+Gdv2I = E du^2 + 2F du dv + G dv^2

其中：E=r⃗u⋅r⃗uE = \vec{r}_u \cdot \vec{r}_u, F=r⃗u⋅r⃗vF = \vec{r}_u \cdot \vec{r}_v, G=r⃗v⋅r⃗vG = \vec{r}_v \cdot \vec{r}_v

直观理解：这就像给生活在曲面上的”二维蚂蚁”一把尺子，它能测量曲面上的距离、角度和面积。

第二把尺子：外在的弯曲

Section titled “第二把尺子：外在的弯曲”

第二基本形式测量曲面在三维空间中的弯曲：

II=Ldu2+2Mdudv+Ndv2II = L du^2 + 2M du dv + N dv^2

其中：L=r⃗uu⋅n⃗L = \vec{r}_{uu} \cdot \vec{n}, M=r⃗uv⋅n⃗M = \vec{r}_{uv} \cdot \vec{n}, N=r⃗vv⋅n⃗N = \vec{r}_{vv} \cdot \vec{n}

直观理解：这测量的是法向量 n⃗\vec{n} 的变化速度，只有”三维神”才能感知到。

两种关键曲率

Section titled “两种关键曲率”
  • 高斯曲率：K=LN−M2EG−F2K = \frac{LN - M^2}{EG - F^2}（内在性质，弯曲纸张不会改变）
  • 平均曲率：H=EN−2FM+GL2(EG−F2)H = \frac{EN - 2FM + GL}{2(EG - F^2)}（外在性质，描述在空间中的弯曲）

曲率公式的来历

Section titled “曲率公式的来历”

  • 高斯曲率公式的由来：

    • 高斯在1827年发现，曲面上任意一点的高斯曲率等于该点两个主曲率的乘积：K=κ1κ2K = \kappa_1 \kappa_2
    • 主曲率是形状算子 S=I−1IIS = I^{-1}II 的特征值，而特征值的乘积等于矩阵的行列式
    • 因此：K=det⁡(S)=det⁡(I−1II)=det⁡(II)det⁡(I)=LN−M2EG−F2K = \det(S) = \det(I^{-1}II) = \frac{\det(II)}{\det(I)} = \frac{LN - M^2}{EG - F^2}
    • 这个公式的神奇之处在于：虽然定义用到了嵌入空间，但结果只依赖于第一基本形式！

  • 平均曲率公式的由来：

    • 平均曲率是两个主曲率的算术平均：H=κ1+κ22H = \frac{\kappa_1 + \kappa_2}{2}
    • 主曲率之和等于形状算子的迹（对角线元素之和）：κ1+κ2=tr(S)\kappa_1 + \kappa_2 = \text{tr}(S)
    • 形状算子的迹可以通过矩阵运算计算：tr(S)=tr(I−1II)\text{tr}(S) = \text{tr}(I^{-1}II)
    • 经过计算得到：H=12tr(S)=EN−2FM+GL2(EG−F2)H = \frac{1}{2}\text{tr}(S) = \frac{EN - 2FM + GL}{2(EG - F^2)}
    • 平均曲率描述了曲面”向外弯曲”的程度，是变分几何中的重要概念

主曲率：理解弯曲的钥匙

Section titled “主曲率：理解弯曲的钥匙”

主曲率的直观理解

Section titled “主曲率的直观理解”
  • 在曲面上任意一点，存在两个特殊方向——主方向
  • 沿这两个方向，曲面的弯曲达到极值：κ1\kappa_1（最大弯曲）和 κ2\kappa_2（最小弯曲）
  • 这两个主曲率完全决定了该点的弯曲性质

欧拉公式：预测任意方向的弯曲

Section titled “欧拉公式：预测任意方向的弯曲”

κn(θ)=κ1cos⁡2θ+κ2sin⁡2θ\boxed{\kappa_n(\theta) = \kappa_1 \cos^2\theta + \kappa_2 \sin^2\theta}

含义：知道了主曲率，就能预测曲面在任何方向 θ\theta 上的弯曲程度。

悬链面：自然界的完美设计

Section titled “悬链面：自然界的完美设计”

什么是悬链面？

Section titled “什么是悬链面？”
  • 定义：平均曲率 H=0H = 0 的曲面，即 κ1+κ2=0\kappa_1 + \kappa_2 = 0
  • 物理意义：表面张力最小的形状，肥皂膜的自然选择
  • 历史：1744年欧拉发现的第一个极小曲面

悬链面的数学描述

Section titled “悬链面的数学描述”

参数方程（取 a=1a = 1 简化计算）：

r⃗(u,v)=(cosh⁡vcos⁡u,cosh⁡vsin⁡u,v)\vec{r}(u,v) = (\cosh v \cos u, \cosh v \sin u, v)

其中 u∈[0,2π]u \in [0, 2\pi], v∈Rv \in \mathbb{R}。

悬链面曲率的完整计算

Section titled “悬链面曲率的完整计算”

步骤1：计算偏导数

Section titled “步骤1：计算偏导数”

r⃗u=(−cosh⁡vsin⁡u,cosh⁡vcos⁡u,0)r⃗v=(sinh⁡vcos⁡u,sinh⁡vsin⁡u,1)r⃗uu=(−cosh⁡vcos⁡u,−cosh⁡vsin⁡u,0)r⃗uv=(−sinh⁡vsin⁡u,sinh⁡vcos⁡u,0)r⃗vv=(cosh⁡vcos⁡u,cosh⁡vsin⁡u,0)\begin{aligned} \vec{r}_u &= (-\cosh v \sin u, \cosh v \cos u, 0) \\ \vec{r}_v &= (\sinh v \cos u, \sinh v \sin u, 1) \\ \vec{r}_{uu} &= (-\cosh v \cos u, -\cosh v \sin u, 0) \\ \vec{r}_{uv} &= (-\sinh v \sin u, \sinh v \cos u, 0) \\ \vec{r}_{vv} &= (\cosh v \cos u, \cosh v \sin u, 0) \end{aligned}

步骤2：第一基本形式

Section titled “步骤2：第一基本形式”

E=r⃗u⋅r⃗u=cosh⁡2vF=r⃗u⋅r⃗v=0G=r⃗v⋅r⃗v=sinh⁡2v+1=cosh⁡2v\begin{aligned} E &= \vec{r}_u \cdot \vec{r}_u = \cosh^2 v \\ F &= \vec{r}_u \cdot \vec{r}_v = 0 \\ G &= \vec{r}_v \cdot \vec{r}_v = \sinh^2 v + 1 = \cosh^2 v \end{aligned}

步骤3：法向量

Section titled “步骤3：法向量”

r⃗u×r⃗v=(cosh⁡vcos⁡u,cosh⁡vsin⁡u,−sinh⁡vcosh⁡v)∣r⃗u×r⃗v∣=cosh⁡2vn⃗=(cos⁡u,sin⁡u,−sinh⁡v)cosh⁡v\begin{aligned} \vec{r}_u \times \vec{r}_v &= (\cosh v \cos u, \cosh v \sin u, -\sinh v \cosh v) \\ |\vec{r}_u \times \vec{r}_v| &= \cosh^2 v \\ \vec{n} &= \frac{(\cos u, \sin u, -\sinh v)}{\cosh v} \end{aligned}

步骤4：第二基本形式

Section titled “步骤4：第二基本形式”

L=r⃗uu⋅n⃗=−1M=r⃗uv⋅n⃗=0N=r⃗vv⋅n⃗=1cosh⁡v\begin{aligned} L &= \vec{r}_{uu} \cdot \vec{n} = -1 \\ M &= \vec{r}_{uv} \cdot \vec{n} = 0 \\ N &= \vec{r}_{vv} \cdot \vec{n} = \frac{1}{\cosh v} \end{aligned}

步骤5：曲率结果

Section titled “步骤5：曲率结果”

K=LN−M2EG−F2=(−1)⋅1cosh⁡vcosh⁡4v=−1cosh⁡4v<0H=EN−2FM+GL2(EG−F2)=cosh⁡2v⋅1cosh⁡v+cosh⁡2v⋅(−1)2cosh⁡4v=0κ1=1cosh⁡2v,κ2=−1cosh⁡2v\begin{aligned} K &= \frac{LN - M^2}{EG - F^2} = \frac{(-1) \cdot \frac{1}{\cosh v}}{\cosh^4 v} = \frac{-1}{\cosh^4 v} < 0 \\ H &= \frac{EN - 2FM + GL}{2(EG - F^2)} = \frac{\cosh^2 v \cdot \frac{1}{\cosh v} + \cosh^2 v \cdot (-1)}{2\cosh^4 v} = 0 \\ \kappa_1 &= \frac{1}{\cosh^2 v}, \quad \kappa_2 = -\frac{1}{\cosh^2 v} \end{aligned}

结果解读

Section titled “结果解读”
  • 极小曲面：H=0H = 0, 表面张力最小
  • 马鞍形：K=−1cosh⁡4v<0K = \frac{-1}{\cosh^4 v} < 0, 每点都向相反方向弯曲
  • 完美平衡：κ1=1cosh⁡2v\kappa_1 = \frac{1}{\cosh^2 v}, κ2=−1cosh⁡2v\kappa_2 = -\frac{1}{\cosh^2 v}, 两个主方向弯曲程度相等但方向相反
  • 具体数值：
    • 在 v=0v = 0（腰部最细处）：cosh⁡0=1\cosh 0 = 1, 所以 κ1=1\kappa_1 = 1, κ2=−1\kappa_2 = -1, K=−1K = -1
    • 当 ∣v∣|v| 增大时，cosh⁡v\cosh v 增大，主曲率的绝对值减小，曲面变得”平坦”
    • 这解释了悬链面的形状：中间最”弯”，两端逐渐趋于平坦

悬链面的工程应用

Section titled “悬链面的工程应用”

建筑结构

Section titled “建筑结构”
  • 圣路易斯拱门：世界最高的拱门，采用悬链线设计原理
  • 薄壳屋顶：悬链面形状的屋顶能以最少材料承受最大荷载
  • 优势：应力分布均匀，避免应力集中，提高结构稳定性

现代工程创新

Section titled “现代工程创新”
  • 3D打印：利用悬链面设计轻量化结构
  • 航空航天：机翼设计中借鉴悬链面的优化特性
  • 海洋工程：海上平台支撑结构采用悬链面设计
  • 材料科学：泡沫材料中的气泡壁自然形成悬链面

未来展望

Section titled “未来展望”
  • 智能材料：能够自适应形成悬链面的材料
  • 生物医学：利用悬链面原理设计人工器官支架
  • 建筑革命：更大跨度、更轻重量的建筑结构

总结：数学之美与工程之用

Section titled “总结：数学之美与工程之用”
  • 数学的力量：一个简单的条件（H=0H = 0）产生了如此丰富的几何结构
  • 自然的智慧：肥皂膜”知道”如何找到最优形状
  • 工程的应用：理解曲率成为解决现实问题的强大工具
  • 未来的桥梁：数学与工程的完美结合将创造更美好的世界

课堂思考题

Section titled “课堂思考题”
  1. 为什么肥皂膜总是形成悬链面？这与表面张力有什么关系？
  2. 如果我们改变悬链面参数方程中的常数 aa, 曲率如何变化？
  3. 除了悬链面，你还能想到哪些自然界中的极小曲面？
  4. 在建筑设计中，什么情况下我们会选择悬链面结构？