万有特性曲线绘制软件实现方法研究

《装备制造技术）２０１０年第９期　万有特性曲线绘制软件实现方法研究　蔡隆玉，谢继鹏，杨敏　（南京理工大学紫金学院，江苏南京２１００４６）　摘要：运用科学计算工具ＭＡＴＬＡＢ，分别基于传统负荷特性作图法和曲面到曲线的方法，实现了万有特性曲线的绘制方法，比较了两　种方法的优缺点．探讨了作图法在应用中的问题，及曲面拟合方法的优势后，进一步用ＭＡＴＬＡＢ自有的ＧＵＩ编辑工具，实现了具有　人机交互界面的万有特性绘制软件。　关键词：万有特性，负荷特性，作图法，曲面拟合，ＭＡＴＬＡＢ　ＧＵＩ　中图分类号：ＴＰ３９１．７　文献标识码：Ａ　文章编号：１　６７２—５４５Ｘ（２０１　０）０９—００４７—０３　内燃机工作特性有多种表现形式，其中负荷特性和速度　在同一张图纸上绘出特性曲线；　／ｋＰａ　特性只能用丁表示某一转速或某一节气门（或油门）时内燃机　经济性、动力性指标随负荷或转速变化的规律。而内燃机的工　－　况变化范围很大，为全面反映发动机的性能，便于与工作机械　；　‘匹配，需要绘制出能够表达内燃机多参数特性的万有特性图　曲线图。万有特性曲线是以转速ｎ为横坐标，以平均有效压力　ｐｍｅ（或转矩Ｔｔｑ）为纵坐标的坐标平面内的一系列特性参数的　等值曲线族，如反映燃油经济性的等油耗线和反映动力性的　等功率曲线…。　§　委　＼　０　：　萎量　ｌ一　ｉ　一一　ｆ　传统万有特性曲线的绘制方法有负荷特性法和速度特性　法，是根据负荷特性曲线或速度特性曲线，运用手工作图法得　＝：：　二　３２０　｝　ｌ　２８Ｏ　ｌ　｝　ｉ　２４０　２∞　到万有特性图的。随着计算机运算能力的提高，计算机高级程　序语言的应用，特别是科学计算工具ＭＡＴＬＡＢ的出现，为万有　特性曲线的自动绘制提供了可能。为此，人们提出了各种根据　曲线拟合函数实现万有特曲线绘制的方法，如最小二乘法拟　６ｃ，Ｉ　ｇ，（ｋＷ・’　Ｈ　ｎ／（ｒ／ａｒｉｎ）　图１万有特性负荷特性作图法　（２）根据内燃机工作转速范围，标出万有特性横坐标１１的　合｝１】、多项式插值法［３】、多元线性回归【　等。但还为见运算程序软　件化研究。　标尺，纵坐标平均有效压力Ｐ　（或转矩　）的标尺与整理得到　的负荷特性上的平均有效压力ｐｍｅ（或转矩　）的标尺相同；　本文首先基于传统的万有特性曲线绘制方法，探讨了运　用科学计算工具ＭＡＴＬＡＢ实现负荷特性法绘制万有特性曲线　的具体途径。分析其应用优缺点后，进一步通过二维插值方法　实现了内燃机特性参数三维曲面的拟合，并得到其二维等油　（３）将负荷特性旋转９０。后，置于万有特性纵坐标轴的左　侧，使同样是平均有效压力Ｐ　（或转矩　）的两个坐标对齐；　（４）在负荷特性图上，引若干条等燃油消耗率线与６ｅ线　相交，每条线各有ｌ～２个交点；　（５）从每一个交点引水平线至万有特性上，与负荷特性线　相同转速的位置上，获得若干新交点，并在每一交点上，标注　耗曲线。比较了汽油机、柴油机在曲线拟合时的不同特点及各　自适用的方法。并选择了曲面拟合方法实现了具有人机交互　界面的万有特性绘制软件。　出燃油消耗率的数值；　（６）所有转速下的负荷特性，都经过这样的转换后，依次　１　传统绘制方法及其数字化实现途径　如前所述，传统的万有特性曲线，是根据负荷特性曲线或　将６ｅ值相等的点连成光滑曲线，即可得到万有特性上的等燃　油消耗率６ｅ线。　依据作图法的过程，运用ＭＡＴＬＡＢ编写计算及绘图程序，　可以将手工绘图的过程通过软件自行实现。具体的程序框图　如图２所示。　速度特性曲线通过作图法得到的。以负荷特性法为例（如图１　所示）㈣。　１．１　绘图方法和步骤　本文根据一组柴油机的万有特性测试数据，依作图法得　到了其万有特性图（图３）。　（１）将各种转速下的负荷特性以平均有效压力Ｐ　（或转　矩　）为横坐标，以有效油耗率６ｅ为纵坐标，以同一比例尺，　收稿日期：２０１０－－０６—１４　作者简介：蔡隆玉（１９８２一），女，江西赣州人，助教，硕士学位，主要研究方向：发动机排放控制，工作过程和性能研究。　４７　Ｅｑｕｉｐｍｅｎｔ　Ｍａｎｕｆａｃｔｒｉｎｇ　Ｔｅｃｈｎｏｌｏｇｙ　Ｎｏ．９，２０１０　耗线变形、多条油耗线交叉的可能性更大。因此，绘制万有特　性曲线前，应通过样条插值拟合，观察负荷特性曲线走向，对　其进行直观分析，消除奇异点后，再进一步完成万有特性曲线　的绘制。　（３）由于一维插值的局限性，对于等油耗线存在封闭曲线　的情况，在曲线拟合时存在困难。关键是当封闭曲线出现在介　于两个负荷特性实验测试转速区域间时，难以合理估计其曲　线走向，拟合出合理的封闭等油耗线。　２　三维曲面拟合到二维等油耗线的实现　运用作图法绘制万有特性曲线，取值点较为可靠，但计算　程序较为繁琐，数据读取和油耗线数组相关的处理、计算过程　难以实现完全的程序自运行。以转速　为　坐标，转矩　为　】，坐标，燃油消耗率为ｚ坐标，运用二维样条插值函数　（ｉｎｔｅｒｐ２）将负荷特性测试数据进行三维曲面拟合（如图４所　示）。再用一系列平行与　—ｙ平面的平面（即一系列不同油耗　值的平面）截取三维曲面投影至　—ｌ，平面，即得到等燃油消　耗率线（如图５所示）。　ｆｑ≥蒿１日　博　罂霹嚣　图２作图法程序框图　卿　鲫　瞄　伽　瑚瑚　３蝴　扭矩Ｔ｜ｑ（Ｎ　ｒｎ）　。　１０∞　转速ｎ　０，ｍ　ｃｎ）　图４负荷特性三维曲面　’蝴’　’　’埘　湖’　建　２４００　。　。。　图３作图法得到的某柴油机万有特性图　１．２作图法绘制万有特性图时应注意的问题　笔者在程序调试过程中，还尝试了作图法实现汽油机万　有特性图的绘制，结果并不理想。在这一过程中，总结出运用　作图法绘制万有特性图时应注意的问题：　（１）由于原始数据为离散的测试结果点，对于内燃机而　；　言，油耗率值与转速、转矩之间无特定的函数关系。要通过有　限的测试数据，获得的更多的发动机工况点，选择正确的插值　蓦　方式，适当的求根初始值是关键。为使插值拟合结果接近测试　结果，可选用三次样条插值方式。而求根计算在选取初始值　时，则应首先观察负荷特性曲线的拟合结果，预估根的大小，　选取根附近的较小区间中的一个值作为求根初始值，以免遗　漏可能根或求根无解返回随机值，干扰后续等油耗线的拟合。　（２）运用作图法时，测试值数据量及其转速、扭矩的分　布，对拟合结果有很大影响。如果测试数据不理想，造成等油　４８　图５　等高线法获得的万有特性曲线　《装备制造技术）２ｏｌｏ年第９期　比较图３、图５可知，通过传统作图法得到的曲线，由于拟　软件运行后，首先，浏览按一定顺序初步整理好的测试数　合时采用了三次样条拟合，所得拟合曲线较为光滑，拟合误差　据文件；为提高计算效率，直接输入发动机负荷特性实验的各　较小。而二维插值拟合后的结果，其曲面较为光顺，等高线投　转速变化范围；最后，点击“显示曲线”按钮，程序将自动完成　影所得曲线却并不十分光滑。但两种不同方法所得万有特性　一系列的计算和绘图过程。万有特性［ｆｆ１线即在图形窗口呈现　等油耗线的变化趋势相同，工况范同较为一致。说明Ｆｈ曲面到　（如图７所示）。　曲线的方法，也能正确反映发动机的特性。而这一计算过程　中，可充分利用数组特点管理计算结果，实现数据的自动处　４　结束语　理。因此采用曲面到曲线的拟合方式更适用于软件。　（１）运用ＭＡＴＬＡＢ工具，采用作图法和曲面到曲线的万　３　万有特性曲线绘制软件的实现　有特性曲线绘制方法，各有其优缺点。作图法取点后曲线趋势　更直观可辨，而曲面拟合方法二维插值结果更为合理，而且程　ＭＡＴＬＡＢ具备编辑图形用户界面（Ｇｒａｐｈｉｃａｌ　Ｕｓｅｒ　Ｉｎｔｅｒｆａｃｅ，　序自动实现更为方便。　ＧＵＩ）的功能。在命令窗口中输入ｇｕｉｄｅ，即出现用户界面编辑　（２）本文所采用的绘制程序，主要为某一柴油机的万有特　窗口蜘。在窗口中，加入文本、按钮、图形的对象，即得到人机交　性测试数据而调试所得。为使运算及绘图程序更具普遍性，应　互的基本界面（如图６所示）。随后，完成各文本、按钮等对象　利用更多不同机型发动机测试结果进行调试。考虑到不同发　的ｃａｌｌ—ｂａｃｋ响应设置。在“显示曲线”部分，加入将曲面到曲　动机工作情况变化范围大，也可以针对不同类型发动机总结　线的拟合程序。　及应用不同的绘制程序。　（３）ＭＡＴＬＡＢ　ＧＵＩ的实现简单、直观，但要得到最终的可　执行程序，还需要进一步的转换。通过ＭＡＴＬＡＢ　ＧＵＩ得到的可　执行程序，还不能完全脱离ＭＡＴＬＡＢ工作环境，可移植性较　差。如将该绘图程序、软件应用于测功机厂家作为配套软件，　可在配套软件中加入ＭＡＴＬＡＢ的工作环境安装包。而对于发　动机生产厂家，如要脱离ＭＡＴＬＡＢ环境进行计算，则需要进一　步开发出基于Ｃ语言环境的计算软件。　图６软件初始界面　参考文献：　［１］周龙保，刘巽俊，高宗英．内燃机学［ＭＩ．北京：机械工业出版社，　２ｏｏ８．　ｆ２］周广猛，郝志刚，刘瑞林，陈东，管金发，张春海．基于ＭＡＴＬＡＢ　语言的发动机万有特性曲线绘制方法【Ｊ】．洛阳工学院学报，２００９，　（１１０）：３４—３６，４８．　【３】杜慧勇，李民，刘建新，赵致和．ＭＡＴＬＡＢ语言在内燃机万有特性　研究中的应用［Ｊ］．洛阳工学院学报，２００１，２２（４）：３２—３４．　【４］马向平，李春燕，骆清国，王宪成，张志远．一种发动机万有特性曲　面拟合的新方法【Ｊ】．装甲兵工学院学报，２００６，２０（１）：５２—５４．　【５】吴建华，常绿，韩同群，周炜．汽车发动机原理【Ｍ］．北京：机械工　业出版社，２００５．　图７软件运行结果界面　［６］罗华飞．ＭＡＴＬＡＢ　ＧＵＩ设计学习手ｉＥ［Ｍ］．北京：北京航空航天大　学出版社，２００９．　Ｒｅｓｅａｒｃｈ　ｏｎ　Ｓｏｆｔｗａｒｅ　ｆｏｒ　Ｕｎｉｖｅｒｓａｌ　Ｃｈａｒａｃｔｅｒｉｓｔｉｃ　Ｃｕｒｖｅ　Ｍａｐｐｉｎｇ　ＣＡＩ　Ｌｏｎｇ—ｙｕ，ＸＩＥ　Ｊｉ－ｐｅｎｇ，ＹＡＮＧ　Ｍｉｎ　（ＮａｎｊｉｎｇＵｎｉｖｅｒｓｉｔｙ　ｏｆＳｃｉｅｎｃｅ　ａｎｄＴｅｃｈｎｏｌｏｇｙＺｉｊｉｎＣｏｌｌｅｇｅ，Ｎａｎｊｉｎｇ　２１００４６，Ｃｈｉｎａ）　Ａｂｓｔｒａｃｔ：Ｉｎ　ｔｈｉｓ　ｐａｐｅｒ，ｓｃｉｅｎｔｉｆｉｃ　ｃｏｍｐｕｔｉｎｇ　ｔｏｏｌｓ　ＭＡＴＬＡＢ　ｉｓ　ｕｓｅｄ　ｔｏ　ｃｏｍｐ￣ｅｎｇｉｎｅ　ｕｎｉｖｅｒｓａｌ　ｃｈａｒａｃｔｅｒｉｓｔｉｃ　ｍａｐｐｉｎｇ．Ｔｒａｄｉｔｉｏｎａｌ　ｍａｐ－－　ｐｉｎｇ　ｍｅｔｈｏｄ　ｂａｓｅｄ　ｏｎ　ｌｏａｄ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ａｎｄ　ｔｈｅ　ｓｕｒｆａｃｅ　ｔｏ—ｃｕｎ．ｅ　ｍｅｔｈｏｄ　ａｌｅ　ｒｅｓｐｅｃｔｉｖｅｌｙ　ａｐｐｌｉｅｄ　ｆｏｒ　ｄａｔａ　ｃｏｍｐｕｔｅ　ａｎｄ　ｃｕｒｖｉｎｇ．　Ａｄｖａｎｔａｇｅｓ　ａｎｄ　ｄｉｓａｄｖａｎｔａｇｅｓ　ｏｆ　ｔｈｅ　ｔｗｏ　ｍｅｔｈｏｄｓ　ａｇｅ　ｃｏｍｐａｒｅｄ；ｔｈｅ　ｐｍｂｌｅｍｓ　ｏｆ　ｔｒａｄｉｔｉｏｎａｌ　ｍａｐｐｉｎｇ　ａｎｄ　ｔｈｅ　ａｄｖａｎｔａｇｅｓ　ｏｆ　ｓｕｒｆａｃｅ　ｆｉｔｔｉｎｇ　ｍｅｔｈｏｄｓ　ａｒｅ　ｂｏｔｈ　ｄｉｓｃｕｓｓｅｄ．Ｆｕｒｔｈｅｒ，ＭＡＴＬＡＢ　ｏｗｎ　ＧＵＩ　ｅｄｉｔｉｎｇ　ｔｏｌ　ｉｓ　ｉｍｐｌｅｍｅｎｔｅｄ　ｔｏ　ｄｅｖｅｌｏｐ　ｄｒａｗｉｎｇ　ｓｏｆｔｗａｒｅ　ｗｉｔｈ　ｔｈｅ　ｈｕｍａｎ—ｃｏｍｐｕｔｅｒ　ｉｎｔｅｒａｃｔｉｏｎ　ｉｎｔｅｒｆａｃｅ．　Ｋｅｙ　ｗｏ￣ｓ：ｕｎｉｖｅｒｓａｌ　ｃｈａｒａｃｔｅｒｉｓｔｉｃ；ｌｏａｄ　ｃｈａｒａｃｔｅｒｉｓｔｉｃ；ｍａｐｐｉｎｇ；ｓｕｒｆａｃｅ　ｆｉｔｔｉｎｇ；ＭＡＴＬＡＢ；ＧＵＩ　４９