发明名称 一种车用发动机转矩估计方法
摘要 本发明公开了一种车用发动机转矩估计方法,包括以下步骤:先离线标定放热速率中的参数和转矩损失中的系数;再在线进行转矩估计,转矩估计是先确定发动机的指示转矩Ti,再确定发动机的损失转矩T损,指示转矩与发动机损失转矩之差即为发动机的有效转矩估计值。本发明根据放热速率能计算气缸内压力,不需要气缸内压力传感器和扭矩传感器,降低了成本,克服了传感器安装困难、有效使用期短的弊端;本发明根据放热速率估计发动机转矩值,不仅能在稳态状态下估计出发动机的转矩值,而且在瞬态状态下也能估计出发动机的转矩值;同时,在满足实时性的条件下,能够获得更高的精度要求。
申请公布号 CN102518521B 申请公布日期 2013.12.25
申请号 CN201110445023.X 申请日期 2011.12.27
申请人 大连理工大学 发明人 连静;李琳辉;周雅夫;韩虎;沈晓勇;吕天赠;化玉伟
分类号 F02D29/02(2006.01)I 主分类号 F02D29/02(2006.01)I
代理机构 大连东方专利代理有限责任公司 21212 代理人 李洪福
主权项 1.一种车用发动机转矩估计方法,其特征在于:包括以下步骤:A、离线参数标定根据发动机实验对有关的参数进行标定,为在线估计提供准备工作;A1、标定放热速率中的参数使用韦伯函数估计发动机的放热速率,韦伯函数的形式为:<img file="FDA00003585205800011.GIF" wi="1891" he="218" />式中,x为进入气缸内已燃的燃油量与总的燃油量的百分比;<img file="FDA00003585205800012.GIF" wi="69" he="77" />为发动机的曲轴转角,单位为°;<img file="FDA00003585205800013.GIF" wi="78" he="80" />为燃烧起始角,单位为°;<img file="FDA00003585205800014.GIF" wi="91" he="90" />为以曲轴转角表示的燃烧持续期,单位为°;m为燃烧品质指数,无量纲;Q为每循环实际燃油燃烧的放热量,单位为kJ;Q<sub>fuel</sub>为每循环进入气缸总的燃油完全燃烧放出的热量,单位为kJ;所述的循环是发动机工作的一个完整的周期;Q<sub>fuel</sub>=m<sub>f</sub>×H<sub>μ</sub>    (2)式中,m<sub>f</sub>为每循环的喷油量,单位为kg;H<sub>μ</sub>为燃油的低热值,单位为kJ/kg;每循环的喷油量m<sub>f</sub>由公式(3)确定:m<sub>f</sub>=δ·p<sub>f</sub>·t<sub>f</sub>(3)式中,δ为喷油量系数,为定值,单位为kg/kP<sub>a</sub>·s;p<sub>f</sub>为喷油压力,单位为kP<sub>a</sub>;t<sub>f</sub>为喷油时间,单位为s;p<sub>f</sub>、t<sub>f</sub>由发动机控制单元得到;公式(1)中需要确定的参数有m、<img file="FDA00003585205800015.GIF" wi="125" he="84" /><img file="FDA00003585205800016.GIF" wi="125" he="93" />其中,对于汽油机,燃烧起始角<img file="FDA00003585205800017.GIF" wi="84" he="80" />为点火提前角,对于柴油机,燃烧起始角<img file="FDA00003585205800018.GIF" wi="88" he="83" />为喷油提前角,无论是汽油机还是柴油机,<img file="FDA00003585205800019.GIF" wi="82" he="85" />均可通过发动机的控制单元得到;标定参数m、<img file="FDA000035852058000110.GIF" wi="95" he="91" />的具体步骤如下:A11、根据发动机实验,测出气缸内指示压力p和曲轴转角<img file="FDA000035852058000111.GIF" wi="98" he="82" />绘制成发动机<img file="FDA000035852058000112.GIF" wi="162" he="77" />示功图;A12、根据公式(2)计算出Q<sub>fuel</sub>,然后根据下式(4)计算Q;接着,根据<img file="FDA00003585205800021.GIF" wi="202" he="162" />得到发动机燃烧特性曲线<img file="FDA00003585205800022.GIF" wi="268" he="86" /><img file="FDA00003585205800023.GIF" wi="1905" he="195" />式中,γ为可燃混合气比热容比,无量纲;V为气缸工作容积,单位为L;<img file="FDA00003585205800024.GIF" wi="61" he="93" />为气缸工作容积的变化率;<maths num="0001"><![CDATA[<math><mrow><mi>V</mi><mo>=</mo><mfrac><mi>&pi;</mi><mn>4</mn></mfrac><msup><mi>D</mi><mn>2</mn></msup><mi>S</mi><mo>&times;</mo><msup><mn>10</mn><mrow><mo>-</mo><mn>6</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mover><mi>V</mi><mo>.</mo></mover><mo>=</mo><mfrac><mi>&pi;</mi><mn>4</mn></mfrac><msup><mi>D</mi><mn>2</mn></msup><mo>&CenterDot;</mo><mfrac><mrow><mi>S</mi><mo>&CenterDot;</mo><msub><mi>&omega;</mi><mi>e</mi></msub></mrow><mn>30</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,ω<sub>e</sub>为发动机的转速,单位为rad/s;D为气缸直径,单位为mm;S为活塞行程,单位为mm;A13、对公式(1)两边取对数,得到如下方程:<img file="FDA00003585205800027.GIF" wi="1910" he="154" />由步骤A12得到的发动机燃烧特性曲线通过拟合直线得到参数m和<img file="FDA00003585205800028.GIF" wi="115" he="89" />A2、标定转矩损失中的系数发动机的转矩损失T<sub>损</sub>与发动机的转速、进排气歧管压力有关,其估计的表达式为:T<sub>损</sub>=d<sub>0</sub>+d<sub>1</sub>·ω<sub>e</sub>+d<sub>2</sub>·ω<sub>e</sub><sup>2</sup>+d<sub>3</sub>(p<sub>2</sub>-p<sub>1</sub>)+d<sub>4</sub>(p<sub>2</sub>-p<sub>1</sub>)<sup>2</sup>    (8)式中,d<sub>i</sub>(i=0,1,...,4)为待拟合的系数;p<sub>1</sub>为进气歧管内气体压力,单位为kP<sub>a</sub>;p<sub>2</sub>为排气歧管内气体压力,单位为kP<sub>a</sub>;确定系数d<sub>i</sub>的具体步骤如下:A21、让发动机在给定的工况下运转,当冷却水和机油温度达到正常值时,立即切断供油,同时将电力测功机转换为电动机,以给定转速ω<sub>e</sub>倒拖发动机,测出此时的进气歧管压力p<sub>1</sub>和排气歧管压力p<sub>2</sub>,此时的电力测功机转矩值为T<sub>损</sub>;在运转过程中维持冷却水和机油温度不变;A22、用步骤A21同样的方法获得一系列的发动机转速ω<sub>e</sub>、进气歧管压力p<sub>1</sub>、排气歧管压力p<sub>2</sub>、电力测功机的转矩值T<sub>损</sub>;A23、将步骤A21-A22获得的实验数据用matlab软件进行拟合,得到d<sub>i</sub>的值;B、在线转矩估计通过计算发动机指示转矩T<sub>i</sub>和转矩损失T<sub>损</sub>,从而得到发动机的有效转矩估计值T<sub>e</sub>;B1、确定发动机的指示转矩T<sub>i</sub>根据公式(1)得到放热速率计算公式,其表达式为:<img file="FDA00003585205800031.GIF" wi="1858" he="218" />式中的参数在步骤A中已标定完,根据获得的喷油压力p<sub>f</sub>、喷油时间t<sub>f</sub>,发动机转速ω<sub>e</sub>得到放热速率<img file="FDA00003585205800032.GIF" wi="110" he="140" />根据公式(4)就可以计算发动机缸内指示压力p;根据公式(10)计算出发动机的指示转矩T<sub>i</sub>;<maths num="0003"><![CDATA[<math><mrow><msub><mi>T</mi><mi>i</mi></msub><mo>=</mo><mi>p</mi><mfrac><mrow><mi>i</mi><mo>&times;</mo><mi>V</mi></mrow><mrow><mi>&pi;</mi><mo>&times;</mo><mi>&tau;</mi></mrow></mfrac><mo>&times;</mo><msup><mn>10</mn><mn>3</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,τ为发动机的行程数;i为发动机的气缸数;B2、确定发动机的损失转矩T<sub>损</sub>根据由发动机控制单元读取的进气歧管压力p<sub>1</sub>、排气歧管压力p<sub>2</sub>、发动机转速ω<sub>e</sub>,利用公式(8)得到发动机的损失转矩T<sub>损</sub>;B3、确定发动机的有效转矩发动机的有效转矩估计值为指示转矩与发动机损失转矩之差,即T<sub>e</sub>=T<sub>i</sub>-T<sub>损</sub>    (11)所述的输入量p<sub>f</sub>、t<sub>f</sub>、ω<sub>e</sub>、p<sub>1</sub>、p<sub>2</sub>均通过读取发动机控制单元得到,从而得到T<sub>i</sub>和T<sub>损</sub>的值,进而得到发动机的有效转矩估计值T<sub>e</sub>。
地址 116024 辽宁省大连市甘井子区凌工路2号