CN108416153A - 自由态与夹紧态下薄壁曲面件加工轨迹映射方法 - Google Patents
自由态与夹紧态下薄壁曲面件加工轨迹映射方法 Download PDFInfo
- Publication number
- CN108416153A CN108416153A CN201810222492.7A CN201810222492A CN108416153A CN 108416153 A CN108416153 A CN 108416153A CN 201810222492 A CN201810222492 A CN 201810222492A CN 108416153 A CN108416153 A CN 108416153A
- Authority
- CN
- China
- Prior art keywords
- cutter
- state
- contact point
- face piece
- curved face
- 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
Links
Classifications
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
Abstract
本发明自由态与夹紧态下薄壁曲面件加工轨迹映射方法属于薄壁曲面零件精密加工技术领域,涉及一种自由态与夹紧态下薄壁曲面件加工轨迹映射方法。该方法通过建立自由态与夹紧态下薄壁曲面件加工轨迹映射关系数学模型,根据薄壁曲面件原始设计几何模型获得自由态加工信息,利用中性层壳体有限元分析获得薄壁曲面件夹紧态加工信息。通过刀触点网格映射建立薄壁曲面件自由态与夹紧态下刀触点对应关系,最终生成夹紧态下薄壁曲面件刀具加工轨迹,由此确保装夹力释放后薄壁曲面件满足加工精度要求。该方法通过建立自由态与夹紧态下薄壁曲面件加工轨迹映射关系,再规划夹紧态下薄壁曲面加工轨迹,提高了加工质量,适用于薄壁曲面件精密加工。
Description
技术领域
本发明属于薄壁曲面零件精密加工技术领域,涉及一种自由态与夹紧态下薄壁曲面件加工轨迹映射方法。
背景技术
随着航空航天国防装备以及医疗民生等领域的飞速发展,薄壁曲面零件以其质量轻、比强度高等特点受到广泛关注。但是由于薄壁曲面件刚性低,在加工前装夹力作用下易产生弹性变形,从而破坏刀具和工件间理想的相对位置。若按照原始设计出的薄壁曲面件规划刀具轨迹进行加工,待加工结束撤去夹具,工件将回弹变形造成加工误差,最终致使加工精度无法满足要求。研究表明,20%至60%的加工误差是由装夹引起的,但是目前对工件-夹具系统装夹变形造成的加工误差控制上尚缺乏高效系统的方法,对装夹误差的研究主要集中在两个方面:一是最小装夹力的计算,二是最佳装夹顺序的规划。然而,薄壁曲面件加工完成装夹力释放后由于回弹变形及由其导致的加工误差的研究较少。鉴于装夹变形诱发的夹紧态下薄壁曲面件与原始设计存在差异,易致使按原始设计出的薄壁曲面件规划刀具加工轨迹进行加工难以满足加工精度要求,由此迫切需要研究一种自由态与夹紧态下薄壁曲面件加工轨迹映射方法,为实现基于原始设计的夹紧态复杂薄壁曲面件加工轨迹再规划奠定基础,以保证装夹力释放产生回弹变形后复杂薄壁曲面件满足精度要求。
陈蔚芳等人专利公开号CN104111625B的“一种薄壁异形件装夹变形主动加工方法”,该专利通过对薄壁异形件的夹紧力进行优化设置,然后找出工件受力点范围,接着对工件进行辅助支撑的方式以达到变形最小的目的,并最终通过增减刀补的方式补偿变形。然而,该方法主要集中在优化夹紧力以减小工件变形上,并且增减刀补的补偿方法对具有变形不均匀、变形方向不确定的薄壁曲面件加工并不适用。王军等人的文献“薄壁壳体件装夹变形机理有限元分析与控制”,兵工学报,2011,32(8),1008-1013,该文献采用有限元分析方法,对装夹过程中的薄壁壳体件装夹方案进行了优选,模拟了在集中载荷与均布载荷作用下,装夹位置、装夹顺序及加载方式对薄壁壳体零件变形的影响。然而,该方法仅从最优装夹布局方面减小工件变形,并未考虑装夹力释放后产生的回弹变形,具有一定的局限性。
发明内容
本发明针对现有技术缺陷,发明了一种自由态与夹紧态下薄壁曲面件加工轨迹映射方法。该方法建立了自由态与夹紧态下薄壁曲面件加工轨迹映射关系数学模型,根据薄壁曲面件原始设计几何模型经由CAM后处理获得自由态加工信息,利用中性层壳体有限元分析获得薄壁曲面件夹紧态加工信息,通过刀触点网格映射建立薄壁曲面件自由态与夹紧态下刀触点对应关系,最终生成夹紧态下薄壁曲面件刀具加工轨迹,该方法考虑装夹力释放后产生的回弹变形,再规划夹紧态下薄壁曲面加工轨迹,由此确保装夹力释放后薄壁曲面件满足加工精度要求,提高了加工质量。
本发明采用的技术方案是一种自由态与夹紧态下薄壁曲面件加工轨迹映射方法,其特征在于,该方法首先建立了自由态与夹紧态下薄壁曲面件加工轨迹映射数学模型;利用有限元数值分析方法,结合薄壁曲面件几何特征和加工过程中装夹布局,求解自由态与夹紧态下的有限元网格、单元、节点信息,通过加工刀触点网格映射方法获得夹紧态下薄壁曲面件刀触点;在此基础上进行最小二乘曲面局部拟合求解夹紧态下薄壁曲面件刀触点法矢,结合自由态下薄壁曲面件加工信息,进行夹紧态刀位点再生成及规划处理,最终减少薄壁曲面件加工完成后装夹回弹变形引起的加工误差;方法的具体步骤如下:
步骤1,自由态与夹紧态下薄壁曲面件加工轨迹映射模型建立
首先,利用UG软件对薄壁圆筒件进行自由态几何建模,在薄壁圆筒件表面加工图案参数方程为:首先,利用UG软件对薄壁圆筒件进行自由态几何建模,在薄壁圆筒件表面加工图案参数方程为:
式中,r表示极坐标系下极径大小,r=10+(3×sin(2.5θ))2,θ表示极坐标系下极角大小;xt、yt、zt分别表示直角坐标系下x、y、z坐标值。
然后利用CAM软件,根据薄壁曲面件几何特征及加工要求,对其加工轨迹规划,后处理获得刀位点、刀触点坐标信息。以薄壁曲面件刀触面为研究对象,定义自由态与夹紧态下薄壁曲面件加工刀触点间映射关系为f。
自由态时刀触点为Di=[xi,yi,zi,1],则自由态刀位点为:
Wi=Di+R·en (2)
其中,R为刀具半径,en为自由态刀触点单位法矢。
夹紧态时刀触点为:
则夹紧态刀位点为:
其中,为夹紧态刀触点单位法矢。
根据式(1)-(4),建立自由态与夹紧态下薄壁曲面件刀位点映射关系如下:
由此,夹紧态下薄壁曲面件加工轨迹生成的关键在于,其一求解自由态与夹紧态下薄壁曲面件加工刀触点间映射关系f,其二将夹紧态刀触点沿法矢生成夹紧态下薄壁曲面件刀位点。
步骤2,自由态与夹紧态下薄壁曲面件加工刀触点间映射关系求解
采用有限元方法计算夹紧态薄壁曲面件加工刀触点,夹紧态薄壁曲面件加工刀触点坐标的求解过程实质上是自由态下刀触点网格映射的过程,在计算出自由态下刀触点网格节点及其所在的薄壁曲面件有限元网格单元后,确定出自由态刀触点相对于薄壁曲面件有限元网格单元的位置关系,然后根据位置关系计算出自由态刀触点随薄壁曲面件有限元网格单元的变形而形成的夹紧态刀触点坐标。
有限元分析中,变形前后薄壁曲面件有限元网格单元和节点之间存在一一对应关系;有如下映射法则:映射前后刀触点所在的有限元网格单元相同,映射前后刀触点在各自有限元网格单元中的位置相似;
自由态薄壁曲面件加工刀触点所在的有限元网格单元位置分三种情况:刀触点在有限元网格单元内;刀触点在有限元网格单元边线上;刀触点为有限元网格节点。
根据自由态薄壁曲面件加工刀触点所在位置,分别对每种情况下的夹紧态加工刀触点坐标计算如下:
1)刀触点在有限元网格单元内
根据有限元网格单元类型确定刀触点在有限元网格单元中的位置,采用三角形单元或四边形单元;对于三角形单元,采用面积坐标表示有限元网格单元内点的位置;
自由态刀触点e在有向三角形单元Δijk中,它的面积坐标表示为:
其中,i,j,k为自由态有限元三角形网格节点,按逆时针排列约定;顶点按逆时针方向排列的三角形面积为正,顶点按顺时针方向排列的三角形面积为负,S为三角形单元Δijk的面积,Sejk,Seki,Seij分别为有向三角形单元Δejk,Δeki,Δeij的面积。
有限元三角形单元Δijk所对应的有限元映射单元为有向三角形Δi*j*k*,由于Li,Lj,Lk不变,则夹紧态刀触点e*的坐标为:
其中,i*,j*,k*为夹紧态有限元三角形网格节点,S*为有向三角形单元Δi*j*k*的面积,分别为有向三角形单元Δe*j*k*,Δe*k*i*,Δe*i*j*的面积。
定义有向三角形第一个顶点A1坐标为(x1,y1,z1),第二个顶点A2坐标为(x2,y2,z2),第三个顶点A3坐标为(x3,y3,z3),则该三角形面积S有向可用如下行列式表示:
其内任意一点A(x,y,z)与另外两个顶点A2,A3构成的三角形面积为:
由此,自由态刀触点e的坐标为(xe,ye,ze,1),与之相关联的节点有限元三角形网格节点i,j,k的坐标为(xi,yi,zi,1),(xj,yj,zj,1),(xk,yk,zk,1),则自由态下刀触点面积坐标为:
由于面积坐标与直角坐标换算关系为:
由此可得:
式中,夹紧态刀触点e*的坐标为与之相关联的节点有限元三角形网格节点i*,j*,k*的坐标为
薄壁曲面件夹紧态下,与自由态刀触点相关联的节点有限元三角形网格节点i,j,k的变形矩阵Ti,Tj,Tk为:
式中,(Δxi,Δyi,Δzi),(Δxj,Δyj,Δzj),(Δxk,Δyk,Δzk)分别为节点i,j,k的变形量。
则有:
e*=[Li Lj Lk]·[i j k]T·[Ti Tj Tk] (13)
其中,[Li Lj Lk]=h(e,i,j,k),h表示矩阵中数值Li,Lj,Lk与自由态下刀触点e、有限元三角形网格节点i,j,k间的函数关系。
由此,求得夹紧态与自由态下薄壁曲面件加工刀触点间映射关系f为:
f=h(e,i,j,k)·[i j k]T·[Ti Tj Tk]·e-1 (14)
2)刀触点在有限元网格单元边线上
自由态刀触点e落在有限元网格单元边线节点i,j之间,夹紧态有限元网格单元边线节点i*,j*为对应的变形节点,e和i,j之间存在一个比例系数λ:
其中,(xe,ye,ze,1)为自由态刀触点e的坐标,(xi,yi,zi,1),(xj,yj,zj,1)分别为有限元网格单元节点i,j的坐标。
由于比例系数λ不变,可得夹紧态刀触点e*的坐标为:
即:
其中,为夹紧态刀触点e*的坐标,分别为有限元网格单元节点i*,j*的坐标。
薄壁曲面件夹紧态下,与自由态刀触点相关联的节点有限元三角形网格节点i,j的变形矩阵分别为Ti,Tj,E为单位矩阵:
式中,(Δxi,Δyi,Δzi),(Δxj,Δyj,Δzj)分别为节点i,j的变形量。
式(17)可化为:
化简得:
由此,求得夹紧态与自由态下薄壁曲面件加工刀触点间映射关系f为:
3)刀触点在网格单元节点上
自由态刀触点e落在有限元网格节点i上,夹紧态有限元网格节点i*为对应的变形节点,按照变形前后节点间的一一对应关系,可得夹紧态下刀触点e*的坐标:
e*=i*=i·T=e·Ti (21)
式中,为节点i的变形量。
由此,夹紧态与自由态下薄壁曲面件加工刀触点间映射关系f为:
f=T (22)
步骤3,夹紧态下薄壁曲面件刀触点法矢计算
刀触点网格节点中某一节点P的几何欧拉距离最短的m个数据点为P的m近邻或P的m邻域,m邻域中m值不是随便选取的,若取较大值可以提高邻域点的搜索精度,但是大大增加了搜索时间,若取较小的值,虽然提高了效率,但邻域点的数量太少会令后续拟合误差增大,无法准确反映出真实的局部几何形状特征,兼顾精度和效率取m=22-32,进行最小二乘曲面局部拟合求解夹紧态下薄壁曲面件刀触点法矢。
利用最小二乘法原理对夹紧态下薄壁曲面件加工刀触点附近局部区域进行曲面拟合,对于任意刀触点,它的法向量可以用该点及其邻域点拟合出的最小二乘曲面在该点的法矢替代。拟合的曲面方程为:
其中,a0,a1,a2,a3,a4,a5为拟合系数。利用夹紧态刀触点拟合计算出式(23)的曲面方程,则使函数值W最小,W为:
为此,应满足:
即满足线性方程组:
解方程组获得aq的值,将拟合曲面方程写成参数方程形式r(x,y,z):
求曲面r(x,y)的偏微分记作曲面各个刀触点处单位法矢为:
步骤4,夹紧态下薄壁曲面件加工轨迹再规划
由自由态下薄壁曲面件刀触点坐标,联合自由态与夹紧态下薄壁曲面件加工刀触点间映射关系f以及求解得到的夹紧态下薄壁曲面件刀触点法矢e* n,获得夹紧态下薄壁曲面件加工刀位点。结合自由态下刀轨连接信息,通过后处理,输出机床能够识别的加工文件,实现减少薄壁曲面件加工完成后装夹回弹变形引起的加工误差。
本发明的显著效果和益处是该方法针对薄壁曲面件实际加工过程中,由装夹力释放后产生回弹变形的问题,通过建立自由态与夹紧态下薄壁曲面件加工轨迹映射关系,再规划夹紧态下薄壁曲面加工轨迹,提高了加工质量,该方法适用于薄壁曲面件精密加工。
附图说明
图1—方法整体流程图。
图2—薄壁曲面件装夹测力装置。图中:1-压块,2-传感器,3-支架1,4-支架2,5-工件,6-底座1,7-底座2。
图3—自由态与夹紧态下薄壁曲面件加工刀触点间映射关系流程图。
图4—自由态与夹紧态下薄壁曲面件加工刀触点局部网格映射结果图;图中:X、Y表示坐标位置,□表示自由态有限元网格节点,+表示夹紧态有限元网格节点,○表示自由态刀触点,*表示夹紧态刀触点,1-有限元网格单元。
图5—未经加工轨迹映射的薄壁圆筒件表面图案点加工深度分布图;其中,Z值为薄壁圆筒件表面图案点加工深度,X值为测量点位置。
图6—经加工轨迹映射的薄壁圆筒件表面图案点加工深度分布图;其中,Z值为薄壁圆筒件表面图案点加工深度,X值为测量点位置。
具体实施方式
结合技术方案与附图详细说明本发明的具体实施方式。
薄壁曲面件加工时,由于刚性低,导致装夹力释放后产生回弹变形,影响薄壁曲面件加工质量。为了减小薄壁曲面件加工后装夹力释放后产生回弹变形导致的加工误差,发明了一种自由态与夹紧态下薄壁曲面件加工轨迹映射方法,整体流程如附图1所示。
以规格尺寸为98×1×200(外径×壁厚×长度,单位mm)的7075铝合金薄壁圆筒件加工为例,借助UG软件、ANSYS软件以及MATLAB软件,详细说明本发明实施过程。加工中,主轴转速5000r/min,切深0.3mm,进给速度200mm/min,残余高度0.005mm。
首先,利用UG软件对薄壁圆筒件进行自由态几何建模,在薄壁圆筒件表面加工图案参数方程为:xt=rcosθ,yt=rsinθ,zt=0。其中,r表示极坐标系下极径大小,其值为r=10+(3×sin(2.5θ))2;θ表示极坐标系下极角大小;xt、yt、zt分别表示直角坐标系下x、y、z坐标值。
附图2为薄壁圆筒件装夹测力装置,夹紧力大小分别为164.1N和130.3N,根据实际装夹布局情况,利用ANSYS软件对薄壁圆筒件中性层壳体分析。薄壁
圆筒件材料弹性模量E=71.7GPa,泊松比μ=0.33,定义网格单元类型为壳单元shell181,单元网格尺寸为4mm,后处理得到薄壁圆筒件有限元网格、单元、节点三者间的拓扑关系,以及薄壁圆筒件装夹变形前后有限元网格节点的坐标值和有限元网格节点的变形位移量。
然后,利用MATLAB软件根据附图3自由态与夹紧态下薄壁曲面件加工刀触点间映射关系流程,结合上述过程得到的自由态刀触点信息,以及薄壁圆筒件装夹变形前后有限元网格节点的坐标值、有限元网格节点的变形位移量以及有限元网格、单元、节点三者间的拓扑关系,根据式(14)、(20)、(22)计算生成如附图4所示的自由态与夹紧态下薄壁曲面件加工刀触点网格映射结果,即夹紧态刀触点。在此基础上,利用最小二乘法原理对夹紧态刀触点附近局部区域进行曲面拟合,自由态与夹紧态下薄壁曲面件加工刀触点网格映射结果。在此基础上,利用最小二乘法原理对夹紧态刀触点附近局部区域进行曲面拟合,兼顾精度和效率取m=22,根据式(23)-(28)计算各刀触点处法矢;再根据式(5),获得夹紧态下薄壁曲面件加工刀位点。最后,结合自由态刀具加工轨迹信息,通过后处理,输出机床能够识别的加工文件,实现薄壁曲面件加工完成后装夹回弹变形引起的加工误差避免目的。
为验证方法有效性,进行薄壁圆筒件经加工轨迹映射后与未经加工轨迹映射对比实验,并以三坐标测量机测量薄壁曲面件表面图案加工深度Z值变化作为加工质量好坏的判断标准。实验结果表明,经加工轨迹映射后,薄壁圆筒件表面图案加工深度变化均匀如附图5所示,未经加工轨迹映射的薄壁圆筒件表面图案加工深度变化明显,如附图6所示。未经加工轨迹映射的图案加工深度变化量为296.4μm,测量数据标准差为0.0189,深度变化较大,加工质量较差;经加工轨迹映射后的图案加工深度变化量为69.7μm,相比提高了76.5%,测量数据标准差为0.0786,深度变化较小,加工质量明显提高。由此判定结果与实验结果较好吻合,说明利用本发明的自由态与夹紧态下薄壁曲面件加工轨迹映射方法,可达到薄壁曲面件加工完成后装夹回弹变形引起的加工误差避免目的,提高加工质量,对工程实际中产生装夹回弹变形的薄壁曲面件加工具有重要的指导作用。
Claims (1)
1.一种自由态与夹紧态下薄壁曲面件加工轨迹映射方法,其特征在于,该方法首先建立自由态与夹紧态下薄壁曲面件加工轨迹映射数学模型;利用有限元数值分析方法,结合薄壁曲面件几何特征和加工过程中装夹布局,求解自由态与夹紧态下的有限元网格、单元、节点信息,通过加工刀触点网格映射方法获得夹紧态下薄壁曲面件刀触点;在此基础上进行最小二乘曲面局部拟合求解夹紧态下薄壁曲面件刀触点法矢,结合自由态下薄壁曲面件加工信息,进行夹紧态刀位点再生成及规划处理,最终减小薄壁曲面件加工完成后装夹回弹变形引起的加工误差;方法的具体步骤如下:
步骤1,自由态与夹紧态下薄壁曲面件加工轨迹映射模型建立
首先,利用UG软件对薄壁圆筒件进行自由态几何建模,在薄壁圆筒件表面加工图案参数方程为:
式中,r为极坐标系下极径大小,r=10+(3×sin(2.5θ))2,θ为极坐标系下极角大小;xt、yt、zt分别表示直角坐标系下x、y、z坐标值;
然后利用CAM软件,根据薄壁曲面件几何特征及加工要求,对其加工轨迹规划,后处理获得刀位点、刀触点坐标信息;以薄壁曲面件刀触面为研究对象,定义自由态与夹紧态下薄壁曲面件加工刀触点间映射关系为f;
自由态时刀触点为Di=[xi,yi,zi,1],则自由态刀位点为:
Wi=Di+R·en(2)
其中,R为刀具半径,en为自由态刀触点单位法矢;
夹紧态时刀触点为:
则夹紧态刀位点为:
其中,为夹紧态刀触点单位法矢;
根据式(1)-(4),建立自由态与夹紧态下薄壁曲面件刀位点映射关系如下:
由此,夹紧态下薄壁曲面件加工轨迹生成的关键在于,其一求解自由态与夹紧态下薄壁曲面件加工刀触点间映射关系f,其二将夹紧态刀触点沿法矢生成夹紧态下薄壁曲面件刀位点;
步骤2,自由态与夹紧态下薄壁曲面件加工刀触点间映射关系求解
采用有限元方法计算夹紧态薄壁曲面件加工刀触点,薄壁曲面件变形前后有限元网格单元和节点之间存在一一对应关系;映射法则如下:映射前后刀触点所在的有限元网格单元相同,映射前后刀触点在各自有限元网格单元中的位置相似;
自由态薄壁曲面件加工刀触点所在的有限元网格单元位置分三种情况:刀触点在有限元网格单元内;刀触点在有限元网格单元边线上;刀触点为有限元网格节点;
根据自由态薄壁曲面件加工刀触点所在位置,分别对每种情况下的夹紧态加工刀触点坐标计算如下:
1)刀触点在有限元网格单元内
根据有限元网格单元类型确定刀触点在有限元网格单元中的位置;采用三角形单元或四边形单元,面积坐标表示有限元网格单元内点的位置;
自由态刀触点e在有向三角形单元Δijk中,它的面积坐标表示为:
其中,i,j,k为自由态有限元三角形网格节点,按逆时针排列,约定顶点按逆时针方向排列的三角形面积为正,顶点按顺时针方向排列的三角形面积为负,S为三角形单元Δijk的面积,Sejk,Seki,Seij分别为有向三角形单元Δejk,Δeki,Δeij的面积;
有限元三角形单元Δijk所对应的有限元映射单元为有向三角形Δi*j*k*,由于Li,Lj,Lk不变,则夹紧态刀触点e*的坐标为:
其中,i*,j*,k*为夹紧态有限元三角形网格节点,S*为有向三角形单元Δi*j*k*的面积,分别为有向三角形单元Δe*j*k*,Δe*k*i*,Δe*i*j*的面积;
定义有向三角形第一个顶点A1坐标为(x1,y1,z1),第二个顶点A2坐标为(x2,y2,z2),第三个顶点A3坐标为(x3,y3,z3),则该三角形面积S有向可用如下行列式表示:
其内任意一点A(x,y,z)与另外两个顶点A2,A3构成的三角形面积为:
由此,自由态刀触点e的坐标为(xe,ye,ze,1),与之相关联的节点有限元三角形网格节点i,j,k的坐标为(xi,yi,zi,1),(xj,yj,zj,1),(xk,yk,zk,1),则自由态下刀触点面积坐标为:
由于面积坐标与直角坐标换算关系为:
由此可得:
式中,夹紧态刀触点e*的坐标为与之相关联的节点有限元三角形网格节点i*,j*,k*的坐标为
薄壁曲面件夹紧态下,与自由态刀触点相关联的节点有限元三角形网格节点i,j,k的变形矩阵Ti,Tj,Tk为:
式中,(Δxi,Δyi,Δzi),(Δxj,Δyj,Δzj),(Δxk,Δyk,Δzk)分别为节点i,j,k的变形量;则有:
e*=[Li Lj Lk]·[i j k]T·[Ti Tj Tk] (13)
其中,[Li Lj Lk]=h(e,i,j,k),h表示矩阵中数值Li,Lj,Lk与自由态下刀触点e、有限元三角形网格节点i,j,k间的函数关系;
由此,求得夹紧态与自由态下薄壁曲面件加工刀触点间映射关系f为:
f=h(e,i,j,k)·[i j k]T·[Ti Tj Tk]·e-1 (14)
2)刀触点在有限元网格单元边线上
自由态刀触点e落在有限元网格单元边线节点i,j之间,夹紧态有限元网格单元边线节点i*,j*为对应的变形节点,e和i,j之间存在一个比例系数λ:
其中,(xe,ye,ze,1)为自由态刀触点e的坐标,(xi,yi,zi,1),(xj,yj,zj,1)分别为有限元网格单元节点i,j的坐标;
由于比例系数λ不变,得到夹紧态刀触点e*的坐标为:
即:
其中,为夹紧态刀触点e*的坐标,分别为有限元网格单元节点i*,j*的坐标;
薄壁曲面件夹紧态下,与自由态刀触点相关联的节点有限元三角形网格节点i,j的变形矩阵分别为Ti,Tj,E为单位矩阵:
式中,(Δxi,Δyi,Δzi),(Δxj,Δyj,Δzj)分别为节点i,j的变形量;
式(17)可化为:
化简得:
由此,求得夹紧态与自由态下薄壁曲面件加工刀触点间映射关系f为:
3)刀触点在网格单元节点上
自由态刀触点e落在有限元网格节点i上,夹紧态有限元网格节点i*为对应的变形节点,按照变形前后节点间的一一对应关系,得到夹紧态下刀触点e*的坐标:
e*=i*=i·T=e·Ti(21)
式中,(Δxi,Δyi,Δzi)为节点i的变形量;
由此,夹紧态与自由态下薄壁曲面件加工刀触点间映射关系f为:
f=T(22)
步骤3,夹紧态下薄壁曲面件刀触点法矢计算
刀触点网格节点中某一节点P的几何欧拉距离最短的m个数据点为P的m近邻或P的m邻域,兼顾精度和效率取m=22-32,进行最小二乘曲面局部拟合求解夹紧态下薄壁曲面件刀触点法矢;
利用最小二乘法原理对夹紧态下薄壁曲面件加工刀触点附近局部区域进行曲面拟合,对于任意刀触点,它的法向量可以用该点及其邻域点拟合出的最小二乘曲面在该点的法矢替代;拟合的曲面方程为:
其中,a0,a1,a2,a3,a4,a5为拟合系数;利用夹紧态刀触点p=1,2,3,...,m拟合计算出式(23)的曲面方程,则使函数值W最小,W为:
为此,应满足:
即满足线性方程组:
解方程组获得aq的值,将拟合曲面方程写成参数方程形式r(x,y,z):
求曲面r(x,y)的偏微分记作曲面各个刀触点处单位法矢为:
步骤4,夹紧态下薄壁曲面件加工轨迹再规划
由自由态下薄壁曲面件刀触点坐标,联合自由态与夹紧态下薄壁曲面件加工刀触点间映射关系f以及求解得到的夹紧态下薄壁曲面件刀触点法矢获得夹紧态下薄壁曲面件加工刀位点;结合自由态下刀轨连接信息,通过后处理,输出机床能够识别的加工文件,实现减少薄壁曲面件加工完成后装夹回弹变形引起的加工误差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810222492.7A CN108416153B (zh) | 2018-03-19 | 2018-03-19 | 自由态与夹紧态下薄壁曲面件加工轨迹映射方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810222492.7A CN108416153B (zh) | 2018-03-19 | 2018-03-19 | 自由态与夹紧态下薄壁曲面件加工轨迹映射方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108416153A true CN108416153A (zh) | 2018-08-17 |
CN108416153B CN108416153B (zh) | 2019-11-29 |
Family
ID=63132034
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810222492.7A Active CN108416153B (zh) | 2018-03-19 | 2018-03-19 | 自由态与夹紧态下薄壁曲面件加工轨迹映射方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108416153B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109448122A (zh) * | 2018-10-10 | 2019-03-08 | 天津大学 | 对包含结构面的非连续体进行变形量测的方法及应用 |
CN110497088A (zh) * | 2019-08-31 | 2019-11-26 | 大连理工大学 | 基于曲面映射的柔性共形天线激光加工误差控制方法 |
CN111581874A (zh) * | 2020-06-28 | 2020-08-25 | 大连理工大学 | 一种薄壁曲面表层图案激光加工轨迹求解方法 |
CN114131426A (zh) * | 2021-11-09 | 2022-03-04 | 中国人民解放军国防科技大学 | 基于快速伺服刀具的弱刚度反射镜加工方法、系统及介质 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020002420A1 (en) * | 1999-09-20 | 2002-01-03 | Junichi Hirai | Numerically controlled curved surface machining unit |
CN101840453A (zh) * | 2010-04-22 | 2010-09-22 | 西北工业大学 | 薄壁曲面结构的有限元网格生成方法 |
CN102081693A (zh) * | 2011-01-11 | 2011-06-01 | 大连理工大学 | 薄壁复杂曲面零件三维加工预处理方法 |
CN103558808A (zh) * | 2013-09-28 | 2014-02-05 | 大连理工大学 | 复杂曲面五轴数控加工刀矢的运动学控制方法 |
CN103592891A (zh) * | 2013-09-28 | 2014-02-19 | 大连理工大学 | 运动学约束的复杂曲面五轴数控加工刀矢光顺方法 |
CN104615084A (zh) * | 2015-01-23 | 2015-05-13 | 大连理工大学 | 加工进给速度优化的刀轨曲线轮廓误差补偿方法 |
CN105739432A (zh) * | 2016-03-17 | 2016-07-06 | 浙江大学 | 基于改进型Butterfly细分的网格自由曲面环形刀具轨迹规划方法 |
CN106125666A (zh) * | 2016-07-28 | 2016-11-16 | 大连理工大学 | 以切削力波动为约束的曲面加工刀具轨迹规划方法 |
CN107728577A (zh) * | 2017-09-28 | 2018-02-23 | 大连理工大学 | 基于薄壁曲面加工变形的瞬时切削量规划方法 |
-
2018
- 2018-03-19 CN CN201810222492.7A patent/CN108416153B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020002420A1 (en) * | 1999-09-20 | 2002-01-03 | Junichi Hirai | Numerically controlled curved surface machining unit |
CN101840453A (zh) * | 2010-04-22 | 2010-09-22 | 西北工业大学 | 薄壁曲面结构的有限元网格生成方法 |
CN102081693A (zh) * | 2011-01-11 | 2011-06-01 | 大连理工大学 | 薄壁复杂曲面零件三维加工预处理方法 |
CN103558808A (zh) * | 2013-09-28 | 2014-02-05 | 大连理工大学 | 复杂曲面五轴数控加工刀矢的运动学控制方法 |
CN103592891A (zh) * | 2013-09-28 | 2014-02-19 | 大连理工大学 | 运动学约束的复杂曲面五轴数控加工刀矢光顺方法 |
CN104615084A (zh) * | 2015-01-23 | 2015-05-13 | 大连理工大学 | 加工进给速度优化的刀轨曲线轮廓误差补偿方法 |
CN105739432A (zh) * | 2016-03-17 | 2016-07-06 | 浙江大学 | 基于改进型Butterfly细分的网格自由曲面环形刀具轨迹规划方法 |
CN106125666A (zh) * | 2016-07-28 | 2016-11-16 | 大连理工大学 | 以切削力波动为约束的曲面加工刀具轨迹规划方法 |
CN107728577A (zh) * | 2017-09-28 | 2018-02-23 | 大连理工大学 | 基于薄壁曲面加工变形的瞬时切削量规划方法 |
Non-Patent Citations (1)
Title |
---|
张振金等: "自由曲面数控加工刀位轨迹与加工精度的研究", 《抚顺石油学院学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109448122A (zh) * | 2018-10-10 | 2019-03-08 | 天津大学 | 对包含结构面的非连续体进行变形量测的方法及应用 |
CN109448122B (zh) * | 2018-10-10 | 2023-02-28 | 天津大学 | 对包含结构面的非连续体进行变形量测的方法及应用 |
CN110497088A (zh) * | 2019-08-31 | 2019-11-26 | 大连理工大学 | 基于曲面映射的柔性共形天线激光加工误差控制方法 |
CN111581874A (zh) * | 2020-06-28 | 2020-08-25 | 大连理工大学 | 一种薄壁曲面表层图案激光加工轨迹求解方法 |
CN114131426A (zh) * | 2021-11-09 | 2022-03-04 | 中国人民解放军国防科技大学 | 基于快速伺服刀具的弱刚度反射镜加工方法、系统及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN108416153B (zh) | 2019-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108416153B (zh) | 自由态与夹紧态下薄壁曲面件加工轨迹映射方法 | |
CN104289748B (zh) | 一种大型薄壁蒙皮自适应等壁厚铣削系统及其加工方法 | |
CN104252566B (zh) | 一种箱体结构的简化及装夹变形仿真分析方法 | |
CN101482403A (zh) | 一种基于有限元方法的确定薄板冲压件测量点云在回弹前位置的方法 | |
Liu et al. | Constant scallop-height tool path generation for three-axis discrete data points machining | |
CN105108215A (zh) | 一种自由曲面微细铣削让刀误差预测及补偿方法 | |
CN108508848A (zh) | 一种基于插补数据的铣削加工轮廓误差的评估方法 | |
CN110837715B (zh) | 一种基于逆向工程技术的复杂曲面加工误差补偿方法 | |
Sun et al. | Application of FEM to tool design for electrochemical machining freeform surface | |
Wu et al. | Adaptive location of repaired blade for multi-axis milling | |
JP2007058508A (ja) | 形状測定に基づく構造解析システム及び解析方法 | |
Jia et al. | Feed speed scheduling method for parts with rapidly varied geometric feature based on drive constraint of NC machine tool | |
Huang et al. | A systematic approach for online minimizing volume difference of multiple chambers in machining processes based on high-definition metrology | |
Sun et al. | Iso-planar feed vector-fields-based streamline tool path generation for five-axis compound surface machining with torus-end cutters | |
Cai et al. | An improved cutting force prediction model in the milling process with a multi-blade face milling cutter based on FEM and NURBS | |
JP3344649B2 (ja) | 自由曲面の形状誤差評価方法 | |
Lu et al. | Tool path generation for turbine blades machining with twin tool | |
CN110497088B (zh) | 基于曲面映射的柔性共形天线激光加工误差控制方法 | |
Jiang et al. | A non-uniform allowance allocation method based on interim state stiffness of machining features for NC programming of structural parts | |
Kim et al. | Material pixel-based process planning for layered manufacturing of heterogeneous objects | |
CN109725595A (zh) | 工件的加工路径的补偿方法、加工方法及工件 | |
CN105373672A (zh) | 一种四边形单元实时加密-减密方法 | |
JP2007193552A (ja) | 面モデルの作成装置と作成方法 | |
CN109635348B (zh) | 一种空间复杂曲面磨削温度数值仿真分析方法 | |
CN105160133A (zh) | 基于stl三角片顶点等距偏置的渐进成形支撑体生成法 |
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 |