发明名称 一种双框架磁悬浮CMG框架系统的高精度摩擦补偿控制方法
摘要 本发明公开了一种双框架磁悬浮CMG框架系统的高精度摩擦补偿控制方法,首先建立双框架伺服系统动力学模型,对框架系统的耦合力矩进行前馈补偿控制,通过时滞估计器(TDE)估计框架系统的非线性摩擦力矩,然后采用TDC进行非线性摩擦力矩的反馈补偿,引入两自由度IMC以实现框架系统跟踪性能和扰动抑制能力的独立控制,同时增强TDC的鲁棒性。本发明增强了框架系统对非线性摩擦力矩的鲁棒性,适用于双框架磁悬浮控制力矩陀螺框架伺服系统的高精度摩擦补偿控制。
申请公布号 CN105159083A 申请公布日期 2015.12.16
申请号 CN201510561163.1 申请日期 2015.09.06
申请人 北京航空航天大学 发明人 崔培玲;杨珊;李海涛;房建成;闫斌;宁欣
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 北京科迪生专利代理有限责任公司 11251 代理人 杨学明;顾炜
主权项 一种双框架磁悬浮CMG框架系统的高精度摩擦补偿控制方法,其特征在于包括以下步骤:步骤(1)、建立双框架伺服系统的动力学模型根据几何约束关系,应用欧拉动力学方程及矢量叠加原理,推导双框架磁悬浮控制力矩陀螺内、外框架系统的动力学模型如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mrow><msub><mi>T</mi><mrow><mi>g</mi><mi>x</mi></mrow></msub><mo>=</mo><msub><mi>K</mi><mrow><mi>i</mi><mi>g</mi><mi>x</mi></mrow></msub><msub><mi>i</mi><mrow><mi>g</mi><mi>x</mi></mrow></msub><mo>=</mo><msub><mi>J</mi><mrow><mi>g</mi><mi>x</mi></mrow></msub><msub><mover><mi>&theta;</mi><mo>&CenterDot;&CenterDot;</mo></mover><mi>g</mi></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>J</mi><mrow><mi>g</mi><mi>y</mi></mrow></msub><mo>-</mo><msub><mi>J</mi><mrow><mi>g</mi><mi>z</mi></mrow></msub><mo>)</mo></mrow><msubsup><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>j</mi><mn>2</mn></msubsup><msub><mi>sin&theta;</mi><mi>g</mi></msub><msub><mi>cos&theta;</mi><mi>g</mi></msub><mo>+</mo><msub><mi>J</mi><mrow><mi>r</mi><mi>r</mi></mrow></msub><mo>&lsqb;</mo><mfrac><msqrt><mn>2</mn></msqrt><mn>2</mn></mfrac><mrow><mo>(</mo><mover><mi>&alpha;</mi><mo>&CenterDot;&CenterDot;</mo></mover><mo>-</mo><mover><mi>&beta;</mi><mo>&CenterDot;&CenterDot;</mo></mover><mo>)</mo></mrow><mo>+</mo><msub><mover><mi>&theta;</mi><mo>&CenterDot;&CenterDot;</mo></mover><mi>g</mi></msub><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><msub><mi>H</mi><mrow><mi>r</mi><mi>z</mi></mrow></msub><mo>&lsqb;</mo><mfrac><msqrt><mn>2</mn></msqrt><mn>2</mn></mfrac><mrow><mo>(</mo><mover><mi>&alpha;</mi><mo>&CenterDot;</mo></mover><mo>+</mo><mover><mi>&beta;</mi><mo>&CenterDot;</mo></mover><mo>)</mo></mrow><mo>+</mo><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>j</mi></msub><msub><mi>cos&theta;</mi><mi>g</mi></msub><mo>&rsqb;</mo><mo>+</mo><msub><mi>T</mi><mrow><mi>f</mi><mi>x</mi></mrow></msub><mo>+</mo><msub><mi>w</mi><mi>x</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>T</mi><mrow><mi>j</mi><mi>y</mi></mrow></msub><mo>=</mo><msub><mi>K</mi><mrow><mi>i</mi><mi>j</mi><mi>y</mi></mrow></msub><msub><mi>i</mi><mrow><mi>j</mi><mi>y</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>J</mi><mrow><mi>j</mi><mi>y</mi></mrow></msub><mo>+</mo><msub><mi>J</mi><mrow><mi>g</mi><mi>y</mi></mrow></msub><msup><mi>cos</mi><mn>2</mn></msup><msub><mi>&theta;</mi><mi>g</mi></msub><mo>+</mo><msub><mi>J</mi><mrow><mi>g</mi><mi>z</mi></mrow></msub><msup><mi>sin</mi><mn>2</mn></msup><msub><mi>&theta;</mi><mi>g</mi></msub><mo>+</mo><msub><mi>J</mi><mrow><mi>r</mi><mi>r</mi></mrow></msub><msup><mi>cos</mi><mn>2</mn></msup><msub><mi>&theta;</mi><mi>g</mi></msub><mo>)</mo></mrow><msub><mover><mi>&theta;</mi><mo>&CenterDot;&CenterDot;</mo></mover><mi>j</mi></msub><mo>+</mo><mfrac><msqrt><mn>2</mn></msqrt><mn>2</mn></mfrac><msub><mi>J</mi><mrow><mi>r</mi><mi>r</mi></mrow></msub><mrow><mo>(</mo><mover><mi>&alpha;</mi><mo>&CenterDot;&CenterDot;</mo></mover><mo>+</mo><mover><mi>&beta;</mi><mo>&CenterDot;&CenterDot;</mo></mover><mo>)</mo></mrow><msub><mi>cos&theta;</mi><mi>g</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><msub><mi>H</mi><mrow><mi>r</mi><mi>z</mi></mrow></msub><mo>&lsqb;</mo><mfrac><msqrt><mn>2</mn></msqrt><mn>2</mn></mfrac><mrow><mo>(</mo><mover><mi>&alpha;</mi><mo>&CenterDot;</mo></mover><mo>-</mo><mover><mi>&beta;</mi><mo>&CenterDot;</mo></mover><mo>)</mo></mrow><mo>+</mo><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>g</mi></msub><mo>&rsqb;</mo><msub><mi>cos&theta;</mi><mi>g</mi></msub><mo>-</mo><mrow><mo>(</mo><msub><mi>J</mi><mrow><mi>r</mi><mi>r</mi></mrow></msub><mo>+</mo><mn>2</mn><msub><mi>J</mi><mrow><mi>g</mi><mi>y</mi></mrow></msub><mo>-</mo><mn>2</mn><msub><mi>J</mi><mrow><mi>g</mi><mi>z</mi></mrow></msub><mo>)</mo></mrow><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>g</mi></msub><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>j</mi></msub><msub><mi>sin&theta;</mi><mi>g</mi></msub><msub><mi>cos&theta;</mi><mi>g</mi></msub><mo>+</mo><msub><mi>T</mi><mrow><mi>f</mi><mi>y</mi></mrow></msub><mo>+</mo><msub><mi>w</mi><mi>y</mi></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000796257580000011.GIF" wi="1438" he="499" /></maths>其中,<img file="FDA0000796257580000012.GIF" wi="53" he="76" />为外框轴相对惯性系的转动角速率;<img file="FDA0000796257580000013.GIF" wi="49" he="80" />为内框轴相对外框系的转动角速率,θ<sub>g</sub>为内框轴相对外框系转动角位置,<img file="FDA0000796257580000014.GIF" wi="54" he="83" />为外框轴相对惯性系的转动角加速率;<img file="FDA0000796257580000015.GIF" wi="54" he="83" />为内框轴相对外框系的转动角加速率,<img file="FDA0000796257580000016.GIF" wi="140" he="77" />分别为高速转子x、y方向的扭转速度,<img file="FDA0000796257580000017.GIF" wi="143" he="77" />分别为高速转子x、y方向的扭转加速度,H<sub>rz</sub>为高速转子的角动量,T<sub>gx</sub>和T<sub>jy</sub>分别为内、外框架电机的输出力矩,K<sub>igx</sub>和K<sub>ijy</sub>分别为内、外框架电机的力矩系数,i<sub>gx</sub>和i<sub>jy</sub>分别为内、外框架电机绕组电流;T<sub>fx</sub>为作用在内框架转动轴的非线性摩擦力矩,T<sub>fy</sub>为作用在外框架转动轴的非线性摩擦力矩,J<sub>jy</sub>为外框架输出力矩方向的转动惯量;J<sub>gx</sub>、J<sub>gy</sub>、J<sub>gz</sub>分别为内框架对内框坐标系相应各轴的转动惯量;J<sub>rr</sub>为高速转子径向的转动惯量,w<sub>x</sub>和w<sub>y</sub>分别为内、外框架系统的未知扰动力矩;由于转子运动被限制在保护间隙内,而且高速转子转速<img file="FDA0000796257580000018.GIF" wi="351" he="77" />J<sub>rr</sub>&lt;J<sub>rz</sub>,J<sub>rz</sub>为高速转子轴向的转动惯量,因而忽略转子径向运动的影响,得到内、外框架系统的简化动力学模型如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>T</mi><mrow><mi>g</mi><mi>x</mi></mrow></msub><mo>=</mo><msub><mi>K</mi><mrow><mi>i</mi><mi>g</mi><mi>x</mi></mrow></msub><msub><mi>i</mi><mrow><mi>g</mi><mi>x</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>J</mi><mrow><mi>g</mi><mi>x</mi></mrow></msub><mo>+</mo><msub><mi>J</mi><mrow><mi>r</mi><mi>r</mi></mrow></msub><mo>)</mo></mrow><msub><mover><mi>&theta;</mi><mo>&CenterDot;&CenterDot;</mo></mover><mi>g</mi></msub><mo>+</mo><msub><mi>H</mi><mrow><mi>r</mi><mi>z</mi></mrow></msub><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>j</mi></msub><msub><mi>cos&theta;</mi><mi>g</mi></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>J</mi><mrow><mi>g</mi><mi>y</mi></mrow></msub><mo>-</mo><msub><mi>J</mi><mrow><mi>g</mi><mi>z</mi></mrow></msub><mo>)</mo></mrow><msubsup><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>j</mi><mn>2</mn></msubsup><msub><mi>sin&theta;</mi><mi>g</mi></msub><msub><mi>cos&theta;</mi><mi>g</mi></msub><mo>+</mo><msub><mi>T</mi><mrow><mi>f</mi><mi>x</mi></mrow></msub><mo>+</mo><msub><mi>w</mi><mi>x</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>T</mi><mrow><mi>j</mi><mi>y</mi></mrow></msub><mo>=</mo><msub><mi>K</mi><mrow><mi>i</mi><mi>j</mi><mi>y</mi></mrow></msub><msub><mi>i</mi><mrow><mi>j</mi><mi>y</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>J</mi><mrow><mi>j</mi><mi>y</mi></mrow></msub><mo>+</mo><msub><mi>J</mi><mrow><mi>g</mi><mi>y</mi></mrow></msub><msup><mi>cos</mi><mn>2</mn></msup><msub><mi>&theta;</mi><mi>g</mi></msub><mo>+</mo><msub><mi>J</mi><mrow><mi>g</mi><mi>z</mi></mrow></msub><msup><mi>sin</mi><mn>2</mn></msup><msub><mi>&theta;</mi><mi>g</mi></msub><mo>+</mo><msub><mi>J</mi><mrow><mi>r</mi><mi>r</mi></mrow></msub><msup><mi>cos</mi><mn>2</mn></msup><msub><mi>&theta;</mi><mi>g</mi></msub><mo>)</mo></mrow><msub><mover><mi>&theta;</mi><mo>&CenterDot;&CenterDot;</mo></mover><mi>j</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mrow><mo>(</mo><msub><mi>J</mi><mrow><mi>r</mi><mi>r</mi></mrow></msub><mo>+</mo><mn>2</mn><msub><mi>J</mi><mrow><mi>g</mi><mi>y</mi></mrow></msub><mo>-</mo><mn>2</mn><msub><mi>J</mi><mrow><mi>g</mi><mi>z</mi></mrow></msub><mo>)</mo></mrow><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>g</mi></msub><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>j</mi></msub><msub><mi>sin&theta;</mi><mi>g</mi></msub><msub><mi>cos&theta;</mi><mi>g</mi></msub><mo>-</mo><msub><mi>H</mi><mrow><mi>r</mi><mi>z</mi></mrow></msub><msub><mover><mi>&theta;</mi><mo>&CenterDot;</mo></mover><mi>g</mi></msub><msub><mi>cos&theta;</mi><mi>g</mi></msub><mo>+</mo><msub><mi>T</mi><mrow><mi>f</mi><mi>y</mi></mrow></msub><mo>+</mo><msub><mi>w</mi><mi>y</mi></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000796257580000019.GIF" wi="1380" he="234" /></maths>由上式可知,影响框架系统伺服精度的主要因素是陀螺耦合力矩和非线性摩擦力矩,陀螺耦合力矩可以根据动力学模型采用计算力矩的方法进行前馈或者反馈补偿;由于陀螺效应,框架摩擦力矩受到陀螺耦合力矩的影响,与框架角速率和角位置都有关系,框架伺服系统的非线性摩擦力矩模型更加复杂,不易进行建模和参数辨识,因而基于摩擦模型的补偿方法无法实现框架伺服系统的高精度控制;步骤(2)、根据所述步骤(1)中的内、外框架系统的动力学模型,首先对框架耦合力矩进行前馈补偿控制,由于框架系统非线性摩擦力矩模型未知,应用TDE(time delay estimation)对非线性摩擦力矩进行估计,并采用TDC(time delay control)对非线性摩擦进行反馈补偿;步骤(3)、根据所述步骤(2)中非线性摩擦补偿后的系统,由于TDE存在估计误差,TDC对突变扰动力矩的鲁棒性不强,采用两自由度IMC(internal model control)对TDE引入的估计误差进行补偿,同时实现了框架系统跟踪性能和扰动抑制性能的独立控制;步骤(4)、根据所述步骤(2)和(3)中的TDC与两自由度IMC,将二者相结合,并化为经典反馈控制形式,实现框架系统跟踪性能和扰动抑制性能的独立控制及框架系统高精度角速率跟踪。
地址 100191 北京市海淀区学院路37号
您可能感兴趣的专利