CN103019092A - 一种机械传动***定位平台的预测控制方法 - Google Patents
一种机械传动***定位平台的预测控制方法 Download PDFInfo
- Publication number
- CN103019092A CN103019092A CN2012105942984A CN201210594298A CN103019092A CN 103019092 A CN103019092 A CN 103019092A CN 2012105942984 A CN2012105942984 A CN 2012105942984A CN 201210594298 A CN201210594298 A CN 201210594298A CN 103019092 A CN103019092 A CN 103019092A
- Authority
- CN
- China
- Prior art keywords
- model
- sigma
- mechanical transmission
- gradient
- output
- 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
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种机械传动***定位平台的预测控制方法,该机械传动***带有滚珠丝杠,所述的机械传动***定位平台包括的饲服电机的输出转角作为***的输入信号u(k),经过齿轮箱变速及丝杠的机械传动,再带动工作台负载工作,工作台的位移y(k),即为整个***的输出,给出控制误差ε,根据***模型,***的工作区间4→1→2和2→3→4都能估计出来,在第k步,能计算出在非光滑点,J(k)关于u(k)的Clarke次梯度,i.e.,其中j∈Jk,Jk={1,2,…,t},|Jk|是Jk的元素的个数,|Jk|≤t1给定的有界的自然数,t=t+1,如果t≤t1,那么Jk={1,…,t};如果t>t1,那么Jk=Jk-1∪{t}\{t-t1};或者在光滑点,J(k)关于u(k)的梯度 其中我们令或者那么我们就能求出***的准输入u1(k)。
Description
技术领域
本发明属于机械传动控制技术领域,特别涉及一种机械传动***定位平台的预测控制方法。
背景技术
许多常见的精密加工***所需的设备,譬如高精度万能铣床、三坐标测量仪等,往往采用含有滚珠丝杠装置的机械传动机构的工作平台,这些机械传动机构由于配合或磨损造成的死区、间隙等非线性特性,而这些特性对***的控制性能影响是不能忽略的,它会降低***的控制精度,有时会引起***的抖动,甚至会出现***的不稳定。因此,现存的不少技术大都是对这些非线性特性进行补偿,由于这些非线性特性建模的时候会有比较大的误差,并且随着设备的磨损、老化,模型的误差会越来越大,因此会存在着较大的补偿误差,这在超精密加工设备中体现得更为明显,能否找到一种不必对这些非线性特性进行直接补偿,同时又能消除其影响,获得满意控制性能的实时控制方法,成为提高超精密***控制精度和性能的关键和难点。
发明内容
本发明提出一种不必对这类***中滚珠丝杠产生的间隙这一类非光滑的非线性特性进行直接补偿、却能消除间隙影响的非光滑预测控制方法。
本发明的技术方案是,一种机械传动***定位平台的预测控制方法,该机械传动***带有滚珠丝杠,所述的机械传动***定位平台包括的饲服电机的输出转角作为***的输入信号u(k),经过齿轮箱变速及丝杠的机械传动,再带动工作台负载工作,工作台的位移y(k),即为整个***的输出,其中,
所述定位平台的间隙用来描述齿轮箱与丝杠的机械传动的非线性特性,
L1(·)是一个线性的动态子模型,用来描述工作平台的负载,
机械传动***的非线性特性的输出,即x(k),且不能直接测量,
所述定位平台的输入和输出为u(k)和y(k)
所述的间隙的特性描述为,
L1(·)模型可以表述为:
由式(1)和(2)组成机械传动***的模型,则有线性子***表示成
其中ε(k)为模型误差,设其为零均值的白噪声,根据Diophantine方程
其中 和 相应的d步超前预测模型为
其中 和
采用式(5)来作为预测控制的目标函数,
其中,r(k+d)是参考轨迹,是预测模型的输出,λ是一个非负的加权系数,J(k)是Lipschitz连续的,在该优化函数的非光滑点,用次微分来替代传统意义的梯度,J(k)是伪凸的,在非光滑点存在唯一的最优值,在该点的广义梯度为0,
在间隙的非光滑点,J(k)关于u(k)的Clarke次梯度可以描述为:
则所述的预测控制方法包括步骤,
第二步,根据***模型,***的工作区间4→1→2和2→3→4都能估计出来,在第k步,我们能计算出,在非光滑点,J(k)关于u(k)的Clarke次梯度,i.e.,其中j∈Jk,Jk={1,2,…,t},|Jk|是Jk的元素的个数,|Jk|≤t1给定的有界的自然数,t=t+1,如果t≤t1,那么Jk={1,…,t};如果t>t1,那么Jk=Jk-1∪{t}\{t-t1};或者在光滑点,J(k)关于u(k)的梯度
第三步, 其中 我们令 或者 那么我们就能求出***的准输入u1(k);
第六步,k=k+1,转入第二步。
本发明是针对带有滚珠丝杠一类机械传动装置的定位平台的一种预测控制方法。在机械传动中滚珠丝杠由于配合及磨损所产生的间隙现象不能直接测量,如果不对间隙的影响进行直接补偿的情况下,用一种非光滑的预测控制策略消除间隙对机械传动的不利影响。用非光滑的优化技术来设计该非光滑***的预测控制器,解决了非光滑点的优化预测控制问题,不必对定位平台中滚珠丝杠产生的间隙进行直接补偿,避免了其他控制方法中需要对间隙求逆的困难。
附图说明
图1本发明的含有滚珠丝杠的机械传动机构工作平台的实际物理结构示意图。
图2本发明的采用滚动丝杠的机械传动机构定位平台模型的结构图。
具体实施方式
含有机械传动机构工作平台的预测模型,例如,含有滚珠丝杠的机械传动机构工作平台的实际物理结构,可以用如图1所示,
其中饲服电机的输出转角作为***的输入信号u(k),经过齿轮箱变速及丝杠的机械传动,再带动工作台(负载)工作,工作台的位移y(k),即为整个***的输出。
图1所示的***可以进一步用图2所示的模型结构图表示,其中间隙用来描述齿轮箱与丝杠的机械传动的非线性特性,L1(·)是一个线性的动态子模型,用来描述工作平台的负载,机械传动的非线性特性的输出,即x(k)不能直接测量,只有整个***的输入和输出,即:u(k)和y(k)可以直接测量。图2中,间隙的特性可以描述为:
L1(·)模型可以表述为:
由此,式(1)和(2)就组成了机械传动机构的模型。则线性子***可以表示成
其中ε(k)为模型误差,设其为零均值的白噪声。根据Diophantine方程
其中 和 相应的d步超前预测模型为
其中 和
2.采用丝杠机械传动机构工作平台的预测控制
我们用式(5)来作为预测控制的目标函数。
J(k)是Lipschitz连续的,因此,在该优化函数的非光滑点,我们可以用次微分来替代传统意义的梯度,又因为J(k)是伪凸的,所以在非光滑点存在唯一的最优值,并且在该点的广义梯度为0.
因此,在间隙的非光滑点,J(k)关于u(k)的Clarke次梯度可以描述为:
该预测的算法为:
第二步,根据***模型,***的工作区间4→1→2和2→3→4都能估计出来,在第k步,我们能计算出,在非光滑点,J(k)关于u(k)的Clarke次梯度,i.e.,其中j∈Jk,Jk={1,2,…,t},|Jk|是Jk的元素的个数,|Jk|≤t1给定的有界的自然数,t=t+1,如果t≤t1,那么Jk={1,…,t};如果t>t1,那么Jk=Jk-1∪{t}\{t-t1};或者在光滑点,J(k)关于u(k)的梯度
第三步, 其中 我们令 或者 那么我们就能求出***的准输入u1(k)。
第六步,k=k+1,转入第二步。
可以采用角编码器作为传感器,对其测量值进行解码和差分后获得含有机械平台的速度,并传送到数字信号处理器(DSP)中。DSP计算控制算法输出。经过放大器放大后驱动饲服电机。
Claims (1)
1.一种机械传动***定位平台的预测控制方法,该机械传动***带有滚珠丝杠,其特征在于,所述的机械传动***定位平台包括的饲服电机的输出转角作为***的输入信号u(k),经过齿轮箱变速及丝杠的机械传动,再带动工作台负载工作,工作台的位移y(k),即为整个***的输出,其中,
所述定位平台的间隙用来描述齿轮箱与丝杠的机械传动的非线性特性,
L1(·)是一个线性的动态子模型,用来描述工作平台的负载,
机械传动***的非线性特性的输出,即x(k),且不能直接测量,
所述定位平台的输入和输出为u(k)和y(k)
所述的间隙的特性描述为,
L1(·)模型表述为:
由式(1)和(2)组成机械传动***的模型,则有线性子***表示成
其中ε(k)为模型误差,设其为零均值的白噪声,根据Diophantine方程
其中 和 相应的d步超前预测模型为
其中 和
采用式(5)来作为预测控制的目标函数,
其中,r(k+d)是参考轨迹,是预测模型的输出,λ是一个非负的加权系数,
J(k)是Lipschitz连续的,在该优化函数的非光滑点,用次微分来替代传统意义的梯度,J(k)是伪凸的,在非光滑点存在唯一的最优值,在该点的广义梯度为0,
在间隙的非光滑点,J(k)关于u(k)的Clarke次梯度可以描述为:
则所述的预测控制方法包括步骤,
第二步,根据***模型,***的工作区间4→1→2和2→3→4都能估计出来,在第k步,我们能计算出,在非光滑点,J(k)关于u(k)的Clarke次梯度,i.e.,其中j∈Jk,Jk={1,2,…,t},|Jk|是Jk的元素的个数,|Jk|≤t1给定的有界的自然数,t t+1,如果t≤t1,那么Jk={1,…,t};如果t>t1,那么Jk=Jk-1∪{t}\{t-t1};或者在光滑点,J(k)关于u(k)的梯度
第三步, 其中 我们令 或者 那么我们就能求出***的准输入u1(k);
第六步,k=k+1,转入第二步。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210594298.4A CN103019092B (zh) | 2012-12-31 | 2012-12-31 | 一种机械传动***定位平台的预测控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210594298.4A CN103019092B (zh) | 2012-12-31 | 2012-12-31 | 一种机械传动***定位平台的预测控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103019092A true CN103019092A (zh) | 2013-04-03 |
CN103019092B CN103019092B (zh) | 2016-07-06 |
Family
ID=47967812
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210594298.4A Expired - Fee Related CN103019092B (zh) | 2012-12-31 | 2012-12-31 | 一种机械传动***定位平台的预测控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103019092B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103558760A (zh) * | 2013-11-04 | 2014-02-05 | 上海师范大学 | 一种称重式塑料原料配料机送料过程的预测控制方法 |
CN105372074A (zh) * | 2015-10-29 | 2016-03-02 | 重庆长安汽车股份有限公司 | 一种发动机非线性阶次分量提取分析方法 |
CN111208818A (zh) * | 2020-01-07 | 2020-05-29 | 电子科技大学 | 基于视觉时空特征的智能车预测控制方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101276207A (zh) * | 2008-05-07 | 2008-10-01 | 浙江大学 | 基于哈默斯坦模型的多变量非线性***预测函数控制方法 |
US20120259600A1 (en) * | 2011-04-11 | 2012-10-11 | King Fahd University Of Petroleum And Minerals | Method of identifying hammerstein models with known nonlinearity structures using particle swarm optimization |
-
2012
- 2012-12-31 CN CN201210594298.4A patent/CN103019092B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101276207A (zh) * | 2008-05-07 | 2008-10-01 | 浙江大学 | 基于哈默斯坦模型的多变量非线性***预测函数控制方法 |
US20120259600A1 (en) * | 2011-04-11 | 2012-10-11 | King Fahd University Of Petroleum And Minerals | Method of identifying hammerstein models with known nonlinearity structures using particle swarm optimization |
Non-Patent Citations (3)
Title |
---|
孙西等: "《具有间隙非线性***的自适应控制》", 《控制理论与应用》, vol. 8, no. 1, 30 March 1991 (1991-03-30) * |
李冬伍等: "《含间隙非线性的模块化***辨识及鲁棒逆补偿控制》", 《第三十一届中国控制会议论文集A卷》, 27 July 2012 (2012-07-27) * |
董瑞丽: "《非光滑三明治***的辨识和控制研究》", 《上海交通大学博士论文2009》, 31 December 2009 (2009-12-31) * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103558760A (zh) * | 2013-11-04 | 2014-02-05 | 上海师范大学 | 一种称重式塑料原料配料机送料过程的预测控制方法 |
CN105372074A (zh) * | 2015-10-29 | 2016-03-02 | 重庆长安汽车股份有限公司 | 一种发动机非线性阶次分量提取分析方法 |
CN105372074B (zh) * | 2015-10-29 | 2017-11-07 | 重庆长安汽车股份有限公司 | 一种发动机非线性阶次分量提取分析方法 |
CN111208818A (zh) * | 2020-01-07 | 2020-05-29 | 电子科技大学 | 基于视觉时空特征的智能车预测控制方法 |
CN111208818B (zh) * | 2020-01-07 | 2023-03-07 | 电子科技大学 | 基于视觉时空特征的智能车预测控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103019092B (zh) | 2016-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104075890B (zh) | 一种伺服电机及谐波减速机综合测试平台 | |
US8731886B2 (en) | Simulator for estimating life of robot speed reducer | |
CN103472730B (zh) | 一种基于谐波减速器扭转刚度迟滞模型的双框架控制力矩陀螺高精度框架速率伺服*** | |
CN108983703B (zh) | 超精密运动***前馈控制器参数整定方法 | |
CN103753351B (zh) | 一种全闭环伺服***的数控机床反向间隙监测方法 | |
CN103092075A (zh) | 一种含有间隙特性的机械传动***模型参数确定方法 | |
CN102430683B (zh) | 精密伺服锻压机的误差实时在线补偿方法 | |
CN103019092A (zh) | 一种机械传动***定位平台的预测控制方法 | |
CN105404237A (zh) | 一种基于空间网格补偿方式的数控机床空间误差建模方法 | |
CN101941304B (zh) | 压力机械控制装置 | |
CN105149645A (zh) | 一种用于机器人制孔***锪窝深度控制的装置及制孔方法 | |
JP5657633B2 (ja) | 移動体が反転するときの位置誤差を補正するサーボ制御装置 | |
CN102637014A (zh) | 数控机床加工过程机电主传动***能量效率获取的方法 | |
CN105305913A (zh) | 一种新型的用于滚珠丝杠进给***的抗扰跟随控制器 | |
CN108020415A (zh) | 一种谐波减速器静态扭转刚度的测量装置 | |
JP5453606B2 (ja) | 波動歯車装置を備えたアクチュエータの位置決め制御方法 | |
CN101751002A (zh) | 一种用于大型数控机床的温度补偿***及补偿方法 | |
US20150045940A1 (en) | Servo control device and servo control method | |
CN110077028A (zh) | 伺服压力机全闭环非线性预测控制方法与*** | |
CN110398359A (zh) | 一种机械传动链综合误差的动态检测方法及装置 | |
CN108241289B (zh) | 一种用于三自由度加载***的力-位移混合控制方法 | |
CN101260967A (zh) | Xyy精密定位平台的校准方法 | |
CN101629801A (zh) | 数控磨床导轨热误差确定方法 | |
CN102658503A (zh) | 一种基于内置传感器的数控机床进给***的模态测试方法 | |
Yong et al. | Dynamic characteristic optimization of ball screw feed drive in machine tool based on modal extraction of state space model |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160706 Termination date: 20181231 |
|
CF01 | Termination of patent right due to non-payment of annual fee |