CN102938008A - 一种人体建模中曲面拟合的进化算法 - Google Patents

一种人体建模中曲面拟合的进化算法 Download PDF

Info

Publication number
CN102938008A
CN102938008A CN2012104761587A CN201210476158A CN102938008A CN 102938008 A CN102938008 A CN 102938008A CN 2012104761587 A CN2012104761587 A CN 2012104761587A CN 201210476158 A CN201210476158 A CN 201210476158A CN 102938008 A CN102938008 A CN 102938008A
Authority
CN
China
Prior art keywords
curved surface
fitting
reference mark
evolutionary algorithm
fitness
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
CN2012104761587A
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.)
Dalian Polytechnic University
Original Assignee
Dalian Polytechnic University
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 Dalian Polytechnic University filed Critical Dalian Polytechnic University
Priority to CN2012104761587A priority Critical patent/CN102938008A/zh
Publication of CN102938008A publication Critical patent/CN102938008A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Image Generation (AREA)

Abstract

本发明公开了一种人体建模中曲面拟合的进化算法,属于计算机应用技术领域。本发明的原理及有益效果:进化算法是一种新兴的搜索寻优技术,提供了一种求解复杂系统的通用框架。本发明将基于进化算法的B样条曲面拟合应用到人体建模中,首先利用进化算法的B样条曲线拟合方法确定曲线边界线,从而得到基面,再按照曲面拟合的基本原理反求出曲面控制点,依据进化算法对控制点适应度大小的判断进行优化选择,最终实现人体建模中的曲面拟合。该方法利用进化算法的鲁棒性提高计算速度和效率,简化求取过程,降低拟合过程产生的误差,进而证明了该方法的优越性。

Description

一种人体建模中曲面拟合的进化算法
技术领域
本发明涉及一种人体建模中曲面拟合进化算法,属于计算机应用技术领域。
背景技术
曲面拟合是计算机图形学中的重要研究内容之一,在工程领域中有着广泛的应用。随着人们对计算机图形显示的真实性、实时性和交互性等方面要求的日益提高,曲面拟合的表示方法被不断地开拓创新。样条曲面是实体造型,样条曲面是在虚拟现实等CAD/CAM领域中普遍使用的几何造型工具。目前CAD系统中的曲面主要是B样条曲面和NURBS曲面,其中,B样条曲面是通过逼近一组控制点来生成,并具有两个最大的优点:(1)B样条多项式次数可独立于控制点数目;(2)B样条允许局部控制曲面。B样条曲面拟合方法中基于四边域曲面的几何重构方法是曲面重建中常用的方法。在CAD中,从点云中重建出被测物体表面的三角网格模型后,再将三角网格模型转换为用B样条表示的曲面。一般针对散乱点数据进行四边域曲面拟合,先根据散乱的数据构造一张基面,将数据点三角网格参数化,然后最小化测量数据到目标曲面的误差的平方和来判断所求目标曲面的精度,但一般的求取过程运算量较大,并且没有目的性。
发明内容
为了克服上述的不足,本发明提出一种人体建模中曲面拟合的进化算法。
本发明采用如下技术方案:
在计算过程中,曲面控制顶点的反求是通过两个阶段曲线反算得到的。
第一步:根据复杂曲面上的呈矩形阵列的数据点pi,j,(i=0,1,...,r;j=0,1,...,s),反求出控制点di,j(i=0,1,...,m;j=0,1,...,s),对控制顶点和数据点横、纵向进行平均规范积累弦长参数化从而得到节点矢量,这样就确定了一个B样条拟合曲面,该曲面逼近于已知的复杂曲面,但由于选取的数据点的数量、位置有可能存在噪声和摆动,所以计算出的拟合曲面在某些点处会过分偏离原有曲面。这时根据误差判断条件,实时调整偏离点处对应的控制顶点,便可修改拟合曲面的局部性状,达到更好的拟合效果。由于双三次B样条曲面具有良好的局部特性,工程实践中常常使用双三次B样条曲面来拟合复杂曲面,所以本发明采用双三次B样条曲面进行解释并运用四边域曲面拟合方法对其拟合。
四边界曲面拟合的具体步骤如下所示是:
1、边界预处理
2、空间散乱点的参数化
3、三角网格的建立
4、建立数学模型
第二步:为了得到散乱数据点的逼近曲面模型,散乱点的个数必须大于所求曲面的未知控制顶点数(l>m×n),这样才能保证未知数变量小于或等于方程数,目标函数为:
I r = α Σ k = 0 j [ S r ( u k r - 1 , v k r - 1 ) - P k ] 2 + β [ S uu r ( u k r - 1 , v k r - 1 ) 2 + S vv r ( u k r - 1 , v k r - 1 ) 2 + 2 S uv r ( u k r - 1 , v k r - 1 ) 2 ] - - - ( 1.1 )
其中:
Sr(u,0)=p(u,0);Sr(u,1)=p(u,1)
Sr(0,v)=p(0,v);Sr(1,v)=p(1,v)    (1.2)
Sv r(u,0)=pv(u,0);Sv r(u,1)=pv(u,1)
Su r(0,v)=pu(0,v);Su r(1,v)=pu(1,v)
r表示曲面的迭代次数;sr(u,v)是经过r次迭代后得到的曲面,并且 s r ( u , v ) = Σ i = 0 m Σ j = 0 n N i , 3 ( u ) N j , 3 ( v ) V i , j r ;
Figure GDA00002444850300032
是曲面Sr(u,v)的u,v两向量的二阶导数和混合偏导数;Pk(k=0,1,...l)为散乱数据点;
Figure GDA00002444850300033
是散乱点Pk在r-1次迭代曲面上的参数值;α为逼近权,β为光顺权,且α+β=1。当α值越大时曲线越靠近逼近点,则β较小,曲线越不光顺;当β越大时,线越光顺,则α变小曲线到数据点的距离则越大,所以应调节好α和β平衡曲线光顺和精度的关系。
由以上叙述可知目标曲面的数学模型公式(1.1)是带有复杂约束条件的二次优化问题,我们的目的就是运用寻优手段使目标函数Ir既能满足约束,又能达到最小,由此就可求得第r次几何迭代目标曲面的控制顶点为
Figure GDA00002444850300034
(i=0,1,…,m;j=0,1,…,n)。
第三步:数学模型的求解即目标函数Ir的求解
对于公式(1.1)的目标函数Ir的优化我们分为两步进行。由于约束条件只涉及到边界曲线和跨边界切矢曲线,而曲面的边界和跨边界切矢可以由曲面的边界和次边界控制顶点完全确定,因此我们可以先处理这些约束条件来得到目标曲面的边界和次边界控制顶点,根据公式(1.1)就可以得到我们所需要的拟合曲面。
曲面重建首先将离散点参数化,设计为三角网格进而初步拟合成第一层曲面(基面)所示。采用Delaunay三角化方法对散乱点进行三角网格拟合,以人体模型数据为例,优化整体曲面的边界线和基曲面。
然后利用公式(1.1)、(1.2)运用最小二乘法反求出曲面控制点,将控制点进行编码产生初始群体,通过适应度函数f=1/|(F(x,y,z)-z)2+(y(x,y)-y)2+(x-xi)2-ε|求得最大值,也就是最终的曲面控制点以便拟合出基面。
因为人体模型是复杂样件模型,并不能用一张曲面表达,因此将曲面细化,就是从粗糙到精细的B样条控制网格确定,最后使B样条曲面逐步逼近给定的散乱点云,而且由粗到细的过程仅在误差超过预定的区域进行,具体流程如图3所示。
首先对反求出的控制点进行编码,对控制点数据分块得到相互连接的矩形区域,为保证曲面的连续性,将每块曲面的左上角和右下角数值记录,以便确定细化曲面的范围,对分块后曲面内的数据建立初始群体,判断它的迭代次数是否超出给定范围(初次分块迭代次数为零),如果超出则进入到上层计算重复分块动作,如果没有则接着从初始群体(含有N个基因)中随机选取第0代群体(含有n个),对其拟合并求出适应度,如果适应度达到给出的标准则停止循环,直接得到拟合曲面,如果不符合适应度的要求则根据适应度的大小对数据进行遗传和变异操作,判断新生代基因个数,如果小于初始群体个数则从余下的数据中随机选取使新生代群体达到n个,从而得到第一代群体,再进入到循环中直到满足停机条件得到最优的拟合曲面。
本发明的原理及有益效果:进化算法是一种新兴的搜索寻优技术,进化算法和传统方法的最大不同,就是它不要求所研究的问题连续或可导,能够很快地得到所求的最优解,适用于处理传统搜索方法难以解决的复杂和非线性问题。进化算法的整体搜索策略和优化搜索方法在计算工程中不依赖于梯度信息或是其他辅助知识,而只需要对搜索反应有影响的目标函数和相应的适应度函数,所以进化算法提供了一种求解复杂体统的通用框架。本发明将基于进化算法的B样条曲面拟合应用到人体建模中,首先利用进化算法的B样条曲线拟合方法确定曲线边界线,从而得到基面,再按照曲面拟合的基本原理反求出曲面控制点,依据进化算法对控制点适应度大小的判断进行优化选择,最终实现人体建模中的曲面拟合。该方法利用进化算法的鲁棒性提高计算速度和效率,简化求取过程,降低拟合过程产生的误差,进而证明了该方法的优越性。
附图说明
图1三角网格化。
图2基面。
图3具体流程图。
图4.人体模型的曲面拟合。
具体实施方式
下面结合附图对本发明做进一步说明,本发明的一个实施例:
在计算过程中,曲面控制顶点的反求是通过两个阶段曲线反算得到的。
第一步:根据复杂曲面上的呈矩形阵列的数据点pi,j,(i=0,1,...,r;j=0,1,...,s),反求出控制点di,j(i=0,1,...,m;j=0,1,...,s),对控制顶点和数据点横、纵向进行平均规范积累弦长参数化从而得到节点矢量,这样就确定了一个B样条拟合曲面,该曲面逼近于已知的复杂曲面,但由于选取的数据点的数量、位置有可能存在噪声和摆动,所以计算出的拟合曲面在某些点处会过分偏离原有曲面。这时根据误差判断条件,实时调整偏离点处对应的控制顶点,便可修改拟合曲面的局部性状,达到更好的拟合效果。由于双三次B样条曲面具有良好的局部特性,工程实践中常常使用双三次B样条曲面来拟合复杂曲面,所以本发明采用双三次B样条曲面进行解释并运用四边域曲面拟合方法对其拟合。
四边界曲面拟合的具体步骤如下所示是:
1、边界预处理
2、空间散乱点的参数化
3、三角网格的建立
4、建立数学模型
第二步:为了得到散乱数据点的逼近曲面模型,散乱点的个数必须大于所求曲面的未知控制顶点数(l>m×n),这样才能保证未知数变量小于或等于方程数,目标函数为:
I r = α Σ k = 0 j [ S r ( u k r - 1 , v k r - 1 ) - P k ] 2 + β [ S uu r ( u k r - 1 , v k r - 1 ) 2 + S vv r ( u k r - 1 , v k r - 1 ) 2 + 2 S uv r ( u k r - 1 , v k r - 1 ) 2 ] - - - ( 1.1 )
其中:
Sr(u,0)=p(u,0);Sr(u,1)=p(u,1)
Sr(0,v)=p(0,v);Sr(1,v)=p(1,v)    (1.2)
Sv r(u,0)=pv(u,0);Sv r(u,1)=pv(u,1)
Su r(0,v)=pu(0,v);Su r(1,v)=pu(1,v)
r表示曲面的迭代次数;sr(u,v)是经过r次迭代后得到的曲面,并且 s r ( u , v ) = Σ i = 0 m Σ j = 0 n N i , 3 ( u ) N j , 3 ( v ) V i , j r ;
Figure GDA00002444850300063
是曲面Sr(u,v)的u,v两向量的二阶导数和混合偏导数;Pk(k=0,1,..l)为散乱数据点;
Figure GDA00002444850300064
是散乱点Pk在r-1次迭代曲面上的参数值;α为逼近权,β为光顺权,且α+β=1。当α值越大时曲线越靠近逼近点,则β较小,曲线越不光顺;当β越大时,线越光顺,则α变小曲线到数据点的距离则越大,所以应调节好α和β平衡曲线光顺和精度的关系。
由以上叙述可知目标曲面的数学模型公式(1.1)是带有复杂约束条件的二次优化问题,我们的目的就是运用寻优手段使目标函数Ir既能满足约束,又能达到最小,由此就可求得第r次几何迭代目标曲面的控制顶点为
Figure GDA00002444850300065
(i=0,1,…,m;j=0,1,…,n)。
第三步:数学模型的求解
对于公式(1.1)的目标函数Ir的优化我们分为两步进行。由于约束条件只涉及到边界曲线和跨边界切矢曲线,而曲面的边界和跨边界切矢可以由曲面的边界和次边界控制顶点完全确定,因此我们可以先处理这些约束条件来得到目标曲面的边界和次边界控制顶点,根据公式(1.1)就可以得到我们所需要的拟合曲面。
曲面重建首先将离散点参数化,设计为三角网格进而初步拟合成第一层曲面(基面)如图1所示。采用Delaunay三角化方法对散乱点进行三角网格拟合,以人体模型数据为例,优化整体曲面的边界线和基曲面。
然后利用公式(1.1)、(1.2)运用最小二乘法反求出曲面控制点,将控制点进行编码产生初始群体,通过适应度函数f=1/|(F(x,y,z)-z)2+(y(x,y)-y)2+(x-xi)2-ε|求得最大值,也就是最终的曲面控制点以便拟合出基面,如图2所示。
因为人体模型是复杂样件模型,并不能用一张曲面表达,因此将曲面细化,就是从粗糙到精细的B样条控制网格确定,最后使B样条曲面逐步逼近给定的散乱点云,而且由粗到细的过程仅在误差超过预定的区域进行,具体流程如图3所示。
首先对反求出的控制点进行编码,对控制点数据分块得到相互连接的矩形区域,为保证曲面的连续性,将每块曲面的左上角和右下角数值记录,以便确定细化曲面的范围,对分块后曲面内的数据建立初始群体,判断它的迭代次数是否超出给定范围(初次分块迭代次数为零),如果超出则进入到上层计算重复分块动作,如果没有则接着从初始群体(含有N个基因)中随机选取第0代群体(含有n个),对其拟合并求出适应度,如果适应度达到给出的标准则停止循环,直接得到拟合曲面,如果不符合适应度的要求则根据适应度的大小对数据进行遗传和变异操作,判断新生代基因个数,如果小于初始群体个数则从余下的数据中随机选取使新生代群体达到n个,从而得到第一代群体,再进入到循环中直到满足停机条件得到最优的拟合曲面。
以双三次B样条曲面为例,方程为:
p ( u , v ) = Σ i = 1 m Σ j = 1 n d i , j N i , 3 ( u ) N j , 3 ( v ) - - - ( 1.3 )
其中节点向量U和V,边界的控制点
Figure GDA00002444850300082
由第三章内容得出,曲面数据点
Figure GDA00002444850300083
及将其参数化为
Figure GDA00002444850300084
为解方程要求控制点di,j(i=2,3,…,m-1;j=2,3,…,n-1)在最小二乘意义下逼近
Figure GDA00002444850300085
即:
f = Σ k L [ Q k - p ( u k , v k ) ] 2 = Σ k L [ Q k - Σ i = 1 m Σ j = 1 n d i , j N i , 3 ( u k ) N j , 3 ( v k ) ] 2 = min - - - ( 1.4 )
f = Σ k L [ Q k - Σ i = 1 m d i , 1 N i , 3 ( u k ) N 1,3 ( v k ) - Σ i m d i , n N i , 3 ( u k ) N n , 3 ( v k ) - Σ j = 2 n - 1 d i , n N i , 3 ( u k ) N j , 3 ( v k )
- Σ j = 2 n - 1 d i , n N m , 3 ( u k ) N j , 3 ( v k ) - Σ i = 2 m - 1 Σ j = 2 n - 1 d i , j N i , 3 ( u k ) N j , 3 ( v k ) ] 2 - - - ( 1.5 )
令:
r k = Q k - Σ i = 1 m d i , 1 N i , 3 ( u k ) N 1,3 ( v k ) - Σ i = 1 m d i , n N i , 3 ( u k ) N n , 3 ( v k )
- Σ j = 2 n - 1 d i , n N 1,3 ( u k ) N j , 3 ( v k ) - Σ j = 2 n - 1 d i , n N m , 3 ( u k ) N j , 3 ( v k ) - - - ( 1.6 )
d(i-2)×(n-2)+j=di,j,t(i-2)×(n-2)+j(wk)=Ni,3(uk)Nj,3(vk)
(i=2,3,...,m-1;j=2,3,...,n-1)    (1.7)
可得:
f = Σ k = 1 L [ r k - Σ i = 2 m - 1 Σ j = 2 n - 1 d i , j N i , 3 ( u k ) N j , 3 ( v k ) ] 2 = Σ k = 1 L [ r k - Σ s = 1 ( m - 2 ) × ( n - 2 ) d s t s ( w k ) ] 2 - - - ( 1.8 )
式中:
Figure GDA000024448503000812
D = d 1 . . . d ( m - 2 ) × ( n - 2 ) ( m - 2 ) × ( n - 2 ) × s R = r 1 . . . r L L × s
公式1.4则可改写成f=(R-FD)T(R-FD)=min欲使目标函数f最小,应使它关于控制顶点D的导数为零,则可解得。
D=F-R    (1.9)
假设f的值为α,对每一块曲面重新拟合分别得到f1,f2,f3,f4,使他们分别与α比较,以fn-a=min,(n=1,2,3,4)为适应度进行判断,对该曲面内的控制点进行遗传、变异操作,从而选取最优的控制点拟合第n块曲面,当迭代次数超过30时或fn-a>0,(n=1,2,3,4)则将第n块曲面平分为4块,在平分曲面时,应将该区面的左上坐标和右下坐标记录,并以这两点作为面片的边界点之一,这样在拟合结束时较容易拼接曲面片并保持拟合后的整体曲面的连续性,重复以上步骤直到得到最优曲面,如图4所示。
基于进化算法B样条曲面拟合是一种自适应细化B样条曲面拟合法,提高了拟合精度,并且具有目的性,一旦达到早期设定的精度则立即停止迭代输出曲面结果,减少了数据的计算量,本发明将其应用在人体建模,并证明了它的优越性。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。

Claims (1)

1.一种人体建模中曲面拟合的进化算法,其特征在于:包括以下几个步骤
第一步:根据复杂曲面上的呈矩形阵列的数据点pi,j,(i=0,1,...,r;j=0,1,...,s),反求出控制点di,j(i=0,1,...,m;j=0,1,...,s),对控制顶点和数据点横、纵向进行平均规范积累弦长参数化从而得到节点矢量,这样就确定了一个B样条拟合曲面,该曲面逼近于已知的复杂曲面,但由于选取的数据点的数量、位置有可能存在噪声和摆动,所以计算出的拟合曲面在某些点处会过分偏离原有曲面;这时根据误差判断条件,实时调整偏离点处对应的控制顶点,便可修改拟合曲面的局部性状,达到更好的拟合效果,由于双三次B样条曲面具有良好的局部特性,工程实践中常常使用双三次B样条曲面来拟合复杂曲面;采用双三次B样条曲面进行解释并运用四边域曲面拟合方法对其拟合;
所述四边域曲面拟合方法步骤如下:
1、边界预处理,
2、空间散乱点的参数化,
3、三角网格的建立,
4、建立数学模型;
第二步:建立目标函数为:
I r = α Σ k = 0 j [ S r ( u k r - 1 , v k r - 1 ) - P k ] 2 + β [ S uu r ( u k r - 1 , v k r - 1 ) 2 + S vv r ( u k r - 1 , v k r - 1 ) 2 + 2 S uv r ( u k r - 1 , v k r - 1 ) 2 ] - - - ( 1.1 )
所述:
Sr(u,0)=p(u,0);Sr(u,1)=p(u,1)
Sr(0,v)=p(0,v);Sr(1,v)=p(1,v)
Sv r(u,0)=pv(u,0);Sv r(u,1)=pv(u,1)
Su r(0,v)=pu(0,v);Su r(1,v)=pu(1,v)
散乱点的个数必须大于所求曲面的未知控制顶点数(l>m×n),r表示曲面的迭代次数;sr(u,v)是经过r次迭代后得到的曲面,并且
Figure FDA00002444850200021
Figure FDA00002444850200022
是曲面Sr(u,v)的u,v两向量的二阶导数和混合偏导数;Pk(k=0,1,...l)为散乱数据点;是散乱点Pk在r-1次迭代曲面上的参数值;α为逼近权,β为光顺权,且α+β=1;当α值越大时曲线越靠近逼近点,则β较小,曲线越不光顺;当β越大时,线越光顺,则α变小曲线到数据点的距离则越大,所以应调节好α和β平衡曲线光顺和精度的关系;所述可知目标曲面的目标函数是带有复杂约束条件的二次优化问题使目标函数Ir既能满足约束,又能达到最小,由此就可求得第r次几何迭代目标曲面的控制顶点为
Figure FDA00002444850200024
(i=0,1,…,m;j=0,1,…,n);
第三步:目标函数Ir的求解
目标函数Ir的优化分为两步;
1.采用Delaunay三角化方法对散乱点进行三角网格拟合,以人体模型数据为例,优化整体曲面的边界线和基曲面;
2.利用公目标函数运用最小二乘法反求出曲面控制点,将控制点进行编码产生初始群体,通过适应度函数f=1/|(F(x,y,z)-z)2+(y(x,y)-y)2+(x-xi)2-ε|求得最大值,即最终的曲面控制点,以便拟合出基面;
所述拟合出基面的过程:首先对反求出的控制点进行编码,对控制点数据分块得到相互连接的矩形区域,为保证曲面的连续性,将每块曲面的左上角和右下角数值记录,以便确定细化曲面的范围,对分块后曲面内的数据建立初始群体,判断它的迭代次数是否超出给定范围,即初次分块迭代次数为零;如果超出则进入到上层计算重复分块动作,如果没有则接着从初始群体中随机选取第0代群体,所述第0代群体含有n个,所述初始群体含有N个基因,第0代群体含有n个;对其拟合并求出适应度,如果适应度达到给出的标准则停止循环,直接得到拟合曲面,如果不符合适应度的要求则根据适应度的大小对数据进行遗传和变异操作,判断新生代基因个数,如果小于初始群体个数则从余下的数据中随机选取使新生代群体达到n个,从而得到第一代群体,再进入到循环中直到满足停机条件得到最优的拟合曲面。
CN2012104761587A 2012-11-21 2012-11-21 一种人体建模中曲面拟合的进化算法 Pending CN102938008A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2012104761587A CN102938008A (zh) 2012-11-21 2012-11-21 一种人体建模中曲面拟合的进化算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2012104761587A CN102938008A (zh) 2012-11-21 2012-11-21 一种人体建模中曲面拟合的进化算法

Publications (1)

Publication Number Publication Date
CN102938008A true CN102938008A (zh) 2013-02-20

Family

ID=47696905

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012104761587A Pending CN102938008A (zh) 2012-11-21 2012-11-21 一种人体建模中曲面拟合的进化算法

Country Status (1)

Country Link
CN (1) CN102938008A (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103676786A (zh) * 2013-12-23 2014-03-26 北京航空航天大学 一种基于加速度原理的曲线光顺方法
CN105160700A (zh) * 2015-06-18 2015-12-16 上海工程技术大学 一种用于三维模型重建的截面曲线重构方法
CN107798732A (zh) * 2017-10-27 2018-03-13 中国工程物理研究院应用电子学研究所 一种自由曲面形态控制方法
CN109035363A (zh) * 2018-07-05 2018-12-18 凌云光技术集团有限责任公司 一种快速迭代的线圆最优拟合方法
CN109657362A (zh) * 2018-12-22 2019-04-19 上海杰达齿科制作有限公司 修复体瓷材料层的缩放方法及加工工艺
CN109961517A (zh) * 2019-03-01 2019-07-02 浙江大学 一种用于参数曲面拟合的三角网格重参数化方法
CN110211201A (zh) * 2019-05-07 2019-09-06 北京航空航天大学 一种b样条曲面重建方法
CN110688768A (zh) * 2019-10-11 2020-01-14 广东工业大学 一种新能源汽车空调压缩机型线优化方法、装置和设备
CN114341942A (zh) * 2019-09-02 2022-04-12 马自达汽车株式会社 曲面的拟合处理方法、拟合处理装置和拟合处理程序以及存储有该拟合处理程序的计算机可读存储介质

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
JINGSHENG LEI等: "The Application of Evolutionary Algorithm in B-Spline Curved Surface Fitting", 《ARITFICIAL INTELLIGENCE AND COMPUTATIONAL INTELLIGENCE》 *
刘佳琪等: "进化算法在最小二乘法曲线拟合中的应用", 《计算机应用》 *
施美雅: "基于NURBS的服装人体建模技术", 《电脑知识与技术》 *
黄海峤等: "基于可展曲面的3D服装原型建模与服装样板生成", 《东华大学学报(自然科学版)》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103676786B (zh) * 2013-12-23 2016-05-25 北京航空航天大学 一种基于加速度原理的曲线光顺方法
CN103676786A (zh) * 2013-12-23 2014-03-26 北京航空航天大学 一种基于加速度原理的曲线光顺方法
CN105160700A (zh) * 2015-06-18 2015-12-16 上海工程技术大学 一种用于三维模型重建的截面曲线重构方法
CN105160700B (zh) * 2015-06-18 2018-04-24 上海工程技术大学 一种用于三维模型重建的截面曲线重构方法
CN107798732B (zh) * 2017-10-27 2021-03-12 中国工程物理研究院应用电子学研究所 一种自由曲面形态控制方法
CN107798732A (zh) * 2017-10-27 2018-03-13 中国工程物理研究院应用电子学研究所 一种自由曲面形态控制方法
CN109035363A (zh) * 2018-07-05 2018-12-18 凌云光技术集团有限责任公司 一种快速迭代的线圆最优拟合方法
CN109035363B (zh) * 2018-07-05 2023-11-10 凌云光技术股份有限公司 一种快速迭代的线圆最优拟合方法
CN109657362A (zh) * 2018-12-22 2019-04-19 上海杰达齿科制作有限公司 修复体瓷材料层的缩放方法及加工工艺
CN109961517B (zh) * 2019-03-01 2021-05-07 浙江大学 一种用于参数曲面拟合的三角网格重参数化方法
CN109961517A (zh) * 2019-03-01 2019-07-02 浙江大学 一种用于参数曲面拟合的三角网格重参数化方法
CN110211201A (zh) * 2019-05-07 2019-09-06 北京航空航天大学 一种b样条曲面重建方法
CN114341942A (zh) * 2019-09-02 2022-04-12 马自达汽车株式会社 曲面的拟合处理方法、拟合处理装置和拟合处理程序以及存储有该拟合处理程序的计算机可读存储介质
CN110688768A (zh) * 2019-10-11 2020-01-14 广东工业大学 一种新能源汽车空调压缩机型线优化方法、装置和设备

Similar Documents

Publication Publication Date Title
CN102938008A (zh) 一种人体建模中曲面拟合的进化算法
Ribeiro et al. An airfoil optimization technique for wind turbines
Bates et al. A new method for moving–boundary hydrodynamic problems in shallow water
Belibassakis et al. A BEM-isogeometric method for the ship wave-resistance problem
Peraire et al. Unstructured finite element mesh generation and adaptive procedures for CFD
Jiang et al. A sharp-interface immersed smoothed finite element method for interactions between incompressible flows and large deformation solids
US11629694B2 (en) Wind turbine model based control and estimation with accurate online models
Jia et al. Extended isogeometric analysis for material interface problems
Becker et al. An advanced NURBS fitting procedure for post-processing of grid-based shape optimizations
Gao et al. Grid generation on free-form surface using guide line advancing and surface flattening method
Berger Cut cells: Meshes and solvers
Hou et al. Multislope MUSCL method applied to solve shallow water equations
Wang et al. A triangular grid generation and optimization framework for the design of free-form gridshells
Xu et al. Wing-body junction optimisation with CAD-based parametrisation including a moving intersection
López et al. Isogeometric structural shape optimization using automatic sensitivity analysis
Isola et al. Finite-volume solution of two-dimensional compressible flows over dynamic adaptive grids
Masters et al. Progressive subdivision curves for aerodynamic shape optimisation
Kim et al. A B-spline based higher order panel method for analysis of steady flow around marine propellers
CN105427362A (zh) 一种快速aiap形状插值算法
Fanxi et al. A fast and automatic full-potential finite volume solver on Cartesian grids for unconventional configurations
Sheng et al. Multiblock approach for calculating incompressible fluid flows on unstructured grids
Li et al. Multi-objective optimization of turbine blade profiles based on multi-agent reinforcement learning
Bandala et al. Implementation of varied particle container for Smoothed Particle Hydrodynamics-based aggregation for unmanned aerial vehicle quadrotor swarm
Masters et al. Three-dimensional subdivision parameterisation for aerodynamic shape optimisation
Wei et al. Hexahedral mesh smoothing via local element regularization and global mesh optimization

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C05 Deemed withdrawal (patent law before 1993)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20130220