基于MATLAB的发动机万有特性曲线绘制方法

2009年第2期(总第110期)内燃机与动力装置 I.C.E&Powerplant2009年4月

=设计研究>

基于MATLAB的发动机万有特性曲线绘制方法

周广猛,郝志刚,刘瑞林,陈 东,管金发,张春海

1

2

1

3

1

4

(1.军事交通学院汽车工程系,天津 300161;2.军事交通学院训练部,天津 300161;3.军事交通学院基础部,天津 300161;4.兰州军区军械汽车技工训练大队,陕西西安 710111) 摘要:利用MATLAB数学运算能力,处理发动机性能试验数据,方法简单可靠,绘制后的万有曲线直观明了,把等燃油消耗率曲线、外特性曲线和等功率曲线较好地拟合在同一张图上,拟合程度较高。

关键词:MATLAB;万有特性曲线;绘图

中图分类号:TK402 文献标识码:A 文章编号:1673-6397(2009)02-0034-03

UniversalCharacteristicsCurvePlottingMethodbasedonMATLAB

ZHOUGuang-meng,HAOZhi-gang,LIURui-lin,CHENDong,GUANJin-fa,ZHANGChun-hai

1

2

1

3

1

4

(1.AutomobileEngineeringDepartment,AcademyofMilitaryTransportation,Tianjin300161,China;

2.TrainingDepartment,AcademyofMilitaryTransportation,Tianjin300161,China;

3.GeneralCourseDepartment,AcademyofMilitaryTransportation,Tianjin300161,China;4.OrdnanceMechanicTrainingBrigade,LanZhouTheater,Xi.an710111,China))

Abstract:TakingadvantageofMATLABmathematicoperation,datafromenginecharacteristictestwasprocessed,themethodissimpleandcredible,Theuniversalcharacteristicscurveplottedisintuitionisticandperspicuous,andwasingoodfitwithdatagotintest.

KeyWords:MATLAB;UniversalCharacteristicsCurve;Plot

引 言

为了能全面反映发动机的性能,把发动机的多个参数画在一张图上而形成的多参数的特性曲线叫做发动机的万有特性曲线

[1]

单,方法较为成熟,研究的重点多集中在等燃油消耗率曲线的拟合与绘制上,MATLAB中提供了二元插值函数、三次插值、V4插值等多种插值方法,但要求数据间隔要足够小,而且对周围节点的精度要求高,个别数据点的误差有可能造成数据畸变。由洛阳凯迈机电有限公司开发的FST2D发动机控制系统利用MATLAB的计算引擎,可以利用发动机试验数据绘制万有特性曲线,但对试验点的密度和准确度仍然要求较高。而文中所采用的方法能够很好地解决这个问题。其它方法如神经网络拟合方法需要重新训练网格,比较繁琐;从一元样条非张量积形式推广到薄板样条形式可较好光滑曲线,但易出现多个插值点;本文中利用的最小二乘法原理,采用多元线性回归的方法,但随着拟合条件增加,也有产生奇点的可能。但在试验点采集受限等条件限制下仍

[6]

[5]

[4]

,传统用作图法制取万

有特性曲线是将不同转速下的负荷特性曲线绘制在同一张坐标图上,形成曲线簇,然后从曲线簇上把等油耗点逐一投影到万有特性图上,并圆滑地连接成等油耗曲线,再做出等功率曲线,画出外特性曲线,进而得到发动机的万有特性曲线,这种万有特性曲线的手工绘制方法费时费力,难以保证数据和图形的精度

[2]

,而MATLAB软件具有强大的矩阵计算和

[3]

数据可视化能力,为万有特性曲线的绘制提供了一种新的方法。国内开展了利用MATLAB进行万有特性曲线绘制的研究,由于外特性曲线拟合较为简

:(),男,,,

2009年第2期周广猛,等: 基于MATLAB的发动机万有特性曲线绘制方法

[7]

#35#

不失为一个较好的手段。

坐标的万有特性曲线运用最广,绘制该万有特性曲线需要画出等燃油消耗率曲线、等功率曲线和边界线,具体绘制方法如表1所示。

1 利用MATLAB绘制万有特性曲线的过程

以转速为横坐标、平均有效压力(或扭矩)为纵

表1 发动机万有特性曲线绘制方法

线型等燃油消耗率曲线边界线等功率曲线

绘制方法

先建立燃油消耗率与转速和扭矩的关系模型,绘出三维曲面图,再生成二维的等值线图。根据发动机外特性曲线的数据绘制边界线。公式Pe=TtqnP9550绘制,绘出三维曲面图,再生成二维等值线图

采用(建立)的模型be=f(Ttq,n)

拟合方式最小二乘法原理,多元线性回归样条插值拟合样条插值拟合

MATLAB命令(参数)meshgridmeshcontourNaNinterp1plotinterp1NaN

Ttq=f(n)Ttq=f(n)

发动机负荷特性试验的数据并不是矩阵格式,需要首先将它们转化为绘制MATLAB所识别的、绘

[8]

制三维图形的矩阵格式,才能绘制出转速n(rPmin)、有效扭矩Ttq(N#m)和燃油消耗率

be(gP(kW#h))的三维曲面图,再利用MATLAB语言里的contour语句生成二维的等值线图(相当于用不同的Ttq-n平面去截be、n、Ttq的三维曲面图),进而得到发动机的等油耗曲线,而等油耗曲线模型的建立实际上是建立燃油消耗率be与转速n和平均有效压力Pme的函数,采用多元线性回归的方法进行曲面拟合。

边界线的绘制实际上是把由外特性试验得到的bel

be2,

b

=

11,1

n1n2,nN

Ttq1Ttq2,TtqN

n1n2,nN

222

Ttq、n数据绘制在平面上,形成Ttq-n关系曲线,构成发动机万有特性曲线的边界线,此时油量调节机构固定在标定循环供油量位置,万有特性曲线上的转速和扭矩不可能超出边界线的范围。

根据公式Pe=TtqnP9550,选择不同的Ttq、n,利用interp1进行拟合,得到函数Pe=f(Ttq,n),绘制得到三维曲线,利用

contour命令得到的双曲线即是等功率曲线。

2 模型的建立

2.1 等油耗消耗率曲线模型的建立

根据多元线性回归理论,

建立模型如下:n1n1,n1

lll

n1Ttq1n2Ttq2

,NNTtqN

Ttq1Ttq2,TtqN

22

2

,,,

n1Ttq1n2Ttq2

,nNTtqN

l-1l-1

l-1

,,

Ttq1Ttq2,

ll

l

a0@

a1,ak-+

e0e1,e(1)

[3,9,10]

,T 可表示为:B=G@A+E,式中k为多项式的项数,l为多项式的最高次幂,其中k=(l+1)(l+2)P2,(1,2,,N),为不同的试验点,A=(a0,a1,,ak-1)为模型中的待定系数,E=(e0,e1,,eN)为随机误差,又称为残差。

2.2 边界线模型的建立

MATLAB中一维插值函数interp1(),提供了三种插值方法可选-linear.(线性的,此选项是默认的,它在两个点之间简单地采用直线拟合,故效果并不光滑),-cubic.(三次的)和-spline.(样条型)等,本论文采用拟合效果较好的-spline.型插值方法建立模型。

2.3 等功率曲线模型的建立

根据公式Pe=TtqnP9550,建立等功率曲线的模[10]

3 曲线的绘制

根据某发动机外特性和负荷特性得到的数据来绘制该发动机的万有特性曲线,外特性和负荷特性数据如表2、表3所示。

表2 外特性试验数据

转速n(rPmin)扭矩Ttq(N#m)转速n(rPmin)扭矩Ttq(N#m)

[1**********]509.5

[1**********]

[1**********]

[1**********]448.4

2102528.82700408.3

N

2199522.82802357.4

2

492.2471.2

由式(1),根据最小二乘法原理,J=i=E0ei=E#E值为最小,此时

T

9J

|=0。进而得到B==A

G@A,从理论上讲系数矩阵G的阶数l越大越能更好地实现等燃油消耗率曲线的拟合,然而l变大时,

#36#内燃机与动力装置2009年4月

函数可能出现病态,文章折衷选取了二次函数进行最小二乘拟合,此时拟合效果较好,得到be=f(Ttq,n)的函数,而由公式Pe=Ttq#nP9550得到Pe=f(Ttq,n),绘制燃油消耗率和功率的三维模型如

图1所示,利用contour命令绘制等油耗率曲线和等功率曲线的二维图,并利用外特性数据采用样条型插值方法绘制边界线,最终得到图2所示的万有特性曲线。

表3 负荷特性试验数据

n

Pr#min-

1

1400

TtqPN#m399.8354.1

318.5278.1236.2203.6185.3157.2117.280.8

2200

TtqPN#m420.7379.6334.6291.6244.4202.8157.5116.074.137.8

beg#(kW#h)-1

234.7259.8235.5237.6242.8292.3277.9308.7396.2605.9

TtpPN#m404.6360.5322.7283.0243.3205.5162.1124.786.852.4

beg#(kW#h)-1

222.8220.4232.4228.5227.8232.6248.5245.9272.4329.7

TtpPN#m409.1365.7328.3284.1243.7203.2164.3123.983.539.7

2400

beg#(kW#h)-1

174.2242.2252.1287.4253.6263.6290.6316.8378.0518.8

TtqPN#m378.0344.7310.3264.3226.1186.8154.2115.376.334.1

1600

beg#(kW#h)-1

222.0221.7235.4226.5230.5236.8249.1276.1407.9487.0

TtqPN#m408.3368.3328.3289.0244.4208.8167.7132.189.546.1

2600

beg#(kW#h)-1

256.9253.7253.5260.0303.8280.7300.6346.6435.6812.9

TtqPN#m315.6275.5242.5210.3178.5145.6118.672.652.822.4

1800

beg#(kW#h)-1

226.0225.3226.4233.9242.1283.3253.9271.4323.5468.6

TtqPN#m425.6380.3332.7290.9244.4205.1160.2114.568.830.7

2800

beg#(kW#h)-1

257.9295.3282.4288.7301.9329.7357.0475.4580.31080.12000

beg#(kW#h)-1

206.5231.1231.1233.0242.0244.9265.0299.8398.0596.8

图1 等燃油消耗率曲线和等功率曲线的三维拟合图

图2 万有特性曲线

(下转第48页)

#48#内燃机与动力装置2009年4月

性消声器进行消声实验,并且记录下相应数据;在保持同一工况转速的情况下,接通电源,有源消声部分开始工作,再次记录相应数据。根据摩托车发动机排气噪声的频谱特征及有源消声的特点,本实验系统主要针对基频噪声及其一次谐波噪声进行控制,从实验结果看其消声效果还是比较理想,由于消声器特有的消声方式,对高频部分也具备一定消声效果。基于扬声器参数的限制,下限暂定为125Hz。

表1 降噪的实验结果

fPHzLP1PdBLP2PdB$LPPdB

125

826616

250786414

500665115

100065605

200061592

400055541

(1)针对内燃机排气噪声的特点,提出了基于模拟电子技术的有源消声系统,该系统把有源消声和无源消声相结合,在不影响中高频消声效果的前提下,对低频段噪声消声效果明显。

(2)虽然受传声器、扬声器等器件好坏和非线性的影响,但在实验消声过程中采取了较好的匹配措施,信号处理部分采用纯模拟电路,实时性好,消声效果未受到很大影响,可以实现宽带消声。参考文献:

[1]葛维晶,王伟生,张宝欢.摩托车发动机排气消声器的研究与设计[J].小型内燃机,1994,23(3):37-41.

[2]冯振东,宋传学.车内噪声主动控制系统初探.汽车工程,1991,13(4):202-207.

[3]吴斌.发动机排气自适应有源消声关键技术的研究[D]:[学位论文].北京:北京工业大学机械工程与应用电子技术学院,2001,13-37.

[4]KrauseP,WeltensH.Advanceddesignofautomotiveexhaustsilencersystems[C].SAEPaper922088,1992.

[5]KimHeung-Seob,HongJin-Seok,SohnDong-Goo,etal.DevelopmentofanActiveMufflerSystemforReducingExhaustNoiseandFlowRestrictioninaHeavyVehicle[J].NoiseControlEngineeringJournal,1999,47(2):57-63.

表1中LP1指的是原有抗性消声器的消声效果,LP2指的是带有有源消声部分的总的消声效果,$LP指的是有源消声部分所达到的消声效果。在理想条件下,各个频率点的消声指数应为无穷大,但是由于幅值和相位的偏差,传声器、扬声器和电路中一些元器件的延时和非线性,对噪声消除能力都有不同程度的影响,使得末端还保留一定的残余噪声,另外消声器管道两端实际存在的声反射对实验结果也有一定影响。

2 结束语

(上接第36页)

4 结 论

(1)MATLAB强大的数学运算能力,方便实用的绘图功能为发动机万有特性曲线的绘制提供了一个很好的途径,利用MATLAB绘制万有曲线,作图方法简单,拟合程度高,提高了数据处理的速度和精度,为发动机动力性、经济性的研究,发动机参数的标定,发动机与传动系的匹配提供了更可靠的依据;

(2)与其它利用MATLAB绘制万有特性的方法相比,文章所述的方法能够很好地把等油耗曲线、等功率曲线及外特性曲线拟合在一张图上,图像更加直观可靠;

(3)该绘图方法可以利用MATLAB图形用户界面设计技术进行进一步处理,进而完成该绘图方法的图形界面设计。同时可以进一步探索利用更好的曲线拟合方法,或采取分段拟合方法,进一步提高拟合程度。参考文献:

[1]孙凤英,汽车性能与使用技术[M].北京:机械工业出版社,

2002.

[2]易雪梅,吴伶.用MATLAB语言绘制发动机万有特性的两种方法[J].北京汽车,2005(5):33-35.

[3]李金辉,徐立友.基于MATLAB语言的发动机特性研究[J].汽车科技,2005(3):40-42.

[4]倪计民,陈源.MATLAB软件在汽车动力传动系匹配计算中的应用[J],上海汽车,2003:24-27.

[5]洛阳凯迈机电有限公司,FST2D系统使用说明书[Z].洛阳.[6]李小华,罗福强,汤东.多项式插值法绘制发动机万有特性曲线[J].农业工程学报,2004,20(5):138-141.

[7]董敬.汽车拖拉机发动机[M].北京:机械工业出版社,2004,149.

[8]苏金明,阮沈勇.MATLAB实用教程[M].北京:电子工业出版社,2005.

[9]CurtisF.Gerald,PatrickO.Wheatley.AppliedNumericalAnalysis[M].Beijing:HigherEducationPress,2006.

[10]关志伟,杨玲,施继红.基于MATLAB语言的发动机万有特性研究[J].吉林农业大学学报,2003,25(3):339-342.[11]薛定宇,陈阳泉.基于MATLABPSimulink的系统仿真与应用[M].北京:清华大学出版社.2002.

2009年第2期(总第110期)内燃机与动力装置 I.C.E&Powerplant2009年4月

=设计研究>

基于MATLAB的发动机万有特性曲线绘制方法

周广猛,郝志刚,刘瑞林,陈 东,管金发,张春海

1

2

1

3

1

4

(1.军事交通学院汽车工程系,天津 300161;2.军事交通学院训练部,天津 300161;3.军事交通学院基础部,天津 300161;4.兰州军区军械汽车技工训练大队,陕西西安 710111) 摘要:利用MATLAB数学运算能力,处理发动机性能试验数据,方法简单可靠,绘制后的万有曲线直观明了,把等燃油消耗率曲线、外特性曲线和等功率曲线较好地拟合在同一张图上,拟合程度较高。

关键词:MATLAB;万有特性曲线;绘图

中图分类号:TK402 文献标识码:A 文章编号:1673-6397(2009)02-0034-03

UniversalCharacteristicsCurvePlottingMethodbasedonMATLAB

ZHOUGuang-meng,HAOZhi-gang,LIURui-lin,CHENDong,GUANJin-fa,ZHANGChun-hai

1

2

1

3

1

4

(1.AutomobileEngineeringDepartment,AcademyofMilitaryTransportation,Tianjin300161,China;

2.TrainingDepartment,AcademyofMilitaryTransportation,Tianjin300161,China;

3.GeneralCourseDepartment,AcademyofMilitaryTransportation,Tianjin300161,China;4.OrdnanceMechanicTrainingBrigade,LanZhouTheater,Xi.an710111,China))

Abstract:TakingadvantageofMATLABmathematicoperation,datafromenginecharacteristictestwasprocessed,themethodissimpleandcredible,Theuniversalcharacteristicscurveplottedisintuitionisticandperspicuous,andwasingoodfitwithdatagotintest.

KeyWords:MATLAB;UniversalCharacteristicsCurve;Plot

引 言

为了能全面反映发动机的性能,把发动机的多个参数画在一张图上而形成的多参数的特性曲线叫做发动机的万有特性曲线

[1]

单,方法较为成熟,研究的重点多集中在等燃油消耗率曲线的拟合与绘制上,MATLAB中提供了二元插值函数、三次插值、V4插值等多种插值方法,但要求数据间隔要足够小,而且对周围节点的精度要求高,个别数据点的误差有可能造成数据畸变。由洛阳凯迈机电有限公司开发的FST2D发动机控制系统利用MATLAB的计算引擎,可以利用发动机试验数据绘制万有特性曲线,但对试验点的密度和准确度仍然要求较高。而文中所采用的方法能够很好地解决这个问题。其它方法如神经网络拟合方法需要重新训练网格,比较繁琐;从一元样条非张量积形式推广到薄板样条形式可较好光滑曲线,但易出现多个插值点;本文中利用的最小二乘法原理,采用多元线性回归的方法,但随着拟合条件增加,也有产生奇点的可能。但在试验点采集受限等条件限制下仍

[6]

[5]

[4]

,传统用作图法制取万

有特性曲线是将不同转速下的负荷特性曲线绘制在同一张坐标图上,形成曲线簇,然后从曲线簇上把等油耗点逐一投影到万有特性图上,并圆滑地连接成等油耗曲线,再做出等功率曲线,画出外特性曲线,进而得到发动机的万有特性曲线,这种万有特性曲线的手工绘制方法费时费力,难以保证数据和图形的精度

[2]

,而MATLAB软件具有强大的矩阵计算和

[3]

数据可视化能力,为万有特性曲线的绘制提供了一种新的方法。国内开展了利用MATLAB进行万有特性曲线绘制的研究,由于外特性曲线拟合较为简

:(),男,,,

2009年第2期周广猛,等: 基于MATLAB的发动机万有特性曲线绘制方法

[7]

#35#

不失为一个较好的手段。

坐标的万有特性曲线运用最广,绘制该万有特性曲线需要画出等燃油消耗率曲线、等功率曲线和边界线,具体绘制方法如表1所示。

1 利用MATLAB绘制万有特性曲线的过程

以转速为横坐标、平均有效压力(或扭矩)为纵

表1 发动机万有特性曲线绘制方法

线型等燃油消耗率曲线边界线等功率曲线

绘制方法

先建立燃油消耗率与转速和扭矩的关系模型,绘出三维曲面图,再生成二维的等值线图。根据发动机外特性曲线的数据绘制边界线。公式Pe=TtqnP9550绘制,绘出三维曲面图,再生成二维等值线图

采用(建立)的模型be=f(Ttq,n)

拟合方式最小二乘法原理,多元线性回归样条插值拟合样条插值拟合

MATLAB命令(参数)meshgridmeshcontourNaNinterp1plotinterp1NaN

Ttq=f(n)Ttq=f(n)

发动机负荷特性试验的数据并不是矩阵格式,需要首先将它们转化为绘制MATLAB所识别的、绘

[8]

制三维图形的矩阵格式,才能绘制出转速n(rPmin)、有效扭矩Ttq(N#m)和燃油消耗率

be(gP(kW#h))的三维曲面图,再利用MATLAB语言里的contour语句生成二维的等值线图(相当于用不同的Ttq-n平面去截be、n、Ttq的三维曲面图),进而得到发动机的等油耗曲线,而等油耗曲线模型的建立实际上是建立燃油消耗率be与转速n和平均有效压力Pme的函数,采用多元线性回归的方法进行曲面拟合。

边界线的绘制实际上是把由外特性试验得到的bel

be2,

b

=

11,1

n1n2,nN

Ttq1Ttq2,TtqN

n1n2,nN

222

Ttq、n数据绘制在平面上,形成Ttq-n关系曲线,构成发动机万有特性曲线的边界线,此时油量调节机构固定在标定循环供油量位置,万有特性曲线上的转速和扭矩不可能超出边界线的范围。

根据公式Pe=TtqnP9550,选择不同的Ttq、n,利用interp1进行拟合,得到函数Pe=f(Ttq,n),绘制得到三维曲线,利用

contour命令得到的双曲线即是等功率曲线。

2 模型的建立

2.1 等油耗消耗率曲线模型的建立

根据多元线性回归理论,

建立模型如下:n1n1,n1

lll

n1Ttq1n2Ttq2

,NNTtqN

Ttq1Ttq2,TtqN

22

2

,,,

n1Ttq1n2Ttq2

,nNTtqN

l-1l-1

l-1

,,

Ttq1Ttq2,

ll

l

a0@

a1,ak-+

e0e1,e(1)

[3,9,10]

,T 可表示为:B=G@A+E,式中k为多项式的项数,l为多项式的最高次幂,其中k=(l+1)(l+2)P2,(1,2,,N),为不同的试验点,A=(a0,a1,,ak-1)为模型中的待定系数,E=(e0,e1,,eN)为随机误差,又称为残差。

2.2 边界线模型的建立

MATLAB中一维插值函数interp1(),提供了三种插值方法可选-linear.(线性的,此选项是默认的,它在两个点之间简单地采用直线拟合,故效果并不光滑),-cubic.(三次的)和-spline.(样条型)等,本论文采用拟合效果较好的-spline.型插值方法建立模型。

2.3 等功率曲线模型的建立

根据公式Pe=TtqnP9550,建立等功率曲线的模[10]

3 曲线的绘制

根据某发动机外特性和负荷特性得到的数据来绘制该发动机的万有特性曲线,外特性和负荷特性数据如表2、表3所示。

表2 外特性试验数据

转速n(rPmin)扭矩Ttq(N#m)转速n(rPmin)扭矩Ttq(N#m)

[1**********]509.5

[1**********]

[1**********]

[1**********]448.4

2102528.82700408.3

N

2199522.82802357.4

2

492.2471.2

由式(1),根据最小二乘法原理,J=i=E0ei=E#E值为最小,此时

T

9J

|=0。进而得到B==A

G@A,从理论上讲系数矩阵G的阶数l越大越能更好地实现等燃油消耗率曲线的拟合,然而l变大时,

#36#内燃机与动力装置2009年4月

函数可能出现病态,文章折衷选取了二次函数进行最小二乘拟合,此时拟合效果较好,得到be=f(Ttq,n)的函数,而由公式Pe=Ttq#nP9550得到Pe=f(Ttq,n),绘制燃油消耗率和功率的三维模型如

图1所示,利用contour命令绘制等油耗率曲线和等功率曲线的二维图,并利用外特性数据采用样条型插值方法绘制边界线,最终得到图2所示的万有特性曲线。

表3 负荷特性试验数据

n

Pr#min-

1

1400

TtqPN#m399.8354.1

318.5278.1236.2203.6185.3157.2117.280.8

2200

TtqPN#m420.7379.6334.6291.6244.4202.8157.5116.074.137.8

beg#(kW#h)-1

234.7259.8235.5237.6242.8292.3277.9308.7396.2605.9

TtpPN#m404.6360.5322.7283.0243.3205.5162.1124.786.852.4

beg#(kW#h)-1

222.8220.4232.4228.5227.8232.6248.5245.9272.4329.7

TtpPN#m409.1365.7328.3284.1243.7203.2164.3123.983.539.7

2400

beg#(kW#h)-1

174.2242.2252.1287.4253.6263.6290.6316.8378.0518.8

TtqPN#m378.0344.7310.3264.3226.1186.8154.2115.376.334.1

1600

beg#(kW#h)-1

222.0221.7235.4226.5230.5236.8249.1276.1407.9487.0

TtqPN#m408.3368.3328.3289.0244.4208.8167.7132.189.546.1

2600

beg#(kW#h)-1

256.9253.7253.5260.0303.8280.7300.6346.6435.6812.9

TtqPN#m315.6275.5242.5210.3178.5145.6118.672.652.822.4

1800

beg#(kW#h)-1

226.0225.3226.4233.9242.1283.3253.9271.4323.5468.6

TtqPN#m425.6380.3332.7290.9244.4205.1160.2114.568.830.7

2800

beg#(kW#h)-1

257.9295.3282.4288.7301.9329.7357.0475.4580.31080.12000

beg#(kW#h)-1

206.5231.1231.1233.0242.0244.9265.0299.8398.0596.8

图1 等燃油消耗率曲线和等功率曲线的三维拟合图

图2 万有特性曲线

(下转第48页)

#48#内燃机与动力装置2009年4月

性消声器进行消声实验,并且记录下相应数据;在保持同一工况转速的情况下,接通电源,有源消声部分开始工作,再次记录相应数据。根据摩托车发动机排气噪声的频谱特征及有源消声的特点,本实验系统主要针对基频噪声及其一次谐波噪声进行控制,从实验结果看其消声效果还是比较理想,由于消声器特有的消声方式,对高频部分也具备一定消声效果。基于扬声器参数的限制,下限暂定为125Hz。

表1 降噪的实验结果

fPHzLP1PdBLP2PdB$LPPdB

125

826616

250786414

500665115

100065605

200061592

400055541

(1)针对内燃机排气噪声的特点,提出了基于模拟电子技术的有源消声系统,该系统把有源消声和无源消声相结合,在不影响中高频消声效果的前提下,对低频段噪声消声效果明显。

(2)虽然受传声器、扬声器等器件好坏和非线性的影响,但在实验消声过程中采取了较好的匹配措施,信号处理部分采用纯模拟电路,实时性好,消声效果未受到很大影响,可以实现宽带消声。参考文献:

[1]葛维晶,王伟生,张宝欢.摩托车发动机排气消声器的研究与设计[J].小型内燃机,1994,23(3):37-41.

[2]冯振东,宋传学.车内噪声主动控制系统初探.汽车工程,1991,13(4):202-207.

[3]吴斌.发动机排气自适应有源消声关键技术的研究[D]:[学位论文].北京:北京工业大学机械工程与应用电子技术学院,2001,13-37.

[4]KrauseP,WeltensH.Advanceddesignofautomotiveexhaustsilencersystems[C].SAEPaper922088,1992.

[5]KimHeung-Seob,HongJin-Seok,SohnDong-Goo,etal.DevelopmentofanActiveMufflerSystemforReducingExhaustNoiseandFlowRestrictioninaHeavyVehicle[J].NoiseControlEngineeringJournal,1999,47(2):57-63.

表1中LP1指的是原有抗性消声器的消声效果,LP2指的是带有有源消声部分的总的消声效果,$LP指的是有源消声部分所达到的消声效果。在理想条件下,各个频率点的消声指数应为无穷大,但是由于幅值和相位的偏差,传声器、扬声器和电路中一些元器件的延时和非线性,对噪声消除能力都有不同程度的影响,使得末端还保留一定的残余噪声,另外消声器管道两端实际存在的声反射对实验结果也有一定影响。

2 结束语

(上接第36页)

4 结 论

(1)MATLAB强大的数学运算能力,方便实用的绘图功能为发动机万有特性曲线的绘制提供了一个很好的途径,利用MATLAB绘制万有曲线,作图方法简单,拟合程度高,提高了数据处理的速度和精度,为发动机动力性、经济性的研究,发动机参数的标定,发动机与传动系的匹配提供了更可靠的依据;

(2)与其它利用MATLAB绘制万有特性的方法相比,文章所述的方法能够很好地把等油耗曲线、等功率曲线及外特性曲线拟合在一张图上,图像更加直观可靠;

(3)该绘图方法可以利用MATLAB图形用户界面设计技术进行进一步处理,进而完成该绘图方法的图形界面设计。同时可以进一步探索利用更好的曲线拟合方法,或采取分段拟合方法,进一步提高拟合程度。参考文献:

[1]孙凤英,汽车性能与使用技术[M].北京:机械工业出版社,

2002.

[2]易雪梅,吴伶.用MATLAB语言绘制发动机万有特性的两种方法[J].北京汽车,2005(5):33-35.

[3]李金辉,徐立友.基于MATLAB语言的发动机特性研究[J].汽车科技,2005(3):40-42.

[4]倪计民,陈源.MATLAB软件在汽车动力传动系匹配计算中的应用[J],上海汽车,2003:24-27.

[5]洛阳凯迈机电有限公司,FST2D系统使用说明书[Z].洛阳.[6]李小华,罗福强,汤东.多项式插值法绘制发动机万有特性曲线[J].农业工程学报,2004,20(5):138-141.

[7]董敬.汽车拖拉机发动机[M].北京:机械工业出版社,2004,149.

[8]苏金明,阮沈勇.MATLAB实用教程[M].北京:电子工业出版社,2005.

[9]CurtisF.Gerald,PatrickO.Wheatley.AppliedNumericalAnalysis[M].Beijing:HigherEducationPress,2006.

[10]关志伟,杨玲,施继红.基于MATLAB语言的发动机万有特性研究[J].吉林农业大学学报,2003,25(3):339-342.[11]薛定宇,陈阳泉.基于MATLABPSimulink的系统仿真与应用[M].北京:清华大学出版社.2002.