基于NURBS曲面插值的船体曲面重构..doc

基于NURBS曲面插值的船体曲面重构 钱 宏1, 刘 敏2, 贺 庆1, 刘朕明1, 荣焕宗1 （1.中国船舶及海洋工程设计研究院，上海 ；2.中国科学技术大学，） A 0前言 在船体型线设计中，船体线框和曲面设计方法[1]是以B样条曲线为基础，先用B样条曲线插值算法，根据给定的型值点生成横剖线、水线、纵剖线、首尾轮廓线和折角线等，创建船体线框模型。然后用NURBS曲面片覆盖此线框模型重构船体曲面，最后生成船体曲面模型。由于该方法符合设计船舶型线的习惯，船舶设计师容易掌握，在国内还有很大市场。 重构船体曲面的方法主要有孔斯（Coons）曲面、蒙皮曲面和插值曲面方法。孔斯曲面方法[2]采用超限插值法构造满足给定边界跨界矢条件的曲面，使构造出来的曲面片光滑拼合Tribon M3[3]和芬兰NAPA[4]等软件属于这类方法，它们占领着中国市场。蒙皮曲面方法[2]是在关键位置处布置若干条描述曲面形状的截面曲线，生成一张通过这些曲线的曲面这种直观闫秋莲杨启采用一种改进的蒙皮法利用广义站改进了船体曲面的首尾控制网格利用节点删除优化了相容性处理后的节点数目具有工程应用价值采用插值法构造给定的曲面这种NURBS曲面插值方法，直接插值船体线框模型的网格型值点，生成船体NURBS曲面，成功地处理了带切矢的点和重点等复杂情况，实现了用较少的曲面片（几至十几块）重构带有球首、球尾、平边线、平底线、部分折角线和方尾的复杂船体曲面。由于NURBS曲面由较少的曲面片组成，所以可以比较方便地对它进行修改，得到满足设计要求的、光顺的船体曲面。 B样条（B-spline）曲线曲面定义，节点矢量，，，称p次（degree）或p+1阶（order）参数曲线 (1) 为p次B样条曲线。式中，m = n + p + 1，为p次B样条基函数， 为控制顶点。p次B样条曲线的定义域为。为便于统一处理，通常将B样条曲线的定义域取成规范参数域。B样条基函数用deBoor-Cox-Mansfield递推公式计算。 本文只讨论非周期B样条曲线，简称为B样条曲线[2]。其节点矢量形式如下： (2) 1.2 B样条曲面定义 给定(n+1) ?(m+1)个空间矢量，节点矢量与，称p?q次参数曲面 (3) 为p?q次B样条曲面。式中，和为p和q次B样条基函数，为控制顶点。 本文只讨论非周期B样条曲面，简称为B样条曲面[2]。其节点矢量形式如下： (4) B样条曲线插值 为了表示一组数据点（型值点）的位置和切矢，将三维空间的的点附加上切矢控制条件，用ida表示左导数条件，idf表示右导数条件，根据它们可以确定该点的左右切矢。称它们为型值点的特征，可取如下值，切矢长度取数据点的弦长。 ORDINARY：无给定的左右切矢，称为普通点。 DIRLINE：以两点所在直线方向为切方向。若切方向为0，切矢条件退化为重顶点条件。 DIRCURVE：以此点出发的连续三个数据点插值出的抛物线在此点处的切方向为切方向。若不够三个数据点，退化成DIRLINE；若切方向为0，切矢条件退化为重顶点条件。作为端点切矢的默认值。 DIRX：以DIRLINE在X轴的投影为切方向。 DIRY：以DIRLINE在Y轴的投影为切方向。 DIRZ：以DIRLINE在Z轴的投影为切方向。 DIRMIDCURVE：以此点和左右两个点插值出的抛物线在此点处的切方向为切方向。用于光滑拼接。若三个数据点中有且仅有两个相等，退化成DIRLINE；若切方向为0，切矢条件退化为重顶点条件。 INVERSE：表示该方向的切矢与其相反方向切矢相等。用于180o的光滑拼接。 ZERO：切矢取为0。 当型值点中某点有左或右切矢时，用带有2重节点的B样条曲线插值；当型值点中某点有左右切矢时，用带有3重节点的B样条曲线插值，解法详见文献[8]。当型值点中有重型值点时，把重型值点处理成重控制顶点，解法与文献[8]类似，由于篇幅限制，不再详述。一般，两重点会产生明显的凸起，但是，如果在该处加上方向相反的左右切矢条件，就可使之比较平滑；三重点及更高重会产生折线。 B样条曲面插值 为了表示一组网格形式的数据点（型值点）的位置和切矢，将三维空间的一个点附加上切矢控制条件，用ida_u和idf_u分别表示沿着u方向的左、右导数条件，ida_v和idf_v分别表示沿着v方向的左、右导数条件，根据它们可以确定该点的u/v方向左右切矢，称它们为网格型值点的特征，它们的取值见第2节，切矢长度取子网格内对应等参数线的弦长。 设切矢条件总是一行行或一列列地出现，这种有切矢控制的线称之为特征线（边界也可以认为是特征线）。如果某一行/列只有部分数据有切矢控制，根据DIRMIDCURVE规则将该行缺失的特征补全。 当网格型值点中某列/行有左或右切矢时，在u/v方向节点