CN103577636B - 一种基于曲率特征的叶片曲面统一离散方法 - Google Patents
一种基于曲率特征的叶片曲面统一离散方法 Download PDFInfo
- Publication number
- CN103577636B CN103577636B CN201310499607.4A CN201310499607A CN103577636B CN 103577636 B CN103577636 B CN 103577636B CN 201310499607 A CN201310499607 A CN 201310499607A CN 103577636 B CN103577636 B CN 103577636B
- Authority
- CN
- China
- Prior art keywords
- line
- curvature
- prime
- blade
- parameter
- 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.)
- Active
Links
Landscapes
- Structures Of Non-Positive Displacement Pumps (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
一种基于曲率特征的叶片曲面统一离散方法,步骤为:一选定截面线r,确定圆弧段中心点B1,B2;二从B1,B2点求出纵向曲率线L2,L6;三从B1,B2两侧找到4点V1,V2,V3,V4;四在叶片两端求出经过该4点的纵向曲率线L1,L3,L5,L7;五L1,L2,L3作为等参数线,对L1,L3中间区域参数化,确定边界线l1,l3;L5,L6,L7作为等参数线,对L5,L7中间区域参数化,确定边界线l5,l7;六分别确定截面线r上圆弧中心点M1,M2;七确定经过M1,M2纵向曲率线L4,L8;八以L2,L4,L6,L8和l1,l3,l5,l7为新参数域的等参数线,重新参数化;九根据新参数对曲面进行非均匀离散。
Description
技术领域
本发明提供一种基于曲率特征的叶片曲面统一离散方法,它是一种用于宽行数控加工中用于刀轨规划的工件曲面的离散方法,它涉及一种宽行数控加工中离散刀位点的选取方法,属于曲面数控加工技术领域。
背景技术
等参数线法是数控加工中常用的刀轨规划方法之一,该方法主要是让刀具沿着工件曲面的一个参数方向运动。等参数线法依赖于曲面的参数化。
在常用的CAM***中,曲面造型一般采用的是弧长参数化,该种参数化方法没有考虑到曲面的弯曲情况,因此可能会产生参数线扭曲或者波动的情况,不利于实际加工。而基于曲率线网的参数化方法能够很好的兼顾曲面的局部弯曲特征,但是在实际的复杂曲面时,曲面的局部造型误差会影响到曲率线的全局走势,从而很难得到可信的曲率线网格。
在实际的宽行加工编程***中,一般先对工件曲面进行离散。传统的曲面离散方法中,均匀离散方法是直接在选取均匀参数网格,非均匀离散方法则主要是通过曲面的局部特征(如弧角变化)来调整离散网格的密度。上述传统的曲面离散方法都是基于曲面既有参数进行离散的,没有解决曲面参数线扭曲的所带来的问题。
一般叶片曲面的参数线扭曲一般发生在进排气边附近,而进排气边附近的圆弧过渡的中心位置处的曲率线走势因其两个主曲率差距大而一般比较稳定,受曲面造型的误差影响较小,而且该处的曲率线一般不会发生扭曲和翻越叶盆叶背的情况,具有叶盆叶背“分界线”的特点。本专利申请考虑了曲面参数线扭曲的情况,结合叶片曲面的上述特征,提出一种基于有限条曲率线的重新参数化方法。并且给出在重新参数化的基础上对曲面进行离散的过程。
发明内容
1、目的:本发明的目的是提供一种基于曲率特征的叶片曲面统一离散方法,用以解决叶片进排气边附近参数线扭曲造成的刀轨波动问题。
2、技术方案:本发明主要是通过以下技术方案来实现。
本发明一种基于曲率特征的叶片曲面统一离散方法,它包括以下步骤:
步骤一选定一个有代表性的截面线r,确定其上进排气边圆弧段的中心点B1,B2位置。具体过程如下:
选定叶片的某一条等参数线S(u,v0)作为叶片的一条截型线(如图2中标示有点的截型线所示),计算该截型线的曲率分布
其中L(u,v0)=Suu(u,v0)·n(u,v0),E(u,v0)=Su(u,v0)·Su(u,v0),n(u,v0)为曲面在点(u,v0)处的单位法矢。如图4所示,函数k(u)随着u的从小到大的会由两个突变,不妨设给定阈值系数F,定义当k(u)≥F×k0时,u处于突变位置之内。分别记录第一个突变的起始位置U1s和终止位置U1e,第二个突变的起始位置U2s和终止位置U2e。定义运算
令
则U2和U6即为圆弧过渡的中心位置。
本专利申请中如果不加说明,都有v0=0;std(a)运算保证了结果落在参数的定义域内;Add(a,b)则是考虑到突变区域跨越参数域边界的情况;k0实际上表示截型线的平均曲率;圆弧过渡的中心位置一般对应着截型线的最弯处;U2和U6的确定的也可以通过选取叶片多条截型线,分别通过上述步骤计算,再取平均值得到;U2和U6确定了图3中的点B1,B2。
步骤二从B1,B2点出发向叶片两端求出经过点B1,B2的两条纵向曲率线L2,L6。具体过程如下:
令u2,1=U2,u6,1=U6,v1=0.
如果己知曲率线上第i个点,首先确定第二条曲率线下一点的参数值。计算叶片曲面在点(u2,i,vi)处与V方向相近的主方向D2,i,则下一点P2,i+1=S(u2,i,vi)+step×D2,i,其中step为给定的步长,其大小影响V方向的离散精度。再计算点P2,i+1在曲面上对应的参数值(u2,i+1,vi+1)。当vi+1≥1时,上述步骤终止。
对于j=6,确定叶片曲面在点(uj,i,vi)处与V方向相近的主方向dj,i,则
从而第j条曲率线第i+1个点的参数值为(uj,i+1,vi+1)。
对于i=2,6,记Ui对应的曲率线为Li。
说明:曲率线是通过一系列给定步长的离散点表示;也可以选取其他曲率线确定vi+1取值;计算过程保证了两条曲率线在第i个点的V参数值是一致的;主方向D2,i为三维空间中的向量,而dj,i则为曲面上一方向的比值表示。
步骤三在代表截面线r上,从B1,B2两点出发,分别在该点两侧进排气边圆弧范围内找到4个指定点V1,V2,V3,V4。
给定距离值dis1,dis2,计算U1,U3,U5,U7,使得
则U1,U3,U5,U7即为所需求的圆弧上V1,V2,V3,V4的位置。
在求解积分时出现U1,U3,U4,U6与步骤一中的U1s,U1e,U2s,U2e意义不一样。实际叶片曲面在圆弧过渡边缘处的曲率走势一般会比较紊乱,所以需要通过参数dis1,dis2来调整U1,U3,U4,U6的取值,使得经过点(U1,v0),(U3,v0),(U4,v0),(U6,v0)处的四条与叶片造型方向一致的曲率线稳定可信。
步骤四从V1,V2,V3,V4出发向叶片两端求出分别经过V1,V2,V3,V4的四条纵向曲率线L1,L3,L5,L7。
u1,1=U1,u3,1=U3,u5,1=U5,u7,1=U7。
己知曲率线上第i个点,对于j=1,3,5,7,确定叶片曲面在点(uj,i,vi)处与V方向相近的主方向dj,i,则
从而第j条曲率线第i+1个点的参数值为(uj,i+1,vi+1)。
对于i=1,3,5,7,记Ui对应的曲率线为Li.
步骤五把L1,L2,L3作为等参数线,对L1,L3中间的区域进行重新参数化,确定重新参数化后的参数域边界线l1,l3;把L5,L6,L7作为等参数线,对L5,L7中间的区域进行重新参数化,确定重新参数化后的参数域边界线l5,l7。
设叶片进排气边的平均过渡圆弧弧长为ArcL,
定义
令
从而把L1,L3,L5,L7往远离进排气边方向(如图6中的箭头所示)平移至新的参数域边界。记平移后得到的曲线分别为l1,l3,l5,l7。
由于造型的误差,实际叶片曲面在进排气边过渡圆弧起末端区域的曲率线走势会比较紊乱,一般不能够直接求解通过该区域的曲率线,步骤四是把该区域外侧的曲率线通过距离比例来移动,从而得到经过该区域的近似曲率线;l1,l3,l5,l7通过图3中的点P1,P2,P3,P4。
步骤六分别确定截面线r上叶盆叶背圆弧的中心点M1,M2。
令
则U4,U8记录了点M1,M2的位置。
步骤七确定经过M1,M2的两条纵向曲率线或等参数线L4,L8。
L4,L8如果为曲率线,则求解过程如下:
己知曲率线上第i个点,对于j=4,8,确定叶片曲面在点(uj,i,vi)处与V方向相近的主方向dj,i,则
从而第j条曲率线第i+1个点的参数值为(uj,i+1,vi+1)。
对于i=4,8,记Ui对应的曲率线为Li.
如果L4,L8为等参数线,则对于j=4,8,uj,i+1=Uj。
步骤八以L2,L4,L6,L8和l1,l3,l5,l7为新参数域下的等参数线,对曲面进行重新参数化。
设重新参数化的后的曲面表达式为S′(u′,v′)。
将八条线看做是新的等U′参数线,首先确定第j条线在新参数化下对应的U′坐标u′j。令u′1=0,则
u′j=u′j-1+M(Uj,Uj-1),j=2,…,8.
记u′9=1。
给定一个U′V′坐标(u′,vi),其对应的UV坐标为
其中j满足u′j≤u′≤u′j+1.
八条线L2,L4,L6,L8和l1,l3,l5,l7实际上是逼近图2中叶片曲面上与V方向一致的八根线,这八根曲率线根据曲面的弯曲特征把叶片曲面分割了八块相对平稳的区域,从而可以保证每块区域中的参数线都不会出现过大的扭曲。
步骤九根据新参数对曲面进行非均匀离散。
V′方向的离散已经由步骤二确定,故只需讨论U′方向的离散。给定U′方向的离散弓高为GgU,最大离散间隔为Offset,记V′方向的离散点数为vnum。如果当前U′参数位置为U′i,对于1≤j≤vnum,计算点(U′i,vj)在UV坐标下的参数,从而确定该点处的沿U′方向的曲率uQLij。令
其中uLenj表示曲面上的曲线V=vj的弧长。
令offseti=min{minjOSi,j,Offset},则下一个离散位置为U′i+1=U′i+offseti。
说明:最大离散间隔Offset限制了离散的最大步长,以防止在离散叶背叶盆等较平区域时出现步长过大的情况;离散弓高的意义如图5所示;U方向的离散密度受到曲面局部曲率的影响,故为非均匀离散;V方向的离散由三维空间中的步长确定,可以近似认为为均匀离散。
叶片截型线一般由叶盆部分、叶背部分和两端小的过渡圆弧组成。其中叶片曲面造型圆弧过渡的中心位置是小过渡圆弧的中点。如图1所示。
其中纵向指的是叶片造型方向,是指叶片截型线的排列的参数方向,一般为V方向。如图2所示。
3、优点和功效:
本发明提供了一种基于曲率特征的叶片曲面统一离散方法,该方法考虑了叶片的造型特征,以八条曲线为基准对把叶片曲面划分成八个边界与曲率线走势一致的平坦的区域,再分别对每个区域进行重新参数化,使得每个区域的边界都为新参数域下的等参数线,解决了叶片在进排气边附近参数线扭曲严重的情况,从而实现叶片进排气加工的刀轨光顺稳定,一方面降低了叶片加工对机床加速度和机床联动的要求,另一方面提高了叶片进排气边加工的表面质量。另外,本发明在叶片较弯曲的方向采用非均匀的离散方法,在保证离散精度的基础上大大节省了存储空间,减少了计算量。同时本发明具有兼容性强、实现成本低的优点。
附图说明
图1:叶片截型线构成示意图。可以看出,叶片截型线一般由叶盆部分、叶背部分和两端小的过渡圆弧组成。
图2:叶片构造方向说明。叶片曲面一般是通过数条排列好的截型线生成,称与叶片截型线的排列的方向一致的参数方向为叶片造型方向,也称为叶片纵向。为方便起见,本专利申请说明书中一般认为叶片造型方向为V方向。
图3:叶片截型线分区点示意图。点P1,P2,P3,P4,M1,M2,B1,B2把截型线分成了八段,这八段截型线自然对应着八块叶片曲面区域。其中点P1,P2,P3,P4把截型线分成了曲率变化均匀的四段,这四段线对应的叶片曲面区域内部曲率变化较为均匀,对重新参数化后的参数线性质起了主要作用。
图4:叶片截型线的曲率分布示意图。叶片截型线曲率函数有两个明显的峰值,分别对应着叶片截型线的两个过渡圆弧的位置。
图5:离散弓高示意图。如果当前点为P1,则可以把2L作为离散步长,得到下一个离散点的坐标。因为要满足离散弓高的要求,此时P1点处的离散步长不能够大于2L。
图6:叶片曲面在新的参数域下的分区示意图。可以看出,曲率线L2,L4,L6,L8和近似曲率线l1,l3,l5,l7为新参数域下的等参数线,而且该八条等参数线自然把叶片划分为八块区域。
图7:叶片曲面统一离散流程。
图中代号、符号说明如下:
M1,M2-截型线上对应叶盆叶背区域的中点。
B1,B2-截型线上进排气边过渡圆弧的中点
V1,V2,V3,V4-在截型线上,从B1出发,在指定范围内确定点V1,V2;从B2出发,在指定范围内确定点V3,V4。
P1,P2,P3,P4-在截型线上,V1,V2,V3,V4按照指定规则平移得到点P1,P2,P3,P4。
U2,U6-对应点B1,B2的U方向参数值。
U1,U3,U5,U7-对应点V1,V2,V3,V4的U方向参数值。
U4,U8-对应点M1,M2的U方向参数值。
Pi-曲面上某一点。
rij-Pi点处沿着离散方向的rij曲率。
GgU-离散弓高。
L-Pi点离散步长的一半。
L1,L2,L3,L4,L5,L6,L7,L8-分别对应曲面上经过点V1,B1,V2,M1,V3,B2,V4,M2的曲率线。
l1,l3,l5,l7-分别对应曲面上经过曲线L1,L3,L5,L7以特定规则移动后得到的经过点P1,P2,P3,P4的曲线。
具体实施方式
本发明提出一种基于曲率特征的叶片曲面统一离散方法。叶片截型线一般由叶盆部分、叶背部分和两端小的过渡圆弧组成。其中叶片曲面造型圆弧过渡的中心位置是小过渡圆弧的中点。如图1所示。本发明的流程如图7所示。
具体实施方式如下:
不妨设叶片的造型参数方向为V方向,叶片曲面的参数方程为S(u,v),其中0≤u,v≤1。
步骤一选定一个有代表性的截面线r,确定其上进排气边圆弧段的中心点B1,B2位置。具体过程如下:
选定叶片的某一条等参数线S(u,v0)作为叶片的一条截型线(如图2中标示有点的截型线所示),计算该截型线的曲率分布
其中L(u,v0)=Suu(u,v0)·n(u,v0),E(u,v0)=Su(u,v0)·Su(u,v0),n(u,v0)为曲面在点(u,v0)处的单位法矢。如图4所示,函数k(u)随着u的从小到大的会由两个突变,不妨设给定阈值系数F,定义当k(u)≥F×k0时,u处于突变位置之内。分别记录第一个突变的起始位置U1s和终止位置U1e,第二个突变的起始位置U2s和终止位置U2e。定义运算
令
则U2和U6即为圆弧过渡的中心位置。
说明:为了便于说明,本专利申请中如果不加说明,都有v0=0;std(a)运算保证了结果落在参数的定义域内;Add(a,b)则是考虑到突变区域跨越参数域边界的情况;k0实际上表示截型线的平均曲率;圆弧过渡的中心位置一般对应着截型线的最弯处;U2和U6的确定的也可以通过选取叶片多条截型线,分别通过上述步骤计算,再取平均值得到;U2和U6确定了图3中的点B1,B2。
步骤二从B1,B2点出发向叶片两端求出经过点B1,B2的两条纵向曲率线L2,L6。
具体过程如下:
令u2,1=U2,u6,1=U6,v1=0.
如果己知曲率线上第i个点,首先确定第二条曲率线下一点的参数值。计算叶片曲面在点(u2,i,vi)处与V方向相近的主方向D2,i,则下一点P2,i+1=S(u2,i,vi)+step×D2,i,其中step为给定的步长,其大小影响V方向的离散精度。再计算点P2,i+1在曲面上对应的参数值(u2,i+1,vi+1)。当vi+1≥1时,上述步骤终止。
对于j=6,确定叶片曲面在点(uj,i,vi)处与V方向相近的主方向dj,i,则
从而第j条曲率线第i+1个点的参数值为(uj,i+1,vi+1)。
对于i=2,6,记Ui对应的曲率线为Li。
说明:曲率线是通过一系列给定步长的离散点表示;也可以选取其他曲率线确定vi+1取值;计算过程保证了两条曲率线在第i个点的V参数值是一致的;主方向D2,i为三维空间中的向量,而dj,i则为曲面上一方向的比值表示。
步骤三在代表截面线r上,从B1,B2两点出发,分别在该点两侧进排气边圆弧范围内找到4个指定点V1,V2,V3,V4。
给定距离值dis1,dis2,计算U1,U3,U5,U7,使得
则U1,U3,U5,U7即为所需求的圆弧上V1,V2,V3,V4的位置。
说明:在求解积分时出现U1,U3,U4,U6与步骤一中的U1s,U1e,U2s,U2e意义不一样。实际叶片曲面在圆弧过渡边缘处的曲率走势一般会比较紊乱,所以需要通过参数dis1,dis2来调整U1,U3,U4,U6的取值,使得经过点(U1,v0),(U3,v0),(U4,v0),(U6,v0)处的四条与叶片造型方向一致的曲率线稳定可信。
步骤四从V1,V2,V3,V4出发向叶片两端求出分别经过V1,V2,V3,V4的四条纵向曲率线L1,L3,L5,L7。
u1,1=U1,u3,1=U3,u5,1=U5,u7,1=U7。
己知曲率线上第i个点,对于j=1,3,5,7,确定叶片曲面在点(uj,i,vi)处与V方向相近的主方向dj,i,则
从而第j条曲率线第i+1个点的参数值为(uj,i+1,vi+1)。
对于i=1,3,5,7,记Ui对应的曲率线为Li.
步骤五把L1,L2,L3作为等参数线,对L1,L3中间的区域进行重新参数化,确定重新参数化后的参数域边界线l1,l3;把L5,L6,L7作为等参数线,对L5,L7中间的区域进行重新参数化,确定重新参数化后的参数域边界线l5,l7。
设叶片进排气边的平均过渡圆弧弧长为ArcL,
定义
令
从而把L1,L3,L5,L7往远离进排气边方向(如图6中的箭头所示)平移至新的参数域边界。记平移后得到的曲线分别为l1,l3,l5,l7。
说明:由于造型的误差,实际叶片曲面在进排气边过渡圆弧起末端区域的曲率线走势会比较紊乱,一般不能够直接求解通过该区域的曲率线,步骤四是把该区域外侧的曲率线通过距离比例来移动,从而得到经过该区域的近似曲率线;l1,l3,l5,l7通过图3中的点P1,P2,P3,P4。
步骤六分别确定截面线r上叶盆叶背圆弧的中心点M1,M2。
令
则U4,U8记录了点M1,M2的位置。
步骤七确定经过M1,M2的两条纵向曲率线或等参数线L4,L8。
L4,L8如果为曲率线,则求解过程如下:
己知曲率线上第i个点,对于j=4,8,确定叶片曲面在点(uj,i,vi)处与V方向相近的主方向dj,i,则
从而第j条曲率线第i+1个点的参数值为(uj,i+1,vi+1)。
对于i=4,8,记Ui对应的曲率线为Li.
如果L4,L8为等参数线,则对于j=4,8,uj,i+1=Uj。
步骤八以L2,L4,L6,L8和l1,l3,l5,l7为新参数域下的等参数线,对曲面进行重新参数化。
设重新参数化的后的曲面表达式为S′(u′,v′)。
将八条线看做是新的等U′参数线,首先确定第j条线在新参数化下对应的U′坐标u′j。令u′1=0,则
u′j=u′j-1+M(Uj,Uj-1),j=2,…,8.
记u′9=1。
给定一个U′V′坐标(u′,vi),其对应的UV坐标为
其中j满足u′j≤u′≤u′j+1.
说明:八条线L2,L4,L6,L8和l1,l3,l5,l7实际上是逼近图2中叶片曲面上与V方向一致的八根线,这八根曲率线根据曲面的弯曲特征把叶片曲面分割了八块相对平稳的区域,从而可以保证每块区域中的参数线都不会出现过大的扭曲。
步骤九根据新参数对曲面进行非均匀离散。
V′方向的离散已经由步骤二确定,故只需讨论U′方向的离散。给定U′方向的离散弓高为GgU,最大离散间隔为Offset,记V′方向的离散点数为vnum。如果当前U′参数位置为U′i,对于1≤j≤vnum,计算点(U′i,vj)在UV坐标下的参数,从而确定该点处的沿U′方向的曲率uQLi,j,令
其中uLenj表示曲面上的曲线V=vj的弧长。
令offseti=min{minjOSi,j,Offset},则下一个离散位置为U′i+1=U′i+offseti。
说明:最大离散间隔Offset限制了离散的最大步长,以防止在离散叶背叶盆等较平区域时出现步长过大的情况;离散弓高的意义如图5所示;U方向的离散密度受到曲面局部曲率的影响,故为非均匀离散;V方向的离散由三维空间中的步长确定,可以近似认为为均匀离散。
本发明的实现途径为编写程序,通过CAM***的二次开发接口获取工件曲面的信息,通过上述过程对曲面进行曲率特征划分离散,离散后的曲面可以进一步通过二次开发程序进行处理。
Claims (1)
1.一种基于曲率特征的叶片曲面统一离散方法,其特征在于:它包括以下步骤:
步骤一选定一个有代表性的截面线r,确定其上进排气边圆弧段的中心点B1,B2位置;具体过程如下:
选定叶片的某一条等参数线S(u,v0)作为叶片的一条截面线,计算该截面线的曲率分布
其中L(u,v0)=S(u,v0)·n(u,v0),E(u,v0)=S(u,v0)·S(u,v0),n(u,v0)为曲面在点(u,v0)处的单位法矢,函数k(u)随着u的从小到大的会有两个突变,设给定阈值系数F,定义当k(u)≥F×k0时,u处于突变位置之内;分别记录第一个突变的起始位置U1s和终止位置U1e,第二个突变的起始位置U2s和终止位置U2e;定义运算
令
则U2和U6即为圆弧过渡的中心位置;
如果不加说明,都有v0=0;std(a)运算保证了结果落在参数的定义域内;Add(a,b)则是考虑到突变区域跨越参数域边界的情况;k0实际上表示截面线的平均曲率;圆弧过渡的中心位置对应着截面线的最弯处;U2和U6的确定通过选取叶片多条截面线,分别通过上述步骤计算,再取平均值得到;U2和U6确定了点B1,B2;
步骤二从B1,B2点出发向叶片两端求出经过点B1,B2的两条纵向曲率线L2,L6;具体过程如下:
令u2,1=U2,u6,1=U6,v1=0
如果已知曲率线上第i个点,首先确定第二条曲率线下一点的参数值;计算叶片曲面在点(u2,i,vi)处与V方向相近的主方向三维空间中的向量D2,i,则下一点P2,i+1=S(u2,i,vi)+step×D2,i,其中step为给定的步长,其大小影响V方向的离散精度;再计算点P2,i+1在曲面上对应的参数值(u2,i+1,vi+1),当vi+1≥1时,上述步骤终止;
对于j=6,确定叶片曲面在点(uj,i,vi)处与V方向相近的主方向曲面上一方向的比值dj,i,则
从而第j条曲率线第i+1个点的参数值为(uj,i+1,vi+1);
对于i=2,6,记Ui对应的曲率线为Li;
步骤三在代表截面线r上,从B1,B2两点出发,分别在B1和B2点两侧进排气边圆弧范围内找到4个指定点V1,V2,V3,V4;
给定距离值dis1,dis2,计算U1,U3,U5,U7,使得
则U1,U3,U5,U7即为所需求的圆弧上V1,V2,V3,V4的位置;
在求解积分时出现U1,U3,U4,U6与步骤一中的U1s,U1e,U2s,U2e意义不一样,实际叶片曲面在圆弧过渡边缘处的曲率走势会比较紊乱,所以需要通过参数dis1,dis2来调整U1,U3,U4,U6的取值,使得经过点(U1,v0),(U3,v0),(U4,v0),(U6,v0)处的四条与叶片造型方向一致的曲率线稳定可信;
步骤四从V1,V2,V3,V4出发向叶片两端求出分别经过V1,V2,V3,V4的四条纵向曲率线L1,L3,L5,L7;
u1,1=U1,u3,1=U3,u5,1=U5,u7,1=U7;
已知曲率线上第i个点,对于j=1,3,5,7,确定叶片曲面在点(uj,i,vi)处与V方向相近的主方向曲面上一方向的比值dj,i,则
从而第j条曲率线第i+1个点的参数值为(uj,i+1,vi+1);对于i=1,3,5,7,记Ui对应的曲率线为Li;
步骤五把L1,L2,L3作为等参数线,对L1,L3中间的区域进行重新参数化,确定重新参数化后的参数域边界线l1,l3;把L5,L6,L7作为等参数线,对L5,L7中间的区域进行重新参数化,确定重新参数化后的参数域边界线l5,l7;
设叶片进排气边的平均过渡圆弧弧长为ArcL,
定义
令
从而把L1,L3,L5,L7往远离进排气边方向平移至新的参数域;记平移后得到的曲线分别为l1,l3,l5,l7;由于造型的误差,实际叶片曲面在进排气边过渡圆弧起末端区域的曲率线走势会比较紊乱,不能够直接求解通过该区域的曲率线,步骤四是把该区域外侧的曲率线通过距离比例来移动,从而得到经过该区域的近似曲率线;
步骤六分别确定截面线r上叶盆叶背圆弧的中心点M1,M2;
令则U4,U8记录了点M1,M2的位置;
步骤七确定经过M1,M2的两条纵向曲率线或等参数线L4,L8;
L4,L8如果为曲率线,则求解过程如下:已知曲率线上第i个点,对于j=4,8,确定叶片曲面在点(uj,i,vi)处与V方向相近的主方向曲面上一方向的比值dj,i,则
从而第j条曲率线第i+1个点的参数值为(uj,i+1,vi+1);对于i=4,8,记Ui对应的曲率线为Li,如果L4,L8为等参数线,则对于j=4,8,uj,i+1=Uj;
步骤八以L2,L4,L6,L8和l1,l3,l5,l7为新参数域下的等参数线,对曲面进行重新参数化;
设重新参数化的后的曲面表达式为S′(u′,v′),将八条线看做是新的等U’参数线,首先确定第j条线在新参数化下对应的U’坐标u′j,令u′1=0,则
u′j=u′j-1+M(Uj,Uj-1),j=2,…,8.
记u′9=1;
给定一个U′V′坐标(u′,vi),其对应的UV坐标为
其中j满足u′j≤u′≤u′j+1,八条线L2,L4,L6,L8和l1,l3,l5,l7实际上是逼近叶片曲面上与V方向一致的八根线,这八根曲率线根据曲面的弯曲特征把叶片曲面分割了八块相对平稳的区域,从而保证每块区域中的参数线都不会出现过大的扭曲;
步骤九根据新参数对曲面进行非均匀离散;
V’方向的离散已经由步骤二确定,故只需讨论U’方向的离散;给定U’方向的离散弓高为GgU,最大离散间隔为Offset,记V’方向的离散点数为vnum;如果当前U’参数位置为U′i,对于1≤j≤vnum,计算点(U′i,vj)在UV坐标下的参数,从而确定(U′i,vj)处的沿U’方向的曲率uQLi,j,令
其中uLenj表示曲面上的曲线V=vj的弧长;
令offseti=min{minjOSi,j,Offset},则下一个离散位置为U′i+1=U′i+offseti。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310499607.4A CN103577636B (zh) | 2013-10-22 | 2013-10-22 | 一种基于曲率特征的叶片曲面统一离散方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310499607.4A CN103577636B (zh) | 2013-10-22 | 2013-10-22 | 一种基于曲率特征的叶片曲面统一离散方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103577636A CN103577636A (zh) | 2014-02-12 |
CN103577636B true CN103577636B (zh) | 2016-10-19 |
Family
ID=50049404
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310499607.4A Active CN103577636B (zh) | 2013-10-22 | 2013-10-22 | 一种基于曲率特征的叶片曲面统一离散方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103577636B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104439987A (zh) * | 2014-11-09 | 2015-03-25 | 沈阳黎明航空发动机(集团)有限责任公司 | 一种保证叶片进排气边圆弧形状完整性的加工方法 |
CN106372321B (zh) * | 2016-08-31 | 2018-10-02 | 中车戚墅堰机车车辆工艺研究所有限公司 | 一种摆线齿轮齿形曲线的变步长离散方法 |
CN106558056B (zh) * | 2016-11-15 | 2019-03-05 | 西北工业大学 | 基于曲率半径自适应阀值的叶片截面线前后缘分割方法 |
CN107480392B (zh) * | 2017-08-25 | 2020-06-09 | 北京航空航天大学 | 一种基于椭圆非均匀变形的叶片造型方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1714844A1 (en) * | 2004-02-10 | 2006-10-25 | Asmo Co., Ltd. | Distribution pressure setting method for wiper blade, and wiper blade |
WO2008036964A2 (en) * | 2006-09-22 | 2008-03-27 | Federal-Mogul Corporation | Contour flat blade packaging design |
CN102222138A (zh) * | 2011-05-31 | 2011-10-19 | 北京航空航天大学 | 一种基于曲面截型线分划的最短距离线对获取方法 |
CN103163836A (zh) * | 2011-12-12 | 2013-06-19 | 沈阳高精数控技术有限公司 | 直纹曲面刀具轨迹规划方法 |
-
2013
- 2013-10-22 CN CN201310499607.4A patent/CN103577636B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1714844A1 (en) * | 2004-02-10 | 2006-10-25 | Asmo Co., Ltd. | Distribution pressure setting method for wiper blade, and wiper blade |
WO2008036964A2 (en) * | 2006-09-22 | 2008-03-27 | Federal-Mogul Corporation | Contour flat blade packaging design |
CN102222138A (zh) * | 2011-05-31 | 2011-10-19 | 北京航空航天大学 | 一种基于曲面截型线分划的最短距离线对获取方法 |
CN103163836A (zh) * | 2011-12-12 | 2013-06-19 | 沈阳高精数控技术有限公司 | 直纹曲面刀具轨迹规划方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103577636A (zh) | 2014-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103577636B (zh) | 一种基于曲率特征的叶片曲面统一离散方法 | |
CN102637216B (zh) | 一种复杂组合曲面的数控侧铣加工刀轨生成方法 | |
CN104881540B (zh) | 一种涡轮气冷叶片叶身内型精确壁厚控制建模方法 | |
CN104076742B (zh) | 一种连续多类型曲线段混合插补三次多项式进给速度前瞻规划方法 | |
CN102945019B (zh) | 一种平面螺旋状和环状铣削轨迹生成方法 | |
CN106502202B (zh) | 一种球头铣刀与导向叶片接触区域的半解析建模方法 | |
CN104708824A (zh) | 一种保留模型特征的3d打印自适应切片方法 | |
CN103631198B (zh) | 基于G2连续Bézier曲线的刀具轨迹压缩方法 | |
CN105629876B (zh) | 一种三角网格模型数控加工的多刀具优化组合选择方法 | |
CN106125666A (zh) | 以切削力波动为约束的曲面加工刀具轨迹规划方法 | |
CN106354098B (zh) | 一种nurbs组合曲面上刀具加工轨迹生成方法 | |
CN103411574A (zh) | 航空发动机叶片型面三坐标测量方法 | |
CN106125672B (zh) | 一种复杂曲面零件高效加工方法 | |
CN107369127A (zh) | 一种沿行切路径进给的三维类摆线抛光轨迹生成方法 | |
CN102794488A (zh) | 类直纹面整体叶轮曲面的侧铣加工方法 | |
CN105955195A (zh) | 一种基于铣削力预测的曲面加工轨迹生成方法 | |
CN107037779B (zh) | 非均匀公差下的自由曲面nc加工刀具轨迹优化方法 | |
CN108073138A (zh) | 适用于高速高精加工的椭圆弧平滑压缩插补算法 | |
CN108469263A (zh) | 一种基于曲率进行形点优化的方法及*** | |
CN103994195A (zh) | 基于儒科夫斯基型线的液力变矩器叶片造型方法 | |
CN103577635A (zh) | 一种叶片型面数据的拟合方法 | |
CN104268317A (zh) | 机械零件圆角结构形状优化方法 | |
CN114912159A (zh) | 一种轨道交通线路平面几何线形的拟合方法 | |
CN103116679B (zh) | 航空叶片椭圆形前后缘工艺模型生成方法 | |
CN105574221A (zh) | 一种改进cst翼型参数化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |