发明名称 一种微耕机汽油发动机空燃比优化控制方法
摘要 一种微耕机汽油发动机空燃比自动调节方法,通过点击控制计算机组态界面上的按钮启动事先编写好的空燃比优化控制的程序,以微耕机汽油发动机空燃比动力学模型为基础,计算与空燃比目标值对应的状态向量目标值和发动机喷油量目标值,滚动优化控制技术,设计微耕机发动机空燃比优化控制器。本发明的最重要的特征是在微耕机汽油发动机偏离空燃比目标值时能实现空燃比的自动控制,使汽油发动机在稳态和瞬态运行时都能尽可能将空燃比保持在理论目标值处,提高了微耕机汽油发动机燃油效率。
申请公布号 CN106499534A 申请公布日期 2017.03.15
申请号 CN201610855700.8 申请日期 2016.09.28
申请人 浙江工业大学 发明人 何德峰;宋秀兰;余世明
分类号 F02D41/24(2006.01)I;F02D41/26(2006.01)I;F02D41/30(2006.01)I 主分类号 F02D41/24(2006.01)I
代理机构 杭州斯可睿专利事务所有限公司 33241 代理人 王利强
主权项 一种微耕机汽油发动机空燃比优化控制方法,其特征在于:所述控制方法包括如下步骤:1)、建立微耕机汽油发动机的空燃比与喷油量的连续时间二阶动力学模型,参见式(1):<maths num="0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>A</mi><mi>x</mi><mrow><mo>(</mo><mi>t</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><mi>y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>q</mi><mi>a</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mfrac></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122213980000011.GIF" wi="1070" he="205" /></maths>其中,符号t是模型式(1)的时间变量;u(t)为微耕机汽油发动机在t时刻的喷油量,单位g/s;y(t)为微耕机汽油发动机在t时刻的空燃比;列向量x(t)=[x<sub>1</sub>(t);x<sub>2</sub>(t)]称为状态向量,x<sub>1</sub>(t)是液化在进气歧管壁上的燃油蒸发质量流量,单位g/s,x<sub>2</sub>(t)是未液化部分燃油蒸发质量流量,单位g/s;对角矩阵A∈R<sup>2×2</sup>和列向量B∈R<sup>2×1</sup>分别称为状态矩阵和控制矩阵;2)、根据微耕机汽油发动机的空燃比控制的目标值为y<sub>r</sub>和当前时刻的空气流量q<sub>a</sub>,定义状态向量的目标值x<sub>r</sub>,参见式(2):<maths num="0002"><math><![CDATA[<mrow><msub><mi>y</mi><mi>r</mi></msub><mo>=</mo><mfrac><msub><mi>q</mi><mi>a</mi></msub><mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><msub><mi>x</mi><mi>r</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122213980000012.GIF" wi="982" he="118" /></maths>3)、定义汽油发动机喷油量的目标值u<sub>r</sub>,参见式(3): u<sub>r</sub>=‑(B<sup>T</sup>B)<sup>‑1</sup>B<sup>T</sup>Ax<sub>r</sub>                     (3)其中,符号T表示矩阵的转置;4)、结合式(2)和式(3),定义符号z=x‑x<sub>r</sub>和v=u‑u<sub>r</sub>,代入式(1),整理得到一个等价的空燃比连续时间动力学模型,参见式(4):<maths num="0003"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>A</mi><mi>z</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mi>B</mi><mi>v</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><mfrac><mrow><msub><mi>q</mi><mi>a</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mrow><mo>(</mo><mi>z</mi><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>x</mi><mi>r</mi></msub><mo>)</mo></mrow></mfrac></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122213980000013.GIF" wi="1086" he="205" /></maths>5)、结合等价模型式(4),定义一个输入函数,参见式(5):<maths num="0004"><math><![CDATA[<mrow><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>u</mi><mi>r</mi></msub><mo>-</mo><mfrac><mrow><msup><mi>z</mi><mi>T</mi></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>A</mi><mi>z</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&theta;</mi><mn>1</mn></msub><msqrt><mrow><msup><mrow><mo>&lsqb;</mo><msup><mi>z</mi><mi>T</mi></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>A</mi><mi>z</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>4</mn><msub><mi>&theta;</mi><mn>2</mn></msub><msup><mrow><mo>&lsqb;</mo><msup><mi>B</mi><mi>T</mi></msup><mi>z</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>4</mn></msup></mrow></msqrt></mrow><mrow><msup><mi>B</mi><mi>T</mi></msup><mi>z</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122213980000014.GIF" wi="1478" he="151" /></maths>其中,θ<sub>1</sub>和θ<sub>2</sub>是未知参数,符号T表示向量的转;6)、结合原始的模型式(1),在t时刻定义一个跟踪偏差和喷油量的性能函数,参见式(6):<maths num="0005"><math><![CDATA[<mrow><mi>J</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><mrow><mi>s</mi><mo>=</mo><mi>t</mi></mrow><mrow><mi>t</mi><mo>+</mo><msub><mi>t</mi><mi>N</mi></msub></mrow></msubsup><mo>&lsqb;</mo><msup><mrow><mo>(</mo><mi>x</mi><mo>(</mo><mi>s</mi><mo>)</mo><mo>-</mo><msub><mi>x</mi><mi>r</mi></msub><mo>)</mo></mrow><mi>T</mi></msup><mi>Q</mi><mrow><mo>(</mo><mi>x</mi><mo>(</mo><mi>s</mi><mo>)</mo><mo>-</mo><msub><mi>x</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>+</mo><mi>R</mi><msup><mrow><mo>(</mo><mi>u</mi><mo>(</mo><mi>s</mi><mo>)</mo><mo>-</mo><msub><mi>u</mi><mi>r</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo><mi>d</mi><mi>s</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122213980000015.GIF" wi="1421" he="94" /></maths>其中,t<sub>N</sub>是预测时间窗口,Q是跟踪偏差加权矩阵,R是喷油量加权矩阵,符号s是积分变量;结合输入函数(5),在线求解优化控制问题,参见式(7):<maths num="0006"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><munder><mi>min</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mn>1</mn></msub><mo>,</mo><msub><mi>&theta;</mi><mn>2</mn></msub><mo>)</mo><mo>&Element;</mo><mi>D</mi></mrow></munder><mi>J</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mrow></mtd><mtd><mrow><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mi>A</mi><mi>x</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>+</mo><mi>B</mi><mi>u</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow></mrow></mtd><mtd><mrow><mi>u</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>u</mi><mi>r</mi></msub><mo>-</mo><mfrac><mrow><msup><mi>z</mi><mi>T</mi></msup><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mi>A</mi><mi>z</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&theta;</mi><mn>1</mn></msub><msqrt><mrow><msup><mrow><mo>&lsqb;</mo><msup><mi>z</mi><mi>T</mi></msup><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mi>A</mi><mi>z</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>4</mn><msub><mi>&theta;</mi><mn>2</mn></msub><msup><mrow><mo>&lsqb;</mo><msup><mi>B</mi><mi>T</mi></msup><mi>z</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>4</mn></msup></mrow></msqrt></mrow><mrow><msup><mi>B</mi><mi>T</mi></msup><mi>z</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mrow></mrow></mtd><mtd><mrow><mi>t</mi><mo>&le;</mo><mi>s</mi><mo>&le;</mo><mi>t</mi><mo>+</mo><msub><mi>t</mi><mi>N</mi></msub></mrow></mtd></mtr></mtable></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122213980000016.GIF" wi="1486" he="398" /></maths>其中,D为参数θ<sub>1</sub>和θ<sub>2</sub>的取值区间,符号“s.t.”表示“约束”;利用数值求解算法,得到参数θ<sub>1</sub>和θ<sub>2</sub>的一组最优解θ<sub>1</sub><sup>*</sup>和θ<sub>2</sub><sup>*</sup>,代入式(5),定义发动机在t时刻喷油量的一个最优控制器,参见式(8):<maths num="0007"><math><![CDATA[<mrow><mi>u</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>u</mi><mi>r</mi></msub><mo>-</mo><mfrac><mrow><msup><mi>z</mi><mi>T</mi></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>A</mi><mi>z</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&theta;</mi><mn>1</mn><mo>*</mo></msubsup><msqrt><mrow><msup><mrow><mo>&lsqb;</mo><msup><mi>z</mi><mi>T</mi></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>A</mi><mi>z</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>4</mn><msubsup><mi>&theta;</mi><mn>2</mn><mo>*</mo></msubsup><msup><mrow><mo>&lsqb;</mo><msup><mi>B</mi><mi>T</mi></msup><mi>z</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>4</mn></msup></mrow></msqrt></mrow><mrow><msup><mi>B</mi><mi>T</mi></msup><mi>z</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122213980000021.GIF" wi="1388" he="154" /></maths>7)、在线测量微耕机汽油发动机的实际空燃比值和空气流量,根据式(7)在线求解参数θ<sub>1</sub>和θ<sub>2</sub>的最优值,获得实时计算微耕机汽油发动机的喷油量,发动机自动调节控制系统根据喷油量的计算结果实时调整进入气缸内的喷油量,如此周而复始,实现微耕机汽油发动机空燃比对目标值y<sub>r</sub>的自动调节。
地址 310014 浙江省杭州市下城区潮王路18号浙江工业大学