CN107085413A - 基于ac型五轴数控机床球头刀加工刀轴矢量光顺方法 - Google Patents

基于ac型五轴数控机床球头刀加工刀轴矢量光顺方法 Download PDF

Info

Publication number
CN107085413A
CN107085413A CN201710280231.6A CN201710280231A CN107085413A CN 107085413 A CN107085413 A CN 107085413A CN 201710280231 A CN201710280231 A CN 201710280231A CN 107085413 A CN107085413 A CN 107085413A
Authority
CN
China
Prior art keywords
msub
mrow
mtd
mtr
beta
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.)
Granted
Application number
CN201710280231.6A
Other languages
English (en)
Other versions
CN107085413B (zh
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.)
Shandong University of Technology
Original Assignee
Shandong University of Technology
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 Shandong University of Technology filed Critical Shandong University of Technology
Priority to CN201710280231.6A priority Critical patent/CN107085413B/zh
Publication of CN107085413A publication Critical patent/CN107085413A/zh
Application granted granted Critical
Publication of CN107085413B publication Critical patent/CN107085413B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/408Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by data handling or data format, e.g. reading, buffering or conversion of data
    • G05B19/4086Coordinate conversions; Other special calculations
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/33Director till display
    • G05B2219/33263Conversion, transformation of coordinates, cartesian or polar

Abstract

本发明公开了一种基于AC型五轴数控机床球头刀加工刀轴矢量光顺方法,包括如下步骤:建立球头刀刀轴矢量与刀位设计变量之间的关系方程;建立球头刀刀轴矢量与五轴数控机床回转轴A和C之间的运动变换方程;建立球头刀刀位设计变量和五轴数控机床回转轴A和C之间的关系方程;确定球头刀刀轴矢量光顺的设计变量、目标函数和约束条件,建立基于AC型五轴数控机床球头刀加工刀轴矢量光顺数学模型;确定上述刀轴矢量光顺数学模型的求解方法。该方法能够避免机床回转轴的急剧变化,使机床回转轴运动更加平稳和光滑,大幅降低机床回转轴的角速度和角加速度,从而提高曲面的加工质量和加工效率,具有较强的实际应用价值。

Description

基于AC型五轴数控机床球头刀加工刀轴矢量光顺方法
技术领域
本发明涉及一种五轴数控机床刀轴矢量光顺方法,尤其涉及基于AC型五轴数控机床球头刀加工刀轴矢量光顺方法,属于五轴数控加工技术领域。
背景技术
当使用球头刀五轴加工复杂曲面时,由于曲面几何性质较差,比如曲面的法向量、主方向、曲率等,容易导致所生成的球头刀刀轴矢量发生突变和大幅波动。即使采用最简单刀具定位方法(如Sturz法)五轴加工上述曲面区域,也会引起刀轴矢量的剧烈变化,从而影响五轴机床运动的平稳性、超出机床进给轴的伺服能力和增大加工过程中的非线性误差等。因此获得五轴数控加工中球头刀光滑刀轴矢量成为曲面加工技术的重要研究方向。为获得光滑的刀轴矢量,国内外学者在刀轴矢量优化方面开展了大量研究工作,并提出许多五轴加工刀轴矢量优化方法,主要集中于两个方面:一是在工件坐标系中仅考虑几何约束的刀轴矢量优化方法;二是在工件坐标系中考虑几何约束和运动学约束的刀轴矢量优化方法。
现有技术一,文献(周波,赵吉宾与刘伟军,复杂曲面五轴数控加工刀轴矢量优化方法研究.机械工程学报,2013(07):184-192)提出一种复杂曲面五轴数控加工刀轴矢量优化方法。首先在非干涉域内有效插入限定的加工点位,保证刀轴矢量的整体优化;同时在干涉域内采用改进的C-Space法,生成刀轴矢量光顺可行域。
现有技术二,文献(王晶等,复杂曲面零件五轴加工刀轴整体优化方法.航空学报,2013(06):1452-1462)提出一种基于临界约束的五轴刀轴矢量整体优化方法。首先构造了给定切触点处所有可行摆刀平面,并在摆刀平面内根据临界约束计算出临界刀轴矢量,在获得临界刀轴矢量的基础上,对其进行平面映射,建立了刀轴摆动的初始可行域;其次,通过对初始可行域进行均匀离散,根据离散点之间相对位置关系构造邻接矩阵,并结合最短路径搜索算法获得了初始参考刀轴,从而构造了新的刀轴摆动可行域;最后,建立当前切削行内无干涉且相邻刀轴变化最小的刀轴矢量优化模型,实现自由曲面五轴加工无干涉刀轴矢量的光滑控制。上述现有技术至少存在以下缺点:
上述方法主要以刀具与工件/机床之间不发生过切和全局干涉等作为约束条件,以工件坐标系中刀轴矢量变化最小或光滑变化作为目标函数进行刀轴矢量优化,从而获得工件坐标系中无干涉和光滑的刀轴矢量。而五轴数控机床因结构差异较大,尽管工件坐标系中刀轴矢量光滑变化,却很难保证机床坐标系中五轴机床各坐标轴尤其是回转轴都能光滑运动而无突变现象发生,从而影响机床运动的平稳性、超出机床进给轴的伺服能力和增大加工过程中的非线性误差等。因此,有必要在工件坐标系中进行刀轴矢量优化时需考虑机床坐标系中机床回转轴的变化情况,或者在机床坐标系中直接对机床回转轴进行光顺优化。
发明内容
为克服现有球头刀五轴加工复杂曲面刀轴矢量发生突变和大幅波动的问题,本发明提供一种基于AC型五轴数控机床球头刀加工刀轴矢量光顺方法。
为了实现上述目的,本发明采用的技术方案是这样的:一种基于AC型五轴数控机床球头刀加工刀轴矢量光顺方法,包括以下步骤:
A、建立球头刀刀轴矢量与刀位设计变量之间的关系方程;
B、建立球头刀刀轴矢量与五轴数控机床回转轴A和C之间的运动变换方程;
C、建立球头刀刀位设计变量和五轴数控机床回转轴A和C之间的关系方程;
D、确定球头刀刀轴矢量光顺的设计变量、目标函数和约束条件,建立基于AC型五轴数控机床球头刀加工刀轴矢量光顺数学模型;
E、确定步骤D中刀轴矢量光顺数学模型的求解方法。
作为优选:所述步骤A具体为:
(1)在刀触点处建立刀具局部坐标系,推导出局部坐标系OLXLYLZL中球头刀刀轴矢量与刀位设计变量之间的关系方程:
式中,θ为局部坐标系中球头刀绕YL轴转动的前倾角,ψ为局部坐标系中球头刀具绕XL轴转动的侧倾角,
在局部坐标系OLXLYLZL中球头刀刀触点处刀位点径矢
(2)建立工件坐标系中球头刀刀轴矢量与刀位设计变量之间的关系方程:
式中,e1=(x1,y1,z1)T,e2=(x2,y2,z2)T,e3=(x3,y3,z3)T分别为局部坐标系OLXLYLZL各坐标轴在工件坐标系OwXwYwZw中的单位矢量;
在工件坐标系OwXwYwZw中球头刀刀触点处刀位点径矢
作为优选:所述步骤B中,根据AC型五轴数控机床具体结构、机床运动链和机床各坐标系之间关系,建立球头刀刀轴矢量与五轴数控机床回转轴A和C之间关系方程:
作为优选:所述步骤C中,联立式(3)和(5),则推导出球头刀刀位设计变量和机床回转轴A和C之间的关系方程:
作为优选:所述步骤D中,以机床坐标系中机床回转轴A和C作为设计变量,以曲面每行刀轨所有刀触点{Pi,i=1,…,N}处复合角加速度的平方和作为目标函数,以机床回转轴A和C的角度、角速度和角加速度限制范围作为约束条件,建立基于AC型五轴数控机床球头刀加工刀轴矢量光顺数学模型:
式中,N为给定行刀轨上刀触点数目,β1和β2分别代表机床回转轴A和C转角,ω1和ω2分别代表机床回转轴A和C角速度,α1和α2分别代表机床回转轴A和C加速度,分别代表机床回转角β1和β2运动范围,分别代表机床回转角β1和β2角速度限制范围,分别代表机床回转角β1和β2角加速度限制范围,曲面上任意刀触点Pi处的复合角加速度定义为:
式中,为刀触点Pi处的刀轴矢量,t代表时间。
作为优选:所述步骤E具体为:
(1)利用现有球头刀五轴刀位优化方法生成给定曲面采样刀触点处的刀位,然后获得上述采样刀触点处的初始刀轴矢量;
(2)利用式(5)计算出各采样刀触点处所对应的机床回转角A和C,再利用三次样条插值函数分别将各采样刀触点处的机床回转角A和C进行三次样条拟合;
(3)利用步骤(2)所得到的三次样条插值拟合函数计算每行刀轨刀触曲线上各刀触点处机床回转角A和C,然后利用式(6)计算各刀触点处球头刀刀位设计变量θ和ψ,再利用式(4)和(3)计算出各刀触点处球头刀刀位点径矢和刀轴矢量,直到求解出该行刀轨刀触曲线上所有刀触点的刀位点径矢和刀轴矢量。
本发明的有益效果是,该方法能够避免机床回转轴突变和不光顺,使机床回转轴运动更加平稳和光滑,大幅减少机床回转轴的角速度和角加速度,从而提高曲面的加工质量和加工效率。
附图说明
图1为基于AC型五轴数控机床球头刀加工刀轴矢量光顺方法流程图
图2为球头刀初步定位;
图3为AC′型五轴数控机床结构示意图;
图4为AC′型五轴数控机床中的坐标系;
图5为基于AC型五轴数控机床球头刀加工刀轴矢量光顺数学模型求解流程图。
具体实施方式
本发明的一种基于AC型五轴数控机床球头刀加工刀轴矢量光顺方法,其基本流程如图1所示,较佳的具体实施方式是,包括:
步骤A、建立球头刀刀轴矢量与刀位设计变量之间的关系方程。所述步骤A具体为:
(1)局部坐标系中球头刀刀轴矢量与刀位设计变量之间的关系方程
如图1所示,设球头刀加工复杂曲面S:r(u,v),Pcc(ucc,vcc)为曲面上任意一点,ncc为该点的单位法矢,OwXwYwZw为工件坐标系。令r为球头刀的半径,OL为局部坐标系的原点,其中和rPcc分别为点OL和Pcc对应的径矢。分别以OL为原点建立点Pcc处的局部坐标系OLXLYLZL和刀具坐标系OtXtYtZt,且点Ot与OL重合。假设初始状态时OtXtYtZt与OLXLYLZL的各坐标轴方向一致,则刀具在点Pcc处被初步定位。在局部坐标系OLXLYLZL中,刀具存在两个自由度:一是绕YL轴旋转的前倾角θ,二是绕XL轴旋转的侧倾角ψ,以上两角即为刀位设计变量。
通过调整上述两个角度可以获得不同的刀位,可以得到局部坐标系OLXLYLZL中刀触点Pcc处的刀轴矢量和刀位点径矢
式中,
由式(1)可得
由式(2)可得
(2)工件坐标系中球头刀具刀轴矢量与刀位设计变量之间的关系方程假设局部坐标系OLXLYLZL各坐标轴在工件坐标系OwXwYwZw中的矢量分别为e1=(x1,y1,z1)T,e2=(x2,y2,z2)T,e3=(x3,y3,z3)T,则在工件坐标系OwXwYwZw中刀轴矢量和刀位点径矢为
由式(5)可得到球头刀刀位设计变量(即前倾角θ和侧倾角ψ)和刀轴矢量之间的关系方程
步骤B、建立球头刀刀轴矢量与五轴数控机床回转轴A和C之间的运动变换方程。所述步骤B具体为:
根据回转轴所在位置不同,AC型五轴数控机床又可以分为三种类型:1)AC型(A和C轴都位于主轴侧),2)AC′型(A轴位于主轴侧和C轴位于转台侧),3)A′C′型(A和C轴都位于转台侧)。为了简单起见,下文将以AC′型五轴数控机床作为研究对象进行阐述,其推导公式适用于所有AC型五轴数控机床,如图2所示。图3所示为AC′型五轴数控机床中各坐标系之间关系,一般默认工件坐标系OwXwYwZw与机床坐标系OmXmYmZm的坐标轴方向一致,数控系统可以将两者进行关联。由此可得,球头刀刀轴矢量与五轴数控机床回转轴A和C之间关系方程:
式中,
重新整理式(8)可得
步骤C、建立球头刀刀位设计变量和五轴数控机床回转轴A和C之间的关系方程。所述步骤C具体为:
将式(9)代入式(7)便得到球头刀刀位设计变量和五轴数控机床回转轴A和C之间的关系方程:
步骤D、确定球头刀刀轴矢量光顺的设计变量、目标函数和约束条件,建立基于AC型五轴数控机床球头刀加工刀轴矢量光顺数学模型。所述步骤D具体为:
(1)五轴加工中复合角加速度的定义
在机床坐标系中引入复合角加速度的概念,在设计曲面上任意刀触点Pi处复合角加速度的定义为
式中,β1和β2分别代表机床回转轴A和C转角,为刀触点Pi处的刀轴矢量,t代表时间,ω1和ω2分别代表机床回转轴A和C角速度,α1和α2分别代表机床回转轴A和C加速度。
(2)构建基于AC型五轴数控机床球头刀加工刀轴矢量光顺数学模型
以机床坐标系中机床回转轴A和C作为设计变量,以曲面每行刀轨所有刀触点{Pi,i=1,…,N}处复合角加速度的平方和作为目标函数,以机床回转轴A和C的角度、角速度和角加速度限制范围作为约束条件,建立基于AC型五轴数控机床球头刀加工刀轴矢量光顺数学模型:
式中,分别代表机床回转角β1和β2运动范围,分别代表机床回转角β1和β2角速度设置范围,分别代表机床回转角β1和β2角加速度设置范围。
步骤E、确定步骤D中刀轴矢量光顺数学模型的求解方法。所述步骤E具体为:提出一种求解步骤D中刀轴矢量光顺的数学模型的方法,求解过程如图4所示,其具体过程如下:
(1)利用现有球头刀五轴刀位优化方法生成给定曲面采样刀触点处的刀位,然后获得上述采样刀触点处的初始刀轴矢量。
首先,通过分析加工曲面的特点,初步给定球头刀刀位设计变量的数值,即前倾角θ和侧倾角ψ,假设给定刀轨刀触曲线上采样刀触点数为N,实际刀触点数为M,给定步长公差为h。然后,在给定刀轨刀触曲线上利用等参数离散法生成采样刀触点{CCi,i=1,…,N};最后,利用Sturz法计算所有采样刀触点CCi处刀轴矢量{Tai,i=1,…,N}和刀位点径矢{Tpi,i=1,…,N}。
(2)利用式(10)计算出各采样刀触点CCi处所对应的机床回转角Ai和Ci,再利用三次样条插值函数将各采样刀触点CCi处的机床回转角Ai和Ci进行三次样条拟合,从而得到机床回转角Ai和Ci的三次样条插值拟合函数F(u,A)和F(u,C),其中u为采样刀触点的参数变量。
(3)对于给定刀轨刀触曲线,由给定M或h并利用等参数离散法或等弦高误差离散法生成实际刀触点{PCi,i=1,…,M},再利用步骤(2)中三次样条插值拟合函数F(u,A)和F(u,C)计算该行刀轨刀触曲线上实际刀触点PCi处所对应的机床回转角Ai和Ci;然后利用式(10)得到实际刀触点PCi处机床回转角Ai和Ci所对应的球头刀刀位设计变量θi和ψi,进而由式(6)和(7)得到该刀触点PCi处球头刀的刀位点径矢Tpi和刀轴矢量Tai,直到i=M求解出该行刀触曲线上所有刀触点的刀位点和刀轴矢量。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。

Claims (6)

1.一种基于AC型五轴数控机床球头刀加工刀轴矢量光顺方法,其特征在于,包括如下步骤:
A、建立球头刀刀轴矢量与刀位设计变量之间的关系方程;
B、建立球头刀刀轴矢量与五轴数控机床回转轴A和C之间的运动变换方程;
C、建立球头刀刀位设计变量和五轴数控机床回转轴A和C之间的关系方程;
D、确定球头刀刀轴矢量光顺的设计变量、目标函数和约束条件,建立基于AC型五轴数控机床球头刀加工刀轴矢量光顺数学模型;
E、确定步骤D中刀轴矢量光顺数学模型的求解方法。
2.根据权利要求1所述的基于AC型五轴数控机床球头刀加工刀轴矢量光顺方法,其特征在于,所述步骤A具体为:
(1)在刀触点处建立刀具局部坐标系,推导出局部坐标系OLXLYLZL中球头刀刀轴矢量与刀位设计变量之间的关系方程:
<mrow> <msubsup> <mi>T</mi> <mrow> <mi>a</mi> <mi>x</mi> <mi>i</mi> <mi>s</mi> </mrow> <mrow> <mi>L</mi> <mi>C</mi> <mi>S</mi> </mrow> </msubsup> <mo>=</mo> <mi>R</mi> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mi>L</mi> </msub> <mo>,</mo> <mi>&amp;psi;</mi> <mo>)</mo> </mrow> <mi>R</mi> <mrow> <mo>(</mo> <msub> <mi>Y</mi> <mi>L</mi> </msub> <mo>,</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <mn>0</mn> <mo>,</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>T</mi> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
式中,θ为局部坐标系中球头刀绕YL轴转动的前倾角,ψ为局部坐标系中球头刀具绕XL轴转动的侧倾角,
在局部坐标系OLXLYLZL中球头刀刀触点处刀位点径矢
<mrow> <msubsup> <mi>T</mi> <mrow> <mi>p</mi> <mi>o</mi> <mi>s</mi> </mrow> <mrow> <mi>L</mi> <mi>C</mi> <mi>S</mi> </mrow> </msubsup> <mo>=</mo> <mi>R</mi> <mrow> <mo>(</mo> <msub> <mi>X</mi> <mi>L</mi> </msub> <mo>,</mo> <mi>&amp;psi;</mi> <mo>)</mo> </mrow> <mi>R</mi> <mrow> <mo>(</mo> <msub> <mi>Y</mi> <mi>L</mi> </msub> <mo>,</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <mn>0</mn> <mo>,</mo> <mn>0</mn> <mo>,</mo> <mo>-</mo> <mi>r</mi> <mo>)</mo> </mrow> <mi>T</mi> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
(2)建立工件坐标系中球头刀刀轴矢量与刀位设计变量之间的关系方程:
<mrow> <msubsup> <mi>T</mi> <mrow> <mi>a</mi> <mi>x</mi> <mi>i</mi> <mi>s</mi> </mrow> <mrow> <mi>W</mi> <mi>C</mi> <mi>S</mi> </mrow> </msubsup> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msub> <mi>e</mi> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mi>e</mi> <mn>2</mn> </msub> </mtd> <mtd> <msub> <mi>e</mi> <mn>3</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <msubsup> <mi>T</mi> <mrow> <mi>a</mi> <mi>x</mi> <mi>i</mi> <mi>s</mi> </mrow> <mrow> <mi>L</mi> <mi>C</mi> <mi>S</mi> </mrow> </msubsup> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mn>1</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;theta;</mi> <mo>-</mo> <msub> <mi>x</mi> <mn>2</mn> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;psi;</mi> <mo>+</mo> <msub> <mi>x</mi> <mn>3</mn> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;psi;</mi> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mn>1</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;theta;</mi> <mo>-</mo> <msub> <mi>y</mi> <mn>2</mn> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;psi;</mi> <mo>+</mo> <msub> <mi>y</mi> <mn>3</mn> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;psi;</mi> </mtd> </mtr> <mtr> <mtd> <msub> <mi>z</mi> <mn>1</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;theta;</mi> <mo>-</mo> <msub> <mi>z</mi> <mn>2</mn> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;psi;</mi> <mo>+</mo> <msub> <mi>z</mi> <mn>3</mn> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;psi;</mi> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
式中,e1=(x1,y1,z1)T,e2=(x2,y2,z2)T,e3=(x3,y3,z3)T分别为局部坐标系OLXLYLZL各坐标轴在工件坐标系OwXwYwZw中的单位矢量;
在工件坐标系OwXwYwZw中球头刀刀触点处刀位点径矢
<mrow> <msubsup> <mi>T</mi> <mrow> <mi>p</mi> <mi>o</mi> <mi>s</mi> </mrow> <mrow> <mi>W</mi> <mi>C</mi> <mi>S</mi> </mrow> </msubsup> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msub> <mi>e</mi> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mi>e</mi> <mn>2</mn> </msub> </mtd> <mtd> <msub> <mi>e</mi> <mn>3</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <msubsup> <mi>T</mi> <mrow> <mi>p</mi> <mi>o</mi> <mi>s</mi> </mrow> <mrow> <mi>L</mi> <mi>C</mi> <mi>S</mi> </mrow> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
3.根据权利要求1所述的基于AC型五轴数控机床球头刀加工刀轴矢量光顺方法,其特征在于,所述步骤B中,根据AC型五轴数控机床具体结构、机床运动链和机床各坐标系之间关系,建立球头刀刀轴矢量与五轴数控机床回转轴A和C之间关系方程:
<mrow> <msubsup> <mi>T</mi> <mrow> <mi>a</mi> <mi>x</mi> <mi>i</mi> <mi>s</mi> </mrow> <mrow> <mi>W</mi> <mi>C</mi> <mi>S</mi> </mrow> </msubsup> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mi>sin</mi> <mi> </mi> <mi>A</mi> <mi> </mi> <mi>sin</mi> <mi> </mi> <mi>C</mi> <mo>,</mo> <mo>-</mo> <mi>sin</mi> <mi> </mi> <mi>A</mi> <mi> </mi> <mi>cos</mi> <mi> </mi> <mi>C</mi> <mo>,</mo> <mi>cos</mi> <mi> </mi> <mi>A</mi> <mo>)</mo> </mrow> <mi>T</mi> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
4.根据权利要求1所述的基于AC型五轴数控机床球头刀加工刀轴矢量光顺方法,其特征在于,所述步骤C中,联立式(3)和(5),则推导出球头刀刀位设计变量和机床回转轴A和C之间的关系方程:
<mrow> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mi>sin</mi> <mi> </mi> <mi>A</mi> <mi> </mi> <mi>sin</mi> <mi> </mi> <mi>C</mi> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <mi>sin</mi> <mi> </mi> <mi>A</mi> <mi> </mi> <mi>cos</mi> <mi> </mi> <mi>C</mi> </mtd> </mtr> <mtr> <mtd> <mi>cos</mi> <mi> </mi> <mi>A</mi> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mn>1</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;theta;</mi> <mo>-</mo> <msub> <mi>x</mi> <mn>2</mn> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;psi;</mi> <mo>+</mo> <msub> <mi>x</mi> <mn>3</mn> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;psi;</mi> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mn>1</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;theta;</mi> <mo>-</mo> <msub> <mi>y</mi> <mn>2</mn> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;psi;</mi> <mo>+</mo> <msub> <mi>y</mi> <mn>3</mn> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;psi;</mi> </mtd> </mtr> <mtr> <mtd> <msub> <mi>z</mi> <mn>1</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;theta;</mi> <mo>-</mo> <msub> <mi>z</mi> <mn>2</mn> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;psi;</mi> <mo>+</mo> <msub> <mi>z</mi> <mn>3</mn> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;psi;</mi> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
5.根据权利要求1所述的基于AC型五轴数控机床球头刀加工刀轴矢量光顺方法,其特征在于,所述步骤D中,以机床坐标系中机床回转轴A和C作为设计变量,以曲面每行刀轨所有刀触点{Pi,i=1,…,N}处复合角加速度的平方和作为目标函数,以机床回转轴A和C的角度、角速度和角加速度限制范围作为约束条件,建立基于AC型五轴数控机床球头刀加工刀轴矢量光顺数学模型:
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>min</mi> <mi>&amp;Gamma;</mi> <mo>=</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mo>|</mo> <mo>|</mo> <msub> <mi>&amp;alpha;</mi> <msub> <mi>P</mi> <mi>i</mi> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>s</mi> <mo>.</mo> <mi>t</mi> <mo>.</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mover> <mi>&amp;beta;</mi> <mo>^</mo> </mover> <mrow> <mn>1</mn> <mi>lim</mi> <mi>i</mi> <mi>t</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mover> <mi>&amp;beta;</mi> <mo>^</mo> </mover> <mrow> <mn>2</mn> <mi>lim</mi> <mi>i</mi> <mi>t</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;omega;</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mover> <mi>&amp;omega;</mi> <mo>^</mo> </mover> <mrow> <mn>1</mn> <mi>lim</mi> <mi>i</mi> <mi>t</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;omega;</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mover> <mi>&amp;omega;</mi> <mo>^</mo> </mover> <mrow> <mn>2</mn> <mi>lim</mi> <mi>i</mi> <mi>t</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mover> <mi>&amp;alpha;</mi> <mo>^</mo> </mover> <mrow> <mn>1</mn> <mi>lim</mi> <mi>i</mi> <mi>t</mi> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>=</mo> <msub> <mover> <mi>&amp;alpha;</mi> <mo>^</mo> </mover> <mrow> <mn>2</mn> <mi>lim</mi> <mi>i</mi> <mi>t</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
式中,N为给定行刀轨上刀触点数目,β1和β2分别代表机床回转轴A和C转角,ω1和ω2分别代表机床回转轴A和C角速度,α1和α2分别代表机床回转轴A和C加速度,分别代表机床回转角β1和β2运动范围,分别代表机床回转角β1和β2角速度限制范围,分别代表机床回转角β1和β2角加速度限制范围,曲面上任意刀触点Pi处的复合角加速度定义为:
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&amp;alpha;</mi> <msub> <mi>P</mi> <mi>i</mi> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msup> <mi>d</mi> <mn>2</mn> </msup> <msub> <mi>T</mi> <msub> <mi>P</mi> <mi>i</mi> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <msup> <mi>dt</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>=</mo> <mi>d</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>T</mi> <msub> <mi>P</mi> <mi>i</mi> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> </mrow> </mfrac> <msub> <mi>&amp;omega;</mi> <mn>1</mn> </msub> <mo>+</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>T</mi> <msub> <mi>P</mi> <mi>i</mi> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> </mrow> </mfrac> <msub> <mi>&amp;omega;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>/</mo> <mi>d</mi> <mi>t</mi> <mo>=</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mn>2</mn> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msub> <mi>T</mi> <msub> <mi>P</mi> <mi>i</mi> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> </mrow> </mfrac> <msub> <mi>&amp;omega;</mi> <mn>1</mn> </msub> <msub> <mi>&amp;omega;</mi> <mn>2</mn> </msub> <mo>+</mo> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msub> <mi>T</mi> <msub> <mi>P</mi> <mi>i</mi> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> </mrow> </mfrac> <msup> <msub> <mi>&amp;omega;</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> <mo>+</mo> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msub> <mi>T</mi> <msub> <mi>P</mi> <mi>i</mi> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> <mn>2</mn> </msup> </mrow> </mfrac> <msup> <msub> <mi>&amp;omega;</mi> <mn>2</mn> </msub> <mn>2</mn> </msup> <mo>+</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>T</mi> <msub> <mi>P</mi> <mi>i</mi> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> </mrow> </mfrac> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> <mo>+</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <msub> <mi>T</mi> <msub> <mi>P</mi> <mi>i</mi> </msub> </msub> <mrow> <mo>(</mo> <msub> <mi>&amp;beta;</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>&amp;beta;</mi> <mn>2</mn> </msub> </mrow> </mfrac> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
式中,为刀触点Pi处的刀轴矢量,t代表时间。
6.根据权利要求1所述的基于AC型五轴数控机床球头刀加工刀轴矢量光顺方法,其特征在于,所述步骤E具体为:
(1)利用现有球头刀五轴刀位优化方法生成给定曲面采样刀触点处的刀位,然后获得上述采样刀触点处的初始刀轴矢量;
(2)利用式(5)计算出各采样刀触点处所对应的机床回转角A和C,再利用三次样条插值函数分别将各采样刀触点处的机床回转角A和C进行三次样条拟合;
(3)利用步骤(2)所得到的三次样条插值拟合函数计算每行刀轨刀触曲线上各刀触点处机床回转角A和C,然后利用式(6)计算各刀触点处球头刀刀位设计变量θ和ψ,再利用式(4)和(3)计算出各刀触点处球头刀刀位点径矢和刀轴矢量,直到求解出该行刀轨刀触曲线上所有刀触点的刀位点径矢和刀轴矢量。
CN201710280231.6A 2017-04-26 2017-04-26 基于ac型五轴数控机床球头刀加工刀轴矢量光顺方法 Active CN107085413B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710280231.6A CN107085413B (zh) 2017-04-26 2017-04-26 基于ac型五轴数控机床球头刀加工刀轴矢量光顺方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710280231.6A CN107085413B (zh) 2017-04-26 2017-04-26 基于ac型五轴数控机床球头刀加工刀轴矢量光顺方法

Publications (2)

Publication Number Publication Date
CN107085413A true CN107085413A (zh) 2017-08-22
CN107085413B CN107085413B (zh) 2019-07-16

Family

ID=59612474

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710280231.6A Active CN107085413B (zh) 2017-04-26 2017-04-26 基于ac型五轴数控机床球头刀加工刀轴矢量光顺方法

Country Status (1)

Country Link
CN (1) CN107085413B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108415374A (zh) * 2018-02-01 2018-08-17 大连理工大学 基于机床旋转进给轴运动学特性的刀轴矢量光顺方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101063880A (zh) * 2007-06-07 2007-10-31 北京航空航天大学 鼓型刀具多坐标宽行加工的刀位优化方法
JP2013117924A (ja) * 2011-12-05 2013-06-13 Fanuc Ltd 工具背面での切削点指令により加工を行う数値制御装置
CN103558808A (zh) * 2013-09-28 2014-02-05 大连理工大学 复杂曲面五轴数控加工刀矢的运动学控制方法
CN103592891A (zh) * 2013-09-28 2014-02-19 大连理工大学 运动学约束的复杂曲面五轴数控加工刀矢光顺方法
CN103869757A (zh) * 2014-03-26 2014-06-18 大连理工大学 复杂曲面五轴数控加工刀矢的动力学控制方法
CN104102171A (zh) * 2014-07-04 2014-10-15 山东理工大学 一种球头刀多轴加工刀轴矢量优化方法
CN105676777A (zh) * 2016-01-15 2016-06-15 广东工业大学 面向三角网格曲面的双向插值数控代码生成方法及其应用

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101063880A (zh) * 2007-06-07 2007-10-31 北京航空航天大学 鼓型刀具多坐标宽行加工的刀位优化方法
JP2013117924A (ja) * 2011-12-05 2013-06-13 Fanuc Ltd 工具背面での切削点指令により加工を行う数値制御装置
CN103558808A (zh) * 2013-09-28 2014-02-05 大连理工大学 复杂曲面五轴数控加工刀矢的运动学控制方法
CN103592891A (zh) * 2013-09-28 2014-02-19 大连理工大学 运动学约束的复杂曲面五轴数控加工刀矢光顺方法
CN103869757A (zh) * 2014-03-26 2014-06-18 大连理工大学 复杂曲面五轴数控加工刀矢的动力学控制方法
CN104102171A (zh) * 2014-07-04 2014-10-15 山东理工大学 一种球头刀多轴加工刀轴矢量优化方法
CN105676777A (zh) * 2016-01-15 2016-06-15 广东工业大学 面向三角网格曲面的双向插值数控代码生成方法及其应用

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
徐汝锋: "《基于机床运动学约束球头刀多轴加工刀轴》", 《机械工程学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108415374A (zh) * 2018-02-01 2018-08-17 大连理工大学 基于机床旋转进给轴运动学特性的刀轴矢量光顺方法
CN108415374B (zh) * 2018-02-01 2019-08-13 大连理工大学 基于机床旋转进给轴运动学特性的刀轴矢量光顺方法

Also Published As

Publication number Publication date
CN107085413B (zh) 2019-07-16

Similar Documents

Publication Publication Date Title
CN105425727B (zh) 五轴侧铣加工刀具路径光顺方法
CN106950916A (zh) 基于ab型五轴数控机床环形刀加工刀轴矢量光顺方法
CN104460516B (zh) 一种基于后置处理五轴刀具半径补偿方法
CN102608952B (zh) 对采用球头刀具的五轴联动机床平滑加工路径的方法
CN104460515B (zh) 一种基于后置处理五轴刀具长度补偿方法
CN107065769A (zh) 基于ab型五轴数控机床球头刀加工刀轴矢量光顺方法
CN106933190A (zh) 基于bc型五轴数控机床环形刀加工刀轴矢量光顺方法
CN104384586A (zh) 四轴数控铣机床加工整体叶轮的方法
Zhu et al. Formulating the swept envelope of rotary cutter undergoing general spatial motion for multi-axis NC machining
CN103218475A (zh) 一种基于在机测评的复杂空间型面误差反馈补偿方法
CN110297462A (zh) 一种考虑机床几何误差影响的磨齿精度预测建模方法
CN106873522A (zh) 一种非轴对称扫掠曲面的数控车削刀具轨迹规划方法
CN105549535A (zh) 五轴摇篮式数控机床未变形切屑的三维几何建模方法
CN106896782A (zh) 基于bc型五轴数控机床球头刀加工刀轴矢量光顺方法
CN106312850A (zh) 螺杆磨削刀具的设计方法
CN107065777A (zh) 基于ba型五轴数控机床环形刀加工刀轴矢量光顺方法
CN105252232B (zh) 一种加工汽轮机动叶片叶身曲面的方法
CN107045328A (zh) 基于ba型五轴数控机床球头刀加工刀轴矢量光顺方法
CN107748538A (zh) 基于bc型五轴数控机床后处理五轴鼓型刀半径补偿方法
CN112883505B (zh) 考虑刀具工件相对振动的超精密端面车削表面建模方法
CN106933189A (zh) 基于ac型五轴数控机床环形刀加工刀轴矢量光顺方法
CN107085413A (zh) 基于ac型五轴数控机床球头刀加工刀轴矢量光顺方法
CN107357256A (zh) 基于ac型五轴数控机床后处理五轴鼓型刀半径补偿方法
CN113065205A (zh) 一种采用平行砂轮磨削圆弧头后刀面的轨迹求解方法
Deng et al. The theoretical and experimental research of machining equal base circle bevel gears with a pot-shaped milling cutter

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