发明名称 基于坐标变换和参数调整直接驱动数控平台鲁棒控制方法
摘要 一种基于坐标变换和参数调整直接驱动数控平台鲁棒控制方法,首先确定电机动子的初始相位;然后把轮廓误差进行坐标变换,并参考调整,得到位置误差,即进行位置偏差计算,判断是否进行位置调节;最后执行鲁棒控制算法,输出控制量,驱动数控平台。本发明方法采用的控制系统包括电压调整电路、整流滤波单元、IPM逆变单元、DSP、霍尔传感器、光栅尺、电流采样电路、位置采样电路、IPM隔离驱动保护电路。本发明针对直接驱动数控平台,提出了一个鲁棒轮廓控制器。基于一个完整的坐标系转换和参数调节函数,能被应用到任何光滑的轮廓曲线,对于设备建模误差和干扰提供鲁棒控制系统稳定性。
申请公布号 CN102637011A 申请公布日期 2012.08.15
申请号 CN201110390966.7 申请日期 2011.11.30
申请人 沈阳工业大学 发明人 王丽梅;郑浩;赵希梅;孙宜标;刘春芳
分类号 G05B19/18(2006.01)I 主分类号 G05B19/18(2006.01)I
代理机构 沈阳东大专利代理有限公司 21109 代理人 梁焱
主权项 1.一种基于坐标变换和参数调整直接驱动数控平台鲁棒控制方法,其特征在于:具体步骤如下:步骤1:确定电机动子的初始相位;通过动子位置采样电路和电流采样电路分别采集电机动子位置、速度和电流信息;步骤2:对轮廓误差进行坐标变换,并进行参数调整,得到位置误差,即进行位置偏差计算,判断是否进行位置调节,是则进行步骤3,否则进行电流调节;步骤2.1:对轮廓误差进行坐标变换假设:(a1)理想轨迹r<sub>i</sub>和它的导数<img file="FDA0000114480170000011.GIF" wi="39" he="59" />和<img file="FDA0000114480170000012.GIF" wi="36" he="51" />可得;(a2)局部坐标系∑<sub>l</sub>的倾角θ和它的导数<img file="FDA0000114480170000013.GIF" wi="39" he="59" />和<img file="FDA0000114480170000014.GIF" wi="40" he="51" />可得;(a3)进给驱动系统位置x和它的导数是可测量的;(a4)对于2设备参数值,都是可得的;根据上述假设确定如下控制,包括驱动轴电机的输入电压v、矩阵H、E和I<sub>e</sub>:<img file="FDA0000114480170000015.GIF" wi="1151" he="69" /><maths num="0001"><![CDATA[<math><mrow><mi>H</mi><mo>=</mo><mi>diag</mi><mo>{</mo><mfrac><mrow><msub><mi>m</mi><mi>i</mi></msub><mo>+</mo><msub><mi>n</mi><mi>i</mi></msub><msup><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><msub><mi>p</mi><mi>i</mi></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><msub><mi>p</mi><mi>i</mi></msub></mfrac></mfrac><mo>}</mo><mo>,</mo></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mi>E</mi><mo>=</mo><mi>diag</mi><mo>{</mo><mfrac><mrow><msub><mi>c</mi><mi>i</mi></msub><mo>+</mo><msub><mi>d</mi><mi>i</mi></msub><msup><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><msub><mi>p</mi><mi>i</mi></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow><mfrac><mrow><mn>2</mn><mi>&pi;</mi><msub><mi>k</mi><mi>i</mi></msub></mrow><msub><mi>p</mi><mi>i</mi></msub></mfrac></mfrac><mo>}</mo><mo>,</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mi>I</mi><mi>e</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mo>-</mo><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中m<sub>i</sub>(>0),c<sub>i</sub>(≥0)和f<sub>i</sub>分别是驱动轴i表质量、粘性摩擦系数和驱动力,其中K<sub>v</sub>和K<sub>p</sub>是速度和位置反馈增益矩阵;它们被假设只有正元素的对角矩阵,下面的等式成立:<img file="FDA0000114480170000019.GIF" wi="1596" he="69" />由于矩阵H和R是非奇异的,下面的关系对于实现等式(2)需要被满足:<img file="FDA00001144801700000110.GIF" wi="1616" he="69" />从上面的等式中,通过适当分配等式(3)中反馈增益矩阵K<sub>v</sub>和K<sub>p</sub>得出结论,当t趋近于无穷时,e<sub>l</sub>趋近于0是可得的,如果对于e<sub>l2</sub>反馈增益被设置大于e<sub>l1</sub>的,对于理想曲线,轮廓误差能被减小得比切线跟踪误差快一些,n<sub>i</sub>为驱动进给驱动系统的电机惯性;k<sub>i</sub>为转矩-电压转换率;p<sub>i</sub>是滚珠丝杠的节;k<sub>i</sub>为转矩-电压转换率;步骤2.2:参数调整,具体如下:移动理想位置和局部坐标系,来减少正交误差和实际轮廓误差的不同,因为一个跟踪误差沿着l<sub>1</sub>存在,在误差e<sub>l2</sub>和e<sub>c</sub>之间可能存在一个不可接受的差异,为减少这个差异在理想轮廓曲线到实际位置x估计最近的位置,用r<sub>a</sub>表示;一个局部坐标系被产生,它的方向与最近位置的坐标系相似,它的原点接近理想轮廓曲线的理想位置,r<sub>n</sub>和∑<sub>n</sub>;用新定义系∑<sub>n</sub>取代等式(1)中旋转矩阵R和跟踪误差e<sub>w</sub>和e<sub>l</sub>,产生e<sub>l2</sub>和e<sub>c</sub>的不同,下面说明产生系∑<sub>n</sub>的方法:假设沿着l<sub>1</sub>的跟踪误差有一个负值,因为为减少e<sub>l1</sub>分配的控制器增益正常比为e<sub>l2</sub>分配的小一些,也假设理想位置r和沿着理想轮廓曲线最近点r<sub>a</sub>之间的距离与跟踪误差e<sub>l1</sub>的大小近似等值,此外,沿着这段的理想速度接近常量,于是,要求通过这段的时间t<sub>d</sub>能被如下估计:<img file="FDA0000114480170000021.GIF" wi="1604" he="143" />坐标系∑<sub>a</sub>中原点r<sub>a</sub>和倾角θ<sub>a</sub>,是在理想轮廓曲线到x的最近位置,能被如下估计:r<sub>a</sub>=r(t-t<sub>d</sub>),θ<sub>a</sub>=θ(t-t<sub>d</sub>)                                        (5)其中r()和θ()表明时间的函数;修改的理想位置r<sub>n</sub>表示成:r<sub>n</sub>=r+R<sub>a</sub>d<sub>r</sub>,d<sub>r</sub>=[0,-d<sub>ra2</sub>]<sup>T</sup>                                       (6)其中R<sub>a</sub>是∑<sub>a</sub>到∑<sub>w</sub>的旋转矩阵,是θ<sub>a</sub>的一个函数,控制输入等式(1)被下式代替:<img file="FDA0000114480170000022.GIF" wi="1475" he="83" />其中e<sub>wn</sub>=x-r<sub>n</sub>;于是,获得下面的控制动力,e<sub>n1</sub>和e<sub>n2</sub>收敛速度能被独立地调整:<img file="FDA0000114480170000023.GIF" wi="449" he="73" />步骤3:执行鲁棒控制算法,输出控制量v,驱动数控平台;实际加工中,包括非线性摩擦和切割力以及设备建模错误的干扰的确存在,因此,控制器等式应该扩展考虑干扰和设备建模误差,含有有界干扰矢量w的进给驱动动力学,<img file="FDA0000114480170000031.GIF" wi="1622" he="55" />提出下面的控制器:<img file="FDA0000114480170000032.GIF" wi="1367" he="69" /><maths num="0004"><![CDATA[<math><mrow><mover><mi>H</mi><mo>^</mo></mover><mo>=</mo><mi>diag</mi><mo>{</mo><mfrac><mrow><msub><mover><mi>m</mi><mo>^</mo></mover><mi>i</mi></msub><mo>+</mo><msub><mover><mi>n</mi><mo>^</mo></mover><mi>i</mi></msub><msup><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><msub><mover><mi>p</mi><mo>^</mo></mover><mi>i</mi></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow><mfrac><mrow><mn>2</mn><mi>&pi;</mi><msub><mover><mi>k</mi><mo>^</mo></mover><mi>i</mi></msub></mrow><msub><mover><mi>p</mi><mo>^</mo></mover><mi>i</mi></msub></mfrac></mfrac><mo>}</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><mover><mi>E</mi><mo>^</mo></mover><mo>=</mo><mi>diag</mi><mo>{</mo><mfrac><mrow><msub><mover><mi>c</mi><mo>^</mo></mover><mi>i</mi></msub><mo>+</mo><msub><mover><mi>d</mi><mo>^</mo></mover><mi>i</mi></msub><msup><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><msub><mover><mi>p</mi><mo>^</mo></mover><mi>i</mi></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow><mfrac><mrow><mn>2</mn><mi>&pi;</mi><msub><mover><mi>k</mi><mo>^</mo></mover><mi>i</mi></msub></mrow><msub><mover><mi>p</mi><mo>^</mo></mover><mi>i</mi></msub></mfrac></mfrac><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中<img file="FDA0000114480170000035.GIF" wi="22" he="39" />表示z的标称值,信号<img file="FDA0000114480170000036.GIF" wi="287" he="65" />式中Λ是一个对角矩阵和只有正元素的常量,信号e<sub>v</sub>被采用作为后续稳定性分析的代替<img file="FDA0000114480170000037.GIF" wi="44" he="51" />的速度信号,符号δv是实现鲁棒控制的一个输入矢量;通过近似分配等式反馈增益矩阵K<sub>v</sub>、K<sub>p</sub>和Λ,当t趋近于无穷时能得到e<sub>n</sub>趋近于零;此外,沿着∑<sub>n</sub>每个轴的误差收敛速度能被独立调整;鲁棒控制得到如下等式:<maths num="0006"><![CDATA[<math><mrow><mi>&delta;v</mi><mo>=</mo><mo>-</mo><mfrac><mrow><mi>&eta;</mi><msub><mi>R</mi><mi>a</mi></msub><msub><mi>e</mi><mi>v</mi></msub></mrow><mi>&rho;</mi></mfrac><mo>,</mo></mrow></math>]]></maths><img file="FDA0000114480170000039.GIF" wi="1276" he="152" />其中||a||是a的欧几里得范数,ε是小的正常数;为了表明提出系统的鲁棒稳定性,利亚普诺夫函数候选被应用;得出结论e<sub>v</sub>和e<sub>n</sub>是一致最终有界,于是,从<img file="FDA00001144801700000310.GIF" wi="287" he="65" /><img file="FDA00001144801700000311.GIF" wi="44" he="51" />也是一致最终有界。
地址 110870 辽宁省沈阳市经济技术开发区沈辽西路111号