曲面上边界向XY平面投影：不可拉伸薄片的圆投影形状求解

这是个挺有意思的几何问题，核心抓住薄片不可拉伸/收缩这个等距约束就好办了，我来一步步拆解怎么求解：

1. 明确核心约束：等距映射 #

首先，薄片的内部坐标$(u,v)$到曲面$(x,y,z(x,y))$的映射必须是等距映射——因为薄片在平面内任意方向都不能拉伸或收缩，意味着薄片上任意两点的弧长，和它们映射到曲面上对应点的弧长完全相等。

• 薄片上的线元（微小弧长的平方）是：du² + dv²
• 曲面的第一基本形式（曲面上的线元）是：E dx² + 2F dx dy + G dy²，其中：
  • $E = 1 + \left(\frac{\partial z}{\partial x}\right)^2$
  • $F = \left(\frac{\partial z}{\partial x}\right)\left(\frac{\partial z}{\partial y}\right)$
  • $G = 1 + \left(\frac{\partial z}{\partial y}\right)^2$
• 等距约束的数学表达就是：du² + dv² = E dx² + 2F dx dy + G dy²

2. 建立圆的方程关联 #

假设薄片上的圆在内部坐标系$(u,v)$下的标准方程是：

$(u - u_0)^2 + (v - v_0)^2 = r^2$
（其中$(u_0,v_0)$是圆心，$r$是圆的半径）

我们的目标是找到这个圆映射到曲面上后，在$xy$平面的投影曲线$(x(\theta), y(\theta))$，需要同时满足圆的方程和等距约束。

3. 具体求解方法 #

方法一：参数化求解 #

我们可以先给薄片上的圆做参数化：

• 令$u = u_0 + r\cos\theta$，$v = v_0 + r\sin\theta$，其中$\theta \in [0, 2\pi)$

接下来利用等距约束的弧长相等条件：薄片上圆的微小弧长是r dθ，对应曲面上曲线的微小弧长必须等于它，即：
$$\sqrt{E \left(\frac{dx}{d\theta}\right)^2 + 2F \left(\frac{dx}{d\theta}\right)\left(\frac{dy}{d\theta}\right) + G \left(\frac{dy}{d\theta}\right)^2} = r$$

这是一个关于$x(\theta)$和$y(\theta)$的微分方程，结合初始条件（比如圆心对应的投影点$(x_0,y_0)$），就可以解出投影曲线的参数方程。

方法二：等距映射反推 #

如果能先求出从薄片$(u,v)$到曲面$(x,y,z(x,y))$的等距映射函数$x = x(u,v)$，$y = y(u,v)$，那么直接把圆的方程$(u - u_0)^2 + (v - v_0)^2 = r^2$代入，消去参数$u,v$就能得到$xy$平面上的投影曲线方程。

求解等距映射可以通过雅可比矩阵的约束：映射的雅可比矩阵$J = \begin{pmatrix} \partial x/\partial u & \partial x/\partial v \ \partial y/\partial u & \partial y/\partial v \end{pmatrix}$需要满足：
$$J^T \begin{pmatrix} E & F \ F & G \end{pmatrix} J = \begin{pmatrix} 1 & 0 \ 0 & 1 \end{pmatrix}$$

4. 特殊情况简化 #

• 如果曲面是平面（$z$为常数），等距映射就是刚性平移/旋转，投影自然还是圆，符合直觉。
• 如果曲面是可展曲面（比如圆柱面、圆锥面），因为可展曲面可以等距展开成平面，我们可以先把曲面展开成平面，把薄片上的圆匹配到展开平面上，再映射回曲面，最后投影到$xy$平面，这个过程会比一般曲面简单很多。

内容的提问来源于stack exchange，提问作者user3433489