CN115660224B - 一种滚柱直线导轨副的预紧拖动力预测方法 - Google Patents
一种滚柱直线导轨副的预紧拖动力预测方法 Download PDFInfo
- Publication number
- CN115660224B CN115660224B CN202211587676.6A CN202211587676A CN115660224B CN 115660224 B CN115660224 B CN 115660224B CN 202211587676 A CN202211587676 A CN 202211587676A CN 115660224 B CN115660224 B CN 115660224B
- Authority
- CN
- China
- Prior art keywords
- tightening
- guide rail
- linear guide
- dragging force
- roller linear
- 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 91
- 238000012549 training Methods 0.000 claims abstract description 41
- 238000004458 analytical method Methods 0.000 claims abstract description 9
- 238000012360 testing method Methods 0.000 claims description 27
- 238000009826 distribution Methods 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000001914 filtration Methods 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 8
- 238000001228 spectrum Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 5
- 238000002474 experimental method Methods 0.000 claims description 4
- 101100056815 Methanocaldococcus jannaschii (strain ATCC 43067 / DSM 2661 / JAL-1 / JCM 10045 / NBRC 100440) artD gene Proteins 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 230000006870 function Effects 0.000 description 22
- 230000015556 catabolic process Effects 0.000 description 11
- 238000006731 degradation reaction Methods 0.000 description 11
- 238000004364 calculation method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000012876 topography Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000002790 cross-validation Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 239000002243 precursor Substances 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 239000004519 grease Substances 0.000 description 1
- 230000001050 lubricating effect Effects 0.000 description 1
- 238000005461 lubrication Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000004901 spalling Methods 0.000 description 1
Images
Classifications
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种滚柱直线导轨副的预紧拖动力预测方法,该方法为:采集滚柱直线导轨副中导轨滚道多个区域对应的表面二维轮廓,同时测量滚柱直线导轨副预紧拖动力变化;使用几何参数法、单重分形法、递归分析法提取所得表面二维轮廓的四个特征参数;通过灰色关联度法计算所采集预紧拖动力与四个特征参数的关联度,对比关联度大小得到与预紧拖动力相关性最高的特征参数作为最优特征参数;将得到的最优特征参数和预紧拖动力进行多项式拟合,提取拟合曲线中的数值点构建训练集;建立高斯回归过程模型,将训练集导入模型进行训练得到最终高斯回归过程模型,用于滚柱直线导轨副的预紧拖动力预测。本发明适用于小样本、实用性好、精确率高。
Description
技术领域
本发明涉及预紧拖动力的监测技术领域,特别是一种滚柱直线导轨副的预紧拖动力预测方法。
背景技术
滚柱直线导轨副是广泛应用于机床和精密测量仪器的关键传动部件,具有承载能力高、传动精度高、寿命长等特点。由于在运行过程中摩擦摩损,滚柱直线导轨副的滚动体与滚道之间的过盈量会逐渐减小,直接导致预紧力减小进而影响其传动精度以及数控机床的加工精度。滚柱直线导轨副的预紧力无法直接测量,一般通过测量其预紧拖动力(摩擦力)来间接计算,因此检测滚柱直线导轨副在服役过程中的预紧拖动力退化情况十分重要。
目前关于预紧拖动力的监测多采用基于振动信号的特征提取的与识别方法,但是该方法易受到外界噪声的干扰。此外,由于滚柱直线导轨副的磨损实验需要大型加载设备作为支撑,以滚道形貌磨损过程表征滚柱直线导轨副预紧拖动力退化的研究较少,相关的研究主要集中在对材料磨损规律的探究,存在样本需求量大、实用性较差等缺点。
发明内容
本发明的目的在于提供一种有效的滚柱直线导轨副的预紧拖动力预测方法,适用于小样本的情况,且实用性好、精确率高。
实现本发明目的的技术解决方案为:一种滚柱直线导轨副的预紧拖动力预测方法,包括以下步骤:
步骤1,通过泰勒霍普森轮廓仪,采集滚柱直线导轨副中导轨滚道M个区域对应的表面二维轮廓,并对表面二维轮廓进行去除形状以及高斯滤波处理,同时测量滚柱直线导轨副预紧拖动力变化;M为大于5的自然数;
步骤2,使用几何参数法、单重分形法、递归分析法提取步骤1所得表面二维轮廓的特征参数,特征参数包括轮廓算数偏差R a 、最大峰谷距R t 、分形维数D、递归率R r ;
步骤3,基于所得到的四个特征参数,通过灰色关联度法计算步骤1所采集预紧拖动力与四个特征参数的关联度,对比关联度大小得到与预紧拖动力相关性最高的特征参数作为最优特征参数;
步骤4,将得到的前M-1个区域的最优特征参数和预紧拖动力进行M-1次多项式拟合,提取拟合曲线中的数值点,最优特征参数作为样本特征,预紧拖动力数值作为标签,构建训练集;
步骤5,建立高斯回归过程模型,将训练集导入高斯回归过程模型进行训练;
步骤6,将第M个区域的最优特征参数导入训练好的高斯回归过程模型,预测得到预紧拖动力数值,并与真实值进行比较,验证高斯回归过程模型准确性;
步骤7,将验证后的高斯回归过程模型用于滚柱直线导轨副的预紧拖动力预测。
本发明与现有技术相比,其显著优点为:(1)基于滚道形貌特征参数进行预紧拖动力预测,通多对小样本数据进行合理扩充,并基于高斯过程回归模型,可以实现对滚柱直线导轨副预紧拖动力的准确预测;(2)不仅可以克服传统根据电流信号、振动信号等手段识别检测滚柱直线导轨副预紧拖动力时存在的噪声干扰的缺点,而且适用于小样本的情况,实用性好、精确率高。
附图说明
图1为本发明滚柱直线导轨副的预紧拖动力预测方法的流程图。
图2为本发明实施例提供的滚柱直线导轨副加载运行试验台结构示意图。
图3为本发明实施例提供的滚柱直线导轨副滚道形貌采集装置结构示意图。
图4为本发明实施例提供的滚柱直线导轨副预紧拖动力测量试验台结构示意图。
图5为本发明实施例提供的基于滚道形貌特征参数的滚柱直线导轨副预紧拖动力预测方法的流程示意图。
图6为本发明实施例提供的滚柱直线导轨副疲劳剥落示意图。
图7为本发明实施例提供的四个特征参数与预紧拖动力退化率关联度示意图。
图8为本发明实施例提供的高斯回归过程模型训练示意图。
图9为本发明实施例提供的高斯回归模型预测示意图。
具体实施方式
由于滚柱直线导轨副的预紧力退化情况与滚道表面形貌参数的变化情况存在较强的映射关系,本发明提出了一种基于滚道表面形貌特征参数的滚柱直线导轨副预紧拖动力预测方法,用于解决滚柱直线导轨副预紧拖动力预测的技术问题。
结合图1,本发明一种滚柱直线导轨副的预紧拖动力预测方法,包括以下步骤:
步骤1,通过泰勒霍普森轮廓仪,采集滚柱直线导轨副中导轨滚道M个区域对应的表面二维轮廓,并对表面二维轮廓进行去除形状以及高斯滤波处理,同时测量滚柱直线导轨副预紧拖动力变化;M为大于5的自然数;
步骤2,使用几何参数法、单重分形法、递归分析法提取步骤1所得表面二维轮廓的特征参数,特征参数包括轮廓算数偏差R a 、最大峰谷距R t 、分形维数D、递归率R r ;
步骤3,基于所得到的四个特征参数,通过灰色关联度法计算步骤1所采集预紧拖动力与四个特征参数的关联度,对比关联度大小得到与预紧拖动力相关性最高的特征参数作为最优特征参数;
步骤4,将得到的前M-1个区域的最优特征参数和预紧拖动力进行M-1次多项式拟合,提取拟合曲线中的数值点,最优特征参数作为样本特征,预紧拖动力数值作为标签,构建训练集;
步骤5,建立高斯回归过程模型,将训练集导入高斯回归过程模型进行训练;
步骤6,将第M个区域的最优特征参数导入训练好的高斯回归过程模型,预测得到预紧拖动力数值,并与真实值进行比较,验证高斯回归过程模型准确性;
步骤7,将验证后的高斯回归过程模型用于滚柱直线导轨副的预紧拖动力预测。
作为一种具体示例,步骤1中M=6,首先在滚柱直线导轨副的导轨滚道上选取6个测量区域,每个测量区域选取3个位置进行标记;对滚柱直线导轨副进行摩擦磨损试验,每运转5km停机,使用泰勒霍普森轮廓仪采集滚柱直线导轨副中导轨滚道6个区域对应的18个表面二维轮廓,并对表面二维轮廓进行去除形状以及高斯滤波处理,同时利用预紧拖动力试验台测量滚柱直线导轨副预紧拖动力,得到预紧拖动力变化曲线。
作为一种具体示例,使用几何参数法提取步骤1所得表面二维轮廓的轮廓算数偏差R a 、最大峰谷距R t ,具体如下:
其中,L为采样长度,z(x)为测量轮廓高度函数,z max 为表面轮廓高度的最大值,z min 为表面轮廓高度最小值;z i 为第i个采样点的轮廓高度,i=1,2,…N,N为采样点数。
作为一种具体示例,使用单重分形法提取步骤1所得表面二维轮廓的分形维数D,具体如下:
W-M函数被用来模拟和表征表面二维轮廓,表达式为:
测量轮廓高度函数z(x)的功率谱函数S(ω)为:
定义z(x)的增量方差为结构函数,如下式所示:
其中,S(t)表示表面轮廓高度函数的结构函数,[·]2表示差方的算数平均值,t表示任意两数据之间的时间间隔;ω表示频率,即粗糙度波长的倒数;z(x)为测量轮廓高度函数,z(x+t)表示在z(x)时间间隔t后的值,S(ω)为功率谱密度函数;c为尺度系数,反映表面轮廓的粗糙程度;
对时间间隔t取不同值时,根据测量得到的z(x)计算得出对应的S(t),然后在双对数坐标中做出logS(t)-logt的关系曲线,再对这个曲线的最佳线性段进行拟合,拟合直线的斜率记为k,截距为b,则表面二维轮廓的分形维数D为
D=(4-k)/2
表面二维轮廓的特征粗糙度G为
G=exp[b/2(D-1)]
作为一种具体示例,步骤2中,使用递归分析法提取步骤1所得表面二维轮廓的递归率R r ,所用公式为:
其中,R r 是递归率,z i 和z j 都是轮廓高度,i,j=1,2,…,N,|·|为无穷大范数,r ij 表示任意两点之间的距离,R i,j 是矩阵的一个元素,ε为阈值,取0.5σ,σ为数据集合的标准差,θ(·)为Heaviside阶跃函数。
作为一种具体示例,步骤3中通过灰色关联度法计算步骤1所采集预紧拖动力与四个特征参数的关联度,具体如下:
关联度计算如下:
预紧拖动力退化率计算公式如下:
作为一种具体示例,步骤5中高斯回归过程模型如下:
对于给定训练集
简化回归模型为:
取均值函数为零,使用试验获得的真实数据进行训练后,得到输出值y的先验分布
将构建好的训练集导入高斯回归过程模型中进行训练,得到训练好的高斯回归过程模型。
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
实施例1
下面首先结合图2对滚柱直线导轨副的加载运行试验进行介绍,加载运行试验台主体由床身、龙门加载装置、齿轮2、齿条3、滚柱直线导轨副7、伺服电机1以及导轨转接板8组成,龙门加载装置由手轮6、垫块4、力传感器5组成。将两根滚柱直线导轨副7通过螺栓安装在加载运行试验台的转接板上;将设计好的压块利用螺栓固定在滑块上,使龙门加载装置中的压头对准压块的中心位置进行加载;以小速度试运行并注入润滑脂(glp-500)充分润滑;正式试验参数设置为:外加载荷33.48kN(即60%C,C表示被测滚柱直线导轨副的基本额定动载荷)、运行速度为50m/min;试验每进行5km测量一次滚道表面二维轮廓和预紧拖动力。
接下来介绍滚道表面二维轮廓和预紧拖动力测量试验,滚柱直线导轨副导轨滚道表面二维轮廓采集装置主要包括泰勒霍普森表面轮廓仪11、两个V型块9、一台数据采集设备14组成,其中V型块9、导轨滚道10、泰勒霍普森表面轮廓仪11、触针12、数据线13、数据采集设备14如图3所示;
滚道表面二维轮廓试验与显微图片拍摄前,都需要选取滚柱直线导轨副导轨滚道上的六个固定测量区域进行标记;将滚柱直线导轨副置于两个V形块时需要保持水平;表面轮廓仪每次使用前需要使用6μm的标准样片进行校准;在导轨的跑合段选取6个区域,定位方式为从没有油嘴的一端开始第9、13、18、21、26、30个孔所对应的导轨滚道区域,区域面积为5mm×10mm,每个区域选取三个固定位置,一共得到9组表面二维轮廓,对其进行滤波得到粗糙度轮廓。
接着进行滚柱直线导轨副的预紧拖动力测量试验,如图4所示,预紧拖动力测量试验台主要由床身17、移动平台18、齿轮齿条、测力计15、垫块16组成。将被测滚柱直线导轨副7利用螺栓和垫片固定在移动平台18上,导轨应平行于测力计15的移动方向和测量方向;测量行程与加载运行试验行程保持一致,测试速度应在10mm/s-20mm/s范围内,本次实验设定为11.6mm/s(0.7m/min);数据处理时,去掉测量长度前5%以及后5%的数据,正向拖动、反向拖动各测量3次,取平均值作为导轨副的预紧拖动力。
结合图1,本实施例提供了滚柱直线导轨副的预紧拖动力预测方法,包括以下步骤:
首先在滚柱直线导轨副的导轨滚道上选取六个测量区域,每个测量区域选取三个位置进行标记。对滚柱直线导轨副进行摩擦磨损试验,每运行5km停机,使用泰勒霍普森轮廓仪采集六个区域位置滚柱直线导轨副导轨滚道表面轮廓曲线,并对该轮廓进行去除形状以及高斯滤波处理,得到滚道表面二维轮廓的粗糙度曲线,同时通过预紧拖动力试验台测量其预紧拖动力,得到预紧拖动力变化曲线。
然后通过传统几何参数法对滚道表面二维轮廓进行特征提取,针对粗糙表面二维轮廓,最常用的几何特征参数就是轮廓算数偏差和轮廓最大峰谷距,两个特征参数可以通过下式定义:
其中,L为采样长度,z(x)为测量轮廓高度函数,z max 为表面轮廓高度的最大值,z min 为表面轮廓高度最小值;z i 为第i个采样点的轮廓高度,i=1,2,…N,N为采样点数。
通过单重分形法提取分形维数。用传统的几何参数对粗糙表面进行表征时,特征参数依赖于测量仪器的分辨率和采样长度,而分形理论具有尺度无关特性,因此引入分形理论对二维粗糙表面进行表征。分形维数的计算方法主要有盒维数法、功率谱法、概率密度法、结构函数法、均方根法等。
本发明使用结构函数法计算分形维数,W-M函数常被用来模拟和表征粗糙表面轮廓,表达式为:
通过任何测量技术测量表面形貌时,最低频率与采样长度有关,最高频率取决于测量仪器或滤波器的分辨率。
测量轮廓高度函数z(x)的功率谱函数S(ω)为:
定义z(x)的增量方差为结构函数,如下式所示:
其中,S(t)表示表面轮廓高度函数的结构函数,[·]2表示差方的算数平均值,t表示任意两数据之间的时间间隔;ω表示频率,即粗糙度波长的倒数;z(x)为测量轮廓高度函数,z(x+t)表示在z(x)时间间隔t后的值,S(ω)为功率谱密度函数;c为尺度系数,反映表面轮廓的粗糙程度;
对时间间隔t取不同值时,根据测量得到的z(x)计算得出对应的S(t),然后在双对数坐标中做出logS(t)-logt的关系曲线,再对这个曲线的最佳线性段进行拟合,拟合直线的斜率记为k,截距为b,则表面二维轮廓的分形维数D为
D=(4-k)/2
表面二维轮廓的特征粗糙度G为
G=exp[b/2(D-1)]
使用递归分析法提取所得表面二维轮廓的递归率R r ,递归分析法是非线性时间序列的一个重要方法,粗糙表面轮廓可以看作是一个时间序列z(x),因此递归分析可以用来描述粗糙表面内部的结构变化。
该方法所用公式为:
其中,R r 是递归率,z i 和z j 都是轮廓高度,i,j=1,2,…,N,|·|为无穷大范数,r ij 表示任意两点之间的距离,R i,j 是矩阵的一个元素,ε为阈值,取0.5σ,σ为数据集合的标准差,θ(·)为Heaviside阶跃函数。
接着通过灰色关联度法计算预紧拖动力退化率与四个特征参数的关联度。灰色关联度是对事物或***因素间描述其关系强弱、大小和次序的一种度量方法。若表征这些因素的时间序列曲线形状彼此相似,则关联程度大;反之,则关联程度小。对比关联度大小可以得到与预紧拖动力相关性最高的特征参数。
具体过程为:
进行关联度计算:
预紧拖动力退化率计算公式如下:
根据得到的最优特征参数,利用五次多项式拟合,提取拟合曲线的数据作为特征集,所对应的预紧拖动力数据为标签,构建训练集。
建立预紧拖动力与最优特征参数之间的高斯过程回归模型。高斯过程回归方法是基于贝叶斯框架对数据进行回归分析的非参数机器学习算法,在处理小样本、高维和线性不可分问题方面具有明显的优势。模型具体建立过程如下:
对于给定训练集
简化回归模型为:
为了方便计算,取均值函数为零,使用试验获得的真实数据进行训练后,得到输出值y的先验分布
基于高斯回归模型,将训练集数据导入模型进行计算,并使用5倍交叉验证防止模型过拟合。
将滚道上第六个区域的最优特征参数导入训练好的模型进行测试,实现滚柱直线导轨副预紧拖动力的预测,整体的流程如图5所示。
选取试验滚柱直线导轨副陕西汉江机床厂生产的DZ35型导轨,主要参数如表1所示。
表1滚柱直线导轨副参数表
在运行试验中途不断测量滚柱直线导轨副的滚道表面二维轮廓数据和预紧拖动力数据,直至滚柱之直线导轨副表面出现点蚀剥落,如图6所示,此时滚柱直线导轨副总运行里程为50km。将所提取的六个区域表面二维轮廓通过高斯滤波得到其粗糙度,随后通过前文提到的特征提取方法得到四个特征参数随运行里程的变化。
接着计算得到预紧拖动力的退化率,计算得到预紧拖动力退化率和轮廓算数偏差R a 、最大峰谷距R t 、分形维数D、递归率R r 的关联度为0.645、0.657、0.718,0.722如图7所示。因此递归率R r 为最优特征参数。
下一步通过五次多项式拟合预紧拖动力和递归率随运行里程的变化,提取五个区域对应的302个插值数据作为样本训练集。
接着通过MATLAB建立高斯回归过程模型将得到的样本训练集导入进行训练,选用5倍交叉验证以防止过拟合,回归模型评价指标方根误差RAME和决定系数R 2分别为0.00935和100%因此预紧拖动力回归模型具有较高的拟合度,模型训练结果如图8所示。
最后第六个区域的三组滚道表面二维轮廓递归率导入训练好的模型进行预测,三组预测结果的精度分别为93.75%,98.5%和98.8%,说明模型的预测精度很高,能够依据递归率准确地预测预紧拖动力,如图9所示。
Claims (5)
1.一种滚柱直线导轨副的预紧拖动力预测方法,其特征在于,包括以下步骤:
步骤1,通过泰勒霍普森轮廓仪,采集滚柱直线导轨副中导轨滚道M个区域对应的表面二维轮廓,并对表面二维轮廓进行去除形状以及高斯滤波处理,同时测量滚柱直线导轨副预紧拖动力变化;M为大于5的自然数;
步骤2,使用几何参数法、单重分形法、递归分析法提取步骤1所得表面二维轮廓的特征参数,特征参数包括轮廓算数偏差R a 、最大峰谷距R t 、分形维数D、递归率R r ;
步骤3,基于所得到的四个特征参数,通过灰色关联度法计算步骤1所采集预紧拖动力与四个特征参数的关联度,对比关联度大小得到与预紧拖动力相关性最高的特征参数作为最优特征参数;
步骤4,将得到的前M-1个区域的最优特征参数和预紧拖动力进行M-1次多项式拟合,提取拟合曲线中的数值点,最优特征参数作为样本特征,预紧拖动力数值作为标签,构建训练集;
步骤5,建立高斯回归过程模型,将训练集导入高斯回归过程模型进行训练;
步骤6,将第M个区域的最优特征参数导入训练好的高斯回归过程模型,预测得到预紧拖动力数值,并与真实值进行比较,验证高斯回归过程模型准确性;
步骤7,将验证后的高斯回归过程模型用于滚柱直线导轨副的预紧拖动力预测;
步骤2中,使用递归分析法提取步骤1所得表面二维轮廓的递归率R r ,所用公式为:
其中,R r 是递归率,z i 和z j 都是轮廓高度,i,j=1,2,…,N,|·|为无穷大范数,r ij 表示任意两点之间的距离,R i,j 是矩阵的一个元素,ε为阈值,取0.5σ,σ为数据集合的标准差,θ(·)为Heaviside阶跃函数;
步骤3中通过灰色关联度法计算步骤1所采集预紧拖动力与四个特征参数的关联度,具体如下:
关联度计算如下:
2.根据权利要求1所述的滚柱直线导轨副的预紧拖动力预测方法,其特征在于,步骤1中M=6,首先在滚柱直线导轨副的导轨滚道上选取6个测量区域,每个测量区域选取3个位置进行标记;对滚柱直线导轨副进行摩擦磨损试验,每运转5km停机,使用泰勒霍普森轮廓仪采集滚柱直线导轨副中导轨滚道6个区域对应的18个表面二维轮廓,并对表面二维轮廓进行去除形状以及高斯滤波处理,同时利用预紧拖动力试验台测量滚柱直线导轨副预紧拖动力,得到预紧拖动力变化曲线。
4.根据权利要求3所述的滚柱直线导轨副的预紧拖动力预测方法,其特征在于,步骤2中,使用单重分形法提取步骤1所得表面二维轮廓的分形维数D,具体如下:
W-M函数被用来模拟和表征表面二维轮廓,表达式为:
测量轮廓高度函数z(x)的功率谱函数S(ω)为:
定义z(x)的增量方差为结构函数,如下式所示:
其中,S(t)表示表面轮廓高度函数的结构函数,[·]2表示差方的算数平均值,t表示任意两数据之间的时间间隔;ω表示频率,即粗糙度波长的倒数;z(x)为测量轮廓高度函数,z(x+t)表示在z(x)时间间隔t后的值,S(ω)为功率谱密度函数;c为尺度系数,反映表面轮廓的粗糙程度;
对时间间隔t取不同值时,根据测量得到的z(x)计算得出对应的S(t),然后在双对数坐标中做出logS(t)-logt的关系曲线,再对这个曲线的最佳线性段进行拟合,拟合直线的斜率记为k,截距为b,则表面二维轮廓的分形维数D为
D=(4-k)/2
表面二维轮廓的特征粗糙度G为
G=exp[b/2(D-1)]。
5.根据权利要求4所述的滚柱直线导轨副的预紧拖动力预测方法,其特征在于,步骤5中高斯回归过程模型如下:
对于给定训练集
简化回归模型为:
取均值函数为零,使用试验获得的真实数据进行训练后,得到输出值y的先验分布
将构建好的训练集导入高斯回归过程模型中进行训练,得到训练好的高斯回归过程模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211587676.6A CN115660224B (zh) | 2022-12-12 | 2022-12-12 | 一种滚柱直线导轨副的预紧拖动力预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211587676.6A CN115660224B (zh) | 2022-12-12 | 2022-12-12 | 一种滚柱直线导轨副的预紧拖动力预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115660224A CN115660224A (zh) | 2023-01-31 |
CN115660224B true CN115660224B (zh) | 2023-03-17 |
Family
ID=85019995
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211587676.6A Active CN115660224B (zh) | 2022-12-12 | 2022-12-12 | 一种滚柱直线导轨副的预紧拖动力预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115660224B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110095280A (zh) * | 2019-04-29 | 2019-08-06 | 南京理工大学 | 一种滚动直线导轨副综合磨损系数测试方法 |
CN113011041A (zh) * | 2021-04-01 | 2021-06-22 | 南京理工大学 | 考虑微观接触特性的滚动直线导轨副预紧力衰退计算方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115272694A (zh) * | 2022-07-06 | 2022-11-01 | 南京理工大学 | 基于滚道表面轮廓的滚珠丝杠副剩余使用寿命预测方法 |
CN115205512A (zh) * | 2022-08-12 | 2022-10-18 | 南京理工大学 | 基于滚道表面轮廓的滚珠丝杠副磨损状态识别方法 |
-
2022
- 2022-12-12 CN CN202211587676.6A patent/CN115660224B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110095280A (zh) * | 2019-04-29 | 2019-08-06 | 南京理工大学 | 一种滚动直线导轨副综合磨损系数测试方法 |
CN113011041A (zh) * | 2021-04-01 | 2021-06-22 | 南京理工大学 | 考虑微观接触特性的滚动直线导轨副预紧力衰退计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115660224A (zh) | 2023-01-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107562696B (zh) | 轮胎产品质量在线检测与控制方法 | |
CN101246369A (zh) | 一种汽车零部件尺寸质量控制***和方法 | |
CN113468479B (zh) | 一种基于数据驱动的冷连轧工业过程监测与异常检测方法 | |
CN116070527B (zh) | 基于退化模型的铣削刀具剩余寿命预测方法 | |
CN113486457A (zh) | 一种压铸件缺陷预测及诊断*** | |
CN115272694A (zh) | 基于滚道表面轮廓的滚珠丝杠副剩余使用寿命预测方法 | |
CN115741235A (zh) | 基于五轴加工中心刀具的磨损预测与健康管理方法 | |
CN116523129A (zh) | 基于xgboost与分形维数的路面抗滑性能预测方法 | |
CN115660224B (zh) | 一种滚柱直线导轨副的预紧拖动力预测方法 | |
CN111639715B (zh) | 基于ls-svm的汽车仪表组装质量预测方法与*** | |
CN113158558B (zh) | 一种高速铁路路基连续压实分析方法、装置及分析仪 | |
CN117781882A (zh) | 一种工业五金件检测mes数据采集分析*** | |
CN109648397B (zh) | 基于刀齿刃带宽度与拉削负载相关性的拉刀寿命预测方法 | |
CN109614758B (zh) | 具有空间相关性的圆形形状误差的监控方法 | |
CN113780852B (zh) | 一种板带轧制过程质量缺陷的诊断方法 | |
CN113973403B (zh) | 基于结构离散测点拓扑的温致应变场重分布智能感知方法 | |
CN101629407A (zh) | 一种路面结构强度预测方法 | |
CN115985411A (zh) | 基于高斯过程回归模型的高分子熔融指数软测量方法 | |
CN110967184B (zh) | 基于振动信号分布特征识别的变速箱故障检测方法和*** | |
CN111368376B (zh) | 基于高斯过程模型的圆形轮廓误差监控方法 | |
CN116362447A (zh) | 一种基于层次分析法的钻井施工难度评价方法 | |
CN112069621A (zh) | 基于线性可靠度指标的滚动轴承剩余使用寿命的预测方法 | |
Liu et al. | Research on intelligent measurement system of workpiece outer dimension | |
CN115468679B (zh) | 一种基于大数据和数字化工厂的产品检测方法 | |
CN117557574B (zh) | 基于图像处理的材料参数检测方法及*** |
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 |