基于四边形网格参数细分的平面与自由曲面求交算法

本资料包含pdf文件1个，下载需要1积分

- 直以来，平面与曲面间求解交线的问题都是计算机辅助设计与制造(computer-aided design/manufae。

turing，CAD/CAM)编程领域中最普遍的工程问题。在曲面造型与裁剪、加工刀具轨迹计算、加工几何图形验证等实际应用中，常常需要对平面与曲面进行求交运算 。国内外在求交计算方面做了大量的研究，大致包括牛顿迭代法H 、曲面离散法 、区间算术法 、光线跟踪法随 以及近年来出现的等值线法n 、三角网格方法n 等。牛顿迭代法通常是在给定初值后将求交问题转化为求线性方程组或常微分方程解的问题，但收敛与否以及收敛的速度与迭代初值的选取有很大关系。曲面离散算法是用小平面片近似逼近曲面的- 种几何求交算法，有较高的可靠性，但缺点是运算量较大、效率低、精度低。区间算术法的提出则只是为了解决直线与隐曲面求交的问题。

鉴于以上各种求交算法各自的缺点，本研究将首先提出-种四边形参数曲面片模型，在此基础上可将平面与自由曲面的求交问题简化为直线段与平面的位置关收稿日期：2013-04-27基金项目：国家自然科学基金资助项目(51275485)；河南省高衅技创新人才支持计划资助项目(13HASTIT036)作者简介：李慧莹(1981-)，女，黑龙江齐齐哈尔人，主要从事多轴数控加工方面的研究.E-mail：lihuiying.1hy###163.corn通信联系人：陈良骥，男，副教授，硕士生导师.E-mail：chenjiaxiclj###163.corn第8期 李慧莹，等：基于四边形网格参数细分的平面与自由曲面求交算法 .957 。

系问题。与现行常见的求交方法相比，本研究算法的特点是计算稳定可靠、精度高而且具有较强的普适性。

1 自由曲面的NURBS表示- 张笛卡尔空间自由曲面可表示为如下非均匀有理B样条形式 n ：S ， ， y(u， ， ， ]T∑∑ ∽ (1)l-------- 0≤u. ≤ 1∑∑ ∽式中： -三维控制点；wi -P 对应权重；(rt1)，(m1)- 向和 向控制点的数目； ， ( ， -沿t向的p次和 向的g次B样条基函数。

由自由曲面的NURBS定义可知，-张自由曲面对应于参数平面内-正方形区域(， )IO≤ ≤1，0≤ ≤1，笛卡尔空间与参数空间的映射关系如图1所示。

图 1 笛卡尔空间与参数空间的映射关系2 空间点、线与平面的关系2.1 空间点 与平面的关系已知平面单位法矢为A，B，c、平面常数为D，平面与曲面相交如图2所示。平面为所有点 ( ，y，z)的集合V1F( )F( ，y， ) ByCzD0。该平面将空间分为3个部分，即 F( )0。因此，空间任-点 与平面有如下关系：图2 平面与曲面相交若点 在平面上 ，则 F(V)-O；反之，若 使得F( )≠0，则 定不在平面上。

2.2 空间线段 与平面的关系空间线段 与平面的关系可由线段两个端点和 分别与平面的关系来确定。具体位置关系可以有以下4种情形：①若 和 分布在平面相同-侧，则线段 与平面不相交，此时F )、F )皆非零且具有相同的符号；②若 和 分布在平面相异- 侧，则线段 I/2与平面交于-点，此时F )、F )皆非零且具有相异的符号；③若 (或 I，2)在平面上而 (或 )不在平面上，则线段与平面交于 (或)；④若 和 都在平面上，则线段位于平面上。

将以上线段与平面之间的关系归纳如下：F )F )>0 与平面不相交F )F )<0 与平面相交于-点)o；F )≠0 与平面交于F )F )0F )≠0；F )0 与平面交于IF )0；F )0 I/2在平面上3 四边形参数面片如图1所示，位于笛卡尔空间的曲面的每个空间四边形都对应于参数平面内的-个小长方形。四边形参数面片模型如图3所示 ，本研究将每个长方形分别按照-定顺序标记出4个角点(1、2、3、4)和4条边(I、Ⅱ、Ⅲ、Ⅳ)，并定义这样的长方形为四边形参数面片。参数平面内，平面与曲面交线的参数变化曲线13 ∽ ，如图4所示。本研究用前-汹中介绍的两种关系可以判定出各个空间四边形每个边与平面的关系，进而可以判断出四边形参数面片与 (u)的位置关系。曲面所有四边形参数面片与 ∽间的位置关系分别会有如图3所示的a、b、c.d、e、厂、g等7种状态。

4Av血 。囟囟图3 四边形参数面片模型本研究对这7种状态分别做如下定义 ：①空间四边形的每个边均在平面的相同-侧为口态；②空间四边形的仅有-个顶点在平面上为b态；③空间四边形的两个相邻边与平面相交为C态；④空间四边形的-个顶点在平面上、-条边与平面相交为d态；⑤空间四边形的两个对边与平面相交为e态；⑥空间四边形的两个对角顶点在平面上为厂态；⑦空间四边形的某边位于平面上时为g态。

· 958 · 机 电 工 程 第30卷4 交线链表生成对图1中的每个小四边形，顺次计算各角点坐标，应用前面所介绍的两种关系和几种状态的定义方法 ，可以得到如图4fi)i：示的-张状态表。

4 口 口 口、 b口 口 口 b 口口 口 6 6 d口 b 口 口口 口 口 d 口cj ，c 口 口 口 口图4 参数变化曲线与状态表求交线的过程为：①建立-张空的四边形面片结构的链表并对照状态表将所有处于非口态和非b态的四边形按-定顺序(如从下到上顺次取出每-行、对取出的每行按从左至右顺次取出每个四边形)存人该链表中，得到-张初始链表；②对初始链表中所有处于非.厂态以及非g态的各结点(小四边形)按同样的方法再进行细分，并将细分后得到的所有非。态和非b态的四边形按与前面相同的顺序插入该链表中，得到- 张新链表；③重复第②步的操作直到链表中所有的结点状态都为 和g态为止。最后得到的四边形链表能很好地逼近交线且交线必然位于平面内，称该链表为交线链表。以图4表示的状态表来说明该计算过程，求交线链表的过程如图5所示。

CCPCd工 ldq始链表逐层细分 线链表图5 求交线链表的过程将交线链表中各四边形位于平面内的角点作为型值点，用NURBS样条进行拟合，即可得到平面与自由曲面的交线。

5 实例计算本节将以-张双三次NURBS曲面为例对所提出的方法进行验证。该曲面由l6个控制点形成，各控制点权重取为1。为说明问题，先找出曲面上已知的3个点P、P。和P ，通过这三点确定-个平面，求这个平面与自由曲面的交线。

取 PS(0.5，0.5)[30，67.5，451T，P s(o.15，0.35)[9，34.528 3，31.5]T，P s(o.85，0.65)[51，34.321 7，58.5 ，经计算，由点P、P。和P 所确定的平面方程各系数分别为：A0.540 757 6，B0.0，C-0.841 178 5，D21.630 304。利用Vc和OpenGL在计算机上编程实现了本研究所提出的求交算法，从初始链表开始，对链表中每个四边形分别进行10x10细分，共经过6层细分可以得到图5中的交线链表，对其进行曲线拟合，求得的平面与自由曲面的交线如图6所示。

图6 平面与自由曲面的交线以下通过求直线与交线的交点来进-步验证算法的精确性。理论上，P，、P 应该在交线上，但实际计算平面内直线 Pl P：与交线的交点时采用将曲线变折线段的方式进行，因而只能得到它们的近似值户，和户 ，求得的 和 分别为：户。[9.003 41，34.528 283 3，31.502 185 31p [50.998 2，34.321 708 86，58.498 842 1T与 Pl和 P2的距离误差分别为：占。 -P。I0.00405mm (2) I户2-P2I0.00214mm (3)6 结束语本研究所提出的平面与自由曲面求交线方法的思路是：根据笛卡尔空间与参数平面之间的映射关系，将空间四边形与平面的位置关系转化为参数平面内四边形与交线参数变化曲线间的关系，进而求得交线。最后，本研究结合实例进行了计算，研究结果表明，求得的交线能够与理论交线精确符合，完全可以在几何造型、曲面裁剪以及五轴编程刀位计算等诸多方面得到实际应用。

(下转第970页)· 970 · 机 电 工 程 第 3O卷h/ram图3 分辨力比较图藏深度2.0 mm的缺陷，而传统方法在计算盲区范围内则基本无法检测到缺陷，图1及图2也说明，本研究方法在检测过程中，测量误差及分辨力均更优于传统方法。

4 结束语为了解决传统TOFD法对埋藏在近表面的缺陷波无法辩识等近表面盲区问题，本研究提出了-种新的近表面检测方法。通过理论分析和实验比较可知，与传统TOFD检测方法相比，基于变型波的方法能有效地检测到埋藏深度 2.0 mm的缺 陷，并有效地解决TOFD近表面盲区问题；该方法通过改变探头中心距和采用变型波检测方法能有效减小测量误差并提高盲区的分辨力；虽然提出的基于变型波的检测方法能有效解决表面盲区的问题，但是由于在波的传播过程中有纵波直通波、底面反射波及变型波等波型，较为复杂，目前要通过人工方法来首先调整探头中心距和人工判断变型波的位置，而智能判别将是下-步研究的重点。