1、大 连 理 工 大 学 硕 士 学 位 论 文- I -硕 士 学 位 论 文检测自由曲面等参数化取样点的方法与研究Isoperimetric line sampling strategy for the inspection of sculptured surfaces作 者 姓 名： 工 程 领 域： 机械工程 学 号： 指 导 教 师： 完 成 日 期： 大连理工大学Dalian University of Technology检 测 自 由 曲 面 等 参 数 化 取 样 点 的 方 法 与 研 究- II -摘 要随着计算机的应用，运用数学方法定义自由曲线/曲面得到实际应用，大大促进
2、了CAGD 的产生和发展。应用 B 样条提供的方法可以依照构造形状的几何信息来建立对应的曲线/曲面方程。这些数学模型可以在计算机上通过被执行计算和处理，提取到曲线/曲面上大量的点和特征信息。在这过程中，通过计算机的分析和综合，还可以实时显示并交互设计修改所定义的形状所拥有的整体和局部信息。这一技术对数控加工、有限元分析、物理性能计算等起到了重要的推动作用。本文重点结合了基于不规则截面数据线的 B 样条曲线/曲面对自由曲线/曲面重构进行了研究。关键词：B 样条曲线；三坐标测量机；检测；自由曲面大 连 理 工 大 学 硕 士 学 位 论 文- III -Isoperimetric line sam
3、pling strategy for the inspection of sculptured surfaces AbstractContents of the abstract. Times New Roman. Key Words：B-spline curve; CMM; Inspection; Sculptured surfacesFair fitting method 光顺拟合 skinning method 逐层截面扫描数据法 sectional points 截面数据 interpolation 插值 computer-aided geometric design 计算机辅助几何设
4、计 parametric surfaces 参数曲面 B-splines B 样条 interaction techniques 交互技术检 测 自 由 曲 面 等 参 数 化 取 样 点 的 方 法 与 研 究- IV -目 录摘 要 .IAbstract .II1 绪论 .11.1 研究背景及意义.11.2 自由曲面检测的研究现状.31.2.1 B 样条曲线/曲面的研究现状 .31.2.2 点到曲线/曲面最小距离的研究现状 .41.2.3 自由曲面检测的研究现状.51.3 论文主要研究内容及结构.62 B 样条曲线和曲面的基本理论 .82.1 形状数学描述的几种方法比较.82.2 B 样条
5、曲线的计算 .92.2.1 B 样条曲线方程及正算反算 .92.2.2 节点矢量的构造.112.2.3 基函数的计算.162.2.4 B 样条曲线插值 .192.3 B 样条曲面的计算 .212.3.1 B 样条曲面的定义 .212.4 本章小结.223 基于截面测量数据光顺拟合 B 样条曲线 .233.1 三次均匀 B 样条 .233.2 最小二乘逼近.243.3 B 样条曲线的光顺拟合算法 .243.2.1 算法思想概述.243.2.2 数学模型.253.2.3 算法步骤.253.4 数值算例.253.4.1 不同线型的算例.253.4.2 开曲线和闭曲线.263.4.3 周期和非周期曲线
6、.26大 连 理 工 大 学 硕 士 学 位 论 文- V -3.4.4 已知条件不同的曲线.263.5 本章小结.264 基于截面测量数据光顺拟合 B 样条曲面 .274.1 B 样条曲面表达 .274.2 B 样条曲面的光顺拟合算法 .274.2.1 算法思想概述.274.2.2 数学模型.274.2.3 算法步骤与实例分析.274.3 数值算例.274.4 本章小结.275 基于轮廓特征的等参数线取样方法在 CMM 测量中的应用 .285.1 CMM 检测策略 .285.2 数值算例.285.3 自由曲线/曲面的误差/光顺性分析 .285.4 本章小结.286 总结与展望 .306.1
7、论文总结.306.2 工作展望.30结 论 .32参 考 文 献 .33攻读硕士学位期间发表学术论文情况 .35致 谢 .36大连理工大学学位论文版权使用授权书.37大 连 理 工 大 学 硕 士 学 位 论 文- 1 -1 绪论机械制造业中涉及到大量的自由曲面的造型、制造与检测，例如压铸的汽车车盖、叶轮的叶片、飞机机翼等（如 图 1.1），这些内容是计算机辅助几何设计（CAGD）领域的研究重点和关键所在。正因为这些自由曲面在 CAD/CAM 领域，特别是航空和汽车工业，有着如此广泛的应用，曲面造型所面临的有关提高精度、检测效率等新问题也愈发突出，因此，对自由曲面的加工质量进行高精度、高效率检
8、测变得越来越重要。图 0.1 自由曲面造型的工业应用Fig. 1.1 Industrial application of sculptured surfaces1.1 研究背景及意义在实际的生产和生活中，常常需要建立三维物体造型，这些造型被广泛应用于计算机辅助设计与制造、计算机动画和艺术、人体模型、地表地形勘察、自然资源分布、还可以应用在军事指挥和机械产品造型方面。对于规则几何形状，可以通过数学表达式描述；对于不规则、表面呈自由曲面的物体，则可以用分割成的曲面片表示，这些曲线曲面的形状不依赖于坐标系的选取，利用直观简便的计算机辅助设计手段，就可以很轻松地离散成造型易于调整、拼接、拟合的几何形状
9、。常用的曲线（曲面）有Bezier 曲线（曲面）和 B 样条曲线曲面。自由曲面是指很难用简单的代数和几何公式表示的曲面，又称雕塑曲面或自由曲面。在实际工程中，常常由于效率和时间等问题无法测得自由曲面上全部的数据，因此需要通过一定的测量，获得自由曲面上的一些离散点集（这些点被要求能反映该曲面的大致轮廓，称这些点为型值点），再通过这些点集逼近生成光滑曲面。这其中涉及到曲面拟合、重构和拼接等多种方法。检 测 自 由 曲 面 等 参 数 化 取 样 点 的 方 法 与 研 究- 2 -目前对自由曲面的高精度检测设备主要有两种：一种是三坐标测量机（Coordinator Measuring Machin
10、e，简称 CMM），另外一种是激光扫描为基础的“点云”数据收集、曲面造型特征求精。而第一种方法因为具有测量精度高、设备造价低、通用性强的特点被广泛应用于工程中。如 图 1.2 所示，CMM 是一种采用点接触测量坐标的仪器，它的工作原理是：沿着理论目标点的法矢方向靠近加工曲面，获得测量的接触点坐标，认为加工误差就是实际测量点和理论值在该方向上的投影距离。因此，为提高检测效率，提高模拟曲面精度，使用 CMM 进行加工曲面检测时，必须首先对自由曲面进行适当的离散处理，构造能使测量逐步逼近过程能实现最终目的的合理的检测样本。判断这个检测样本是否合理的指标有两个：一是应在尽量反映自由曲面的加工情况下，同
11、时还要兼顾检测效率和检测成本。图 0.2 扫描测量Fig. 1.2 Scanning measurement自由曲面（曲线）的测量实际上是利用测量者采集的离散点去准确表达曲面（曲线）的原始几何轮廓的过程：问题的重点在于如何分布曲面上测点的位置和数量，且高效地表达原形状，这个问题部分：（1)当样本数据数量相同时，哪种测点分布算法能最大程度地表达曲面的原始形状； (2)当测量准确度相同时，哪种测点分布算法可以减少采样点数。虽然数据点取得越密集，插值法越具有收敛性，但在工程实践中，都不希望很麻烦。人们希望能用尽可能少但又足以表达形状的数据点，方便地生成所要求的曲线或曲面。由于自由曲面的设计和加工特点
12、，采用 CMM 对自由曲面构建数学模型进行测量时，得到检测点数量检测点数量的确定除了和工艺能达到的精度和检测设备置信度相大 连 理 工 大 学 硕 士 学 位 论 文- 3 -关，还和曲面的面积大小呈一定正相关。因此，检测点的分布应综合考虑曲率变化和离散后曲面片的相互制约关系。考虑到非均匀有理 B 样条在表达曲面上有一定的特点，可以确定研究方法是对复杂自由曲面应采取分片测量，同时在几何量测量中采用优化搜索的方法，在误差范围内确定最小不合格区域，去掉不合格区域后，再继续搜索，直到达到最小距离的最大值，可以基本认为是实现了用尽量少的点最大限度的还原曲面形状。自由曲面加工中的检测取样方法确定不仅可以
13、实现曲面重构，为后续加工曲面提供曲面模型，还能根据微分算法的改进方法求取加工自由曲面模型的误差，通过修改原理论计算所得的刀位轨迹可以实现曲面加工的误差补偿，实验证明，这种误差补偿具有良好的稳定性和准确性。因此，这种技术未来一定会在我国科研领域以及航空航天、汽车、造船和模具等制造行业带来深远的影响。1.2 自由曲面检测的研究现状1.2.1 B 样条曲线/曲面的研究现状20 世纪 80 年代中期以后，随着 CAGD 成为一门应用广泛的新兴学科，B 样条作为该领域内最有发展前景的方法而变成人们研究的重点。1991 年，NURBS（Non-Uniform Rational B-Splines）方法作为
14、国际标准化组织（ International Standardization Organization, 也就是人们熟知的 ISO）规定的定义产品形状的唯一数学方法，用于表示用计算机处理的产品模型形状的数据表示、设计和交换的工业标准。历史上对非均匀 B 样条曲线和曲面有过突出贡献的是 Piegl 和 Tiller1。他们所著的 The NURBS Book 作为一本经典之作详细地介绍了有关曲线和曲面的基本理论，将Bezier-B 样条 -NURBS 作为学科主线，将理论和应用的研究不断推进。国内刘鼎元教授等 2在 1981 年发表了有关 Bezier 曲线和 B 样条曲线光顺拟合法的研究，提出
15、了便于实现交互设计的权因子。在这基础上，刘鼎元又在文献 3中提出 次的 Bezier 曲面拟4合的算法。Woodword 4在 1988 年提出了一种针对截面测量数据进行 B 样条插值曲面的算法，该算法主要针对的是拥有同样个数的截面数据类型，尤其是封闭的自由曲面类型往往容易出现波动、折皱、光顺性差等情况，因此该种方法适用范围也很受限。综合以上文献可知，目前对自由曲线曲面的研究无论是用 B 样条还是 Bezier，无论是曲线还是曲面拟合，其关键主体步骤都是曲面重构看是否在允许的误差范围内，因此构造数学模型，求取合理的定点进行插值拟合，计算拟合精度是本方向研究的关键。目前的研究工作仍面临着以下问题
16、：（1） 点到曲线/曲面的最小距离的算法仍需探讨；检 测 自 由 曲 面 等 参 数 化 取 样 点 的 方 法 与 研 究- 4 -（2） 如何构造合理的曲线/曲面数学模型进行重构曲线 /曲面；（3） 拟合精度的判断标准。1.2.2 点到曲线/曲面最小距离的研究现状重构自由曲线/曲面，其实就要检测自由曲线/曲面拟合过程中和原模型之间的误差是否在精度范围内的问题，这就势必要涉及到求取点到参数曲线和曲面的最小距离的问题。求一个定点到曲线/曲面的最小距离是微分几何学和 CAD/CAM 中的基本问题，在向量空间中，通常解决点到曲线/曲面的最小距离这一问题运用的方法是过该定点向参数曲线/曲面进行空间投
17、影。它的应用范围很广泛，在曲面求交、数控加工刀具轨迹规划、曲面匹配、轮廓度误差评定等领域中都有着十分重要的地位。因此，研究有效且精确的点到参数曲线/曲面最小距离的方法十分必要。牛顿迭代法是最初被提出来并广泛使用的一种求点到参数曲线/曲面最小距离的方法。Mortenson 在 5中曾经提出点到不同类型曲线/ 曲面的最小距离的度量方程，并采用牛顿- 拉斐尔算法找到多项式的根求解方程。文献 6提出求点到参数曲线/ 曲面最小距离最关键的问题是求解非线性方程组，作者提出通过构造辅助方程的方法可解。之后，文献 7提出了一种计算点到曲面的投影的二阶算法作为改进。牛顿型求解点到曲线/ 曲面最小距离的方法具有计
18、算精度高、速度快、算法思想简单易懂、易于编程的特点，因此该类型在 CAD/CAM 领域被广泛采用。但是，这些利用牛顿型的计算方法虽然有这么多优点，但是由于对初始值很敏感，稳定性和有效性并不是那么理想。在求解过程中要想实现对于最优解的收敛必须使每一步迭代的初始值都是良好的，否则就会在逐步迭代的过程中发生离散或者循环的情况，尤其是当目标解位于临近边界位置或端点位置时。为了解决这一问题，一方面，文献 8提供了一种实用的算法用于找到牛顿算法中良好的初值。另一方面，2001 年 Piegl 和 Tiller9给出了一种求取点到非均匀 B 样条曲面（NURBS）投影的方法，该算法的思想是：将非均匀 B 样条曲面分解成多个四边形曲面片，将定点投影到最近的曲面片上，然后根据这个最近的曲面片还原定点的参数值。之后，Ma 和 Hewitt 在文献 10中提出了将 NURBS 曲线分割成 Bezier 曲线段、将NURBS 曲面分割成 Bezier 曲面片的细分算法，这种方法算法的主要步骤有三步：首先分析出被分割出的曲线段或曲面片的控制顶点与定点的关系，之后根据 Bezier 曲线或曲面的强凸包性搜索出待定的曲线段/曲面片，最后通过分别计算定点到待定曲线段 /曲面片上的待定点的距离，比较得出距离最小的点作为该待定曲线段/曲面片上的牛顿型迭代初始点。
