CN103926875A - 一种滚珠丝杠进给***摩擦补偿方法 - Google Patents
一种滚珠丝杠进给***摩擦补偿方法 Download PDFInfo
- Publication number
- CN103926875A CN103926875A CN201410157347.7A CN201410157347A CN103926875A CN 103926875 A CN103926875 A CN 103926875A CN 201410157347 A CN201410157347 A CN 201410157347A CN 103926875 A CN103926875 A CN 103926875A
- Authority
- CN
- China
- Prior art keywords
- friction
- population
- speed
- individual
- model
- 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
Landscapes
- Feedback Control In General (AREA)
Abstract
本发明公开一种滚珠丝杠进给机构摩擦补偿方法。在分析***实际摩擦现象基础上,提出摩擦力矩在速度区间的分段模型,采用斯特克贝里模型对低速区间摩擦建模,对数模型对高速区间摩擦进行建模。利用混合遗传算法辨识出摩擦模型未知参数,并基于该辨识出的摩擦模型估测出实时速度下的摩擦补偿量,通过补偿量前馈的方式实现了对***摩擦的补偿。本发明提出的分段摩擦模型精确地描述了滚珠丝杠进给***中的摩擦,且设计的混合遗传算法准确地摩擦模型的参数进行辨识,利用基于摩擦模型前馈的补偿方式,提高了滚珠丝杠进给***的运动精度,可广泛应用于各类含滚珠丝杠进给机构的***。
Description
技术领域
本发明涉及一种机床伺服进给***的摩擦补偿方法,提出一种基于分段摩擦模型的摩擦误差补偿方法,属于高精度机床的控制技术领域。
背景技术
随着工业生产对机床加工精度的要求越来越高,消除摩擦对机床加工核心部件—滚珠丝杠机构运动精度的影响已成为提高机床加工精度的重要手段。作为一种干扰因素,摩擦能使高精度进给机构在位置跟踪时出现“削顶”现象,速度跟踪时出现“过零畸变”现象。在位置伺服控制时,摩擦还是导致“死区”和“极限换振荡”的主要原因,因而对于精密加工,超精密加工领域,摩擦干扰的影响不容忽视。
为了消除非线性摩擦对进给***运动性能的影响,国内外学者在研究摩擦补偿控制策略方面做了如下研究:
1)基于库伦摩擦+粘性摩擦模型的补偿法
该方法设计了基于速度和位置信号的摩擦参数观测器,实现了对***的摩擦补偿。此补偿方法优点在于模型简单,但是所选用的摩擦模型属于静态摩擦模型,其无法描述摩擦完整特性,因而在实际应用中存在很大局限。
2)基于力矩传感器反馈控制的补偿法
该思想是在机器人基部安装力矩传感器,通过基部测量的力矩信号计算出各关节上的净力矩,从而构成力矩反馈回路。力矩反馈控制虽然不依赖于模型,但由于传感器价格高,安装困难等原因,其应用亦受到限制。
3)基于模糊控制的摩擦干扰补偿法
该方法针对摩擦的非线性环节,采用模糊聚簇技术进行建模,从大量数据中提取出摩擦环节的模糊模型并进行补偿。但是模糊控制学习能力不强,设计时控制规则过于依赖经验和专家知识,此缺点限制了模糊逻辑控制在实际补偿中应用。
目前在消除摩擦干扰的各种方法中,基于摩擦模型补偿的方法因为其成本较低,实现较为方便,成为机床进给机构中消除摩擦影响最为常用的方法。
发明内容
本发明适用于基于模型的摩擦补偿方式,提出了滚珠丝杠进给机构中摩擦力矩关于速度的分段模型,特别是在***高速运行阶段,通过对实测的摩擦数据进行分析,并结合现有的相关理论,提出了高速阶段对数摩擦模型,能够反映“粘性摩擦增速随速度增加而降低”这一现象,同时利用混合遗传算法对模型中的待辨识参数进行精确辨识,最后利用摩擦补偿量前馈输入的方式对***进行了摩擦误差的补偿,提高了滚珠丝杠进给机构的运动精度。
为了实现上述目的,本发明采用了如下的技术方案:
一种滚珠丝杠进给***摩擦补偿方法,其特征在于步骤如下:
步骤1首先控制滚珠丝杠进给***,使滚珠丝杠分别以从小到大的速度v1、v2、v3、…、v53匀速运转,分别获取速度为v1、v2、v3、…、v53时***的摩擦力矩T1、T2、T3、…、T53,并将速度和摩擦力矩组成***摩擦数据离散样本序列St={(v1、T1)、(v2、T2)、…、(v53、T53)},
步骤2设:vz为低、高速区间最佳分界速度,分别对低速、高速运行时***的摩擦建模,
步骤3使用步骤1得到的***摩擦数据离散样本序列St={(v1、T1)、(v2、T2)、…、(v53、T53)},分别对低速、高速运行时***的摩擦模型中的待辨识参数进行辨识,得到辨识后的低速、高速运行时***的摩擦模型,
步骤4将电机的期望速度Vd代入辨识后的低速、高速运行时***的摩擦模型,得到期望速度为Vd时***所需的摩擦力矩补偿量T(Vd),并求得等效的速度补偿增量Vf,
步骤5将速度补偿增量Vf叠加到期望速度Vd上,获得补偿速度Vdf,并得到所述补偿速度Vdf对应的电机控制电压信号Udf,再将电机控制电压信号Udf转换为电机驱动电压U用于电机驱动。
在|v|≤vz的低速区间,摩擦模型采用斯特克贝里模型;在|v|>vz的高速区间,采用对数摩擦模型,
所述斯特克贝里模型为:
所述对数摩擦模型为:
自变量v为滚珠丝杠的瞬时速度,因变量T为摩擦力矩值;模型中待辨识参量:Tc为库伦摩擦力矩,Ts为最大静摩擦力矩,vs为斯特克贝里速度,σ2为粘滑系数,ξ1为粘性摩擦斜率因子,ξ2为粘性摩擦斜率修正因子,ξ3为速度连接系数,sgn(v)为速度的符号函数,满足:
对低速、高速运行时***摩擦模型的待辨识参数进行辨识的方法为:
步骤3.1令i=4,i为循环指针,
步骤3.2选择vi为分界速度,则低速摩擦数据为Stl={(v1、T1)}、…、(vi、Ti)},高速摩擦数据为Sth={(vi+2、Ti+2)、…、(v53、T53)},基于Stl及Sth,采用混合遗传算法辨识得出斯特克贝里摩擦模型待辨识参数Tc、Ts、vs、σ2的辨识值分别为以及对数摩擦模型待辨识参数ξ1、ξ2、ξ3的辨识值分别为 则分界速度为vi时***摩擦模型式如下:
步骤3.3将St中每个样本点的速度vj代入式(4),j为1到53的自然数,得速度vj时***的摩擦力矩模型估计值计算与vj对应的***实际摩擦力矩Tj与模型估计摩擦力矩误差的平方和Ei,即Ei为分界速度为vi时摩擦模型对实际摩擦数据的拟合误差,Ei表达式为:
步骤3.4如果i≥53,则进入步骤3.5;否则,令i=i+1,返回步骤3.2,
步骤3.5取E4、E5、…、E53最小值时的所选分界速度为最佳分界速度vz,取最佳分界速度为vz时所得模型待辨识参数的辨识值为待辨识参数最终估计值。
待辨识参数Tc、Ts、vs、σ2、ξ1、ξ2、ξ3的混合遗传辨识算的流程如下:
1)编码
首先将所有待辨识参数可行解范围设定为[-10,10],利用16位二进制序列对每个待辨识参数进行编码,编码精度δ为:
编码公式为:
G依次为Tc、Ts、vs、σ2、ξ1、ξ2及ξ3,在w依次取值为1、2、…、7时分别为Tc、Ts、vs、σ2、ξ1、ξ2及ξ3的第k个基因位,分别是以G所表示的Tc、Ts、vs、σ2、ξ1、ξ2或ξ3的基因型,k为1~16的正整数,的数值为0或1,
2)初始化种群
设由斯特克贝里摩擦模型及对数摩擦模型的待辨识参数组合成的种群个体的表现型为:
X=[Tc Ts vs σ2 ξ1 ξ2 ξ3] (8)
由式(7)所示的编码规则,与X对应的种群个体的基因型为
设定种群个体的数目为300,随机产生初始种群其中u为1~300的正整数,u为种群个体的序数,为第0代种群(即初始种群)的第u个种群个体的表现型,
设定种群的最大进化代数200,交叉概率0.6,变异概率0.3,局部搜索次数为5,
3)种群个体适应度值求取
第c代种群Pc为:其中c为0~200的正整数,u为1~300的正整数,为c代种群中第u个个体,设定种群个体的优化目标函数:
上式,Ti为速度为vi时***的实际摩擦力矩,为速度为vi时***的摩擦力矩估计值,其由如下所示的摩擦模型估计所得:
式中,分别为Tc、Ts、vs、σ2、ξ1、ξ2、ξ3在种群个体下的辨识值,
种群个体的适应度函数与目标函数的映射关系,如下:
式中,为c代种群的第u个个体的优化目标函数值,为c代种群的第u个个体的适应度值,
4)选择操作
采用赌盘选择策略确定第c代种群中种群个体的被选择与否,赌盘选择策略如下:
(a)Pc中所有个体的适应度值总和Fc:
式中,Fc为Pc中所有个体适应度值总和,
(b)Pc中第u个个体被选择的概率Ru:
式中,u为1~300的正整数,
(c)Pc中第u个种群个体选择与否判断:
随机生成一个[0,1]范围内的数r,如果r大于选择概率Ru,则该种群个体被选择,否则,该种群个体被舍弃,
5)交叉操作
随机不重复地从种群Pc中选择两个种群个体其中a,b为1~300的不相等的正整数,种群个体的基因型分别为112位,由交叉概率0.6随机选择并确定112位基因型中的一位基因位并将所选基因位位数记录h,将种群个体中的h至第112位基因位的部分与种群个体中的h至第112位基因位的部分进行交叉互换,用交叉后得到的种群个体分别取代原种群个体
6)变异操作
第c代种群Pc中个体的基因型有112位基因位,从第1个基因位开始,产生随机数rand,将此随机数与概率值为0.3的变异概率比较,当rand<0.3时,该基因位值取反,依次重复此操作,直到第112位基因位,
7)局部搜索操作
对种群Pc中个体采用爬山算法实现局部搜索,具体步骤为:随机交换种群个体基因型中两个基因位值,得新的种群个体如果适应度值用取代反之,当前种群个体仍为局部搜索操作的次数为5次,
8)最优种群个体获取
c依次取值为0、1、2、……、200,第200代种群中最优个体所对应的基因型为
9)解码,输出最优解
对最优种群个体的基因型进行解码,解码公式为:
G为Tc、Ts、vs、σ2、ξ1、ξ2的辨识值 分别是以G所表示的 的基因型,在w依次取值为1、2、…、7时分别为 或的第f个基因位,f为1~16的正整数,的值为0或1。
与现有技术相比,本发明具有以下优点:
(1)本发明提出的速度区间分段摩擦模型,解决了以往单一模型无法精确描述整个速度区间摩擦的问题。本发明提出对数模型T=(ln(v+ξ2sgn(v))ξ1+ξ3)sgn(v)对***高速运动时的摩擦进行建模,当ξ1>0且ξ2>0,对数模型所描述的摩擦力矩关于速度的一阶导数dT/dv恒正,二阶导数d2T/dv2恒负,满足“粘滑阶段的摩擦增速率随速度增加而降低”的结论,因而从摩擦模型构建的角度,提高了摩擦误差的补偿精度。
(2)本发明设计的混合遗传算法针对经典遗传算法局部寻优能力弱的缺点,在种群个体经过遗传操作后对其进行多次基于爬山搜索策略的局部搜索操作。因为此局部搜索操作是以选择、交叉及变异等遗传操作后所得较优的种群个体为基础,且每次的局部搜索操作都沿着使此种群较优个体适应度值增加的方向进行,所以经过若干代,且每代多次的局部寻优操作,混合遗传算法最终能够搜索到比遗传算法适应度值更高的种群最优个体。在对非线性摩擦模型未知参数的辨识时,混合遗传算法极大地提高了其参数辨识的精度。
(3)本发明采用的摩擦补偿量前馈输入方式,在无需速度,位置反馈的开环***中,根据分段摩擦模型实时估测指令速度所需的补偿速度Vf,并将其叠加到指令速度中,以实现对***摩擦干扰的补偿,该方法操作简单,能够方便地在伺服控制器中实现。
(4)本发明不仅可以应用于高精度机床进给机构的摩擦补偿,也可以应用于其他高精度、超高精度设备的摩擦补偿。
附图说明
图1滚珠丝杠进给***摩擦补偿框图。
图2控制器框图。
图3执行机构框图。
图4混合遗传算法流程图。
图5不同假设分界点时分段摩擦模型拟合误差图。
图6基于混合遗传算法辨识的分段摩擦模型对***摩擦数据的拟合图。
图7指令速度、有/无摩擦补偿速度输出图。
图8指令位移、有/无摩擦补偿位移输出图。
具体实施方式
本发明为基于模型的摩擦前馈补偿方式,该补偿策略首先是采集***的摩擦数据,同时建立能够精确描述***摩擦的数学模型,在此基础上利用混合遗传算法对所建摩擦模型的未知参数进行离线辨识。然后,在控制器中利用此摩擦模型估计***指令速度下的摩擦干扰,并在指令速度上加上相应的速度补偿量。最后利用补偿后的信号控制进给机构运转。
本发明的具体实施方式的具体步骤如下:
1、指定滚珠丝杠匀速运转时的速度vi,获得此时电机的电枢电流I(vi)及力矩常量kt。
速度为vi时,电机输出转矩与电枢电流满足如下关系:
Mi=kt·I(vi) (1)
电机空载转动时的力矩平衡式:
Jα=Mi-Ti (2)
式中,J为***等效转动惯量。
匀速运动时,加速度α=0,电机驱动力矩Mi近似等于滚珠丝杠摩擦力矩Ti。
指定53个不同速度,初始速度v1为7.5Rpm,即0.785rad.s-1,相邻速度的差值为7.5Rpm,得***速度区间上速度由小到大增加的摩擦样本点序列St={(v1、T1)、(v2、T2)、…、(v53、T53)}。
2、进行摩擦在速度区间的分段建模,低速区间摩擦满足斯特克贝里模型,其表达式如下:
式中,Tc为库伦摩擦力矩,Ts为最大静摩擦力矩,vs为斯特克贝里速度,σ2为粘滑系数。
高速粘滑区间摩擦力矩与速度的关系满足对数模型,其表达式为:
自变量v为瞬时速度,因变量T为***的摩擦力矩,ξ1为粘性摩擦斜率因子,ξ2为粘性摩擦斜率修正因子,ξ3为速度连接系数。sgn(v)为速度的符号函数,满足:
对式(4)求一阶导,得摩擦力矩变化率为:
由式(6)可得摩擦力矩T关于速度v的二阶导倒数如下:
由式(6)、式(7)可得,存在ξ1>0,ξ2>0,使dT/dv恒为正,d2T/dv2恒为负,因而此模型能够描述“高速阶段的粘性摩擦增加趋势随速度增加减缓”这一现象。
综上,当求得***低、高速区间的最佳分解速度vz时,***的摩擦模型如下:
3、假设St中第i个样本点的速度vi为低、高速区间的最佳分界速度,则低速摩擦数据为Stl={(v1、T1)、…、(vi、Ti)},高速摩擦数据Sth={(vi+1、Ti+1)、…、(v53、T53),其中,i为循环指针,为4~53的正整数。
4、待辨识参数Tc、Ts、vs、σ2、ξ1、ξ2、ξ3的混合遗传辨识算的流程如下:
1)待辨识参数编码
首先将所有待辨识参数可行解范围设定为[-10,10],利用16位二进制序列对每个待辨识参数进行编码,编码精度δ为:
编码公式为:
G依次为Tc、Ts、vs、σ2、ξ1、ξ2及ξ3,在w依次取值为1、2、…、7时分别为Tc、Ts、vs、σ2、ξ1、ξ2及ξ3的第k个基因位,分别是以G所表示的Tc、Ts、vs、σ2、ξ1、ξ2或ξ3的基因型,k为1~16的正整数,的数值为0或1,
2)混合遗传种群初始化
设由斯特克贝里摩擦模型及对数摩擦模型的待辨识参数组合成的种群个体的表现型为:
X=[Tc Ts vs σ2 ξ1 ξ2 ξ3] (11)
由式(7)所示的编码规则,与X对应的种群个体的基因型为
设定种群个体的数目为300,随机产生初始种群其中u为1~300的正整数,u为种群个体的序数,为第0代种群(即初始种群)的第u个种群个体的表现型,
设定种群的最大进化代数200,交叉概率0.6,变异概率0.3,局部搜索次数为5,
3)种群个体适应度值求取
第c代种群Pc为:其中c为0~200的正整数,u为1~300的正整数,为c代种群中第u个个体,设定算法的优化目标函数:
上式,Ti为速度为vi时***的实际摩擦力矩,为速度为vi时***的摩擦力矩估计值,此估计值由如下所示的摩擦模型估计所得:
式中,分别为Tc、Ts、vs、σ2、ξ1、ξ2、ξ3在种群个体下的辨识值,
种群个体的适应度函数与目标函数的映射关系如下:
式中,为c代种群的第u个个体的优化目标函数值,为c代种群的第u个个体的适应度值,
4)选择操作
采用赌盘选择策略确定第c代种群中种群个体被选择与否,赌盘选择策略如下:
(b)Pc中所有个体的适应度值总和Fc:
式中,Fc为Pc中所有个体适应度值总和,
(b)Pc中第u个个体被选择的概率Ru:
式中,u为1~300的正整数,
(d)Pc中第u个种群个体选择与否判断:
随机生成一个[0,1]范围内的数r,如果r大于选择概率Ru,则该种群个体被选择,否则,该种群个体被舍弃,
5)交叉操作
随机不重复地从种群Pc中选择两个种群个体其中a,b为1~300的不相等的正整数,种群个体的基因型分别为112位,由交叉概率0.6随机选择并确定112位基因型中的一位基因位并将所选基因位位数记录为h,将种群个体中的h至第112位基因位的部分与种群个体中的h至第112位基因位的部分进行交叉互换,用交叉后得到的种群个体分别取代原种群个体
6)变异操作
第c代种群Pc中个体的基因型有112位基因位,从第1个基因位开始,产生随机数rand,将此随机数与概率值为0.3的变异概率比较,当rand<0.3时,该基因位值取反,依次重复此操作,直到第112位基因位,
7)局部搜索操作
对种群Pc中个体采用爬山算法实现局部搜索,具体步骤为:随机交换种群个体基因型中两个基因位值,得新的种群个体如果适应度值用取代反之,当前种群个体仍为局部搜索操作的次数为5次,
8)最优种群个体获取
c依次取值为0、1、2、……、200,第200代种群中最优个体所对应的基因型为
9)解码,输出最优解
对最优种群个体的基因型进行解码,解码公式为:
G为参数Tc、Ts、vs、σ2、ξ1、ξ2所对应的辨识值 而在w依次取值为1、2、…、7时分别为 或的第f个基因位,f为1~16的正整数,的值为0或1。
混合遗传算法的流程图如附图4所示。
5、将步骤4中各待辨识参数的估计值带入式(8)中,得分界速度为vi时***摩擦模型式如下:
将St中每个样本点的速度vj代入式(18),j为1~53的正整数,得速度vj时***的摩擦力矩模型估计值则vj下***实际摩擦力矩Tj与模型估计摩擦力矩误差的平方和Ei:
6、如果i≥53,则取E4、E5、…、E53序列中的最小值对应的分界速度为最佳分界速度vz,,如附图5所示,此时的参数辨识值即为所求最终摩擦模型待辨识参数的辨识值。否则,令i=i+1,返回3。
最佳分界速度为vz时利用混合遗传算法对***摩擦数据拟合效果如图6所示。
7、采用基于摩擦模型的前馈补偿控制方式消除***的摩擦干扰。在控制器中,设定指令为Vd,由上述辨识所得的摩擦模型估计出摩擦补偿量T(Vd)。
指令速度为Vd时的摩擦补偿量T(Vd)与等效的补偿速度Vf如下:
Vf=KfT(Vd) (20)
上式,Kf为力矩-速度转换系数,Vf为与摩擦补偿量等效的速度补偿量。
8、将速度补偿量Vf叠加到指令速度Vd上,得期望输出速度为Vd时,***所需的实际速度指令Vdf。
Vdf=Vd+Vf (21)
控制器向驱动器输出的与补偿速度对应的补偿后电压控制信号,如下:
Udf=KvVdf (22)
Kv为速度-电压转换系数,Udf为实际电压控制信号。
控制器中驱动放大电路对补偿后的控制信号Udf进行功率放大,得到输入给电机的驱动电压U为:
U=KpUdf (23)
Kp为放大系数,U为经过补偿后的实际电机驱动电压。
在此驱动电压U的作用下,电机驱动滚珠丝杠转动,由此消除摩擦干扰Tf的影响。
附图7给出的是指令速度Vd为20sin(0.4pi·t)mm/s的滚珠丝杠速度的理论输出、无摩擦补偿输出、有摩擦补偿输出结果对比图。附图8给出了此速度下位置输出的理论曲线、无摩擦补偿输出曲线、有摩擦补偿输出曲线结果对比图,由图可知,采用摩擦补偿的速度、效果明显优于未采用摩擦补偿的速度、位移输出。
Claims (4)
1.一种滚珠丝杠进给***的摩擦补偿方法,其特征在于步骤如下:
步骤1首先控制滚珠丝杠进给***,使滚珠丝杠分别以从小到大的速度v1、v2、v3、…、v53匀速运转,分别获取速度为v1、v2、v3、…、v53时***的摩擦力矩T1、T2、T3、…、T53,并将速度和摩擦力矩组成***摩擦数据离散样本序列St={(v1、T1)、(v2、T2)、…、(v53、T53)},
步骤2设vz为低、高速区间最佳分界速度,分别对低速、高速运行时***的摩擦建模,
步骤3使用步骤1得到的***摩擦数据离散样本序列St={(v1、T1)、(v2、T2)、…、(v53、T53)},分别对低速、高速运行时***的摩擦模型中的待辨识参数进行辨识,得到辨识后的低速、高速运行时***的摩擦模型,
步骤4将电机的期望速度Vd代入辨识后的低速、高速运行时***的摩擦模型,得到期望速度为Vd时***所需的摩擦力矩补偿量T(Vd),并求得等效的速度补偿增量Vf,
步骤5将速度补偿增量Vf叠加到期望速度Vd上,获得补偿速度Vdf,并得到所述补偿速度Vdf对应的电机控制电压信号Udf,再将电机控制电压信号Udf转换为电机驱动电压U用于电机驱动。
2.根据权利要求1所述的滚珠丝杠进给***的摩擦补偿方法,其特征在于:在|v|≤vz的低速区间,摩擦模型采用斯特克贝里模型;在|v|>vz的高速区间,采用对数摩擦模型,
所述斯特克贝里模型为:
所述对数摩擦模型为:
自变量v为滚珠丝杠的瞬时速度,因变量T为摩擦力矩值;模型中待辨识参量:Tc为库伦摩擦力矩,Ts为最大静摩擦力矩,vs为斯特克贝里速度,σ2为粘滑系数,ξ1为粘性摩擦斜率因子,ξ2为粘性摩擦斜率修正因子,ξ3为速度连接系数,sgn(v)为速度的符号函数,满足:
。
3.根据权利要求2所述的滚珠丝杠进给***的摩擦补偿方法,其特征在于:对低速、高速运行时***摩擦模型的待辨识参数进行辨识的方法为:
步骤3.1令i=4,i为循环指针,
步骤3.2选择vi为分界速度,则低速摩擦数据为Stl={(v1、T1)}、…、(vi、Ti)},高速摩擦数据为Sth={(vi+2、Ti+2)、…、(v53、T53)},基于Stl及Sth,采用混合遗传算法辨识得出斯特克贝里摩擦模型待辨识参数Tc、Ts、vs、σ2的辨识值分别为以及对数摩擦模型待辨识参数ξ1、ξ2、ξ3的辨识值分别为 则分界速度为vi时***摩擦模型式如下:
步骤3.3将St中每个样本点的速度vj代入式(4),j为1到53的自然数,得速度vj时***的摩擦力矩模型估计值计算与vj对应的***实际摩擦力矩Tj与模型估计摩擦力矩误差的平方和Ei,Ei即为分界速度为vi时摩擦模型对实际摩擦数据的拟合误差,Ei表达式为:
步骤3.4如果i≥53,则进入步骤3.5;否则,令i=i+1,返回步骤3.2,
步骤3.5取E4、E5、…、E53最小值时的分界速度为最佳分界速度vz,取最佳分界速度为vz时所得模型待辨识参数的辨识值为待辨识参数最终估计值。
4.根据权利要求3所述的滚珠丝杠进给***的摩擦补偿方法,其特征在于:对参数Tc、Ts、vs、σ2、ξ1、ξ2、ξ3进行辨识的混合遗传算法步骤如下:
1)待辨识参数编码
首先将所有待辨识参数可行解范围设定为[-10,10],利用16位二进制序列 对每个待辨识参数进行编码,编码精度δ为:
编码公式为:
G依次为Tc、Ts、vs、σ2、ξ1、ξ2及ξ3,在w依次取值为1、2、…、7时分别为Tc、Ts、vs、σ2、ξ1、ξ2及ξ3的第k个基因位,分别是以G所表示的Tc、Ts、vs、σ2、ξ1、ξ2或ξ3的基因型,k为1~16的正整数,的数值为0或1,
2)混合遗传种群初始化
设由斯特克贝里摩擦模型及对数摩擦模型的待辨识参数组合成的种群个体的表现型为:
X=[Tc Ts vs σ2 ξ1 ξ2 ξ3] (8)
由式(7)所示的编码规则,与X对应的种群个体的基因型为
设定种群个体的数目为300,随机产生初始种群其中u为1~300的正整数,u为种群个体的序数,为第0代种群(即初始种群)的第u个种群个体的表现型,
设定种群的最大进化代数200,交叉概率0.6,变异概率0.3,局部搜索次数为5,
3)种群个体适应度值求取
第c代种群Pc为:其中c为0~200的正整数,u为1~300的正整数,为c代种群中第u个个体,设定算法的优化目标函数:
上式,Ti为速度为vi时***的实际摩擦力矩,为速度为vi时***的摩擦力矩估计值,此估计值由如下所示的摩擦模型估计所得:
式中,分别为Tc、Ts、vs、σ2、ξ1、ξ2、ξ3在种群个体下的辨识值,
种群个体的适应度函数与目标函数的映射关系如下:
式中,为c代种群的第u个个体的优化目标函数值,为c代种群的第u个个体的适应度值,
4)选择操作
采用赌盘选择策略确定第c代种群中种群个体 被选择与否,赌盘选择策略如下:
(a)Pc中所有个体的适应度值总和Fc:
式中,Fc为Pc中所有个体适应度值总和,
(b)Pc中第u个个体被选择的概率Ru:
式中,u为1~300的正整数,
(c)Pc中第u个种群个体选择与否判断:
随机生成一个[0,1]范围内的数r,如果r大于选择概率Ru,则该种群个体被选择,否则,该种群个体被舍弃,
5)交叉操作
随机不重复地从种群Pc中选择两个种群个体其中a,b为1~300的不相等的正整数,种群个体的基因型分别为112位,由交叉概率0.6随机选择并确定112位基因型中的一位基因位并将所选基因位位数记录为h,将种群个体中的h至第112位基因位的部分与种群个体中的h至第112位基因位的部分进行交叉互换,用交叉后得到的种群个体分别取代原种群个体
6)变异操作
第c代种群Pc中个体的基因型有112位基因位,从第1个基因位开始,产生随机数rand,将此随机数与概率值为0.3的变异概率比较,当rand<0.3时,该基因位值取反,依次重复此操作,直到第112位基因位,
7)局部搜索操作
对种群Pc中个体采用爬山算法实现局部搜索,具体步骤为:随机交换种群个体基因型中两个基因位值,得新的种群个体如果适应度值 用取代反之,当前种群个体仍为局部搜索操作的次数为5次,
8)最优种群个体获取
c依次取值为0、1、2、……、200,第200代种群中最优个体所对应的基因型为
9)解码,输出最优解
对最优种群个体的基因型进行解码,解码公式为:
G为参数Tc、Ts、vs、σ2、ξ1、ξ2所对应的辨识值 而在w依次取值为1、2、…、7时分别为或的第f个基因位,f为1~16的正整数,的值为0或1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410157347.7A CN103926875B (zh) | 2014-04-18 | 2014-04-18 | 一种滚珠丝杠进给***摩擦补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410157347.7A CN103926875B (zh) | 2014-04-18 | 2014-04-18 | 一种滚珠丝杠进给***摩擦补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103926875A true CN103926875A (zh) | 2014-07-16 |
CN103926875B CN103926875B (zh) | 2016-08-31 |
Family
ID=51145134
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410157347.7A Active CN103926875B (zh) | 2014-04-18 | 2014-04-18 | 一种滚珠丝杠进给***摩擦补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103926875B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104950806A (zh) * | 2015-06-18 | 2015-09-30 | 南京航空航天大学 | 一种基于gmdh数据挖掘算法的机床进给***前馈控制方法 |
CN105022345A (zh) * | 2015-07-21 | 2015-11-04 | 哈尔滨工业大学 | 基于时间库伦摩擦模型的伺服***摩擦补偿方法 |
CN105700470A (zh) * | 2016-02-01 | 2016-06-22 | 华中科技大学 | 一种用于减小机床伺服进给***跟踪误差的方法 |
CN106596087A (zh) * | 2016-12-12 | 2017-04-26 | 山东大学 | 一种新型可实现全组件摩擦辨识的装置和方法 |
CN106695458A (zh) * | 2017-02-22 | 2017-05-24 | 华中科技大学 | 一种数控机床进给***的摩擦力建模方法 |
CN109520734A (zh) * | 2018-10-28 | 2019-03-26 | 北京工业大学 | 精密减速器几何传动误差最佳测量转速的确定方法 |
CN109556778A (zh) * | 2018-12-06 | 2019-04-02 | 北京工业大学 | 一种基于电流的滚珠丝杠副预紧力在线监测方法 |
CN109660154A (zh) * | 2019-02-01 | 2019-04-19 | 东南大学 | 一种伺服柔性负载新型速度控制方法 |
CN109725659A (zh) * | 2017-10-27 | 2019-05-07 | 北京精密机电控制设备研究所 | 一种面向伺服***负载匹配的球窝喷管负载特性辨识方法 |
CN111788040A (zh) * | 2018-12-28 | 2020-10-16 | 深圳配天智能技术研究院有限公司 | 机器人的动力学参数辨识方法、机器人和存储装置 |
CN112152522A (zh) * | 2020-09-29 | 2020-12-29 | 中国船舶重工集团公司第七二四研究所 | 一种电机全速范围无位置传感器交互多模型控制方法 |
CN112830403A (zh) * | 2021-03-15 | 2021-05-25 | 西安丰树电子科技发展有限公司 | 塔式起重机吊重精确检测的补偿方法 |
CN113110305A (zh) * | 2021-04-30 | 2021-07-13 | 华中科技大学 | 一种机电***的摩擦力建模方法及其应用 |
WO2022047935A1 (zh) * | 2020-09-01 | 2022-03-10 | 瑞声声学科技(深圳)有限公司 | 电机摩擦力补偿方法、***及其计算机可读存储介质 |
CN114765439A (zh) * | 2022-03-16 | 2022-07-19 | 北京建筑大学 | 一种基于摩擦补偿的连续回转马达***简化动态柔性变结构pd复合控制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070088446A1 (en) * | 2005-10-13 | 2007-04-19 | Honeywell International Inc. | Apparatus and method for stiction compensation in a process control system |
CN102023612A (zh) * | 2010-12-01 | 2011-04-20 | 西安交通大学 | 一种数控机床伺服***摩擦误差补偿方法 |
CN102033509A (zh) * | 2011-01-18 | 2011-04-27 | 上海交通大学 | 数控机床误差虚拟补偿*** |
-
2014
- 2014-04-18 CN CN201410157347.7A patent/CN103926875B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070088446A1 (en) * | 2005-10-13 | 2007-04-19 | Honeywell International Inc. | Apparatus and method for stiction compensation in a process control system |
CN102023612A (zh) * | 2010-12-01 | 2011-04-20 | 西安交通大学 | 一种数控机床伺服***摩擦误差补偿方法 |
CN102033509A (zh) * | 2011-01-18 | 2011-04-27 | 上海交通大学 | 数控机床误差虚拟补偿*** |
Non-Patent Citations (2)
Title |
---|
刘红: "基于遗传算法的摩擦模型参数辨识研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》, 15 April 2008 (2008-04-15) * |
向红标: "开放式伺服***的摩擦建模与补偿研究", 《中国博士学位论文全文数据库信息科技辑》, 15 July 2011 (2011-07-15) * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104950806A (zh) * | 2015-06-18 | 2015-09-30 | 南京航空航天大学 | 一种基于gmdh数据挖掘算法的机床进给***前馈控制方法 |
CN104950806B (zh) * | 2015-06-18 | 2017-07-28 | 南京航空航天大学 | 一种基于gmdh数据挖掘算法的机床进给***前馈控制方法 |
CN105022345A (zh) * | 2015-07-21 | 2015-11-04 | 哈尔滨工业大学 | 基于时间库伦摩擦模型的伺服***摩擦补偿方法 |
CN105700470A (zh) * | 2016-02-01 | 2016-06-22 | 华中科技大学 | 一种用于减小机床伺服进给***跟踪误差的方法 |
CN105700470B (zh) * | 2016-02-01 | 2017-12-29 | 华中科技大学 | 一种用于减小机床伺服进给***跟踪误差的方法 |
CN106596087A (zh) * | 2016-12-12 | 2017-04-26 | 山东大学 | 一种新型可实现全组件摩擦辨识的装置和方法 |
CN106695458A (zh) * | 2017-02-22 | 2017-05-24 | 华中科技大学 | 一种数控机床进给***的摩擦力建模方法 |
CN109725659A (zh) * | 2017-10-27 | 2019-05-07 | 北京精密机电控制设备研究所 | 一种面向伺服***负载匹配的球窝喷管负载特性辨识方法 |
CN109725659B (zh) * | 2017-10-27 | 2021-12-17 | 北京精密机电控制设备研究所 | 一种面向伺服***负载匹配的球窝喷管负载特性辨识方法 |
CN109520734A (zh) * | 2018-10-28 | 2019-03-26 | 北京工业大学 | 精密减速器几何传动误差最佳测量转速的确定方法 |
CN109556778A (zh) * | 2018-12-06 | 2019-04-02 | 北京工业大学 | 一种基于电流的滚珠丝杠副预紧力在线监测方法 |
CN111788040B (zh) * | 2018-12-28 | 2023-06-06 | 深圳配天智能技术研究院有限公司 | 机器人的动力学参数辨识方法、机器人和存储装置 |
CN111788040A (zh) * | 2018-12-28 | 2020-10-16 | 深圳配天智能技术研究院有限公司 | 机器人的动力学参数辨识方法、机器人和存储装置 |
CN109660154A (zh) * | 2019-02-01 | 2019-04-19 | 东南大学 | 一种伺服柔性负载新型速度控制方法 |
WO2022047935A1 (zh) * | 2020-09-01 | 2022-03-10 | 瑞声声学科技(深圳)有限公司 | 电机摩擦力补偿方法、***及其计算机可读存储介质 |
CN112152522A (zh) * | 2020-09-29 | 2020-12-29 | 中国船舶重工集团公司第七二四研究所 | 一种电机全速范围无位置传感器交互多模型控制方法 |
CN112152522B (zh) * | 2020-09-29 | 2023-09-26 | 中国船舶集团有限公司第七二四研究所 | 一种电机全速范围无位置传感器交互多模型控制方法 |
CN112830403A (zh) * | 2021-03-15 | 2021-05-25 | 西安丰树电子科技发展有限公司 | 塔式起重机吊重精确检测的补偿方法 |
CN112830403B (zh) * | 2021-03-15 | 2022-11-04 | 西安丰树电子科技发展有限公司 | 塔式起重机吊重精确检测的补偿方法 |
CN113110305A (zh) * | 2021-04-30 | 2021-07-13 | 华中科技大学 | 一种机电***的摩擦力建模方法及其应用 |
CN114765439B (zh) * | 2022-03-16 | 2023-03-10 | 北京建筑大学 | 一种基于摩擦补偿的连续回转马达***简化动态柔性变结构pd复合控制方法 |
CN114765439A (zh) * | 2022-03-16 | 2022-07-19 | 北京建筑大学 | 一种基于摩擦补偿的连续回转马达***简化动态柔性变结构pd复合控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103926875B (zh) | 2016-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103926875A (zh) | 一种滚珠丝杠进给***摩擦补偿方法 | |
Maiti et al. | An adaptive speed sensorless induction motor drive with artificial neural network for stability enhancement | |
CN105915121B (zh) | 一种采用遗传算法优化的伺服***惯量辨识方法 | |
CN109946975A (zh) | 一种未知伺服***的强化学习最优跟踪控制方法 | |
CN102385342B (zh) | 虚拟轴机床并联机构运动控制的自适应动态滑模控制方法 | |
CN102129259B (zh) | 用于砂尘环境试验风洞的基于神经网络pi的智能温度控制***和方法 | |
CN103406909B (zh) | 一种机械臂***的跟踪控制设备及方法 | |
CN105351596A (zh) | 比例电磁阀复合控制方法 | |
CN113379164B (zh) | 基于深度自注意力网络的负荷预测方法及*** | |
Wenjing | Parameter identification of LuGre friction model in servo system based on improved particle swarm optimization algorithm | |
CN103279039A (zh) | 一种机器人神经网络式计算力矩控制器训练平台及训练方法 | |
CN101520857A (zh) | 一种基于神经网络的永磁球形电动机逆运动学求解方法 | |
CN104391497A (zh) | 一种关节伺服***参数辨识和控制参数在线优化方法 | |
CN105159069A (zh) | 一种压电陶瓷作动器的位移控制方法 | |
CN106059418A (zh) | 一种永磁同步直线电机神经网络自适应轨迹跟踪控制方法 | |
CN104038128B (zh) | 基于nurbs的直线电机推力波动控制方法 | |
CN113268040A (zh) | 一种数控机床伺服进给***控制参数的优选方法 | |
Liu et al. | Modeling of hysteresis in piezoelectric actuator based on adaptive filter | |
CN113885321B (zh) | 基于记忆相关Koopman的双模超声电机死区模糊补偿与线性预测控制方法及*** | |
CN112433475B (zh) | 基于scn***偏移量辨识的坐垫机器人限时学习控制方法 | |
CN104022701B (zh) | 一种永磁同步直线电机牛顿法内模速度控制方法 | |
Chu et al. | Classifier-based approximator for friction compensation in high accelerated positioning system | |
Zhang et al. | Adaptive fault-tolerant sliding mode control for high-speed trains with actuator faults under strong winds | |
CN108959787A (zh) | 考虑实际工况的宏宏双驱动***的热变形预测方法及*** | |
CN109884890B (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |