CN112417537B - 一种基于车削加工将几何误差可视化的表面形貌仿真方法 - Google Patents
一种基于车削加工将几何误差可视化的表面形貌仿真方法 Download PDFInfo
- Publication number
- CN112417537B CN112417537B CN202011484847.3A CN202011484847A CN112417537B CN 112417537 B CN112417537 B CN 112417537B CN 202011484847 A CN202011484847 A CN 202011484847A CN 112417537 B CN112417537 B CN 112417537B
- Authority
- CN
- China
- Prior art keywords
- axis
- error
- steps
- carrying
- treatment
- 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 101
- 238000004088 simulation Methods 0.000 title claims abstract description 56
- 238000007514 turning Methods 0.000 title claims abstract description 24
- 238000012545 processing Methods 0.000 claims abstract description 20
- 238000003754 machining Methods 0.000 claims abstract description 15
- 230000000007 visual effect Effects 0.000 claims abstract description 5
- 238000013178 mathematical model Methods 0.000 claims abstract description 3
- 238000005520 cutting process Methods 0.000 claims description 45
- 238000005094 computer simulation Methods 0.000 claims description 8
- 238000005096 rolling process Methods 0.000 claims description 6
- 229910003460 diamond Inorganic materials 0.000 claims description 4
- 239000010432 diamond Substances 0.000 claims description 4
- 238000009826 distribution Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 10
- 238000012800 visualization Methods 0.000 abstract description 4
- 238000013461 design Methods 0.000 abstract description 3
- 238000005457 optimization Methods 0.000 abstract description 3
- 238000012876 topography Methods 0.000 description 34
- 230000000694 effects Effects 0.000 description 15
- 230000008859 change Effects 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000002493 microarray Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000002301 combined effect Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 3
- 230000000737 periodic effect Effects 0.000 description 3
- 238000007516 diamond turning Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000000877 morphologic effect Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003746 surface roughness Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Numerical Control (AREA)
Abstract
本发明属于超精密加工领域,公开了一种基于车削加工将几何误差可视化的表面形貌仿真方法,包括以下步骤:(a)建立特征曲面的数学模型;(b)结合加工条件,分析在无几何误差下特征曲面经仿真实验得到的表面形貌,得到理想情况下特征曲面的三维表面形貌模型;(c)在一个或几个运动轴存在几何误差的情况下,分析在预先设定的几何误差值下特征曲面经仿真实验得到的表面形貌,得到可视化结果,进行比较后即可得到运动轴的几何误差对于工件形貌的影响结果。本发明通过对仿真方法整体处理流程设计等进行改进,能够清晰地表征车削加工过程中几何误差与工件三维表面形貌的关系,便于后续基于几何误差可视化的工艺参数优化。
Description
技术领域
本发明属于超精密加工领域,更具体地,涉及一种基于车削加工将几何误差可视化的表面形貌仿真方法,即一种考虑刀具工件相对振动的超精密端面车削表面形貌仿真建模方法,利用该方法能够仿真得到几何误差可视化的工件三维表面。
背景技术
随着国防、航空航天、医疗、光学元器件等行业的发展,对于一些关键元件有了越来越高的要求,主要体现在元件的加工精度、表面粗糙度等方面。这些元件通常来自于超精密机床的加工。
衡量超精密机床很重要的一个指标是加工精度。加工精度决定了工件的表面质量好坏、尺寸精度高低。加工精度的测量有很多方法,其中工件的三维表面形貌可以直接反映加工精度的好坏。工件的三维表面形貌主要由工件与刀具之间的相对振动、位置误差以及车刀的刀齿误差影响,这些误差的存在,会使刀具在实际加工过程中的刀具轨迹与理论轨迹不同,造成工件的尺寸精度降低的结果。补偿几何误差后,提高的几何精度也可以反映到工件的三维表面形貌中。
目前对于机床的误差辨识方法大多是借助球杆仪等设备。用球杆仪对机床刀具的插补运动轨迹进行测量,再根据机床运动模型,从插补运动轨迹中得到各个机床误差。这种非切削误差辨识方法不能反映真实的加工状态下,机床产生的误差对工件加工的影响结果。对于机床使用者和机床厂商来说,在加工状态下进行误差辨识更有实际意义。
发明内容
针对现有技术的以上缺陷或改进需求,本发明的目的在于提供一种基于车削加工将几何误差可视化的表面形貌仿真方法,其中通过对仿真方法整体处理流程设计等进行改进,紧密围绕车削,并对特征曲面进行进一步优选控制,能够清晰地表征车削加工过程中几何误差与工件三维表面形貌的关系,既可以分析单轴几何误差对工件表面形貌的影响,又可以分析多轴几何误差共同影响下的工件表面形貌,研究几何误差对尺寸精度的影响规律,便于后续基于几何误差可视化的工艺参数优化。
为实现上述目的,按照本发明,提供了一种考虑刀具工件相对振动的超精密端面车削表面形貌仿真建模方法,其特征在于,包括以下步骤:
(a)建立加工面的特征曲面的数学模型;
(b)结合加工条件,分析在无几何误差下特征曲面经仿真实验得到的表面形貌,得到理想情况下特征曲面的三维表面形貌模型;
(c)假设单个或多个运动轴存在几何误差,分析在预先设定的几何误差值下特征曲面经仿真实验得到的表面形貌,得到可视化结果,并与所述步骤(b)得到的理想情况下特征曲面的三维表面形貌进行比较,从而得到单轴或多轴几何误差对于工件形貌的影响结果;其中,存在几何误差的运动轴选自车床的平动轴与转动轴中的一个或多个;其中,
平动轴i轴的几何误差选自以下至少一种:沿着i轴运动,在X轴方向上产生的直线度误差δxi;沿着i轴运动,在Y轴方向上产生的直线度误差δyi;沿着i轴运动,在Z轴方向上产生的直线度误差δzi;沿着i轴运动,在X轴方向上产生的滚动误差θxi;沿着i轴运动,在Y轴方向上产生的颠摆误差θyi;沿着i轴运动,在Z轴方向上产生的偏摆误差θzi;
转动轴j轴几何误差选自以下至少一种:随着j轴转动,在X轴方向上产生的跳动误差δxj;随着j轴转动,在Y轴方向上产生的跳动误差δyj;随着j轴转动,在C轴方向上产生的窜动误差δzj;随着j轴转动,绕X轴的转角误差θxj;随着j轴转动,绕Y轴的转角误差θyj;随着j轴转动,绕自身的转角误差θzj。
作为本发明的进一步优选,当车床为三轴车床时,所述平动轴为车床的X轴、车床的Z轴中的至少一者,所述转动轴为车床的转动轴C轴;其中,
车床X轴的几何误差选自以下至少一种:沿着X轴运动,在X轴方向上产生的定位误差δxx;沿着X轴运动,在Y轴方向上产生的直线度误差δyx;沿着X轴运动,在Z轴方向上产生的直线度误差δzx;沿着X轴运动,在X轴方向上产生的滚动误差θxx;沿着X轴运动,在Y轴方向上产生的颠摆误差θyx;沿着X轴运动,在Z轴方向上产生的偏摆误差θzx;
车床Z轴的几何误差选自以下至少一种:沿着Z轴运动,在X轴方向上产生的直线度误差δxz;沿着Z轴运动,在Y轴方向上产生的直线度误差δyz;沿着Z轴运动,在Z轴方向上产生的定位误差δzz;沿着Z轴运动,在X轴方向上产生的滚动误差θxz;沿着Z轴运动,在Y轴方向上产生的颠摆误差θyz;沿着Z轴运动,在Z轴方向上产生的偏摆误差θzz;
车床转动轴C轴的几何误差选自以下至少一种:随着C轴转动,在X轴方向上产生的跳动误差δxc;随着C轴转动,在Y轴方向上产生的跳动误差δyc;随着C轴转动,在C轴方向上产生的窜动误差δzc;随着C轴转动,绕X轴的转角误差θxc;随着C轴转动,绕Y轴的转角误差θyc;随着C轴转动,绕自身的转角误差θzc。
作为本发明的进一步优选,所述步骤(a)中,所述特征曲面同时包含环状阵列和平面结构;优选的,所述特征曲面是在初始圆形区域的平面结构上设置阵列单元沿圆环均匀分布的多个凹凸交错分布的凹凸阵列。
作为本发明的进一步优选,记特征曲面的初始圆形区域的半径为R,则所述凹凸阵列分布在的圆环区域内。
作为本发明的进一步优选,所述步骤(c)中,当所述多个存在几何误差的运动轴同时包含车床的X轴与Z轴时,预先设定的几何误差值满足:
δxx=-δxz,δyx=-δyz,δzx=-δzz。
作为本发明的进一步优选,所述仿真实验具体是通过先对特征曲面模型进行离散化处理,再求解出仿真条件下特征曲面上每个离散点对应的工件轮廓高度值,最后输出工件的三维表面形貌。
作为本发明的进一步优选,所述仿真实验具体包括以下步骤:
(1)以特征曲面的初始圆形区域的圆心为原点,建立笛卡尔直角坐标系,同时按预先设定的角度分辨率Δθ以及预先设定的长度分辨率Δl分别沿圆周方向和半径方向对仿真区域进行分割,使仿真区域沿圆周方向的分割数目Nc、仿真区域沿径向方向的分割数目Nr两者均为整数,实现离散化处理,得到N=Nc·Nr个离散点;其中,Rw为工件半径;
接着,从仿真区域外侧开始,沿内侧对离散的圆周由1至Nr进行编号;同时,从X轴开始,沿顺时针方向对每个圆周上的离散点由1至Nc进行编号,则,第i个圆周上的第j个点P(i,j)在仿真区域内的极坐标为:
式中,i=1~Nr,j=1~Nc;
(2)计算离散点处的切削区域内所有切削刃参考点的半径rk;记第i个参考点的半径为ri,则在此离散点处位于切削区域内所有的参考点的半径满足:其中,Lp为切削区域的长度;
记切削深度为ap,切削刃的半径为Rt,则Lp满足:
(3)计算切削区域内的刀具切削刃轮廓在离散点处的轮廓高度值:
记第k个参考点的坐标为(rk,zk),其中zk为该点在Z轴方向上的坐标,特征曲面的曲面方程为f,则,切削刃圆弧的圆心坐标为(Rk,Zk),且满足:
设离散点P(i,j)处的半径为r,经过第k个刀具轮廓切削后,轮廓高度为hk,则hk满足:
如此求出离散点P(i,j)在每个刀具轮廓切削时的轮廓高度,计算出其中最小的轮廓高度;
最后,记工件表面在P(i,j)处的轮廓高度理论值为h(i,j),再在高度理论值与切削深度之和h(i,j)+ap与最小轮廓高度中取出最小值,该最小值即为切削区域内的刀具切削刃轮廓在离散点处的轮廓高度值。
作为本发明的进一步优选,所述车床为超精密金刚石车床。
通过本发明所构思的以上技术方案,与现有技术相比,本发明中基于车削加工的将几何误差可视化的表面形貌仿真方法,能够清晰地表征几何误差与工件三维表面形貌的关系,进而可实现基于几何误差可视化的工艺参数优化。本方法填补了对于超精密车削加工的工件表面形貌研究很少,不能直观地体现几何误差影响的缺陷,弥补了现有技术的不足。
本发明仿真建模方法所采用的特征曲面优选同时包含环状阵列和平面结构,记工件半径为R,工件表面上在的区域内为内部的平面结构,该平面结构的形貌特征可以直观反映C轴几何误差对表面形貌的影响。工件表面上位于/>的圆环区域内,沿圆周方向凹凸相间地排列着若干个微型阵列,通过这若干个微型阵列的位置变化能够直观反映出导轨几何误差对表面形貌的影响,对于显示车床的几何误差影响结果有很明显的作用。
基于本发明仿真建模方法,通过无几何误差的工件表面形貌图和考虑了几何误差的工件表面形貌图的对比,可以清晰地看出各项几何误差在工件表面的影响作用,尤其适用于对诸如超精密端面车削表面形貌的仿真建模。
附图说明
图1为本发明中超精密车削的三维表面形貌仿真分析过程示意图;其中,图1中的(a)为用一基面截刀具轨迹,图1中的(b)为将仿真区域离散化,图1中的(c)为基面内截面的切削示意图,图1中的(d)为仿真点轮廓高度和刀具位置间的关系。
图2(a)和图2(b)均为特征曲面的三维形状及其俯视图。其中,图2(a)为特征曲面的理想曲面形状,图2(b)为无几何误差下车床加工得到的特征曲面三维表面形貌图。
图3为本发明中分别考虑X轴上几何误差的表面形貌图;其中,图3中的(a)考虑δxx的影响,图3中的(b)考虑δyx的影响,图3中的(c)考虑δzx的影响,图3中的(d)考虑θxx的影响,图3中的(e)考虑θyx的影响,图3中的(f)考虑θzx的影响。
图4为本发明中X轴上的六项几何误差对三维表面形貌的影响示意图。
图5为本发明中分别考虑Z轴上几何误差的表面形貌图;其中,图5中的(a)考虑δxz的影响,图5中的(b)考虑δyz的影响,图5中的(c)考虑δzz的影响,图5中的(d)考虑θxz的影响,图5中的(e)考虑θyz的影响,图5中的(f)考虑θzz的影响。
图6为Z轴上的六项几何误差对三维表面形貌的影响。
图7为分别考虑C轴上几何误差的表面形貌图;其中,图7中的(a)考虑δxc的影响,图7中的(b)考虑δyc的影响,图7中的(c)考虑δzc的影响,图7中的(d)考虑θxc的影响,图7中的(e)考虑θyc的影响,图7中的(f)考虑θzc的影响。
图8为本发明中C轴上的六项几何误差对三维表面形貌的影响。
图9为与X轴和Z轴有关的12项几何误差对三维表面形貌的共同影响。
图10为与X轴和C轴有关的12项几何误差对三维表面形貌的共同影响。
图11为与Z轴和C轴有关的12项几何误差对三维表面形貌的共同影响。
图12为与X轴、Z轴、C轴有关的18项几何误差对三维表面形貌的共同影响。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
总的来说,本发明中基于车削加工将几何误差可视化的表面形貌仿真方法,包括以下步骤:
步骤a,建立加工面的特征曲面模型;
步骤b,结合加工条件,建立理想情况下特征曲面的三维表面形貌模型(加工条件,可以参照常规的车削工艺,例如可包含车削过程中的主轴参数与X轴的进给速度。由于加工条件的改变会影响工件的表面形貌,因此本发明预先设定了一定的加工条件,排除了加工条件对表面形貌的干扰);
步骤c,对各运动轴分别分析几何误差对工件表面形貌的影响规律。对各个运动轴分别进行仿真实验。在对每个轴的仿真分析中,分别施加对此轴产生影响的几何误差作用,通过分析在每个几何误差下特征曲面的表面形貌,可以得出该几何误差对工件形貌的影响结果;
步骤d,分析多轴几何误差共同影响下的工件表面形貌,研究几何误差对尺寸精度的影响规律。对运动轴两两耦合以及所有运动轴共同耦合进行仿真实验,可以得到几何误差对于工件加工精度的影响规律。
具体的,三轴机床为例,所述步骤a中,可以包括以下内容:
工件半径为R,工件表面上在的区域内为内部的平面结构,该平面结构的形貌特征可以直观反映C轴几何误差对表面形貌的影响。工件表面上位于/>的圆环区域内,沿圆周方向凹凸相间地排列着8个微型阵列,通过这8个微型阵列的位置变化能够直观反映出导轨几何误差对表面形貌的影响。
所述步骤b中,可以包括以下内容:
1.对仿真区域进行离散。
(1)在基于超精密金刚石车削加工的特征曲面三维表面形貌建模方法中,由于车削加工的毛坯件通常为圆柱形棒料,所以取仿真区域为圆形区域,在极坐标系内将仿真区域等间隔离散成一系列离散点。特征曲面的车削过程仿真建模就是求出特征曲面上每个离散点对应的工件轮廓高度值。
(2)首先沿圆周方向对仿真区域进行离散,设角度分辨率为Δθ,则仿真区域沿圆周方向的分割数目Nc为:
(3)然后沿径向方向对仿真区域进行离散,设长度分辨率为Δl,仿真区域半径也就是工件半径为Rw,则仿真区域沿径向方向的分割数目Nr为:
为了方便,一般情况下可预先设定合适的Δθ和Δl使Nc和Nr为整数,则仿真区域内所有离散点的数目N为:
N=Nc·Nr
(4)过特征曲面圆心作笛卡尔直角坐标系,从仿真区域外侧开始,沿内侧对离散的圆周进行编号(1~Nr);从X轴开始,沿顺时针方向对圆周上的离散点进行编号(1~Nc),则第i个圆周上的第j个点P(i,j)在仿真区域内的极坐标为:
式中,i=1~Nr,j=1~Nc。
2.计算离散点处的切削区域内所有切削刃参考点的半径rk。如图1中的(a)所示,用一个与X轴正方向成θ的截面截特征曲面。假设第i个参考点的半径为ri,则在此离散点处位于切削区域内所有的参考点的半径满足:
式中的Lp为切削区域的长度,假设切削深度为ap,切削刃的半径为Rt,则Lp满足:
3.计算切削区域内的刀具切削刃轮廓在离散点处的轮廓高度值。
用一个与X轴正方向成θ的截面截加工面,得到图1中的(d)所示的截面。假设第k个参考点的坐标为(rk,zk),特征曲面的曲面方程为f(r,θ)。金刚石车刀在加工过程中通常沿曲面法向方向进行刀补。曲面在参考点处的法向量为切削刃圆弧的圆心坐标为(Rk,Zk),其满足:
设离散点P(i,j)处的半径为r,经过第k个刀具轮廓切削后,轮廓高度为hk,则hk满足:
使用上式可以求出离散点P(i,j)在每个刀具轮廓切削时的轮廓高度,计算出其中最小的轮廓高度。记工件表面在P(i,j)处的轮廓高度理论值为h(i,j)(该轮廓高度理论值,即,根据特征曲面的曲面方程得到的值),再在高度理论值与切削深度之和h(i,j)+ap与最小轮廓高度中取出最小值,即为切削区域内的刀具切削刃轮廓在离散点处的轮廓高度值。
4.计算完全部离散点后,输出工件的三维表面形貌。如图2(b)所示。
以下为具体实施例:
实施例1
本实施例提供了一种基于车削加工的将几何误差可视化的表面形貌仿真方法,以NanoformX型三轴超精密金刚石车床为例,具体仿真方法实现步骤如下所示:
步骤1,设计特征曲面。
用超精密金刚石车床对工件的端面进行车削,工件半径为0.5mm。特征曲面为:在半径0.25mm到0.5mm的环状区域内,均匀地分布8个环状微型阵列,阵列凹凸相间。在半径小于0.25mm的区域内为平面。通过环状阵列的位置变化,可以分析出X、Y轴几何误差对表面形貌的影响,而平面区域的表面形貌特征变化可以反映出回转轴C轴的几何误差对工件加工的影响。主轴转速为100r/min,进给速度为5mm/min,切削深度为5μm,刀具切削刃半径为0.1mm。则在无几何误差的情况下,特征曲面的曲面方程为(式中c代表回转轴C轴的旋转角度):
特征曲面的三维表面形貌如图2(a)所示。
步骤2,对各运动轴相关的几何误差分别进行分析。
在对超精密车床建模过程中可以得到21项几何误差与刀具和工件的相对位置有关。在对几何误差积累过程的分析中,垂直度误差不是单个机床部件的几何误差。因此不考虑垂直度误差,只分析各运动轴的各项几何误差。
(1)X轴相关的几何误差对工件形貌的影响
与X轴相关的几何误差一共有6项,将设计7组实验分别与理想情况下的工件表面形貌做对比。如下表所示:
考虑X轴相关几何误差的特征曲面的三维表面形貌特征如图3、图4所示。
图3中的(a)-(f)展示了仿真参数为a-f的结果。从表面形貌特征图中可以直观地看出各几何误差对加工精度的影响:δxx影响外部的环状阵列沿径向移动,δyx使环状阵列沿周向移动,δzx改变整个特征曲面的轮廓高度,θyx使特征曲面的表面形貌产生倾斜的变化,而θxx和θzx对表面形貌没有很大的影响。图4是仿真实验g的表面形貌,反映了所有X轴相关的几何误差对工件表面的共同作用。分别用与X轴成45°、270°的平面截取特征曲面曲面,并截取半径为0.375mm的周向曲截面。在两个径向方向的截面结果中,可以看出δxx对径向的影响,δzx对纵向的轮廓高度的影响,以及θyx对轮廓的倾斜程度的影响。从周向的曲截面结果中可以看出δyx对加工轨迹周向移动的影响。
(2)Z轴相关的几何误差对工件形貌的影响
与Z轴相关的几何误差一共有6项,将设计7组实验分别与理想情况下的工件表面形貌做对比。如下表所示:
考虑Z轴相关几何误差的特征曲面的三维表面形貌特征如图5、图6所示。
图5展示了仿真参数为a-f的结果。从表面形貌特征图中可以直观地看出各几何误差对加工精度的影响:δxz影响外部的环状阵列沿径向移动,δyz使环状阵列沿周向移动,δzz改变整个特征曲面的轮廓高度,从图5中可以看出θxz、θyz、θzz对加工精度的影响很小。图6是仿真实验g的表面形貌,反映了所有Z轴相关的几何误差对工件表面的共同作用。分别用与X轴成45°、270°的平面截取特征曲面曲面,并截取半径为0.375mm的周向曲截面。在两个径向方向的截面结果中,可以看出δxz对径向的影响,δzz对纵向的轮廓高度的影响,从周向的曲截面结果中可以看出δyz对加工轨迹周向移动的影响。
(3)C轴相关的几何误差对工件形貌的影响
C轴是一个转动轴,与其相关的几何误差是随着角度的变化而变化的,且产生的影响呈周期性。本实验将与C轴有关的几何误差设计为三角函数形式,随着旋转角度θ的变化,值周期性变化。这样的几何误差一共有6项,将设计7组实验分别与理想情况下的工件表面形貌做对比。如下表所示(表中的C表示旋转角度θ,单位是rad):
考虑C轴相关几何误差的特征曲面的三维表面形貌特征如图7、图8所示。
图7展示了仿真参数为a-f的结果。从表面形貌特征图中可以直观地看出各几何误差对加工精度的影响:从图7中的(a)、(b)、(f)中可以看出:δxc、δyc、θzc影响环状阵列沿径向偏移,偏移程度不同。从图7中的(c)中可以看出存在δzc时,特征曲面的表面形貌产生四条脊线,脊线轮廓高度相同。图7中的(d)、(e)显示θyc与θxc的影响是:环状阵列区域产生轮廓高度的变化,平面区域出现五条脊线,表面形貌产生倾斜趋势。图8是仿真实验g的表面形貌,反映了所有C轴相关的几何误差对工件表面的共同作用。分别用与X轴成45°、270°的平面截取特征曲面曲面,并截取半径为0.375mm的周向曲截面。在两个径向方向的截面结果中,径向方向上可以看出δxc、δyc、θzc对环状阵列的影响,轮廓高度上可以看出δzc的影响,倾斜趋势上可以看出θxc和θyc的影响。周向表面形貌受六个几何误差的影响较小。
(4)多轴几何误差对工件形貌的影响
将设置四组实验讨论多轴几何误差对工件的影响:分别是X轴与Z轴,X轴与C轴,Z轴与C轴,X轴、Z轴、C轴的耦合。从前面的分析可以看出X轴和Z轴的线性几何误差产生的影响作用方向相反,因此设置X轴的线性误差为负值,大小与Z轴的线性几何误差相同。实验的具体设置如下表所示(表中的C表示旋转角度θ,单位是rad):
从图9中可以看出X轴与Z轴的几何误差耦合作用是将特征曲面的表面形貌沿径向、周向方向偏移,同时影响整个特征曲面的轮廓高度并产生了倾斜的趋势。从图10和11中可以看出X轴和Z轴的几何误差主要影响特征曲面的径向和周向,C轴会使特征曲面的表面产生周期性的纹理。综上所述,X轴与Z轴的线性几何误差类似,都是使工件形貌沿径向、周向偏移,轮廓高度也会改变,但X轴与Z轴产生的影响方向相反,关系是线性叠加的。C轴主要对特征曲面产生周期性的表面特征,只与C轴的旋转角度有关,和X轴、Z轴的运动无关。
可见,本实施例实现了基于车削加工的将几何误差可视化的表面形貌仿真。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种考虑刀具工件相对振动的超精密端面车削表面形貌仿真建模方法,其特征在于,包括以下步骤:
(a)建立加工面的特征曲面的数学模型;
(b)结合加工条件,分析在无几何误差下特征曲面经仿真实验得到的表面形貌,得到理想情况下特征曲面的三维表面形貌模型;
(c)假设单个或多个运动轴存在几何误差,分析在预先设定的几何误差值下特征曲面经仿真实验得到的表面形貌,得到可视化结果,并与所述步骤(b)得到的理想情况下特征曲面的三维表面形貌进行比较,从而得到单轴或多轴几何误差对于工件形貌的影响结果;其中,存在几何误差的运动轴选自车床的平动轴与转动轴中的一个或多个;其中,
平动轴i轴的几何误差选自以下至少一种:沿着i轴运动,在X轴方向上产生的直线度误差δxi;沿着i轴运动,在Y轴方向上产生的直线度误差δyi;沿着i轴运动,在Z轴方向上产生的直线度误差δzi;沿着i轴运动,在X轴方向上产生的滚动误差θxi;沿着i轴运动,在Y轴方向上产生的颠摆误差θyi;沿着i轴运动,在Z轴方向上产生的偏摆误差θzi;
转动轴j轴几何误差选自以下至少一种:随着j轴转动,在X轴方向上产生的跳动误差δxj;随着j轴转动,在Y轴方向上产生的跳动误差δyj;随着j轴转动,在C轴方向上产生的窜动误差δzj;随着j轴转动,绕X轴的转角误差θxj;随着j轴转动,绕Y轴的转角误差θyj;随着j轴转动,绕自身的转角误差θzj。
2.如权利要求1所述方法,其特征在于,当车床为三轴车床时,所述平动轴为车床的X轴、车床的Z轴中的至少一者,所述转动轴为车床的转动轴C轴;其中,
车床X轴的几何误差选自以下至少一种:沿着X轴运动,在X轴方向上产生的定位误差δxx;沿着X轴运动,在Y轴方向上产生的直线度误差δyx;沿着X轴运动,在Z轴方向上产生的直线度误差δzx;沿着X轴运动,在X轴方向上产生的滚动误差θxx;沿着X轴运动,在Y轴方向上产生的颠摆误差θyx;沿着X轴运动,在Z轴方向上产生的偏摆误差θzx;
车床Z轴的几何误差选自以下至少一种:沿着Z轴运动,在X轴方向上产生的直线度误差δxz;沿着Z轴运动,在Y轴方向上产生的直线度误差δyz;沿着Z轴运动,在Z轴方向上产生的定位误差δzz;沿着Z轴运动,在X轴方向上产生的滚动误差θxz;沿着Z轴运动,在Y轴方向上产生的颠摆误差θyz;沿着Z轴运动,在Z轴方向上产生的偏摆误差θzz;
车床转动轴C轴的几何误差选自以下至少一种:随着C轴转动,在X轴方向上产生的跳动误差δxc;随着C轴转动,在Y轴方向上产生的跳动误差δyc;随着C轴转动,在C轴方向上产生的窜动误差δzc;随着C轴转动,绕X轴的转角误差θxc;随着C轴转动,绕Y轴的转角误差θyc;随着C轴转动,绕自身的转角误差θzc。
3.如权利要求1所述方法,其特征在于,所述步骤(a)中,所述特征曲面同时包含环状阵列和平面结构;所述特征曲面是在初始圆形区域的平面结构上设置阵列单元沿圆环均匀分布的多个凹凸交错分布的凹凸阵列。
4.如权利要求3所述方法,其特征在于,记特征曲面的初始圆形区域的半径为R,则所述凹凸阵列分布在的圆环区域内。
5.如权利要求1所述方法,其特征在于,所述步骤(c)中,当所述多个存在几何误差的运动轴同时包含车床的X轴与Z轴时,预先设定的几何误差值满足:
δxx=-δxz,δyx=-δyz,δzx=-δzz。
6.如权利要求1所述方法,其特征在于,所述仿真实验具体是通过先对特征曲面模型进行离散化处理,再求解出仿真条件下特征曲面上每个离散点对应的工件轮廓高度值,最后输出工件的三维表面形貌。
7.如权利要求6所述方法,其特征在于,所述仿真实验具体包括以下步骤:
(1)以特征曲面的初始圆形区域的圆心为原点,建立笛卡尔直角坐标系,同时按预先设定的角度分辨率Δθ以及预先设定的长度分辨率Δl分别沿圆周方向和半径方向对仿真区域进行分割,使仿真区域沿圆周方向的分割数目Nc、仿真区域沿径向方向的分割数目Nr两者均为整数,实现离散化处理,得到N=Nc·Nr个离散点;其中,Rw为工件半径;
接着,从仿真区域外侧开始,沿内侧对离散的圆周由1至Nr进行编号;同时,从X轴开始,沿顺时针方向对每个圆周上的离散点由1至Nc进行编号,则,第i个圆周上的第j个点P(i,j)在仿真区域内的极坐标为:
式中,i=1~Nr,j=1~Nc;
(2)计算离散点处的切削区域内所有切削刃参考点的半径rk;记第i个参考点的半径为ri,则在此离散点处位于切削区域内所有的参考点的半径满足:其中,Lp为切削区域的长度;
记切削深度为ap,切削刃的半径为Rt,则Lp满足:
(3)计算切削区域内的刀具切削刃轮廓在离散点处的轮廓高度值:
记第k个参考点的坐标为(rk,zk),其中zk为该点在Z轴方向上的坐标,特征曲面的曲面方程为f,则,切削刃圆弧的圆心坐标为(Rk,Zk),且满足:
设离散点P(i,j)处的半径为r,经过第k个刀具轮廓切削后,轮廓高度为hk,则hk满足:
如此求出离散点P(i,j)在每个刀具轮廓切削时的轮廓高度,计算出其中最小的轮廓高度;
最后,记工件表面在P(i,j)处的轮廓高度理论值为h(i,j),再在高度理论值与切削深度之和h(i,j)+ap与最小轮廓高度中取出最小值,该最小值即为切削区域内的刀具切削刃轮廓在离散点处的轮廓高度值。
8.如权利要求1所述方法,其特征在于,所述车床为超精密金刚石车床。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011484847.3A CN112417537B (zh) | 2020-12-16 | 2020-12-16 | 一种基于车削加工将几何误差可视化的表面形貌仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011484847.3A CN112417537B (zh) | 2020-12-16 | 2020-12-16 | 一种基于车削加工将几何误差可视化的表面形貌仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112417537A CN112417537A (zh) | 2021-02-26 |
CN112417537B true CN112417537B (zh) | 2024-02-02 |
Family
ID=74775826
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011484847.3A Active CN112417537B (zh) | 2020-12-16 | 2020-12-16 | 一种基于车削加工将几何误差可视化的表面形貌仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112417537B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114472989A (zh) * | 2022-02-10 | 2022-05-13 | 中国工程物理研究院机械制造工艺研究所 | 表面微结构加工的几何误差补偿方法及加工方法、装置 |
CN114996975B (zh) * | 2022-07-19 | 2022-12-30 | 中国科学院光电技术研究所 | 一种基于超精密车削的光学曲面表面形貌仿真方法 |
CN116728158B (zh) * | 2023-08-09 | 2023-12-08 | 成都飞机工业(集团)有限责任公司 | 一种用于五轴机床R-test检测的误差检测结果可视化方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102490081A (zh) * | 2011-11-14 | 2012-06-13 | 华中科技大学 | 基于球头铣削加工的工件三维表面形貌的仿真方法 |
CN102592035A (zh) * | 2012-03-20 | 2012-07-18 | 北京航空航天大学 | 一种车铣复合切削加工表面粗糙度及表面形貌仿真预测方法 |
CN106228606A (zh) * | 2016-07-29 | 2016-12-14 | 南京航空航天大学 | 一种三维椭圆振动辅助切削微织构形貌建模方法 |
CN107526875A (zh) * | 2017-07-31 | 2017-12-29 | 电子科技大学 | 一种航空叶片型面加工误差可视化方法 |
CN107942931A (zh) * | 2018-01-09 | 2018-04-20 | 吉林大学 | 一种正弦柱面超精密车削螺旋刀具轨迹生成方法 |
CN108508848A (zh) * | 2018-04-20 | 2018-09-07 | 华中科技大学 | 一种基于插补数据的铣削加工轮廓误差的评估方法 |
-
2020
- 2020-12-16 CN CN202011484847.3A patent/CN112417537B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102490081A (zh) * | 2011-11-14 | 2012-06-13 | 华中科技大学 | 基于球头铣削加工的工件三维表面形貌的仿真方法 |
CN102592035A (zh) * | 2012-03-20 | 2012-07-18 | 北京航空航天大学 | 一种车铣复合切削加工表面粗糙度及表面形貌仿真预测方法 |
CN106228606A (zh) * | 2016-07-29 | 2016-12-14 | 南京航空航天大学 | 一种三维椭圆振动辅助切削微织构形貌建模方法 |
CN107526875A (zh) * | 2017-07-31 | 2017-12-29 | 电子科技大学 | 一种航空叶片型面加工误差可视化方法 |
CN107942931A (zh) * | 2018-01-09 | 2018-04-20 | 吉林大学 | 一种正弦柱面超精密车削螺旋刀具轨迹生成方法 |
CN108508848A (zh) * | 2018-04-20 | 2018-09-07 | 华中科技大学 | 一种基于插补数据的铣削加工轮廓误差的评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112417537A (zh) | 2021-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112417537B (zh) | 一种基于车削加工将几何误差可视化的表面形貌仿真方法 | |
JP5448634B2 (ja) | 機械の誤差同定方法およびプログラム | |
EP3367056B1 (en) | Error identification method for machine tool | |
US7251580B2 (en) | Method for measuring curved surface of workpiece, program and medium thereof | |
Kim et al. | Texture prediction of milled surfaces using texture superposition method | |
CN112580160B (zh) | 一种用于成形磨齿机在机测量***标定方法 | |
CN108568567B (zh) | 一种基于通用四轴数控机床与球头铣刀的螺旋锥齿轮加工方法 | |
Li et al. | Sensitivity analysis of relationship between error motions and machined shape errors in five-axis machining center-Peripheral milling using square-end mill as test case | |
KR20100119494A (ko) | 가공 장치 및 가공 방법 | |
Chang et al. | Five-axis automated measurement by coordinate measuring machine | |
Wu et al. | A novel prediction method of machining accuracy for machine tools based on tolerance | |
CN111283478A (zh) | 一种类回转零件加工自适应定位方法 | |
Zha et al. | An accuracy evolution method applied to five-axis machining of curved surfaces | |
KR101200171B1 (ko) | 5축 cνc 가공기계를 이용한 드럼 캠의 황삭가공을 위한 후처리방법 | |
CN112799299A (zh) | 一种机器人多轴铣削稳定性模型及其构建方法 | |
Yan et al. | A dimension-driven adaptive programming for tool-path planning and post-processing in 5-axis form milling of hyperboloidal-type normal circular-arc gears | |
Chen et al. | A novel geometric error compensation method for improving machining accuracy of spiral bevel gear based on inverse kinematic model | |
KR101199086B1 (ko) | 5축 cνc 가공기계를 이용한 롤러 기어 캠의 가공을 위한 후처리방법 | |
Guo et al. | Study on tooth profile error of cylindrical gears manufactured by flexible free-form milling | |
Xiang et al. | Five-Axis Numerical Control Machining of the Tooth Flank of a Logarithmic Spiral Bevel Gear Pinion | |
CN113427320A (zh) | 一种降低测量不确定度的多轴在机测量规划方法 | |
JP2007203457A (ja) | Nc加工方法およびnc加工機 | |
Zhang et al. | A CNC program module based on polar coordinate system | |
JP2005098752A (ja) | ブローチの形状測定装置 | |
Thasana et al. | A comparison between simulation and experiment of virtual machining in CNC turning machine considering kinematic motion deviations, tool wear and workpiece deflection errors |
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 |