车身曲线曲面数学模型基础.pptx

《车身曲线曲面数学模型基础.pptx》由会员分享，可在线阅读，更多相关《车身曲线曲面数学模型基础.pptx（78页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、4.1 4.1 参数样条曲线及孔斯曲面参数样条曲线及孔斯曲面4.1.1 4.1.1 三次样条曲线三次样条曲线(cubic spline curve)(cubic spline curve)数数学学上上的的样样条条函函数数是是对对绘绘图图用用的的样样条条的的模模拟拟。如将样条简化为弹性细杆如将样条简化为弹性细杆,必定满足欧拉方程：必定满足欧拉方程：M(x)=EIK(x)(4.1.1-1)(4.1.1-1)其其中中M（x）是是弯弯矩矩，E E是是杨杨氏氏系系数数，I I是是截截面面惯惯性性矩矩，K(x)是是样样条条的的曲曲率率。从从(4.1.1-1)(4.1.1-1)式式出出发发，经经数数学学推推

2、导导可可得出如下的得出如下的三次样条函数表达式：三次样条函数表达式：（4.1.1-15)4.1.1-15)三次样条函数三次样条函数S(x)的本质是是：一致通过型值点、二阶连续可导一致通过型值点、二阶连续可导的分段三次多项式函数。的分段三次多项式函数。第1页/共78页M（二阶导数）关系式（二阶导数）关系式在各中间（连接）点一阶导数连续，SS(xi-0)=)=S S(xi+0)，即(式中:hi=xi-xi-1)各项乘以 ，得：令:则有则有Mi-1+2Mi+i Mi+1=di (i=1,2,n-1)(4.1.1-16)当当i取取值值1,2,n-1时时，可可得得到到n-1个个形形如如(4.1.1-16

3、)的的M关关系系式式。但未知数二阶导数Mi却有n n+1+1个,即M0,M1,Mn 。要唯一定解,必须再附加两个方程。通常按实际问题的具体情况,在样条两端,即P P0 0和P Pn n处给出约束条件，常用的边界条件有：第2页/共78页常用的边界条件有：1给定两端的斜率m0=y0 和mn=yn以x=x0,i=1代入式(4.1.1-14)，得(4.1.1-17)以x=xn,i=n代入式(4.1.1-14)，得(4.1.1-18)式(4.1.1-16)和这两个附加方程合在一起得到有确定解的线性方程组。写成矩阵形式为(4.1.1-19)。2.给定两端的二阶导数M0=y0,Mn=yn这可以写成：2M0+

4、0M1=2y00Mn-1+2Mn=2 yn此 时 式(4.1.1-19)中 的0=0,d0=2y0，n=0,dn=2 yn。如果y0=yn=0,则称为自然插值三次样条函数。3.如果取0=-2,d0=0，n=-2,dn=0则M0=M1,Mn-1=Mn，这就是抛物端边界条件。第3页/共78页三次样条函数的解法（4.1.1-16)由端点条件补充两个方程后，得出如下线性方程组：（4.1.1-19)式中:可以用“追赶法”(参看附录A)求解(4.1.1-19)式,解出Mi(i=1,2,n)代入(4.1.1-15),便可由(4.1.1-15)式计算出样条曲线上的一系列插值点。第4页/共78页4.1.2三次参

5、数样条曲线在大挠度情况下,三次样条函数的光顺性可能变坏。用三次样条函数表示的插值曲线,依赖于座标系的选择,不具有几何不变性。有时旋转座标轴也不可能满足小挠度条件在这些情况下,最常用的处理办法之一是将曲线参数化,即将曲线上点的座标分别用某种参数表示：(4.1.2-1)其中t为参数，常取曲线内在的量弧长作为参数，它与坐标系无关。若将t取作弧长s，则x和y作为分量，dxds和dyds都不会大于1，在(x,s),(y,s)平面上各构造一个三次样条函数：(4.1.2-2)曲线上的点比较密时，弦长之和近似于弧长，因此可取累加弦长作为三次参数样条曲线的参数。第5页/共78页设给定个点Pi（xi,yi）,i=

6、0,1,n,两相邻点之间的弦长为：(4.1.2-3)记:这里ti的几何意义是累加弦长，它近似等于弧长参数。在每一个节点pi 都有一个确定的ti与之对应。当然pi 的每一个坐标xi或yi也与ti一一对应。这就相当于给定了两组点（xi,ti）和（yi,ti）,i=0n。对于每一组点，都可按4.1.1节所述方法构造一个三次样条函数。这种曲线称为累加弦长三次参数样条曲线。在平面曲线的情况下，构造三次参数样条曲线相当于构造两遍三次样条曲线。这在工程上是经常使用的方法。第6页/共78页4.1.3弗格森曲弗格森曲线线下面讨论参数样条曲线中的某一段，并用端点及端点的导数来表达出这段曲线的方程。设参数为u，第i

7、段曲线对应的参数范围为，在0,1区间上对应于两个端点型值点的函数值及一阶导数值分别为r(0),r(1),r(0),r(1)。则插值函数为(4.1.3-1)那么(4.1.3-2)将四个已知条件代入以上两式，可解得四个系数a0,a1,a2,a3,再将求得的系数代回上式则得曲线段的方程为：(4.1.3-3)式中(4.1.3-4)我 们 称F0(u),F1(u),G0(u),G1(u)为 埃 尔 米 特(Hermite)基函数。第7页/共78页由式(4.1.3-4)可见，F0与F1专门控制端点的函数值对曲线形态的影响，G0和G1专门控制端点的一阶导数对曲线形态的影响。或者说，F0和G0控制左端的影响，

8、F1和G1控制右端的影响。由式(4.1.3-3)确定的曲线可以进一步整理为矩阵形式：（4.1.3-5）该曲线也叫弗格森曲线。第8页/共78页4.1.4 4.1.4 孔斯孔斯(Coons)(Coons)曲曲面面 是S.A.Coons(S.A.Coons(是美国波音公司搞实际设计的专家)提出的一种适用于CAGDCAGD的构作自由型曲面的方法。孔斯曲面的基本思想是：把所要描述的曲面看作由若干曲面片光滑拼接而成,每个曲面片一般用四条边界曲线来定义。且尽量用简缩符号来表达。双三次孔斯曲面的表达形式为：uw=UMCMTWT (4.1.4-1)式中 U=1 u u2 u3 MTWT分别是M和W的转置。第9页

9、/共78页 =（4.1.4.2)C称为角点信息矩阵。Coons曲面主要用于曲面设计，但用于曲面拟合也有较好的效果。（如，用零扭矢Coons曲面拟合飞机局部复杂曲面）。Coons曲面方法的主要缺点：1）角点扭矢与曲面内部形状的联系，难以掌握；2）确定扭矢是个难题；3）双三次Coons曲面的拼接只能达到一阶连续；4）角点信息多，占内存大，计算量也大。这些缺点限制了Coons曲面的广泛应用。第10页/共78页4.2贝齐尔(Bezier)方法 法法国国雷雷诺诺汽汽车车公公司司的的贝贝齐齐尔尔(P(PBezier)Bezier)于于19621962年年着着手手研研究究一一种种以以逼逼近近为为基基础础的的

10、构构造造曲曲线线与与曲曲面面的的方方法法,并并以以这这种种方方法法建建立立了了一一种种自自由由型型曲曲线线与与曲曲面面设设计计系系统统UNISURFUNISURF系系统统。该该系系统统于于19721972年年投投入入使使用用,至至今今已已有有很很大大的的发发展展。应应用用该该方方法法的的还还有有英英国国剑剑桥桥大大学学的的MultiobjectMultiobject实实用用设设计计系系统统等等。贝贝齐齐尔尔方方法法有有许许多多优优良良的的性性质质，已已成成为为 自自由由型型曲曲线线曲曲面面造造型型先先进进的的数数学学方方法之一。法之一。4.2.1贝齐尔曲线贝齐尔曲线 设设 有有n+1个个 控控

11、 制制 顶顶 点点P0,P1,Pn，借借 助助 于于 一一 组组Bernstein基基函函数数Bi,n(u)=C(n,i)ui(1-u)n-i(i=0,1,2,n),可以定义一条曲线可以定义一条曲线:(4.2.1-1)该曲线称为该曲线称为n次贝齐尔曲线次贝齐尔曲线。其中其中：C(n，i)是组合数，是组合数，n是贝齐尔曲线次数，是贝齐尔曲线次数，i是顶点标号是顶点标号,u为曲线的参数为曲线的参数。第11页/共78页贝齐尔曲线的分量的形式：P(u)=x(u),y(u),z(u)(4.2.1-2)Pi=(xi,yi,zi,)(4.2.1-3)(4.2.1-4)为了了解控制顶点Pi(i=0,1,2,n

12、)如何通过基函数Bi,n(u)控制曲线的形状，先考察n3的情形。由定义可得：(4.2.1-6)Bi,3(u)(i=0,1,2,3)的图形见图4.2.1-1。第12页/共78页控制顶点Pi通过基函数Bi,n(u)控制曲线的形状 B0,3(u)在u=0附近影响最大，B1,3(u)在u=13附近，B2,3(u)在u=23附近影响最大，而B3,3(u)在附近影响最大。因为u=0时P(0)=P0,故P0决定曲线的起点。U=1时P(1)=P3，故P3决定曲线的终点。P1则主要通过B1,3(u)影响曲线在u=13附近的形状；P2则主要通过B2,3(u)影响曲线在u=23附近的形状。图4.2.1-2是P0,P

13、1,P2,P3控制曲线形状的几个例子。第13页/共78页贝齐尔曲线的性质1 1端点性质。根据贝齐尔曲线的定义可以证明：P P(0)=(0)=P P0 0，P P(1)=(1)=P Pn n，表明P0和Pn是贝齐尔曲线的起点和终点。还可以证明：P P(0)=(0)=nana1 1 ,P P(1)=(1)=nanan n 式中ai=Pi-Pi-1(i=1,2,n)代表贝齐尔多边形的边向量，说明贝齐尔曲线分别以a1和an为起端和终端的切向，类似地也有：P”(0)=n(n-1)(a2-a1)P”(1)=n(n-1)(an-an-1)2 2对称性。保持贝齐尔曲线诸顶点的位置不变，只把次序完全颠倒过来，新

14、的顶点序列记为P*i=Pn-i(i=0,1,2,n)生成贝齐尔曲线仍然是原曲线，第14页/共78页 3 3凸包性。贝齐尔曲线上各点一定落在特征多边 形在的凸包之中。4 4几何不变性。由于贝齐尔曲线曲线的形状与特征多边形各顶点P Pi i有关，它不依赖于座标系的选择。5 5交互能力。控制多边形P P0 0,P P1 1,P Pn n大致地勾画出BezierBezier曲线P P(u u)的形状，要改变P P(u u)的形状，只需改变P P0 0,P P1 1,P Pn n 的位置，把控制多边形作为曲线输入和人机交互的手段，既直观又简便。6 6保凸性。如果平面上的凸控制多边形能导致所生成的曲线为凸

15、曲线，则称这个生成曲线的方法具有保凸性。BezierBezier曲线具有这种保凸性质。7.7.变差缩减性。如果BezierBezier曲线P P(u u)的控制多边形是一平面图形,则该平面内的任意直线与P P(u u)的交点个数不多于该直线与控制多边形的交点的个数。这一性质被称为变差缩减性。此性质反映了BezierBezier曲线比控制多边形所在的折线更光顺。第15页/共78页8.曲线的可分割性。为求出曲线上任意一点，贝齐尔曾给出一种有趣的几何分割作图法。设Pk(k=0,1,2,n)是控制多边形顶点,在每边上取点 ，重复上面的过程 (4.2.1-14)可以证明，如图是入=1/2的分割情况。入可

16、取几个不同的值，得到曲线上几个不同的点，再根据贝齐尔曲线端点性质，便可绘出曲线。（画图讲解）第16页/共78页4.2.2三次贝齐尔曲线三次贝齐尔曲线 在产品外形设计中，C2阶连续的三次贝齐尔曲线已相当理想。高次贝齐尔曲线的许多问题还有待于理论上的解决。下面给出实际中常用的三次贝齐尔曲线的表达式。根据贝齐尔曲线的表达式并设n=3，则得三次贝齐尔曲线的表达式：（4.2.2-1）或写成矩阵形式：（4.2.2-2）这里的四个基函数分别为:第17页/共78页贝齐尔曲线的计算1.正算根据给定的特征多边形顶点Pi，构造贝齐尔曲线表达式（按照（4.2.2-1)式),并计算曲线上的一系列点。（写成分量形式表达式

17、进行计算）。2.反算用给定的曲线上的型值指点xi，yi及对应的参数ui反求出顶点Pi(运用贝齐尔曲线公式进行逆运算)，再构造贝齐尔曲线表达式并计算曲线上的一系列点。第18页/共78页4.2.3贝齐尔样条贝齐尔样条曲线曲线贝齐尔曲线是一整段n次参数曲线，不是样条，但可以把多段贝齐尔曲线光滑连接起来构造样条曲线，为了使合成后的整条曲线达到一定的连续性，连接处要满足特定的条件。设已经给定两条贝齐尔曲线L(n次)和L*(m次)，它 们 的 特 征 顶 点 分 别 为Pi(i=0,1,2,n)和Qi(i=0,1,2,m)，特 征 多 边 形 的 边 向 量 分 别 为ai(i=1,2,n)和bi(i=1

18、,2,m)，见图4.2.3-1。两条曲线L和L*达到C1连续的充耍条件是:(1)P(1)=Q(0)(4.2.3-1)；(2)Q/(0)=P/(1)0(4.2.3-2)由此可知两条贝齐尔曲线L与L*达到C1连续的充要条件是L的终点同L*的起点重合,且b1=an 0(4.2.3-3)对于两条空间曲线L和L*如果要达到C2连续要求，那么除了上述(4.2.3-2)条件以外,还应使连接点处有连续变化的曲率K和单位主法矢N。第19页/共78页因为在连接处已有相同的单位切矢T，又由于单位副法矢B=TN。因此在连接处还应满足(3)、(4)两条件：(3)副法矢同向；(4)曲率相等。从贝齐尔曲线的端点性质可推得，

19、L在终点的副法矢和L*在起点的副法矢分别是:(1)=n2(n-1)(an-1 an)(4.2.3-4)(0)=m2(m-1)(b1 b2)(4.2.3-5)由于要求副法线向量同向，则边向量an-1,an,b1,b2共面。又 知an-1,an是 两 个 线 性 无 关 的 量，再 考 虑 到 式(4.2.3-3)，则有(3)b2=-an-1+an (4.2.3-6)如果只考虑条件(3)，和可以是任意常数。若考虑条件(4)，则还有一定的限制。由曲线L和L*在连接点的曲率可推得:，(4)（4.2.3-7）综上所述可得出如下结论:1.L和L*拼接达到一阶连续,要满足条件（1）和（2），即L的末端点与L

20、*的首端点重合且斜率相等；2.L和L*拼接达到二阶连续,要满足条件（1）、（2）、（3）、(4)，即除了满足（1）、（2）外，在连接点处L和L*的副法矢同向且曲率相等。第20页/共78页4.2.4贝齐尔曲面贝齐尔曲面利用控制顶点和基函数生成曲线的方法很容易推广来生成曲面。现在，我们考虑(n+1)(m+1)个排成网格的控制顶点Pij(i=0,1,n;j=0,1,m)利用基函数Bi,n(u),Bi,m(w)就可以生成一块曲面:（4.2.4-1)该曲面称为nm次的贝齐尔曲面。显然,固定w,对u而言是一簇贝齐尔曲线;固定u,对w而言也是一簇贝齐尔曲线。可以认为,贝齐尔曲面是由贝齐尔曲线交织而成的曲面。

21、也就是说可以利用贝齐尔曲线的网格来绘制或显示贝齐尔曲面。第21页/共78页注意到 时,；时，；令 ，得，即Pi,0(i=0,1,n)恰好是P(u,0)的控制顶点。令 ，得 ，即Pi，m(i=0,1,n)恰好是P(u,1)的控制顶点。同理P0,j(j=0,1,m)恰好是P(0,w)的控制顶点，Pn,j(j=0,1,m)恰好是P(1,w)的控制顶点。而其余的Pi,j并不是P(u,w0)、P(u0,w)的控制顶点。此外，只有四个顶点P0,0,P0,m,Pn,0,Pn,m,与贝齐尔曲面的4个角点重合,并且在那里相切。第22页/共78页一个复杂的曲面往往不能用单一的贝齐尔曲面来实现。于是要用几块贝齐尔曲

22、面拼接起来,这时就要注意一定的连续性。对于3 3次贝齐尔曲面的情形，要4 4个控制点阵。根据定义用矩阵表示为 (4.2.4-2)其中 第23页/共78页曲面拼接条件：如果有二块33次贝齐尔曲面要进行拼接，如图4.2.4-2，当P(1)(1,w)=P(2)(0,w)对所有0w1成立时，拼接处连续。显然，这只要控制顶点满足P3,i(1)=P0,i(2)(i=0,1,2,3)(4.2.45)就行。而为了在拼接处满足C1连续，还要满足Pu(1)(1,w)Pw(1)(1,w)=Pu(2)(0,w)Pw(2)(0,w)(4.2.4-6)其中为常数，0w 1。这时除条件P3,i(1)=P0,i(2)(i=0

23、,1,2,3)外，最简的充分条件是:P3,i(1)-P2,i(1)=(P1,i(2)-P0,i(2)(i=0,1,2,3)(4.2.4-7)图4.2.4-2表示的两拼接曲面达到了C1级连续。图4.2.4-2第24页/共78页4.3 4.3 均匀均匀B B样条曲线样条曲线 4.3.1 B4.3.1 B样条方法引论 。以上讨论的贝齐尔方法，以上讨论的贝齐尔方法，因为它的许多优良性质，诸如直观性、凸包性、变差减小，因为它的许多优良性质，诸如直观性、凸包性、变差减小，可分割性等，已经成为从事机电产品几何外形设计的得力可分割性等，已经成为从事机电产品几何外形设计的得力工具。但使用中也发现，工具。但使用中

24、也发现，由于贝齐尔曲线曲面是采取单一由于贝齐尔曲线曲面是采取单一参数多项式的整体表示，它难以构造形状复杂的曲线曲面，参数多项式的整体表示，它难以构造形状复杂的曲线曲面，而不得不借助于拼接。拼接虽然灵活，但使用起来却不够而不得不借助于拼接。拼接虽然灵活，但使用起来却不够方便。方便。正由于贝齐尔曲线曲面的整体性，不可能作局部修正由于贝齐尔曲线曲面的整体性，不可能作局部修改，移动一个控制顶点的影响将波及整条曲线或整张曲面。改，移动一个控制顶点的影响将波及整条曲线或整张曲面。在在1972197219761976年年期期间间，里里森森费费尔尔德德(Riesnfeld)(Riesnfeld)、戈戈登登(G

25、ordon)(Gordon)、福福雷雷斯斯特特(Forrest)(Forrest)等等人人推推广广了了贝贝齐齐尔尔曲曲线线，改改用用B B样样条条基基代代替替贝贝齐齐尔尔曲曲线线的的伯伯恩恩斯斯坦坦基基，用用这这种种方方法法构构作作的的曲曲线线曲曲面面叫叫B B样样条条曲曲线线曲曲面面。B B样样条条曲曲线线与与曲曲面面也也具具备备良良好好的的性性质质。它它继继承承了了贝贝齐齐尔尔曲曲线线的的直直观观性性等等优优良良属属性性，又又克克服服了了贝贝齐齐尔尔方方法法的的不不足足之之处处。B B样样条条曲曲线线与与特特征多边形相当接近，便于局部修改。征多边形相当接近，便于局部修改。首先介绍工程上常用

26、的三次首先介绍工程上常用的三次B B样条曲线，然后再对样条曲线，然后再对B B样条曲线作进一步讨论。样条曲线作进一步讨论。第25页/共78页4.3.2 4.3.2 三次B B样条曲线段 B B样条曲线也是即分段又连续，先讨论各分段的特性，再解决各分段间的连续性问题。(1)(1)三次B B样条基。B B样条基函数可以由多种方法推导：如差商定义、德布尔一考克斯的递推定义、考虑曲线段之间连续性要求的几何定义等。由于推导的途径不一，B B样条基函数的表达式各有不同，但实质是完全一致的。现在直接引出工程上经常应用的三次B B样条基函数的矩阵表达式：(4.3.2-1)(4.3.2-1)Nj,4,4(u)(

27、)(j=0,1,2,3)=0,1,2,3)是一组重要的基函数，利用它和四个相邻顶点线性组合，可构成三次B B样条曲线段。第26页/共78页(2)(2)三次三次B B样条曲线样条曲线段段 (4.3.2-2)(4.3.2-2)它的端点具有如下的一些性质：(4.3.2-3)(4.3.2-3)根据上述端点几何性质，三次B B样条曲线段的形状就大体确定了。第27页/共78页4.3.3 4.3.3 三次三次B B样条曲样条曲线线当特征多边形的顶点超过四点时，其上每增加一个当特征多边形的顶点超过四点时，其上每增加一个顶点。则相应地在样条上增加一段曲线顶点。则相应地在样条上增加一段曲线,下图表示下图表示B B

28、特征多特征多边形及其对应的边形及其对应的B B样条曲线样条曲线,多边形中每四个相邻的顶点按多边形中每四个相邻的顶点按公式公式（4.3.2-24.3.2-2）定义一段曲线。定义一段曲线。现在从曲线连续、光现在从曲线连续、光滑的要求出发，推导滑的要求出发，推导出三次出三次B B样条曲线方样条曲线方程。推导过程的几何程。推导过程的几何意义十分明显。意义十分明显。已已知知n+2+2个个按按顺顺序序排排列列的的顶顶点点矢矢量量V V0 0,V V1 1,V Vn+n+1 1 (见见上上图图)，设设N N0,40,4(u u),),N N1,41,4(u u),),N N2,42,4(u u),),N N

29、3,43,4(u u),),分分别别为为u u的的三次多项式。三次多项式。顺次以相邻的四顶点顺次以相邻的四顶点V Vi i,V Vi+i+1 1 ,V Vi+i+2 2,V Vi+i+3 3作为一组作为一组,共得到共得到(n-1)-1)个线性组合：个线性组合：(4.3.3-1)(4.3.3-1)第28页/共78页这些线性组合在连接点处要求直到二阶连续，即 ri(1)=ri+1(0)，ri(1)=ri+1 (0)，ri”(1)=ri+1”(0)，由此可推导出三次B样条基函数：（见115-116页）（4.3.3-6)(4.3.3-6)式正是前述的三次B样条基。在在(4.3.3-1)(4.3.3-1

30、)式式中中，特特征征多多边边形形顶顶点点V Vi+ji+j和和三三次次B B样样条条基基函函数数N Nj j,4,4(u u)()(j j=0,1,2,3)=0,1,2,3)线线性性组组合合得得到到r ri i(u u)。当当参参数数u u从从0 0变变化化到到1 1时时，上上式式描描绘绘出出第第i i段段曲曲线线。各各段段曲曲线线在在连连接接点点处处保保持持C C2 2连连续续。由由于于N Nj j,4,4(u u)是是三三次次B B样样条条基基，故故上上述述曲曲线叫三次线叫三次B B样条曲线。它在工程上的应用最普遍。样条曲线。它在工程上的应用最普遍。第29页/共78页将将(4.3.3-6)

31、(4.3.3-6)代入代入(4.3.3-1)(4.3.3-1),写成矩阵形式得三次写成矩阵形式得三次B B样条样条曲线公式为：曲线公式为：（4.3.3-7)此式为三次此式为三次B B样条曲线的计算公式，请记住并会计算。样条曲线的计算公式，请记住并会计算。第30页/共78页4.3.4B样条曲线的几何性样条曲线的几何性质质除了图4.3.2-1所示的曲线端点性质外，B样条曲线还具备另一些性质。下面以三次B样条曲线为例加以说明。但这些性质对任意次B样条曲线也都成立。(1)直观性B B样样条条曲曲线线的的形形状状决决定定于于B B特特征征多多边边形形，而且曲线和多边形相当逼近。而且曲线和多边形相当逼近。

32、(2)(2)局局部部性性 由由于于三三次次B B样样条条曲曲线线段段ri(u)仅仅由由四四个个顶顶点点矢矢量量确确定定,而而与与其其它它顶顶点点矢矢量量无无关关。所所以以改改变变特特征征多多边边形形的的某某一一顶顶点点矢矢量量,只只对对相相邻邻的的四四段段曲曲线线段段产产生生影影响响,而而对对其其它它曲曲线线段段不不会会引引起起变变化化。B B样样条条所所具具备备的的局局部部性性在在B B样样条条曲曲线的几何性质中占有很重要的地位。线的几何性质中占有很重要的地位。(3)(3)凸凸包包性性 所所谓谓空空间间点点的的凸凸包包是是指指连连接接各各空空间间点点所所围围成的空间区域。曲线段必成的空间区域

33、。曲线段必 在在V Vi iV Vi+i+1 1V Vi i+2+2V Vi i+3+3所张所张 成的凸包内。而整条成的凸包内。而整条 B B样条曲线必定落在样条曲线必定落在 这种由相继的四个多这种由相继的四个多 边形顶点所组成的凸包的并集之中。边形顶点所组成的凸包的并集之中。第31页/共78页(4)对称性把特征多边形的顶点V0,V1,Vn反序排成Vn Vn-1,V0按此顶点矢量序列用式(4.3.3-7)构成曲线，就将沿相反的方向描画出同一条曲线。三次B样条曲线的几种退化情况（亦即构造特殊B样条曲线的技巧）：(1)三顶点Vi,Vi+1,Vi+2共线m是线ViVi+2的中点，见图,则三次B样条曲

34、线段的起点ri(0)所在位置为利用端点性质容易证明，利用端点性质容易证明，r ri i(0)(0)处处的的曲曲率率K K=0=0。曲曲线线起起点点ri(0)与与直直线线V Vi iV Vi i+2+2相切，利用此性质可设计出需要的拐点。相切，利用此性质可设计出需要的拐点。第32页/共78页(2)四顶点ViVi+1Vi+2Vi+3共线四个相邻顶点共线，构作的B样条曲线退化为直线，见图4.3.4-3。ri(0)和ri(1)的位置同样可由作图法求出。(3)两顶点Vi+1,Vi+2重合相当于三顶点共线，三次B样条曲线的端点ri(0)满足关系式而且端点曲率K=0,见图4.3.4-4。第33页/共78页(

35、4)三顶点重合为了构作含有尖点的为了构作含有尖点的B B样条曲线，可以取三次重样条曲线，可以取三次重合顶点，即把一个顶点重复取三次。合顶点，即把一个顶点重复取三次。如图如图4.3.4-54.3.4-5所示，由顶点所示，由顶点V Vi i,V Vi+i+1,1,V Vi i+2,+2,V Vi+i+3,3,V Vi i+4+4可以定义两段三次可以定义两段三次B B样样条曲线，具有条曲线，具有C C2 2连续条件。若把连续条件。若把V Vi i+2+2看作为重复的两个顶看作为重复的两个顶点，则由顶点点，则由顶点V Vi i,V Vi+i+1,1,V Vi i+2,+2,V Vi i+2,+2,V

36、Vi+i+3,3,V Vi i+4+4可以构作三段可以构作三段B B样样条曲线。若把条曲线。若把V Vi i+2+2看作为重复三次的顶点，则由顶点看作为重复三次的顶点，则由顶点i i,V Vi+i+1,1,V Vi i+2,+2,V Vi i+2,+2,V Vi i+2+2V Vi+i+3,3,V Vi i+4+4可以定义四段三次可以定义四段三次B B样条曲线。样条曲线。三重顶点三重顶点V Vi i+2+2处，曲线的斜率不连续形成尖点，而且两处，曲线的斜率不连续形成尖点，而且两侧含有直线段和。在尖点处侧含有直线段和。在尖点处斜率尽管不连续。然而对于斜率尽管不连续。然而对于参数曲线来说参数曲线来

37、说,确实是达到确实是达到了了C C2 2连续。因为在三重顶点连续。因为在三重顶点处的一阶和二阶导矢都退化处的一阶和二阶导矢都退化为零。为零。第34页/共78页上述退化情况表明上述退化情况表明(构造特殊构造特殊B B样条曲线的技样条曲线的技巧巧)：(1)如果想构造一段直线如果想构造一段直线,只要使四个顶点共线就可只要使四个顶点共线就可以了。以了。(2)为为了了使使样样条条曲曲线线和和特特征征多多边边形形相相切切，可可以以采采用用三顶点共线或两重顶点的技巧。三顶点共线或两重顶点的技巧。(3)要要使使样样条条曲曲线线通通过过某某一一顶顶点点，即即在在曲曲线线上上使使之之形成一个尖点，可以运用三重顶点

38、的技巧。形成一个尖点，可以运用三重顶点的技巧。在在具具有有复复杂杂外外形形的的机机电电产产品品设设计计中中，常常常常要要用用到到上上述述技技巧。巧。第35页/共78页4.3.5三次三次B样条曲线的算样条曲线的算法法从已知B特征多边形顶点Vi计算三次B样条曲线的结点Pi以及曲线上的任意点，是逼近问题，称为正算，而从已知型值点列Pi反推多边形顶点Vi，是应用于插值的反问题，称为反算。(1)曲线的正算给定特征多边形顶点Vi，构作B样条曲线，按式(4.3.3-7)计算曲线上结点以及任意点的位置矢量，这是不用赘述的。在数控绘图和数控加工中往往需要对参数进行等间隔插值，这时不必将参数代入式(4.3.3-7

39、)逐点计算，而可采用差分运算，使计算可以高速进行。差分法的算法及源程序见附录B。第36页/共78页(2)(2)曲线的反算曲线的反算 在在机机电电产产品品外外形形设设计计实实际际中中，常常常常是是给给出出曲曲线线上上一一批批型型值值点点，希希望望用用B B样样条条曲曲线线来来拟拟合合这这些些点点。然然后后求求出出其其它它需需要要的的插插值值点点。这这时时首首先先要要求求出出B B样样条条曲曲线线特特征征多多边边形形顶顶点点，才才能能构构造造曲曲线线，并并对对曲曲线线进行插值计算。进行插值计算。设 已 知(n n+1)+1)个 有 序 型 值 点 列P Pi i(i i=0,1,=0,1,n n)

40、。求特征多边形顶点位置矢量V Vi i(i i=-1,0,=-1,0,n n+1)+1)。从(4.3.2-3)(4.3.2-3)式的第一、第二式可看出，反算问题归结为下列线性代数方程组的求解：(4.3.5-5)(4.3.5-5)如果补充两个适当的端点条件，方程组就有唯一解。工程中常见的情况有如下三种端点条件：第37页/共78页(1)两端点给出切矢量补充条件为：(4.3.5-6)其中P0和Pn在实际设计中可以依据抛物线端点条件给出。(4.3.5-6)经推导整理后与(4.3.5-5)构成三对角线性方程组：(4.3.5-9)用“追赶法”求解方程组(4.3.5-9)并将结果代入式(4.3.5-6)中求

41、出Vi(i=-1,0,n+1)。全部未知数求解完毕。关于“追赶法”求解方程的原理和程序参见本书附录A。第38页/共78页(2)自由端点条件一般可取(4.3.5-10)由(4.3.5-10)、(4.3.5-5)构成三对角线性方程组(4.3.5-11)可由“追赶法”求解此方程组第39页/共78页(3)(3)封闭曲线封闭曲线为使曲线起点和终点光滑连接,应考虑多生成一小段曲线PnP0将原曲线光滑封闭。考虑Pn点的连续性,PnP0应由顶点Vn-1VnVn+1Vx生成;考虑P0点的连续性，PnP0应由顶点VyV-1V0V1生成,又知PnP0是同一条曲线,所以必有：即要曲线封闭,必须:Vn+1=V0,Vn=

42、V-1(4.3.5-12)由(4.3.5-12)和(4.3.5-5)联立求解，但它不能构成三对角线性方程组，因此不能用解三对角方程的追赶法求解，只能用一般的线性方程组解法求解Vi。方程的形式为(4.3.5-13)第40页/共78页4.3.6 4.3.6 三次参数曲线段的三种等价表三次参数曲线段的三种等价表示示三次参数曲线可以用不同的方法构造,列成附表4.3.6-1。第41页/共78页三次参数曲线段的三种等价表示三次参数曲线段的三种等价表示(续续)三种构作方法有其内在联系。从几何角度分析(见图4.3.6-1所示),B0和B3是曲线段r(u)的始点r(0)和终点r(1),r(0)和r(1)是曲线段

43、始点和终点的切矢量。根据贝齐尔曲线的性质可以验证,从B0开始,沿r(0)的方向截取其模长的13,得B1点；从B3开始,沿r(1)的反方向截取其模长的13,得B2点,则B0B1B2B3即为r(u)曲线段的贝齐尔特征多边形。再根据B样条的端点性质,将线段B1B2向两侧各延长自身的长度,分别得V1和V2。用线段连接V1B0,并延长两倍到d0，再用线段连接V2d0，并延长自身的长度到点V0，在V2和B3方面作对称的操作，得到V3,则V0 V1V2 V3即为r(u)曲线段的B特征多边形。第42页/共78页三者的几何关系是明显的三者的几何关系是明显的,三种几何表示方法可以相互转换三种几何表示方法可以相互转

44、换。由一种几何表示方法很容易用作图法找出另外两种等价的几何表示。三次参数曲线段可以用不同方法构造，表示形式有所不同,但它们之间有其内在联系。从代数上看，这三种表示式可以统一在矩阵形式之下：r（u）=U Mj BjT 0 u 1 （j=1，2，3）其中:U是参数矢量,Mj是基函数阵,BjT 是顶点信息(见前表)。从上述等价表示中可以看出，贝齐尔方法和B样条方法用特征多边形表示曲线，比一般参数曲线更加直观。学习了三次参数曲线段的等价表示,同学们要会解决如下问题：1）给定一曲线段的B特征多边形顶点,如何求出该曲线段的贝齐尔多边形顶点？2）给定一曲线段的贝齐尔多边形顶点，如何求出该曲线段的B特征多边形

45、顶点？第43页/共78页4.3.7二次二次B样条样条曲线曲线在实际应用中，用得最多的是三次B样条曲线，其次就是二次。现在考察二次均匀B样条曲线，其公式为：曲线段的两端点是二次B特征多边形两边的中点，并且以两边为其端点切线(图4.3.7-1)。一次B样条曲线就是B特征多边形本身。对于同一特征多边形而言，随着曲线次数的增高，曲线拉紧，离特征多边形越来越远。三三次次B B样样条条能能保保持持C C2 2连连续续，对对特特征征多多边边形形又又相相当当逼逼近近，所所以以最最常常用用。二二次次B B样样条条曲曲线线由由于于简简单单，与与特特征征多多边边形形更更加加逼逼近近，尽尽管管只只能保持能保持C C1

46、 1连续，在工程中也经常应用。连续，在工程中也经常应用。第44页/共78页4.4 4.4 非均匀非均匀B B样条曲线样条曲线 上节对工程上常用的均匀B B样条曲线作了较全面的介绍。它最重要的特征是：（1 1）参数轴采用均匀分布的节点；（2 2）各段B B样条曲线均采用相同的基函数；（3 3）计算简单。但当型值点或顶点间距相差较大时，应用均匀B B样条曲线效果不好。因因此此，机机电电产产品品设设计计过过程程有有时时要要用用到到非非均均匀匀B B样样条条曲曲线线，它它可可以以在在基基函函数数里里用用节节点点的的不不均均匀匀分分割割来来适适应应特特征征顶顶点点或或型型值值点点严严重重不不均均匀匀的的

47、情情况况。虽虽然然计计算算量量增增加加，但但造造型型功功能能和和效效果果得得以以提提高高，特特别别是是可可以以应应用重节点技术，更增加了用重节点技术，更增加了B B样条曲线的应用灵活性。样条曲线的应用灵活性。第45页/共78页4.4.1B样条基函数样条基函数首先明确几个概念：节点是参数轴上的分割点。将参数轴等距分割,得到均匀节点。若非等距分割,则得到非均匀节点。前面提到过的型值点,是指曲线段间的连接点。前者是曲线上的点，后者是参数轴上的点。而顶点则是确定曲线形状的控制点。定义：在区间a,b上，取分割a=x0 x1xn=b为节点，构造B样条基函数。仅在区间xi x xi+M内其值不为零的M阶(M

48、-1次)B样条基函数Ni,M(x)称为在xi,xi+M上具有局部支集性。M为阶数(为大于或等于1的整数)，基函数Ni,M(x)由下列递推关系给出：式中约定0/0=0，Ni,M(x)称为a,b上第i个M阶（M1次）B样条基函数。第46页/共78页（4.4.1-1）递推公式的几何意义可以归结为“移位”,“升阶”和“线性组合”。Ni,M1（x）“移位”得Ni+1,M1（x）（对非均匀节点而言则指Ni,M1（x）的下一个基函数为Ni+1,M1（x）；M1阶基函数乘以x的线性函数，就是“升阶”；再将其“线性组合”，得到M阶B样条基。由上述递推关系可知，M阶B样条基是一个只在M个子区间xi,xi+M上非零

49、的分段M1次多项式，它具有直到M2阶连续导数。现在具体地讨论0次到3次B样条基当M1，零次(一阶)B样条基为:Ni,1（x）的 图 形 如 图(4.4.1-1)所 示。在 区 间a,b上，它只在一个子区间xi,xi+1上非零，且为常数1(即为零次多项式)。在其他子区间上均为零。Ni,1（x）称为平台函数。第47页/共78页M=2，一次(二阶)B样条基：由(4.4.1-2)的Ni,1（x）的“移位”得将Ni,1（x）和Ni+1,1（x）代入递推公式，得Ni,1（x）和Ni+1,1（x）及由递推所得到的Ni,2（x）见图4.4.1-2。人们形象地称Ni,2（x）为屋顶函数。Ni,2（x）只在两个子

50、区间xi,xi+1,xi+1,xi+2上非零，且各段均为x的一次多项式。第48页/共78页M=3M=3，二次(三阶)B)B样条基：将Ni,2（x）的表达式(4.4.1-3)的下标i换成i+1,即得到Ni+1,2（x）的表达式。再将Ni,2（x）和Ni+1,2（x）代入递推公式，即可得到Ni,3（x）。二次B样条基为:(4.4.1-3)Ni,2（x）、Ni+1,2（x）以及递推所得到的Ni,3（x）见图4.4.1-3。Ni,3（x）在三个子区间上非零,且为分段的二次多项式,形象地被称为钟形函数。第49页/共78页M=4，三次B样条基Ni,4（x）Ni,4（x）是工程上经常用到的一种B样条基，它只