发明名称 一种四轮独立驱动电动汽车的横摆角速度控制方法
摘要 本发明公开了一种四轮独立驱动电动汽车的横摆角速度控制方法,通过设计跟踪微分器获得电动汽车横摆角速度设定值的跟踪信号和该跟踪信号的微分。同时,设计扩张状态观测器估计电动汽车系统的状态参数和不确定扰动。接着,根据前两步得到的值计算出控制系统的误差信号和误差微分信号,并将这两个信号带入非线性组合函数,得到误差反馈控制量,再用扩张状态观测器所得到的电动汽车系统不确定扰动的估计值对误差反馈控制量进行补偿,得到附加横摆力矩值,最后,通过这个附加横摆力矩值对各个车轮转矩进行分配,给汽车四个电机输入分配的指令转矩,从而控制电动汽车的横摆角速度。
申请公布号 CN104590253B 申请公布日期 2017.04.05
申请号 CN201410781886.8 申请日期 2014.12.16
申请人 电子科技大学 发明人 辛晓帅;陈锐;邹见效;徐红兵
分类号 B60W30/045(2012.01)I;B60W10/08(2006.01)I 主分类号 B60W30/045(2012.01)I
代理机构 成都行之专利代理事务所(普通合伙) 51220 代理人 温利平
主权项 一种四轮独立驱动电动汽车的横摆角速度控制方法,其特征在于,包括以下步骤:(1)、设计跟踪微分器u<sup>*</sup>(t)<maths num="0001"><math><![CDATA[<mrow><msup><mi>u</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>x</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>f</mi><mi>h</mi><mi>a</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>-</mo><mi>v</mi><mo>(</mo><mi>t</mi><mo>)</mo><mo>,</mo><msub><mi>x</mi><mn>2</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>,</mo><mi>r</mi><mo>,</mo><mi>h</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001116494710000011.GIF" wi="878" he="155" /></maths>其中,v(t)为电动汽车横摆角速度的设定值,x<sub>1</sub>(t)是对设定值v(t)的跟踪信号,x<sub>2</sub>(t)为x<sub>1</sub>(t)的微分,h为积分步长,r为决定跟踪速度的速度因子,fhan(x<sub>1</sub>',x'<sub>2</sub>,r,h<sub>0</sub>)是最速控制综合函数,<img file="FDA0001116494710000012.GIF" wi="798" he="214" />其中,<img file="FDA0001116494710000013.GIF" wi="733" he="566" />(2)、设计扩展状态观测器u'(t)<maths num="0002"><math><![CDATA[<mrow><msup><mi>u</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>E</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>Z</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>Z</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>Z</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&beta;</mi><mn>01</mn></msub><mi>E</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>Z</mi><mo>&CenterDot;</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>Z</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&beta;</mi><mn>02</mn></msub><mi>f</mi><mi>a</mi><mi>l</mi><mrow><mo>(</mo><mi>e</mi><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>,</mo><mi>h</mi><mo>)</mo></mrow><mo>+</mo><mi>b</mi><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>Z</mi><mo>&CenterDot;</mo></mover><mn>3</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><msub><mi>&beta;</mi><mn>03</mn></msub><mi>f</mi><mi>a</mi><mi>l</mi><mrow><mo>(</mo><mi>e</mi><mo>,</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><mo>,</mo><mi>h</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001116494710000014.GIF" wi="813" he="407" /></maths>其中,<img file="FDA0001116494710000015.GIF" wi="630" he="223" />在电动汽车系统,以u(t)作为输入,以y(t)作为输出,则该系统可表示为如下的二阶系统形式:<maths num="0003"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mover><mi>X</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>X</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>X</mi><mo>&CenterDot;</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>X</mi><mn>1</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>,</mo><msub><mi>X</mi><mn>2</mn></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow><mo>+</mo><mi>b</mi><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>X</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001116494710000016.GIF" wi="646" he="255" /></maths>其中,u(t)为附加横摆力矩值,y(t)为附加横摆力矩值u(t)对横摆角速度的控制结果;Z<sub>1</sub>(t)、Z<sub>2</sub>(t)是对该二阶系统中状态变量X<sub>1</sub>(t),X<sub>2</sub>(t)的估计值,Z<sub>3</sub>(t)是对该二阶系统中不确定扰动f(X<sub>1</sub>(t),X<sub>2</sub>(t))的估计值,<img file="FDA0001116494710000021.GIF" wi="126" he="71" />i=1,2,3是Z<sub>i</sub>(t)的微分;当积分步长h给定时,扩张状态观测器u'(t)的参数β<sub>01</sub>β<sub>02</sub>β<sub>03</sub>按下列公式确定:<maths num="0004"><math><![CDATA[<mrow><msub><mi>&beta;</mi><mn>01</mn></msub><mo>&ap;</mo><mfrac><mn>1</mn><mi>h</mi></mfrac></mrow>]]></math><img file="FDA0001116494710000022.GIF" wi="153" he="118" /></maths><maths num="0005"><math><![CDATA[<mrow><msub><mi>&beta;</mi><mn>02</mn></msub><mo>=</mo><mfrac><mn>1</mn><mrow><msup><mi>e</mi><mn>0.4762</mn></msup><msup><mi>h</mi><mn>1.4673</mn></msup></mrow></mfrac><mo>&ap;</mo><mfrac><mn>1</mn><mrow><mn>1.6</mn><msup><mi>h</mi><mn>1.5</mn></msup></mrow></mfrac></mrow>]]></math><img file="FDA0001116494710000023.GIF" wi="525" he="118" /></maths><maths num="0006"><math><![CDATA[<mrow><msub><mi>&beta;</mi><mn>03</mn></msub><mo>=</mo><mfrac><mn>1</mn><mrow><msup><mi>e</mi><mn>2.1567</mn></msup><msup><mi>h</mi><mn>2.2093</mn></msup></mrow></mfrac><mo>&ap;</mo><mfrac><mn>1</mn><mrow><mn>8.6</mn><msup><mi>h</mi><mn>2.2</mn></msup></mrow></mfrac></mrow>]]></math><img file="FDA0001116494710000024.GIF" wi="533" he="119" /></maths>(3)、计算误差信号根据跟踪微分器u<sup>*</sup>(t)得到的跟踪信号x<sub>1</sub>(t)、x<sub>2</sub>(t),以及扩展状态观测器u'(t)的估计值Z<sub>1</sub>(t)、Z<sub>2</sub>(t),分别求出控制系统的误差信号e<sub>1</sub>(t)和误差微分信号e<sub>2</sub>(t):e<sub>1</sub>(t)=x<sub>1</sub>(t)‑Z<sub>1</sub>(t)e<sub>2</sub>(t)=x<sub>2</sub>(t)‑Z<sub>2</sub>(t)(4)、误差信号的非线性组合将误差信号e<sub>1</sub>(t)和误差微分信号e<sub>2</sub>(t)代入到非线性组合函数k(),得到误差反馈控制量u<sub>0</sub>(t):u<sub>0</sub>(t)=k(e<sub>1</sub>(t),e<sub>2</sub>(t),p)其中,p为非线性组合函数对应的参数;(5)、通过附加横摆力矩值u(t)控制电动汽车的横摆角速度用估计值Z<sub>3</sub>(t)对误差反馈控制量u<sub>0</sub>(t)进行补偿,得到附加横摆力矩值u(t)<maths num="0007"><math><![CDATA[<mrow><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>u</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><msub><mi>Z</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><msub><mi>b</mi><mn>0</mn></msub></mfrac></mrow>]]></math><img file="FDA0001116494710000025.GIF" wi="382" he="127" /></maths>其中,b<sub>0</sub>是决定补偿强弱的“补偿因子”;接下来,将附加横摆力矩值u(t)在车轮间进行力矩分配,得到四个车轮的指令转矩,再将这4个指令转矩输入给汽车车轮的四个电机,从而控制电动汽车的横摆角速度。
地址 611731 四川省成都市高新区(西区)西源大道2006号