¹浙江长征职业技术学院基础部，浙江 杭州

²内蒙古民族大学数理学院，内蒙古 通辽

收稿日期：2020年10月7日；录用日期：2020年10月21日；发布日期：2020年10月28日

摘要

本文构造了通过给定控制多边形指定顶点且与始末控制多边形的始末边相切的两段连续三次B样条曲线，并以此两段连续的三次B样条曲线为轴线的管道光滑拼接了轴线异面椭圆管道。该方法相较于三段连续的三次均匀B样条曲线为轴线的管道、插值于两端点并与控制多边形的始末边相切的三次有理B样条曲线为轴线的管道、两段连续二次B样条曲线为轴线的管道光滑拼接轴线异面管道，各有其优越性。具有理论意义和应用价值。

关键词

B样条，分段连续，光滑拼接，轴线异面，椭圆管道

Cubic B-Spline Curve of Passing Appointed Vertexes and Its Application

Genzhu Bai^1,2

¹Basis Department, Zhejiang Changzheng Vocational & Technical College, Hangzhou Zhejiang

²College of Mathematics and Physics, Inner Mongolia University for the Nationalities, Tongliao Inner Mongolia

Received: Oct. 7^th, 2020; accepted: Oct. 21^st, 2020; published: Oct. 28^th, 2020

ABSTRACT

In this paper, we construct continuous two-section of cubic B-spline curve which is tangent to the start and end points of the control polygon. And we construct a tube that smoothly blends two ellipse tubes whose axes non-coplanar with the two continuous cubic B-spline curves as its axes. The results show that the method is better than that in which tubes are smoothly blended between smoothly connected three-segment tubes with cubic uniform B-spline curves as its axis, and tubes are smoothly blended by tube whose axis is a cubic rational B-spline curve interpolated at both ends and tangent to the start and end edges of the control polygon, and tubes are smoothly blended between smoothly connected two-segment tubes and quadratic uniform B-spline curves as its axis, which has theoretical significance and application value.

Keywords:B-Spline, Piecewise Continuous, Smooth Blending, Non-Coplaner, Elliptic Tube

Copyright © 2020 by author(s) and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

1. 引言

轴线异面管道光滑拼接问题研究是计算几何领域中难点基础问题。我们提出基于轴线光滑拼接的轴线异面管道光滑拼接技术旨在构造更多的自由曲线来光滑拼接轴线异面管道的轴线，进而实现两个轴线异面的管道的光滑拼接。Bézier曲线具有很好的端点性质，但是不具有局部修改性质。人们尝试将具有很好的局部修改性质的B样条曲线也具有像Bézier曲线那样很好的端点性质。主要通过两个途径，一种是对于 $V_0,V_1,\cdots ,V_n$ 为控制多边形顶点的三次B样条曲线在原来控制多边形的始末两段再增加两个顶点 $V_{-1}$ 和 $V_{n+1}$，满足 $V_{-1}=2V_0-V_1$ 和 $V_{n+1}=2V_n-V_{n-1}$，使B样条曲线通过原来控制多边形的起点 $V_0$ 和终点 $V_n$ ；另一种是对于给定的以 $V_0,V_1,\cdots ,V_n$ 为顶点的控制多边形构造过指定顶点的三次B样条曲线，反算B样条曲线的实际控制顶点 $P_0,P_1,\cdots ,P_n$。由于B样条曲线是分段光滑曲线，增加顶点会增加分段曲线的段数。构造插值于两端点的B样条曲线，并要求分段光滑曲线的段数较少，可能由于约束条件较少，构造出来的B样条曲线的光顺性不够理想。我们在文献 [1] 中用增加顶点的方法构造了三段连续的三次B样条曲线，用于轴线光滑拼接的轴线异面管道的拼接。在文献 [2] 中，对于给定的控制顶点，尝试用一段三次B样条曲线光滑拼接轴线异面管道的轴线，并将其应该用于轴线异面管道光滑拼接，发现光顺性较差，改用一段三次有理B样条曲线，并将其应该用于轴线异面圆管道光滑拼接。在选取权因子调整轴线的形状时，两端点附近的切线发生急剧变化，致使不能应用于轴线异面椭圆管道的拼接中。而在文献 [3] 中，对于给定的一段三次B样条曲线的控制顶点 $V_0,V_1,V_2,V_3$，利用二次B样条曲线与控制多边形的两边相切并通过控制多边形中点这一特性，构造了两段连续的二次B样条曲线，并将其应用于轴线异面管道光滑拼接。由于二次B样条曲线是平面曲线，而轴线异面管道的轴线不在同一平面内，造成连接点处的切线不连续，在构造光滑拼接的椭圆管道时需要反复调整椭圆管道的长半轴和短半轴才能实现，且曲线只能过控制多边形的边 $V_1{V}_2$ 的中点。本文对于给定

控制多边形顶点 $V_0,V_1,V_2,V_3$，在 $V_1$ 和 $V_2$ 中间插入一个顶点 $V_{1\frac{1}{2}}$，构造两段连续的三次B样条曲线过指定顶点 $V_0,V_{1\frac{1}{2}}$ 和 $V_3$，与始末两条边相切，并将其应用于轴线异面椭圆管道拼接中。

2. 过指定顶点的三次B样条曲线

在轴线异面管道拼接中，假设管道轴线 $P_0{P}_1$ 和 $P_2{P}_3$ 是异面直线， $P_0{P}_1$ 位于空间直角坐标系XOY平面内且与Y轴平行， $P_2{P}_3$ 位于YOZ平面内且与Y轴相交，连接 $P_1$ 和 $P_2$ 两点构成以 $P_0,P_1,P_2,P_3$ 为控制顶点

的控制多边形。为了构造两段连续的三次B样条曲线，在 $P_1$ 和 $P_2$ 中间插入一个顶点 $P_{1\frac{1}{2}}$ 构成以 $P_0,P_1,P_{1\frac{1}{2}},P_2,P_3$ 为顶点的控制多边形。

定义1设

$r\left(t\right)={\displaystyle \underset{i=0}{\overset{3}{\sum }}{N}_{i,3}\left(t\right)V_i}$

是三次均匀B样条函数。其中 $N_{0,3}\left(t\right),N_{1,3}\left(t\right),N_{2,3}\left(t\right),N_{3,3}\left(t\right)$ 为三次均匀B样条基， $V_0,V_1,V_2,V_3$ 为特征多边形的顶点。

令 ${P^{\prime }}_0=P_0,{P^{\prime }}_1=P_1,{P^{\prime }}_2=P_{1\frac{1}{2}},{P^{\prime }}_3=P_2,{P^{\prime }}_4=P_3,{V^{\prime }}_0=V_0,{V^{\prime }}_1=V_1,{V^{\prime }}_2=V_{1\frac{1}{2}},{V^{\prime }}_3=V_2,{V^{\prime }}_4=V_3$。下面构造插值于 ${P^{\prime }}_0,{P^{\prime }}_2,{P^{\prime }}_4$，并在 ${P^{\prime }}_0$ 处与 ${P^{\prime }}_0{P^{\prime }}_1$ 相切，在 ${P^{\prime }}_4$ 处与 ${P^{\prime }}_3{P^{\prime }}_4$ 相切的两段连续三次B样条曲线，需满足：

$\{\begin{array}{l}\frac{1}{6}{V^{\prime }}_0+\frac{2}{3}{V^{\prime }}_1+\frac{1}{6}{V^{\prime }}_2={P^{\prime }}_0,\\ \frac{1}{6}{V^{\prime }}_1+\frac{2}{3}{V^{\prime }}_2+\frac{1}{6}{V^{\prime }}_3={P^{\prime }}_2,\\ \frac{1}{6}{V^{\prime }}_2+\frac{2}{3}{V^{\prime }}_3+\frac{1}{6}{V^{\prime }}_4={P^{\prime }}_4,\\ -\frac{1}{2}{V^{\prime }}_0+\frac{1}{2}{V^{\prime }}_2={P^{\prime }}_1-{P^{\prime }}_0,\\ -\frac{1}{2}{V^{\prime }}_2+\frac{1}{2}{V^{\prime }}_4={P^{\prime }}_4-{P^{\prime }}_3.\end{array}$

解以 ${V^{\prime }}_0,{V^{\prime }}_1,{V^{\prime }}_2,{V^{\prime }}_3,{V^{\prime }}_4$ 为变量的线性方程组，得

$\{\begin{array}{l}{V^{\prime }}_0=\frac{5}{3}{P^{\prime }}_0-\frac{13}{6}{P^{\prime }}_1+2{P^{\prime }}_2-\frac{1}{6}{P^{\prime }}_3-\frac{1}{3}{P^{\prime }}_4,\\ {V^{\prime }}_1=\frac{7}{6}{P^{\prime }}_0+\frac{7}{12}{P^{\prime }}_1-{P^{\prime }}_2+\frac{1}{12}{P^{\prime }}_3+\frac{1}{6}{P^{\prime }}_4,\\ {V^{\prime }}_2=-\frac{1}{3}{P^{\prime }}_0-\frac{1}{6}{P^{\prime }}_1+2{P^{\prime }}_2-\frac{1}{6}{P^{\prime }}_3-\frac{1}{3}{P^{\prime }}_4,\\ {V^{\prime }}_3=\frac{1}{6}{P^{\prime }}_0+\frac{1}{12}{P^{\prime }}_1-{P^{\prime }}_2+\frac{7}{12}{P^{\prime }}_3+\frac{7}{6}{P^{\prime }}_4,\\ {V^{\prime }}_4=-\frac{1}{3}{P^{\prime }}_0-\frac{1}{6}{P^{\prime }}_1+2{P^{\prime }}_2-\frac{13}{6}{P^{\prime }}_3+\frac{5}{3}{P^{\prime }}_4.\end{array}$

于是，所求两段连续的三次B样条曲线的表达式为

$r_1\left(t\right)={\displaystyle \underset{i=0}{\overset{3}{\sum }}{B}_{i3}\left(t\right){V^{\prime }}_i},r_2\left(t\right)={\displaystyle \underset{i=0}{\overset{3}{\sum }}{B}_{i3}\left(t\right){V^{\prime }}_{i+1}}.$

例1设 $P_0=\left(5,-3,0\right),P_1=\left(0,-3,0\right),P_2=\left(0,5,0\right),P_3=\left(0,5,3\right)$，适当选取 ${V^{\prime }}_2$，调整三次B样条曲线的形状。

1) 取 ${P^{\prime }}_2=\left(\frac{5}{2},\frac{5}{2},0\right)$，得以 ${P^{\prime }}_0,{P^{\prime }}_1,{P^{\prime }}_2,{P^{\prime }}_3,{P^{\prime }}_4$ 为控制顶点，并过 ${P^{\prime }}_0,{P^{\prime }}_2,{P^{\prime }}_4$，且在 $P_0$ 处与 $P_0{P}_1$ 相切，在 $P_3$ 处与 $P_2{P}_3$ 相切的两段连续三次B样条曲线，其参数表达式

$r_1\left(t\right)=\{\begin{array}{l}{x}_1=\frac{5}{4}{t}^3-\frac{15}{4}{t}^2+5,\\ y_1=-\frac{11}{4}{t}^3+\frac{21}{4}{t}^2+3t-3,\\ z_1=\frac{3}{2}{t}^3-\frac{2}{3}{t}^2.\end{array}{r}_2\left(t\right)=\{\begin{array}{l}{x}_2=\frac{5}{4}{t}^3-\frac{15}{4}t+\frac{5}{2},\\ y_2=\frac{1}{4}{t}^3-3t^2+\frac{21}{4}t+\frac{5}{2},\\ z_2=-\frac{3}{2}{t}^3+3t^2+\frac{3}{2}t.\end{array}t\in \left[0,1\right]$

如图1。

2) 取 ${P^{\prime }}_2=\left(\frac{5}{2},\frac{5}{2},\frac{\text{1}}{\text{5}}\right)$，得以 ${P^{\prime }}_0,{P^{\prime }}_1,{P^{\prime }}_2,{P^{\prime }}_3,{P^{\prime }}_4$ 为控制顶点，过点 ${P^{\prime }}_0,{P^{\prime }}_2,{P^{\prime }}_4$，并在 ${P^{\prime }}_0$ 处与 ${P^{\prime }}_0{P^{\prime }}_1$ 相切，在 ${P^{\prime }}_4$ 处与 ${P^{\prime }}_3{P^{\prime }}_4$ 相切的两段连续三次B样条曲线，其参数表达式为

$r_1\left(t\right)=\{\begin{array}{l}{x}_1=\frac{5}{4}{t}^3-\frac{15}{4}{t}^2+5,\\ y_1=-\frac{11}{4}{t}^3+\frac{21}{4}{t}^2+3t-3,\\ z_1=\frac{\text{11}}{\text{10}}{t}^3-\frac{\text{9}}{\text{10}}{t}^2.\end{array}{r}_2\left(t\right)=\{\begin{array}{l}{x}_2=\frac{5}{4}{t}^3-\frac{15}{4}t+\frac{5}{2},\\ y_2=\frac{1}{4}{t}^3-3t^2+\frac{21}{4}t+\frac{5}{2},\\ z_2=-\frac{11}{10}{t}^3+\frac{12}{5}{t}^2+\frac{3}{2}t+\frac{1}{5}.\end{array}t\in \left[0,1\right]$

如图1。

3) 取 ${P^{\prime }}_2=\left(\frac{5}{2},\frac{5}{2},\frac{\text{2}}{\text{5}}\right)$，得以 ${P^{\prime }}_0,{P^{\prime }}_1,{P^{\prime }}_2,{P^{\prime }}_3,{P^{\prime }}_4$ 为控制顶点，过点 ${P^{\prime }}_0,{P^{\prime }}_2,{P^{\prime }}_4$，并在 ${P^{\prime }}_0$ 处与 ${P^{\prime }}_0{P^{\prime }}_1$ 相切，在 ${P^{\prime }}_4$ 处与 ${P^{\prime }}_3{P^{\prime }}_4$ 相切的两段连续三次B样条曲线，其参数表达式为

$r_1\left(t\right)=\{\begin{array}{l}{x}_1=\frac{5}{4}{t}^3-\frac{15}{4}{t}^2+5,\\ y_1=-\frac{11}{4}{t}^3+\frac{21}{4}{t}^2+3t-3,\\ z_1=\frac{\text{7}}{\text{10}}{t}^3-\frac{\text{3}}{\text{10}}{t}^2.\end{array}{r}_2\left(t\right)=\{\begin{array}{l}{x}_2=\frac{5}{4}{t}^3-\frac{15}{4}t+\frac{5}{2},\\ y_2=\frac{1}{4}{t}^3-3t^2+\frac{21}{4}t+\frac{5}{2},\\ z_2=-\frac{\text{7}}{10}{t}^3+\frac{\text{9}}{5}{t}^2+\frac{3}{2}t+\frac{\text{2}}{5}.\end{array}t\in \left[0,1\right]$

如图1。

将 ${V^{\prime }}_2$ 上述三种取值时得到的三条两段连续的三次B样条曲线放在一起，如图2。

${V^{\prime }}_2$ 可以取多种值，利用B样条曲线的局部修改性质可得满足应用条件的轴线应用于轴线异面管道拼接中。

图1. 为 ${V^{\prime }}_2$ 取三种不同值时两段连续的三次B样条曲线与轴线连接情况

图2. 为 ${V^{\prime }}_2$ 取三种不同值时三条两段连续的三次B样条曲线放在一起的情况

3. 过指定顶点的三次B样条曲线及其在轴线异面椭圆管道拼接中的应用

我们利用三次B样条曲线的局部修改性质，构造了过指定顶点的两段连续的三次B样条曲线，曲线通过控制多边形的两端点并与始末两边相切，而曲线的形状可由第三个指定顶点的位置来调控。这样，我们可将其应用于轴线光滑拼接的轴线异面管道拼接中 [4] [5]。

例2设 ${\Phi }_1\left(t,\phi \right)$ 和 ${\Phi }_2\left(t,\phi \right)$ 分别是长半轴和短半轴分别相同的两个椭圆管道，其参数表达式为

${\Phi }_{\text{1}}\left(t,\phi \right):\{\begin{array}{l}x\left(t,\phi \right)=5+0.5\cdot \cos \phi ,\\ y\left(t,\phi \right)=-3+t,\\ z\left(t,\phi \right)=0+0.8\cdot \sin \phi .\end{array}t\in \left[0,1\right],\phi \in \left[0,2\pi \right]$

和

${\Phi }_2\left(t,\phi \right):\{\begin{array}{l}x\left(t,\phi \right)=0+0.8\cdot \cos \phi ,\\ y\left(t,\phi \right)=5+0.5\cdot \sin \phi ,\\ z\left(t,\phi \right)=3+t.\end{array}t\in \left[0,1\right],\phi \in \left[0,2\pi \right]$

其中 ${\gamma }_1\left(t\right)=\{\begin{array}{l}{x}_1(t)=5,\\ y_1(t)=-3+t,\\ z_1(t)=0.\end{array}{\gamma }_2\left(t\right)=\{\begin{array}{l}{x}_2(t)=0,\\ y_2(t)=5,\\ z_2(t)=3+t.\end{array}t\in \left[0,1\right]$ 是轴线的参数表达式，取三次B样条曲线的控制多边形的顶点为 ${P^{\prime }}_0=\left(5,-3,0\right),{P^{\prime }}_1=\left(0,-3,0\right),{P^{\prime }}_2=\left(\frac{5}{2},\frac{5}{2},\frac{\text{2}}{\text{5}}\right),{P^{\prime }}_3=\left(0,5,0\right),{P^{\prime }}_4=\left(0,5,3\right)$，则两段连续的三次B样

条曲线的表达式为例1(3)。构造基于两段连续的轴线光滑拼接的两段椭圆管道，其参数表达式为

${\Psi }_1\left(t,\phi \right)=\{\begin{array}{l}{x}_1\left(t,\phi \right)=\frac{5}{4}{t}^3-\frac{15}{4}{t}^2+5+a{B}_{11}\left(t\right)\cos \phi +b{N}_{11}\left(t\right)\sin \phi ,\\ y_1\left(t,\phi \right)=-\frac{11}{4}{t}^3+\frac{21}{4}{t}^2+3t-3+a{B}_{12}\left(t\right)\cos \phi +b{N}_{12}\left(t\right)\sin \phi ,\\ z_1\left(t,\phi \right)=\frac{\text{7}}{\text{10}}{t}^3-\frac{\text{3}}{\text{10}}{t}^2+a{B}_{13}\left(t\right)\cos \phi +b{N}_{13}\left(t\right)\sin \phi .\end{array}t\in \left[0,1\right],\phi \in \left[0,2\pi \right]$

和

${\Psi }_2\left(t,\phi \right)=\{\begin{array}{l}{x}_2=\frac{5}{4}{t}^3-\frac{15}{4}t+\frac{5}{2}a{B}_{21}\left(t\right)\cos \phi +b{N}_{21}\left(t\right)\sin \phi ,\\ y_2=\frac{1}{4}{t}^3-3t^2+\frac{21}{4}t+\frac{5}{2}a{B}_{22}\left(t\right)\cos \phi +b{N}_{22}\left(t\right)\sin \phi ,\\ z_2=-\frac{\text{7}}{10}{t}^3+\frac{\text{9}}{5}{t}^2+\frac{3}{2}t+\frac{\text{2}}{5}a{B}_{23}\left(t\right)\cos \phi +b{N}_{23}\left(t\right)\sin \phi .\end{array}t\in \left[0,1\right],\phi \in \left[0,2\pi \right]$

其中 $N_i\left(t\right)=\left(N_{i1}\left(t\right),N_{i2}\left(t\right),N_{i3}\left(t\right)\right)$ 和 $B_i\left(t\right)=\left(B_{i1}\left(t\right),B_{i2}\left(t\right),B_{13}\left(t\right)\right),i=1,2$ 分别是连续两段B样条曲线的主法矢和副法矢与 ${\Phi }_1\left(t,\phi \right)$ 和 ${\Phi }_2\left(t,\phi \right)$ 的拼接效果如图3。

图3. $a=0.5$,$b=0.8$ 时，构造的两段光滑拼接的椭圆管道与给定轴线异面椭圆管道拼接效果图

椭圆管道拼接不同于圆管道。在构造椭圆管道时，轴线上每一点处与轴线垂直的平面上的自然坐标系 $\left[r_i\left(t\right),N_i\left(t\right),B_i\left(t\right)\right]$ 上构造椭圆，椭圆的长半轴和短半轴的方向与主法矢和副法矢的方向相关。例2中, ${\Phi }_1\left(t,\phi \right)$ 与 ${\Psi }_1\left(t,\phi \right)$ 的长半轴与短半轴是相对应的，而 ${\Phi }_2\left(t,\phi \right)$ 与 ${\Psi }_2\left(t,\phi \right)$ 的长半轴与短半轴是不相对应的。在 $r_2\left(t\right)$ 每一点 $t\in \left[0,1\right]$ 处对 $\left[r_2\left(t\right),N_2\left(t\right),B_2\left(t\right)\right]$ 坐标轴旋转变换

$\{\begin{array}{l}{e}_1\left(t\right)=-B_2\left(t\right)\sin \theta +N_2\left(t\right)\cos \theta ,\\ e_2\left(t\right)=B_2\left(t\right)\cos \theta +N_2\left(t\right)\sin \theta .\end{array}t\in \left[0,1\right]$

其中 $\theta ={\theta }_1+\left({\theta }_2-{\theta }_1\right)t$，${\theta }_1$ 为 $N_2\left(0\right)$ (或 $B_2\left(0\right)$ )与 ${\Phi }_2\left(0\right)$ 长半轴(或短半轴)夹角， ${\theta }_2$ 为 $N_2\left(1\right)$ (或 $B_2\left(1\right)$ )与 ${\Phi }_2\left(1\right)$ 长半轴(或短半轴)夹角。调整后的拼接效果如图4。

图4. $a=0.5$,$b=0.8$ 时，调整后两段光滑拼接的椭圆管道与给定轴线异面椭圆管道光滑拼接效果图

当 $a_1\ne a_2,b_1\ne b_2$ 时，还可以构造两段光滑拼接的椭圆管道光滑拼接轴线异面的椭圆管道 [5]。

例3设 ${\Phi }_1\left(t,\phi \right)$ 和 ${\Phi }_2\left(t,\phi \right)$ 分别是长半轴和短半轴分别不同的两个椭圆管道，其参数表达式为

${\Phi }_{\text{1}}\left(t,\phi \right):\{\begin{array}{l}x\left(t,\phi \right)=5+0.\text{48}\cdot \cos \phi ,\\ y\left(t,\phi \right)=-3+t,\\ z\left(t,\phi \right)=0+0.8\cdot \sin \phi .\end{array}t\in \left[0,1\right],\phi \in \left[0,2\pi \right]$

和

${\Phi }_2\left(t,\phi \right):\{\begin{array}{l}x\left(t,\phi \right)=0+\text{1}\text{.2}\cdot \cos \phi ,\\ y\left(t,\phi \right)=5+0.\text{72}\cdot \sin \phi ,\\ z\left(t,\phi \right)=3+t.\end{array}t\in \left[0,1\right],\phi \in \left[0,2\pi \right]$

光滑拼接两个异面轴线的两段连续三次B样条曲线如例1 (3)。构造光滑拼接两个轴线异面椭圆管道某一母线的两段连续三次B样条曲线为

${r^{\prime }}_1\left(t\right)=\{\begin{array}{l}{x^{\prime }}_1=1.05t^3-3.15t^2+4.2,\\ {y^{\prime }}_1=-2.45t^3+4.35t^2+3t-3\cdot 1,\\ {z^{\prime }}_1=\frac{7}{10}{t}^3-\frac{3}{10}{t}^2.\end{array}{r^{\prime }}_2\left(t\right)=\{\begin{array}{l}{x^{\prime }}_2=\text{1}\text{.05}{t}^3+2.1t^2-3.15t,\\ {y^{\prime }}_2=0.55t^3+3t^2+4.35t+1.\text{9},\\ {z^{\prime }}_2=-\frac{7}{10}{t}^3+\frac{9}{5}{t}^2+\frac{3}{2}t+\frac{2}{5}.\end{array}t\in \left[0,1\right]$

光滑拼接两个轴线异面椭圆管道的两段椭圆管道的参数表达式为

${\Psi }_i\left(t,\phi \right):\{\begin{array}{l}{x}_i\left(t,\phi \right)=x_i\left(t\right)+a_i\alpha \left(t\right)e_{i1}\left(t\right)\cos \phi +b_i\alpha \left(t\right){e^{\prime }}_{i1}\left(t\right)\sin \phi ,\\ y_i\left(t,\phi \right)=y_i\left(t\right)+a_i\alpha \left(t\right)e_{i2}\left(t\right)\cos \phi +b_i\alpha \left(t\right){e^{\prime }}_{i2}\left(t\right)\sin \phi ,\\ z_i\left(t,\phi \right)=z_i\left(t\right)+a_i\alpha \left(t\right)e_{i3}\left(t\right)\cos \phi +b_i\alpha \left(s\right){e^{\prime }}_{i3}\left(t\right)\sin \phi .\end{array}i=1,2,t\in \left[0,1\right],\phi \in \left[0,2\pi \right]$

其中 ${\alpha }_i\left(t\right)=\sqrt{{\left({x^{\prime }}_i\left(t\right)-x_i\left(t\right)\right)}^2+{\left({y^{\prime }}_i\left(t\right)-y_i\left(t\right)\right)}^2+{\left({z^{\prime }}_i\left(t\right)-z_i\left(t\right)\right)}^2},e_i=\left(e_{i1},e_{i2},e_{i3}\right)$ 和 ${e^{\prime }}_i=\left({e^{\prime }}_{i1},{e^{\prime }}_{i2},{e^{\prime }}_{i3}\right),i=1,2,t\in \left[0,1\right]$，为轴线上互为垂直的单位向量， $a_1=0.48$，$b_1=0.8$，$a_2=0.72$，$b_1=1.2$。其拼接效果如图5。

图5. 两个长半轴和短半轴互不相同的轴线异面椭圆管道光滑拼接效果图

图5中， ${\Phi }_1\left(t,\phi \right),{\Phi }_2\left(t,\phi \right)$ 分别为

${\Phi }_{\text{1}}\left(t,\phi \right):\{\begin{array}{l}x\left(t,\phi \right)=5+0.\text{48}\cdot \cos \phi ,\\ y\left(t,\phi \right)=-3+t,\\ z\left(t,\phi \right)=0+0.8\cdot \cos \phi .\end{array}t\in \left[0,1\right],\phi \in \left[0,2\pi \right]$

和

${\Phi }_2\left(t,\phi \right)=\{\begin{array}{l}{x}_2=\text{0}+a{B}_{21}\left(1\right)\cos \phi +b{N}_{21}\left(1\right)\sin \phi ,\\ y_2=\text{5}+a{B}_{22}\left(1\right)\cos \phi +b{N}_{22}\left(1\right)\sin \phi ,\\ z_2=\text{3+}t+a{B}_{23}\left(1\right)\cos \phi +b{N}_{23}\left(1\right)\sin \phi .\end{array}t\in \left[0,1\right],\phi \in \left[0,2\pi \right]$

4. 结束语

用指定顶点的两段连续三次B样条曲线为轴线的两段光滑拼接的椭圆管道拼接两个轴线异面的椭圆管道，与三段连续的三次B样条曲线为轴线的三段光滑拼接的管道拼接两个轴线异面管道相比较段数少，与插值于两端点的三次有理B样条曲线相比较能够构造椭圆管道拼接，而与二次B样条曲线为轴线的管道光滑拼接轴线异面管道相比较能够调整轴线的位置。具有理论意义和较好的应用价值。

基金项目

国家自然科学基金项目资助(11561052)，浙江省教育厅一般项目资助Y (201636628)。

文章引用

白根柱. 过指定顶点的三次B样条曲线及其应用
Cubic B-Spline Curve of Passing Appointed Vertexes and Its Application[J]. 应用数学进展, 2020, 09(10): 1821-1828. https://doi.org/10.12677/AAM.2020.910211

参考文献