[发明专利]一种基于延伸的产品外形空间曲线的拼接方法

摘要

一种基于延伸的产品外形空间曲线拼接的CAD方法，属于曲线曲面的CAD领域，其特征在于，初始化阶段输入要进行拼接操作的两条空间曲线C₁和D；对其中的曲线C₁相继实行三次延伸操作使得与曲线D实现相连通；调整延伸曲线的两个控制顶点，实现两条曲线的二阶几何连续。基于该方法在不添加第三条曲线的情况下，既填补了两条曲线间原有的缝隙，又不改变曲线的原有部分，实现了一种新的曲线拼接效果。

主权项

1.一种基于延伸的产品外形空间曲线拼接的CAD方法，其特征在于，所述方法是在计算机上依次按照如下步骤实现的：步骤(1)，初始化：输入要实行拼接操作的空间曲线C₁(t)和空间曲线D(s)，每条曲线采用B样条的表示方法，对于曲线C₁其形式为：C1(t)=Σi=0I1-1P1,iNi,c1(t)]]>其中t是空间曲线C₁(t)的参数坐标，I₁是曲线C₁(t)所含的控制顶点个数，P_1，i是曲线C₁(t)中序号为i的控制顶点的函数值，i＝0，1，…，I₁-1，采用笛卡尔坐标表示，是曲线C₁(t)中序号为i的B样条基函数，该B样条基函数定义在如下节点向量序列之上：其中k₁是曲线C₁的幂，为设定值，等为节点向量序列中的节点，为设定值，上述的B样条基函数采用Matlab系统样条工具库中B样条基函数的定义方式，对于曲线D其形式为：D(s)=Σj=0J-1QjNj,D(s)]]>其中s是空间曲线D(s)的参数坐标，J是曲线D(s)所含的控制顶点个数，Q_j是曲线D(s)中序号为j的控制顶点的函数值，j＝0，1，…，J-1，采用笛卡尔坐标表示，N_j，D是曲线D(s)中序号为j的B样条基函数，该B样条基函数定义在如下节点向量序列之上：其中k_D是曲线D的幂，为设定值，…，W_J-1等为节点向量序列中的节点，为设定值，上述的B样条基函数采用Matlab系统样条工具库中B样条基函数的定义方式；步骤(2)，对曲线C₁(t)进行向曲线D(s)方向的延伸，其步骤如下：步骤(2.1)在所述的曲线C₁(t)和曲线D(s)之间按以下方式生成R₁和R₂两个空间点，R1=0.67P1,I1-1+0.33Q0,]]>R2=0.33P1,I1-1+0.67Q0,]]>其中是曲线C₁(t)上最末一个控制顶点的函数值，Q₀是曲线D(s)上最初一个控制顶点的函数值，步骤(2.2)把曲线C₁延伸到所述的点R₁，并且设延伸后的曲线为C₂(t)，其表示形式为C2(t)=Σi=0I2-1P2,iNi,c2(t)]]>其中I₂是包含曲线C₁(t)在内的曲线C₂(t)中所含的控制顶点个数，且有I₂＝I₁+1，P_2，i是曲线C₂(t)中序号为i的控制顶点的函数值，i＝0，1，…，I₁，是曲线C₂(t)中序号为i的B样条基函数，所述的定义在如下节点向量序列之上：其中a=1+||P1,I1-1-R1||Σi=0I1-k1-1||C1(uk1+i+1)-C1(uk1+i)||,]]>||·||表示空间的欧氏距离，步骤(2.3)按下述步骤计算曲线C₂(t)中序号为i的控制顶点的值P_2，i步骤(2.3.1)设置初值P~i0=P1,i,]]>i＝I₁-k₁，I₁-k₁+1，…，I₁-1步骤(2.3.2)按下式递推计算空间点P~ir=P~ir-1i=I1-k1,···,I1-r-1P~ir=P~ir-1-(1-λ)P~i-1rλi=I1-r,I1-r+1···,I1-1r=1,2,···,k1-1]]>其中λ=uI1-uiui+r+1-ui]]>步骤(2.3.3)得出P_2，i的设置结果P2,i=P1,ii=0,1,···,I1-k1-1P~ik1-1i=I1-k1,I1-k1+1,···,I1-1R1i=I1,]]>步骤(2.4)按照步骤(2.2)所述的方式将所述曲线C₂(t)延伸到点R₂，并且设延伸后的包含曲线C₂(t)的曲线为C₃(t)，步骤(2.5)再按照步骤(2.2)所述的方式将所述曲线C₃(t)延伸到所述曲线D(s)的第0个控制顶点Q₀，并且设延伸后的包含曲线C₃(t)的曲线为C₄(t)：C4(t)=Σi=0I4-1P4,iNi,c4(t)]]>其中I₄是曲线C₄(t)所含的控制顶点个数且有I₄＝I₁+3，是曲线C₄(t)中序号为i的B样条基函数，i＝0，1，…，I₁+3，定义在如下归一化形式的节点向量序列之上：步骤(3)，修改曲线C₄(t)的两个控制顶点和，使得修改后的曲线C₄(t)与曲线D(s)在其交接处即t＝1和s＝0处实现2阶几何连续，和和计算方法如下：P4,I1+1=1NI1+1,C4′(1)[α(Q0N0,D′(0)+Q1N1,D′(0))-P4,I1+2NI1+2,C4′(1)]]]>P4,I1=1NI1,C4′′(1)[α2(Q0N0,D′′(0)+Q1N1,D′′(0)+Q2N2,D′′(0))-P4,I1+2NI1+2,C4′′(1)-P4,I1+1NI1+1,C4′′(1)]]]>其中α为用户指定的值，取值范围为(0，1]之间的正值，N′和N″为相应的B样条基函数的1阶导数以及2阶导数。