一种三角网格模型的分治加工方法
【专利摘要】本发明提供一种三角网格模型的分治加工方法，属于CAM【技术领域】，其特征在于：选取R半径范围邻域点以精确计算三角面片模型的微分几何信息；以加工子区域的特征表述为生长原则将三角面片模型分割成不同特征的子区域，进行小面积区域和误判区域的优化合并消除过分割现象，对子区域进行边界光顺以减少锯齿状边界；对不同类型加工子区域采用不同的轨迹策略，生成等残留高度轨迹时，对凸子区域与凹子区域采用环切初始轨迹生成方法，对鞍形子区域采用行切初始轨迹生成方法，刀具轨迹投影偏置扩展时，进行子区域边界处轨迹的处理，得到合理完整的子区域加工轨迹。该方法能够由复杂三角面片模型有效地生成兼顾加工效率和加工质量的数控加工刀具轨迹。
【专利说明】一种三角网格模型的分治加工方法
【技术领域】
[0001]本发明涉及CAM【技术领域】，尤其涉及一种三角网格模型的分治加工方法。
【背景技术】
[0002]随着反求工程、快速成形与计算机辅助设计等技术的发展，三角网格模型在近十几年里得到了广泛应用，与其相关的轨迹规划研究和网格分割研究也成为研究热点。从三角网格模型直接生成刀具轨迹，可以避免曲面重构过程精度的损失，也可消除利用3D打印方法加工金属零件对成本及材料的限制。迄今为止，针对参数曲面模型的三轴数控加工刀具轨迹规划研究相对较多；而基于网格模型的数控加工刀具轨迹规划研究还较少，且现有研究基本上是将整个网格模型看成单一特征，采用一种轨迹生成策略生成整个加工表面的刀具加工轨迹。单一轨迹生成策略为了保证加工精度，一般要采用较小的轨迹参数，因此采用分治加工策略来兼顾加工效率和加工质量。
[0003]所谓三角网格模型的分治加工，即针对复杂三角网格模型不同区域的几何特点，分别采用不同的轨迹生成策略来生成数控加工轨迹，使得在各子区域内的加工轨迹具有相同轨迹类型、在误差允许条件下具有相同或相近的、且尽可能大的轨迹行距和步距，以便达到兼顾加工效率和加工质量的目的。然而，目前三角网格模型分治加工的研究几乎没有深入开展，只有郁萌在其硕士学位论文“复杂三角网格模型数控加工特征提取及刀具轨迹生成研究”中研究了三角网格模型局部干涉特征和狭槽特征的识别分割及其刀具轨迹生成，并没有对三角网格模型的其它特征进行研究。三角网格模型分治加工研究包括两个重要的方面，即三角网格模型的加工子区域分割以及各子区域的轨迹生成。在三角网格模型的分割方面，由于网格模型的复杂性、多样性与分割目的不同，到目前为止还没有一种公认通用的分割准则。根据不同的分割目的，网格模型可以分割成具有几何意义区域、语义意义区域与加工意义区域等，如Wang J等在学术期刊《Computers&Graphics》2011，35 (3), P661-667 发表的学术论文 “Surface feature based mesh segmentation，，、Xiao D 等在学术期刊《Computers&Graphics》2011, 35(3)，P685-691 发表的学术论文“CADmesh model segmentation by clustering”等就是将模型分割为有几何意义区域。但这些分割都不是针对分治加工，因而没有考虑分割子区域的刀具轨迹是否具有相似性的问题。分治加工的另外一个方面是由分割子区域生成有效的刀具轨迹，而刀具轨迹的不同对加工精度和加工效率有重要的影响。目前，针对三角网格模型的刀具轨迹生成方法主要有如黄常标等在学术期刊《机床与液压》2009, 37 (12) P1-4发表的学术论文“三角网格模型等高线刀具轨迹生成研究”中采用的截面线法[5]、陈晓兵等在学术期刊《中国机械工程》2013，24 (08)，P1047-1051发表的学术论文“三角网格曲面等参数线刀轨生成算法”采用的参数化法、以及徐金亭等在学术期刊《机械工程学报》2010，46(11),P193-198发表的学术论文“基于网格曲面模型的等残留刀位轨迹生成方法”和陈晓兵等在学术期刊《计算机辅助设计与图形学学报》.2009, 21 (12)，P1800-1804发表的学术论文“三角网格表面等残留高度刀轨生成算法”等等采用的等残留高度法。截面线法只适用于简单模型，对于复杂的三角网格模型无法生成等残留刀具轨迹；而参数化法是将模型映射到平面后再生成轨迹，然后再逆映射到模型表面，计算过程较为复杂。等残留高度法根据最大允许残留高度直接计算各刀触点行距，能够生成更加精确的加工轨迹。但是目前的等残留高度法研究，主要集中在基于整体模型的刀具轨迹偏置生成方法，都是以模型边界作为初始轨迹，这并不能很好解决分治加工中不同类型子模型轨迹生成问题。
[0004]综上所述，目前还没有成熟方法来有效实现三角网格模型的分治加工，现有三角网格模型的分割方法、轨迹生成方法等相关技术也都无法直接或有效地应用于实现三角网格模型分治加工。
【发明内容】
[0005]本发明的目的在于提供一种三角网格模型的分治加工方法，可以将复杂三角网格曲面模型分解为多个子区域，使得各子区域内能采用相同或相似的轨迹，并能够根据子区域特点生成合适的刀具轨迹，最终生成兼顾加工效率和加工质量的加工轨迹。
[0006]为了解决上述的技术问题，本发明提供了一种三角网格模型的分治加工方法，包括以下几个步骤:
[0007](I)读取三角网格模型进入内存，去除冗余顶点，建立基于半边结构的拓扑信息；
[0008](2)选取R半径范围邻域点后，采用局部二次曲面拟合法来精确计算三角网格模型的微分几何信息，包括主曲率和主方向；
[0009](3)根据三角网格模型的形状与微分几何信息将其分成六类自定义的加工特征子区域，分别为陡壁子区 域、平缓子区域、局部干涉子区域、凹子区域、凸子区域和鞍形子区域；并以最大主曲率和最小主曲率构成平面坐标系来实现所述加工特征子区域的特征表述；
[0010](4)采用以所述加工特征子区域的特征表述为生长原则的区域生长法，将三角网格模型分割成不同特征的加工特征子区域；对分割后产生过分割现象，进行小面积区域和误判区域的优化合并处理，消除面积很小的区域和误判的过渡区域；再对呈现出较多锯齿状边界的所述加工特征子区域进行边界光顺；
[0011](5)对不同类型加工子区域采用不同的刀具轨迹策略来生成刀具轨迹:
[0012]陡壁子区域采用等高线切法；
[0013]平缓子区域采用等行距行切法；
[0014]凹子区域与凸子区域采用等残留高度环切法；
[0015]鞍形子区域采用等残留高度行切法；
[0016](6)将所述初始刀具轨迹用投影偏置扩展的方法生成整个子区域的加工轨迹，并对所述整个子区域的加工轨迹的边界进行处理，从而得到合理完整的子区域加工轨迹。
[0017]作为优选:所述的R半径范围邻域点具体是指:
[0018]以任意面片F的中心点O为起始点，在三角网格模型表面上与所述中心点O的弧线距离小于等于R的所有顶点的集合。
[0019]作为优选:所述加工特征子区域的特征表述具体是指:
[0020]以最大主曲率K1为横坐标、最小主曲率K2为纵坐标建立主曲率坐标系,所述主曲率坐标系内由于!^ SK2，所有的坐标点都在直线K1=K2上或其右侧；根据不同情况设定Z轴分量阈值V。、平缓曲率阈值K。、凹区域曲率阈值1(,、凸区域曲率阈值Kn等阈值后，所述六类自定义的加工特征子区域的特征表述具体是:
[0021]①陡壁子区域:面片法矢的Z轴分量K V。，其中V。为Z轴分量阈值；
[0022]②平缓子区域=K1与K2的平方和? + O C，其中K。为平缓曲率阈值；
[0023]③局部干涉子区域=K1X).9KT，其中Kt为球头铣刀顶端的曲率，0.9为避免局部干涉的安全系数；
[0024]④凹子区域且要求Κ2>Κτ ,或Κ2〈Κτ且K2M^K1，其中Kτ为凹区域曲率阈值；
[0025]⑤凸子区域:⑷+ K;-)>K:,且要求Κ,Κη，或Κ,Κη且Κ^-0.2K2，其中Kn为凸区域曲率阈值；
[0026]⑥鞍形子区域=K1XK K2〈0且不属于上述各区域。
[0027]作为优选:所述小面积区域和误判区域的优化合并处理的具体做法是:
[0028]①找出需要被合并的区域记为R1，其中，
[0029]小面积区域的判别方法是:计算区域内所有面片面积的总和Α，并判断其是否超过给定面积阈值Α。；若Α〈Α。，则判定其为小面积区域；
[0030]误判区域的判别方法是:先找出凸特征区域或凹特征区域，并找到所述凸特征区域或凹特征区域中的最值点；判断所述最值点是否为所述凸特征区域或凹特征区域的极值点；若不是，则将所述凸特征区域或凹特征区域判定为误判区域；所述极值点的判断方法是:凹特征区域的最小值点或凸特征区域的最大值点是否在区域边界上或接近边界，若是，则该最值点就不是极值点；
[0031]②搜索R1的相邻区域:沿着R1的边界搜索其相邻区域并计算合并参数，若对应相邻区域的合并参数最大，则将该区域作为合并区域R2，其中所述合并参数为:
[0032]小区域的合并参数是其边界垂直方向的法曲率半径的平均值；
[0033]误判区域合并参数包括两个方面，其一是相邻边界的长度；其二为特征类型要求，即若R1为凹特征，则R2为凸特征；若R1为凸特征，则R2为凹特征。
[0034]③将所述需要被合并的区域R1合并到最合适的所述合并区域R2,分为两部分:
[0035]一是边界的合并，将两区域的相邻边界去除，然后把所述需要被合并的区域Rl的剩余边界作为所述合并区域R2边界的一部分；
[0036]二是面片的合并，将所述需要被合并的区域R1的所有面片归为所述合并区域民的一部分。
[0037]作为优选:所述对呈现出较多锯齿状边界的所述加工特征子区域进行边界光顺的具体做法是:
[0038]以边界分叉点为端点，将区域边界分成若干条边界线段；所述边界分叉点是指与其相邻区域超过两个的边界点；若存在无所述边界分叉点的区域边界，则将整条所述区域边界作为一条边界线段；然后以所述边界线段为处理路径，逐条进行光顺处理，具体步骤是: [0039]①取一条所述边界线段中的一个边界点P1，搜索所述边界点P1的一阶邻域面片，并区分属于不同区域的面片；[0040]②计算以所述边界点P1为顶点的边界线两个夹角α与β，其中α〈 β，当比值β /α>λ时，则判定该局部边界不光顺，其中λ为夹角比值阈值；
[0041]③若局部边界被判定为不光顺，则把与所述边界点P1相邻的两条所述边界线段从所述区域边界中删除，同时将与所述边界点P1相对的边界线段E作为新的区域边界；
[0042]④将处理对象转至下一边界点P2,重复步骤①-③，直到所有边界点都处理完毕为止。
[0043]作为优选:所述生成等残留高度环切的初始刀具轨迹的具体做法是:
[0044]①找出所述特征子区域的极值点Pl ；
[0045]②按照最大允许残留高度计算行距的方法计算所述极值点匕周围的加工刀触点在某一邻边E1方向行距L，将所述极值点匕沿所述邻边E1向外的单位矢量ei方向偏置L/2距离得到刀触点Pi=P1+!^.θ!；
[0046]③取θι与所述邻边E1顺时针方向相邻邻边向外方向的单位矢量e2之和的角平分线方向e3，即e3=ei+e2 ;将e3单位化，用②中方法求e3方向的相应刀触点；
[0047]④重复步骤②-③直到所述极值点匕在所有邻边方向都计算完毕，将上述求得刀触点依次放入刀触点容器中，得到环切初始刀具轨迹。
[0048]作为优选:所述生成等残留高度行切的初始刀具轨迹的具体做法是: [0049]①求所述特征子区域的包围盒,过所述包围盒的中心点沿刀轴矢量方向做一条射线，求该射线与三角网格曲面的交点P，则得到以所述交点P为初始轨迹的一个刀触点，将所述交点P放入刀触点容器中；
[0050]②将所述交点P作为当前处理的刀触点，求出所述交点P的最小主方向T2,过所述交点P以叉积T2XVt为法矢做平面，其中Vt为刀轴矢量。将π与所述交点P所在面片求交，得到两交点PpP2 ;
[0051]③将所述交点P1插入刀触点容器末端，把所述交点P1作为当前处理刀触点，求所述交点P1的最小主方向T2，过所述交点P1以叉积T2XVt为法矢做平面JI，用JI与所述交点P1的相邻面片F1求交，将求得的非P1交点P3为下一刀触点；
[0052]④将所述交点P3插入刀触点容器末端，把所述交点P3作为当前处理刀触点，重复步骤3，直到求得的刀触点超出所述特征子区域的包围盒；
[0053]⑤将所述交点P2插入刀触点容器前端，用步骤③-④中的方法求刀触点，将求得刀触点插入刀触点容器前端，直到求得的刀触点超出所述特征子区域的包围盒；最终得到行切初始刀具轨迹。
[0054]作为优选:所述对整个子区域的加工轨迹的边界进行处理的具体做法是:
[0055]①当所述环切初始刀具轨迹的扩展达到边界时，先将完整的环切刀具轨迹复制备份到一个清空的交点数组中；然后处理所述环切刀具轨迹，删除其中位于子区域外的刀触点，保留处于区域内的刀触点作为有效的环切刀具轨迹；后续的环切刀具轨迹则基于备份的所述完整的所述环切刀具轨迹来扩展生成；
[0056]②对于所述行切刀具轨迹，采用偏置扩展法，由上一条备份的所述行切刀具轨迹获得一条首末点超出所述特征子区域包围盒的未处理的行切刀具轨迹；将所述未处理的行切刀具轨迹复制备份到一个清空的交点数组中；将所述未处理的行切刀具轨迹中位于所述特征子区域边界以外的顶点删除，使每条行切刀具轨迹的首末点处于所述特征子区域边界。
[0057]本发明与现有相关技术相比，具有以下优点:
[0058](I)本发明提供的三角网格模型分治加工方法能够有效实现三角网格模型的分治加工，生成兼顾加工效率和加工质量的刀具轨迹；
[0059](2)基于R半径范围邻域点来进行局部二次曲面拟合，能够避免因模型网格疏密程度不同而导致邻域点所覆盖的区域大小相差较大的问题，从而保证曲率计算的准确性；
[0060](3)采用以加工子区域的特征表述为生长原则的区域生长法，并在分割后进行小面积区域和误判区域的优化合并处理以及边界光顺处理，可以有效实现用于分治加工的三角网格模型分割；
[0061](4)不同特征采用不同的刀具轨迹策略，特别是针对凸子区域与凹子区域、鞍形子区域，采用不同方法生成等残余高度轨迹的刀具轨迹，并对边界处的刀具轨迹进行了特殊处理，可以有效实现三角网格模型分割后的各类型子区域的刀具轨迹的生成。
【专利附图】
【附图说明】
[0062]图1是本发明程序流程图；
[0063]图2是本发明R半径范围邻域点示意图；
[0064]图3是本发明加工子区域的特征表述示意图；
[0065]图4是本发明采用的区域生长法示意图；
[0066]图5是本发明采用的区域生长法的实现流程；
[0067]图6是本发明误判区域示意图；
[0068]图7是本发明子区域边界光顺示意图；
[0069]图8是本发明进行子区域边界光顺后的示意图；
[0070]图9是本发明用行切法生成初始刀具轨迹中极值点偏置示意图；
[0071]图10是本发明用行切法生成初始刀具轨迹中偏置后刀触点示意图；
[0072]图11是本发明采用的最大允许残留高度计算行距的方法示意图；
[0073]图12是本发明用行切法生成初始刀具轨迹示意图；
[0074]图13是本发明采用的初始刀具轨迹投影偏置扩展方法示意图；
[0075]图14是本发明环切边界刀具轨迹处理示意图；
[0076]图15是本发明行切边界刀具轨迹处理示意图；
[0077]图16是本发明实施实例三角网格模型的立体显示线框图；
[0078]图17是本发明实施实例进行加工子区域分割后的显示效果；
[0079]图18是本发明实施实例加工子区域优化合并后的显示效果；
[0080]图19是本发明实施实例加工子区域边界光顺后的显示效果；
[0081]图20是本发明实施实例各子区域轨迹生成后的顶视图效果；
[0082]图21是本发明实施实例各子区域轨迹生成后的立体显示效果；
[0083]图22是本发明实施实例凸子区域加工轨迹的局部放大显示效果；
[0084]图23是本发明实施实例鞍形子区域加工轨迹的局部放大显示效果；
[0085]图24是本发明实施实例平缓子区域加工轨迹的局部放大显示效果；
[0086]图25是本发明实施实例陡壁子区域加工轨迹的局部放大显示效果。【具体实施方式】
[0087]下文结合附图和实施例对本发明做进一步说明。
[0088]如图1所示，本发明提供的一种三角网格模型的分治加工方法，具体采用以下几个步骤:
[0089]步骤一、读取三角网格模型进入内存，去除冗余顶点，建立基于半边结构的拓扑信息；
[0090]从硬盘中将三角网格模型，如STL模型，读入内存时，采用C++标准模板库中的关联容器类，如set、map、multimap等,进行冗余顶点的去除,建立无重复顶点的vector类型点表和面片表；在此基础上，同样采用关联容器类来创建半边数据结构，其方法是:
[0091]①从面片表中取一个三角面片，分别创建半边结构的三个对象，并依次处理三个半边对象；
[0092]②对于面片的每一个半边对象，则就建立新的半边映射，并搜索多映射Multimap中键值等于新半边的终点的半边映射；
[0093]③若步骤②中的搜索结果为O个半边映射，则把新半边映射存入多映射Multimap中，以待下次被搜索，返回重复步骤②； [0094]④若步骤②中的搜索结果为η (η>=1)个半边映射，则在这η个半边映射中搜索，查找是否有半边对象终点等于新半边对象的起点；
[0095]⑤若步骤④中有搜索到对应的半边对象，则将该半边对象的伙伴半边指针指向新半边对象，同样将新半边的伙伴半边指针指向该半边对象；
[0096]⑥若步骤④中没有搜索到对应的半边映射，则将新半边映射存入Multimap中，以待下次被搜索，返回重复步骤②；处理完一个面片后返回步骤①，当模型中面片全部处理完后，则就结束伙伴半边对应关系建立。
[0097]步骤二、选取R半径范围邻域点后，采用局部二次曲面拟合法来精确计算三角网格模型的微分几何信息，包括主曲率和主方向；
[0098]所述的R半径范围邻域点是指:以面片任意F的中心点O为起始点，在三角网格模型表面上与点O的弧线距离小于等于R的所有顶点的集合。
[0099]如图2所示，以面片F中心点O为圆心、以指定R为半径的虚线圆内顶点为R半径范围邻域点，用于拟合二次曲面以计算三角网格模型的微分几何信息，具体的局部拟合的二次曲面为:
[0100]S(u, v) = (u, v, h(u, v))
[0101]式中h(u, v)=au2+buv+cv2+eu+fv+g, (u, v, h)为拟合邻域点(记为 Pk, k=3,4,5，……，n)的局部坐标值。利用最小二乘法求得二次曲面参数a-g后，再计算出曲面S(u，v)在顶点处的第一、第二基本量，进而利用Weingarten变换矩阵求解中心点O的最大主曲率K1与最小主曲率K2、以及与两个主曲率相对应的主方向。
[0102]步骤三、根据三角网格模型的形状与微分几何信息将其分成六类自定义的加工特征子区域，分别为陡壁子区域、平缓子区域、局部干涉子区域、凹子区域、凸子区域和鞍形子区域；并以最大主曲率和最小主曲率构成平面坐标系来实现所述加工特征子区域的特征表述；[0103]所述步骤3中的加工特征子区域的特征表述，具体是指:
[0104]如图3所示，以最大主曲率K1为横坐标、最小主曲率K2为纵坐标建立主曲率坐标系，坐标系内由于K1 SK2，所有的坐标点都在直线K1=K2上或其右侧；根据不同情况设定Z轴分量阈值V。、平缓曲率阈值K。、凹区域曲率阈值1(,、凸区域曲率阈值Kn等阈值后，六类自定义加工特征子区域的特征表述具体是:
[0105]①陡壁子区域:面片法矢的Z轴分量K V。，其中V。为Z轴分量阈值；
[0106]②平缓子区域=K1与K2的平方和μ?2，其中K。为平缓曲率阈值；
[0107]③局部干涉子区域=K1X).9ΚΤ，其中Kt为球头铣刀顶端的曲率，0.9为避免局部干涉的安全系数；
[0108]④凹子区域+，且要求1(2>1^，或Κ2〈Κτ且K2XJK1，其中Kτ为凹区域曲率阈值；
[0109]⑤凸子区域+ K;)> K;,且要求K1Iη，或?ρκη且&〈-().2Κ2，其中Kη为凸区域曲率阈值；
[0110]⑥鞍形子区域=K1XKKyo且不属于上述各区域。
[0111]步骤四、采用以加工特征子区域的特征表述为生长原则的区域生长法，将三角面片模型分割成不同的加工特征子区域；对分割后产生过分割现象，进行小面积区域和误判区域的优化合并处理，消除面积很小的区域和误判的过渡区域；再对呈现出较多锯齿状边界的所述加工特征子区域进行边界光顺；
[0112]所述步骤4中以加工特征子区域的特征表述为生长原则的区域生长法，如图4所示示意图，图4中阴影部分与非阴影部分分别代表区域生长过程中的已生长区域与未生长区域。其中B1与B2为区域边界半边，F为待生长面片,E1与E2为F不与已生长区域相邻的两条半边；算法的流程图如图5所示，具体做法是:
[0113]①首次区域生长时，以模型面表中的第一个三角面片作为初始种子面片；后续种子面片通过搜索之前已生长区域的边界，找到还未被划分到任何区域的三角面片作为下一个生长区域的种子面片。将种子面片插入到新子区域对象的面片数组中，以种子面片所属的半边为初始区域边界，并对区域编号。
[0114]②通过生长半边B1寻找待生长面片F，并判断待生长面片F是否符合生长原则；若符合则转至步骤③，不符合则转至步骤④。
[0115]③将待生长面片F插入到面片数组中，将B1从子区边界边数组中删除，并将E1与E2插入到子区边界边数组中，将生长半边转到下条边界半边B2，转至步骤5。
[0116]④将当前边界半边B1标记为区域不再生长的最终边界半边，并将生长半边转到下条半边B2,转至步骤⑤。
[0117]⑤判断区域边界半边是否已全部遍历，并标记为最终边界半边；若是则结束该区域的生长，转至步骤⑥，否则转至步骤②。
[0118]⑥判断三角网格模型中所有面片是否已全部遍历并插入到某一子区域；若是则结束三角网格模型的分割，否则转至步骤①。
[0119]步骤4中所述小面积区域和误判区域的优化合并处理，其中为小面积区域为面积小于设定值的区域；而误判区域的类型一般是凸特征区域或凹特征区域，其特点为形状狭长、不存在极值点且与凹特征区域或凸特征区域相邻。图6所示为误判区域中凸特征区域的截面图，根据局部的微分几何信息，则在区域生长过程中会将图中粗线所代表的区域判定为凹特征区域，而实际上其应归为凸特征区域的一部分。优化合并处理的具体做法是:
[0120]①找出需要被合并的区域记为R1，其中，
[0121]小面积区域的判别方法是:计算区域内所有面片面积的总和A，并判断其是否超过给定面积阈值A。；若A〈A。，则判定其为小面积区域；
[0122]误判区域的判别方法是:先找出凸特征区域或凹特征区域，并找到区域中的最值点Pm，如图6中所示；然后，判断最值点是否为区域的极值点Prait，如图6中所示；若最值点Pm不是极值点Prait，则将其判定为误判区域；极值点Prart的判断方法是:凹特征区域的最小值点或凸特征区域的最大值点是否在区域边界上或接近边界，若是，则该最值点Pm就不是
极值点Prart ；
[0123]②搜索R1的相邻区域，找到最合适的合并区域记为R2,沿着R1的边界搜索其相邻区域并计算合并参数，若对应相邻区域的合并参数最大，则将该区域作为合并区域R2，其中合并参数为:
[0124]小区域的合并参数是边界垂直方向的法曲率半径的平均值；
[0125]误判区域合并参数包括两个方面，其一是相邻边界的长度；其二为特征类型要求，即若R1为凹特征，则R2为凸特征；若R1为凸特征，则R2为凹特征。
[0126]③将需要被合并的区域R1合并到最合适的合并区域R2，分为两部分:
[0127]一是边界的合并，将两区域的相邻边界去除，然后把R1的剩余边界作为R2边界的一部分；
[0128]二是面片的合并，将R1的所有面片归为R2的一部分。
[0129]所述步骤4中的锯齿状边界的子区域边界光顺，以边界分叉点(与其相邻区域超过两个的边界点)为端点，将区域边界分成若干条边界线段；若存在无边界分叉点的区域边界，则将整条边界作为一条边界线段；然后以这些边界线段为处理路径，逐条进行光顺处理，具体做法是:
[0130]①取边界线段中的一个边界顶点P1，搜索P1的一阶邻域面片，并区分属于不同区域的面片。
[0131]②计算WP1为顶点的边界线夹角，如图7中所示的α与β，其中α〈β。当比值β/α>λ时，则判定该局部边界不光顺，其中λ为夹角比值阈值。
[0132]③若局部边界被判定为不光顺，则把与P1相邻的两边界边从区域边界中删除，同时将与P1相对的边E作为新的区域边界，则光顺后的边界如图8所示。
[0133]④将处理对象转至下一边界顶点P2，重复步骤①-③，直到所有边界点都处理完毕为止。
[0134]步骤五、对不同类型加工子区域采用不同的刀具轨迹策略来生成初始刀具轨迹:陡壁子区域采用等高线切法；平缓子区域采用等行距行切法；凹子区域与凸子区域采用等残留高度环切法；鞍形子区域采用等残留高度行切法。
[0135]参考图9、图10，所述对凸子区域与凹子区域采用环切法生成初始刀具轨迹的具体做法是:
[0136]①找出特征子区域的极值点匕；[0137]②按照最大允许残留高度计算行距的方法计算1\周围加工刀触点在某一邻边E1S向行距L，将沿邻边E1向外的单位矢量(记为心方向偏置L/2距离得到刀触点Pi=PL+L/2.θι
[0138]③取θι与E1顺时针方向相邻邻边向外方向的单位矢量(记为e2)之和为角平分线方向(记为e3)，即e3=ei+e2 ;将e3单位化，用②中方法求e3方向的相应刀触点；
[0139]④重复步骤②-③直到匕在所有邻边方向都计算完毕，将上述求得刀触点依次放入刀触点容器中，得到环切初始刀具轨迹。
[0140]图11所示为最大允许残留高度计算行距的方法，具体做法是:
[0141]①计算刀触点Pm处沿行距方向的法曲率K，根据K计算出刀触点Pm与在下一条刀具轨迹的相应刀触点Pn的距离PmPn,则PmPn在投影平面上的投影长度计算如下:
[0142]
【权利要求】
1.一种三角网格模型的分治加工方法，其特征在于:包括以下几个步骤: (1)读取三角网格模型进入内存，去除冗余顶点，建立基于半边结构的拓扑信息； (2)选取R半径范围邻域点后，采用局部二次曲面拟合法来精确计算三角网格模型的微分几何信息，包括主曲率和主方向； (3)根据三角网格模型的形状与微分几何信息将其分成六类自定义的加工特征子区域，分别为陡壁子区域、平缓子区域、局部干涉子区域、凹子区域、凸子区域和鞍形子区域；并以最大主曲率和最小主曲率构成平面坐标系来实现所述加工特征子区域的特征表述； (4)采用以所述加工特征子区域的特征表述为生长原则的区域生长法，将三角网格模型分割成不同特征的加工特征子区域；对分割后产生过分割现象，进行小面积区域和误判区域的优化合并处理，消除面积很小的区域和误判的过渡区域；再对呈现出较多锯齿状边界的所述加工特征子区域进行边界光顺； (5)对不同类型加工子区域采用不同的刀具轨迹策略来生成刀具轨迹: 陡壁子区域采用等高线切法； 平缓子区域采用等行距行切法； 凹子区域与凸子区域采用等残留高度环切法； 鞍形子区域采用等残留高度行切法； (6)将所述初始刀具轨迹用投影偏置扩展的方法生成整个子区域的加工轨迹，并对所述整个子区域的加工轨迹的边界进行处理，从而得到合理完整的子区域加工轨迹。`
2.根据权利要求1所述的一种三角网格模型的分治加工方法，其特征在于:所述的R半径范围邻域点具体是指: 以任意面片F的中心点O为起始点，在三角网格模型表面上与所述中心点O的弧线距离小于等于R的所有顶点的集合。
3.根据权利要求1所述的一种三角网格模型的分治加工方法，其特征在于:所述加工特征子区域的特征表述具体是指: 以最大主曲率K1为横坐标、最小主曲率K2为纵坐标建立主曲率坐标系,所述主曲率坐标系内由于K1 SK2，所有的坐标点都在直线K1=K2上或其右侧；根据不同情况设定Z轴分量阈值V。、平缓曲率阈值K。、凹区域曲率阈值1(,、凸区域曲率阈值Kn等阈值后，所述六类自定义的加工特征子区域的特征表述具体是: ①陡壁子区域:面片法矢的Z轴分量KV。，其中V。为Z轴分量阈值； ②平缓子区域=K1与K2的平方和?+(2X，其中K。为平缓曲率阈值； ③局部干涉子区域=K1X).9KT，其中Kt为球头铣刀顶端的曲率，0.9为避免局部干涉的安全系数； ④凹子区域:(κ；+ κ:)>κ；，且要求Κ2>Κτ，或Κ2〈Κτ且Κ2>-0.21，其中K τ为凹区域曲率阈值； ⑤凸子区域且要求K1In,或K1Xn且!^-。.2K2，其中Kn为凸区域曲率阈值； ⑥鞍形子区域=K1XKK2〈0且不属于上述各区域。
4.根据权利要求1中所述的一种三角网格模型的分治加工方法，其特征在于:所述小面积区域和误判区域的优化合并处理的具体做法是: ①找出需要被合并的区域记为R1，其中， 小面积区域的判别方法是:计算区域内所有面片面积的总和A，并判断其是否超过给定面积阈值A。；若A〈A。，则判定其为小面积区域； 误判区域的判别方法是:先找出凸特征区域或凹特征区域，并找到所述凸特征区域或凹特征区域中的最值点；判断所述最值点是否为所述凸特征区域或凹特征区域的极值点；若不是，则将所述凸特征区域或凹特征区域判定为误判区域；所述极值点的判断方法是:凹特征区域的最小值点或凸特征区域的最大值点是否在区域边界上或接近边界，若是，则该最值点就不是极值点； ②搜索R1的相邻区域:沿着R1的边界搜索其相邻区域并计算合并参数，若对应相邻区域的合并参数最大，则将该区域作为合并区域R2，其中所述合并参数为: 小区域的合并参数是其边界垂直方向的法曲率半径的平均值； 误判区域合并参数包括两个方面，其一是相邻边界的长度；其二为特征类型要求，即若R1为凹特征，则R2为凸特征；若R1为凸特征，则R2为凹特征。 ③将所述需要被合并的区域R1合并到最合适的所述合并区域R2，分为两部分: 一是边界的合并，将两区域的相邻边界去除，然后把所述需要被合并的区域R1的剩余边界作为所述合并区域R2边界的一部分； 二是面片的合并，将所述需要被合并的区域R1的所有面片归为所述合并区域民的一部分。
5.根据权利要求1中所述的一种三角网格模型的分治加工方法，其特征在于:所述对呈现出较多锯齿状边界的所述加工特征子区域进行边界光顺的具体做法是: 以边界分叉点为端点，将区域边界分成若干条边界线段；所述边界分叉点是指与其相邻区域超过两个的边界点；若存在无所述边界分叉点的区域边界，则将整条所述区域边界作为一条边界线段；然后以所述边界线段为处理路径，逐条进行光顺处理，具体步骤是: ①取一条所述边界线段中的一个边界点P1，搜索所述边界点P1的一阶邻域面片，并区分属于不同区域的面片； ②计算以所述边界点P1为顶点的边界线两个夹角α与β，其中α〈β，当比值β/α>λ时，则判定该局部边界不光顺，其中λ为夹角比值阈值； ③若局部边界被判定为不光顺，则把与所述边界点P1相邻的两条所述边界线段从所述区域边界中删除，同时将与所述边界点P1相对的边界线段E作为新的区域边界； ④将处理对象转至下一边界点P2,重复步骤①-③，直到所有边界点都处理完毕为止。
6.根据权利要求1中所述的一种三角网格模型的分治加工方法，其特征在于:所述生成等残留高度环切的初始刀具轨迹的具体做法是: ①找出所述特征子区域的极值点匕; ②按照最大允许残留高度计算行距的方法计算所述极值点匕周围的加工刀触点在某一邻边E1方向行距L，将所述极值点匕沿所述邻边E1向外的单位矢量ei方向偏置L/2距离得到刀触点Pi=PL+L/2.e!； ③取ei与所述邻边E1顺时针方向相邻邻边向外方向的单位矢量e2之和的角平分线方向e3，即e3=ei+e2 ;将e3单位化，用②中方法求e3方向的相应刀触点； ④重复步骤②-③直到所述极值点匕在所有邻边方向都计算完毕，将上述求得刀触点依次放入刀触点容器中，得到环切初始刀具轨迹。
7.根据权利要求1中所述的一种三角网格模型的分治加工方法，其特征在于:所述生成等残留高度行切的初始刀具轨迹的具体做法是: ①求所述特征子区域的包围盒，过所述包围盒的中心点沿刀轴矢量方向做一条射线，求该射线与三角网格曲面的交点P，则得到以所述交点P为初始刀具轨迹的一个刀触点，将所述交点P放入刀触点容器中； ②将所述交点P作为当前处理的刀触点，求出所述交点P的最小主方向T2,过所述交点P以叉积T2XVt为法矢做平面，其中Vt为刀轴矢量。将与所述交点P所在面片求交，得到两交点PpP2 ； ③将所述交点P1插入刀触点容器末端，把所述交点P1作为当前处理刀触点，求所述交点P1的最小主方向T2，过所述交点P1以叉积T2XVt为法矢做平面，用JI与所述交点P1的相邻面片F1求交，将求得的非P1交点P3为下一刀触点； ④将所述交点P3插入刀触点容器末端，把所述交点P3作为当前处理刀触点，重复步骤3，直到求得的刀触点超出所述特征子区域的包围盒； ⑤将所述交点P2插入刀触点容器前端，用步骤③-④中的方法求刀触点，将求得刀触点插入刀触点容器前端，直到求得的刀触点超出所述特征子区域的包围盒；最终得到行切初始刀具轨迹。
8.根据权利要求1中所`述的一种三角网格模型的分治加工方法，其特征在于:所述对整个子区域的加工轨迹的边界进行处理的具体做法是: ①当所述环切初始刀具轨迹的扩展达到边界时，先将完整的环切刀具轨迹复制备份到一个清空的交点数组中；然后处理所述环切刀具轨迹，删除其中位于子区域外的刀触点，保留处于区域内的刀触点作为有效的环切刀具轨迹；后续的环切刀具轨迹则基于备份的所述完整的所述环切刀具轨迹来扩展生成； ②对于所述行切刀具轨迹，采用偏置扩展法，由上一条备份的所述行切刀具轨迹获得一条首末点超出所述特征子区域包围盒的未处理的行切刀具轨迹；将所述未处理的行切刀具轨迹复制备份到一个清空的交点数组中；将所述未处理的行切刀具轨迹中位于所述特征子区域边界以外的顶点删除，使每条行切刀具轨迹的首末点处于所述特征子区域边界。
【文档编号】G05B19/19GK103885385SQ201410081684
【公开日】2014年6月25日 申请日期:2014年3月7日 优先权日:2014年3月7日
【发明者】黄常标, 江开勇, 林俊峰, 刘斌, 林俊义, 路平 申请人:华侨大学