发明名称 机电伺服系统基于极限环振荡的快速摩擦辨识方法
摘要 本发明涉及一种机电伺服系统基于极限环振荡的快速摩擦辨识方法,首先将机电伺服系统的线性模型与参数未知的摩擦模型分别用数学表达式进行描述,然后通过调整两个继电器的增益使系统产生自振荡,测得周期运动的数据,最后利用谐波平衡条件得到的原理公式进行数学运算,计算出摩擦模型的参数。本发明能够快速准确地辨识机电伺服系统的摩擦模型,可以为伺服控制系统的设计提供数学模型,提高机电伺服系统的设计效率。由于这种方法运算量小,不需要精心设计激励信号,方便快捷易于实现,适用于高精度的机电伺服控制系统设计。
申请公布号 CN106505922A 申请公布日期 2017.03.15
申请号 CN201610929493.6 申请日期 2016.10.31
申请人 中国航空工业集团公司洛阳电光设备研究所 发明人 张新勇;冉承平;刘昇;杨光;杨渊强
分类号 H02P21/05(2006.01)I;H02P21/14(2016.01)I 主分类号 H02P21/05(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 一种机电伺服系统基于极限环振荡的快速摩擦辨识方法,其特征在于:设定继电器1的动作电压幅值d<sub>1</sub>和继电器2的动作电压幅值d<sub>2</sub>,选择的动作电压幅值使得伺服系统产生振动频率ω的低频震荡;所述动作电压幅值d<sub>2</sub>使得振动幅值A<sub>v</sub>保证不会超出伺服系统的运动范围;摩擦辨识步骤如下:步骤1:任意改变继电器1的动作电压幅值d<sub>1</sub>,令继电器2的动作电压幅值d<sub>2</sub>不变,进行两次改变d<sub>1</sub>为d<sub>11</sub>和d<sub>12</sub>并测试继电反馈信号,得到静摩擦力矩T<sub>s</sub>的估计值为:<maths num="0001"><math><![CDATA[<mrow><msub><mover><mi>T</mi><mo>^</mo></mover><mi>s</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>&pi;A</mi><mn>1</mn></msub><msub><mi>A</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>R</mi><mrow><mi>e</mi><mi>q</mi><mn>2</mn></mrow></msub><msub><mi>R</mi><mrow><mi>g</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>R</mi><mrow><mi>e</mi><mi>q</mi><mn>1</mn></mrow></msub><msub><mi>R</mi><mrow><mi>g</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow><mo>+</mo><mn>4</mn><mrow><mo>(</mo><msub><mi>d</mi><mn>12</mn></msub><msub><mi>A</mi><mn>1</mn></msub><msub><mi>R</mi><mrow><mi>g</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>d</mi><mn>11</mn></msub><msub><mi>A</mi><mn>2</mn></msub><msub><mi>R</mi><mrow><mi>g</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow></mrow><mrow><mn>4</mn><msub><mi>R</mi><mrow><mi>g</mi><mn>1</mn></mrow></msub><msub><mi>R</mi><mrow><mi>g</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><msub><mi>M</mi><mn>1</mn></msub><msub><mi>A</mi><mn>2</mn></msub><mo>-</mo><msub><mi>M</mi><mn>2</mn></msub><msub><mi>A</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0001136949500000011.GIF" wi="1014" he="135" /></maths>同时还可以得到在低速模式下的等效阻尼系数f<sub>0</sub>的估计值为:<maths num="0002"><math><![CDATA[<mrow><msub><mover><mi>f</mi><mo>^</mo></mover><mn>0</mn></msub><mo>=</mo><mfrac><mrow><mn>4</mn><msub><mi>d</mi><mn>11</mn></msub></mrow><mrow><msub><mi>&pi;A</mi><mn>1</mn></msub><msub><mi>R</mi><mrow><mi>g</mi><mn>1</mn></mrow></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>4</mn><msub><mover><mi>T</mi><mo>^</mo></mover><mi>s</mi></msub><msub><mi>M</mi><mn>1</mn></msub></mrow><mrow><msub><mi>&pi;A</mi><mn>1</mn></msub></mrow></mfrac><mo>+</mo><mfrac><msub><mi>R</mi><mrow><mi>e</mi><mi>q</mi><mn>1</mn></mrow></msub><msub><mi>R</mi><mrow><mi>g</mi><mn>1</mn></mrow></msub></mfrac></mrow>]]></math><img file="FDA0001136949500000012.GIF" wi="497" he="135" /></maths>f<sub>0</sub>的表达式为f<sub>0</sub>=(T<sub>s</sub>‑T<sub>c</sub>)/v<sub>s</sub>其中,A<sub>1</sub>、A<sub>2</sub>、ω<sub>1</sub>和ω<sub>2</sub>分别为两次实验中产生周期振动的幅值和频率,相对应的d<sub>11</sub>和d<sub>12</sub>分别为两次实验中所设定的继电器1的增益。式中其他参数具体为:<img file="FDA0001136949500000013.GIF" wi="435" he="79" />为传感器2频响函数的实部;<img file="FDA0001136949500000014.GIF" wi="454" he="80" />为机电伺服控制系统线性单元传递函数的实部;<img file="FDA0001136949500000015.GIF" wi="301" he="93" />为传感器传递函数的模,其中i=1,2;步骤2:保持继电器1的动作电压幅值d<sub>1</sub>,增大继电器2的动作电压幅值d<sub>2</sub>产生幅值A<sub>v</sub>>v<sub>s</sub>的周期振动,再任意改变继电器1的动作电压幅值为d<sub>13</sub>和d<sub>14</sub>,进行两次改变并测试继电反馈信号,得到粘性阻尼系数v<sub>s</sub>和等效参数T<sub>0</sub>的估计值为:<maths num="0003"><math><![CDATA[<mrow><msub><mover><mi>T</mi><mo>^</mo></mover><mn>0</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>&pi;A</mi><mn>1</mn></msub><msub><mi>A</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>R</mi><mrow><mi>e</mi><mi>q</mi><mn>2</mn></mrow></msub><msub><mi>R</mi><mrow><mi>g</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>R</mi><mrow><mi>e</mi><mi>q</mi><mn>1</mn></mrow></msub><msub><mi>R</mi><mrow><mi>g</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow><mo>+</mo><mn>4</mn><mrow><mo>(</mo><msub><mi>d</mi><mn>12</mn></msub><msub><mi>A</mi><mn>1</mn></msub><msub><mi>R</mi><mrow><mi>g</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>d</mi><mn>11</mn></msub><msub><mi>A</mi><mn>2</mn></msub><msub><mi>R</mi><mrow><mi>g</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow></mrow><mrow><mn>4</mn><msub><mi>R</mi><mrow><mi>g</mi><mn>1</mn></mrow></msub><msub><mi>R</mi><mrow><mi>g</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><msub><mi>M</mi><mn>1</mn></msub><msub><mi>A</mi><mn>2</mn></msub><mo>-</mo><msub><mi>M</mi><mn>2</mn></msub><msub><mi>A</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0001136949500000016.GIF" wi="1028" he="135" /></maths><maths num="0004"><math><![CDATA[<mrow><msub><mover><mi>f</mi><mo>^</mo></mover><mi>v</mi></msub><mo>=</mo><mo>-</mo><mfrac><mrow><mn>4</mn><msub><mi>d</mi><mn>11</mn></msub></mrow><mrow><msub><mi>&pi;A</mi><mn>1</mn></msub><msub><mi>R</mi><mrow><mi>g</mi><mn>1</mn></mrow></msub></mrow></mfrac><mo>-</mo><mfrac><mrow><mn>4</mn><msub><mover><mi>T</mi><mo>^</mo></mover><mn>0</mn></msub><msub><mi>M</mi><mn>1</mn></msub></mrow><mrow><msub><mi>&pi;A</mi><mn>1</mn></msub></mrow></mfrac><mo>-</mo><mfrac><msub><mi>R</mi><mrow><mi>e</mi><mi>q</mi><mn>1</mn></mrow></msub><msub><mi>R</mi><mrow><mi>g</mi><mn>1</mn></mrow></msub></mfrac></mrow>]]></math><img file="FDA0001136949500000017.GIF" wi="542" he="135" /></maths>通过谐波平衡条件改写得来库仑摩擦力矩T<sub>c</sub>及边界润滑速度v<sub>s</sub>,求解<maths num="0005"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>J</mi><mrow><mo>(</mo><msub><mover><mi>v</mi><mo>^</mo></mover><mi>s</mi></msub><mo>,</mo><msub><mover><mi>T</mi><mo>^</mo></mover><mi>c</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>{</mo><msup><mrow><mo>&lsqb;</mo><mi>Re</mi><mrow><mo>(</mo><msub><mi>N</mi><mrow><mi>e</mi><mi>q</mi></mrow></msub><mo>(</mo><mrow><msub><mi>A</mi><mi>i</mi></msub><mo>,</mo><msub><mi>&omega;</mi><mi>i</mi></msub><mo>,</mo><msub><mover><mi>v</mi><mo>^</mo></mover><mi>s</mi></msub><mo>,</mo><msub><mover><mi>T</mi><mo>^</mo></mover><mi>c</mi></msub></mrow><mo>)</mo><mo>)</mo></mrow><mo>+</mo><mi>Re</mi><mrow><mo>(</mo><mn>1</mn><mo>/</mo><msub><mover><mi>G</mi><mo>^</mo></mover><mrow><mi>e</mi><mi>q</mi></mrow></msub><mo>(</mo><mrow><msub><mi>j&omega;</mi><mi>i</mi></msub></mrow><mo>)</mo><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><msup><mrow><mo>&lsqb;</mo><mi>Im</mi><mrow><mo>(</mo><msub><mi>N</mi><mrow><mi>e</mi><mi>q</mi></mrow></msub><mo>(</mo><mrow><msub><mi>A</mi><mi>i</mi></msub><mo>,</mo><msub><mi>&omega;</mi><mi>i</mi></msub><mo>,</mo><msub><mover><mi>v</mi><mo>^</mo></mover><mi>s</mi></msub><mo>,</mo><msub><mover><mi>T</mi><mo>^</mo></mover><mi>c</mi></msub></mrow><mo>)</mo><mo>)</mo></mrow><mo>+</mo><mi>Im</mi><mrow><mo>(</mo><mn>1</mn><mo>/</mo><msub><mover><mi>G</mi><mo>^</mo></mover><mrow><mi>e</mi><mi>q</mi></mrow></msub><mo>(</mo><mrow><msub><mi>j&omega;</mi><mi>i</mi></msub></mrow><mo>)</mo><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup><mo>}</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001136949500000021.GIF" wi="1149" he="222" /></maths>在满足<img file="FDA0001136949500000022.GIF" wi="395" he="71" />和<img file="FDA0001136949500000023.GIF" wi="221" he="71" />时,得到库仑摩擦力矩T<sub>c</sub>及边界润滑速度v<sub>s</sub>。
地址 471099 河南省洛阳市凯旋西路25号