CN101384810B - 用于内燃机的停止位置控制装置 - Google Patents

Abstract

一种发动机模型(60)被构造成包括曲轴外围运动方程计算块(62)、摩擦模型(64)、进气压强估计模型(66)、气缸内压强估计模型(68)、燃烧波形计算块(70)、大气压强修正项计算块(72)和大气温度修正项计算块(74)。当燃烧被切断时，使用涉及由发动机模型(60)计算出的曲柄停止位置的估计值和涉及曲柄停止位置的测量数据来获知并更新摩擦模型(64)。

Description

用于内燃机的停止位置控制装置

技术领域

本发明通常涉及一种用于控制内燃机的停止位置的装置。更确切地说，本发明涉及一种用于内燃机的停止位置控制装置，该装置应用一种适于在相关汽车临时停车后自动地停止内燃机和重新启动内燃机的控制方案。

背景技术

公开号为2004-293444(专利文件1)的日本专利申请公开了一种作为传统装置的发动机起动器，其在相关汽车临时停车后执行自动地停止内燃机和重新启动内燃机的控制(经济运行控制)。为了平稳地重新启动内燃机，这种传统装置旨在通过控制即将对车辆停止燃料供应时的发动机转速来优化在内燃机自动停止期间活塞的停止位置(曲轴的停止位置)。

同样地，公开号为2005-16505(专利文件2)的日本专利申请公开了一种混合动力汽车的车辆驱动控制装置，其使用内燃机和电动机来驱动车辆。这种传统装置使用曲轴转角(实际的曲轴转角)的检测值来修正使具有被停止于期望停止位置的曲轴的内燃机停止所需的目标转速(目标燃烧切断发动机转速)。专利文件2描述了这样的事实：即使内燃机的内部摩擦或润滑油和冷却油的温度或粘度由于修正目标转速而出现了变化，所述装置仍然能够使内燃机停止在期望的停止位置。

包括上述文件在内，申请人认为以下文件是本发明的相关技术。

【专利文件1】公开号为2004-293444的日本专利申请。

【专利文件2】公开号为2005-16505的日本专利申请。

【专利文件3】公开号为平02-305342的日本专利申请。

【专利文件3】公开号为昭58-18535的日本专利申请。

发明内容

当内燃机自动停止时，曲轴的停止位置可能偏离期望位置。需要考虑诸如上述摩擦变化的多种因素。然而，在专利文件1中，没有考虑到这种引起曲轴的停止位置偏离期望位置的多个因素。另外，根据专利文件2，尽管考虑了偏离的原因，但由于是在将各种因素一起考虑的情况下来调节期望的停止位置，所以对曲轴停止位置的控制在精度上和收敛性上都被认为是低的。

本发明旨在解决上述问题，并且本发明的目的是提供一种装置，该装置应用一种适于自动地停止和重新启动内燃机的控制方案，并且该装置能够精确地估计曲轴在内燃机中的停止位置。

通过一种用于内燃机的停止位置控制装置实现上述目的，所述停止位置控制装置包括计算内燃机摩擦的摩擦模型。所述装置基于计算出的摩擦来修正曲轴的停止位置。

在本发明的第二个方案中，所述装置可适用于通过控制内燃机的燃烧即将停止时的燃烧切断发动机转速来控制曲轴的停止位置。可以设置基于包括所述摩擦的必要参数来计算曲轴的停止位置的估计值的曲柄停止位置计算部件。所述曲柄停止位置计算部件可以基于由摩擦模型计算出的摩擦来修正曲轴的停止位置的估计值。

本发明的第三个方案包括用于判断曲轴的停止位置的曲柄停止位置判断部件。如果判断为曲轴的停止位置处于必要误差范围之外，则所述装置基于曲轴的停止位置的估计值和曲轴的停止位置的测量值之间的偏差来获知摩擦模型并且更新摩擦模型。

在本发明的第四个方案中，所述装置适用于通过控制内燃机的燃烧即将停止时的燃烧切断发动机转速来控制曲轴的停止位置。可以设置用于基于包括摩擦的必要参数来计算燃烧切断发动机转速的目标值的目标燃烧切断发动机转速计算部件。所述目标燃烧切断发动机转速计算部件可以基于由摩擦模型计算出的摩擦来修正所述目标值。

本发明的第五个方案包括温差判断部件，其用于判断冷却水温度和润滑油温度之间的差值。可以设置用来计算与所述温差相关的摩擦误差的摩擦误差计算部件。如果判断为冷却水温度和润滑油温度彼此不同，则所述装置可以基于摩擦误差来修正摩擦模型。

在本发明的第六个方案中，如果判断为冷却水温度和润滑油温度彼此不同，则所述装置可以禁止对曲轴的停止位置的控制。

根据本发明的第一个方案，由于对于影响曲柄停止位置精度的各个因素而言都考虑了所述因素对曲柄停止位置精度的影响，所以能够精确地控制曲轴的停止位置。

根据本发明的第二个方案，由于能够估计出曲轴停止位置的数据，并且对于影响曲柄停止位置精度的各个因素而言都考虑了所述因素对曲柄停止位置精度的影响，所以能够实现一种执行高效算术运算的停止位置控制装置。

根据本发明的第三个方案，由于只对曲轴停止位置具有显著影响的摩擦的影响能够被获知并更新，所以能够实现对估计曲轴停止位置的精度的改进和高效算术运算。

根据本发明的第四个方案，由于能够获得停止燃烧所需的发动机转速(以下称作“燃烧切断发动机转速”)的目标值，并且对于影响曲柄停止位置精度的各个因素而言都考虑了所述因素对曲轴停止位置精度的影响，所以能够实现一种可执行高效算术运算的停止位置控制装置。

根据本发明的第五个方案，由于考虑到了摩擦受冷却水温度和润滑油温度影响的事实，所以能够更精确估计曲轴的停止位置。

根据本发明的第六个方案，由于如果判断为未精确估计摩擦则禁止对曲柄停止位置的控制，所以能够获得曲轴停止位置的高可控性。

附图说明

图1为说明将用于内燃机的停止位置控制装置应用于作为本发明的第一实施例的内燃机的构造的图解。

图2为显示配置在图1的ECU中的发动机模型的构造的方块图。

图3为显示分配给曲轴的外围组成元件的附图标记和记号的图解。

图4示出了图2的摩擦模型包括用于获取摩擦转矩TRQf的摩擦设定表(friction map)的实例。

图5是说明用作对获取气缸内压强P的记录的改进的方法的图解。

图6是使用发动机模型，为了实现曲柄停止位置的控制而由根据第一实施例的ECU执行的程序的流程图。

图7为显示图1的ECU存储的曲轴的目标停止区域的图解。

图8为说明在步骤110中执行的用于将PID修正数据并入摩擦模型的摩擦设定表上的各个映射点的值中的方法的图解。

图9为说明使用反馈控制的目标燃烧切断发动机转速的计算方法。

图10为说明在图6所示程序的步骤114中进行的处理的时间图。

图11为由根据本发明的第二实施例的ECU执行的程序的流程图。

图12为由根据本发明的第三实施例的ECU执行的程序的流程图。

具体实施方式

第一实施方式

【根据第一实施方式的系统构造】

图1为说明将用于内燃机的停止位置控制装置应用于作为本发明的第一实施例的内燃机10的构造的图解。本实施方式的系统包括内燃机10。并且，本实施方式假定内燃机10为直列四缸发动机。内燃机10的各个气缸都包含活塞12。活塞12经由连杆14与曲轴16结合。内燃机10的各个气缸还具有形成于活塞12上方的燃烧室18。进气通道20和排气通道22与燃烧室18连通。

进气通道20包括节气门24。节气门24为电控节气门，其能够独立于加速器的打开位置来控制自身的打开位置。检测节气门的打开位置TA的节气门位置传感器26靠近节气门24设置。用于将燃料喷入进气口的燃料喷射阀28相对于节气门24设置在下游。在内燃机的气缸盖上，用于各个气缸的火花塞30以火花塞30从燃烧室18的上部凸出并伸入燃烧室18中的方式被安装。进气门32和排气门34分别设置在进气口和排气口处。进气门32使燃烧室18和进气通道20之间连续或中断，而排气门34使燃烧室18和排气通道22之间连续或中断。

进气门32和排气门34分别由进气可变气门正时(VVT)控制机构36和排气可变气门正时(VVT)控制机构38驱动。进气可变气门正时(VVT)控制机构36和排气可变气门正时(VVT)控制机构38分别能够与曲轴16的旋转同步地打开和关闭进气门32和排气门34。VVT控制机构还能够改变进气门和排气门的各气门开启特性(诸如气门开启正时，工作角和气门升程量)。

内燃机10具有靠近曲轴16的曲轴转角传感器40。每当曲轴16旋转一个所需角度时，曲轴转角传感器40就会反转一个高水平(“Hi”)输出状态和一个低水平(“Lo”)输出状态。能够根据曲轴转角传感器40的输出来检测曲轴16的旋转位置及其转速(即，发动机转速，“Ne”)。内燃机10还具有靠近进气凸轮轴的凸轮转角传感器42。凸轮转角传感器42实质上具有与曲轴转角传感器40相同的构造。能够根据凸轮转角传感器42的输出来检测进气凸轮轴的旋转位置(提前值)和其它因素。

图1所示的系统还包括电子控制模块(ECU)50。除了上述的各种传感器之外，用于检测排气通道22的内部排气空/燃比的空燃比传感器52和用于检测内燃机10中的冷却水温度的水温传感器54电性地连接到ECU 50上。此外，上述各种执行器电性地连接到ECU 50上。ECU50能够基于各种传感器的输出且基于使用实际上配置在ECU 50中的发动机模型60所计算出的结果来控制内燃机10的运行状态。

【发动机模型概述】

图2为显示配置在图1的ECU 50中的发动机模型60的构造的方块图。如图2所示，发动机模型60包括曲轴外围运动方程计算块62、摩擦模型64、进气压强估计模型66、气缸内压强估计模型68、燃烧波形计算块70、大气压强修正项计算块72和大气温度修正项计算块74。下面给出对这些构件的详细结构描述。

(1)关于曲轴外围运动方程计算块的描述

曲轴外围运动方程计算块62计算曲轴转角θ和发动机转速“Ne”的估计值，即，曲柄角转速dθ/dt。曲轴外围运动方程计算块62接收来自气缸内压强估计模型68或燃烧波形计算块70的内燃机10的气缸内压强P。在计算开始期间内，曲轴外围运动方程计算块62还接收初始曲轴转角“θ0”和初始发动机转速“Ne0”。

为了分别补偿实际曲轴转角θ和实际发动机转速“Ne”的偏差，估计曲轴转角θ和估计发动机转速“Ne”都受到图2中所示的PID控制器76的反馈控制，估计曲轴转角θ和估计发动机转速“Ne”二者都是由曲轴外围运动方程计算块62计算出来的。通过摩擦模型64将涉及内燃机10的内部摩擦的影响并入到由曲轴外围运动方程计算块62计算出的结果中。

接下来，将在下面描述在曲轴外围运动方程计算块62中执行的计算的进一步的细节。

图3为显示分配给曲轴16的外围组成元件的附图标记和记号的图解。如图3所示，本实施例假定承受气缸内压强P的活塞12的上端具有表面面积A。同样地，连杆14的长度被称作L，而曲柄的旋转半径为“r”。另外，在连接连杆14的活塞安装点与曲轴16的轴心的虚线(气缸轴线)和连杆14的轴线之间形成的角被称作f(下面称作“连杆转角f”)。此外，在气缸轴线和曲柄销17的轴线之间形成的角被称作θ。

就气缸间的曲轴转角而言，由于具有四个气缸的内燃机10具有180°CA的相差，气缸间的曲轴转角的关系可以定义为方程式(1a)。另外，由于气缸的曲柄角转速dθ/dt是通过将各个气缸的曲轴转角θ相对于时间进行微分而获得的，所以各个角能够表达为方程式(1b)。

θ₁＝θ，θ₂＝θ+π，θ₃＝θ+2π，θ₄＝θ+3π    ........(1a)

$\stackrel{\·}{\mathrm{\θ}}={\stackrel{\·}{\mathrm{\θ}}}_1,$ $\stackrel{\·}{\mathrm{\θ}}={\stackrel{\·}{\mathrm{\θ}}}_2,$ $\stackrel{\·}{\mathrm{\θ}}={\stackrel{\·}{\mathrm{\θ}}}_3,$ $\stackrel{\·}{\mathrm{\θ}}={\stackrel{\·}{\mathrm{\θ}}}_4$ ........(1b)

$(\stackrel{\·}{\mathrm{\θ}}=\frac{\mathrm{d\θ}}{\mathrm{dt}})$

在以上方程式(1a)和方程式(1b)中分配给曲轴转角θ和曲柄角转速dθ/dt的参考标记1至4依照内燃机10的所需爆发次序对应于各个气缸中的燃烧次序。同样地，在以下的表达式中，参考标记1至4可以由“i”来代表。

另外，在图3所示的活塞曲柄机构中，曲轴转角θi和连杆转角fi具有由下面的方程式(2)表达的关系。

$\sin \left({\mathrm{\φ}}_i\right)=\frac{r}{L}\sin \left({\mathrm{\θ}}_i\right),\cos \left({\mathrm{\φ}}_i\right)=\sqrt{1-\left(\frac{r}{L}\right){\sin }^2\left({\mathrm{\θ}}_i\right)},$

${\stackrel{\·}{X}}_i=r\·\sin \left({\mathrm{\θ}}_i\right)\{1+\frac{\frac{r}{L}\cos \left({\mathrm{\θ}}_i\right)}{\sqrt{1-{\left(\frac{r}{L}\right)}^2{\sin }^2\left({\mathrm{\θ}}_i\right)}}\}{\stackrel{\·}{\mathrm{\θ}}}_i---\left(2\right)$

$(\stackrel{\·}{X}i=\mathrm{dXi}/\mathrm{dt})$

其中“dXi/dt”表示活塞的速度。

此外，曲轴外围的总动能可以表达为方程式(3)。推导方程式(3)使得方程式(3)的各项中的不同参数由1/2(dθ/dt)²的系数来表示。这里，由此获得的系数表示曲轴转角θ的函数f(θ)。

$T=\frac{1}{2}(I_k+I_{f1}+I_{\mathrm{mi}}){\stackrel{\·}{\mathrm{\θ}}}^2+\underset{i=1}{\overset{4}{\mathrm{\Σ}}}\frac{1}{2}(m_p+m_c){\stackrel{\·}{X}}_i^2+\underset{i=1}{\overset{4}{\mathrm{\Σ}}}\frac{1}{2}{I}_c{\stackrel{\·}{\mathrm{\φ}}}_i^2$

$=\frac{1}{2}[(I_k+I_{f1}+I_{\mathrm{mi}})+(m_p+m_c)r^2\·\underset{i=1}{\overset{4}{\mathrm{\Σ}}}{\sin }^2\left({\mathrm{\θ}}_i\right).$

${\{1+\frac{\frac{r}{L}\cos \left({\mathrm{\θ}}_i\right)}{\sqrt{1-{\left(\frac{r}{L}\right)}^2{\sin }^2\left({\mathrm{\θ}}_i\right)}}\}}^2+I_c{\left(\frac{r}{L}\right)}^2\·\underset{i=1}{\overset{4}{\mathrm{\Σ}}}\frac{{\cos }^2\left({\mathrm{\θ}}_i\right)}{1-{\left(\frac{r}{L}\right)}^2{\sin }^2\left({\mathrm{\θ}}_i\right)}]\·{\stackrel{\·}{\mathrm{\θ}}}^2$

$=\frac{1}{2}\·f\left(\mathrm{\θ}\right)\·{\stackrel{\·}{\mathrm{\θ}}}^2$ ........(3)

方程式(3)右侧的第一项相当于与曲轴16的旋转运动有关的动能，右侧的第二项相当于与活塞12和连杆14的直线运动有关的动能，而右侧的第三项相当于与连杆14的旋转运动有关的动能。同样地，方程式(3)中的“Ik”表示曲轴16围绕其旋转轴的惯性力矩，而“If1”为飞轮围绕其旋转轴的惯性力矩。另外，“Imi”表示包括传动系统(即，变速装置、驱动轴、轮胎以及与内燃机10结合的其它元件)的旋转构件围绕旋转构件的旋转轴的惯性力矩，而“Ic”为连杆14的惯性力矩。此外，“mp”表示活塞12的位移，而“mc”表示连杆14的位移。

接下来，拉格朗日算子L被定义为系统的总动能T和势能U之间的偏差，如方程式(4a)。如果作用在曲轴16上的输入转矩定为TRQ，则运动的拉格朗日方程式能够用来表达拉格朗日算子L、曲轴转角θ和输入转矩TRQ之间的关系，如方程式(4b)。

L＝T-U                  ........(4a)

$\frac{d}{\mathrm{dt}}\frac{\∂L}{\∂\stackrel{\·}{\mathrm{\θ}}}-\frac{\∂L}{\∂\mathrm{\θ}}=\mathrm{TRQ}$ ........(4b)

$\frac{\∂L}{\∂\stackrel{\·}{\mathrm{\θ}}}=f\left(\mathrm{\θ}\right)\stackrel{\·}{\mathrm{\θ}},$ $\frac{d}{\mathrm{dt}}\frac{\∂L}{\∂\stackrel{\·}{\mathrm{\θ}}}=\frac{d}{\mathrm{dt}}\frac{\∂f\left(\mathrm{\θ}\right)}{\∂\mathrm{\θ}}{\stackrel{\·}{\mathrm{\θ}}}^2+f\left(\mathrm{\θ}\right)\stackrel{\·\·}{\mathrm{\θ}}$ ........(4c)

$\frac{\∂L}{\∂\mathrm{\θ}}=\frac{1}{2}\frac{\∂f\left(\mathrm{\θ}\right)}{\∂\mathrm{\θ}}{\stackrel{\·}{\mathrm{\θ}}}^2$ ........(4d)

$\stackrel{\·}{\·\·}\frac{d}{\mathrm{dt}}\frac{\∂L}{\∂\stackrel{\·}{\mathrm{\θ}}}-\frac{\∂L}{\∂\mathrm{\θ}}=\mathrm{TRQ}\⇔f\left(\mathrm{\θ}\right)\stackrel{\·\·}{\mathrm{\θ}}+\frac{1}{2}\frac{\∂f\left(\mathrm{\θ}\right)}{\∂\mathrm{\θ}}{\stackrel{\·}{\mathrm{\θ}}}^2=\mathrm{TRQ}$ ........(4e)

在方程式(4a)中，由于势能U的影响相比动能T的影响显得微不足道，所以前者的影响可以忽略。通过对方程式(3)相对于曲柄角转速(dθ/dt)进行偏微分而获得的值进行时间微分，方程式(4b)左侧的第一项可以表达为曲轴转角θ的函数，如方程式4C。同样地，通过对方程式(3)相对于曲轴转角θ进行偏微分，方程式4B左侧的第二项可以表达为曲轴转角θ的函数，如方程式(4d)。

因此，方程式(4b)可以改写为方程式(4e)，由此能够得到曲轴转角θ和输入转矩TRQ之间的关系。同样地，输入转矩TRQ在这里被定义为具有三个参数的值，如方程式(5)。

TRQ＝TRQ_e-TRQ_L-TRQ_f                 ......(5)

在上述方程式(5)中，TRQe为发动机中产生的转矩。更确切地说，TRQe是来自承受气体压强(气缸内压强P)的活塞12并作用在曲轴16上的转矩。同样地，TRQL为负荷转矩，其作为根据装有内燃机10的车辆的具体特性而变化的已知值存储在ECU 50中。另外，TRQf为摩擦转矩，也就是与诸如活塞12和曲轴16的滑动构件中的总摩擦损失对应的转矩。摩擦转矩TRQf是从摩擦模型64得到的值。

接下来，可以依照方程式(6a)至(6c)来计算发动机转矩TRQe。也就是说，首先，由于气缸内压强P可以表达为作用在活塞12上端的力PA的连杆14的轴向线性方向分力，所以力Fc作用在连杆14上。方程式(6a)表达了力Fc。如图3所示，由于形成于连杆14的轴线和曲柄销17的移动路线的切线之间的角a为{π/2-(Φ+θ)}，由于气缸内压强P而作用在曲柄销17的移动路线的切线方向上的力Fk可以使用施加在连杆14上的力Fc来表达。方程式(6b)表达了力Fk。发动机转矩TRQe是作用在曲柄销17的移动路线的切线方向上的力Fk和曲柄的旋转半径“r”的乘积。因此，方程式(6a)和(6b)可以用来将Fk表达为方程式(6c)。

F_c＝P·Acos(Φ)     ........(6a)

F_k＝F_csin(Φ+θ)    ........(6b)

$\stackrel{\·}{\·\·}{\mathrm{TRQ}}_e=F_k\·r=P\·A\·r\·\cos \left(\mathrm{\φ}\right)\sin (\mathrm{\φ}+\mathrm{\θ})$

$=P\·A\·r\·\left[\right\{1-{\left(\frac{r}{L}\right)}^2{\sin }^2\left(\mathrm{\θ}\right)\}+\frac{r}{L}\cos \left(\mathrm{\θ}\right)]\sin \left(\mathrm{\θ}\right)$ ........(6c)

根据曲轴外围运动方程计算块62的构造，通过利用气缸内压强估计模型68或燃烧波形计算块70来获取气缸内压强P，可以依照方程式(6c)和(5)获得输入转矩TRQ。这使得能够通过解答方程式(4e)来获得曲轴转角θ和曲柄角转速dθ/dt。

(2)关于摩擦模型的描述

图4示出了图2的摩擦模型64包括用于获取摩擦转矩TRQf的摩擦设定表的实例。在图4的曲线图中，摩擦转矩TRQf是根据发动机转速“Ne”和发动机冷却水温度的之间关系定义的。曲线图中显示的这种特性是通过预先进行的实验等确定的，并且随着发动机冷却水温度的下降，摩擦转矩TRQf趋向于增加。

为了降低ECU 50的计算负荷，摩擦模型64包含了上述的曲线图(摩擦设定表)的情况已经在本实施例的上文中进行了描述。然而，摩擦模型的构造不局限于以上描述或受以上描述的限制，并且可以使用诸如方程式(7)的关系表达式。方程式(7)被配置为使得摩擦转矩TRQf成为以内燃机10的发动机转速“Ne”和润滑油的动力粘度υ为参数的函数。

TRQ_f＝C₁·Ne²+C₂·υ+C₃             ......(7)

其中，C1、C2、C3是实验所用的或在实验期间所用的各个系数。

(3)关于进气压强估计模型的描述

进气压强估计模型66具有用于估计进气压强的进气压强设定表(未示出)。该进气压强设定表是根据进气压强、节气门角度“TA”、发动机转速“Ne”以及进气门和排气门的气门正时“VVT”之间的关系定义的。根据这种进气压强估计模型的构造，在最小化ECU 50的计算负荷的同时，能够获取进气压强。在详细计算进气压强时，可以使用估计流经节气门24的空气流速的节气门模型和估计流经进气门32的空气的流速(即，气缸内进气流速)的气门模型，以取代使用上述的这种进气压强设定表，来配置进气压强估计模型。

(4)气缸内压强估计模型

气缸内压强估计模型68用来计算在不执行燃烧的情形下的气缸内压强P。在气缸内压强估计模型68中，在内燃机10的各个行程中的气缸内压强P是使用下面示出的方程式(8a)至(8d)来计算的。也就是说，如方程式(8a)，首先从由进气压强估计模型66的进气压强设定表得来的气缸内压强设定表数据“Pmap”中获得在进气行程期间的气缸内压强P。

进气行程    P＝P_map(Ne，VVT，TA)          ......(8a)

压缩行程 $P={\left(\frac{V_{\mathrm{bdc}}}{V}\right)}^{\mathrm{\κ}}\·P_{\mathrm{map}}$ ......(8b)

膨胀冲程 $P={\left(\frac{V_{\mathrm{tdc}}}{V}\right)}^{\mathrm{\κ}}\·P_C$ ......(8c)

排气行程    P＝P_ex≈P_air                  ......(8d)

接下来，依照气体的可逆绝热变化的方程式，可以将压缩行程期间的气缸内压强P表达为方程式(8b)。

在以上方程式(8b)中，V_BDC表示当活塞12出现在进气下止点时存在的行程容积V，而“k”表示比热比(specific heat ratio)。

如同压缩行程，在膨胀冲程期间的气缸内压强P可以表达为方程式(8c)。

在以上方程式(8c)中，V_TDC表示当活塞12出现在压缩上止点时存在的行程容积V，而Pc表示在压缩行程结束期间获得的气缸内压强。

同样地，如方程式(8d)所示，假定在排气行程期间的气缸内压强P是排气通道22中的内部压强“Pex”。压强“Pex”可以看作大体等于大气压强“Pair”。因此，估计模型使用大气压强“Pair”作为在排气行程期间获得的气缸内压强P。

(5)关于燃烧波形计算块的描述

燃烧波形计算块70用于计算在从压缩行程的中间阶段至膨胀冲程的中间阶段的燃烧期间内的气缸内压强(燃烧压强)P。燃烧波形计算块70使用方程式(9a)计算燃烧压强P的估计值，所述方程式(9a)是基于Weibe函数和下面描述的方程式(10)的关系表达式。

$\frac{\mathrm{dQ}}{\mathrm{d\θ}}=a\·\frac{k\·Q}{{\mathrm{\θ}}_p}\·(m+1)\·{\left(\frac{\mathrm{\θ}-{\mathrm{\θ}}_b}{{\mathrm{\θ}}_p}\right)}^m\·\exp \{-a\·{\left(\frac{\mathrm{\θ}-{\mathrm{\θ}}_b}{{\mathrm{\θ}}_p}\right)}^{m+1}\}$ ........(9a)

$=-a\·(m+1)\·{\left(\frac{\mathrm{\θ}-{\mathrm{\θ}}_b}{{\mathrm{\θ}}_p}\right)}^m\·\exp \{-a\·{\left(\frac{\mathrm{\θ}-{\mathrm{\θ}}_b}{{\mathrm{\θ}}_p}\right)}^{m+1}\}$ ........(9b)

上述的方程式(9a)可以改写为如下：

$\frac{\mathrm{dQ}}{\mathrm{d\θ}}=-\frac{k\·Q}{{\mathrm{\θ}}_p}\·\frac{d}{\mathrm{d\θ}}\left(\exp \right\{-a\·{\left(\frac{\mathrm{\θ}-{\mathrm{\θ}}_b}{{\mathrm{\θ}}_p}\right)}^{m+1}\left\}\right)=-\frac{k\·Q}{{\mathrm{\θ}}_p}\·\frac{\mathrm{dg}\left(\mathrm{\θ}\right)}{\mathrm{d\θ}}$

$\⇔\frac{\mathrm{dQ}}{Q}\·\frac{1}{\mathrm{d\θ}}=-\frac{k}{{\mathrm{\θ}}_p}\·\frac{\mathrm{dg}\left(\mathrm{\θ}\right)}{\mathrm{d\θ}}$ ........(9c)

在方程式(9c)的两侧对θ进行积分，生成如下结果：

$\∫\frac{1}{Q}\·\frac{\mathrm{dQ}}{\mathrm{d\θ}}\mathrm{d\θ}=-\frac{k}{{\mathrm{\θ}}_p}\·\∫\frac{\mathrm{dg}\left(\mathrm{\θ}\right)}{\mathrm{d\θ}}\mathrm{d\θ}$

$\⇔\∫\frac{1}{Q}\mathrm{d\θ}=-\frac{k}{{\mathrm{\θ}}_p}\·\∫\mathrm{dg}\left(\mathrm{\θ}\right)$

$\⇔\log Q+C_2=-\frac{k}{{\mathrm{\θ}}_p}\·g\left(\mathrm{\θ}\right)+C_1$

$\log Q=-\frac{k}{{\mathrm{\θ}}_p}\·g\left(\mathrm{\θ}\right)+C$

(其中C＝C1-C2：C、C1、C2都是积分常数)

$Q=\exp (C-\frac{k}{{\mathrm{\θ}}_p}\·g\left(\mathrm{\θ}\right))=\exp [C-\frac{k}{{\mathrm{\θ}}_p}\·\exp \{-a\·{\left(\frac{\mathrm{\θ}-{\mathrm{\θ}}_b}{{\mathrm{\θ}}_p}\right)}^{m+1}\}]$

........(9d)

更确切地说，燃烧波形计算块70首先使用方程式(9a)来计算与当前曲轴转角θ相关的热释放率dQ/dθ。

上述方程式(9a)中的参数“m”、“k”、“θb”和“a”分别表示形状系数、燃烧效率、点火延迟周期和燃烧速度。这里，燃烧速度“a”为固定值6.9。各个参数都使用预先适用的值。同样地，Q表示热释放量。

为了使用方程式(9a)计算热释放率dQ/dθ，需要计算热释放量Q。热释放量Q能够通过解微分方程式(9a)来计算。因此，在方程式(9b)中，与方程式(9a)中的Weibe函数等价的部分首先由“g(θ)”替代。这使得方程式(9a)能够表达为方程式(9c)。接下来，在方程式(9c)的两侧对曲轴转角θ进行积分后，推导方程式(9c)使得热释放量Q能够表达为方程式(9d)。之后，通过将依照方程式(9d)计算出的热释放量Q赋值给方程式(9a)，能够计算出热释放率dQ/dθ。

通过使用基于能量守恒定律的关系表达式，能够将热释放率dQ/dθ和气缸内压强(燃烧压强)P表达为方程式(10)。因而，在将依照方程式(9a)计算出的热释放率dQ/dθ赋值给方程式(10)后，可以通过解方程式(10)来计算燃烧压强P。

$\frac{\mathrm{dQ}}{\mathrm{d\θ}}=\frac{1}{\mathrm{\κ}-1}\·(V\·\frac{\mathrm{dP}}{\mathrm{d\θ}}+\mathrm{\κ}\·P\·\frac{\mathrm{dV}}{\mathrm{d\θ}})$ ......(10)

使用上述的气缸内压强估计模型68和燃烧波形计算块70使得，无论是否进行燃烧，都能够获取到内燃机10的气缸内压强P的记录。也就是说，通过使用气缸内压强估计模型68来计算非燃烧状态下的气缸内压强P，以及通过使用燃烧波形计算块70来计算燃烧状态下的气缸内压强P，都能够获取到内燃机10的气缸内压强P的记录。

获取内燃机10的气缸内压强P的记录的方法不局限于如上所述的方法，并且可以是例如参照图5所示的方法。

图5是说明用作上述的改进的方法的图解。所述方法不是通过使用上述的方程式(9a)和(10)来计算用于每个所需曲轴转角θ的燃烧压强P。替代性地，所述方法预先使用方程式(9a)和(10)来计算如图5A所示的燃烧模式，即，只是由燃烧引起的气缸内压强P的波形的变化(换句话说，是由燃烧引起的压强增加)。

更确切地说，是保存了一个设定表，其中点火延迟周期、燃烧周期和“ΔPmax”(燃烧状态下的最大压强“Pmax”和非燃烧状态下的最大压强“Pmax0”之间的偏差)是用于确定上述燃烧模式的三个参数。所述设定表定义了发动机转速“Ne”、空气填充速率KL、进气门和排气门的气门正时VVT和点火正时之间的关系。为了计算与由燃烧引起的压强增加所对应的波形，作为通过合并诸如二次函数的简化函数而近似波形，与所述近似波形相关的不同系数与上述发动机转速“Ne”形成映射关系。接下来，如图5B所示，参照这种映射得到的燃烧诱导的压强增加的波形被添加到由气缸内压强估计模型68计算出的气缸内压强P的值上。这样，就得到了燃烧压强P。

(6)关于大气压强修正项计算块的描述

大气压强修正项计算块72包括估计吸入到气缸中的气缸充气量“Mc”的模型(下面，该模型被称为“空气模型”)。该空气模型按照以下方程式(11)来计算气缸充气量“Mc”。

Mc＝aPm-b            ......(11)

在以上的方程式(11)中，符号“a”和“b”表示根据诸如发动机转速“Ne”和气门正时VVT的特定运行参数的适合的系数。同样地，“Pm”是进气压强。例如，可以使用由进气压强估计模型66计算出的值。

大气压强修正项计算块72同样地包括估计吸入到气缸中的燃料量“Fc”的模型(下面，该模型被称为“燃料模型”)。考虑到燃料在从燃料喷射阀28喷出之后的状态(behavior)，即，考虑到部分喷出燃料粘着在例如进气口的内壁上以及所述粘着的燃料随后气化的情况，如果在第k个循环的燃料喷射开始时粘着在壁部上的燃料量为“fw(k)”，而在第k个循环中的实际燃料喷射量为“fi(k)”，则以下的方程式(12a)和(12b)可以用来表达在第k个循环后粘着在壁部的燃料量，如“fw(k+1)”，以及在第k个循环期间进入气缸内的燃料量，如“fc”。

f_w(k+1)＝P(k)·f_w(k)+R(k)·f_i(k)             ........(12a)

f_c(k)＝(1-P(k))·f_w(k)(1-R(k))·f_i(k)        ........(12b)

在以上的方程式(12)中，P表示粘着率，确切地说，该粘着率是粘着在例如进气口的内壁上的燃料量相对于燃料喷射量“fi”的比率。同样地，R表示残留率，确切地说，该残留率是粘着燃料量相对于在进气行程之后粘着在例如壁部上的燃料量的比率。使用方程式(12)能够以粘着率P和残留率R为参数来计算每个循环的上述燃料量“fc”。

因此，能够使用空气模型和燃料模型的计算结果来计算估计的空/燃(A/F)比。接下来，大气压强修正项计算块72计算在估计A/F比和A/F比的测量值之间的稳定状态偏差，所述A/F比的测量值是在考虑从燃料喷射完成至燃烧的燃料到达空燃比传感器52的燃料传输延迟时测量到的。上述稳定状态偏差是气缸充气量“Mc”的误差。如果稳定状态偏差很大，则大气压强修正项计算块72认为大气压强不合适，并且计算大气压强修正系数“kairp”。更确切地说，空气模型逆向计算进气压强“Pm”，然后基于计算出的进气压强“Pm”来计算大气压强修正系数“kairp”作为与标准大气压强“Pa0”相关的修正比。将计算出的大气压强修正系数“kairp”用于上述的进气压强估计模型66和气缸内压强估计模型68以分别修正进气压强“Pmap”和排气压强(大气压强“Pair”)。

(7)关于大气温度修正项计算块的描述

通过将排气行程中存在的行程容积V、残留气体质量“m”(从存在于排气上止点处的余隙容积Vc计算出来)、残留气体(已燃气体)的气体常数R和测量到的大气温度“Tair”赋值给理想气体的状态方程式，大气温度修正项计算块74计算气缸内压强“Pth”。之后，计算在计算出的气缸内压强“Pth”和由气缸内压强估计模型68计算出的气缸内压强P之间的偏差。如果偏差很大，则从上述偏差来计算修正系数。将所述计算出的修正系数用于上述进气压强估计模型66以修正进气压强“Pmap”。

【使用发动机模型的曲轴停止位置控制】

当装有内燃机的车辆即将临时停车时，可以执行用于自动停止并重新启动内燃机的控制(经济运行控制)。对于由内燃机和电动机驱动的混合动力汽车来说，也可以在车辆系统的加电(powerup)期间内(同样地，可以在车辆的运行状态期间)执行这种控制。在本说明书中，在上述后一种情况下执行的控制同样被称为“经济运行控制”。

为了在上述经济运行控制期间平稳地起动内燃机，要求在内燃机自动停止期间的曲轴16的实际停止位置(活塞12的实际停止位置)应该被精确地控制在期望的停止位置上。本实施例的系统假定将上述的发动机模型60用作估计在经济运行控制期间的曲轴16的停止位置的停止位置估计模型。根据发动机模型60，在内燃机10的自动停止期间内的曲轴16的停止位置可以通过计算曲柄角转速dθ/dt为零时的曲轴转角θ的值来获得。另外，可以通过执行发动机模型60的逆向计算来计算在内燃机10以曲轴16的期望停止位置被停止时的燃烧切断发动机转速。在此说明书中，曲轴16的停止位置可以简称为“曲柄停止位置”或“曲柄的停止位置”。

已知的设计技术为：当内燃机自动停止时，用于停止燃料供应所需的发动机转速(燃烧切断发动机转速)被控制，使得曲柄的停止位置被设定在目标位置上。使用这种技术时，曲柄的停止位置可能偏差其目标停止位置。在这种情况下，这种偏差很可能是由诸如发动机内部滑动部件的摩擦转矩的变化的多个因素所引起的。如果引起停止位置误差的这些原因得不到检查或澄清并且考虑所有这些原因来获知曲柄的目标停止位置的话，则这种获知的精度和会聚性会降低。另外，当通过控制燃烧切断发动机转速来控制曲柄的停止位置时，摩擦的影响会变得比其它原因的影响更明显。

因此，在本实施例的系统中，如果在曲柄的实际停止位置和其目标停止位置之间存在明显偏差，则只有摩擦转矩被用作待获知的参数，并且基于曲柄停止位置的上述偏差来获知包含在前述发动机模型60中的摩擦模型64。在本实施例的系统中，其它原因都适合使用传感器输出。因此，只有对曲柄停止位置的估计产生显著影响的摩擦转矩的作用能够被获知，并且系统不受除摩擦转矩之外的任何原因的影响。因此，能够提高对曲柄的停止位置的估计精度。

接下来，参照图6详细描述根据第一实施例的处理。

图6是使用发动机模型，为了实现曲柄停止位置的控制而由根据第一实施例的ECU 50执行的程序的流程图。

(涉及步骤100的处理)

在图6的程序中，首先在步骤100中判断经济运行模式标记是否设定为开启状态。如果ECU 50判断为所需的经济运行控制执行条件成立，更确切地说，如果车速为零并且内燃机10处于怠速状态，则ECT将经济运行模式标记设为开启状态。

(涉及步骤102的处理)

如果在上述步骤100中判断为经济运行模式标记未设定为开启状态，则在执行步骤102时获知燃烧压强P。在步骤102中，获取一个循环中的燃烧压强P的最大值“Pmax”，并将通过对几百个循环的“Pmax”值平均后而得到的值进行保存。将作为几百个循环的“Pmax”值中的最小值和最大值的最小燃烧压强和最大燃烧压强也进行保存。本实施例假定所述平均值、最小值和最大值是在当内燃机10处于怠速状态时，例如以10rpm的间隔，以所需的定时进行保存的。由于在内燃机10的怠速状态期间的燃烧程度的变化会相对明显，因此保存最小和最大燃烧压强使得能够在随后的步骤中将所述这些值作为燃烧不一致性时的信息。

在步骤102中，在用于由燃烧波形计算块70计算燃烧压强P的上述方程式(9a)中，Weibe函数的诸如“θb”和“θp”的不同的参数被修改(适用)，从而在通过基于曲轴转角传感器40的输出而测量获得的发动机转速“Ne”值和由曲轴外围运动方程计算块62算术估计出的发动机转速“Ne”值之间的没有偏差。

更确切地说，将用于补偿估计“Ne”值和测量“Ne”值之间偏差所需的估计发动机转速“Ne”值赋给上述方程式(4e)，从而能够计算为获得所需的估计“Ne”值而必需的输入转矩TRQ。之后，能够使用计算出的输入转矩TRQ及上述方程式(5)和(6)来计算为获得所需的估计“Ne”值而需要的燃烧压强P。如上所述，修改方程式(9a)中的参数以便能够获得燃烧压强P。使用此技术，使得通过持续地反馈测量“Ne”值和估计“Ne”值之间的误差，能够获得燃烧波形计算块70的高计算精度。因此，燃烧压强P的估计值和测量值能够精确地匹配。

尽管在上述步骤102的处理中，如果在测量发动机转速“Ne”值和估计发动机转速“Ne”值之间存在误差，则方程式(9a)中的参数被判断为不正确，并因此对参数进行修改。用于从发动机转速“Ne”数据误差来获知燃烧压强P的技术不限于上述技术，并且可以是例如下文所述的技术。也就是说，不是对方程式(9a)中的参数进行修改，而是可以基于方程式(9a)和通过在方程式(10)的右侧增加修正项而建立的方程式来计算燃烧压强P，所述修正项是曲轴转角θ的函数(即，用于合并发动机转速“Ne”数据误差的修正项)。

(涉及步骤104的处理)

在图6的程序中，如果在步骤100中判断为经济运行模式标记开启，则在此判断之后为步骤104，在步骤104中，判断实际发动机转速“Ne”是否高于当前设定的用于控制曲柄的停止位置的燃烧切断发动机转速“Ne0”。

(涉及步骤106的处理)

结果，如果没有判断为燃烧切断发动机转速“Ne0”小于实际发动机转速“Ne”，即，如果判断为实际发动机转速“Ne”已经达到燃烧切断发动机转速“Ne0”，则在步骤106中计算曲柄停止位置的估计值。

更确切地说，在怠速状态期间在步骤102的处理中获得的平均燃烧压强值P、进气压强“Pmax”、曲轴转角“θ0”和发动机转速(燃烧切断发动机转速)“Ne0”(＝曲柄角转速dθ/dt)是作为步骤106中的初始值的被首先输入的。然后，使用曲轴外围运动方程计算块62依次计算曲轴转角θ和曲柄角转速dθ/dt的各自的估计值。下面，将使用方程式(13)和(14)描述此计算技术的细节。使用此技术在上述图2所示的各个箭头方向上求解发动机模型60，这在此说明书中被称为“正向模型计算”。

首先，将表达式 $(\∂f\left(\mathrm{\θ}\right)/\∂\mathrm{\θ}=h\left(\mathrm{\θ}\right))$ 代入由方程式(4e)表达的曲柄外围运动方程式中。然后，在将方程式(5)赋值给方程式(4e)中的输入转矩TRQ之后，使方程式(4e)离散化以获得下面所示的方程式(13)。

$\{\mathrm{\θ}(k+2)-\mathrm{\θ}(k+1)\}-\{\mathrm{\θ}(k+1)-\mathrm{\θ}\left(k\right)\}$

$=[{\mathrm{TRQ}}_e\left(\mathrm{\θ}\left(k\right)\right)-{\mathrm{TRQ}}_{\mathrm{fr}}(\mathrm{\θ}(k+1)-\mathrm{\θ}\left(k\right))$

$-\frac{1}{2}h\left(\mathrm{\θ}\left(k\right)\right)\·{(\mathrm{\θ}(k+1)-\mathrm{\θ}\left(k\right))}^2]/f\left(\mathrm{\θ}\left(k\right)\right)$ ........(13)

然后，曲轴转角“θ0”、曲柄角转速dθ/dt和其它值作为基于方程式(13)的正向模型计算期间的初始计算值被给出。随后，继续更新步进计数(step count)“k”并因此顺序计算曲轴转角θ和曲柄角转速dθ/dt的各自的估计值。通过将步进计数“k＝1”赋值给方程式(13)，方程式(13)可以表达为方程式(14a)。

如果k＝1，那么：

$\{\mathrm{\θ}\left(2\right)-\mathrm{\θ}\left(1\right)\}-\{\mathrm{\θ}\left(1\right)-\mathrm{\θ}\left(0\right)\}$

$=[{\mathrm{TRQ}}_e\left(\mathrm{\θ}\left(0\right)\right)-{\mathrm{TRQ}}_{\mathrm{fr}}(\mathrm{\θ}\left(1\right)-\mathrm{\θ}\left(0\right))$

$-\frac{1}{2}h\left(\mathrm{\θ}\left(0\right)\right)\·{(\mathrm{\θ}\left(1\right)-\mathrm{\θ}\left(0\right))}^2]/f\left(\mathrm{\θ}\left(0\right)\right)$ ........(14a)

其中，由于 $\stackrel{\·}{\mathrm{\θ}}\left(1\right)=\{\mathrm{\θ}\left(2\right)-\mathrm{\θ}\left(1\right)\},$ $\stackrel{\·}{\mathrm{\θ}}\left(0\right)=\{\mathrm{\θ}\left(1\right)-\mathrm{\θ}\left(0\right)\},$

接下来：

$\stackrel{\·}{\mathrm{\θ}}\left(1\right)-\stackrel{\·}{\mathrm{\θ}}\left(0\right)$

$=[{\mathrm{TRQ}}_e\left(\mathrm{\θ}\left(0\right)\right)-{\mathrm{TRQ}}_{\mathrm{fr}}\left(\stackrel{\·}{\mathrm{\θ}}\left(0\right)\right)-\frac{1}{2}h\left(\mathrm{\θ}\left(0\right)\right)\·{\left(\stackrel{\·}{\mathrm{\θ}}\left(0\right)\right)}^2]/f\left(\mathrm{\θ}\left(0\right)\right)$

........(14b)

$\⇔\stackrel{\·}{\mathrm{\θ}}\left(1\right)=[{\mathrm{TRQ}}_e\left(\mathrm{\θ}\left(0\right)\right)-{\mathrm{TRQ}}_{\mathrm{fr}}\left(\stackrel{\·}{\mathrm{\θ}}\left(0\right)\right)-\frac{1}{2}h\left(\mathrm{\θ}\left(0\right)\right)\·{\left(\stackrel{\·}{\mathrm{\θ}}\left(0\right)\right)}^2]$

$/f\left(\mathrm{\θ}\left(0\right)\right)+\mathrm{\θ}\left(0\right)$

$\⇔\stackrel{\·}{\mathrm{\θ}}\left(1\right)=[{\mathrm{TRQ}}_e\left(\mathrm{\θ}\left(0\right)\right)-{\mathrm{TRQ}}_{\mathrm{fr}}\left(\stackrel{\·}{\mathrm{\θ}}\left(0\right)\right)-\frac{1}{2}h\left(\mathrm{\θ}\left(0\right)\right)\·{\left(\stackrel{\·}{\mathrm{\θ}}\left(0\right)\right)}^2]$

$/f\left(\mathrm{\θ}\left(0\right)\right)+{\stackrel{\·}{\mathrm{\θ}}}_0$ ........(14c)

$\mathrm{\θ}\left(1\right)=\mathrm{\θ}\left(0\right)+\stackrel{\·}{\mathrm{\θ}}\left(0\right)=\mathrm{\θ}\left(0\right)+{\stackrel{\·}{\mathrm{\θ}}}_0$ ........(14d)

用相关的曲柄角转速dθ/dt替换上述方程式(14a)中的曲轴转角“θ(k)”的部分，生成方程式(14b)。接下来，通过使用曲轴转角“θ0”和最终作为初始值被输入的曲柄角转速dθ0/dt展开方程式(14b)使其能够表达步进计数“k＝1”的曲柄角转速dθ(1)/dt，如方程式(14c)。另外，对上述方程式(14c)积分使其能够表达步进计数“k＝1”的曲轴转角“θ(1)”，如方程式(14d)。

接下来，当上述处理重复进行直至步进计数“k”达到N时，即，曲柄角转速已经变为dθ(N)/dt＝0，则计算曲柄角转速dθ(N)/dt＝0和曲轴转角“θ(N)”。这意味着，可以使用以上处理来计算在内燃机10停止期间的发动机转速“Ne”(＝0)和曲柄停止位置。

在步骤106的处理中，由于以下原因，将燃烧压强P的平均值用作初始值。即，在从判断为实际发动机转速“Ne”等于燃烧切断发动机转速时至燃烧实际上被切断时的期间内，因为与内燃机10的硬件性能相关的原因，需要等同于大约90°CA的相差的时间。在此时间内，发动机转速“Ne”的下降程度基于燃烧切断开始前存在的燃烧压强P的值而变化。更确切地说，如果在燃烧切断开始前存在的燃烧压强P高，则发动机转速“Ne”的下降慢。相反，如果燃烧压强P低，则发动机转速“Ne”下降快。根据步骤106的处理，通过将燃烧压强P的平均值作为初始值，燃烧切断开始前存在的燃烧压强P的影响可以反映在曲柄停止位置的估计值上。

在步骤106中，当燃烧压强P的最大值和最小值为燃烧切断开始前存在的燃烧压强P时，同样计算曲柄停止的估计值。最小和最大燃烧压强在步骤102的处理中已经保存。更确切地说，使用与保存的最小燃烧压强和最大燃烧压强对应的输入转矩TRQ值作为初始值，来计算曲柄停止位置。根据此处理，由于同样计算了在最小和最大燃烧压强处的停止位置，所以能够获知由于内燃机10的燃烧状态的改变而引起的曲柄停止位置的变化。

步骤106中的上述处理可以由以下述的处理替代。在步骤106中，为了考虑燃烧状态的变化，将在怠速期间预先获取的燃烧压强P的平均值作为初始值，并且将与最小和最大燃烧压强相关的输入转矩TRQ值作为初始值。然而，除了使用这些值以外，在所需的循环期间内不断地计算燃烧压强P的平均值。当判断为实际发动机转速“Ne”已经达到燃烧切断发动机转速，则将在燃烧切断即将开始之前存在的上述平均值作为初始值。

(涉及步骤108的处理)

接下来，在图6的程序中，在步骤108中判断在步骤106的处理中计算出的估计曲柄停止位置值和由曲轴转角传感器40检测到的测量曲柄停止位置值是否大于所需的极限水平。如果判断为所述误差大于所需的极限水平，则在步骤110中执行摩擦模型64的适应性获知然后对其更新。

在步骤108中，对用于判断摩擦模型64的适应性获知是否应该开始的初始判断值同样进行修改，从而使得上述极限水平能够根据内燃机10的特定运行状态或其特定的劣化状态(deterioration state)而变化。为此，ECU 50保存表示曲轴16的这种目标停止区域的数据，该数据被限定为如图7所示。为了在执行经济运行控制期间适当地重新启动内燃机10，曲轴16(活塞12)需要被控制在适当的停止位置上。因此，如图7所示，在压缩冲程期间在达到上止点之前的90°CA-180°CA的范围内的区域由ECU 50设定为曲轴16的目标停止区域。

如果判断为在步骤106的处理中计算出的处于最小燃烧压强和最大燃烧压强时的两个曲柄停止位置中的至少一个偏离上述目标停止区域的两侧的必要范围(例如，如图7所示的目标停止区域的+10％的范围)，则上述极限水平改变。更确切地说，最小燃烧压强时的曲柄停止位置和最大燃烧压强时的曲柄停止位置被保存，然后根据曲柄的所述停止位置之间的特定偏差来改变上述极限水平。更确切地说，如果偏差明显，则表明发动机模型60的计算精度下降。因此，通过增加摩擦模型64的获知频率来降低极限水平，以提高发动机模型60的计算精度。通过考虑内燃机10的燃烧水平的变化来获知摩擦转矩的影响，上述技术的使用能够提高曲柄的停止位置的估计精度。

另外，用于判断在所述程序中摩擦模型64是否需要获知的技术并不局限于步骤108的处理方法。例如，如果将使用上述方法计算出的曲柄停止位置的估计值作为对步骤106的处理的改进，即，使用采用在燃烧切断即将开始前获得的平均燃烧压强值P作为初始值的方法，当判断为曲柄的估计停止位置不处于必要范围(例如，排除目标停止位置的10％的部分的范围，如图7所示)内时，可以在步骤110中进行摩擦模型64的适应性获知然后对其进行更新。

此外，在所述程序中用于使极限水平根据内燃机10的特定运行状态或劣化状态而改变的技术不限于步骤108的处理方法。例如，在根据偏差来判断此劣化水平之后，可以根据内燃机10内的燃烧的特定劣化水平来改变极限水平，所述偏差是在使用采用恰好在燃烧切断前获得的平均燃烧压强值P作为初始值的方法估计出的曲柄停止位置和在最小或最大燃烧压强时的曲柄停止位置之间的偏差。更确切地说，可以将在最小燃烧压强时的曲柄停止位置和在最大燃烧压强时的曲柄停止位置之间存在明显偏差的状态判断为内燃机10的燃烧水平发生显著变化的状态。因此，如果偏差明显，可以判断为内燃机10的燃烧已经劣化。

因而，如果对应的偏差明显，即，如果内燃机10的燃烧的劣化水平明显，则可以增加用于开始摩擦模型64的获知的所述极限水平。使用所述技术可以通过降低摩擦模型64的获知频率来防止错误的获知。另外，尽管使用所述技术需要单独地获知燃烧压强P的任何变化(即，燃烧水平变化的影响)，但只需要在此步骤的阶段中更新摩擦模型64即可。例如，可以在本文后面描述的步骤114中进行燃烧压强P的变化的获知。下面，将描述所述获知方法。

在获知燃烧压强P的变化时，燃烧压强P的改变与曲柄停止位置的偏差之间的关系是基于在最小和最大燃烧压强时的曲柄停止位置上计算出的数据而保存的。这使得可以根据燃烧压强P的变化判断出曲柄停止位置的偏差很可能发生。基于曲柄停止位置的所述偏差，在下述的步骤114中修改在计算燃烧切断发动机转速期间使用的目标停止位置。根据这种获知，能够考虑燃烧水平的变化来精确地计算出在获得期望的目标停止位置时的燃烧切断发动机转速。

(涉及步骤110的处理)

如上所述，在图6的程序中，如果在步骤108中判断为以上误差超出了所需的极限水平，则在步骤110中执行摩擦模型64的适应性获知然后对其进行更新。

曲柄停止位置的计算需要由摩擦模型64计算出的摩擦转矩TRQf。所述计算出的摩擦转矩TRQf值依赖于发动机转速“Ne”。当燃烧停止并且内燃机10停止时，发动机转速“Ne”立即改变并且经过在摩擦模型64的摩擦设定表上的不同的映射点。因而，当发动机转速“Ne”发生每一个变化时，很难立即获知摩擦模型64。

因此，在步骤110中，使用下述方法获知并更新摩擦模型64。顺便要说的是，当获知摩擦模型64时，在由发动机模型60计算出的曲轴转角θ和测量出的曲轴转角θ之间的偏差(误差)，以及在由发动机模型60计算出的曲柄角转速dθ/dt和测量出的曲柄角转速dθ/dt之间的偏差(误差)，也都可以通过PID控制器76的定期反馈控制来补偿。在步骤110中，所述补偿是作为用于获知摩擦模型64的预备操作而进行的。更确切地说，通过使所需的反馈增益与所述偏差相乘而获得的修正数据(下面，所述修正数据被称为“PID修正数据”)是由PID控制器76计算出来的，然后所述PID修正数据被并入曲轴外围运动方程计算块62的输入数据中。进行所述处理是为了防止摩擦转矩误差影响发动机转速“Ne”。例如，如果摩擦转矩TRQr的偏差在发动机转速“Ne”为600rpm时出现，则在发动机转速“Ne”降低的同时，会持续保存所述摩擦转矩误差。因此，摩擦转矩误差会影响处于发动机转速“Ne”低于600rpm的区域的发动机。然而，上述处理防止了这种影响。

由PID控制器76进行的上述处理能够提高发动机模型60的计算精度，以使在曲轴转角θ和曲柄角转速dθ/dt中出现的估计误差最小化。然而，摩擦模型64的精度不会提高。因此，在步骤110中，使用PID修正数据进行获知摩擦模型64的以下处理：

首先，每当摩擦转矩经过已经被作为摩擦模型64的摩擦设定表上的映射点而赋值的发动机转速“Ne”的不同的设定点(setting)时，更确切地说，例如，如果燃烧在600rpm时被切断，则当发动机转速“Ne”下降至诸如500rpm和400rpm的转速范围时，在适当的时刻计算并保存PID修正数据。然而，在此时，并没有获得与各个映射点完全对应的PID修正数据。作为替代，为了获得与类似噪声的行为无关的修正数据，保存在邻近的或不同的映射点上经过必要次数而计算出来的平均PID修正值，或者保存在各个映射点的附近的时间积分PID修正数据，作为在各映射点附近的PID修正数据。所有与模型计算开始(燃烧切断的开始)发生的偏差都被保存在PID修正数据的第I项(积分项)分量中。因此，为了在各映射点的附近计算PID修正数据，则从此次计算出的PID修正数据中减去直至前次已经被保存的所有偏差。

图8为说明在步骤110中执行的用于将PID修正数据并入到摩擦模型64的摩擦设定表上的各个映射点的值中的方法的图解。图8示出了摩擦转矩TRQf和数据采样点数N之间的关系。这里，在冷却水温度和发动机转速“Ne”之间的关系中定义的摩擦设定表上的各个映射值被首先绘制在图8所示的图解上。接下来，使用最小二乘方的方法计算相对于绘制点的回归线(回归曲线)。所述回归线是在冷却水温度“Thw”和发动机转速“Ne”之间的关系下定义的，此线可以表达为方程式(15)。

TRQ_f＝aNe+bT_hw+c             ......(15)

其中“a”、“b”、“c”是被确定为使各个绘制点至回归线的距离最小的系数。

接下来，ECU 50使用摩擦模型64为保存有PID修正数据的各个发动机转速“Ne”区域来计算当前冷却水温度“Thw”下的摩擦转矩TRQf。之后，由摩擦模型64已经计算出的摩擦转矩TRQf值被绘制在图8的曲线图上。紧接着在添加完这些新的绘制点之后，再次确定所述方程式(15)中的系数“a”、“b”和“c”，然后由此获得的新回归线被用于计算对于各个映射点的加权系数。

加权系数用来乘以上述PID修正数据，如图8所示，各个系数根据各个映射点与新回归线之间的特定距离而改变。更确切地说，当距离变短时，各加权系数变为趋近于1的值。相反，当距离变长时，加权系数变为趋近于0的值。接下来，使用通过将PID修正数据乘以加权系数而获得的数据来修正各映射点的值。选定这些加权系数来进行修正使得合并的PID修正数据在更加远离回归线的映射点处减小。

这样，为各个发动机转速区域而获得的上述PID修正数据被并入到各映射点的值中。即使当前水温值是邻近映射的中间值，也能够使用上述的方法将适当的PID修正数据赋予考虑了其与邻近映射的关系各映射点。在步骤110中，通过上述处理来获知并更新摩擦模型64。

(涉及步骤112的处理)

在图6的程序中，在步骤112中随后获知曲柄停止位置的变化。依照独立于此程序的程序处理流程，ECU 50持续地保存在不同燃烧切断发动机转速下的曲柄停止位置的测量变化。在步骤112中，如果判断为在当前燃烧切断发动机转速下的曲柄停止位置的测量变化已经增加，这表明从摩擦模型64获得的摩擦转矩信息出现增加的误差，因而，用于判断摩擦模型64的适应性获知是否应该开始的值被改变。更确切地说，如果判断为对应的变化已经增加，则在步骤108中的极限水平被更新为更小的值，以提高发动机模型60的计算精度。

(涉及步骤114的处理)

接下来，使用在步骤110的处理中被更新的摩擦模型64，在步骤114中对用于将曲柄的停止位置修改为其期望的目标停止位置的目标燃烧切断发动机转速进行计算并更新。通过执行步骤114的此处理，能够获取目标燃烧切断发动机转速，所述目标燃烧切断发动机转速合并了在上述一系列处理步骤中已经获知的摩擦的影响。

例如，三种方法可用于实现步骤114的处理。第一种可能的方法是逆向模型计算(inverse model computation)，其用于在与正向模型计算(forward model computation)完全相反的方向上求解上述发动机模型60。所述逆向模型计算涉及到内燃机10即将停止前发生的行为。在这种状态中，内燃机10内部的影响实质上是没有变化的并且发动机转速也没有明显的起伏。因而，通过利用曲轴16的目标曲柄停止位置(曲轴转角)作为初始值来对发动机模型60进行逆向模型计算，能够计算出目标燃烧切断发动机转速(等同于在正向模型计算中的上述初始曲柄角转速dθ/dt)。

然而，由于发动机模型60在正向模型计算中通常是作为二阶滞后系统运算的，所以与发动机模型60相反的模型是作为二阶微分系统运算的。因此，由于逆向模型计算系统是发散系统，不容易保证稳定的计算。因而，第二种可能的方法是使用参照图9描述的反馈。在图9的方法中，计算由上述发动机模型(停止位置估计模型)60计算出的曲柄停止位置和相关的曲柄的目标停止位置之间的偏差，并且所述偏差由PID控制器76并入到燃烧切断发动机转速(初始曲柄角转速(dθ/dt)₀)中，所述燃烧切断发动机转速是发动机模型60的输入值。使用图9所示的方法多次重复所述处理能够获取没有上述偏差的燃烧切断发动机转速，即，与曲柄目标停止位置对应的目标燃烧切断发动机转速。当发动机模型60在步骤110的处理期间被获知后，此时的发动机模型60使用摩擦模型64。因此，通过合并获知的摩擦转矩的影响，使用所述方法使得与曲柄目标停止位置对应的目标燃烧切断发动机转速的计算精度非常高。

上述第二种方法的优点是，总能够获得稳定的计算结果并且因此能够执行可靠的计算。然而，同时也具有缺点：ECU 50的计算负荷趋于增加并且计算的收敛性下降。这使得很难执行频繁的计算。因此，下面将描述一种作为第三种可能方法的使用发动机模型60的正向模型计算结果的方法，所述发动机模型60使用在步骤110的处理中被获知并更新的摩擦模型64。

图10为用于说明在图6所示的程序的步骤114中进行的处理的时间图。更确切地说，图10A示出了根据燃烧切断发动机转速的特定值而减小的发动机转速“Ne”；图10B示出了依照燃烧切断发动机转速的值而改变的曲柄停止位置；而图10C示出了在燃烧切断发动机转速和与其相关的曲柄停止位置之间的关系。

确切地说，在步骤114中，使用燃烧切断发动机转速“Ne0”(初始曲柄角转速(dθ/dt)₀)和在此程序开始期间在步骤106中输入的初始曲轴转角θ0作为初始值，首先利用在步骤110的处理中被更新的摩擦模型64执行正向模型计算。因此，计算在对应(dθ/dt)₀的燃烧切断发动机转速下的曲柄停止位置“θx0”。除了“θx0”之外，通过进行类似的正向模型计算，在步骤114中分别计算在燃烧切断发动机转速(初始曲柄角转速)(dθ/dt)₀+10和(dθ/dt)₀-10下的曲柄停止位置“θx1”和“θx2”。当燃烧切断发动机转速对应于(dθ/dt)₀+10和(dθ/dt)₀-10时同时输入的初始曲轴转角为θ0，这对于dθ0/dt是一样的。另外，可以根据燃烧的特定劣化水平和/或摩擦增加量(依赖于诸如润滑油劣化的因素)改变前面描述中的±10rpm的范围。

如图10A所示，如果将不同的值赋予燃烧切断发动机转速，这会引起使发动机转速“Ne”下降至零所需的时间的不同。另外，如图10B所示，如果为了获得(dθ/dt)₀±10rpm而赋予不同的燃烧切断发动机转速值，则最终的曲柄停止位置会在θx1至θx2的范围内变化。更确切地说，燃烧切断发动机转速的增加使曲柄停止位置向更靠前的位置转移。因此，进行上述处理能够获得一个范围，在该范围内，当为了获得(dθ/dt)₀±10rpm而改变燃烧切断发动机转速时，曲柄停止位置会改变。

接下来在步骤114中，从以上述方式已经获得的曲柄停止位置变化范围相对于燃烧切断发动机转速变化范围的关系来计算目标燃烧切断发动机转速，使得曲柄停止位置与其期望的停止位置匹配。更确切地说，如图10C所示，逼近线(approximated line)是从燃烧切断发动机转速和由上述三种正向模型计算获得的曲柄位置停止数据之间的关系计算出来的。获得此逼近线能够计算与曲轴16的期望目标曲柄停止位置相对应的目标燃烧切断发动机转速。

同样地，下述的方法用于降低在三种正向模型计算期间的ECU 50的计算负荷。更确切地说，与其他项的大小相比较，上述方程式(13)中的“f(θ(k))”可以用作实质上固定的值。另外，尽管受初始曲轴转角θ(0)的影响，但是“h(θ(k))”不受初始曲柄角转速((dθ/dt)₀)的影响。此外，“TRQe(θ(k))”并不依赖于那些初始值而是依赖于直至计算完成一直存在的步进计数“k”。

由于方程式(13)中的各个上述参数具有上述特性，所以对上述三个初始值中的每一个只计算方程式(13)右侧的“TRQf(θ)”。严格来说，“TRQe(θ(k))”、“h(θ(k))”和“f(θ(k))”的每一个都依赖于曲轴转角θ。然而，在方程式(13)中，如果用作平均值的话，“TRQe(θ(k))/f(θ(k))”项和“1/2·h(θ(k))·(θ(k+1)-θ(k))²/f(θ(k))”项几乎变为固定值。因此，各个参数可以计算为与步进计数相乘的固定值(例如，如果燃烧切断发动机转速为(dθ/dt)₀，则那些参数中的每一个都可以计算为与步进计数N相乘的固定值)。

另外，在三种正向模型计算期间，当在用于(dθ/dt)₀的燃烧切断发动机转速的曲柄角转速dθ/dt变为零时的步进计数为N时，如果燃烧切断发动机转速变为((dθ/dt)₀+10)，则在步进计数N时的曲柄角转速dθ/dt不变为零。当步进计数Np加上步进计数N时，即，当达到(N+Np)的步进计数时，曲柄角转速变为零。相反，如果燃烧切断发动机转速变为((dθ/dt)₀-10)时，则曲柄角转速dθ/dt在(N-Nm)的步进计数时变为零。因此，保存步进计数Np、Nm。因而，如上所述，当“TRQe(θ(k))/f(θ(k))”项和“1/2·h(θ(k))·(θ(k+1)-θ(k))²/f(θ(k))”项各自被计算为与步进计数相乘的固定值时，所述各项可以被计算为与步进计数(N+Np)相乘的用于((dθ/dt)₀+10)的燃烧切断发动机转速的固定值，或者计算为与步进计数(N-Nm)相乘的用于((dθ/dt)₀-10)的燃烧切断发动机转速的固定值。

根据上述方法，可以基于发动机模型60的正向模型计算结果来计算与曲轴16的目标曲柄停止位置相对应的目标燃烧切断发动机转速。

由发动机模型60计算出的曲柄停止位置与曲柄的实际停止位置之间的偏差很可能是由摩擦模型64的影响和发动机模型60中其它部分的影响而引起的。很难仅通过获取偏差来判断曲柄停止位置的上述偏差是否仅仅是由摩擦模型64的影响引起的。换句话说，很难判断所述偏差是由计算摩擦转矩TRQf引起的还是由计算发动机转矩TRQe引起的。

与此相反，根据图6的上述程序，燃烧压强P是在内燃机10内的燃烧期间获知的。之后，摩擦是在完成燃烧切断后在内燃机10的运行即将停止之前获知的。根据上述程序的处理，由于因燃烧压强的改变引起的曲柄停止位置的变化是在执行燃烧期间获知的，因此除了其它与燃烧相关的影响外，只有对曲柄的停止位置产生显著影响的摩擦影响能够在内燃机10的即将停止前被独立地获知。因而，可以更新摩擦模型64。从以上可以理解，根据本实施例中的发动机模型60，由于摩擦、大气压强、大气温度、节气门角度、气门正时VVT和其它影响曲柄停止位置的因素(在本发明中提到的“必要参数”)的影响被适当地模拟，所以对摩擦的适应性获知的精度可以被提高。

另外，根据上述程序的处理，使用诸如利用并入的已被获知并更新的摩擦模型64的逆向模型计算的方法，可直接计算用于使曲柄停止位置与其目标停止位置匹配的燃烧切断发动机转速。因此，可以实现能够执行高效算术运算的发动机模型60。

在如上所述的第一实施例中，根据本发明第二个方案的“曲柄停止位置计算部件”是在ECU 50执行步骤106时实现的。

此外，根据本发明第三个方案的“曲柄停止位置判断部件”是在ECU 50执行步骤108时实现的。

此外，根据本发明第四个方案的“目标燃烧切断发动机转速计算部件”是在ECU 50执行步骤114时实现的。

第二实施例

现在，将参照图11来描述本发明的第二实施例。

通过采用图1所示的硬件配置和图2所示的发动机模型60，并通过允许ECU 50执行图11所示的程序而不是图6所示的程序，实现了第二实施例的系统。

为了提高在执行经济运行控制期间的发动机再起动性(restartability)，通过使用所述发动机模型60，适当地控制了在内燃机10中的曲柄停止位置。由于摩擦是一个非常有影响的因素，为了适当地控制曲柄的停止位置，有必要以高精度来估计内燃机10内部的摩擦(摩擦转矩TRQf)。摩擦量通常依赖于内燃机10的润滑油温度。然而，典型地，摩擦是基于冷却水温度来修正的，如同上述的摩擦模型64。当润滑油温度和冷却水温度之间存在明显差值时，由于很难精确地修正摩擦，因此对曲柄停止位置的估计精度就降低了。

因而，参照图11所示的程序，本实施例使用下述方法。同样地，即使没有用于检测润滑油温度的传感器，本实施例也能够通过使用所述方法来判断润滑油温度和冷却水温度之间的差值。因此，能够使用廉价的系统配置来实现对经济运行控制期间的曲柄停止位置的高估计精度。

【第二实施例中的详细处理】

图11是为了实现上述功能性，ECU 50在本发明的第二实施例中执行的程序的流程图。在图11的程序中，首先在步骤200中判断内燃机10当前是否处于稳定运行状态中。进行所述判断步骤是为了在步骤202的下列处理中获得稳定的计算数据：

如果在步骤200中判断为内燃机10当前处于稳定运行状态，则依照前面的方程式(4e)在发动机模型60中计算发动机转矩TRQe。同样计算用于必要循环数量(例如，几百个循环)的与发动机转矩TRQe相对应的平均值A。所述计算是在步骤202中进行的。

接下来，在步骤204中获取内燃机10的基本转矩(base torque)B。更确切地说，获取到了由内燃机10的进气量、发动机转速“Ne”、空/燃(A/F)比(燃料喷射量)、点火正时和指示转矩(indicated torque)之间的关系定义的设定表。通过从指示转矩中减去摩擦转矩和负载转矩来计算基本转矩B。指示转矩、摩擦转矩和负载转矩是在冷却水温度和润滑油温度彼此相关的状态中预先测量出的数据。

接下来，在步骤206中判断在由发动机模型60计算出的发动机转矩TRQe的平均值A和基本转矩B之间的转矩偏差ΔTRQ的绝对值是否大于所需判断值。因此，如果判断为转矩偏差ΔTRQ的绝对值小于所需判断值，即，如果判断为当前发动机转矩TRQe是一个接近在冷却水温度和润滑油温度相关状态下获得的基本转矩B的值，则这表明润滑油温度和冷却水温度之间的差值不太明显。在这种情况下，允许在步骤208中执行经济运行模式。

如果在步骤206中判断为转矩偏差ΔTRQ的绝对值大于所需判断值，则与上述情况相反，这表明在发动机模型60中获取到的摩擦转矩TRQf存在明显变化，即，润滑油温度与冷却水温度之间的差值相对要明显。在这种情况下，禁止在步骤210中执行经济运行模式。

为了便于描述，以下描述涉及到基于流程图的一个程序中的处理。然而，如果在步骤206中判断为转矩偏差ΔTRQ的绝对值大于所需判断值，则可以进行下列处理。也就是说，如果判断出润滑油温度和冷却水温度之间的差值相对明显，则通过求解方程式(16)在步骤212中计算目标燃烧切断发动机转速，在方程式(16)中，上述转矩偏差(摩擦误差)ΔTRQ被加在方程式(4e)的右侧。在这种情况下，允许在步骤208中执行经济运行模式。这样，可以在获知摩擦误差之后继续执行经济运行模式。

根据图11的上述程序，如果在冷却水温度和润滑油温度之间存在相对明显的差值，则禁止执行经济运行模式。也就是说，如果摩擦转矩TRQf偏差明显并且无法精确估计，则禁止执行经济运行模式。因此，使用发动机模型60可以精确地估计出曲柄的停止位置。

另外，根据在以上程序的描述中阐明的其它处理，如果润滑油温度和冷却水温度之间的差值相对明显，如方程式(16)所示，通过求解方程式来计算目标燃烧切断发动机转速，在所述方程式中，上述转矩偏差ΔTRQ加在从摩擦模型64获得的摩擦转矩TRQf上。通过使用方程式(16)计算目标燃烧切断发动机转速对应于，基于作为摩擦转矩TRQf的误差的ΔTRQ值来修正摩擦模型64然后通过使用修正过的摩擦模型64来计算目标燃烧切断发动机转速，在方程式(16)中，转矩偏差ΔTRQ被加在从摩擦模型64获得的摩擦转矩TRQf上。因此，基于上述程序的处理能够反映如下事实：摩擦受到冷却水温度和润滑油温度的影响。因此，即使当冷却水温度和润滑油温度之间存在相对明显的差值时，也能够精确地估计曲柄的停止位置。

在如上所述的第二实施例中，根据本发明第四个方案的“温差判断部件”是在ECU 50执行步骤206时实现的；而根据本发明第五个方案的“摩擦误差计算部件”是在ECU 50执行步骤212时实现的。

第三实施例

现在，将参照图12来描述本发明的第三实施例。

通过采用图1所示的硬件配置和图2所示的发动机模型60，并通过允许ECU 50执行图12所示的程序而不是图11所示的程序，实现了根据第三实施例的系统。

本实施例的系统与上述第二实施例的系统不同之处仅在于，判断冷却水温度和润滑油温度之间差值的方法。下面，将依照图12所示的程序，主要对所述方法进行描述。

【第三实施例中的详细处理】

图12为ECU 50在本发明的第三实施例中执行的程序的流程图。在图12的程序中，如果在步骤300中判断为内燃机10当前处于稳定运行中，则在步骤302中计算基于运行参数的记录的评价函数“f”。由以下方程式(17)表达的评价函数“f”是一个考虑到由燃烧引起的典型热释放量的计算的函数。

f＝∫KL·K_Ne·K_abyf·K_aop         ......(17)

如以上方程式(17)所表达的，评价函数“f”被赋值为通过在所需周期上进行时间积分而获得的值，通过使修正系数“Kne”与充气率(air-filling rate)KL相乘而获得的值取决于发动机转速“Ne”，修正系数“Kabyf”取决于空/燃(A/F)比，而修正系数“Kaop”取决于点火正时。恰好在内燃机持续以高功率输出长时间运行后，诸如恰好在严苛的条件下驾驶之后，冷却水温度和润滑油温度之间的差值增大。

充气率KL的增大会增加热释放量。另外，如方程式(17)所示，如果使对应发动机转速“Ne”的修正系数“Kne”与充气率KL相乘，则可以获得关于燃烧气体的速率的信息。此外，使取决于空/燃(A/F)比的修正系数“Kabyf”与KL×Kne相乘提供了适当的燃烧温度信息，这取决于空/燃(A/F)比是浓(rich)空/燃比还是稀(lean)空/燃比。此外，使取决于点火正时修正系数“Kaop”与KL×Kne×Kabyf相乘提供了是否能够从表明点火正时是否被延迟的信息中获得发动机已经被撞击的信息。因此，使用考虑了上述影响的评价函数“f”能够判断驾驶所述车辆是多么严苛。定义各个修正系数K的值，以便当车辆驾驶条件变得更加严苛时，评价函数“f”取更大的值。

接下来，在步骤304中判断评价函数“f”的值是否大于所需判断值。因此，如果判断为评价函数“f”小于判断值，即，如果判断为车辆当前没有处于诸如恰好在严苛条件下驾驶后的状态中，并且因此冷却水温度和润滑油温度之间的差值相对小时，则允许在步骤306中执行经济运行模式。

然而，如果在步骤304中判断为评价函数“f”大于判断值，即，如果判断为车辆当前处于诸如恰好在严苛条件下驾驶后的状态中，并且因此冷却水温度和润滑油温度之间的差值相对明显，则禁止在步骤308中执行经济运行模式。

根据图12的上述程序，即使没有用于检测润滑油温度的传感器，也能够判断冷却水温度和润滑油温度之间的差值。

已经基于以下假定对上述第三实施例进行了描述：如果判断为评价函数“f”大于判断值，则在步骤308中禁止执行经济运行模式。然而，与禁止经济运行模式不同，如果进行了上述判断，可以依照转矩偏差ΔTRQ来修正摩擦模型64，然后可以计算目标燃烧切断发动机转速，以便允许执行经济运算模式。所述处理实质上与图10的程序中的步骤112的处理相同。

Claims (5)

1.一种用于内燃机的停止位置控制装置，其适用于通过控制内燃机的燃烧即将停止时的燃烧切断内燃机转速来控制曲轴的停止位置，所述用于内燃机的停止位置控制装置包括：

摩擦模型，其计算所述内燃机的滑动构件的摩擦；

曲柄停止位置计算部件，其基于包括所述摩擦的必要参数来计算所述曲轴的所述停止位置的估计值；以及

曲柄停止位置判断部件，其判断所述曲轴的所述停止位置；

其中，如果判断为所述曲轴的所述停止位置处于必要误差范围之外，则所述装置基于所述曲轴的所述停止位置的估计值和所述曲轴的所述停止位置的测量值之间的偏差来获知所述摩擦模型并且更新所述摩擦模型；并且

其中所述曲柄停止位置计算部件基于由更新后的摩擦模型计算出的摩擦来修正所述曲轴的所述停止位置的所述估计值。

2.如权利要求1所述的用于内燃机的停止位置控制装置，进一步包括：

温差判断部件，其判断冷却水温度和润滑油温度之间的差值；及

摩擦误差计算部件，其计算与所述温差相关的摩擦误差；

其中，如果判断为所述冷却水温度和所述润滑油温度彼此不同，则所述装置基于所述摩擦误差来修正所述摩擦模型。

3.如权利要求2所述的用于内燃机的停止位置控制装置，

其中，如果判断为所述冷却水温度和所述润滑油温度彼此不同，则所述装置禁止对所述曲轴的所述停止位置的控制。

4.一种用于内燃机的停止位置控制装置，其适用于通过控制内燃机的燃烧即将停止时的燃烧切断内燃机转速来控制曲轴的停止位置，所述用于内燃机的停止位置控制装置包括：

摩擦模型，其计算所述内燃机的滑动构件的摩擦；

目标燃烧切断内燃机转速计算部件，其基于包括所述摩擦的必要参数来计算所述燃烧切断内燃机转速的目标值；

温差判断部件，其判断冷却水温度和润滑油温度之间的差值；及

摩擦误差计算部件，其计算与所述温差相关的摩擦误差；

其中，如果判断为所述冷却水温度和所述润滑油温度彼此不同，则所述装置基于所述摩擦误差来修正所述摩擦模型；并且

其中所述目标燃烧切断内燃机转速计算部件基于由修正后的摩擦模型计算出的摩擦来修正所述目标值。

5.如权利要求4所述的用于内燃机的停止位置控制装置，

其中，如果判断为所述冷却水温度和所述润滑油温度彼此不同，则所述装置禁止对所述曲轴的所述停止位置的控制。

Images