3D空间中已知弧端点与拱高，如何求解圆弧圆心？

问题背景 #

已知3D空间中一段圆弧的两个端点坐标、拱高（弦中点到弧中点的距离），且弧始终位于圆心所在平面的上方，需求解该圆弧对应的完整圆的圆心。此前套用网上的2D公式计算结果不符合预期，以下结合示例说明错误点与正确解法。

示例中的错误分析 #

给定示例参数：

• P1=(0,0,500)，P2=(500,0,0)，预期圆心为(0,0,0)，半径500，拱高h=146.42

原计算存在以下错误：

1. 弦长公式书写错误：
   原公式写成SQRT((p1.X + p2.X)^2 + ...)，正确应为坐标差的平方和：

   弦长L = SQRT((p2.X-p1.X)^2 + (p2.Y-p1.Y)^2 + (p2.Z-p1.Z)^2)
   

   本次计算结果碰巧正确，但逻辑错误。
2. 弦中点Zm公式书写错误：
   原公式误用X坐标计算Z中点，正确应为Zm=(p1.Z + p2.Z)/2，本次结果正确但公式错误。
3. 圆心Y坐标代入错误：
   计算时将Ym=0误写为250，导致计算逻辑混乱。
4. 核心错误：套用2D公式到3D场景：
   2D中圆心仅需在弦的垂直方向移动，但3D空间中，圆心位于弦的垂直平分面上，需额外确定圆弧所在平面的法向量，才能锁定唯一方向。

正确求解步骤 #

步骤1：计算基础参数 #

设P1=(x1,y1,z1)，P2=(x2,y2,z2)，拱高h：

• 弦中点M：

  Xm=(x1+x2)/2，Ym=(y1+y2)/2，Zm=(z1+z2)/2
  

• 弦向量V：

  V=(x2-x1, y2-y1, z2-z1)
  

• 弦长L：

  L=√[(x2-x1)² + (y2-y1)² + (z2-z1)²]
  

• 圆弧半径r（由几何关系推导，r²=(L/2)²+(r-h)²）：

  r=(L²)/(8*h) + h/2
  

• 圆心到弦中点的距离d：

  d=r-h
  

步骤2：确定圆弧所在平面的法向量 #

3D空间中需额外指定圆弧所在平面的法向量n（如示例中圆弧在XOZ平面，法向量为(0,1,0)），否则圆心不唯一。

步骤3：计算垂直于弦的方向向量 #

在圆弧所在平面内，计算垂直于弦向量V的单位方向向量u：

1. 计算法向量与弦向量的叉乘，得到平面内垂直于V的向量：

   U = n × V （叉乘顺序决定方向，需根据弧的“上方”要求调整）
   

2. 单位化向量U：

   u = U / ||U|| （||U||为U的模长）
   

步骤4：计算圆心坐标 #

根据弧的方向（位于圆心所在平面上方），圆心C为：

C = M - d*u （或M + d*u，需验证方向是否符合弧的位置要求）

示例验证（正确计算） #

代入示例参数：

• M=(250,0,250)，V=(500,0,-500)，n=(0,1,0)
• U = (0,1,0) × (500,0,-500) = (-500,0,-500)
• ||U||=707.11，u=(-500,0,-500)/707.11≈(-0.7071,0,-0.7071)
• d=500-146.42=353.58
• C=(250,0,250) - 353.58*(-0.7071,0,-0.7071) = (250-250,0,250-250)=(0,0,0)，与预期一致。

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