CN105303017B - 线性移动平台的摩擦力参数鉴别方法 - Google Patents
线性移动平台的摩擦力参数鉴别方法 Download PDFInfo
- Publication number
- CN105303017B CN105303017B CN201410319893.6A CN201410319893A CN105303017B CN 105303017 B CN105303017 B CN 105303017B CN 201410319893 A CN201410319893 A CN 201410319893A CN 105303017 B CN105303017 B CN 105303017B
- Authority
- CN
- China
- Prior art keywords
- linear moving
- moving table
- motor
- friction force
- parametric equations
- 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 abstract description 4
- 230000003068 static effect Effects 0.000 claims abstract description 7
- 238000012850 discrimination method Methods 0.000 claims description 12
- 239000004576 sand Substances 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 238000005259 measurement Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B17/00—Systems involving the use of models or simulators of said systems
- G05B17/02—Systems involving the use of models or simulators of said systems electric
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Automation & Control Theory (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Control Of Position Or Direction (AREA)
- Braking Arrangements (AREA)
Abstract
一种线性移动平台的摩擦力参数鉴别方法,由于不同摩擦力的作用区间取决于两个接触面之间的相对速度,当两个接触面之间的相对速度远大于史崔贝克速度时,此时仅有库仑摩擦力与黏滞摩擦力存在于接触面,所以可以利用此区间的量测扭矩信号而鉴别出库仑摩擦力矩、线性移动平台的摩擦力矩及线性移动平台的等效惯量,当两个接触面之间的相对速度小于史崔贝克速度时,可以配合前述三个已知参数将最大静摩擦力矩与史崔贝克速度鉴别出来。由此,线性移动平台只要进行一次往复运动即可完成所有摩擦力参数的鉴别,以提升实务上的可行性。
Description
技术领域
本发明与线性***有关,尤指一种线性移动平台的摩擦力参数鉴别方法。
背景技术
对于使用滚珠螺杆的自动化设备来说,自动化设备的定位精度主要仰赖于滚珠螺杆的预压力,以消除滚珠螺杆的背隙,同时增加滚珠螺杆的刚性,然而此一预压力将无可避免地增加接触面的摩擦力,使得螺杆轴在速度转换方向的时候会产生象限误差(quadranterror),进而对自动化设备的精密度造成影响。
为了解决前述问题,一般会通过LuGre摩擦力模型来建构出摩擦力矩与速度之间的关系曲线,最后再利用曲线拟合(curve fitting)的方式将相关参数鉴别出来,但是在使用LuGre摩擦力模型的过程中需要进行多次的定速实验,在实务应用上会受到诸多限制而缺乏实用性,另外,在使用曲线拟合方式的过程中会因为未知参数的数量过多而大幅增加鉴别的困难度。
发明内容
本发明的主要目的在于提供一种线性移动平台的摩擦力参数鉴别方法,其不需进行多次定速实验,并且能降低参数鉴别的困难度及提升实务应用的可行性。
为了达成前述目的,本发明的摩擦力参数鉴别方法包含有三个步骤。第一个步骤提供一参数方程式,该参数方程式为
其中的Tm为马达的输出力矩,J为线性移动平台的等效惯量,α为马达的输出轴的角加速度,Tc为库仑摩擦力矩,ω为马达的输出轴的角速度,Ts为最大静摩擦力矩,ωs为史崔贝克速度(Stribeckvelocity),σ2为黏滞摩擦系数;第二个步骤设定在ω远大于ωs的条件下,通过该参数方程式得到J、Tc及σ2,较佳地,该参数方程式可通过弦波速度规划或梯形速度规划取得J、Tc及σ2;第三个步骤设定在ω低于ωs的条件下且配合第二个步骤所得的已知参数,通过该参数方程式而鉴别出Ts及ωs,较佳地,该参数方程式可以通过曲线拟合的方式取得Ts及ωs,或者再转换成一线性方程式得到Ts及ωs,该线性方程式为
p=q-ω2×γ,
其中的p为
ln(Tm-Jα-Tcsgn(ω)-σ2ω),q为ln(Ts-Tc),γ为1/(ωs)2。
由此,本发明的摩擦力参数鉴别方法将该线性移动平台的移动速度区分为高速段区间及低速段区间,如此让该线性移动平台只要进行一次往复运动即可将所有的参数鉴别出来,以有效降低参数鉴别的困难度及显著提升在实务应用上的可行性。
附图说明
为进一步说明本发明的技术内容,以下结合实施例及附图详细说明如后,其中:
图1为本发明的方块图。
图2为本发明的弦波速度规划示意图。
图3为本发明的梯形速度规划示意图。
具体实施方式
请参阅图1,本发明的摩擦力参数鉴别方法包含有步骤a)S1、步骤b)S2,以及步骤c)S3。
步骤a)S1:先通过LuGre摩擦力模型取得一第一方程式,第一方程式为
其中的Tf为线性移动平台的摩擦力矩,Tc为库仑摩擦力矩,ω为马达的输出轴的角速度,Ts为最大静摩擦力矩,ωs为史崔贝克速度(Stribeck velocity),σ2为黏滞摩擦系数,接着同样通过LuGre摩擦力模型取得一第二方程式,第二方程式为Tm=Jα+Tf,其中的Tm为马达的输出力矩,J为线性移动平台的等效惯量,α为马达的输出轴的角加速度,接着将第一、第二方程式结合的后会得到一参数方程式,参数方程式为
步骤b)S2:当ω远大于ωs时,也就是线性移动平台位于高速区段的情况下,这时候的
趋近于0,
所以参数方程式可以简化为Tm=Jα+Tcsgn(ω)+σ2ω,其中的Tm与ω可以直接量测而得,在取得ω的后再经过一次微分可以得到α,这时候有两种方式可以鉴别出J、Tc及σ2。
第一种方式是利用弦波速度规划(如第2图所示)将高速区段的多数个量测信号整理成如下所示的矩阵形式:
并且令Y=AX,
其中的Y为马达的输出力矩所组成的向量,A为马达的输出轴的角加速度及马达的输出轴的角速度所共同组成的矩阵,X为所欲鉴别的参数所组成的向量,此时可以再将前述矩阵改写成下列形式:
并利用最小平方法求出J、Tc及σ2。
第二种方式利用梯形速度规划(如图3所示)定义出ωp、ωn、Tp,以及Tn,其中的ωp为线性移动平台于去程时在定速区段的角速度,且|ωp|>>ωs,ωn为线性移动平台于回程时在定速区段的角速度,且|ωn|>>ωs,Tp为马达在线性移动平台于去程时在定速区段所输出的力矩,Tn为马达在线性移动平台于回程时在定速区段所输出的力矩,由于在定速区段时的角速度α为0,这时候可以将步骤a)的参数方程式改写成如此即可得到
在得到σ2与Tc的后再利用高速区段(ω远大于ωs)的量测信号通过步骤a)的参数方程式求出
步骤c)S3:当ω小于ωs或是接近于ωs时,也就是线性移动平台位于低速区段的情况下,
并不等于0,由于J、Tc及σ2在步骤b)已经鉴别出来,所以在参数方程式内只剩下Ts与ωs为未知参数,这时候有两种方式可以使用:
第一种方式是将未知参数跟步骤b)所取得的已知参数分开来且同时取对数,使步骤a)的参数方程式变成一线性方程式,线性方程式为
p=q-ω2×γ,其中的
p=ln(Tm-Jα-Tcsgn(ω)-σ2ω),q=ln(Ts-Tc),γ=1/(ωs)2,
由于p可以将已知参数代入而求得,ω可以直接量测而得,如此即可求解出q跟γ,在得到q跟γ的后就能够进一步鉴别出Ts与ωs。
第二种方式先将参数方程式改写成下列形式:
接着再以曲线拟合的方式将Ts与ωs鉴别出来,由于此时只有两个未知参数,因而可以大幅降低在曲线拟合过程的困难度。
综上所陈,本发明的摩擦力参数鉴别方法将线性移动平台的移动速度区分为高速段区间及低速段区间,如此让线性移动平台只要进行一次往复运动即可将所有的相关参数鉴别出来,相较于习用技术更能够有效降低参数鉴别的困难度及提升在实务应用上的可行性。
Claims (5)
1.一种线性移动平台的摩擦力参数鉴别方法,包含有下列步骤:
a)提供一参数方程式,该参数方程式如下所示:
其中的Tm为马达的输出力矩,J为线性移动平台的等效惯量,α为马达的输出轴的角加速度,Tc为库仑摩擦力矩,ω为马达的输出轴的角速度,Ts为最大静摩擦力矩,ωs为史崔贝克速度,σ2为黏滞摩擦系数;
b)设定在ω远大于ωs的条件下,这时候的sgn(ω)趋近于0,通过该参数方程式鉴别出J、Tc及σ2;以及,
c)设定在ω小于ωs的条件下且配合步骤b)所取得的已知参数,通过该参数方程式鉴别出Ts及ωs。
2.如权利要求1所述的线性移动平台的摩擦力参数鉴别方法,其中该参数方程式在步骤c)中通过曲线拟合的方式取得Ts及ωs。
3.如权利要求1所述的线性移动平台的摩擦力参数鉴别方法,其中该参数方程式在步骤c)中将未知参数与步骤b)所取得的已知参数区分开来且同时取对数而得到一线性方程式,并通过该线性方程式鉴别出Ts及ωs,该线性方程式为
p=q-ω2×r,其中的p为ln(Tm-Jα-Tcsgn(ω)-σ2ω),q为ln(Ts-Tc),r为1/(ωs)2。
4.如权利要求1所述的线性移动平台的摩擦力参数鉴别方法,其中该参数方程式在步骤b)中利用弦波速度规划鉴别出
J、Tc及σ2,其中的A为Y为
其中的Y为马达的输出力矩所组成的向量,A为马达的输出轴的角加 速度及马达的输出轴的角速度所共同组成的矩阵。
5.如权利要求1所述的线性移动平台的摩擦力参数鉴别方法,其中该参数方程式在步骤b)中利用梯形速度规划鉴别出
J、Tc及σ2,
其中的ωp为该线性移动平台于去程时在定速区段的角速度,|ωp|>>ωs,ωn为该线性移动平台于回程时在定速区段的角速度,|ωn|>>ωs,Tp为马达在该线性移动平台于去程时在定速区段所输出的力矩,Tn为马达在该线性移动平台于回程时在定速区段所输出的力矩。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW103122131 | 2014-06-26 | ||
TW103122131A TW201600846A (zh) | 2014-06-26 | 2014-06-26 | 線性移動平台之摩擦力參數鑑別方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105303017A CN105303017A (zh) | 2016-02-03 |
CN105303017B true CN105303017B (zh) | 2018-04-20 |
Family
ID=53441559
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410319893.6A Active CN105303017B (zh) | 2014-06-26 | 2014-07-07 | 线性移动平台的摩擦力参数鉴别方法 |
Country Status (6)
Country | Link |
---|---|
US (1) | US20150377726A1 (zh) |
JP (1) | JP5860935B2 (zh) |
KR (1) | KR101632940B1 (zh) |
CN (1) | CN105303017B (zh) |
DE (1) | DE102014111605A1 (zh) |
TW (1) | TW201600846A (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE112019002515T5 (de) * | 2018-05-17 | 2021-02-25 | Arcus Technology, Inc. | Zustands- bzw. gesundheitsmanagement eines bewegungssystems unter verwendung vorhandener servoantriebsvariablen |
CN110823338B (zh) * | 2018-08-10 | 2021-05-18 | 河南工业大学 | 基于底面单圈标准差对数模型的粮仓检测方法及*** |
CN110266232A (zh) * | 2019-06-20 | 2019-09-20 | 哈尔滨工业大学 | 用于电机伺服***动态转动惯量辨识的方法和装置、电机伺服*** |
CN113110305B (zh) * | 2021-04-30 | 2022-07-05 | 华中科技大学 | 一种机电***的摩擦力建模方法及其应用 |
TWI778729B (zh) * | 2021-07-30 | 2022-09-21 | 新代科技股份有限公司 | 傳動裝置及其控制方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101122779A (zh) * | 2007-09-13 | 2008-02-13 | 北京航空航天大学 | 一种利用蚁群算法辨识高精度伺服***摩擦参数的方法 |
CN101226068A (zh) * | 2008-02-01 | 2008-07-23 | 西安电子科技大学 | 动态摩擦参数测试***及测试方法 |
CN101515778A (zh) * | 2009-04-03 | 2009-08-26 | 哈尔滨工程大学 | 无刷直流力矩电机位置伺服***摩擦力矩自适应补偿方法 |
CN101846975A (zh) * | 2010-05-28 | 2010-09-29 | 北京理工大学 | 带有动态摩擦补偿的伺服***自适应鲁棒控制器 |
CN102269638A (zh) * | 2011-04-27 | 2011-12-07 | 中国科学院光电技术研究所 | 伺服转台LuGre模型摩擦参数及转动惯量的一体化测量方法 |
CN103344243A (zh) * | 2013-07-02 | 2013-10-09 | 北京航空航天大学 | 一种航空遥感惯性稳定平台摩擦参数辨识方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SG84531A1 (en) * | 1999-05-20 | 2001-11-20 | Univ Singapore | Disturbance attenuation in a precision servomechanism by a frequency-separated acceleration soft sensor |
DE102008021848B4 (de) * | 2008-05-02 | 2017-01-12 | Volkswagen Ag | Verfahren zur Berücksichtigung der statischen und dynamischen Reibung in einem System und System |
CN102208891B (zh) * | 2010-11-18 | 2013-04-24 | 东南大学 | 基于摩擦和扰动补偿的pmsm伺服***控制方法 |
KR101915191B1 (ko) * | 2011-09-19 | 2018-11-05 | 한화지상방산 주식회사 | 회전축 제어 시스템 |
JP5924944B2 (ja) * | 2012-01-10 | 2016-05-25 | 株式会社ミツトヨ | すべり案内装置 |
JP6214948B2 (ja) * | 2013-07-12 | 2017-10-18 | 三菱重工業株式会社 | 摩擦補償装置及び摩擦補償方法並びにサーボ制御装置 |
-
2014
- 2014-06-26 TW TW103122131A patent/TW201600846A/zh unknown
- 2014-07-07 CN CN201410319893.6A patent/CN105303017B/zh active Active
- 2014-07-29 KR KR1020140096265A patent/KR101632940B1/ko active IP Right Grant
- 2014-07-30 JP JP2014154916A patent/JP5860935B2/ja active Active
- 2014-08-14 DE DE102014111605.4A patent/DE102014111605A1/de not_active Withdrawn
- 2014-09-05 US US14/478,940 patent/US20150377726A1/en not_active Abandoned
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101122779A (zh) * | 2007-09-13 | 2008-02-13 | 北京航空航天大学 | 一种利用蚁群算法辨识高精度伺服***摩擦参数的方法 |
CN101226068A (zh) * | 2008-02-01 | 2008-07-23 | 西安电子科技大学 | 动态摩擦参数测试***及测试方法 |
CN101515778A (zh) * | 2009-04-03 | 2009-08-26 | 哈尔滨工程大学 | 无刷直流力矩电机位置伺服***摩擦力矩自适应补偿方法 |
CN101846975A (zh) * | 2010-05-28 | 2010-09-29 | 北京理工大学 | 带有动态摩擦补偿的伺服***自适应鲁棒控制器 |
CN102269638A (zh) * | 2011-04-27 | 2011-12-07 | 中国科学院光电技术研究所 | 伺服转台LuGre模型摩擦参数及转动惯量的一体化测量方法 |
CN103344243A (zh) * | 2013-07-02 | 2013-10-09 | 北京航空航天大学 | 一种航空遥感惯性稳定平台摩擦参数辨识方法 |
Non-Patent Citations (1)
Title |
---|
基于LuGre模型实现精密伺服转台摩擦参数辨识及补偿;于伟 等;《光学精密工程》;20111130;第19卷(第11期);第2738-2739页 * |
Also Published As
Publication number | Publication date |
---|---|
DE102014111605A1 (de) | 2015-12-31 |
TW201600846A (zh) | 2016-01-01 |
JP2016009481A (ja) | 2016-01-18 |
KR20160001557A (ko) | 2016-01-06 |
CN105303017A (zh) | 2016-02-03 |
KR101632940B1 (ko) | 2016-06-23 |
TWI481847B (zh) | 2015-04-21 |
US20150377726A1 (en) | 2015-12-31 |
JP5860935B2 (ja) | 2016-02-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105303017B (zh) | 线性移动平台的摩擦力参数鉴别方法 | |
CN106289708B (zh) | 用于捕获轨迹风洞试验的运动机构的位姿标定方法 | |
CN109249397B (zh) | 一种六自由度机器人动力学参数辨识方法和*** | |
CN103869834B (zh) | 基于经验模态法的三轴气浮台质心智能调节方法 | |
CN108297093A (zh) | 一种机械臂动力学参数的分步辨识方法 | |
CN104535577B (zh) | 一种工件质量损失检测方法 | |
KR20170044987A (ko) | 저크가 제한된 궤적 생성 방법 | |
Fortin-Côté et al. | Improving cable driven parallel robot accuracy through angular position sensors | |
CN109631847A (zh) | 一种基于点云数据的含螺纹目标物位姿的计算方法 | |
CN115586724A (zh) | 一种齿轮巡检机器人***自适应分数阶全局滑模控制方法 | |
US20170261529A1 (en) | Method for identifying friction parameter for linear module | |
CN106813570A (zh) | 基于线结构光扫描的长圆柱形物体三维识别与定位方法 | |
CN111047002A (zh) | 一种基于bp神经网络的校直行程预测方法 | |
WO2023050188A1 (zh) | 一种基于三维建模技术的数字煤场管理方法及*** | |
JP2019046336A (ja) | マルチコプタ及びその制御方法 | |
CN108227493A (zh) | 一种机器人轨迹跟踪方法 | |
CN106643566B (zh) | 一种自动测量轮毂轮辐初始角度的方法 | |
CN107547028A (zh) | 一种电机转子角度检测的方法及装置 | |
CN108120450A (zh) | 一种静止状态的判断方法及装置 | |
CN109932925B (zh) | 多侦察点仿真方法 | |
US20150300842A1 (en) | Device and Method For Determining the Change of Position of a 3D Measuring Head | |
CN107966701B (zh) | 一种雷达倾角测量方法及装置 | |
KR102059305B1 (ko) | 엘리베이터용 패널 제조 공정의 제어 방법 및 장치 | |
CN104657595A (zh) | 一种单颗粒曳力模型系数标定方法 | |
CN104280046B (zh) | 飞行器运动参数测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |