首页 / 专利库 / 物理 / 连续介质力学 / 一种滚珠丝杠副定位精度衰退预测方法

一种滚珠丝杠定位精度衰退预测方法

阅读:758发布:2020-06-23

专利汇可以提供一种滚珠丝杠定位精度衰退预测方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种滚珠 丝杠 副 定位 精度 衰退预测方法,包括:滚珠丝杠副定位精度衰退分析、准连续介质形变影响的滚珠丝杠副应变能分布分析、不同进给工况下的滚珠丝杠副定位精度衰退模型的建立与 预测分析 ,基于准连续介质理论,根据滚珠丝杠副的磨损表面的微观 变形 特征确定其应变能分布,考虑滚珠丝杠副的不同进给工况,建立了滚珠丝杠副的指数磨损模型,根据宏‑微多尺度方法,对滚珠丝杠副的定位精度保持性的有效寿命进行了预测。本发明建立一套实现精度、精度保持性协调评价,融合贯通的精度评价体系,具有承前启后的重要意义。,下面是一种滚珠丝杠定位精度衰退预测方法专利的具体信息内容。

1.一种滚珠丝杠定位精度衰退预测方法,其特征在于,实现该方法的步骤包括如下:
步骤一 滚珠丝杠副定位精度衰退分析
滚珠丝杠副精度的衰退主要是由磨损造成的;滚珠丝杠副存在定位精度误差、度误差以及平行度误差等由于磨损,导致其精度呈现衰退特性;对滚珠丝杠副的误差做出分析:
滚珠丝杠副经过磨损后,直接导致平面内以及垂直面内的定位精度误差增大;
滚珠丝杠副的磨损导致其在X、Y、Z三个方向的精度衰退;假设处于滚道中的滚珠个数为m个,则滚珠在丝杠滚道的弧长ls1-m、螺母滚道中的弧长ln1-m分别为:
其中,βsm是所有滚珠在丝杠滚道上所对应的圆心角,βnm是所有滚珠在螺母滚道上所对应的圆心角,X′(β)、Y′(β)、Z'(β)滚珠球心参数方程X(β)、Y(β)、Z(β)的一阶导数;
设滚珠所在丝杠滚道中的初始位置为xsi,滚珠与丝杠滚道所有接触微凸体沿丝杠螺旋方向回归直线的平均偏差为ebs(xsi);设滚珠所在螺母滚道中的初始位置为xni,滚珠与螺母滚道所有接触微凸体沿螺母螺旋方向回归直线的平均偏差为ebn(xni),则滚珠丝杠副在进给方向上的定位精度测量值gfd(xi)表示为:
设滚珠丝杠副在进给方向上的定位精度误差为Gfd,Gfd表示全体微凸点偏向性分布的正负极大值之和,反映了滚珠丝杠副进给精度的衰退过程与滚道表面形貌的特征、分布情况密切相关;
步骤二 准连续介质形变影响的滚珠丝杠副应变能分布
滚珠丝杠副的磨损直接影响滚道表面形貌特征,考虑磨损的微观过程,根据滚珠丝杠副磨损表面微观形变特性确定其应变能分布状况;采用准连续介质方法,从而建立滚珠丝杠副微观磨损模型;
当接触表面受到预加载荷时,微观上分为局部区域与非局部区域;
处于受载荷变形的核心区,把所有的原子都选作为代表原子,处于变形梯度较小的区域,选取若干原子作为代表原子;
根据准连续介质法,则滚珠丝杠副准连续介质系统的能量表达式为:
其中, 为准连续介质系统耦合局部连续部分与非局部原子部分的静载应能的总和, 为准连续介质系统得到外部负载所做功的总和, 为准连续介质系统局部连续部分的线弹性应力能的总和,由下式计算得到:
式中,Ncn为连续区域代表原子数,Λφ为单元cn的体积, 为单元应变能密度,为准连续介质系统局部连续区域的总变形梯度;
其中, 为准连续介质系统非局部原子部分的原子间相互作用势能的总和,由下式计算得到:
式中,Nnla为非局部原子区域代表原子数,Wnla为权函数,Enla为SW势能;Enla与原子k上的力的关系为:
式中,rrw为原子k相邻范围内的所有原子,其大小取为原子半径的2-3倍;
滚珠丝杠副在服役过程中,由于外界工况的作用导致其发生磨损;根据对准连续介质系统局部与非局部的判定原理,认为滚珠丝杠副磨损过程中,准连续介质系统局部区域的总变形梯度的变化率为定值δφ,通过 计算得到;
ψ为形变位移量,χ为原子区域的比例系数,χψ为非局部原子区域的范围, 为局部连续部分单元的最小变形梯度值, 为非局部原子部分单元的最大变形梯度,对于局部连续部分单元的总变形梯度变化率较小且为定值,设其变化的斜率为tanβ;
其中,Q为非局部原子部分指数方程的系数,nχ为非局部原子部分指数方程的原指数;
当预加载荷卸载后,滚珠丝杠副在准连续介质系统内的形变位移量不能完全恢复;η为虚延长率,ηχψ为虚延长后非局部原子区域的范围,结合二项式定理对指数部分进行简化,可得:
其中,nη为非局部原子部分指数方程的虚延长比例指数,可结合量化理论中晶格弛豫理论求得;
考虑扩散过程的影响,认为原子扩散损失后,非局部原子区域的远端指数曲线回缩,k为扩散系数,kηχψ为扩散原子区域的范围,结合二项式定理对指数部分进行简化,可得:
其中,nk为原子部分指数方程的扩散比例指数;
滚珠丝杠副磨损会导致材料的损失,材料的损失会影响原子比例系数χ以及原子比例指数nχ,将连续区域的材料损失部分用虚延长的指数方程近似;
σψ为损失量,σ为连续部分损失比例,结合二项式定理对指数部分进行简化,可得:
其中,R为指数方程连续部分的常数项系数;
步骤三.不同进给工况下的滚珠丝杠副定位精度衰退模型的建立与预测分析
通过对准连续介质系统内变形梯度分布的分析,设准连续介质系统原子部分垂直于σ的截面积为Aal,连续部分垂直于σ的截面积为λAal,λ为比例系数;基于耦合特性以及扩散影响,考虑磨损导致滚珠丝杠副准连续介质系统内的受损,可得单向准连续介质系统连续部分与原子部分的简化能量表达式:
其中,nt为磨损指数,由下式计算得到:
nt=nχ+nη+nk     (14)
根据准连续介质力学原理和虚功原理,准连续介质系统单元平衡点位移由下式确定:
由式(15)得:
基于二项式定理化简并通分后得:
其中, 变化相对较小,由于连续部分损失比例σ相对较大,但是线弹性区域能量的损失几乎不影响原子区域的指数分布,主要影响常数项系数; 中χ相
对其他数量级;综上所述,基于准连续介质原理以及指数近似假设,滚珠丝杠副在服役过程中,其磨损量随着累计磨损次数呈现指数变化;假设滚珠丝杠副受到平稳接触磨损的次数为C,可得滚珠丝杠副连续介质系统局部磨损量Δdb与接触磨损次数Cb之间的关系式:
其中,Qb为滚珠丝杠副脆性受载情况下的指数部分原子区域磨损系数,Rb为滚珠丝杠副脆性受载情况下的指数部分连续区域磨损系数,Kb为滚珠丝杠副脆性受载情况下的弹性部分磨损系数,nχb为非局部原子部分指数方程的原指数,nηb为非局部原子部分指数方程的虚延长比例指数,nkb为原子部分指数方程的扩散比例指数;
当滚珠丝杠副处于韧性磨损阶段过程中,则不存在原子部分的虚延长率,即不考虑非局部原子部分指数方程的虚延长比例指数nη,可得滚珠丝杠副的磨损深度Δdt与磨损次数Ct之间的关系式:
其中,Qt为滚珠丝杠副韧性受载情况下的指数部分原子区域磨损系数,Kt为滚珠丝杠副韧性受载情况下的弹性部分磨损系数,nχt为非局部原子部分指数方程的原指数,nkt为原子部分指数方程的扩散比例指数;
滚珠丝杠副在服役过程中,丝杠与螺母均会发生磨损;由于机床加工工件大小的不同,导致滚珠丝杠副每次的进给量不同;
假设滚珠丝杠副在服役过程中的进给总数为N,螺母相对于丝杠的最小进给量为fmin,螺母相对于丝杠的最大进给量为fmax,进给量f服从正态分布;引入进给量系数ξ,0≤ξ≤1,ξ是一个随机变量,用来衡量滚珠丝杠副进给量的不同程度;
滚珠丝杠副进给方式影响其磨损量的分布,根据实际服役工况,滚珠丝杠副的进给量系数ξ基本服从正态分布、负偏态分布以及正偏态分布三种状态;正态分布表示实际工况中较小与较大进给距离所占的比重小,负偏态分布表示实际工况偏向于大进给量,正偏态分布表示实际工况偏向于小进给量;
假设滚珠丝杠副在相同的服役时间内按照三种方式的进给总量分别为Ωa、Ωb、Ωc;
滚珠丝杠副在服役过程中,不同的进给方式导致其微凸体接触次数不同,从而导致其磨损特性不同;将每次相互接触微凸体的应变能损失计算过程进行简化,结合式(4)、(18)与(19),得到滚珠丝杠副进给方向上简化的直线度测量点的表达式
sfd(xi)≈Λxiw        (20)
式(20)中,Λ为滚珠丝杠副不同受载情况下的磨损系数,w为滚珠丝杠副不同受载情况下的磨损比例指数,当滚珠丝杠副处于脆性受载时,有:
当滚珠丝杠副处于韧性受载时,有:
随着滚珠丝杠副的服役与进给方式,处于不同位置微凸体的接触次数会呈现不同;微凸体的接触次数不断累计,处于不同位置的微凸体呈现不同的磨损程度;结合式(20),根据滚珠丝杠副的磨损深度与进给方向上直线度的几何关系,滚珠丝杠副的进给精度保持性通用公式表示如下:
Sfd≈Sfd0(xi)tw      (23)
其中,Sfd0为滚珠丝杠副的初始进给精度值,Sfd为滚珠丝杠副服役一段时间后的实时进给精度值,t为滚珠丝杠副的服役时间,w为滚珠丝杠副的实时磨损指数。

说明书全文

一种滚珠丝杠定位精度衰退预测方法

技术领域

[0001] 本发明提供了一种基于宏-微多尺度的机床滚珠丝杠副定位精度衰退预测方法,属于数控机床滚珠丝杠副技术领域。

背景技术

[0002] 滚珠丝杠副在数控机床等设备的传递动与运动方式转化中起着关键性的作用,滚珠丝杠副的定位精度对数控机床的加工精度影响很大。滚珠丝杠副在服役过程中会发生磨损,磨损是导致其定位精度衰退的主要因素,也是导致定位精度衰退的直接因素。一种新型的考虑滚珠丝杠副在不同进给工况下的定位精度衰退预测方法被发明了,基于准连续介质理论,根据滚珠丝杠副的磨损表面的微观变形特征确定其应变能分布,考虑滚珠丝杠副的不同进给工况,建立了滚珠丝杠副的指数磨损模型,根据宏-微多尺度方法,对滚珠丝杠副的定位精度保持性的有效寿命进行了预测。
[0003] 本发明对发掘磨损对精度保持性的作用影响机理,建立一套实现精度、精度保持性协调评价,融合贯通的精度评价体系,具有承前启后的重要意义。

发明内容

[0004] 本发明目的是提供一种滚珠丝杠副定位精度衰退预测方法,包括:滚珠丝杠副定位精度衰退分析、准连续介质形变影响的滚珠丝杠副应变能分布分析、不同进给工况下的滚珠丝杠副定位精度衰退模型的建立与预测分析
[0005] 本发明采用的技术方案为一种滚珠丝杠副定位精度衰退预测方法,实现该方法的步骤包括如下:
[0006] S1、滚珠丝杠副定位精度衰退分析;
[0007] S2、基于准连续介质理论,考虑准连续介质形变影响下,分析滚珠丝杠副磨损面的应变能分布;
[0008] S3、不同进给工况下的滚珠丝杠副定位精度衰退模型的建立与预测分析。附图说明
[0009] 图1滚珠丝杠副磨损状态示意图。
[0010] 图2准连续介质方法示意图。
[0011] 图3准连续介质系统变形梯度示意图。
[0012] 图4准连续介质系统内的耦合变形梯度示意图。
[0013] 图5考虑扩散的准连续介质系统内变形梯度示意图。
[0014] 图6准连续介质系统内的受磨损变形梯度示意图。
[0015] 图7滚珠丝杠副进给示意图。
[0016] 图8ξ的概率密度分布函数。
[0017] 图9进给总量与进给量系数关系。

具体实施方式

[0018] 以下结合附图和实施例对本发明进行详细说明。
[0019] 一种滚珠丝杠副定位精度衰退预测方法,实现该方法的步骤包括如下:
[0020] 步骤一 滚珠丝杠副定位精度衰退分析
[0021] 滚珠丝杠副精度的衰退主要是由磨损造成的。滚珠丝杠副存在定位精度误差、度误差以及平行度误差等由于磨损,导致其精度呈现衰退特性。对滚珠丝杠副的误差做出分析:滚珠丝杠副经过磨损后,直接导致平面内以及垂直面内的定位精度误差增大。滚珠丝杠副由于磨损导致的接触状态如图1所示。
[0022] 滚珠丝杠副的磨损导致其在X、Y、Z三个方向的精度衰退。假设处于滚道中的滚珠个数为m个,则滚珠在丝杠滚道的弧长ls1-m、螺母滚道中的弧长ln1-m分别为:
[0023]
[0024]
[0025] 其中,βsm是所有滚珠在丝杠滚道上所对应的圆心角,βnm是所有滚珠在螺母滚道上所对应的圆心角,X′(β)、Y′(β)、Z'(β)滚珠球心参数方程X(β)、Y(β)、Z(β)的一阶导数。
[0026]
[0027] 设滚珠所在丝杠滚道中的初始位置为xsi,滚珠与丝杠滚道所有接触微凸体沿丝杠螺旋方向回归直线的平均偏差为ebs(xsi);设滚珠所在螺母滚道中的初始位置为xni,滚珠与螺母滚道所有接触微凸体沿螺母螺旋方向回归直线的平均偏差为ebn(xni),则滚珠丝杠副在进给方向上的定位精度测量值gfd(xi)表示为:
[0028]
[0029] 设滚珠丝杠副在进给方向上的定位精度误差为Gfd,Gfd表示全体微凸点偏向性分布的正负极大值之和,反映了滚珠丝杠副进给精度的衰退过程与滚道表面形貌的特征、分布情况密切相关。
[0030] 步骤二 准连续介质形变影响的滚珠丝杠副应变能分布
[0031] 滚珠丝杠副的磨损直接影响滚道表面形貌特征,考虑磨损的微观过程,根据滚珠丝杠副磨损表面微观形变特性确定其应变能分布状况。采用准连续介质方法,从而建立滚珠丝杠副微观磨损模型。
[0032] 当接触表面受到预加载荷时,微观上分为局部区域与非局部区域。如图2所示为准连续介质方法示意图。
[0033] 在图2中,处于受载荷变形的核心区,把所有的原子都选作为代表原子,处于变形梯度较小的区域,选取若干原子作为代表原子。
[0034] 根据准连续介质法,则滚珠丝杠副准连续介质系统的能量表达式为:
[0035]
[0036] 其中, 为准连续介质系统耦合局部连续部分与非局部原子部分的静载应力能的总和, 为准连续介质系统得到外部负载所做功的总和, 为准连续介质系统局部连续部分的线弹性应力能的总和,由下式计算得到:
[0037]
[0038] 式中,Ncn为连续区域代表原子数,Λφ为单元cn的体积, 为单元应变能密度, 为准连续介质系统局部连续区域的总变形梯度。
[0039] 其中, 为准连续介质系统非局部原子部分的原子间相互作用势能的总和,由下式计算得到:
[0040]
[0041] 式中,Nnla为非局部原子区域代表原子数,Wnla为权函数,Enla为SW(Stillinger-weber)势能。Enla与原子k上的力的关系为:
[0042]
[0043] 式中,rrw为原子k相邻范围内的所有原子,其大小一般取为原子半径的2-3倍。
[0044] 滚珠丝杠副在服役过程中,由于外界工况的作用导致其发生磨损。根据对准连续介质系统局部与非局部的判定原理,认为滚珠丝杠副磨损过程中,准连续介质系统局部区域的总变形梯度的变化率为定值δφ,通过 计算得到。假设当滚珠丝杠副受载后,其准连续介质系统内的单元变形梯度分布如图3所示。
[0045] 图3中,ψ为形变位移量,χ为原子区域的比例系数,χψ为非局部原子区域的范围,为局部连续部分单元的最小变形梯度值, 为非局部原子部分单元的最大变形梯度,对于局部连续部分单元的总变形梯度变化率较小且为定值,设其变化的斜率为tanβ。由图3可知:
[0046]
[0047] 其中,Q为非局部原子部分指数方程的系数,nχ为非局部原子部分指数方程的原指数。
[0048] 当预加载荷卸载后,滚珠丝杠副在准连续介质系统内的形变位移量不能完全恢复。根据局部与非局部区域判定定理,由图3可得耦合部分梯度变化率大于局部连续区域的梯度变化率,将原子区域延长至梯度变化率为0处,如图4所示。
[0049] 图4中,η为虚延长率,ηχψ为虚延长后非局部原子区域的范围,结合二项式定理对指数部分进行简化,可得:
[0050]
[0051] 其中,nη为非局部原子部分指数方程的虚延长比例指数,可结合量化理论中晶格弛豫理论求得。
[0052] 考虑扩散过程的影响,认为原子扩散损失后,非局部原子区域的远端指数曲线回缩,如图5所示。
[0053] 图5中,k为扩散系数,kηχψ为扩散原子区域的范围,结合二项式定理对指数部分进行简化,可得:
[0054]
[0055] 其中,nk为原子部分指数方程的扩散比例指数。
[0056] 滚珠丝杠副磨损会导致材料的损失,材料的损失会影响原子比例系数χ以及原子比例指数nχ,将连续区域的材料损失部分用虚延长的指数方程近似,如图6所示。
[0057] 图6中,σψ为损失量,σ为连续部分损失比例,结合二项式定理对指数部分进行简化,可得:
[0058]
[0059] 其中,R为指数方程连续部分的常数项系数。
[0060] 步骤三.不同进给工况下的滚珠丝杠副定位精度衰退模型的建立与预测分析[0061] 通过对准连续介质系统内变形梯度分布的分析,设准连续介质系统原子部分垂直于σ的截面积为Aal,连续部分垂直于σ的截面积为λAal,λ为比例系数。基于耦合特性以及扩散影响,考虑磨损导致滚珠丝杠副准连续介质系统内的受损,可得单向准连续介质系统连续部分与原子部分的简化能量表达式:
[0062]
[0063] 其中,nt为磨损指数,由下式计算得到:
[0064] nt=nχ+nη+nk  (14)
[0065] 根据准连续介质力学原理和虚功原理,准连续介质系统单元平衡点位移由下式确定:
[0066]
[0067] 由式(15)得:
[0068]
[0069] 基于二项式定理化简并通分后得:
[0070]
[0071] 其中, 变化相对较小,由于连续部分损失比例σ相对较大,但是线弹性区域能量的损失几乎不影响原子区域的指数分布,主要影响常数项系数。中χ相对其他数量级。综上所述,基于准连续介质原理以及指数近似假设,滚珠丝杠副在服役过程中,其磨损量随着累计磨损次数呈现指数变化。假设滚珠丝杠副受到平稳接触磨损的次数为C,可得滚珠丝杠副连续介质系统局部磨损量Δdb与接触磨损次数Cb之间的关系式:
[0072]
[0073] 其中,Qb为滚珠丝杠副脆性受载情况下的指数部分原子区域磨损系数,Rb为滚珠丝杠副脆性受载情况下的指数部分连续区域磨损系数,Kb为滚珠丝杠副脆性受载情况下的弹性部分磨损系数,nχb为非局部原子部分指数方程的原指数,nηb为非局部原子部分指数方程的虚延长比例指数,nkb为原子部分指数方程的扩散比例指数。
[0074] 当滚珠丝杠副处于韧性磨损阶段过程中,则不存在原子部分的虚延长率,即不考虑非局部原子部分指数方程的虚延长比例指数nη,可得滚珠丝杠副的磨损深度Δdt与磨损次数Ct之间的关系式:
[0075]
[0076] 其中,Qt为滚珠丝杠副韧性受载情况下的指数部分原子区域磨损系数,Kt为滚珠丝杠副韧性受载情况下的弹性部分磨损系数,nχt为非局部原子部分指数方程的原指数,nkt为原子部分指数方程的扩散比例指数。
[0077] 滚珠丝杠副在服役过程中,丝杠与螺母均会发生磨损。由于机床加工工件大小的不同,导致滚珠丝杠副每次的进给量不同。如图7所示为螺母相对于丝杠的进给量示意图。
[0078] 由图7:假设滚珠丝杠副在服役过程中的进给总数为N,螺母相对于丝杠的最小进给量为fmin,螺母相对于丝杠的最大进给量为fmax,进给量f服从正态分布。引入进给量系数ξ,0≤ξ≤1,ξ是一个随机变量,用来衡量滚珠丝杠副进给量的不同程度。
[0079] 滚珠丝杠副进给方式影响其磨损量的分布,根据实际服役工况,滚珠丝杠副的进给量系数ξ基本服从正态分布(如图8中的(a)所示)、负偏态分布(如图8中的(b)所示)以及正偏态分布(如图8中的(c)所示)三种状态。正态分布表示实际工况中较小与较大进给距离所占的比重小,负偏态分布表示实际工况偏向于大进给量,正偏态分布表示实际工况偏向于小进给量。
[0080] 假设滚珠丝杠副在相同的服役时间内按照三种方式的进给总量分别为Ωa、Ωb、Ωc,则进给总量Ωa与正态分布的进给量系数ξa的关系如图9中的(a)所示,进给总量Ωb与负偏态分布的进给量系数ξb的关系如图9中的(b)所示,进给总量Ωc与正偏态分布的进给量系数ξc的关系如图9中的(c)所示。
[0081] 滚珠丝杠副在服役过程中,结合图9,不同的进给方式导致其微凸体接触次数不同,从而导致其磨损特性不同。将每次相互接触微凸体的应变能损失计算过程进行简化,结合式(4)、(18)与(19),得到滚珠丝杠副进给方向上简化的直线度测量点的表达式[0082] sfd(xi)≈Λxiw  (20)
[0083] 式(20)中,Λ为滚珠丝杠副不同受载情况下的磨损系数,w为滚珠丝杠副不同受载情况下的磨损比例指数,当滚珠丝杠副处于脆性受载时,有:
[0084]
[0085] 当滚珠丝杠副处于韧性受载时,有:
[0086]
[0087] 随着滚珠丝杠副的服役与进给方式,处于不同位置微凸体的接触次数会呈现不同。微凸体的接触次数不断累计,处于不同位置的微凸体呈现不同的磨损程度。结合图1与式(20),根据滚珠丝杠副的磨损深度与进给方向上直线度的几何关系,滚珠丝杠副的进给精度保持性通用公式表示如下:
[0088] Sfd≈Sfd0(xi)tw  (23)
[0089] 其中,Sfd0为滚珠丝杠副的初始进给精度值,Sfd为滚珠丝杠副服役一段时间后的实时进给精度值,t为滚珠丝杠副的服役时间,w为滚珠丝杠副的实时磨损指数。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈