第三章曲线曲面设计1

1、第三章 曲线曲面设计几何造型技术几何造型技术是一项研究在计算机中，如何表达物体模型形状的技术。描述物体的三维模型有三种: 线框模型、曲面模型和实体模型线框模型、曲面模型和实体模型。l线框模型用顶点和棱边来表示物体。l由于没有面的信息，它不能表示表面含有曲面的物体；l它不能明确地定义给定点与物体之间的关系（点在物体内部、外部或表面上）。l表面模型用面的集合来表示物体，而用环来定义面的边界。l表面模型能够满足面面求交、线面消隐、明暗色彩图、数控加工等需要。l但在该模型中，只有一张张面的信息，物体究竟存在于表面的哪一侧，并没有给出明确的定义，无法计算和分析物体的整体性质。如物体的表面积、体积、重心等

2、。l也不能将这个物体作为一个整体去考察它与其它物体相互关联的性质，如是否相交等。l实体模型能完整表示物体的所有形状信息，可以无歧义地确定一个点是在物体外部、内部或表面上。是最高级的模型。l这种模型能够进一步满足物性计算、有限元分析等应用的要求。三维表面模型表示三维物体的信息并不完整，但它能够表达复杂的雕刻曲面，在几何造型中具有重要的地位，对于支持曲面的三维实体模型，表面模型是它的基础几何造型的历史l曲面造型：60年代，法国雷诺汽车公司、Pierre Bzier、汽车外形设计的UNISURF系统。l实体造型：1973英国剑桥大学CAD小组的Build系统、美国罗彻斯特大学的PADL-1系统等。l

3、独立发展起来，又合二为一。l主流：基于线框、曲面、实体、特征统一表示的造型设计系统3.1 参数曲线和曲面3. .1. .1 曲线曲面参数表示曲线曲面参数表示l显式表示:y=f（x）l隐式表示:f（x，y）=0l参数表示:P(t)=x(t), y(t), z(t)显式或隐式表示存在下述问题：（1）与坐标轴相关；（2）会出现斜率为无穷大的情形(如垂线)；（3) 不便于计算机编程。参数表示:曲线上任一点的坐标均表示成给定参数的函数。假定用t表示参数，平面曲线上任一点P可表示为： 空间曲线上任一三维点P可表示为：)(),()(tytxtP)(),(),()(tztytxtP参数表示例子：l直线l圆 1

4、 , 0,)()(121ttPPPtP 1 , 012,11)(222ttttttP参数表示的优点：（1）以满足几何不变性的要求。（2）有更大的自由度来控制曲线、曲面的形状（3）对曲线、曲面进行变换，可对其参数方程直接进行几何变换。（4）便于处理斜率为无穷大的情形，不会因此而中断计算。（5）便于用户把低维空间中曲线、曲面扩展到高维空间去。（6）规格化的参数变量t0, 1，使其相应的几何分量是有界的，而不必用另外的参数去定义边界。（7）易于用矢量和矩阵表示几何分量，简化了计算。3.1.2 位置矢量、切矢量、法矢量、曲率和挠率l曲线上任一点的位置矢量可表示为： P(t)=x(t), y(t), z

5、(t)；位置矢量切向量(切矢量)l切向量(切矢量)l选择弧长s作为参数，则 是单位切矢l根据弧长微分公式有：l于是有 ，即为单位矢量sPdsdPTs0lim2222dzdydxds22222)(/tPdtdzdtdydtdxdtds0)(tPdtds)()(tPtPdsdtdtdPdsdP法矢量 T(S)与 垂直l法矢量l与 平行的法矢称为曲线在该点的主法矢主法矢N Nl矢量积 是第三个单位矢量，它垂直于T和N。把平行于矢量B的法矢称为曲线的副副法矢量法矢量l我们可以推导出：dsdTNTB)()()()()()()()()()(tPtPtPtPtPtPTBNtPtPtPtPB 密切面从切面法平

6、面TBN主法线图3.1.2 曲线的法矢dsdTlT(切矢)、N(主法矢)和B(副法矢)构成了曲线上的活动坐标架lN、B构成的平面称为法平面，N、T构成的平面称为密切平面，B、T构成的平面称为从切平面。密切面从切面法平面TBN主法线图3.1.2 曲线的法矢 l曲率和挠率 即称为曲率曲率，其几何意义是曲线的单位切矢对弧长的转动率曲率k的倒数 称为曲率半径曲率半径。挠率挠率 的绝对值等于副法线方向(或密切平面)对于弧长的转动率.ssTsTTss00limlimss0lim1sslim.对于一般参数t，我们可以推导出曲率和挠率的计算公式如下：3)()()(tPtPtP 2)()()(),(),(tPt

7、PtPtPtP )(ssT)(sTTO图3.1.3 曲率和挠率（a）（b）1N1B1T0N0B0T0B1B3.1.3 插值、拟合、逼近和光顺给定一组有序的数据点Pi，i=0, 1, , n，构造一条曲线顺序通过这些数据点，称为对这些数据点进行插值插值，所构造的曲线称为插值曲线。 l线性插值：假设给定函数f(x)在两个不同点x1和x2的值，用一个线形函数：y=ax+b，近似代替，称为的线性插值函数。l抛物线插值:已知在三个互异点 的函数值为 ，要求构造一个函数 使抛物线 在结点 处与 在 处的值相等cbxaxx2)()(x321,xxx321,yyy) 3 , 2 , 1( ixi)(xfixx

8、yo1y2y)(xfy )(xy1x2xxyo1y2y)(xfy )(xy1x2x3x3y(a)(b)图3.1.4 线性插值和抛物插值拟合：拟合：构造一条曲线使之在某种意义下最接近给定的数据点(但未必通过这些点)，所构造的曲线为拟合拟合曲线。在计算数学中，逼近通常指用一些性质较好的函数近似表示一些性质不好的函数。在计算机图形学中，逼近继承了这方面的含义，因此插值和拟合都可以视为逼近。光顺(Firing)指曲线的拐点不能太多。对平面曲线而言，相对光顺的条件是：la. 具有二阶几何连续性(G2)；lb. 不存在多余拐点和奇异点；lc. 曲率变化较小。3.1.4 参数化过三点P0、P1和P2构造参数

9、表示的插值多项式可以有无数条，这是因为对应地参数t, 在0, 1区间中有无数种取法。即P0、P1和P2可对应不同的参数值，比如， 或 其中每个参数值称为节点(knot)。, 1,21, 0210ttt, 1,21, 0210ttt对于一条插值曲线，型值点 与其参数域 内的节点之间有一种对应关系。对于一组有序的型值点，所确定一种参数分割，称之这组型值点的参数化。nPPP,10,0nttt 参数化常用方法有：l均匀参数化(等距参数化)l节点在参数轴上呈等距分布， +正常数。l累加弦长参数化 这种参数法如实反映了型值点按弦长的分布情况，能够克服型值点按弦长分布不均匀的情况下采用均匀参数化所出现的问题

10、。iitt1 niPtttiii, 2 , 1,0110iiiPPP1l向心参数化法 向心参数化法假设在一段曲线弧上的向心力与曲线切矢从该弧段始端至末端的转角成正比，加上一些简化假设，得到向心参数化法。此法尤其适用于非均匀型值点分布。niPtttiii, 2 , 1,021110 l修正弦长参数化法弦长修正系数Ki=1。从公式可知，与前后邻弦长及相 niPKtttiiii, 2 , 1,0110iiiiiiiiiPPPPPPK112122310,2,min111niiiiPPPPP比，若越小，且与前后邻弦边夹角的外角i-1和 i(不超过时)越大，则修正系数就K i 就越大。参数区间的规格化我们

11、通常将参数区间 规格化为0, 1， ，只需对参数化区间作如下处理： , 0nttnittttnii, 1 , 0, 00 1 , 0, 0ntt3.1.5 参数曲线的代数和几何形式我们以三次参数曲线为例，讨论参数曲线的代数和几何形式。代数形式l上述代数式写成矢量式是zzzzyyyyxxxxatatatatztatatatatyatatatatx012233012233012233)( 1 , 0)()( 1 , 0)(012233tatatatatP几何形式l对三次参数曲线，若用其端点位矢P(0)、P(1)和切矢P(0)、P(1)描述。l将P(0)、P(1)、P(0)和P(1)简记为P0、P1

12、、P0和P1，代入 得 1 , 0)(012233tatatatatP1010310102010022233PPPPaPPPPaPaPa 1 , 0)()2()32() 132()(123023123023tPttPtttPttPtttPl令：可将其简化为：上式是三次Hermite(Ferguson)曲线的几何形式几何形式，几何系数几何系数是P0、P1、P0和P1。 称为调和函数调和函数（或混合函数）1010,GGFF132)(230tttF23132)(tttFttttG2302)(231)(tttG 1 , 0)(11001100tPGPGPFPFtP 3.1.5 Ferguson曲线端点

13、位矢和切矢0P0P1P1P)(tP)(tP0Fto11to111Fto11to110G1G图3.1.6 三次调和函数3.1.6 连续性曲线间连接的光滑度的度量有两种：l函数的可微性：组合参数曲线在连接处具有直到n阶连续导矢，即n阶连续可微，这类光滑度称之为 或n阶参数连续性。l几何连续性：组合曲线在连接处满足不同于 的某一组约束条件，称为具有n阶几何连续性，简记为 。nCnCnG反例：21,3)(213 10,3)(01010010tVVtVVVttVVVt0131)1 (VV 0132)1 (VV 若要求在结合处达到 连续或 连续，即两曲线在结合处位置连续：若要求在结合处达到 连续，就是说两

14、条曲线在结合处在满足 连续的条件下，并有公共的切矢 当a1时， 连续就成为 连续若要求在结合处达到 连续，就是说两条曲线在结合处在满足 连续的条件下，并有公共的曲率矢：0G0C)0() 1 (QP1G0G)0() 1 ()0(PQ1G1C2G1G这个关系可写为： 为任意常数。当 ， 时， 连续就成为 连续。33)0()0()0() 1 () 1 () 1 ( QQQPPP) 1 () 1 ()0(2PPQ102G2C)(tP)(tQ) 0(P) 1 (P) 0(Q) 1 (Q图3.1.7 两条曲线的连续性我们已经看到， 连续保证 连续， 连续能保证 连续，但反过来不行。也就是说 连续的条件比 连续的条件要苛刻。1G1C2G2CnCnG3.1.7 参数曲面基本概念一张定义在矩形域上的参数曲面可以表示为 可记为 1 , 0 1 , 0),(,),().(),(vuvuzzvuyyvuxx),(),(),(),(vuzvuyvuxvuP示意图 P00P00P00P00 xyzw=wjPw(ui,wj)n(ui,wj)Pu(ui,wj)P(ui,wj)u=0w=0u=1