CN105651524B - 汽车试验台架发动机扭振模拟方法 - Google Patents
汽车试验台架发动机扭振模拟方法 Download PDFInfo
- Publication number
- CN105651524B CN105651524B CN201610006317.5A CN201610006317A CN105651524B CN 105651524 B CN105651524 B CN 105651524B CN 201610006317 A CN201610006317 A CN 201610006317A CN 105651524 B CN105651524 B CN 105651524B
- Authority
- CN
- China
- Prior art keywords
- mrow
- mfrac
- msup
- msub
- torque
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M15/00—Testing of engines
- G01M15/04—Testing internal-combustion engines
- G01M15/12—Testing internal-combustion engines by monitoring vibrations
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Combustion & Propulsion (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Engines (AREA)
Abstract
本发明公开了一种汽车试验台架发动机扭振模拟方法，包括如下步骤：1)根据发动机的固有参数，建立发动机“油门开度—转速—扭矩”模型，2)建立惯量扭矩解算模型，获得惯量扭矩TI；3)建立燃烧扭矩解算模型，获得燃烧扭矩TC；4)对平均扭矩TM，惯量扭矩TI，以及燃烧扭矩TC进行拟合，得到电动机的输入扭矩，即总的发动机扭矩Ten。本发明模拟精确度高，实时性高，计算量小，能够满足发动机周期扭矩脉动(扭振)模拟时的高动态控制响应。
Description
技术领域
本发明涉及一种汽车发动机性能试验模拟方法，尤其涉及一种汽车传动系试验台架发动机扭振模拟方法。
背景技术
近年来汽车的NVH问题越来越多的引起工程师的关注，发动机扭矩的周期性波动会引起传动系的扭转振动，从而产生噪声，降低零部件的使用寿命。如果变速器齿轮设计不当，发动机周期扭矩脉动会引起变速器的敲击声，从而影响汽车的NVH品质。利用实物发动机驱动进行NVH试验具有诸多缺点，1)发动机燃油为易燃物，且其工作过程中易产生着火等危险因素，同时其排放的废气也不利于环保。2)NVH试验主要测量被试变速器的噪声问题，而在试验中实物发动机的较高的噪声会影响变速器NVH的测试。3)在传动系开发的初始阶段，通常情况下与其匹配的发动机也在开发测试中，因此很难获得发动机作为驱动在室内台架上对传动系进行NVH等性能测试。综上所述，利用高动态驱动电机代替发动机进行室内传动系台架试验具有安全环保、试验方便等诸多优点。
如果要在室内台架上复现因发动机周期扭转振动引起的NVH问题，就需要一个高精度的扭振模型(模拟活塞连杆的惯量扭矩及各缸的燃烧扭矩)。发动机扭振模拟系统包括一个低惯量高动态特性的驱动电机和一个满足实时控制要求的扭振模型(模拟活塞连杆的惯量扭矩及各缸的燃烧扭矩)。发动机扭矩周期脉动频率很高，如一个4缸2冲程发动机怠速转速为900r/min，则其点火频率为30Hz，又经测试发现一个完整的周期扭振波形至少需要12个点才能完整地再现出来，所以要对转速为900r/min时的发动机周期扭矩振动进行模拟，其通信频率或控制周期至少要达到30×12＝360Hz，如果发动机转速为6000r/min，则其通信频率或控制周期应达到2400Hz，因此发动机模型的计算速度应满足实时控制要求。精确的发动机模型在很多文献中都有论述，如黑箱模型和基于物理结构的模型等，由于这些模型需要大量的发动机开发数据或详细的发动机设计参数，这些参数在传动系试验中不易获得，而且传统的物理模型在模拟发动机燃烧特性时复杂的数学计算很难满足高动态控制响应的要求，因此这些模型不适用于传动系台架的动态模拟。因此迫切需要提供一种满足高动态实时控制要求的发动机扭振模型及模拟方法。
发明内容
针对现有技术存在的上述不足，本发明的目的在于提供一种汽车传动系试验台架发动机扭振模拟方法，精确度高，实时性高，计算量小，能够满足发动机周期扭矩脉动(扭振)模拟时的高动态控制响应。
为了解决上述技术问题，本发明采用的技术方案是这样的：一种汽车试验台架发动机扭振模拟方法，其特征在于：包括如下步骤：
1)根据发动机的固有参数，建立发动机“油门开度—转速—扭矩”模型，通过该“油门开度—转速—扭矩”模型能根据输入的油门开度或转速信号，获得与之对应的扭矩信号，即获得发动机的平均扭矩TM；
2)建立惯量扭矩解算模型，获得惯量扭矩TI：
发动机往复运动部件所产生的往复惯性力为：
式中：m为活塞和连杆的总质量，表示活塞加速度；
由于往复惯性力通过连杆作用在曲柄销上，在曲轴上产生周期性变化的力矩，从而引起轴系的惯量扭振，因此：
式中：x为活塞行程，L为连杆长度，r为曲柄半径，α为曲柄转角，β为曲柄销中心和曲轴旋转中心所在的直线与连杆轴心线之间的夹角，λ为曲柄半径与连杆长度比；
从而得到活塞加速度(近似公式)：
式中ωc为曲柄旋转角速度；
由此得到往复惯性力为：
PI＝-mrωc 2(cosα+λcos2α)；
由往复惯性力所产生的力矩为：
3)建立燃烧扭矩解算模型，获得燃烧扭矩TC：
燃烧扭矩是由气缸内气体燃烧膨胀对活塞引起的力传递到曲轴和飞轮产生的周期性力矩，其燃烧压力为：
式中：Pcomp(αIVC)进气支管压力，Vd(αIVC)气缸容积，V(α)压缩过程中气缸容积(瞬时容积)，k为膨胀系数；
瞬时气缸工作容积为：
式中：为压缩比，Vc为燃烧室容积，D为气缸直径；
曲柄受到的单缸燃烧扭矩为活塞受到的燃烧压力乘以活塞行程随曲柄转角的变化率，即可得到：
式中：P(α)为燃烧压力，Patm大气压力，Ap气缸活塞面积，活塞运动速率；
4)对平均扭矩TM，惯量扭矩TI，以及燃烧扭矩TC进行拟合，得到电动机的输入扭矩，即总的发动机扭矩Ten:
Ten＝TM+TI+TC。
进一步地，当发动机为多缸发动机时，首先通过步骤1)—步骤4)上述方法得到各单缸的扭矩，然后根据缸数和发火顺序进行扭矩叠加，从而获得发动机总的扭矩信号。
与现有技术相比，本发明具有如下优点：
1、本发明的发动机扭振模拟方法，把发动机扭矩分为平均扭矩、惯性扭矩和燃烧扭矩，通过数学模型分别对平均扭矩、惯性扭矩和燃烧扭矩解算，各数学模型计算量小，从而节省了系统运行资源，解决了发动机周期扭振模拟过程中，对扭振模型的实时性和动态响应需求。
2、本发明提供的发动机扭振模拟方法能够模拟发动机在不同缸数、不同曲轴和飞轮转动惯量下的扭振特性，易于移植到不同的控制器中，较好地满足了台架高动态响应的控制需求。
附图说明
图1为发动机扭矩MAP图。
图2发动机扭矩组成原理图。
图3发动机活塞结构示意图。
图4四缸发动机扭矩叠加原理图。
图5汽车传动系试验台架。
图中：1—上位机，2—驱动器，3—永磁同步电机，4—交流异步电机，5—待测变速器，6—编码器，7—转速扭矩传感器。
具体实施方式
下面将结合附图及实施例对本发明作进一步说明。
实施例：本汽车试验台架发动机扭振模拟方法主要应用于室内汽车传动系试验台架，通过驱动电机模拟发动机的周期脉动扭矩(即扭振)，目标是为变速器等传动系部件的性能测试提供一个逼真的发动机工况条件，以达到离合器、变速器等传动系部件的性能测试的目的。
具体实施时，如：对变速器传动性能进行测试，该试验台架如图5所示，其包括上位机1、驱动器2、永磁同步电机3、交流异步电机4、离合器、待测变速器5、编码器6以及转速扭矩传感器7，其中，永磁同步电机3通过离合器与待测变速器5的输入端相连，以驱动待测变速器5；待测变速器5的输出端接交流异步电机4，同时在待测试变速器输入、输出端分别连接转速扭矩传感器7用来采集测试数据，两个电机的转速闭环回馈信号由安装在两个电机上的编码器6获得。为了满足快速响应的闭环控制要求，在进行测试之前，输入发动机的物理尺寸参数和压缩比等燃烧参数，在试验过程中，本汽车试验台架发动机扭振模拟方法根据输入的发动机参数、台架编码器6反馈的实时转速和油门踏板开度信号计算出总的发动机周期扭矩脉动目标曲线，发送给电机驱动器2，由电机驱动器2驱动电机模拟发动机实际工作过程的周期扭矩脉动，为被试变速器等传动系部件提供一个逼真的发动机周期扭矩脉动的工作环境，便于在室内台架上对被试变速器的NVH等性能进行测试。
参见图1—图3，本汽车试验台架发动机扭振模拟方法，包括如下步骤：
1)根据发动机的固有参数，建立发动机“油门开度—转速—扭矩”模型，发动机“油门开度—转速—扭矩”模型的建立为现有技术，通常根据发动机“油门开度—转速—扭矩”MAP图进行建立。通过该“油门开度—转速—扭矩”模型能根据输入的油门开度或转速信号，获得与之对应的扭矩信号，即获得发动机的平均扭矩TM。
2)建立惯量扭矩解算模型，获得惯量扭矩TI：
发动机往复运动部件所产生的往复惯性力为：
式中：m为活塞和连杆的总质量，表示活塞加速度；
由于往复惯性力通过连杆作用在曲柄销上，在曲轴上产生周期性变化的力矩，从而引起轴系的惯量扭振，因此：
式中：x为活塞行程，L为连杆长度，r为曲柄半径，α为曲柄转角，β为曲柄销中心和曲轴旋转中心所在的直线与连杆轴心线之间的夹角，λ为曲柄半径与连杆长度比；
从而得到活塞加速度(近似公式)：
式中ωc为曲柄旋转角速度；
由此得到往复惯性力为：
PI＝-mrωc 2(cosα+λcos2α)；
由往复惯性力所产生的力矩为：
其中，发动机固有的参数(m、L、r、α、β、λ、以及ωc)在实验前根据需要模拟的发动机的参数输入到惯量扭矩解算模型中。
3)建立燃烧扭矩解算模型，获得燃烧扭矩TC：
燃烧扭矩是由气缸内气体燃烧膨胀对活塞引起的力传递到曲轴和飞轮产生的周期性力矩，其燃烧压力为：
式中：Pcomp(αIVC)进气支管压力，Vd(αIVC)气缸容积，V(α)压缩过程中气缸容积(瞬时容积)，k为膨胀系数；
瞬时气缸工作容积为：
式中：为压缩比，Vc为燃烧室容积，D为气缸直径；
曲柄受到的单缸燃烧扭矩为活塞受到的燃烧压力乘以活塞行程随曲柄转角的变化率，即可得到：
式中：P(α)为燃烧压力，Patm大气压力，Ap气缸活塞面积，活塞运动速率；
其中，发动机固有的参数(Pcomp(αIVC)、Vd(αIVC)、k、Patm、以及Ap)在实验前根据需要模拟的发动机的参数输入到惯量扭矩解算模型中。
4)对平均扭矩TM，惯量扭矩TI，以及燃烧扭矩TC进行拟合，得到电动机的输入扭矩，即总的发动机扭矩Ten:
Ten＝TM+TI+TC。
由于上述过程获得的是单杠发动机的扭矩，因此，当发动机为多缸发动机时，首先通过步骤1)—步骤4)上述方法得到各单缸的扭矩，然后根据缸数和发火顺序进行扭矩叠加，从而获得发动机总的扭矩信号，参见图4，为4缸发动机的燃烧扭矩信号叠加。
现有传动系台架上发动机扭振信号一般由基于物理结构的发动机模型提供，由于这些模型需要大量的发动机开发数据或详细的发动机设计参数，这些参数在传动系台架试验中不易获得，而且传统的物理模型在模拟发动机燃烧特性时复杂的数学计算很难满足高动态控制响应的要求，因此这些模型不适用于传动系台架的发动机扭振模拟。而本发明的发动机扭振模拟方法，把发动机扭矩分为平均扭矩、惯性扭矩和燃烧扭矩，从而分别对平均扭矩、惯性扭矩和燃烧扭矩进行数学模型解算，并且各数学模型计算量小，从而有效节省了系统运行资源，解决了发动机周期扭振模拟过程中，对扭振模型的实时性和动态响应需求。
最后需要说明的是，以上实施例仅用以说明本发明的技术方案而非限制技术方案，本领域的普通技术人员应当理解，那些对本发明的技术方案进行修改或者等同替换，而不脱离本技术方案的宗旨和范围，均应涵盖在本发明的权利要求范围当中。
Claims (2)
1.一种汽车试验台架发动机扭振模拟方法，其特征在于：包括如下步骤：
1)根据发动机的固有参数，建立发动机“油门开度—转速—扭矩”模型，通过该“油门开度—转速—扭矩”模型能根据输入的油门开度或转速信号，获得与之对应的扭矩信号，即获得发动机的平均扭矩TM；
2)建立惯量扭矩解算模型，获得惯量扭矩TI：
发动机往复运动部件所产生的往复惯性力为：
<mrow>
<msub>
<mi>P</mi>
<mi>I</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mi>m</mi>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mo>;</mo>
</mrow>
式中：m为活塞和连杆的总质量，表示活塞加速度；
由于往复惯性力通过连杆作用在曲柄销上，在曲轴上产生周期性变化的力矩，从而引起轴系的惯量扭振，因此：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>x</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mo>+</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&beta;</mi>
<mo>+</mo>
<mi>r</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>cos</mi>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&alpha;</mi>
</mrow>
</msqrt>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
式中：x为活塞行程，L为连杆长度，r为曲柄半径，α为曲柄转角，β为曲柄销中心和曲轴旋转中心所在的直线与连杆轴心线之间的夹角，λ为曲柄半径与连杆长度比；
从而得到活塞加速度：
<mrow>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mo>=</mo>
<msup>
<msub>
<mi>r&omega;</mi>
<mi>c</mi>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&alpha;</mi>
<mo>+</mo>
<mi>&lambda;</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mn>2</mn>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
式中ωc为曲柄旋转角速度；
由此得到往复惯性力为：
PI＝-mrωc 2(cosα+λcos2α)；
由往复惯性力所产生的惯量力矩为：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mi>I</mi>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mi>I</mi>
</msub>
<mfrac>
<mrow>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>+</mo>
<mi>&beta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>cos</mi>
<mi>&beta;</mi>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mi>r</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>&ap;</mo>
<mo>-</mo>
<msup>
<mi>mr</mi>
<mn>2</mn>
</msup>
<msup>
<msub>
<mi>&omega;</mi>
<mi>c</mi>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>cos</mi>
<mi>&alpha;</mi>
<mo>+</mo>
<mi>&lambda;</mi>
<mi>cos</mi>
<mn>2</mn>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>sin</mi>
<mi>&alpha;</mi>
<mo>+</mo>
<mfrac>
<mi>&lambda;</mi>
<mn>2</mn>
</mfrac>
<mi>sin</mi>
<mn>2</mn>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msup>
<mi>mr</mi>
<mn>2</mn>
</msup>
<msup>
<msub>
<mi>&omega;</mi>
<mi>c</mi>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>&lambda;</mi>
<mn>4</mn>
</mfrac>
<mi>sin</mi>
<mi>&alpha;</mi>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mi>sin</mi>
<mn>2</mn>
<mi>&alpha;</mi>
<mo>-</mo>
<mfrac>
<mrow>
<mn>3</mn>
<mi>&lambda;</mi>
</mrow>
<mn>4</mn>
</mfrac>
<mi>sin</mi>
<mn>3</mn>
<mi>&alpha;</mi>
<mo>-</mo>
<mfrac>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
<mn>4</mn>
</mfrac>
<mi>sin</mi>
<mn>4</mn>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
3)建立燃烧扭矩解算模型，获得燃烧扭矩TC：
燃烧扭矩是由气缸内气体燃烧膨胀对活塞引起的力传递到曲轴和飞轮产生的周期性力矩，其燃烧压力为：
<mrow>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>m</mi>
<mi>p</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>I</mi>
<mi>V</mi>
<mi>C</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<msub>
<mi>V</mi>
<mi>d</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>I</mi>
<mi>V</mi>
<mi>C</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>k</mi>
</msup>
</mrow>
<mrow>
<mi>V</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mi>k</mi>
</msup>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
式中：Pcomp(αIVC)进气支管压力，Vd(αIVC)气缸容积，V(α)瞬时气缸工作容积，k为膨胀系数；
瞬时气缸工作容积为：
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>&pi;D</mi>
<mn>2</mn>
</msup>
</mrow>
<mn>4</mn>
</mfrac>
<mo>{</mo>
<mfrac>
<mi>x</mi>
<mrow>
<mi>&epsiv;</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mi>x</mi>
<mn>2</mn>
</mfrac>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mi>&lambda;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mi>&pi;</mi>
<mn>180</mn>
</mfrac>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>&lambda;</mi>
</mfrac>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>&lambda;</mi>
<mn>2</mn>
</msup>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>&pi;</mi>
<mn>180</mn>
</mfrac>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
<mo>&rsqb;</mo>
<mo>}</mo>
<mo>;</mo>
</mrow>
式中：为压缩比，Vc为燃烧室容积，D为气缸直径；
曲柄受到的单缸燃烧扭矩为活塞受到的燃烧压力乘以活塞行程随曲柄转角的变化率，即可得到：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mi>C</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>a</mi>
<mi>t</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<msub>
<mi>A</mi>
<mi>p</mi>
</msub>
<mo>&times;</mo>
<mfrac>
<mrow>
<mi>d</mi>
<mi>x</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>&alpha;</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>a</mi>
<mi>t</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>&pi;</mi>
<mfrac>
<msup>
<mi>D</mi>
<mn>2</mn>
</msup>
<mn>4</mn>
</mfrac>
<mo>&lsqb;</mo>
<mi>r</mi>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<msup>
<mi>r</mi>
<mn>2</mn>
</msup>
<mrow>
<mn>2</mn>
<mi>L</mi>
</mrow>
</mfrac>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
式中：P(α)为燃烧压力，Patm大气压力，Ap气缸活塞面积，活塞运动速率；
4)对平均扭矩TM，惯量扭矩TI，以及燃烧扭矩TC进行拟合，得到电动机的输入扭矩，即总的发动机扭矩Ten:
Ten＝TM+TI+TC。
2.根据权利要求1所述的汽车试验台架发动机扭振模拟方法，其特征在于：当发动机为多缸发动机时，首先通过步骤1)—步骤4)上述方法得到各单缸的扭矩，然后根据缸数和发火顺序进行扭矩叠加，从而获得发动机总的扭矩信号。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610006317.5A CN105651524B (zh) | 2016-01-06 | 2016-01-06 | 汽车试验台架发动机扭振模拟方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610006317.5A CN105651524B (zh) | 2016-01-06 | 2016-01-06 | 汽车试验台架发动机扭振模拟方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105651524A CN105651524A (zh) | 2016-06-08 |
| CN105651524B true CN105651524B (zh) | 2018-01-26 |
Family
ID=56490671
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610006317.5A Active CN105651524B (zh) | 2016-01-06 | 2016-01-06 | 汽车试验台架发动机扭振模拟方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105651524B (zh) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106096193A (zh) * | 2016-06-28 | 2016-11-09 | 重庆长安汽车股份有限公司 | 整车怠速、加速仿真发动机激励求解系统及方法 |
| CN106872174A (zh) * | 2017-02-22 | 2017-06-20 | 重庆理工大学 | 汽车变速器台架敲击试验发动机瞬态周期扭矩模拟方法 |
| CN110068463B (zh) * | 2019-03-26 | 2021-05-28 | 中国汽车技术研究中心有限公司 | 一种用于汽车动力总成的分离协同测试方法及系统 |
| CN112628001B (zh) * | 2019-09-24 | 2022-07-05 | 上海汽车集团股份有限公司 | 一种热机怠速的设定方法及装置 |
| CN113758557B (zh) * | 2020-06-05 | 2024-06-11 | 上海汽车集团股份有限公司 | 一种获取车辆结构振动噪声的方法及系统 |
| CN112557026B (zh) * | 2020-12-22 | 2023-03-14 | 中国汽车技术研究中心有限公司 | 模拟驱动方法、模拟负载方法及发动机齿轮试验系统 |
| CN112977394B (zh) * | 2021-02-18 | 2024-01-23 | 精进电动科技股份有限公司 | 一种抑制发动机扭矩脉动的方法和混合动力系统 |
| CN113389645B (zh) * | 2021-07-20 | 2022-11-08 | 中国第一汽车股份有限公司 | 一种混合动力汽车及消除起机抖动的控制方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6336070B1 (en) * | 2000-03-01 | 2002-01-01 | Ford Global Technologies, Inc. | Apparatus and method for engine crankshaft torque ripple control in a hybrid electric vehicle |
| CN102518521A (zh) * | 2011-12-27 | 2012-06-27 | 大连理工大学 | 一种车用发动机转矩估计方法 |
| CN103162966A (zh) * | 2011-12-09 | 2013-06-19 | 上海船舶运输科学研究所 | 一种柴油机推进系统 |
| CN103674571A (zh) * | 2014-01-03 | 2014-03-26 | 重庆市科学技术研究院 | 用于室内台架试验的汽车传动系载荷谱采集、编制方法 |
| CN104179573A (zh) * | 2014-08-06 | 2014-12-03 | 清华大学 | 一种内燃机分缸燃烧状态估计方法 |
-
2016
- 2016-01-06 CN CN201610006317.5A patent/CN105651524B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6336070B1 (en) * | 2000-03-01 | 2002-01-01 | Ford Global Technologies, Inc. | Apparatus and method for engine crankshaft torque ripple control in a hybrid electric vehicle |
| CN103162966A (zh) * | 2011-12-09 | 2013-06-19 | 上海船舶运输科学研究所 | 一种柴油机推进系统 |
| CN102518521A (zh) * | 2011-12-27 | 2012-06-27 | 大连理工大学 | 一种车用发动机转矩估计方法 |
| CN103674571A (zh) * | 2014-01-03 | 2014-03-26 | 重庆市科学技术研究院 | 用于室内台架试验的汽车传动系载荷谱采集、编制方法 |
| CN104179573A (zh) * | 2014-08-06 | 2014-12-03 | 清华大学 | 一种内燃机分缸燃烧状态估计方法 |
Non-Patent Citations (2)
| Title |
|---|
| 汽车传动系统动态载荷谱研究;丁小辉 等;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20140315(第3期);C035-37 * |
| 轿车动力总成悬置系统隔振性能匹配研究;李士保;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20100715(第07期);C035-12 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105651524A (zh) | 2016-06-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105651524B (zh) | 汽车试验台架发动机扭振模拟方法 | |
| CN106872174A (zh) | 汽车变速器台架敲击试验发动机瞬态周期扭矩模拟方法 | |
| CN104111650B (zh) | 一种双燃料发动机控制器仿真测试系统 | |
| JP2013163436A (ja) | モータによるエンジン変動抑制装置 | |
| JP6588165B2 (ja) | 内燃機関の運転中に使用される燃料の組成を求めるための方法 | |
| CN107476891A (zh) | 一种用于增程器中的发电机补偿扭矩的控制方法 | |
| CN107220487A (zh) | 一种柴油机轴系扭振计算方法 | |
| CN104179573B (zh) | 一种内燃机分缸燃烧状态估计方法 | |
| Xu et al. | Misfire fault diagnosis of range extender based on harmonic analysis | |
| CN103900820A (zh) | 发动机负荷检测装置以及发动机负荷检测方法 | |
| Li et al. | A test technology of a vehicle driveline test bench with electric drive dynamometer for dynamic emulation | |
| CN112818520B (zh) | 一种基于振动信号的fplg控制策略仿真模型构建方法及系统 | |
| Giakoumis et al. | Study of crankshaft torsional deformation under steady-state and transient operation of turbocharged diesel engines | |
| JP2007127458A (ja) | 動力伝達機構の試験装置 | |
| Roberts et al. | Detailed Powertrain Dynamics Modelling in Dymola-Modelica | |
| Campos et al. | Cylinder deactivation strategy for fuel consumption reduction | |
| Chen et al. | Development of engine model using modulization method for EMS verification through MIL and HIL | |
| Cakir | Experimental dynamic analysis of the piston assembly of a running single-cylinder diesel engine | |
| CN115099112B (zh) | 一种考虑轴系扭振与缸压耦合的内燃机表面振动预测方法 | |
| Bey et al. | Downsizing concept with two-cylinder CNG engine | |
| Togai et al. | A Reduced Order Powertrain Model Concept for Model Based Development Process with a Driver Agent | |
| Feng et al. | A transient dynamic model for HEV engine and its implementation for fuzzy-PID governor | |
| Johansson et al. | Analysis of engine cold start simulation in GT-power | |
| Winke et al. | Dynamic Simulation of Hybrid Powertrains using Different Combustion Engine Models | |
| Pohjalainen et al. | Novel crank mechanism increasing engine efficiency and reducing CO2 emissions |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |