发明名称 基于分散控制的混合动力汽车节能预测控制方法
摘要 本发明公开了一种基于分散控制的混合动力汽车节能预测控制方法,包括以下步骤:从全球定位系统、车间通信系统、车路通信系统和智能交通系统获取实时自车、前车和后车交通信息作为系统输入;建立混合动力汽车队列行驶分散控制数学模型作为预测未来车辆状态的依据;定义混合动力汽车分散控制队列行驶最优控制问题,提供求解最优控制量的函数方程;实时反馈最优控制,求解最优控制量,运用行星齿轮机构为电子无极变速器,发动机工作于最佳工作点,运用道路交通信息,预测前车和后车行驶状态,在线调整混合动力汽车能量流动,达到节能减排的目标,大大减少了计算时间,提高了车辆的实时控制特性。
申请公布号 CN105083276A 申请公布日期 2015.11.25
申请号 CN201510410042.7 申请日期 2015.07.13
申请人 河南理工大学 发明人 余开江;许孝卓;谭兴国;刘巍;荆鹏辉;谢贝贝;胡治国;王莉;王允建
分类号 B60W20/00(2006.01)I;B60W10/06(2006.01)I;B60W10/08(2006.01)I;B60W10/26(2006.01)I 主分类号 B60W20/00(2006.01)I
代理机构 郑州科维专利代理有限公司 41102 代理人 亢志民;张欣棠
主权项 一种基于分散控制的混合动力汽车节能预测控制方法,其特征在于:包括以下步骤:步骤1)信息采集:由全球定位系统采集前车、后车和自车的位置信息,作为实时车辆状态反馈,由车载雷达测速装置采集前方车辆速度,后方车辆速度,用于跟踪控制,由智能交通系统、车路通信系统和车间通信系统采集交通信号信息,实时路况信息以及自车、后车和前车速度,加速度信息,用于智能交通控制,由卡尔曼滤波器利用采集的蓄电池信息对蓄电池荷电状态进行测定;步骤2)车辆建模:行星齿轮式混联混合动力汽车包含5大动态部件,它们是发动机、蓄电池、2个发电电动一体机和车轮,行星齿轮作为动力分配装置既有速度耦合器的作用,又有电子无极变速器作用,根据车辆机械耦合和电子耦合关系,列写系统动力学方程,对动力学方程解耦,获得系统的状态空间模型,如式(1)所示:x=[p<sub>p</sub> v<sub>p</sub> p<sub>2</sub> v<sub>2</sub> SOC<sub>2</sub> p<sub>f</sub> v<sub>f</sub>]u=[u<sub>2</sub> P<sub>batt2</sub>]<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>u</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000758306110000011.GIF" wi="222" he="62" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>u</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><msub><mi>v</mi><mi>p</mi></msub></mtd></mtr><mtr><mtd><msub><mi>a</mi><mi>p</mi></msub></mtd></mtr><mtr><mtd><msub><mi>v</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mrow><msub><mi>u</mi><mn>2</mn></msub><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>&rho;C</mi><mrow><mi>D</mi><mn>2</mn></mrow></msub><msub><mi>A</mi><mn>2</mn></msub><msup><msub><mi>v</mi><mn>2</mn></msub><mn>2</mn></msup><mo>/</mo><msub><mi>m</mi><mn>2</mn></msub><mo>-</mo><mn>9.8</mn><mi>&mu;</mi><mo>-</mo><mn>9.8</mn><msub><mi>sin&theta;</mi><mn>2</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mfrac><mrow><msub><mi>V</mi><mrow><mi>max</mi><mn>2</mn></mrow></msub><mo>-</mo><msqrt><mrow><msubsup><mi>V</mi><mrow><mi>max</mi><mn>2</mn></mrow><mn>2</mn></msubsup><mo>-</mo><mn>4</mn><msub><mi>P</mi><mrow><mi>b</mi><mi>a</mi><mi>t</mi><mi>t</mi><mn>2</mn></mrow></msub><msub><mi>R</mi><mrow><mi>b</mi><mi>a</mi><mi>t</mi><mi>t</mi><mn>2</mn></mrow></msub></mrow></msqrt></mrow><mrow><mn>2</mn><msub><mi>R</mi><mrow><mi>b</mi><mi>a</mi><mi>t</mi><mi>t</mi><mn>2</mn></mrow></msub><msub><mi>Q</mi><mrow><mi>b</mi><mi>a</mi><mi>t</mi><mi>t</mi><mn>2</mn></mrow></msub></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><msub><mi>v</mi><mi>f</mi></msub></mtd></mtr><mtr><mtd><msub><mi>a</mi><mi>f</mi></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000758306110000012.GIF" wi="966" he="646" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>a</mi><mi>p</mi></msub><mrow><mo>(</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>a</mi><mi>p</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>&beta;</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>v</mi><mi>&rho;</mi></msub><mo>(</mo><mi>&tau;</mi><mo>)</mo><mo>-</mo><msub><mi>v</mi><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub><mo>)</mo></mrow></mrow></msup><mo>)</mo><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><msub><mi>&beta;</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>v</mi><mi>&rho;</mi></msub><mo>(</mo><mi>&tau;</mi><mo>)</mo><mo>-</mo><msub><mi>v</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>)</mo></mrow></mrow></msup><mo>)</mo></mrow></mfrac></mrow>]]></math><img file="FDA0000758306110000013.GIF" wi="781" he="133" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>a</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>a</mi><mi>f</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>&beta;</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>v</mi><mi>f</mi></msub><mo>(</mo><mi>&tau;</mi><mo>)</mo><mo>-</mo><msub><mi>v</mi><mi>min</mi></msub><mo>)</mo></mrow></mrow></msup><mo>)</mo><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><msub><mi>&beta;</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>v</mi><mi>f</mi></msub><mo>(</mo><mi>&tau;</mi><mo>)</mo><mo>-</mo><msub><mi>v</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>)</mo></mrow></mrow></msup><mo>)</mo></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000758306110000014.GIF" wi="1576" he="144" /></maths>式中,x为状态量,u为控制量,参数p<sub>p</sub>、v<sub>p</sub>、p<sub>f</sub>和v<sub>f</sub>为前车位置、前车速度、后车位置和后车速度,参数p<sub>2</sub>、v<sub>2</sub>和SOC<sub>2</sub>为自车的位置、速度和蓄电池荷电状态,参数u<sub>2</sub>和P<sub>batt2</sub>为自车的驱动加速度和自车蓄电池的充放电功率,参数ρ、C<sub>D2</sub>、A<sub>2</sub>、m<sub>2</sub>、g、μ和θ<sub>2</sub>是空气密度、自车空气阻力系数、自车迎风面积、自车质量、重力加速度、滚动阻力系数和自车道路坡度,V<sub>OC</sub>、R<sub>batt</sub>和Q<sub>batt</sub>是蓄电池开路电压、内阻和容量,预测区间内由于车辆的惯性、假设前方车辆加速度一定,如果前行车速度大于最大值或者小于一定值,则前行车加速度为0,如果前方遭遇交通信号灯红灯,则假定一辆速度为0的前行车停在交通信号灯位置处,车辆的启动和停止速度模式采用实验曲线,运用实际驾驶员的特性测取;车辆的燃油经济性评价采用威兰氏线性模型,如式(2)所示:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mover><mi>m</mi><mo>&CenterDot;</mo></mover><mi>f</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mover><mi>m</mi><mo>&CenterDot;</mo></mover><mi>f</mi></msub><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>r</mi><mi>e</mi><mi>q</mi></mrow></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>-</mo><msub><mi>P</mi><mrow><mi>b</mi><mi>a</mi><mi>t</mi><mi>t</mi></mrow></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow><mo>&ap;</mo><msub><mi>c</mi><mi>f</mi></msub><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>r</mi><mi>e</mi><mi>q</mi></mrow></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>-</mo><msub><mi>P</mi><mrow><mi>b</mi><mi>a</mi><mi>t</mi><mi>t</mi></mrow></msub><mo>(</mo><mi>t</mi><mo>)</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000758306110000021.GIF" wi="1598" he="92" /></maths>式中m<sub>f</sub>为燃油消耗率,参数P<sub>req</sub>为车辆需求功率,c<sub>f</sub>为常数参数;步骤3)公式化控制策略:基于分散控制的混合动力汽车能量管理模型预测最优控制策略的步骤为:首先检测自车、后车和前车状态,包括位置、速度和加速度信息,其次运用所建立的数学模型和公式化控制策略求解最优控制问题,最后应用所求得的最优控制序列的第一个控制量于系统,由于模型预测控制为区间最优控制,其求得的最优控制量是数量为预测区间除以采样间隔的序列,最优控制序列的第一个控制量与实际状态最接近,所以采用它来作为实际的控制量;最优控制问题定义如式(3)所示:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>min</mi><mi>i</mi><mi>m</mi><mi>i</mi><mi>z</mi><mi>e</mi><mi> </mi><mi>J</mi><mo>=</mo><msubsup><mo>&Integral;</mo><mi>t</mi><mrow><mi>t</mi><mo>+</mo><mi>T</mi></mrow></msubsup><mi>L</mi><mrow><mo>(</mo><mi>x</mi><mo>(</mo><mrow><mi>&tau;</mi><mo>|</mo><mi>t</mi></mrow><mo>)</mo><mo>,</mo><mi>u</mi><mo>(</mo><mrow><mi>&tau;</mi><mo>|</mo><mi>t</mi></mrow><mo>)</mo><mo>)</mo></mrow><mi>d</mi><mi>&tau;</mi></mrow>]]></math><img file="FDA0000758306110000022.GIF" wi="793" he="108" /></maths>subject to P<sub>batt2min</sub>≤P<sub>batt2</sub>(τ|t)≤P<sub>batt2max</sub>  (3)u<sub>2min</sub>≤u<sub>2</sub>(τ|t)≤u<sub>2max</sub>式中T为预测区间,参数P<sub>batt2min</sub>、P<sub>batt2max</sub>、u<sub>2max</sub>和u<sub>2min</sub>为控制量约束;评价函数定义如式(4)所示:L=w<sub>x</sub>L<sub>x</sub>+w<sub>y</sub>L<sub>y</sub>+w<sub>z</sub>L<sub>z</sub>+w<sub>d</sub>L<sub>d</sub>+w<sub>e</sub>L<sub>e</sub>+w<sub>f</sub>L<sub>f</sub>+w<sub>r</sub>L<sub>r</sub>+w<sub>s</sub>L<sub>s</sub><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>L</mi><mi>x</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><msub><mi>w</mi><mn>2</mn></msub><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>&rho;C</mi><mrow><mi>D</mi><mn>2</mn></mrow></msub><msub><mi>A</mi><mn>2</mn></msub><msup><msub><mi>v</mi><mn>2</mn></msub><mn>2</mn></msup><mo>/</mo><msub><mi>m</mi><mn>1</mn></msub><mo>-</mo><mn>9.8</mn><mo>*</mo><mi>&mu;</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>/</mo><mn>2</mn></mrow>]]></math><img file="FDA0000758306110000031.GIF" wi="835" he="128" /></maths>L<sub>y</sub>=(v<sub>2</sub>‑v<sub>d</sub>)<sup>2</sup>/2<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>L</mi><mi>z</mi></msub><mo>=</mo><mn>0.0874</mn><mo>*</mo><mrow><mo>(</mo><msub><mi>m</mi><mn>2</mn></msub><mo>*</mo><msub><mi>w</mi><mn>2</mn></msub><mo>*</mo><msub><mi>v</mi><mn>2</mn></msub><mo>/</mo><mn>1000</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>b</mi><mi>a</mi><mi>t</mi><mi>t</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow><mo>/</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mo>-</mo><mn>0.5</mn><mo>*</mo><mrow><mo>(</mo><msub><mi>m</mi><mn>2</mn></msub><mo>*</mo><msub><mi>w</mi><mn>2</mn></msub><mo>*</mo><msub><mi>v</mi><mn>2</mn></msub><mo>/</mo><mn>1000</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>b</mi><mi>a</mi><mi>t</mi><mi>t</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow></mrow></msup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000758306110000032.GIF" wi="1320" he="85" /></maths>L<sub>d</sub>=(SOC<sub>2</sub>‑SOC<sub>d</sub>)<sup>2</sup>  (4)L<sub>e</sub>=(m<sub>2</sub>*w<sub>2</sub>*v<sub>2</sub>/1000‑P<sub>batt2</sub>)<sup>2</sup>/2L<sub>f</sub>=(‑ln[SOC<sub>2</sub>‑0.6]‑ln[0.8‑SOC<sub>2</sub>])L<sub>r</sub>=‑ln(d<sub>f</sub>‑d<sub>d</sub>)‑ln(d<sub>p</sub>‑d<sub>d</sub>)L<sub>s</sub>=(a<sub>p</sub>)<sup>2</sup>/2+(a<sub>f</sub>)<sup>2</sup>/2式中SOC<sub>d</sub>是目标蓄电池荷电状态,v<sub>d</sub>是车辆目标速度,它取值为车辆最优等速燃油经济性速度,w<sub>x</sub>、w<sub>y</sub>、w<sub>z</sub>、w<sub>d</sub>、w<sub>e</sub>、w<sub>f</sub>、w<sub>r</sub>和w<sub>s</sub>是权重系数,d<sub>d</sub>为最低车辆间距,评价函数设置使其在最低车辆间距以上浮动,从而增加控制自由度,提高车辆燃油经济性,障碍函数用于处理系统状态约束,d<sub>f</sub>为后车和自车间距,d<sub>p</sub>为前车和自车间距;步骤4)在线最优控制:为保证系统的实时最优性能,运用基于哈密顿方程的数值快速求解方法来求解上述最优控制问题,运用极小值原理将最优控制问题转化为两点边值问题,在处理哈密顿函数相关的微分方程组和代数方程组时采用部分空间法求解,这是一种GMRES解法;在每个采样时刻,首先测取前车位置、自车位置、后车位置、前车速度、自车速度、后车速度、前车加速度、自车加速度、后车加速度和自车蓄电池荷电状态实时状态信号,其次利用全球定位系统、车间通信系统、车路通信系统和智能交通系统预测未来一定区间车辆及周围环境的状态,包括前车加速度、后车加速度、前车与自车间距、后车与自车间距,再次根据建立的车辆模型和最优控制问题,利用上述数值快速解法求解预测区间内的最优控制序列,应用预测区间内的最优控制序列的第一个控制量于车辆,之后在下一个采样时刻,将预测区间向前推进一步,如此循环往复,实现在线最优控制。
地址 454000 河南省焦作市高新区世纪路2001号河南理工大学