基于MATLAB的发动机万有特性曲线绘制方法图文精.docx

1、基于MATLAB的发动机万有特性曲线绘制方法图文精=设计研究 基于 MATLAB 的发动机万有特性曲线绘制方法周广猛 1, 郝志刚 2, 刘瑞林 1, 陈 东 3, 管金发 1, 张春海4(1. 军事交通学院 汽车工程系 , 天津 300161; 2. 军事交通学院 训练部 , 天津 300161; 3. 军事交通学院 基础部 , 天津 300161; 4. 兰州军区 军械汽车技工训练大队 , 陕西 西安 710111 摘要 :利用 MATLAB 数学运算能力 , 处理发动机性能试验数据 , 方法简单可靠 , 绘制后的万 有曲线直观明了 , 把等燃油消耗率曲线 、 外特性曲线和等功率曲线较好地

2、拟合在同一张图上 , 拟合程度较高 。关键词 :MATLAB; 万有特性曲线 ; 绘图中图分类号 :TK402 文献标识码 :A 文章编号 :1673-6397(2009 02-0034-03Universal Characteristics C urve Plotting Method based on MATLABZ HO U Guang-me ng 1, H AO Zhi-gang 2, LI U Rui-lin 1, C HE N Dong 3, GUA N Jin-fa 1, Z HANG C hun-hai4(1. Automobile Engineering Department

3、, Academy of Military Transportation, Tianjin 300161, China;2. Training Department, Acade my of Military Transportation, Tianjin 300161, China;3. General Course Department, Academy of Military Transporta tion, Tianjin 300161, China; 4. Ordnance Mechanic Training Brigade, Lan Zhou Theater, Xi . an 71

4、0111, China Abstract :Taking advantage of MATLAB mathematic operation, data from engine characteristic test was processed, the method is simple and credible, The universal characteristics curve plotted is intuitionistic and perspicuous, and was in good fit with data got in test.Key W ords :MATLAB; U

5、niversal Characteristics Curve; Plot:( , 男 , , , 引 言为了能全面反映发动机的性能 , 把发动机的多 个参数画在一张图上而形成的多参数的特性曲线叫 做发动机的万有特性曲线1, 传统用作图法制取万有特性曲线是将不同转速下的负荷特性曲线绘制在 同一张坐标图上 , 形成曲线簇 , 然后从曲线簇上把等 油耗点逐一投影到万有特性图上 , 并圆滑地连接成 等油耗曲线 , 再做出等功率曲线 , 画出外特性曲线 , 进而得到发动机的万有特性曲线 , 这种万有特性曲 线的手工绘制方法费时费力 , 难以保证数据和图形 的精度2, 而 MATLAB 软件具有强大的矩阵

6、计算和数据可视化能力 3, 为万有特性曲线的绘制提供了 一种新的方法。国内开展了利用 MATLAB 进行万有 特性曲线绘制的研究 , 由于外特性曲线拟合较为简单 , 方法较为成熟 , 研究的重点多集中在等燃油消耗 率曲线的拟合与绘制上 , MATLAB 中提供了二元插 值函数、 三次插值、 V4插值等多种插值方法 , 但要求 数据间隔要足够小 , 而且对周围节点的精度要求高 , 个别数据点的误差有可能造成数据畸变 4。由洛阳 凯迈机电有限 公司开发的 FST2D 发动 机控制系统 利用 MATLAB 的计算引擎 , 可以利用发动机试验数 据绘制万有特性曲线 , 但对试验点的密度和准确度 仍然要

7、求较高 5。而文中所采用的方法能够很好地 解决这个问题。其它方法如神经网络拟合方法需要 重新训练网格 , 比较繁琐 ; 从一元样条非张量积形式 推广到薄板样条形式可较好光滑曲线 , 但易出现多 个插值点 6; 本文中利用的最小二乘法原理 , 采用多 元线性回归的方法 , 但随着拟合条件增加 , 也有产生 奇点的可能。但在试验点采集受限等条件限制下仍2009年第 2期 (总第 110期 内 燃机与动力装置 I. C. E &Powerplan t 2009年 4月不失为一个较好的手段。1利用 MATLAB 绘制万有特性曲线的过程 以转速为横坐标、 平均有效压力 (或扭矩 为纵 坐标的万有特性曲线

8、运用最广 7, 绘制该万有特性 曲线需要画出等燃油消耗率曲线、 等功率曲线和边 界线 , 具体绘制方法如表 1所示。表 1发动机万有特性曲线绘制方法线型 绘制方法 采用 (建立 的模型 拟合方式 MATLAB 命令 (参数 等 燃 油 消 耗 率 曲线 先建 立 燃 油消 耗 率 与转 速 和 扭矩的 关 系 模 型 , 绘 出 三 维 曲 面图 , 再生成二维的 等值线图。be=f (Ttq, n最小 二乘 法原 理 , 多元 线性回归meshgrid mesh contour NaN边界线 根据 发 动 机外 特 性 曲线 的 数 据绘制边界线。Ttq=f (n 样 条插值拟合 inter

9、p1plot等功率曲线 公式 Pe=Ttqn P 9550绘制 , 绘出三维曲面图 , 再生成二维等值线图Ttq=f (n 样 条插值拟合 interp1NaN发动机负荷特性试验的数据并不是矩阵格式 , 需要首先将它们转化为绘制 MATLAB 所识别的、 绘 制 三 维 图 形 的 矩 阵 格 式 8, 才 能 绘 制 出 转 速n (r P min 、 有 效 扭 矩 T tq (N #m 和 燃 油 消 耗 率b e (gP (kW #h 的三维曲面图 , 再利用 MATLAB 语言 里的 contour 语句生成二维的等值线图 (相当于用不 同的 T tq -n 平面去截 b e 、 n

10、 、 T tq 的三维曲面图 , 进而 得到发动机的等油耗曲线 , 而等油耗曲线模型的建 立实际上是建立燃油消耗率 b e 与转速 n 和平均有 效压力 P me 的函数 , 采用多元线性回归的方法进行曲 面拟合。边界线的绘制实际上是把由外特性试验得到的 T tq 、 n 数据绘制在平面上 , 形成 T tq -n 关系曲线 , 构 成发动机万有特性曲线的边界线 , 此时油量调节机 构固定在标定循环供油量位置 , 万有特性曲线上的 转速和扭矩不可能超出边界线的范围。根据公式 P e =T tq n P 9550, 选择不同的 T tq 、 n , 利 用 interp1进行拟合 , 得到函数

11、P e =f (T tq , n , 绘制 得到三维曲线 , 利用 contour 命令得到的双曲线即是 等功率曲线。2模型的建立2. 1等油耗消耗率曲线模型的建立根据多元线性回归理论 , 建立模型如下 :b elb e2, b=1n 1T tq 1n 21n 1T tq1T 2tq1, n l 1n l -11T tq 1, T l tq 11n 2T tq 2n 22n 2T tq2T 2tq2, n l 1n l -12T tq 2, T l tq 2, , , , , , , , ,1n N T tqN n 2N N N T tq N T 2tq N , n l 1n l -1N T

12、tqN , T la 0a 1,a k -+e 0e 1,e(1 3, 9, 10可表示为 :B =G A +E , 式中 k 为多项式的项数 , l 为多项式的最高次幂 , 其中 k =(l +1 (l +2 P 2, (1, 2, , N , 为不同的试验点 , A =(a 0, a 1, ,a k-1 为模型中的待定系数 , E =(e 0, e 1, , e N 为随机误差 , 又称为残差。2. 2边界线模型的建立MATLAB 中一维插值 函数 interp1( , 提供 了三种插值方法可选 -linear . (线性的 , 此选项是默认的 ,它在两个点之间简单地采用直线拟合 , 故效

13、果并不光滑 , -cubic . (三次的 和 -spline . (样条型 等 , 本论文采用拟合效果较好的 -spline . 型插值方法 建立模型 10。2. 3等功率曲线模型的建立根据公式 P e =T tq n P 9550, 建立等功率曲线的模3曲线的绘制根据某发动机外特性和负荷特性得到的数据来 绘制该发动机的万有特性曲线 , 外特性和负荷特性 数据如表 2、 表 3所示。表 2外特性试验数据转速 n (r P mi n 140315971797198621022199扭矩 Ttq(N#m 474497515526528. 8522. 8转速 n (r P mi n 2303240

14、02507259827002802扭矩 T tq (N#m 509. 5492. 2471. 2448. 4408. 3357. 4由 式 (1 , 根据 最 小二 乘法 原 理 , J =E N i =0e 2i = E T #E 值为最 小 , 此时9J|=A =0。进 而得 到 B = G A , 从理论上讲系数矩阵 G 的阶数 l 越大越能更 好地实现等燃油消耗率曲线的拟合 , 然而 l 变大时 , # 35 #2009年第 2期 周广猛 , 等 :基于 MATLAB 的发动机万有特性曲线绘制方法函数可能出现病态 , 文章折衷选取了二次函数进行 最小 二 乘 拟 合 , 此 时 拟 合

15、 效 果 较 好 , 得 到 b e =f (T tq , n 的函数 , 而由公式 P e =T tq #n P 9550得到 P e =f (T tq , n , 绘制燃油消耗率和功率的三维模型如图 1所示 , 利用 contour 命令绘制等油耗率曲线和等 功率曲线的二维图 , 并利用外特性数据采用样条型 插值方法绘制边界线 , 最终得到图 2所示的万有特 性曲线。表 3 负荷特性试验数据n P r #min - 11400160018002000T tq P N #m b e g #(k W #h -1T tp P N #m b e g #(kW #h -1T tq P N #m b

16、e g #(kW #h -1T tq P N #m b e g #(k W #h -1399. 8222. 8409. 1222. 0408. 3226. 0425. 6206. 5354. 1220. 4365. 7221. 7368. 3225. 3380. 3231. 1318. 5232. 4328. 3235. 4328. 3226. 4332. 7231. 1278. 1228. 5284. 1226. 5289. 0233. 9290. 9233. 0236. 2227. 8243. 7230. 5244. 4242. 1244. 4242. 0203. 6232. 6203.

17、2236. 8208. 8283. 3205. 1244. 9185. 3248. 5164. 3249. 1167. 7253. 9160. 2265. 0157. 2245. 9123. 9276. 1132. 1271. 4114. 5299. 8117. 2272. 483. 5407. 989. 5323. 568. 8398. 080. 8329. 739. 7487. 046. 1468. 630. 7596. 82200240026002800T tq P N #m b e g #(k W #h -1T tp P N #m b e g #(kW #h -1T tq P N #m

18、 b e g #(kW #h -1T tq P N #m b e g #(k W #h -1420. 7234. 7404. 6174. 2378. 0256. 9315. 6257. 9379. 6259. 8360. 5242. 2344. 7253. 7275. 5295. 3334. 6235. 5322. 7252. 1310. 3253. 5242. 5282. 4291. 6237. 6283. 0287. 4264. 3260. 0210. 3288. 7244. 4242. 8243. 3253. 6226. 1303. 8178. 5301. 9202. 8292. 320

19、5. 5263. 6186. 8280. 7145. 6329. 7157. 5277. 9162. 1290. 6154. 2300. 6118. 6357. 0116. 0308. 7124. 7316. 8115. 3346. 672. 6475. 474. 1396. 286. 8378. 076. 3435. 652. 8580. 337. 8605. 952. 4518. 834. 1812. 922. 41080. 1图 1 等燃油 消耗率曲线和等功率曲线的三维拟合图图 2 万有特性曲线(下 转第 48页 #36#内燃机与动力装置 2009年 4月性 消声器进行消声实验 , 并且

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

21、 LP1指的是原有抗性消声器的消声效果 , LP2指的是带 有有源消声部分的总的消声效果 , $LP 指 的 是有源消声部分所达到的消声效果。在理想条件 下 , 各个频率点的消 声指数应为无穷大 , 但是由于幅值 和 相位的偏差 , 传声 器、 扬声器和电路中一些元器件的 延 时和非线性 , 对噪 声消除能力都有不同程度的影响 , 使得末端还保留一定的残余噪声 , 另外消声器管道两 端 实际存在的声反射对实 验结果也有一定影响。2 结束语(1针对内燃机排气噪声的特点 , 提出了基于模拟 电子技术的有源消声系统 , 该系统把有源消声和无源 消声相结合 , 在不影响中高频消声效果的前提下 , 对低

22、 频段噪声消声 效果明显。(2虽然受传声器、 扬声器等器件好坏和非线性的 影响 , 但 在实验消声过程中采取了较好的匹配措施 , 信 号处理部分采 用纯模拟电路 , 实时性好 , 消声效果未受 到很大影响 , 可以实现宽带消声 。 参考文献 :1葛维晶 , 王伟生 , 张宝欢 . 摩托车发动机排气消声器的研究 与设计 J.小型内燃机 , 1994, 23(3 :37-41.2冯振东 , 宋 传学 . 车内噪声主动控 制系统初探 . 汽 车工程 , 1991, 13(4 :202-207.3吴斌 . 发动机排气自适应有源消声关键技术的研究 D:学 位论文 .北京 :北京工业 大学机械工程与应 用

23、电子技术学 院 , 2001, 13-37.4Krause P, Welten s H. Ad vanced design of automotive exhaust silencer systemsC.SAE Paper 922088, 1992.5Kim Heung -Seob, Hong Jin-Seok, Sohn Dong -Goo, et al. Develop men t of an Acti ve Muffler System for Reduci ng E xhaust Noise and Flow Restriction in a Heavy Vehicle J.Noise

24、 Control Engineering Journal,1999, 47(2 :57-63.(上接第 36页 4 结 论(1MATLAB 强大 的数学运算能力 , 方便实用的绘 图 功能为发动机万有特性 曲线的绘制提供了一个很好 的 途径 , 利用 MA TLAB 绘制万有曲线 , 作图方法简单 , 拟 合程度高 , 提高了 数据处理的速度和精度 , 为发动机 动 力性、 经济性的研究 , 发动机参数的标定 , 发动机与 传 动系的匹配提供了更可 靠的依据 ;(2与其它利用 MA TLAB 绘制万有特性的方法相 比 , 文章所述的方法 能够很好地把等油耗曲线、 等功率 曲线及外特性曲线拟合在

25、一张图上 , 图像更加直观可 靠 ;(3该绘图方法可以利用 MA TLAB 图形用户界面 设计技术进行进一步处理 , 进而完成该绘图方法的图 形 界面设计。同时可以 进一步探索利用更好的曲线拟 合 方法 , 或采取分段 拟合方法 , 进一 步提高拟合程度。 参 考文献 :1孙凤英 , 汽车性能与使用技术 M.北京 :机械工业出 版社 ,2002.2易雪梅 , 吴伶 . 用 MA TLAB 语言绘制发动机万有特性的两种 方法 J. 北京汽车 ,2005(5 :33-35.3李金辉 , 徐立友 . 基于 MA TLAB 语言的发动机特性研究 J.汽 车科技 , 2005(3 :40-42.4倪计民

26、 , 陈源 . MA TLAB 软件在汽车动力传动系匹配计算中的 应用 J, 上海汽车 ,2003:24-27.5洛阳凯迈机电有限公司 , FS T2D 系统使用说明书 Z.洛阳 . 6李小华 , 罗福强 , 汤东 . 多项式插值法绘制发动机万有特性曲 线 J.农业工程学报 , 2004,20(5:138-141.7董敬 . 汽车拖拉机发动机 M.北京 :机械工业出版社 , 2004, 149.8苏金明 , 阮沈勇 .MA TLAB 实用教程 M.北京 :电子工业出版 社 , 2005.9Curtis F.Gerald, Patrick O. Wheatley. Applied Nu merical Analysi s M.Beijing:HigherEducation Press, 2006.10关志伟 , 杨玲 , 施继红 . 基于 MA TLAB 语言的发动机万有特 性研究 J.吉林农业大学学报 , 2003,25(3 :339-342. 11薛定宇 , 陈阳泉 . 基于 MA TLAB P Simulink 的系统仿真与应用 M. 北京 :清华大学出版社 . 2002.#48#内燃机与动力装置 2009年 4月