CN111597646B - 一种五轴铣削中刀具与工件曲面求交方法与系统 - Google Patents
一种五轴铣削中刀具与工件曲面求交方法与系统 Download PDFInfo
- Publication number
- CN111597646B CN111597646B CN202010262166.6A CN202010262166A CN111597646B CN 111597646 B CN111597646 B CN 111597646B CN 202010262166 A CN202010262166 A CN 202010262166A CN 111597646 B CN111597646 B CN 111597646B
- Authority
- CN
- China
- Prior art keywords
- discrete
- intersection
- layer
- workpiece
- points
- 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
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000003801 milling Methods 0.000 title claims abstract description 26
- 238000004590 computer program Methods 0.000 claims description 4
- 238000003860 storage Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 4
- 238000003754 machining Methods 0.000 description 7
- 238000005520 cutting process Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 241001502050 Acis Species 0.000 description 1
- 238000000627 alternating current impedance spectroscopy Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000011438 discrete method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
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/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/04—Manufacturing
-
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Business, Economics & Management (AREA)
- Manufacturing & Machinery (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Human Resources & Organizations (AREA)
- General Health & Medical Sciences (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Numerical Control (AREA)
Abstract
Description
技术领域
本发明涉及铣削加工制造领域,尤其涉及一种五轴铣削中刀具与工件曲面求交方法。
背景技术
五轴立铣广泛应用于涡轮叶片、螺旋桨叶片、飞机起落架等结构和表面复杂零件的加工,以避免刀具和工件的干涉,提高加工效率和表面质量。在进行上述加工时,切削力是影响加工精度、刀具寿命和加工振动等的重要因素,其预测及准确度对工艺参数的优化具有重要意义。五轴铣削力的建模通常使用机械模型,其中,刀具与工件曲面的求交对判断刀刃微元的切削状态至关重要。
针对刀具与工件曲面的求交方法大致可以分为三类。第一类是基于ACIS或者Parasolid内核的实体模型法,比如2008年,Ferry W等人在J.Manuf.Sci.Eng.上发表文章《Cutter-Workpiece Engagement Calculations by Parallel Slicing for Five-AxisFlank Milling of Jet Engine Impellers》,通过解析几何将三维刀具-工件相交问题简化为更简单的二维平面问题集,减少了所需的数据存储和计算复杂性。第二类是离散法,典型的方法包括z-buffer、Octree、dexel以及最常用的z-map等。第三类是解析法,这种方法主要依赖于利用刀具位置(CL)文件本身逼近工件表面信息,比如2014年,Kiswanto G等人在Computer-Aided Design上发表文章《An analytical method for obtaining cutterworkpiece engagement during a semi-finish in five-axis milling》,针对平底刀和环形刀加工半精加工留下的阶梯面,分别计算下啮合点和上啮合点来计算刀具与工件表面的交线。
本发明通过将刀具和工件曲面离散为点阵来获得交线,具备比实体模型法更高的计算效率和比解析法更广的适用范围。
发明内容
本发明要解决的技术问题在于针对现有技术中计算效率低和使用范围窄的缺陷,提供一种五轴铣削中刀具与工件曲面求交方法和系统。
本发明解决其技术问题所采用的技术方案是:
提供一种五轴铣削中刀具与工件曲面求交方法,包括以下步骤:
S1、将刀具曲面沿轴向和周向离散化,形成的离散轴向高度和周向角为:
S2、对于z1,按照到周向搜索,计算刀具表面相邻离散点(i,j)和(i,j+1)的工件最近平面的法矢量夹角,如果刀具表面相邻离散点(i,j)和(i,j+1)的工件最近平面的法矢量夹角小于阈值,则(i,j)为交点,进而获得z1层的交点;
S3、从z2到沿轴向逐层搜索,第zi层搜索以zi-1层确定的交点的周向位置为基点向两边进行交点判断,直到确定该层的交点为止,再以zi层为参考进行zi+1层交点搜索,直到层离散点处理完毕,进而获得刀具表面的离散点中所有交点;
S4、将搜索到的所有交点,通过曲线拟合方法计算出交线。
接上述技术方案,步骤S2中刀具表面相邻离散点(i,j)和(i,j+1)的工件最近平面是指工件曲面点阵中距离刀具表面离散点(i,j)最近的不共线的三点构成的平面WS(i,j),和距离刀具表面离散点(i,j+1)最近的不共线的三点构成的平面WS(i,j+1),所述最近平面的判断条件为:
h(i,j)<εh (2)
其中,h(i,j)为刀具表面离散点(i,j)到其最近平面WS(i,j)的距离,εh为最近距离阈值。
接上述技术方案,步骤S2中刀具表面相邻离散点(i,j)和(i,j+1)的工件最近平面的法矢量夹角小于阈值是指:
|〈-n(i,j),n(i,j+1)〉|<εm (3)
其中,n(i,j)为刀具表面离散点(i,j)的最近平面WS(i,j)的法矢量,n(i,j+1)为刀具表面离散点(i,j+1)的最近平面WS(i,j+1)的法矢量,εm→+0为容差。
接上述技术方案,步骤S3中曲线拟合方法为最小二乘曲线拟合法。
提供一种五轴铣削中刀具与工件曲面求交系统,包括:
离散化模块,用于将刀具曲面沿轴向和周向离散化,形成的离散轴向高度和周向角为:
初始交点搜索模块,对于z1,按照到周向搜索,计算刀具表面相邻离散点(i,j)和(i,j+1)的工件最近平面的法矢量夹角,如果刀具表面相邻离散点(i,j)和(i,j+1)的工件最近平面的法矢量夹角小于阈值,则(i,j)为交点,进而获得z1层的交点;
逐层交点搜索模块,从z2到沿轴向逐层搜索,第zi层搜索以zi-1层确定的交点的周向位置为基点向两边进行交点判断,直到确定该层的交点为止,再以zi层为参考进行zi+1层交点搜索,直到层离散点处理完毕,进而获得刀具表面的离散点中所有交点;
曲线拟合模块,将搜索到的所有交点,通过曲线拟合方法计算出交线。
提供一种计算机存储介质,其内存储有可被处理器执行的计算机程序,该计算机程序执行五轴铣削中刀具与工件曲面求交方法。
提供一种所述的五轴铣削中刀具与工件曲面求交方法在五轴铣削刀具过程中的应用。
本发明产生的有益效果是:本发明提供的一种五轴铣削中刀具与工件曲面求交方法,通过将刀具曲面沿轴向和周向离散化,通过每层离散点对的工件最近平面的法矢量夹角小于阈值则说明该离散点对的起点为交点,先搜索最小轴向高度对应层的交点,再以该层交点位置为基点进行上一层的搜索,如此迭代直到所有层的交点都搜索完毕,然后对所有交点通过曲线拟合获得交线。本发明的方法计算效率更高和适用范围更广。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是本发明方法的流程图;
图2是本发明刀具曲面离散化示意图;
图3是本发明刀具表面相邻点最近平面及其法矢量示意图;
图4是本发明刀具表面与工件曲面交线示意图;
图5是本发明图4中A区域的放大示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
如图1所示,提供一种五轴铣削中刀具与工件曲面求交方法,包括以下步骤:
S1、将刀具曲面沿轴向和周向离散化,形成的离散轴向高度和周向角为:
S2、对于z1,i=1,按照到周向搜索,计算刀具表面相邻离散点(i,j)和(i,j+1)的工件最近平面的法矢量夹角,如果刀具表面相邻离散点(i,j)和(i,j+1)的工件最近平面的法矢量夹角小于阈值,则(i,j)为交点,进而获得z1层的交点;
S3、从z2到沿轴向逐层搜索,第zi层搜索以zi-1层确定的交点的周向位置为基点向两边进行交点判断,直到确定该层的交点为止,再以zi层为参考进行zi+1层交点搜索,直到层离散点处理完毕,进而获得刀具表面的离散点中所有交点;
S4、将搜索到的所有交点,通过曲线拟合方法计算出交线。
进一步地,如图3所示,步骤S2中刀具表面相邻离散点(i,j)和(i,j+1)的工件最近平面是指工件曲面点阵中距离刀具表面离散点(i,j)最近的不共线的三点构成的平面WS(i,j),和距离刀具表面离散点(i,j+1)最近的不共线的三点构成的平面WS(i,j+1),所述最近平面的判断条件为:
h(i,j)<εh (2)
其中,h(i,j)为刀具表面离散点(i,j)到其最近平面WS(i,j)的距离,εh为最近距离阈值。
进一步地,步骤S2中刀具表面相邻离散点(i,j)和(i,j+1)的工件最近平面的法矢量夹角小于阈值是指:
|〈-n(i,j),n(i,j+1)〉|<εm (3)
其中,n(i,j)为刀具表面离散点(i,j)的最近平面WS(i,j)的法矢量,n(i,j+1)为刀具表面离散点(i,j+1)的最近平面WS(i,j+1)的法矢量,εm→+0为容差。
进一步地,步骤S3中曲线拟合方法为最小二乘曲线拟合法。如图4所示为利用该方法获得刀具与工件曲面交线示意图,其中,TS为刀具表面,WS为工件表面。图5为图4中A区域的放大图。曲面交线的精度随着刀具表面划分密度和工件表面点阵密度的增大而提升,但密度过大会迅速增加计算量,所以实际使用时应当在满足精度要求时尽量减小刀具表面划分密度和工件表面的密度。应当注意的是,在提升求交精度时,应同步增大刀具表面划分密度和工件表面的点阵密度,两者不可差异过大,否则会造成精度下降。
提供一种五轴铣削中刀具与工件曲面求交系统,包括:
离散化模块,用于将刀具曲面沿轴向和周向离散化,形成的离散轴向高度和周向角为:
初始交点搜索模块,对于z1,按照到周向搜索,计算刀具表面相邻离散点(i,j)和(i,j+1)的工件最近平面的法矢量夹角,如果刀具表面相邻离散点(i,j)和(i,j+1)的工件最近平面的法矢量夹角小于阈值,则(i,j)为交点,进而获得z1层的交点;
逐层交点搜索模块,从z2到沿轴向逐层搜索,第zi层搜索以zi-1层确定的交点的周向位置为基点向两边进行交点判断,直到确定该层的交点为止,再以zi层为参考进行zi+1层交点搜索,直到层离散点处理完毕,进而获得刀具表面的离散点中所有交点;
曲线拟合模块,将搜索到的所有交点,通过曲线拟合方法计算出交线。
提供一种计算机存储介质,其内存储有可被处理器执行的计算机程序,该计算机程序执行五轴铣削中刀具与工件曲面求交方法。
提供一种所述的五轴铣削中刀具与工件曲面求交方法在五轴铣削刀具过程中的应用。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (6)
1.一种五轴铣削中刀具与工件曲面求交方法,其特征在于,包括以下步骤:
S1、将刀具曲面沿轴向和周向离散化,形成的离散轴向高度和周向角为:
S2、对于z1,按照到周向搜索,计算刀具表面相邻离散点(i,j)和(i,j+1)的工件最近平面的法矢量夹角,如果刀具表面相邻离散点(i,j)和(i,j+1)的工件最近平面的法矢量夹角小于阈值,则(i,j)为交点,进而获得z1层的交点;
S3、从z2到沿轴向逐层搜索,第zi层搜索以zi-1层确定的交点的周向位置为基点向两边进行交点判断,直到确定该层的交点为止,再以zi层为参考进行zi+1层交点搜索,直到层离散点处理完毕,进而获得刀具表面的离散点中所有交点;
S4、将搜索到的所有交点,通过曲线拟合方法计算出交线;
其中,步骤S2中刀具表面相邻离散点(i,j)和(i,j+1)的工件最近平面是指工件曲面点阵中距离刀具表面离散点(i,j)最近的不共线的三点构成的平面WS(i,j),和距离刀具表面离散点(i,j+1)最近的不共线的三点构成的平面WS(i,j+1),所述最近平面的判断条件为:
h(i,j)<εh (2)
其中,h(i,j)为刀具表面离散点(i,j)到其最近平面WS(i,j)的距离,εh为最近距离阈值。
2.根据权利要求1所述的五轴铣削中刀具与工件曲面求交方法,其特征在于,步骤S2中刀具表面相邻离散点(i,j)和(i,j+1)的工件最近平面的法矢量夹角小于阈值是指:
|<-n(i,j),n(i,j+1)>|<εm (3)
其中,n(i,j)为刀具表面离散点(i,j)的最近平面WS(i,j)的法矢量,n(i,j+1)为刀具表面离散点(i,j+1)的最近平面WS(i,j+1)的法矢量,εm→+0为容差。
3.根据权利要求1所述的五轴铣削中刀具与工件曲面求交方法,其特征在于,步骤S3中曲线拟合方法为最小二乘曲线拟合法。
4.一种五轴铣削中刀具与工件曲面求交系统,其特征在于,包括:
离散化模块,用于将刀具曲面沿轴向和周向离散化,形成的离散轴向高度和周向角为:
初始交点搜索模块,用于对于z1,按照到周向搜索,计算刀具表面相邻离散点(i,j)和(i,j+1)的工件最近平面的法矢量夹角,如果刀具表面相邻离散点(i,j)和(i,j+1)的工件最近平面的法矢量夹角小于阈值,则(i,j)为交点,进而获得z1层的交点;其中,工件最近平面是指工件曲面点阵中距离刀具表面离散点(i,j)最近的不共线的三点构成的平面WS(i,j),和距离刀具表面离散点(i,j+1)最近的不共线的三点构成的平面WS(i,j+1),所述最近平面的判断条件为:
h(i,j)<εh (2)
其中,h(i,j)为刀具表面离散点(i,j)到其最近平面WS(i,j)的距离,εh为最近距离阈值;
逐层交点搜索模块,用于从z2到沿轴向逐层搜索,第zi层搜索以zi-1层确定的交点的周向位置为基点向两边进行交点判断,直到确定该层的交点为止,再以zi层为参考进行zi+1层交点搜索,直到层离散点处理完毕,进而获得刀具表面的离散点中所有交点;
曲线拟合模块,用于将搜索到的所有交点,通过曲线拟合方法计算出交线。
5.一种计算机存储介质,其特征在于,其内存储有可被处理器执行的计算机程序,该计算机程序执行如权利要求1-3中任一项所述的五轴铣削中刀具与工件曲面求交方法。
6.一种如权利要求1-3中任一项所述的五轴铣削中刀具与工件曲面求交方法在五轴铣削刀具过程中的应用方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010262166.6A CN111597646B (zh) | 2020-04-06 | 2020-04-06 | 一种五轴铣削中刀具与工件曲面求交方法与系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010262166.6A CN111597646B (zh) | 2020-04-06 | 2020-04-06 | 一种五轴铣削中刀具与工件曲面求交方法与系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111597646A CN111597646A (zh) | 2020-08-28 |
CN111597646B true CN111597646B (zh) | 2022-05-20 |
Family
ID=72187318
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010262166.6A Active CN111597646B (zh) | 2020-04-06 | 2020-04-06 | 一种五轴铣削中刀具与工件曲面求交方法与系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111597646B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101866162A (zh) * | 2010-06-12 | 2010-10-20 | 北京航空航天大学 | 点到曲面距离计算的邻近三角形方法 |
CN107544433A (zh) * | 2017-09-21 | 2018-01-05 | 大连理工大学 | 五轴数控机床平面加工过程球头铣刀与工件接触区域半解析建模方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI414376B (zh) * | 2011-04-14 | 2013-11-11 | Nat Univ Tsing Hua | 一種五軸曲面側銑加工系統及其刀具路徑規劃方法 |
-
2020
- 2020-04-06 CN CN202010262166.6A patent/CN111597646B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101866162A (zh) * | 2010-06-12 | 2010-10-20 | 北京航空航天大学 | 点到曲面距离计算的邻近三角形方法 |
CN107544433A (zh) * | 2017-09-21 | 2018-01-05 | 大连理工大学 | 五轴数控机床平面加工过程球头铣刀与工件接触区域半解析建模方法 |
Non-Patent Citations (4)
Title |
---|
《Influence mechanism of machining angles on force induced error and their selection in five axis bullnose end milling》;Zerun ZHU等;《Chinese Journal of Aeronautics》;20200314;第3447-3459页 * |
刀具扫描体与曲面法矢求交新算法;王大镇等;《华侨大学学报(自然科学版)》;20050720(第03期);第295-298页 * |
叶轮数控加工曲面法矢与刀具求交研究;王大镇;《陕西工学院学报》;20030630(第02期);第19-21页 * |
复杂曲面加工精度检验中曲面法矢与刀具扫描体求交算法的研究;李吉平等;《计算机辅助设计与图形学学报》;20001201(第12期);第931-935页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111597646A (zh) | 2020-08-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Altintas et al. | Virtual process systems for part machining operations | |
Can et al. | A novel iso-scallop tool-path generation for efficient five-axis machining of free-form surfaces | |
CN109343468B (zh) | 一种基于投影偏置的叶片多轴轨迹生成方法 | |
JP2804829B2 (ja) | 曲面加工方法 | |
CN107451382B (zh) | 高速切削加工工件表面形貌的控制方法 | |
Zhu et al. | High efficiency simulation of five-axis cutting force based on the symbolically solvable cutting contact boundary model | |
Ren et al. | Clean-up tool path generation by contraction tool method for machining complex polyhedral models | |
Scippa et al. | Milled surface generation model for chip thickness detection in peripheral milling | |
Ju et al. | A solid-analytical-based method for extracting cutter-workpiece engagement in sculptured surface milling | |
Huo et al. | CNC tool path generation for freeform surface machining based on preferred feed direction field | |
Wei et al. | Force predictive model for five-axis ball end milling of sculptured surface | |
Li et al. | A voxel model-based process-planning method for five-axis machining of complicated parts | |
CN111324955A (zh) | 一种自由曲面激光切削加工的方法 | |
Zhang et al. | Modeling, simulation, and optimization of five-axis milling processes | |
CN111597646B (zh) | 一种五轴铣削中刀具与工件曲面求交方法与系统 | |
Ning et al. | Parametric design and surface topography analysis of turbine blade processing by turn-milling based on CAM | |
He et al. | Geodesic distance field-based process planning for five-axis machining of complicated parts | |
Si et al. | A solid-discrete-based method for extracting the cutter-workpiece engagement in five-axis flank milling | |
Tsainis et al. | A novel convex hull method for optimum multi-point 5-axis tool positioning for machining of complex sculptured surfaces | |
Ferry | Virtual five-axis flank milling of jet engine impellers | |
CN111061218B (zh) | 一种基于acis的复杂曲面加工切削力预测方法及设备 | |
Rao et al. | Integrated tool positioning and tool path planning for five-axis machining of sculptured surfaces | |
Xu et al. | Non-redundant tool trajectory generation for surface finish machining based on geodesic curvature matching | |
Ichikawa et al. | Estimation Method of Machining Error on Low Rigidity Workpiece for Tool Posture Planning | |
Sutrisno et al. | Development of initial tool orientation method at close bounded area for 5-axis roughing based on faceted models |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |