偏置圆柱体相贯线的空间坐标计算

两个圆柱形管的相贯存在4种形式，正交、偏置、斜交和偏置斜交。

本文仅讨论偏置情况下的计算方法。

建立如下坐标系，圆柱1的轴线为Z轴，圆柱2的轴线与Y轴平行，通过右手法则确定X轴，圆柱2的轴线在XOY平面上。圆柱1半径为a，圆柱2半径为b，圆柱2的轴线相较于Y轴偏移距离为ΔX，其中b > a。

采用极坐标系表示两个圆柱面和法线矢量：

${\sum }_1$： $x_1=a\cos \theta$， $y_1=a\sin \theta$，$z_1=u$；
法向量：$n_{x1}=\cos \theta ,n_{y1}=\sin \theta ,n_{z1}=0$

${\sum }_2$：$x_2=\Delta x+b\cos \phi$，$y_2=v$，$z_2=b\sin \phi$
法向量：$n_{x2}=\cos \phi ,n_{y2}=0,n_{z2}=\sin \phi$

相贯线同时位于两个圆柱面上，即满足
$$x_1=x_2,y_1=y_2,z_1=z_2$$

由于a, b, $\Delta x$为已知，且a < b，则$\theta$的取值范围为（$0-2\pi$），根据需要的点距可以选择不同的步长。

可解出相贯线的空间坐标关于$\theta$的方程如下：

$x=a\cos \theta$
$y=a\sin \theta$
$z=\sqrt{b^2-(x-\Delta x{)}^2}$

用python实现如下：

import math
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

# 圆柱1半径：radius_1;  偏置距离：offset
# 圆柱2半径：rasius_2; （直径大）
# 点距： theta；（按弧度）

class IntersectingLine():
    def __init__(self):
        self.radius_1 = 8.0
        self.radius_2 = 30.0
        self.center_offset = 10.0
        self.step_theta = 0.04
	
	# 计算空间坐标，保存在 ndarray中
    def calculate_coordinates(self, r1, r2, offset, theta):
        shape = ((int)(math.pi * 2 / theta), 3)
        coords = np.zeros(shape, dtype=float)
        for i in range(shape[0]):
            coords[i][0] = r1 * math.cos(i*theta)
            coords[i][1] = r1 * math.sin(i*theta)
            coords[i][2] = math.sqrt(r2*r2 - math.pow((coords[i][0] - offset), 2))
        return coords

    # 根据坐标点绘制三维散点图
    def plot_coordinates(self, coordinates):
        fig = plt.figure()
        ax = Axes3D(fig)
        ax.view_init(elev=45, azim=135)
        ax.scatter3D(coordinates[:,0], coordinates[:,1], coordinates[:,2], cmap='Blues')
        plt.show()

绘制结果如下：