CN115358009A - 交通工具外壳造型中nurbs曲线曲面光顺生成方法 - Google Patents

交通工具外壳造型中nurbs曲线曲面光顺生成方法 Download PDF

Info

Publication number
CN115358009A
CN115358009A CN202211045861.2A CN202211045861A CN115358009A CN 115358009 A CN115358009 A CN 115358009A CN 202211045861 A CN202211045861 A CN 202211045861A CN 115358009 A CN115358009 A CN 115358009A
Authority
CN
China
Prior art keywords
nurbs
control point
jth
fairing
vector
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202211045861.2A
Other languages
English (en)
Inventor
蔺宏伟
蒋旖旎
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202211045861.2A priority Critical patent/CN115358009A/zh
Publication of CN115358009A publication Critical patent/CN115358009A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00—Road transport of goods or passengers
    • Y02T10/10—Internal combustion engine [ICE] based vehicles
    • Y02T10/40—Engine management systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Data Mining & Analysis (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Computer Hardware Design (AREA)
  • Operations Research (AREA)
  • Evolutionary Computation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Computing Systems (AREA)
  • Image Generation (AREA)

Abstract

本发明公开了一种交通工具外壳造型中NURBS曲线曲面光顺生成方法。已知拟合参数模型的模型方程，对交通工具外壳造型进行逆向工程的数据采集获得数据点集，参数化后获得模型方程的节点向量，并选取初始的控制点，由控制点和模型方程的节点向量生成拟合参数模型；根据交通工具外壳造型光顺效果的实际要求对拟合参数模型和控制点进行迭代优化，获得最优的拟合参数模型。本发明通过调整控制点对应的光顺权重实现了对曲线曲面的精准光顺，提高了光顺操作的灵活性与交互性。同时，几何迭代过程保证了用户在每步迭代后可以全局或局部地修改控制点的权值、节点向量或数据参数化等。

Description

交通工具外壳造型中NURBS曲线曲面光顺生成方法
技术领域
本发明涉及一般的三维模型或者二维模型的数据拟合和光顺处理方法，尤其涉及NURBS曲线曲面光顺的渐进迭代逼近方法。
背景技术
给定离散的数据点集，利用参数曲线曲面进行拟合，得到光顺的曲线曲面是许多领域的基本研究问题。在几何设计领域中，通常需要生成光滑的机翼翼型、船体造型和汽车引擎盖等，因为这些几何特性会极大地影响它们的物理性能。近年来，许多学者深入研究了光顺曲线在机器人路径规划中的应用，这使得如何灵活和高效地进行曲线曲面光顺成为了一个热点问题。
传统的光顺方法大致分为两类：全局光顺方法(参见H.Hagen，G.-P.Bonneau，Variational design of smooth rational bézier curves，Computer aided geometricdesign 8，5(1991)，393–399.)和局部光顺方法(参见G.Farin，G.Rein，N.Sapidis，A.Worsey，Fairing cubic b-spline curves，Computer Aided Geometric Design 4，1(1987)，91–103.)。这两类方法的主要特点是能高效的改变模型的能量值，从而获得光滑的曲线曲面。全局光顺法通常全局地调整能量值，而局部光顺法仅对局部区域进行调整。然而这两类方法，通常缺少足够的参数来微调和提升曲线曲面的光顺性。这使得曲线曲面的光顺过程可控性较差，灵活度较低。因此，精细调整光顺效果的能力是光顺曲线曲面的关键因素。
坏点移除法是一种交互式的光顺算法，通过对不满足光顺要求的点进行修改来提高曲线的光顺性(参见J.A.P.Kjellander，Smoothing of cubic parametric splines，Computer aided design 15，3(1983)，175–179.)。随后，Farin等人首次提出了先节点删除再重新插入的方法以提高修改点的效率。然而，当遇到大量不满足光顺要求的点时，这种局部光顺方法只能进行逐点处理，这使得光顺过程变得十分低效。
能量最小化法的基本思想是找到使得曲线曲面能量模型最小的参数模型，从而得到最光滑的曲线曲面。通常，曲线曲面对应的能量值越小，该曲线曲面的光顺度越高、物理性能越好。其中常用的能量包括应变能(参见H.Hagen，G.-P.Bonneau，Variational designof smooth rational bézier curves，Computer aided geometric design 8，5(1991)，393–399.)、张力能(参见H.Pottmann，Smooth curves under tension，Computer aideddesign 22，4(1990)，241–245.)、和曲率半径的变化(参见T.C.Rando，Automatic fairnessin computer-aided geometric design，Ph.D.thesis，aAI9102032(1990).)等。能量法虽然可以对参数模型进行全局优化，但在工业设计和制造中采集得到的数据点规模通常较大。能量法的计算耗时并且不能很好地控制拟合效果。
Vassilev基于上述两种方法，提出了一种自动光顺算法(参见T.I.Vassilev，Fairinterpolation and approximation of b-splines by energy minimization andpoints insertion，Computer aided design 28，9(1996)，753–760.)。与之前的方法不同的是，该方法只在必要的时候才会插入新的点。从而大大减少了问题规模，提高了光顺效率。受这项工作的启发，许多研究者尝试提出更高效地识别坏点和更新能量最小化导出的控制点的方法。
发明内容
本发明的目的在于针对现有技术的不足，提出一种交通工具外壳造型中NURBS曲线曲面光顺生成方法。
本方法利用了对应每个控制点的光顺权重，生成了可全局或局部调整的光顺曲线(曲面)；发明了曲线曲面上差向量在控制点上的聚合方法，即归一化系数，保证了几何迭代光顺算法的收敛性；利用了几何迭代过程中参数便于灵活修改的特性，适用于节点插入，实现更好的光顺权重或数据参数化调整。
本发明采用的技术方案如下：
1)已知拟合参数模型的模型方程，针对交通工具外壳造型进行逆向工程的数据采集获得数据点集，数据点集通常为通过扫描仪扫描获得的点云。对数据点集参数化获得各个数据点的参数，根据各个数据点的参数处理获得模型方程的节点向量，并选取初始的控制点，由初始的控制点和模型方程的节点向量生成初始的拟合参数模型；
2)根据交通工具外壳造型光顺效果的实际要求对拟合参数模型和控制点进行双重迭代优化，获得最优的拟合参数模型，可再用拟合参数模型输入计算机中进行3D打印或者制造。
所述的拟合参数模型为NURBS曲线或NURBS曲面。
所述步骤2)具体为：
2.1)计算拟合参数模型与各个数据点之间的差向量di，i表示数据点的序数，将参数值在一个基函数的支撑域的数据点对应的差向量di分为一组，一组对应设置一个控制点；然后对一组中的所有差向量di进行加权求和，得到该组对应的控制点的拟合向量δj；一组数据点对应设置一个控制点。
拟合参数模型的每个控制点对应一个基函数，基函数的非零域作为支撑域。
2.2)设置光顺泛函
Figure BDA0003822363790000031
根据光顺泛函
Figure BDA0003822363790000032
获得每个控制点的光顺向量ηj；
2.3)根据交通工具外壳造型光顺效果的实际要求设置控制点的光顺权重ωj，并计算控制点的差向量Bj：
Figure BDA0003822363790000033
其中，Bj表示第j个控制点的差向量，上标k为迭代次数；
Figure BDA0003822363790000034
表示第k次迭代下获得的第j个控制点的光顺向量ηj，
Figure BDA0003822363790000035
表示第k次迭代下获得的第j个控制点的拟合向量δj；
2.4)根据步骤2.3)得到的结果，按照以下公式更新控制点的坐标：
Figure BDA0003822363790000036
式中：μj为第j个控制点的归一化系数，以保证迭代收敛；
Figure BDA0003822363790000037
表示第k次迭代的第j个控制点的坐标；
2.5)将更新后的控制点的坐标代入拟合参数模型的模型方程，生成新的拟合参数模型，再计算当前第k次迭代的相对迭代误差
Figure BDA0003822363790000038
Figure BDA0003822363790000039
A＝(I-Ω)NTN+ΩDr
Figure BDA00038223637900000310
式中：k为迭代次数，I为单位矩阵，Ω表示一个以光顺权重ωj为对角元素的对角矩阵，N表示配置矩阵，Q为由数据点构成的矩阵，Dr是由步骤3)中的光顺泛函构成的矩阵，r是对拟合参数模型的模型方程中的基函数求导的阶数，P[k]为第k次迭代后获得的控制点集合，j表示矩阵的第j行；A表示系数矩阵，(j,:)表示取矩阵第j行所有列的元素，||||2表示范数的平方，T表示矩阵的转置；
将当前第k次迭代的相对迭代误差
Figure BDA00038223637900000311
与前一次第k-1次迭代的相对迭代误差
Figure BDA00038223637900000312
进行比较判断：
若之间的差值小于预设的差值阈值εiter时，则迭代终止，进行下一步骤；
若之间的差值不小于预设的差值阈值εiter时，则以新的拟合参数模型回到步骤2.1)进行下一次迭代；
2.6)根据上一步骤中迭代终止得到的最新的拟合参数模型，计算相对拟合误差
Figure BDA0003822363790000041
Figure BDA0003822363790000042
其中，Qi表示第i个数据点，i表示数据点的序号，m表示数据点的总数，ti表示第i个数据点的参数，P[k](ti)表示第k次迭代后第i个数据点对应的参数值在拟合参数模型中的值；
将相对拟合误差
Figure BDA0003822363790000043
和预设的目标拟合精度εfit进行比较判断：
若相对拟合误差
Figure BDA0003822363790000044
未达到目标拟合精度εfit时，则调整拟合参数模型的模型方程的节点向量、调整控制点对应的光顺权重ωj或者调整数据点集的参数化；
a)调整拟合参数模型的模型方程的节点向量：对拟合误差较大的区域，采用节点插入算法，按需插入节点，得到新的控制点。
a)b)调整控制点对应的光顺权重ωj：修改每个控制点对应的光顺权重。ωj的数值越大，光顺效果越好；反之，拟合效果越好。
b)c)调整数据点集的参数化：依照数据点的几何特征，重新参数化。
回到步骤2.1)进行下一次迭代；
若相对拟合误差
Figure BDA0003822363790000045
达到目标拟合精度εfit时，则最新生成的拟合参数模型作为最终结果
由此不断迭代直至生成的拟合参数模型满足所述的目标拟合精度εfit和光顺要求。
所述的拟合参数模型为NURBS曲线，所述步骤1)中的NURBS曲线的模型方程为：
Figure BDA0003822363790000046
式中：
Figure BDA0003822363790000047
为随机选取的初始的第j个控制点，t为NURBS(非均匀有理B样条)曲线的参数，Nj(t)表示第j个基函数，n表示控制点的总数；P[0](t)表示初始的NURBS曲线，其中上标[0]表示初始拟合曲线，t表示曲线参数。
t值不加下标i时，表示一个连续的曲线参数值；加上下标i表示第i个数据点对应的曲线上的第i个数据点ti。t和ti本质表示的都是曲线上的参数。一个是连续值，一个是离散值。
所述的拟合参数模型为NURBS曲线，所述步骤2.1)具体为：
先按照以下公式计算NURBS曲线P[k](t)与每个数据点Qi之间的差向量
Figure BDA0003822363790000048
Figure BDA0003822363790000051
其中，P[k](ti)表示第k次迭代后第i个数据点的参数值对应的拟合参数模型的值，
Figure BDA0003822363790000052
表示第k次迭代后第i个数据点对应的差向量；
然后将参数值在一个基函数的支撑域的数据点对应的差向量
Figure BDA0003822363790000053
分为一组，对一组中所有数据点的差向量取加权和，作为这一组数据点对应的控制点的拟合向量δj [k]，其中加权的权值为第j个基函数Ni(t)在NURBS曲线的参数t处的值；
2.2)设置光顺泛函
Figure BDA0003822363790000054
根据光顺泛函
Figure BDA0003822363790000055
获得每个控制点的光顺向量ηj。
所述的拟合参数模型为NURBS曲线，所述步骤2.2)中，具体按照以下公式获得第j个控制点的光顺向量：
Figure BDA0003822363790000056
Figure BDA0003822363790000057
式中：t为NURBS曲线的参数，范围是[t1,tm]，tm表示第m个参数值，
Figure BDA0003822363790000058
表示利用第j个基函数的r阶导构造的关于第l个基函数r阶导的光顺泛函
Figure BDA0003822363790000059
Nr,j(t)表示对NURBS曲线的第j个基函数求r阶导的结果，Pl [k]表示为第k次迭代下的第l个控制点，Nr,l(t)表示对NURBS曲线的第l个基函数求r阶导的结果，n表示基函数个数，
Figure BDA00038223637900000510
表示第k次迭代后第j个控制点的光顺向量。
所述的拟合参数模型为NURBS曲面，所述步骤1)中，曲面相关的数据点集和对应的参数值按曲面上的两个维度方向的顺序排列，将基函数和控制点也按字典序排列，表示为：
Figure BDA00038223637900000511
Figure BDA00038223637900000512
式中：u,v为NURBS曲面上的两个维度方向参数，n1,n2分别为NURBS曲面上u,v方向的控制点个数；
Figure BDA00038223637900000513
分别表示u方向上第n1个基函数，v方向上第n2个基函数，
Figure BDA00038223637900000514
表示控制网格上第n1行第n2列所对应的控制点坐标，
Figure BDA00038223637900000515
表示将所有控制点按上述规则排序后的第n1*n2个控制点，
Figure BDA00038223637900000516
表示将所有基函数按上述规则排序后的第n1*n2个基函数；
NURBS曲面的模型方程为：
Figure BDA00038223637900000517
式中：
Figure BDA00038223637900000518
为随机选取的初始的第j个控制点，Nj(u,v)表示第j个基函数N(u,v)，(u,v)为基函数N(u,v)的参数值；P[0](u,v)表示初始的NURBS曲面。
按照上述模型方程生成NURBS曲面。
所述的拟合参数模型为NURBS曲面，所述步骤2.1)具体为：
先按照以下公式计算NURBS曲面P[k](u,v)与每个数据点Qi之间的差向量
Figure BDA0003822363790000061
Figure BDA0003822363790000062
其中，
Figure BDA0003822363790000063
表示第k次迭代后第i个数据点对应的差向量，P[k](ui,vi)表示第k次迭代后第i个数据点的参数值对应的拟合参数模型的值；
然后将参数值在一个基函数的支撑域的数据点对应的差向量
Figure BDA0003822363790000064
分为一组，对一组中所有数据点的差向量取加权和，作为该组数据点对应的控制点的拟合向量δj [k]，其中加权的权值为第j个基函数Ni(u,v)在NURBS曲面的参数(u,v)处的值。
所述的拟合参数模型为NURBS曲面，所述步骤2.2)中，选择基于一阶偏导的光顺向量、基于二阶偏导的光顺向量、基于三阶偏导的光顺向量中的一个获得控制点的光顺向量，分别为：
A)第j个控制点的基于一阶偏导的光顺向量计算为：
Figure BDA0003822363790000065
式中：u,v为NURBS曲面P[k](u,v)上的两个维度方向参数，n1,n2分别为NURBS曲面沿两个方向u,v上的控制点个数；Pl [k]表示为第k次迭代下的第l个控制点，l表示控制点的序数，
Figure BDA0003822363790000066
表示第k次迭代后的第j个控制点对应的光顺向量；Nu,l(u,v)、Nv,l(u,v)分别表示NURBS曲面的第j个基函数Nj(u,v)在两个方向u,v上的求一阶偏导的结果，
Figure BDA0003822363790000067
表示u方向第j个一阶偏导光顺泛函，
Figure BDA0003822363790000068
表示v方向第j个一阶偏导光顺泛函；
上述两个光顺泛函计算为：
Figure BDA0003822363790000069
Figure BDA00038223637900000610
其中，um表示u方向第m个参数值，vm表示v方向第m个参数值，m表示数据点的序数；
B)第j个控制点的基于二阶偏导的光顺向量计算为：
Figure BDA0003822363790000071
式中：Nuu,j(u,v)、Nuv,j(u,v)和Nvv,j(u,v)分别表示表示对第j个NURBS基函数Nj(u,v)求二阶偏导的结果，
Figure BDA0003822363790000072
表示第j个第一类二阶偏导光顺泛函，
Figure BDA0003822363790000073
表示第j个第二类二阶偏导光顺泛函，
Figure BDA0003822363790000074
表示第j个第三类二阶偏导光顺泛函；
上述三个光顺泛函计算为：
Figure BDA0003822363790000075
Figure BDA0003822363790000076
Figure BDA0003822363790000077
C)第j个控制点的基于三阶偏导的光顺向量计算为：
Figure BDA0003822363790000078
式中：Nuuu,j(u,v)、Nuuv,j(u,v)、Nuvv,j(u,v)和Nvvv,j(u,v)分别表示对第j个NURBS基函数Nj(u,v)求三阶偏导的结果，
Figure BDA0003822363790000079
表示第j个第一类三阶偏导光顺泛函，
Figure BDA00038223637900000710
表示第j个第二类三阶偏导光顺泛函，
Figure BDA00038223637900000711
表示第j个第三类三阶偏导光顺泛函，
Figure BDA00038223637900000712
表示第j个第四类三阶偏导光顺泛函；
上述四个光顺泛函计算为：
Figure BDA00038223637900000713
Figure BDA00038223637900000714
Figure BDA00038223637900000715
Figure BDA00038223637900000716
具体实施中，在步骤3.3)中控制点的光顺权重ωj对于每个控制点是可以不相同的。特别地，当光顺权重ωj等于一个常数值时，为能量最小化光顺法。
步骤6)中迭代终止后，判断当前光顺效果和拟合误差是否达到阈值，若未满足，则全局或局部地插入新节点、调整光顺权重或数据参数化。然后利用前一轮迭代的曲线(曲面)结果作为下一轮迭代的初始曲线(曲面)，继续迭代直至满足要求。这样可以得到光顺性高并且拟合误差也满足阈值的NURBS曲线(曲面)。
所述的交通工具外壳造型为机翼翼型、船体造型和汽车引擎盖等光滑外壳件，但不限于此。
本发明通过调整控制点对应的光顺权重实现了对曲线曲面的精准光顺，提高了光顺操作的灵活性与交互性。同时，几何迭代过程保证了用户在每步迭代后可以全局或局部地修改控制点的权值、节点向量或数据参数化等。
本发明与背景技术相比较，其优点在于：
1、本发明方法采用丰富的参数，即每个控制点对应一个光顺权重，来精细地调整和提高曲线曲面的光顺性。对于光顺权重的调节可以全局或局部调控光顺的效果，提高了模型光顺过程的可控性。
2、本发明方法是一种高效且收敛的迭代算法，大大减少了计算的内存消耗。当所有的光顺权重相等时，本方法等效于能量最小化光顺法。
3、本发明方法是一种灵活的、可交互的光顺算法，用户可以在每轮迭代后改变控制点的权重、节点向量或数据参数化等，这大幅度提高了光顺过程的灵活度。
附图说明
图1是本发明方法的流程图。
图2是实例1中的机翼截面数据点集，初始的控制点集和初始拟合曲线。
图3是实例1中拟合曲线和拟合曲线的曲率杆图。
图4是实例1中光顺曲线和光顺曲线的曲率杆图。
图5是实例2中的三维牙齿数据点集。
图6是实例2中拟合曲面和拟合曲面的斑马纹图。
图7是实例2中光顺曲面和光顺曲面的斑马纹图。
具体实施方式
下面结合附图和具体实施对本发明作进一步说明。
如图1所示，本发明的实施例及其实施过程如下：
实施例1：
本实施例针对飞机机翼，通过三维扫描采集获得飞机机翼的截面数据点集，如图2所示。
1)数据点的导入与参数化
输入采集得到的数据点集Qi,i＝1,2,...,m，采用弦长参数化方法，给数据点分配对应参数值t1≤t2≤…≤tm。
2)初始样条曲线的构造
从数据点中随机选取初始控制点Pj [0],j＝1,2,...,n，获得初始拟合的NURBS曲线
Figure BDA0003822363790000091
其中，Nj(t)为第j个基函数。如图2所示。
3)拟合向量，光顺向量和差向量的计算
首先，计算NURBS曲线P[k](t)与每个数据点Qi之间的差向量
Figure BDA0003822363790000092
Figure BDA0003822363790000093
其次将参数值在一个基函数的支撑域的数据点对应的差向量
Figure BDA0003822363790000094
分为一组，对一组中所有数据点的差向量取加权和，作为这一组数据点对应的控制点的拟合向量δj [k]，其中加权的权重为第j个基函数Nj(t)在NURBS曲线的参数t处的值。
然后，设置光顺泛函
Figure BDA00038223637900000913
这里选取r＝2。根据光顺泛函
Figure BDA00038223637900000914
获得每个控制点对应的光顺向量ηj [k]：
Figure BDA0003822363790000095
Figure BDA0003822363790000096
式中：t为NURBS曲线的参数，范围是[t1,tm]，tm表示第m个参数值，
Figure BDA0003822363790000097
表示一个光顺泛函，N2,j(t)表示对NURBS曲线的第j个基函数求2阶导的结果，Pl [k]表示为第k次迭代下的第l个控制点，N1,l(t)表示对NURBS曲线的第l个基函数求2阶导的结果，n表示基函数个数，
Figure BDA0003822363790000098
表示第k次迭代后第j个控制点的光顺向量。
4)设置控制点的光顺权重ωj并计算控制点的差向量Bj：
Figure BDA0003822363790000099
5)控制点的更新
根据差向量Bj调整控制点：
Figure BDA00038223637900000910
其中，μj是归一化系数，保证迭代收敛。
Figure BDA00038223637900000911
和
Figure BDA00038223637900000912
分别表示第k次和第k+1次迭代后生成的第j个控制点的坐标。
6)节点，光顺权重或数据参数化的调整
在整个迭代过程达到稳定停止之后，如果当前光顺曲线没有达到用户的光顺要求和拟合精度εfit时，我们可以全局或局部地插入一些节点，调整光顺权重ωj或数据参数化，进一步提高整个网格的自由度，减小拟合误差或改善光顺效果。
A)对拟合误差较大的区域，采用节点插入算法，按需插入节点，生成新的控制点。
B)修改每个控制点对应的光顺权重ωj。ωj的数值越大，光顺效果越好；ωj的数值越小，拟合效果越好。
C)依照数据点的几何特征，重新参数化。
在完成上述修改后，令前一轮迭代的结果作为初始拟合曲线，继续迭代，直至满足用户需求，输出光顺曲线。如图4所示。
对图3中三个方框内对应的光顺权重依次进行调整，得到如图4所示的光顺曲线。曲率杆图中的曲率杆长度反映了曲线上各点曲率值的大小；曲率杆变化的程度反映了曲线曲率的变化程度。如图4所示，光顺曲线的曲率值更小，曲率变化更平滑，即光顺性更好。
实施例2：
本实施例针对牙齿产品，通过三维扫描采集获得牙齿的数据点集，如图5所示。
1)数据点的导入与参数化
输入采集得到的数据点集Qij,i＝1,2,...,m1,j＝1,2,...,m2，采用弦长参数化方法，给数据点分配对应参数值参数值{ui,vj},i＝1,2,...,m1,j＝1,2,...,m2，满足
Figure BDA0003822363790000101
和
Figure BDA0003822363790000102
将数据点Qi,i＝1,2,...,m按曲面上的两个维度方向的顺序排列，对应的参数值为(ui,vi),i＝1,2,...,m。
2)初始样条曲面的构造
从数据点中随机选取初始控制点并将控制点也按字典序排列，表示为：
Figure BDA0003822363790000103
Figure BDA0003822363790000104
NURBS曲面的模型方程为：
Figure BDA0003822363790000105
3)拟合向量，光顺向量和差向量的计算
首先，计算NURBS曲面P[k](u,v)与每个数据点Qi之间的差向量
Figure BDA0003822363790000106
Figure BDA0003822363790000107
其次将参数值在一个基函数的支撑域的数据点对应的差向量
Figure BDA0003822363790000111
分为一组，对一组中所有数据点的差向量取加权和，作为这一组数据点对应的控制点的拟合向量δj [k]，其中加权的权重为第j个基函数Nj(u,v)在NURBS曲面的参数(u,v)处的值。
然后，设置光顺泛函
Figure BDA0003822363790000112
根据光顺泛函
Figure BDA0003822363790000113
获得每个控制点对应的光顺向量ηj [k]，这里采用二阶偏导进行光顺向量的计算。
第j个控制点的基于二阶偏导的光顺向量为：
Figure BDA0003822363790000114
式中：Nuu,j(u,v)、Nuv,j(u,v)和Nvv,j(u,v)分别表示表示对第j个NURBS基函数Nj(u,v)求二阶偏导的结果，
Figure BDA0003822363790000115
表示第j个第一类二阶偏导光顺泛函，
Figure BDA0003822363790000116
表示第j个第二类二阶偏导光顺泛函，
Figure BDA0003822363790000117
表示第j个第三类二阶偏导光顺泛函；
上述三个光顺泛函计算为：
Figure BDA0003822363790000118
Figure BDA0003822363790000119
Figure BDA00038223637900001110
4)设置控制点的光顺权重ωj并计算控制点的差向量Bj：
Figure BDA00038223637900001111
5)控制点的更新
根据差向量Bj调整控制点：
Figure BDA00038223637900001112
其中，μj是归一化系数，保证迭代收敛。
Figure BDA00038223637900001113
和
Figure BDA00038223637900001114
分别表示第k次和第k+1次迭代后生成的第j个控制点的坐标。
6)节点，光顺权重或数据参数化的调整
在整个迭代过程达到稳定停止之后，如果当前光顺曲面没有达到用户的光顺要求和拟合精度εfit时，我们可以全局或局部地插入一些节点，调整光顺权重ωj或数据参数化，进一步提高整个网格的自由度，减小拟合误差或改善光顺效果。
A)对拟合误差较大的区域，采用节点插入算法，按需插入节点，生成新的控制点。
B)修改每个控制点对应的光顺权重ωj。ωj的数值越大，光顺效果越好；ωj的数值越小，拟合效果越好。
C)依照数据点的几何特征，重新参数化。
在完成上述修改后，令前一轮迭代的结果作为初始拟合曲面，继续迭代，直至满足用户需求，输出光顺曲面。
光顺曲面和光顺曲面的斑马纹图如图7所示。图6是牙齿模型的拟合曲面和拟合曲面的斑马纹图。斑马纹的疏密、弯曲能反映曲面的光顺程度，即斑马纹越宽、转折越少、越平缓，曲面的光顺性越好。对比图6与图7，可以说明本发明方法生成的曲面光顺性更好。

Claims (9)

1.一种交通工具外壳造型中NURBS曲线曲面光顺生成方法，其特征在于：
方法包括以下几个步骤：
1)已知拟合参数模型的模型方程，针对交通工具外壳造型进行逆向工程的数据采集获得数据点集，对数据点集参数化获得各个数据点的参数，根据各个数据点的参数处理获得模型方程的节点向量，并选取初始的控制点，由初始的控制点和模型方程的节点向量生成初始的拟合参数模型；
2)根据交通工具外壳造型光顺效果的实际要求对拟合参数模型和控制点进行迭代优化，获得最优的拟合参数模型。
2.根据权利要求1所述的一种交通工具外壳造型中NURBS曲线曲面光顺生成方法，其特征在于：所述的拟合参数模型为NURBS曲线或NURBS曲面。
3.根据权利要求1所述的一种交通工具外壳造型中NURBS曲线曲面光顺生成方法，其特征在于：所述步骤2)具体为：
2.1)计算拟合参数模型与各个数据点之间的差向量di，i表示数据点的序数，将参数值在一个基函数的支撑域的数据点对应的差向量di分为一组，一组对应设置一个控制点；然后对一组中的所有差向量di进行加权求和，得到该组对应的控制点的拟合向量δj；
2.2)设置光顺泛函
Figure FDA0003822363780000011
根据光顺泛函
Figure FDA0003822363780000012
获得每个控制点的光顺向量ηj；
2.3)根据交通工具外壳造型光顺效果的实际要求设置控制点的光顺权重ωj，并计算控制点的差向量Bj：
Figure FDA0003822363780000013
其中，Bj表示第j个控制点的差向量，上标k为迭代次数；
Figure FDA0003822363780000014
表示第k次迭代下获得的第j个控制点的光顺向量ηj，
Figure FDA0003822363780000015
表示第k次迭代下获得的第j个控制点的拟合向量δj；
2.4)根据步骤2.3)得到的结果，按照以下公式更新控制点的坐标：
Figure FDA0003822363780000016
式中：μj为第j个控制点的归一化系数；
Figure FDA0003822363780000017
表示第k次迭代的第j个控制点的坐标；
2.5)将更新后的控制点的坐标代入拟合参数模型的模型方程，生成新的拟合参数模型，再计算当前第k次迭代的相对迭代误差
Figure FDA0003822363780000018
Figure FDA0003822363780000021
A＝(I-Ω)NTN+ΩDr
Figure FDA0003822363780000022
式中：k为迭代次数，I为单位矩阵，Ω表示一个以光顺权重ωj为对角元素的对角矩阵，N表示配置矩阵，Q为由数据点构成的矩阵，Dr是由步骤3)中的光顺泛函构成的矩阵，r是对拟合参数模型的模型方程中的基函数求导的阶数，P[k]为第k次迭代后获得的控制点集合，j表示矩阵的第j行；A表示系数矩阵，(j,:)表示取矩阵第j行所有列的元素，|| ||2表示范数的平方，T表示矩阵的转置；
将当前第k次迭代的相对迭代误差
Figure FDA0003822363780000023
与前一次第k-1次迭代的相对迭代误差
Figure FDA0003822363780000024
进行比较判断：
若之间的差值小于预设的差值阈值εiter时，则迭代终止，进行下一步骤；
若之间的差值不小于预设的差值阈值εiter时，则以新的拟合参数模型回到步骤2.1)进行下一次迭代；
2.6)根据上一步骤中迭代终止得到的拟合参数模型，计算相对拟合误差
Figure FDA0003822363780000025
Figure FDA0003822363780000026
其中，Qi表示第i个数据点，i表示数据点的序号，m表示数据点的总数，ti表示第i个数据点的参数，P[k](ti)表示第k次迭代后第i个数据点对应的参数值在拟合参数模型中的值；
将相对拟合误差
Figure FDA0003822363780000027
和预设的目标拟合精度εfit进行比较判断：
若相对拟合误差
Figure FDA0003822363780000028
未达到目标拟合精度εfit时，则调整拟合参数模型的模型方程的节点向量、调整控制点对应的光顺权重ωj或者调整数据点集的参数化；
回到步骤2.1)进行下一次迭代；
若相对拟合误差
Figure FDA0003822363780000029
达到目标拟合精度εfit时，则最新生成的拟合参数模型作为最终结果。
4.如权利要求3所述的一种交通工具外壳造型中NURBS曲线曲面光顺生成方法，其特征在于：所述的拟合参数模型为NURBS曲线，所述步骤1)中的NURBS曲线的模型方程为：
Figure FDA0003822363780000031
式中：
Figure FDA0003822363780000032
为初始的第j个控制点，t为NURBS曲线的参数，Nj(t)表示第j个基函数，n表示控制点的总数；P[0](t)表示初始的NURBS曲线，其中上标[0]表示初始拟合曲线，t表示曲线参数。
5.如权利要求3所述的一种交通工具外壳造型中NURBS曲线曲面光顺生成方法，其特征在于：所述的拟合参数模型为NURBS曲线，所述步骤2.1)具体为：
先按照以下公式计算NURBS曲线P[k](t)与每个数据点Qi之间的差向量di [k]：
Figure FDA0003822363780000033
其中，P[k](ti)表示第k次迭代后第i个数据点的参数值对应的拟合参数模型的值，
Figure FDA0003822363780000034
表示第k次迭代后第i个数据点对应的差向量；
然后将参数值在一个基函数的支撑域的数据点对应的差向量
Figure FDA0003822363780000035
分为一组，对一组中所有数据点的差向量取加权和，作为这一组数据点对应的控制点的拟合向量
Figure FDA0003822363780000036
其中加权的权值为第j个基函数Ni(t)在NURBS曲线的参数t处的值；
2.2)设置光顺泛函
Figure FDA0003822363780000037
根据光顺泛函
Figure FDA0003822363780000038
获得每个控制点的光顺向量ηj。
6.如权利要求3所述的一种交通工具外壳造型中NURBS曲线曲面光顺生成方法，其特征在于：所述的拟合参数模型为NURBS曲线，所述步骤2.2)中，具体按照以下公式获得第j个控制点的光顺向量：
Figure FDA0003822363780000039
Figure FDA00038223637800000310
式中：t为NURBS曲线的参数，范围是[t1,tm]，tm表示第m个参数值，
Figure FDA00038223637800000311
表示利用第j个基函数的r阶导构造的关于第l个基函数r阶导的光顺泛函
Figure FDA00038223637800000312
Nr,j(t)表示对NURBS曲线的第j个基函数求r阶导的结果，
Figure FDA00038223637800000313
表示为第k次迭代下的第l个控制点，Nr,l(t)表示对NURBS曲线的第l个基函数求r阶导的结果，n表示基函数个数，
Figure FDA00038223637800000314
表示第k次迭代后第j个控制点的光顺向量。
7.如权利要求3所述的一种交通工具外壳造型中NURBS曲线曲面光顺生成方法，其特征在于：所述的拟合参数模型为NURBS曲面，所述步骤1)中，曲面相关的数据点集和对应的参数值按曲面上的两个维度方向的顺序排列，将基函数和控制点也按字典序排列，表示为：
Figure FDA0003822363780000041
Figure FDA0003822363780000042
式中：u,v为NURBS曲面上的两个维度方向参数，n1,n2分别为NURBS曲面上u,v方向的控制点个数；
Figure FDA0003822363780000043
分别表示u方向上第n1个基函数，v方向上第n2个基函数，
Figure FDA0003822363780000044
表示控制网格上第n1行第n2列所对应的控制点坐标，
Figure FDA0003822363780000045
表示将所有控制点按上述规则排序后的第n1*n2个控制点，
Figure FDA0003822363780000046
表示将所有基函数按上述规则排序后的第n1*n2个基函数；
NURBS曲面的模型方程为：
Figure FDA0003822363780000047
式中：
Figure FDA0003822363780000048
为初始的第j个控制点，Nj(u,v)表示第j个基函数N(u,v)，(u,v)为基函数N(u,v)的参数值；P[0](u,v)表示初始的NURBS曲面。
8.如权利要求3所述的一种交通工具外壳造型中NURBS曲线曲面光顺生成方法，其特征在于：所述的拟合参数模型为NURBS曲面，所述步骤2.1)具体为：
先按照以下公式计算NURBS曲面P[k](u,v)与每个数据点Qi之间的差向量
Figure FDA0003822363780000049
Figure FDA00038223637800000410
其中，
Figure FDA00038223637800000411
表示第k次迭代后第i个数据点对应的差向量，P[k](ui,vi)表示第k次迭代后第i个数据点的参数值对应的拟合参数模型的值；
然后将参数值在一个基函数的支撑域的数据点对应的差向量
Figure FDA00038223637800000412
分为一组，对一组中所有数据点的差向量取加权和，作为该组数据点对应的控制点的拟合向量
Figure FDA00038223637800000413
其中加权的权值为第j个基函数Ni(u,v)在NURBS曲面的参数(u,v)处的值。
9.如权利要求3所述的一种交通工具外壳造型中NURBS曲线曲面光顺生成方法，其特征在于：所述的拟合参数模型为NURBS曲面，所述步骤2.2)中，选择基于一阶偏导的光顺向量、基于二阶偏导的光顺向量、基于三阶偏导的光顺向量中的一个获得控制点的光顺向量，分别为：
A)第j个控制点的基于一阶偏导的光顺向量计算为：
Figure FDA00038223637800000414
式中：u,v为NURBS曲面P[k](u,v)上的两个维度方向参数，n1,n2分别为NURBS曲面沿两个方向u,v上的控制点个数；Pl [k]表示为第k次迭代下的第l个控制点，l表示控制点的序数，
Figure FDA0003822363780000051
表示第k次迭代后的第j个控制点对应的光顺向量；Nu,l(u,v)、Nv,l(u,v)分别表示NURBS曲面的第j个基函数Nj(u,v)在两个方向u,v上的求一阶偏导的结果，
Figure FDA0003822363780000052
表示u方向第j个一阶偏导光顺泛函，
Figure FDA0003822363780000053
表示v方向第j个一阶偏导光顺泛函；
上述两个光顺泛函计算为：
Figure FDA0003822363780000054
Figure FDA0003822363780000055
其中，um表示u方向第m个参数值，vm表示v方向第m个参数值，m表示数据点的序数；
B)第j个控制点的基于二阶偏导的光顺向量计算为：
Figure FDA0003822363780000056
式中：Nuu,j(u,v)、Nuv,j(u,v)和Nvv,j(u,v)分别表示表示对第j个NURBS基函数Nj(u,v)求二阶偏导的结果，
Figure FDA0003822363780000057
表示第j个第一类二阶偏导光顺泛函，
Figure FDA0003822363780000058
表示第j个第二类二阶偏导光顺泛函，
Figure FDA0003822363780000059
表示第j个第三类二阶偏导光顺泛函；
上述三个光顺泛函计算为：
Figure FDA00038223637800000510
Figure FDA00038223637800000511
Figure FDA00038223637800000512
C)第j个控制点的基于三阶偏导的光顺向量计算为：
Figure FDA00038223637800000513
式中：Nuuu,j(u,v)、Nuuv,j(u,v)、Nuvv,j(u,v)和Nvvv,j(u,v)分别表示对第j个NURBS基函数Nj(u,v)求三阶偏导的结果，
Figure FDA00038223637800000514
表示第j个第一类三阶偏导光顺泛函，
Figure FDA00038223637800000515
表示第j个第二类三阶偏导光顺泛函，
Figure FDA00038223637800000516
表示第j个第三类三阶偏导光顺泛函，
Figure FDA00038223637800000517
表示第j个第四类三阶偏导光顺泛函；
上述四个光顺泛函计算为：
Figure FDA0003822363780000061
Figure FDA0003822363780000062
Figure FDA0003822363780000063
Figure FDA0003822363780000064
CN202211045861.2A 2022-08-30 2022-08-30 交通工具外壳造型中nurbs曲线曲面光顺生成方法 Pending CN115358009A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211045861.2A CN115358009A (zh) 2022-08-30 2022-08-30 交通工具外壳造型中nurbs曲线曲面光顺生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211045861.2A CN115358009A (zh) 2022-08-30 2022-08-30 交通工具外壳造型中nurbs曲线曲面光顺生成方法

Publications (1)

Publication Number Publication Date
CN115358009A true CN115358009A (zh) 2022-11-18

Family

ID=84004427

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211045861.2A Pending CN115358009A (zh) 2022-08-30 2022-08-30 交通工具外壳造型中nurbs曲线曲面光顺生成方法

Country Status (1)

Country Link
CN (1) CN115358009A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117274425A (zh) * 2023-08-31 2023-12-22 武汉理工大学 一种基于几何特征的样条曲线生成方法、系统及终端
CN117371222A (zh) * 2023-10-23 2024-01-09 内蒙古工业大学 一种基于离散曲率及多级优化点的曲线光顺优化方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117274425A (zh) * 2023-08-31 2023-12-22 武汉理工大学 一种基于几何特征的样条曲线生成方法、系统及终端
CN117274425B (zh) * 2023-08-31 2024-05-24 武汉理工大学 一种基于几何特征的样条曲线生成方法、系统及终端
CN117371222A (zh) * 2023-10-23 2024-01-09 内蒙古工业大学 一种基于离散曲率及多级优化点的曲线光顺优化方法
CN117371222B (zh) * 2023-10-23 2024-04-26 内蒙古工业大学 一种基于离散曲率及多级优化点的曲线光顺优化方法

Similar Documents

Publication Publication Date Title
CN115358009A (zh) 交通工具外壳造型中nurbs曲线曲面光顺生成方法
CN109711048B (zh) 一种考虑气动和结构多类型设计变量的机翼设计方法
Ren et al. Combined reparameterization-based spiral toolpath generation for five-axis sculptured surface machining
CN104898564B (zh) 一种降低三轴联动轮廓误差的方法
Maute et al. Adaptive topology optimization of shell structures
CN108279643A (zh) 一种基于测量点与自适应差分进化算法的工件姿态调整方法
CN110543654B (zh) 激光喷丸成形复杂曲面分布式工艺参数的确定方法
Bohez et al. Adaptive nonlinear tool path optimization for five-axis machining
CN112036041B (zh) 一种3dp工艺中stl模型渗透误差综合补偿方法
CN111709097A (zh) 一种基于零亏格网格曲面连续变形的柔顺机构生成方法
Makhanov et al. On the tool-path optimization of a milling robot
CN112132943B (zh) 一种面向3d打印的过程纹理合成系统及方法
Zheng et al. Interpolation over arbitrary topology meshes using a two-phase subdivision scheme
CN105631158A (zh) 一种基于流线场理论的压气机叶片吸力面建模方法
Makhanov et al. Grid generation as applied to optimize cutting operations of the five-axis milling machine
CN112318499A (zh) 一种异型预制体机器人针刺成形路径规划方法
Wang et al. Tool path generation algorithm based on covariant field theory and cost functional optimization and its applications in blade machining
CN111553034A (zh) 一种基于插值方法的自由变形参数化方法
CN110689620B (zh) 一种多层次优化的网格曲面离散样条曲线设计方法
CN111610751B (zh) 过点集nurbs插值曲线的插值误差多次细分迭代计算方法
CN111079353B (zh) 应用于复杂流体分析的快速均匀网格划分的方法及装置
CN111400798A (zh) 一种建筑3d切片建造方法及3d打印装置
CN117251946A (zh) Nurbs曲面光顺调控的自由曲面抛光驻留时间计算方法及系统
CN117162492A (zh) 一种基于机械臂的曲面3d打印路径规划方法
CN116720268A (zh) 一种周期性描述的叶片叶型全局光顺重建方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination