发明名称 一种基于磁通切换直线电机的轨道交通速度脉动抑制方法
摘要 本发明公开了一种基于磁通切换直线电机的轨道交通速度脉动抑制方法,地铁车辆运行时可分为变速运行和恒速运行两者情况,为达到抑制定位力和速度脉动的目的,当处于变速运行时本文提出了位置x域重复控制和R‑PI调节器结合的方法来抑制定位力的新型控制策略,首先将时变的t域的信号,变换为具有固定周期的x域信号,再对其进行重复控制器的设计;当处于恒速运行时,采用时间t域重复控制和R‑PI调节器结合的方法。同时为了准确获得LFSPM的位置和速度信息,本专利同时提出了位置递推控制算法来取代位置传感器技术。
申请公布号 CN104104297B 申请公布日期 2016.09.07
申请号 CN201410340209.2 申请日期 2014.07.16
申请人 东南大学 发明人 余海涛;孟高军;胡敏强;黄磊
分类号 H02P21/05(2006.01)I;H02P21/18(2016.01)I;H02P25/06(2016.01)I 主分类号 H02P21/05(2006.01)I
代理机构 南京瑞弘专利商标事务所(普通合伙) 32249 代理人 黄成萍
主权项 一种基于磁通切换直线电机的轨道交通速度脉动抑制方法,其特征在于:根据轨道车车辆速度不同而切换x域重复控制器或t域重复控制器,同时与位置递推控制算法相结合,最终达到抑制定位力和速度脉动的目的,同时实现了无位置传感器技术;具体包括如下步骤:(1)监测LFSPM电机的实际运行速度v(t):若v(t)为恒定值,则进入步骤(2),速度环采用t域重复控制器和R‑PI调节器结合的方法进行控制;若v(t)为变速值,则进入步骤(3),速度环采用x域重复控制器和R‑PI调节器结合的方法进行控制;(2)恒速时,速度环的控制方法:恒速时,速度环采用t域重复控制器和R‑PI调节器结合的方法进行控制:R‑PI调节器的传递函数G<sub>R‑PI</sub>(s)为:<maths num="0001"><math><![CDATA[<mrow><msub><mi>G</mi><mrow><mi>R</mi><mo>-</mo><mi>P</mi><mi>I</mi></mrow></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><msub><mi>K</mi><mi>p</mi></msub><mo>+</mo><mfrac><msub><mi>K</mi><mi>i</mi></msub><mi>s</mi></mfrac><mo>+</mo><mfrac><mrow><mn>2</mn><msub><mi>w</mi><mi>c</mi></msub><msub><mi>K</mi><mi>R</mi></msub><mi>s</mi></mrow><mrow><msup><mi>s</mi><mn>2</mn></msup><mo>+</mo><mn>2</mn><msub><mi>w</mi><mi>c</mi></msub><mi>s</mi><mo>+</mo><msubsup><mi>w</mi><mi>o</mi><mn>2</mn></msubsup></mrow></mfrac></mrow>]]></math><img file="FDA0001050515460000011.GIF" wi="692" he="135" /></maths>其中,K<sub>p</sub>为比例系数;K<sub>i</sub>为积分系数;K<sub>R</sub>为谐振系数;w<sub>o</sub>为谐振频率;w<sub>c</sub>为低通截止频率,s为复变量;t域重复控制器:首先将<img file="FDA0001050515460000012.GIF" wi="72" he="63" />与y(s)做差得到e(s),使用二阶低通滤波器Q(s)和时滞环节<img file="FDA0001050515460000013.GIF" wi="90" he="54" />构成重复控制器对e(s)进行重复控制,其中T<sub>0</sub>是e(s)输入重复控制器的周期;e(s)经重复控制器后的输出经过串联补偿器C(s)控制P(s)得到y(s),P(s)为控制对象;最后对y(s)进行拉普拉斯逆变换得到y(t),v<sub>ref</sub>=y(t);(3)变速时,速度环的控制方法:首先对<img file="FDA0001050515460000014.GIF" wi="66" he="68" />进行拉普拉斯变换得到v(s),将v(s)与y(s)做差得到e(s),对e(s)进行拉普拉斯逆变换得到e(t);然后设定有线性算子T与逆算子T<sup>‑1</sup>,计算<img file="FDA0001050515460000015.GIF" wi="283" he="71" />再通过k<sub>RP</sub>调节闭环系统的收敛速度得到<img file="FDA0001050515460000016.GIF" wi="131" he="77" />x为动子齿轴线所在的位置;<img file="FDA0001050515460000017.GIF" wi="878" he="71" />D(x)为时滞环节,Q(x)为低通滤波器;接着对u(t)进行拉普拉斯变换得到u(s),经过串联补偿器C(s)控制P(s)得到y(s),P(s)为控制对象;最后对y(s)进行拉普拉斯逆变换得到y(t),v<sub>ref</sub>=y(t);(4)v<sub>ref</sub>与v(t)的差值作为PI控制器的输入,得到q轴参考电流i<sub>qref</sub>,i<sub>qref</sub>与q轴实际电流i<sub>q</sub>的差值作为PI控制器的输入,得到q轴电压u<sub>q</sub>;d轴参考电流i<sub>dref</sub>=0与d轴实际电流i<sub>d</sub>的差值作为PI控制器的输入,得到d轴参考电压u<sub>d</sub>;对u<sub>q</sub>和u<sub>d</sub>进行dq/αβ变换得到u<sub>α</sub>和u<sub>β</sub>,以u<sub>α</sub>和u<sub>β</sub>作为参考值,对LFSPM电机进行SVPWM控制,检测得到逆变器的输出电流i<sub>a</sub>、i<sub>b</sub>和i<sub>c</sub>;(5)对i<sub>a</sub>、i<sub>b</sub>和i<sub>c</sub>进行abc/αβ变换得到i<sub>α</sub>和i<sub>β</sub>,对i<sub>α</sub>和i<sub>β</sub>进行αβ/dq变换得到i<sub>q</sub>和i<sub>d</sub>;以i<sub>a</sub>、i<sub>b</sub>和i<sub>c</sub>驱动LFSPM电机;(6)通过无位置传感器对LFSPM电机的位置信息进行检测,具体为:设定θ<sub>m</sub>为LFSPM电机的位置角,τ为相邻两个定子之间的极距,则LFSPM电机的电压方程在两相旋转坐标系dq轴坐标下表示为:<maths num="0002"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>u</mi><mi>d</mi></msub></mtd></mtr><mtr><mtd><msub><mi>u</mi><mi>q</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mi>R</mi><mo>+</mo><mi>p</mi><mi>L</mi></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>R</mi><mo>+</mo><mi>p</mi><mi>L</mi></mrow></mtd></mtr></mtable></mfenced><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>i</mi><mi>d</mi></msub></mtd></mtr><mtr><mtd><msub><mi>i</mi><mi>q</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mi>p</mi><mi>&Phi;</mi><mfenced open = "[" close = "]"><mtable><mtr><mtd><mi>c</mi><mi>o</mi><mi>s</mi><mi>&Delta;</mi><msub><mi>&theta;</mi><mi>m</mi></msub></mtd></mtr><mtr><mtd><mrow><msub><mi>sin&Delta;&theta;</mi><mi>m</mi></msub></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001050515460000021.GIF" wi="933" he="151" /></maths>Δθ<sub>m</sub>=θ<sub>m</sub>(n)‑θ<sub>m</sub>(n‑1)其中:R为电枢绕组电阻,L为电枢绕组电感,Φ为永磁体的永磁磁链,p是微分因子;据此得到:<maths num="0003"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><msub><mi>cos&Delta;&theta;</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>cos&Delta;&theta;</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><mfrac><msub><mi>T</mi><mi>s</mi></msub><mi>&Phi;</mi></mfrac><mo>{</mo><msub><mi>u</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>Ri</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mi>L</mi><mfrac><mrow><msub><mi>i</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>i</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><msub><mi>T</mi><mi>s</mi></msub></mfrac><mo>}</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001050515460000022.GIF" wi="1157" he="223" /></maths><maths num="0004"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><msub><mi>sin&Delta;&theta;</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>sin&Delta;&theta;</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><mfrac><msub><mi>T</mi><mi>s</mi></msub><mi>&Phi;</mi></mfrac><mo>{</mo><msub><mi>u</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>Ri</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mi>L</mi><mfrac><mrow><msub><mi>i</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>i</mi><mi>q</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><msub><mi>T</mi><mi>s</mi></msub></mfrac><mo>}</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001050515460000023.GIF" wi="1150" he="223" /></maths>其中:n为采样点,T<sub>s</sub>为采样时间;因此得到:Δθ<sub>m</sub>(n)=tan<sup>‑1</sup>(sinΔθ<sub>m</sub>(n)/cosΔθ<sub>m</sub>(n))据此得到为LFSPM电机的位置角θ<sub>m</sub>(n)为:θ<sub>m</sub>(n)=θ<sub>m</sub>(n‑1)+Δθ<sub>m</sub>(n)(7)对θ<sub>m</sub>(n)作时间微分得到实际运行速度v(t);计算实际位置信号<img file="FDA0001050515460000031.GIF" wi="209" he="119" />将给定位置信号x<sub>ref</sub>与x相减后作为PI控制器的输入,得到<img file="FDA0001050515460000032.GIF" wi="98" he="63" />
地址 211189 江苏省南京市江宁区东南大学路2号