发明名称 电液负载模拟器多余力的舵机指令动态补偿控制方法
摘要 本发明涉及一种电液负载模拟器多余力的舵机指令动态补偿控制方法,包括:对电液负载模拟器系统建模;对电液负载模拟器系统的结构变化形式,将舵机控制系统及加载系统分为第一部分及第二部分;第一部分为舵机做位置控制,舵机受到的负载力矩大小与加载系统力矩传感器采集相等;第二部分为加载系统做力矩闭环控制,加载系统受到的运动干扰与舵机角度传感器采集相等;获取加载系统的输出力矩传递函数;将舵机的运动干扰代入加载系统的输出力矩传递函数式中;令基于舵机指令信号消除多余力的控制输出与生成加载系统的多余力和为零,消除加载系统多余力。本发明能够有效抑制舵机不同频率运动时的多余力。
申请公布号 CN106055753A 申请公布日期 2016.10.26
申请号 CN201610346560.1 申请日期 2016.05.24
申请人 北京航空航天大学 发明人 焦宗夏;李志辉;尚耀星;吴帅
分类号 G06F17/50(2006.01)I;G09B9/08(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京华创博为知识产权代理有限公司 11551 代理人 管莹;张波涛
主权项 一种电液负载模拟器多余力的舵机指令动态补偿控制方法,包括:对电液负载模拟器系统建模,所述电液负载模拟器系统包括加载系统及舵机控制系统相互耦合,所述加载系统包括加载系统伺服阀、加载系统控制器、加载马达、以及力矩传感器;所述舵机控制系统作为位置伺服系统,包括舵机控制器、舵机伺服阀、惯量盘、角度编码器以及舵机马达;对电液负载模拟器系统的结构变化形式,将舵机控制系统及加载系统分为第一部分及第二部分;第一部分为舵机做位置控制,舵机受到的负载力矩大小与加载系统力矩传感器采集相等;第二部分为加载系统做力矩闭环控制,加载系统受到的运动干扰与舵机角度编码器采集相等;获取加载系统的输出力矩传递函数:<maths num="0001"><math><![CDATA[<mrow><mi>T</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>G</mi><mi>s</mi></msub><mo>&CenterDot;</mo><msub><mi>G</mi><mrow><mi>L</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>G</mi><mrow><mi>L</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><msub><mi>u</mi><mi>L</mi></msub></mrow><mrow><msub><mi>G</mi><mrow><mi>L</mi><mn>4</mn></mrow></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mfrac><mrow><msub><mi>G</mi><mi>s</mi></msub><mo>&CenterDot;</mo><msub><mi>G</mi><mrow><mi>L</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>G</mi><mrow><mi>L</mi><mn>3</mn></mrow></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><msub><mi>s&theta;</mi><mi>a</mi></msub></mrow><mrow><msub><mi>G</mi><mrow><mi>L</mi><mn>4</mn></mrow></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000997883080000011.GIF" wi="1453" he="206" /></maths>式中:G<sub>L1</sub>(s)=D<sub>L</sub>K<sub>qL</sub>,G<sub>L2</sub>(s)=J<sub>f</sub>s<sup>2</sup>+B<sub>f</sub>s+G<sub>f</sub>,<maths num="0002"><math><![CDATA[<mrow><msub><mi>G</mi><mrow><mi>L</mi><mn>3</mn></mrow></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mfrac><msub><mi>V</mi><mi>t</mi></msub><mrow><mn>4</mn><msub><mi>&beta;</mi><mi>e</mi></msub></mrow></mfrac><mi>s</mi><mo>+</mo><msub><mi>K</mi><mrow><mi>c</mi><mi>L</mi></mrow></msub><mo>+</mo><msub><mi>K</mi><mrow><mi>t</mi><mi>p</mi><mi>L</mi></mrow></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>J</mi><mi>L</mi></msub><mi>s</mi><mo>+</mo><msub><mi>B</mi><mi>L</mi></msub><mo>)</mo></mrow><mo>+</mo><msup><msub><mi>D</mi><mi>L</mi></msub><mn>2</mn></msup><mo>,</mo></mrow>]]></math><img file="FDA0000997883080000012.GIF" wi="860" he="150" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mi>G</mi><mrow><mi>L</mi><mn>4</mn></mrow></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mfrac><msub><mi>V</mi><mi>t</mi></msub><mrow><mn>4</mn><msub><mi>&beta;</mi><mi>e</mi></msub></mrow></mfrac><mi>s</mi><mo>+</mo><msub><mi>K</mi><mrow><mi>c</mi><mi>L</mi></mrow></msub><mo>+</mo><msub><mi>K</mi><mrow><mi>t</mi><mi>p</mi><mi>L</mi></mrow></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>G</mi><mrow><mi>L</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>G</mi><mi>s</mi></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>G</mi><mrow><mi>L</mi><mn>2</mn></mrow></msub><mo>(</mo><mi>s</mi><mo>)</mo><mo>+</mo><msub><mi>G</mi><mi>s</mi></msub><mo>)</mo></mrow><msub><mi>G</mi><mrow><mi>L</mi><mn>3</mn></mrow></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>s</mi><mo>;</mo></mrow>]]></math><img file="FDA0000997883080000013.GIF" wi="1341" he="135" /></maths>式中:G<sub>s</sub>为力矩传感器与传动轴的综合刚度;u<sub>L</sub>为加载系统伺服阀的控制电压;s为拉普拉斯算子;θ<sub>a</sub>为舵机马达转子的角位移;D<sub>L</sub>为加载液压马达的排量;K<sub>qL</sub>为加载系统伺服阀的流量增益;J<sub>f</sub>为惯性负载的转动惯量;B<sub>f</sub>为加载系统负载粘性阻尼系数;G<sub>f</sub>为力矩传感器与角度传感器的综合刚度;V<sub>L</sub>为加载液压马达的容积;β<sub>e</sub>为液压油的弹性模量;C<sub>slL</sub>为加载液压马达的泄漏系数;KcL为加载系统伺服阀流量压力系数;J<sub>L</sub>为加载液压马达转动惯量;B<sub>L</sub>为加载液压马达粘性阻尼系数;将舵机的运动干扰代入加载系统的输出力矩传递函数式中,得到:<maths num="0004"><math><![CDATA[<mrow><mi>T</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>G</mi><mi>s</mi></msub><msub><mi>G</mi><mrow><mi>L</mi><mn>2</mn></mrow></msub><msub><mi>G</mi><mrow><mi>L</mi><mn>6</mn></mrow></msub></mrow><mrow><msub><mi>G</mi><mrow><mi>L</mi><mn>4</mn></mrow></msub><msub><mi>G</mi><mrow><mi>L</mi><mn>6</mn></mrow></msub><mo>-</mo><mn>1</mn></mrow></mfrac><mo>&lsqb;</mo><msub><mi>G</mi><mrow><mi>L</mi><mn>1</mn></mrow></msub><msub><mi>u</mi><mi>L</mi></msub><mo>-</mo><mfrac><mrow><msub><mi>K</mi><mrow><mi>q</mi><mi>a</mi></mrow></msub><msub><mi>G</mi><mrow><mi>L</mi><mn>3</mn></mrow></msub><msub><mi>D</mi><mi>a</mi></msub><msub><mi>&theta;</mi><mrow><mi>i</mi><mi>n</mi></mrow></msub></mrow><msub><mi>G</mi><mrow><mi>L</mi><mn>5</mn></mrow></msub></mfrac><mo>&rsqb;</mo><mo>,</mo></mrow>]]></math><img file="FDA0000997883080000021.GIF" wi="1078" he="182" /></maths>式中:<img file="FDA0000997883080000022.GIF" wi="918" he="191" /><maths num="0005"><math><![CDATA[<mrow><msub><mi>G</mi><mrow><mi>L</mi><mn>6</mn></mrow></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>G</mi><mrow><mi>L</mi><mn>5</mn></mrow></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>G</mi><mi>s</mi></msub><msub><mi>G</mi><mrow><mi>L</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><msub><mi>G</mi><mrow><mi>L</mi><mn>3</mn></mrow></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mrow><mo>(</mo><mfrac><msub><mi>V</mi><mi>a</mi></msub><mrow><mn>4</mn><msub><mi>&beta;</mi><mi>e</mi></msub></mrow></mfrac><mi>s</mi><mo>+</mo><msub><mi>K</mi><mi>c</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000997883080000023.GIF" wi="910" he="253" /></maths>K<sub>c</sub>=K<sub>ca</sub>+C<sub>sla</sub>;式中:K<sub>qa</sub>为舵机伺服阀的流量增益;D<sub>a</sub>为舵机液压马达排量;θ<sub>in</sub>为舵机控制系统的运动指令;V<sub>a</sub>为舵机液压马达控制容积;K<sub>c</sub>为舵机控制系统的总流量压力系数;K<sub>ca</sub>为舵机伺服阀流量压力系数;C<sub>sla</sub>为舵机液压马达泄漏系数;J<sub>a</sub>为舵机液压马达转动惯量;B<sub>a</sub>为舵机液压马达粘性阻尼系数;基于动态舵机指令前馈补偿控制的加载系统控制输出如下:u<sub>L</sub>=u<sub>pid</sub>+G<sub>d</sub>·θ<sub>in</sub>,式中:u<sub>pid</sub>为加载系统力矩闭环PID控制输出,G<sub>d</sub>为动态舵机指令前馈补偿环节。将式(3)代入式(2)中,可得:<maths num="0006"><math><![CDATA[<mrow><mi>T</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>G</mi><mi>s</mi></msub><msub><mi>G</mi><mrow><mi>L</mi><mn>2</mn></mrow></msub><msub><mi>G</mi><mrow><mi>L</mi><mn>6</mn></mrow></msub></mrow><mrow><msub><mi>G</mi><mrow><mi>L</mi><mn>4</mn></mrow></msub><msub><mi>G</mi><mrow><mi>L</mi><mn>6</mn></mrow></msub><mo>-</mo><mn>1</mn></mrow></mfrac><mrow><mo>&lsqb;</mo><mrow><msub><mi>G</mi><mrow><mi>L</mi><mn>1</mn></mrow></msub><msub><mi>u</mi><mrow><mi>p</mi><mi>i</mi><mi>d</mi></mrow></msub><mo>+</mo><msub><mi>G</mi><mrow><mi>L</mi><mn>1</mn></mrow></msub><msub><mi>G</mi><mrow><mi>d</mi><mi>s</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>&theta;</mi><mrow><mi>i</mi><mi>n</mi></mrow></msub><mo>-</mo><mfrac><mrow><msub><mi>K</mi><mrow><mi>q</mi><mi>a</mi></mrow></msub><msub><mi>G</mi><mrow><mi>L</mi><mn>3</mn></mrow></msub><msub><mi>D</mi><mi>a</mi></msub><msub><mi>&theta;</mi><mrow><mi>i</mi><mi>n</mi></mrow></msub></mrow><msub><mi>G</mi><mrow><mi>L</mi><mn>5</mn></mrow></msub></mfrac></mrow><mo>&rsqb;</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000997883080000024.GIF" wi="1277" he="166" /></maths>上式的等式右边第一项为加载系统的PID控制输出,第二项为基于舵机指令信号消除多余力的控制输出,第三项为包含舵机动态性能的对加载系统的影响,即生成加载系统的多余力;令第二项与第三项之和为零,消除加载系统多余力。
地址 100191 北京市海淀区学院路37号