CN101226068A - 动态摩擦参数测试***及测试方法 - Google Patents
动态摩擦参数测试***及测试方法 Download PDFInfo
- Publication number
- CN101226068A CN101226068A CNA2008100174737A CN200810017473A CN101226068A CN 101226068 A CN101226068 A CN 101226068A CN A2008100174737 A CNA2008100174737 A CN A2008100174737A CN 200810017473 A CN200810017473 A CN 200810017473A CN 101226068 A CN101226068 A CN 101226068A
- Authority
- CN
- China
- Prior art keywords
- friction
- moment
- parameter
- test
- motor
- 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
Images
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种动态摩擦参数的测试***及测试方法,该***在电机(13)上安装绝对编码器(14),在测试件及砝码固定装置(6)与锁紧机构(2)之间同轴安装力矩传感器(4),测试件安装在砝码固定装置的上下层。测试时通过控制器控制电机运行,同步记录采样时刻、摩擦力矩、转速、转角及电机绕组电流参数;设定当前摩擦参数的估计值0、1、c、s、s、2,并求解***运动微分方程,得到每个采样时刻的力矩电机转角估计值(t);分别计算各时刻的力矩电机转角估计值误差和摩擦力矩估计值误差,将该误差参数与所设定的测试参数比较,满足精度要求的输出摩擦参数估计值,否则重新修正。本发明具有测试简便和精度高的优点,适用对各种材料的动态摩擦参数测试。
Description
技术领域
本发明属于测试技术领域,特别是涉及动态摩擦参数的测量,用于指导高性能伺服***的设计。
背景技术
在伺服***中,具有相对运动或相对运动趋势的两个接触面上会产生摩擦力。摩擦降低了***的稳定性和稳态精度,导致了低速爬行现象,对***性能的提高构成了严重的障碍。大量的研究表明摩擦是一种复杂的物理现象,摩擦系数不仅与接触面间的相对运动速度非线性相关,称之为摩擦的静态特性,而且还存在摩擦记忆、增加的静摩擦力等动态效应。Armstrong B,Dupont P,Canudas de Wit C.在《Automatica》,1994年30卷第7期上发表的论文“A Survey of Models,Analysis Tools and Compensation Methods forthe Control of Machines with Friction”对此进行了很好的综述。在伺服***设计时,不考虑或仅考虑最简单的摩擦现象,如静摩擦力或Coulomb摩擦力,将导致跟踪精度下降、控制参数整定困难、低速运行平稳性差等问题。为了更好地描述伺服***中存在的摩擦并对其造成的危害进行补偿,有必要研究精确的摩擦模型及其参数的测试与辨识,即首先对摩擦进行准确的建模,其次精确测量出特定情况下接触表面间的摩擦力,然后依据特定的摩擦模型,采用参数辨识的方法获得模型中各参数的数值。到目前为止,人们已提出了多达30余种摩擦模型,其中最具影响力的模型是Canudas De Wit C等发表于《IEEE Transaction on Automatic Control》1995年40卷第3期的论文“A New Model forControl of Systems with Friction”中提出的LuGre模型,它能全面地描述接触面从静止到运动各个状态摩擦力的变化情况。下式即为LuGre模型的数学表达式,其中T为接触表面正压力一定时的摩擦力或摩擦力矩,z为接触表面刚毛的平均变形量,ωl为相对运动角速度;刚度系数σ0、粘性阻尼系数σ1,Coulomb摩擦力矩FC、静摩擦力矩FS、临界Stribeck角速度ωS、粘性摩擦系数σ2为6个摩擦参数。
为了在工程中应用此模型,就需要测量出实际接触表面间的摩擦,然而,刚毛的平均变形量z是无法测量的,6个模型参数FC、FS、ωs、σ0、σ1、σ2也是高度耦合的,因此不能通过简单的实验获得相应参数的取值,只有采用参数辨识的方法才能得到这6个参数的最佳估计。
目前国内外已有的摩擦参数测试设备有:
(1)英国RAY-RAN公司的FT摩擦测试***
该测试***能自动测量接触表面间的静摩擦系数和动摩擦系数,并能进行数据存储、显示,该静摩擦系数相当于LuGre模型的Fs,该动摩擦系数相当于LuGre模型中的Fc,σ2。由于使用的模型简单,测量数据中既无法体现Stribeck效应,又不能描述摩擦记忆、增加的静摩擦力这些动态摩擦效应。
(2)济南实验机厂的MHK-500摩擦磨损试验机
试验机采用手动加载、人工调节砝码和游码的方法测量接触面间的摩擦系数,即LuGre模型中的Fs和Fc,无法测量速度与摩擦力之间的关系,更不能测量动态摩擦参数,如LuGre模型中的σ0、σ1。目前已有在其基础上通过增加电子装置实现自动测量的研究报道,如俞建卫发表于《密封与润滑》2007,32(2)的论文:“MHK-500环块摩擦磨损试验机智能测控***的研制”,但仍不能测试动态摩擦参数。
目前国内外在LuGre摩擦参数的测量方法主要有二步辨识法,该方法最先由C.Candus等发表于《Internatinal Journal of Adaptive Control and Signal Processing》1997年第11卷的论文“Adaptive friction compensation with partially known dynamic frictionmodel”首先提出,其思路是将LuGre模型中的参数分为两类,即4个静态参数FC、FS、ωS、σ2和2个动态参数σ0、σ1。假定在匀速运动时,刚毛的平均变形量z不发生改变,即 这样上述LuGre模型就可简化为
即将动态的LuGre模型变为静态的Stribeck模型。于是辨识步骤为:
第一步:采用匀速控制方法,分别测量相对运动速度或角速度和对应的摩擦力或摩擦力矩,就可使用曲线拟合的方法确定4个静态参数FC、FS、ωS、σ2。
第二步:在确定静态参数的前提下,辨识动态参数σ0、σ1。
对于动态参数,目前又有以下几种辨识方法:
(1)线性化方法:给***作用一个较小的力u,使其出现预滑动微小转角θ,直接利用测量出的转角,由式 辨识出刚度系数σ0。
将滞滑运动时的摩擦模型在平衡点处线性化,可得:
其中,J为运动部分转动惯量。
此方法将***简化为一个二阶***,使其处于临界阻尼状态,利用表达式 得到σ1。例如C.Candus等发表于《International Journall of AdaptiveControl and Signal Processing》1997年第11卷的论文“Adaptive friction compensationwith partially known dynamic friction model就使用这种方法。
由于这种线性化方法的辨识结果依赖于参数初值的选取,因此辨识精度及收敛性较难保证,而且预滑动微小位移的测量非常困难。
(2)优化方法:驱动负载使其进行低速周期运动,将位置误差作为优化目标、两个动态参数σ0、σ1作为设计变量进行优化计算,最终获得满意的σ0、σ1取值。如DE-PENG LIU 2005年发表于《Proceedings of the Fourth International Conference onMachine Learning and Cybernetics》的论文“Research on Parameter Identification of FrictionModel for Servo System Based on Genetic Algorithms”采用的就是这种方法。
(3)频域辨识法:采用局部线性化的方法将原***简化为二阶线性***,采用频域方法辨识动态参数。使***处于预滑动阶段,给***输入任意噪声激励,就能测出频域响应函数。频率ω接近0的区域内,根据|G(j·0)|=Cm/σ0,直接辨识出刚度系数σ0,依据测量出的频率响应函数,估计出二阶线性***的系数,可辨识出σ1+σ2。例如RonH.A.Hensen等人发表于《IEEE Transactions on Control systems Technology》2002年第10卷第2期的论文“Frequency Domain Identification of Dynamic Friction ModelParameters”描述的就是这种方法。
上述辨识方法均存在如下不足:
1.由于这些方法都要采用匀速实验获得静态参数,而匀速实验要逐点进行,因此工作量大、耗时,同时还要求***在低速运行时不产生爬行现象,不易实现。事实上在伺服***中测试摩擦的目的之一就是获得良好的低速运行平稳性,在未知摩擦参数时就要达到这一目的,难度很大。
2.由于线性化方法和频域方法都对原模型进行了近似处理,因而对于动态参数的测量必然导致其准确性下降。
3.由于原LuGre模型中的动态参数和静态参数本身就是耦合在一起的,因此按“匀速运动对刚毛的平均变形量z不发生改变”的假设,将***解耦也没有理论依据的。
综上所述,目前既没有一种能简捷、有效地辨识LuGre摩擦参数的方法,也没有可用于获取相关数据的仪器设备。
发明的内容
本发明的目的在于避免上述现有技术的不足,提供一种测量动态摩擦参数的***及方法,以便捷地实现对不同材料间动态摩擦参数的高精度测量,提高伺服***的性能。
实现本发明目的的技术思路是:采用高精度力矩传感器和高精度绝对编码器设计摩擦力矩测试装置,采用停滞时间可调的变周期正弦驱动方式,使接触面间摩擦的动、静态特性得以充分体现,并辅以控制器实现摩擦力矩、转速、转角及驱动力矩的实时测量、存储;将测试数据进行自动分析处理,从而一次性得到全部的6个参数的辨识值。
本发明的测试***包括:试件及砝码固定装置、电机、驱动轴、驱动轴、锁紧机构和控制器,其中,电机下面固定有绝对编码器,用于检测试件表面的相对位置与速度;试件及砝码固定装置与锁紧机构之间通过固定轴同轴安装有力矩传感器,用于检测摩擦力矩。
所述的试件及砝码固定装置分上下两部分,下部与驱动轴固定为一体,上部与力矩传感器固连;该下部的表面设有正方形的凹槽,用于安装下测试件,上部安装有上测试件。
所述的电机采用力矩电机,并与绝对编码器同轴安装在驱动轴的底部,直接驱动负载运转。
所述的控制器中设有自动测试和参数自动辨识程序。该控制器分别与绝对编码器、力矩传感器和力矩电机连接,构成闭环控制回路,以控制测试件表面的相对速度、测量实时位置、速度和摩擦力矩。
所述的驱动轴的中部设有二层水平支架,该支架与机架固定为一体。
本发明的测试方法采用一步辨识法,其技术原理是:
当所估计的各参数接逼近真值时,摩擦力矩的误差和***位置响应的误差最小,这样就可将参数辨识问题转化为求解下述优化问题:
find σ0、σ1,FC、FS、ωS、σ2
其中J为***的转动惯量,Tu为驱动力矩,θi、和Tfi分别为第i个采样点的转角预估值、实测转角、摩擦力矩的预估值和摩擦力矩的实测值,λ1为表征累计转角误差要求的权值,其取值范围为0~1,λ2为表征最大转角误差要求的权值,其取值范围为0~1,λ3为表征累计摩擦力矩误差要求的权值,其取值范围为0~1,λ4为为表征最大摩擦力矩误差要求的权值,其取值范围为0~1。
根据上述原理本发明的测试方法,包括如下过程:
(1)安装测试件并设置测试参数;
(2)启动控制器按照设置的测试参数控制力矩电机驱动负载运行,并同步记录采样时刻、摩擦力矩、转速、转角及电机绕组电流参数;
(3)设定当前摩擦参数,即刚度系数σ0、粘性阻尼系数σ1,Coulomb摩擦力矩FC、静摩擦力矩FS、临界Stribeck角速度ωS、粘性摩擦系数σ2的相应估计值并根据所述电机绕组电流通过求解***运动微分方程,得到每个时刻的力矩电机转角估计值
(4)计算每个时刻的力矩电机转角估计值与记录该时刻的转角实测值的误差累计值EA,并找出最大的力矩电机转角误差值EMA;
(5)计算每个时刻的摩擦力矩估计值与记录该时刻的摩擦力矩实测值的误差累计值ET,并找出最大的摩擦力矩误差值EMT;
本发明与现有技术相比具有如下优点:
1).由于本发明的测试***设有控制器,可以实现被测材料的自动测试和参数的自动辨识;同时由于本发明***设有绝对编码器,可有效检测被测材料表面的相对转角与角速度;此外由于本发明***安装有力矩传感器,可实现对摩擦力矩进行高精度测量。
2)由于本发明的测试方法采用一步辨识法,无需进行匀速实验,因而测试简单,使所有摩擦参数一次辨识完成;同时解决了按“匀速运动不改变刚毛的平均变形量z”假设没有理论依据的问题。
仿真结果表明,采用本发明动态摩擦测试方法,6个摩擦参数的平均误差为0.601%。
附图说明
图1是本发明测试***的结构示意图;
图2是本发明控制器原理框图;
图3是本发明一步辨别识方法测试流程图;
图4是本发明一步辨识方法的仿真实验结果图,其中(a)为迭代1000次各摩擦参数的辨识值表示图,(b)为迭代50000次各摩擦参数的辨识值表示图。
以下结合附图对本发明作进一步详细描述。
具体实施方式
参照图1,本发明的测试***包括:机械结构部分由手柄1、锁紧机构2、力矩传感器4、测试件及砝码固定装置6、待测试件7、驱动轴8、机架10、水平支架11、力矩电机13、绝对编码器14和控制器16组成。该测试件及砝码固定装置6分为上下两层,下层6a的表面设有正方形凹槽,上层6b为阶梯轴,力矩传感器上固连有上下法兰盘4a和4b,机架10为整体铸造的长方体框架结构,便于操作,中间设有二层水平支架11以增强刚度,驱动轴8通过轴承9安装于水平支架11上,可以绕其轴线旋转。绝对编码器14与力矩电机13固定,并通过连轴器12同轴安装于驱动轴8的下端,该驱动轴的上端与测试件及砝码固定装置的下层6a固定,测试件及砝码固定装置的上层6b通过下法兰盘4a与力矩传感器4连接,该上层6b的两边放置砝码5,为上下测试件施加载荷,力矩传感器4通过上法兰盘与固定轴3下端连接,固定轴3通过锁紧机构2与机架10连接,固定轴3的上端与手柄1固连。测试件7分别置于测试件及砝码固定装置的上下层,其中下测试件7a置于测试件及砝码固定装置的下层6a的凹槽内,上测试件7通过螺钉15固定在测试件及砝码固定装置的上层6b表面。控制器16通过电缆分别与力矩电机13、绝对编码器14和力矩传感器4连接,构成闭环控制回路,以控制测试件表面的相对速度、测量实时位置、速度和摩擦力矩。
参照图2,控制器16以数字信号处理器DSP为核心,由用户界面、电机控制部分和摩擦力矩检测部分组成,具有运动控制、数据采集、信息处理等功能。
(1)DSP
所述的DSP用于接收用户设置的测试信息、对测试件进行自动测试和参数自动辨识、将测试结果发送到PC机,该DSP中固化有自动测试和参数自动辨识程序,并扩展有64KB的静态RAM,以存储测试数据。本实例的DSP采用TMS320LF2407A型号,但不限于此型号。
(2)用户界面
所述的用户界面主要由320×240点阵LCD显示屏、键盘和USB2.0接口电路组成。
320×240点阵LCD显示屏通用并行I/O端口与DSP通过相连,DSP可将控制命令和字符信息通过并行I/O口发送给LCD显示屏,显示出设定的信息、工作状态和参数辨识结果。
键盘采用16键矩阵式键盘,但不限于16键,通过8根数据线与DSP的通用I/O口相连,DSP通过扫描相应I/O口的状态获取按键信息。通过键盘输入的参数包括:运行速度范围、运行周期、辨识精度、停滞时间范围。
USB2.0接口由接口控制器PHILIPS ISP1581构成,通过***总线与DSP相连,同时USB2.0接口控制器与PC之间采用现有的USB电缆连接。测量辨识完成后,可将测量及辨识结果发送到PC机上进行存储。
(3)电机控制部分
所述的电机控制部分由D/A转换器、伺服放大器、电机绕组电流检测电路、绝对式编码器组成。
D/A转换器采用12位高速D/A芯片MX7541,转换速度0.6μs,其数字输入端通过***总线与DSP相连,其模拟输出端与伺服放大器的模拟输入端相连。
伺服放大器的模拟输入端与D/A转换器的输出端相连,输出端与力矩电机连接。
绕组电流检测电路由精密电阻及A/D转换器组成,精密电阻串联于电机驱动回路中,其两端的压降反映了电机绕组中电流的大小,采用高速12位A/D转换器MAX120,转换速度1.6μs,将此电压信号转换为数字量并通过***总线传送给DSP。
绝对式编码器通过RS485接口与DSP相连,该绝对式编码器用于检测力矩电机的转角。
该电机控制部分可构成电流环、速度环和位置环三闭环位置控制回路,也可工作于开环方式。在进行摩擦参数测试时,使用开环方式工作,DSP根据输入的运行速度范围给出呈正弦变化的控制信号,经D/A变换后送伺服放大器产生相应的驱动电压,驱动电机带动负载跟随参考信号运动。
(4)摩擦力矩检测部分
所述的摩擦力矩检测部分由高精度的静态力矩传感器Staiger Mohilo 0510RD,量程10Nm,精度±0.02Nm,和12位高速A/D转换器MAX120构成。
静态力矩传感器输出的模拟信号经A/D变换后存储于扩展的RAM中。A/D转换器MAX120和静态RAM均通过***总线与DSP相连。
测试时,力矩电机13通过传动轴8驱动下测试件7a运转,上测试件7b保持静止,通过力矩传感器4和控制器就可测量出上下测试件7b和7a间的摩擦力矩。
参照图3,本发明对动态摩擦参数的测量主要包括测试准备、自动测试和参数自动辨识三个过程,具体步骤如下:
第一步:测试准备,包括安装测试件和设置测试参数。
1)安装测试件:首先松开锁紧机构,向上提起手柄,使摩擦测试面脱离,再锁紧手柄,使其不至下落;其次,将测试件分别安装在上下两个测量面上;最后,松开锁紧机构2,使上下试件表面完全接触;
2)设定接触面压力:根据不同的工况,选择适当的砝码加载,锁紧手柄;
3)参数设置:通过键盘输入运行速度范围、参数辨识精度、测试周期数、停滞时间范围参数。
第二步:自动测试
按下键盘上的“测试”按键,力矩电机在DSP控制电路的作用下启动,其转速在要求的运行范围内按正弦曲线变化,同时将每1ms采样周期的摩擦力矩、转速、转角、电机绕组电流这些参数记录在静态RAM中,在完成设定的测试周期后,测试结束。
第三步:参数自动辨识
1)设定当前摩擦参数
2)计算力矩电机转角的估计值
该转角的估计值可通过四阶龙格库塔法对以下运动微分方程进行求解获得
其中,为t时刻摩擦力矩的估计值,依据当前各参数的估计值 按LuGre模型的数学表达式(1)可计算得到;J为已知的***转动惯量,含电机转子、传动轴及试件、砝码固定装置的下部及下测试件;Ka为已知的电机力矩常数;I(t)第为t时刻电机绕组的电流的测量值,I(t)从记录的第i个采样点的实验数据中导出,设tp为采样周期,则t=i·tp;为t时刻的转角。
3)计算累计转角估值误差并选出最大转角估值误差
累计转角误差根据下式计算:
最大转角估值误差根据下式选择:
4)计算累计摩擦力矩估值误差并选出最大摩擦力矩估值误差
最大摩擦力矩估值误差根据下式选择:
5)计算目标函数,判断其是否达到给定的估值精度目标函数按下式计算:
J=λ1EA+λ2EMA+λ3ET+λ4EMT
其中,λ1为表征累计转角误差要求的权值,其取值范围为0~1;λ2为表征最大转角误差要求的权值,其取值范围为0~1;λ3为表征累计摩擦力矩误差要求的权值,其取值范围为0~1;λ4为表征最大摩擦力矩误差要求的权值,其取值范围为0~1。
如果摩擦参数估值达到给定的精度,则在LCD显示屏上输出结果,并将实验测试数据和辨识结果通过USB接口发送到PC机上进行存储。
本发明的效果可通过以下仿真数据进一步说明:
仿真条件:取J=1Kgm2,给定一组摩擦参数作为参数的实际值,以最优个体作为辨识输出,采用遗传算法计算100代、1000代和50000代的辨识值。
仿真结果:如图4和表1所示,由图4(a)可看出遗传算法运行1000代时,各参数的最优个体取值即辨识值已趋近于其真值,由图4(b)可看出遗传算法运行50000代时,各参数的辨识值已满足误差要求。由表1可知遗传算法运行1000代时的6个参数的平均辨识误差为2.612%,最大辨识误差为9.4%;遗传算法运行50000代时,参数σ0、σ1,FS、FC、ωS、σ2的辨识误差分别为0.000%、0.030%、0.596%、0.007%、0.030%、2.940%,这6个参数的最大辨识误差为2.940%,平均辨识误差为0.601%,从而说明了本发明方法的有效性。
表1
参数 | 实际值 | 辨识值(100代) | 辨识值(1000代) | 辨识值(50000代) |
σ0 | 0.5 | 0.48721 | 0.50000 | 0.50000 |
σ1 | 0.4 | 0.4073 | 0.39980 | 0.40012 |
FS | 2.4 | 2.55469 | 2.62620 | 2.38570 |
FC | 0.6 | 0.73354 | 0.59981 | 0.60004 |
ωS | 0.01 | 0.019 | 0.009773 | 0.010003 |
σ2 | 0.005 | 0.00539 | 0.0051964 | 0.004853 |
Claims (10)
1.一种动态摩擦参数测试***,包括测试件及砝码固定装置、电机、驱动轴、锁紧机构和控制器,其特征在于,电机(13)下面固定有绝对编码器(14),用于测量测试件表面的相对转角与角速度;测试件及砝码固定装置(6)与锁紧机构(2)之间通过固定轴(3)同轴安装有力矩传感器(4),用于检测摩擦力矩。
2.根据权利要求1所述的测试***,其特征在于,测试件及砝码固定装置(6)分上下两层(6b和6a),下层(6a)与驱动轴(8)固定为一体,上层(6b)与力矩传感器(4)固连;该下层的表面设有正方形的凹槽,用于安装下测试件(7a),上层安装有上测试件(7b)。
3.根据权利要求1所述的测试***,其特征在于,电机(13)采用力矩电机,并与绝对编码器(14)同轴安装在驱动轴(8)的底部,直接驱动负载运转。
4.根据权利要求1所述的测试***,其特征在于,控制器(16)中设有自动测试和参数自动辨识程序。
5.根据权利要求1或4所述的测试***,其特征在于,控制器(16)分别与绝对编码器(14)、力矩传感器(4)和力矩电机(13)连接,构成闭环控制回路,以控制测试件表面的相对转速、测量实时转角、转速和摩擦力矩。
6.根据权利要求1所述的测试***,其特征在于,驱动轴(8)的中部设有二层水平支架(11),该支架与机架(10)固定为一体。
7.一种动态摩擦参数测试方法,包括如下过程:
(1)安装测试件并设置测试参数;
(2)启动控制器按照设置的测试参数控制力矩电机驱动负载运行,并同步记录采样时刻、摩擦力矩、转速、转角及电机绕组电流参数;
(3)设定当前摩擦参数,即刚度系数σ0、粘性阻尼系数σ1,Coulomb摩擦力矩FC、静摩擦力矩FS、临界Stribeck角速度ωS、粘性摩擦系数σ2的相应估计值并根据所述电机绕组电流通过求解***运动微分方程,得到每个时刻的力矩电机转角估计值
(4)计算每个时刻的力矩电机转角估计值与记录该时刻的转角实测值的误差累计值EA,并找出最大的力矩电机转角误差值EMA;
(5)计算每个时刻的摩擦力矩估计值与记录该时刻的摩擦力矩实测值的误差累计值ET,并找出最大的摩擦力矩误差值EMT;
8.根据权利要求7所述的测试方法,其特征在于步骤(1)所述的设置测试参数,包括力矩电机的运行速度范围、运行周期、停滞时间范围及摩擦参数辨识精度。
10.根据权利要求7所述的测试方法,其特征在于,步骤(6)所述的判断步骤(4)、(5)所得到的参数是否满足步骤(1)所设定的测试参数,按如下过程进行:
(6a).通过如下公式对步骤(4)、(5)所得到的参数进行求和:
J=λ1EA+λ2EMA+λ3ET+λ4EMT
式中,λ1为表征累计转角误差要求的权值,其取值范围为0~1,λ2为表征最大转角误差要求的权值,其取值范围为0~1,λ3为表征累计摩擦力矩误差要求的权值,其取值范围为0~1,λ4为为表征最大摩擦力矩误差要求的权值,其取值范围为0~1;
(6b).用所求和的值J与设定的辨识精度参数进行比较,小于该设定辨识精度参数的作为辨识值输出,反之重新进行修正。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100174737A CN101226068B (zh) | 2008-02-01 | 2008-02-01 | 动态摩擦参数测试***及测试方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100174737A CN101226068B (zh) | 2008-02-01 | 2008-02-01 | 动态摩擦参数测试***及测试方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101226068A true CN101226068A (zh) | 2008-07-23 |
CN101226068B CN101226068B (zh) | 2010-04-07 |
Family
ID=39858190
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008100174737A Expired - Fee Related CN101226068B (zh) | 2008-02-01 | 2008-02-01 | 动态摩擦参数测试***及测试方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101226068B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102269638A (zh) * | 2011-04-27 | 2011-12-07 | 中国科学院光电技术研究所 | 伺服转台LuGre模型摩擦参数及转动惯量的一体化测量方法 |
CN103411715A (zh) * | 2013-07-11 | 2013-11-27 | 河南科技大学 | 一种力矩电机阻力矩测试*** |
CN103975283A (zh) * | 2011-12-08 | 2014-08-06 | 西门子公司 | 确定驱动***的摩擦分量 |
CN105303017A (zh) * | 2014-06-26 | 2016-02-03 | 上银科技股份有限公司 | 线性移动平台的摩擦力参数鉴别方法 |
CN106596087A (zh) * | 2016-12-12 | 2017-04-26 | 山东大学 | 一种新型可实现全组件摩擦辨识的装置和方法 |
CN108195497A (zh) * | 2017-11-29 | 2018-06-22 | 北京控制工程研究所 | 一种动压气浮轴承在轨摩擦力矩测试方法 |
CN110118723A (zh) * | 2019-06-06 | 2019-08-13 | 辽宁工程技术大学 | 一种测试岩石自然断面摩擦系数装置及测试方法 |
CN110208691A (zh) * | 2018-02-28 | 2019-09-06 | 株式会社安川电机 | 一种电机负荷参数检测方法及设备 |
CN110672121A (zh) * | 2019-09-19 | 2020-01-10 | 北京控制工程研究所 | 一种控制力矩陀螺框架动态响应测试方法及*** |
CN111191344A (zh) * | 2019-12-09 | 2020-05-22 | 北京动力机械研究所 | 一种运动副摩擦参数的识别方法 |
CN112721895A (zh) * | 2021-01-26 | 2021-04-30 | 同济大学 | 基于新型摩擦模型的iehb***主缸液压力估计方法 |
CN113029194A (zh) * | 2021-02-26 | 2021-06-25 | 北京控制工程研究所 | 一种cmg组合体测试方法、性能评价方法及测试*** |
CN113110305A (zh) * | 2021-04-30 | 2021-07-13 | 华中科技大学 | 一种机电***的摩擦力建模方法及其应用 |
CN113359478A (zh) * | 2021-07-15 | 2021-09-07 | 广东工业大学 | 单自由度直线运动平台非均匀导轨摩擦参数的辨识方法 |
CN113420455A (zh) * | 2021-07-06 | 2021-09-21 | 同济大学 | 一种基于ehb自身特性的主缸液压力估计方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2604674Y (zh) * | 2003-03-27 | 2004-02-25 | 洛阳轴研科技股份有限公司 | 轴承摩擦力矩测量仪 |
CN2826406Y (zh) * | 2005-10-22 | 2006-10-11 | 中国矿业大学 | 不离位表面测量摩擦磨损试验机 |
CN201004011Y (zh) * | 2007-02-15 | 2008-01-09 | 洛阳轴研科技股份有限公司 | 一种轴系摩擦力矩测量仪 |
-
2008
- 2008-02-01 CN CN2008100174737A patent/CN101226068B/zh not_active Expired - Fee Related
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102269638A (zh) * | 2011-04-27 | 2011-12-07 | 中国科学院光电技术研究所 | 伺服转台LuGre模型摩擦参数及转动惯量的一体化测量方法 |
CN103975283A (zh) * | 2011-12-08 | 2014-08-06 | 西门子公司 | 确定驱动***的摩擦分量 |
CN103975283B (zh) * | 2011-12-08 | 2016-08-24 | 普锐特冶金技术德国有限公司 | 确定驱动***的摩擦分量 |
CN103411715A (zh) * | 2013-07-11 | 2013-11-27 | 河南科技大学 | 一种力矩电机阻力矩测试*** |
CN103411715B (zh) * | 2013-07-11 | 2016-04-06 | 河南科技大学 | 一种力矩电机阻力矩测试*** |
CN105303017B (zh) * | 2014-06-26 | 2018-04-20 | 上银科技股份有限公司 | 线性移动平台的摩擦力参数鉴别方法 |
CN105303017A (zh) * | 2014-06-26 | 2016-02-03 | 上银科技股份有限公司 | 线性移动平台的摩擦力参数鉴别方法 |
CN106596087A (zh) * | 2016-12-12 | 2017-04-26 | 山东大学 | 一种新型可实现全组件摩擦辨识的装置和方法 |
CN106596087B (zh) * | 2016-12-12 | 2019-07-05 | 山东大学 | 一种新型可实现全组件摩擦辨识的装置和方法 |
CN108195497A (zh) * | 2017-11-29 | 2018-06-22 | 北京控制工程研究所 | 一种动压气浮轴承在轨摩擦力矩测试方法 |
CN110208691A (zh) * | 2018-02-28 | 2019-09-06 | 株式会社安川电机 | 一种电机负荷参数检测方法及设备 |
CN110208691B (zh) * | 2018-02-28 | 2021-09-17 | 株式会社安川电机 | 一种电机负荷参数检测方法及设备 |
CN110118723A (zh) * | 2019-06-06 | 2019-08-13 | 辽宁工程技术大学 | 一种测试岩石自然断面摩擦系数装置及测试方法 |
CN110118723B (zh) * | 2019-06-06 | 2021-09-14 | 辽宁工程技术大学 | 一种测试岩石自然断面摩擦系数装置及测试方法 |
CN110672121A (zh) * | 2019-09-19 | 2020-01-10 | 北京控制工程研究所 | 一种控制力矩陀螺框架动态响应测试方法及*** |
CN110672121B (zh) * | 2019-09-19 | 2020-11-20 | 北京控制工程研究所 | 一种控制力矩陀螺框架动态响应测试方法及*** |
CN111191344A (zh) * | 2019-12-09 | 2020-05-22 | 北京动力机械研究所 | 一种运动副摩擦参数的识别方法 |
CN111191344B (zh) * | 2019-12-09 | 2023-09-15 | 北京动力机械研究所 | 一种运动副摩擦参数的识别方法 |
CN112721895A (zh) * | 2021-01-26 | 2021-04-30 | 同济大学 | 基于新型摩擦模型的iehb***主缸液压力估计方法 |
CN112721895B (zh) * | 2021-01-26 | 2022-02-18 | 同济大学 | 基于摩擦模型的iehb***主缸液压力估计方法 |
CN113029194A (zh) * | 2021-02-26 | 2021-06-25 | 北京控制工程研究所 | 一种cmg组合体测试方法、性能评价方法及测试*** |
CN113029194B (zh) * | 2021-02-26 | 2024-03-15 | 北京控制工程研究所 | 一种cmg组合体测试方法、性能评价方法及测试*** |
CN113110305A (zh) * | 2021-04-30 | 2021-07-13 | 华中科技大学 | 一种机电***的摩擦力建模方法及其应用 |
CN113420455A (zh) * | 2021-07-06 | 2021-09-21 | 同济大学 | 一种基于ehb自身特性的主缸液压力估计方法 |
CN113359478A (zh) * | 2021-07-15 | 2021-09-07 | 广东工业大学 | 单自由度直线运动平台非均匀导轨摩擦参数的辨识方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101226068C (zh) | |
CN101226068B (zh) | 2010-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101226068B (zh) | 动态摩擦参数测试***及测试方法 | |
Yang et al. | Dynamic modeling for machine tool thermal error compensation | |
CN104639003B (zh) | 一种交流伺服***的转动惯量辨识方法 | |
CN101821582B (zh) | 用于确定测量机中的动态误差的方法 | |
CN105196294A (zh) | 采用位置测量的可重构机械臂分散控制***及控制方法 | |
Vashisht et al. | Online chatter detection for milling operations using LSTM neural networks assisted by motor current signals of ball screw drives | |
CN103558079A (zh) | 基于并联机构驱动力闭环的多自由度加载方法 | |
CN109739182A (zh) | 一种对冷却***扰动不敏感的主轴热误差补偿方法 | |
CN105880449B (zh) | 模锻压机运行状态的在线预测方法及*** | |
CN103259479A (zh) | 一种永磁同步电机神经网络左逆状态观测方法 | |
CN101691020A (zh) | 一种用于虚拟轴机床刀具运动控制的滑模控制方法 | |
De Moerlooze et al. | On the relationship between normal load and friction force in pre-sliding frictional contacts. Part 2: Experimental investigation | |
CN102288362A (zh) | 一种振动叶片表面非定常压力测试***及测试方法 | |
CN102022966A (zh) | 数控机床导轨丝杆螺母座综合装配精度测量方法及测量装置 | |
CN108225424A (zh) | 机器人关节测试装置及方法 | |
CN111547263B (zh) | 一种飞机操纵***地面动态特性测试方法 | |
Guo et al. | An integrated machine-process-controller model to predict milling surface topography considering vibration suppression | |
Liu et al. | Thermal–mechanical coupling analysis and experimental study on CNC machine tool feed mechanism | |
Zhang et al. | Dynamics analysis and energy consumption modelling based on bond graph: Taking the spindle system as an example | |
Dong et al. | Prediction and compensation of dynamic errors for coordinate measuring machines | |
Tao et al. | High-accurate cutting forces estimation by machine learning with voice coil motor-driven fast tool servo for micro/nano cutting | |
CN101804586A (zh) | 一种虚拟轴机床末端刀具运动位姿的检测方法 | |
CN104569811B (zh) | 马达参数测量方法与马达参数测量*** | |
Zhu et al. | Compilation of program-loading spectrum for milling of a motorized spindle based on cutting force model | |
Burzyński et al. | Dynamics loading by swinging bells—Experimental and numerical investigation of the novel yoke–bell–clapper system with variable geometry |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100407 Termination date: 20140201 |