发动机工作过程计算.docx

发动机工作过程计算.docx

《发动机工作过程计算.docx》由会员分享，可在线阅读，更多相关《发动机工作过程计算.docx（22页珍藏版）》请在冰点文库上搜索。

发动机工作过程计算

发动机工作过程计算

————————————————————————————————作者：

————————————————————————————————日期:

ﻩ

第2组

西华大学实验报告

开课学院及实验室:

　　　　　　　　　实验时间：

20１2年　　４月　７　日

学生姓名

学　号

实　验成绩

课程名称

发动机工作过程的计算

课程代码

实验项目名称

发动机工作过程计算上机实验

项目代　码

指导教师

项　目学　分

实验课考勤

实验工作表现

实验报告

1、实验目的

2、实验设备、仪器及材料

3、实验内容

３.1一般实验（非上机实验）:

3.1.1实验方案设计与选择（设计性、创新性实验需写该项内容）

3.1.2实验原理及实验步骤（实验工作原理或实验的主要操作过程）

3.1.3实验记录（核心代码及调试过程）

3.2上机实验：

3.2．１上机实验的内容及要求

3．2．2算法设计思想与算法实现步骤

3.2.３程序核心代码,程序调试过程中出现的问题及解决方法

3.２.4程序运行的结果

注解:

理工科实验需记录实验过程中的数据、图表、计算、现象观察等,实验过程中出现的问题；其它如在计算机上进行的编程、仿真性或模拟性实验需记录程序核心代码以及程序在调式过程中出现的问题及解决方法;记录程序执行的结果。

4、实验总结

4．１实验结果分析及问题讨论

4.２实验总结心得体会

注解:

实验总结的内容根据不同学科和类型实验要求不一样，一般理工科类的实验需要对实验结果进行分析,并且对实验过程中问题进行讨论;在计算机上进行的编程、仿真性或模拟性实验需要对上机实践结果进行分析，上机的心得体会及改进意见。

其它实验应总结实验过程写出心得体会及改进意见。

说明:

各门实验课程实验报告的格式及内容要求,请按照实验指导书的要求手工书写。

1.实验目的

2.实验仪器,设备及材料

3.实验内容

　MsｇBox＂西华大学2０09级交通与汽车工程学院　热能与动力工程（汽车发动机）3班唐俊

Eｎd　Suｂ

运行结果依次如下:

　进入界面时:

Private　SｕｂComｍaｎd1_Click（）

　　　Ｆoｒm1.Show

Ｆorm２.Ｈide

EndSub

运行出的图像如下：

PrivateSｕbCｏmmaｎd２＿Ｃliｃk（）

MsｇBoｘ＂西华大学200９级交通与汽车工程学院热能与动力工程（汽车发动机）　　３班　唐俊

EnｄSｕb

当进入主界面之后可对发动机类型，是否需要增压,经验参数以及发动机参数进行选择和修改。

选择界面如下:

主程序的详细代码如下：

DimεAsSinglｅ,Ｓ　AｓSinｇle,nAsInteger,nhＡsInteger，αAs　Singlｅ,vc,vh,pa,　n１,　λ,pz,vz,　n2,jｘ,　zy

PriｖａteSubCommand１_Cｌicｋ（）

Pｉctｕre１.Scalｅ（-０.05,－0.5）-（1．05,11）

　　Picture１.Ｌinｅ（0,　1０）-（1,10）　　　　　　　　　　　＇画X轴

　Pｉctuｒe1.Ｌine（0.98,9.9）-（1,１0）

Ｐictｕre1.Line（0.98,　１０．１）－（1,　1０）

　Picture1．CｕｒreｎtX=0.９8：

Picｔuｒe1.ＣurrentY　=1０.2:

　Picture１.Ｐｒint"V/L"

Fｏri=　０To0.9Step0.1

　　Picｔuｒe１.Ｌine（i,９．９）-（i，10）

　　Picture1.Priｎt"0"　＆i

　　Ｎeｘｔ　i

　Picture１.Lｉne　（0,10）-（0，0）　　　　　　　　　　　'画Y轴

　Ｐｉctｕre1.Ｌinｅ（0,0）－（－0.0１,０.3）

　　Ｐicturｅ１.Line　（0．01,０.3）-（０,　0）

　Ｐictｕrｅ1.CurrｅnｔＸ=-０.０4:

Pｉcｔure1．CurrentY　=　-0.4:

Pictｕre１.Print"P／Mpa＂

　Ｆｏri=9To　1Sｔep-1

　Piｃture１.Ｌine（０．01,ｉ）-（０,i）

　　Pｉｃture1.Ｐriｎt10-ｉ

　Ｎｅxti

Picture1.FｏrｅColor=vbRｅd　　　　　　　　　　　'画P0线

　　　Ｐｉcｔure１.Line（0,9.9）-（１,9.9）

　　Picture1.CuｒrentX＝1:

Piｃｔure1.CurｒｅntY　=9.4：

　Piｃｔｕｒｅ1.Priｎｔ"P０＂

　Picture1.ＦｏｒｅＣolｏr　＝vbBluｅ　　　　　　＇画Vc线　'压缩曲线

　　　ε=Vaｌ（Ｔｅxt2.Ｔext）

　ｄ=Tｅxt4.Tｅxｔ

　S=　Text3.Text

　ｒ＝Ｓ/　2

vｈ＝3.14＊d　^２*　S　/　4＊１0^-６

Ｔｅxt１６．Teｘｔ　=vh

ｖc=ｖh/（ε-　１）

Ｐictｕｒe1．CurｒentX=vc+　０.01：

Pｉｃｔure１.ＣurrｅntY=1:

Pｉcture１．Prｉnt　"Vc"

Pｉcｔure1.Ｌiｎｅ（vc,１0）－（vc，　0）

　　Va=vｃ+ｖh

l=Teｘt１5．Ｔeｘt

　　Ｐicture1.ForeColor=　ｖbGreen

　Forφx＝　３.14To　0　Sｔｅp　－0．0001

　　vcx=　vh／2＊（（１　-Cｏs（φx））　＋　（1-Ｃｏs（2＊φｘ））*　Rl/　4）+vｃ

　　　'vｃx=（3.1４　*d　^２＊r　*（（1-Cos（φx）+r　*（1-Cｏｓ（2*φx）　/4　/l）））/　４）＊10^　-８　＋　vc

　　ｐcx=pa*（Va/vｃｘ）^　ｎ1

　Picture１．PＳet（vｃx,１0-　pｃｘ）

Nexｔφｘ　　　　　　　　　　＇燃烧过程曲线

　　Pｉｃture1.Lｉne（vｃｘ,１0-pcx）-（vcx，１0　-ｐcx*　λ）　　　　＇膨胀过程曲线

Ｆorφx＝０Tｏ3.１4Ｓtep０.0０0１

ｖｂx＝vh／2*（（1－Cos（φx））+（1　－Cos（2　＊φx））＊　Ｒl　／4）　+　vｃ

　　　　　'ｖcx=（3.14*ｄ^　２＊　r*（（１－Coｓ（φx）＋r*（１-Cos（2　＊　φx）/4/　l）））/４+vc）　*　10^-6

　pbx=pz*（vc／vbx）　^ｎ2

　Ｐiｃturｅ１．ＰSet（vｂx,　10－　pｂx）

　Nextφｘ

　Ｐiｃture1．Line（ｖbx,10-pｂx）－（ｖbx,１0－pa）

Pｉcture１.ForｅＣolor=　ｖｂBlaｃk

EnｄSub

Pｒiｖａｔe　Ｓｕb　Coｍｍａnd２_Clｉcｋ（）

　Picture1.Cｌｓ

　ε　=Vａｌ（Text2．Texｔ）

　　ｄ=Texｔ4．Teｘt

　S=　Ｔｅxｔ3．Texｔ

ｒ=　S/2

　l=Text15.Text

　ｖh=　3．14*d　^2*S／　4＊１0^-６

　Tｅｘt16．Texｔ＝vｈ　　　　　　　　＇排气过程

　n=Vaｌ（Text５.Ｔeｘt）

　　nh=Vａｌ（Ｔext６.Ｔeｘt）

　If　zy=１Theｎ

　　ｐr　=0.1+0．03＊n/nh

　ＥｌseIfｚy=２Then

pk=Ｔext４5.Text

　　Δpｒ=　２　*10＾-8*ｎ^2

　pｒ=ｐk　+Δpr

ＥndIf

Tｅxt18.Ｔeｘｔ　＝pr

　　α　=　Ｖａl（Ｔexｔ1．Text）

ε　=Val（Text2.Texｔ）

　　Ｉfε=1Then

　　MsgＢox（"压缩比不能为“１”,请重新输入!

＂）

　　Tｅxｔ2.Text=　""

　　End　Ｉf

　Tr　=3５0/（１．2　/Lｏg10（ｎ）　+　0．0０５　＊（ε-　３）+０.01*　（α　-　1））'调用自定义函数就以10为底的对数

Text17．Texｔ=Tｒ　　　　　　　　'进气过程

ΔT=１0　＊　（１１0-　0．0125＊　n）　／（１１0－0.０12５*nh）　　　'1０为ΔTn,选取

φ=０.75

　　δ=０.5

　k=　１．３5　　　　　　　　　　　　　　'a=ｖh/　f

ｐ0=Ｔeｘt7.Teｘt'Pａ=　P0*　（1　-ｎ^2/2＊9.8*ｋ*p０/φ^2*（（ε-　δ）/（ε-1））＾2）　^k　／（k　-１）

　pa　=Ｉnｔ（１0　^4　＊　（ｐ0*（1-（（１０　＊　n^2/520／　1０＾6）/　１0.6　^2／０.7５　＾　2）*　（ε－　δ）^2/（ε-1）＾2）　^（k/　（k　-１））））／1０　^4

Ｔeｘt２0.Ｔext＝ｐa

　　t0=　Tｅxt8.Text

　γ=　Int（1０^　4　*（t0+ΔＴ）　*　ｐr/Tr/（ε　*　pa-pr））/10　^4　'残余废气系数

　　Tｅｘｔ21.Texｔ　=γ

tａ　=　（t0+ΔT+γ*φ＊Ｔr*（ｐａ/pｒ）^　（（n　-　１）/n））/（1　＋γ*　φ）

　　Ｔeｘｔ19.Tｅxｔ＝ta

ηv　＝　ε*pa　*ｔ0/（ε　-　1）／p0/ｔa/（１+　γ）　　　　　　＇充气效率

Text22．Teｘt＝ηv　　　　　　　　　　　　'压缩过程

　n1=1.４6　-　０.05*　nｈ　/n

　Teｘt25.Text=n１

　　pｃ=　pa　＊ε　^　n1

Tｅxt29.Teｘt=pc

tｃ=　ta＊　ε^（ｎ1－1）

　　Tｅxt2８．Text=tc

　'燃烧过程：

　gｃ=Val（Text９.Texｔ）

　gh=Val（Text10.Ｔeｘt）

　ｇo　=Ｖaｌ（Teｘt11.Teｘt）

　　gs=Vａｌ（Ｔｅxt12.Ｔeｘt）

　ｌ0　＝1/0．2１*（gc／12+ｇｈ/　４－go　/　32）　　　　　　'理论空气量

　　ｍt＝Ｔext1３．Tｅxt

Ｉfjx＝1Then

　m1=α　*　ｌ０＋１　/mt　　　　'汽油机新鲜充量

　ElseIｆｊｘ=2Then

　m1=α＊　ｌ０

ＥnｄIｆ

　　　If　α　＞1Then　　　　　　　　　　　＇m2燃烧产物总量

　　　m２=　α＊l0　+　gｈ/　4+　go　/　32

ElseＩｆα<　１　Tｈeｎ

　　m2=０.7９　*α　*l0+ｇc　/1２+　gh/　2

　　Ｅnｄ　Ｉf

　μ0=m2／m1　　　　　　　　　　　'理论分子变更系数

　　μ=（μ0＋　γ）/　（1+　γ）　　　　　　＇实际分子变更系数

Tｅxｔ3４.Text＝μ

　Iｆα<1　Tｈen

　　Δhu＝5８000*（1-α）　　　　　　　　＇化学损失

　End　If

　＇Tｚ的计算

　ζ=Ｔext23.Ｔｅｘt

　Ifjx　＝1Then

　hu=44100

　EｌｓｅIｆjｘ＝　2Thｅｎ

hu　=４２5００

EnｄIｆ

ｑ1　=ζ　*（hu-Δhu）　　　　　　　　　＇加热工质的热量

　cv1　=　４.８15+　0．4１5*10　^　-3　*（tc+t0）　/2

　　Ｓｅlｅct　Caseα

　　Case０.８Tｏ1

　　　　Ifｊx=1　Then

　　　　　aa=（3．7+3.３　＊　α）*10^-4　*　μ

　　　　　　　bb＝（4.4+0.62＊α）*μ

　ｃc=　-（q1／ｍ1　/（1　+γ）+cｖ1*tc）

　tz=　（－bb+　（bb　^2-　4*aa*ｃc）^0.5）　／2/aa／2.５

　　　Ｔext31.Tｅxt　=ｔz

　　　　λ=μ*　tz／tc

　　　Text33．Tｅxt＝λ

　　pz　=　λ　*pｃ

　　　Teｘt32.Text=pｚ

　　EｌｓeIf　ｊｘ＝　2Then

　　λ=Teｘt３3.Teｘt

　　　ａa=μ*　（３.7+3．3*α）*　１0　^－4

　　　　　　　　　ｂb　=　（6.385+0.６２　*α）＊μ

　　　ｃc=　-（ζ＊ｈu/ｍ1/　（1+　γ）+（cv1＋1.98*λ）*tc）

　　ｔz＝（-bb+（ｂb　^2　－4　*aa　*cｃ）^0．5）/　2　/aa　／　2．５

　　　　　Text３1.Ｔeｘｔ=tz

　　pz＝λ*ｐc

　　Text3２．Text=pz

　　　ρ　＝μ＊tz　/　λ　/tc

　　　　Texｔ３０.Teｘｔ=ρ

　　　EｎdIf

　　Ｃasｅ　1Ｔｏ２

　Iｆjｘ　=1　Thｅｎ

　　　　　　ａa　=（３.7+3.3/α）*　10^-4　*μ

　　bb　=（4.8+0.22/α）*μ

　　　　cc＝－（q1/m1/（1+γ）　+cv１*ｔc）

　　　　ｔz=　（-bb　+（bb^2　－4　＊aａ*　cc）＾0.5）/2/aa　／　2．5

　　　Texｔ３１．Text=tz

　　　　λ=μ　＊ｔz／ｔc

　　　　Teｘｔ３3.Ｔexｔ＝　λ

　　　　　pz=λ*ｐc

　　Ｔeｘt３2.Text　=pz

　　　ElsｅIfjx　＝2Theｎ

　λ=Ｔext３３.Tｅxt

　　　　　　　ａa=　（3.7+3．3/　α）*10^　-4＊μ

　　bb=（6.785+0.22　／α）　*μ

　　cｃ　=　-（ζ*hｕ　/m1/（1+γ）+　（cv1　+1．98*λ）*tc）

　tz=（-bb＋　（bb　^2　-　4＊ａa　*cc）　^　0.5）/2　／aａ　/2．5

　　Text３1．Text=tｚ

　　pz＝λ＊pc

　　　Teｘt32.Teｘt　=　pｚ

　　ρ=μ*tz/　λ/ｔｃ

　　　　　Text３0.Ｔeｘt=　ρ

　　　　EndIf

EnｄSelect

　'膨胀过程:

　n2=　１．2+　０.039　*nh／ｎ

　Texｔ26．Tｅｘt＝ｎ2

　'δ=ε/　ρ

　'Ｔext37．Teｘt=δ

　　pb=pz/ε＾n2

Tｅxt36.Ｔext＝pb

　tb　=　ｔz／ε^　（ｎ2-1）

　Tｅxt３５.Ｔext＝tb

＇发动机性能分析

ρ=1

Text3０.Teｘｔ=ρ

δ　=0

　　Ｔext37.Ｔext=　δ

　　　Pii=　ｐｃ＊（λ*　（１-１／ε^　（n2　-　1））　／　（n２　-　1）　-（１-１/　ε^　（n1　-１））／（n1－1））/　（ε-1）

　＇Pｉi=　pa*　ε^ｎ1　*（λ*（ρ－1）　+λ*ρ/　（n2　-1）*（1　-1/　δ^（ｎ2　－1））　-　（１-1/ε^（n1－１）/（n1－　１）））　／　（ε-1）

φi=Text2４．Texｔ

　　pｉ=φi*Piｉ

Text3９.Teｘt=pi

　ηｉ=8．3１4＊m1*t0　*pi　／hu/　ｐ０/ηv

　Teｘt38.Teｘt　＝ηｉ

　gi＝　3.６*1０＾6/　hu／ηi

Text4０.Tｅxt=ｇi

　ｐm=0.０４*（1+n/1０００）

　ηｍ＝1-pm/pi

　Text27.Tｅｘt=ηm

　ηｅ=ηi＊ηm

　Teｘt42.Text　=ηe

　ｐe=ｐｉ　*ηｍ

　Ｔexｔ4１.Ｔｅxt=pe

gｅ　=3．６＊１0　^6／　ｈu　/ηe

Text43．Teｘt　=ｇｅ

i　=Tｅxt1４.Texｔ

nｅ=１0　*i*　ｖh*　pｅ*n　/1２０0

　Teｘｔ44.Ｔext=ne

EnｄSub

Private　SubCommaｎd３_Cｌick（）

　MｓgBox　（"确定退出?

"）

　　Eｎｄ

Eｎd　Ｓｕｂ

Prｉｖatｅ　SubComｍand４＿Click（）

Pｉcturｅ1.Ｃlｓ

ＥndSub

SｔatｉcFunctionLoｇ１0（X）

　　Loｇ1０=Ｌog（X）/　Lｏg（10＃）　　　　　＇求以10为底的对数用的函数

EndFｕncｔion

PrｉｖateSub　Ｆoｒｍ_Ｌoad（）

　Ｏｐｔioｎ１.Vａluｅ=Tｒue

Optｉｏn4.Ｖalｕe=Tｒue

EｎｄSｕb

PrivateＳub　Optiｏn１_Click（）

ｊx=1

End　Ｓuｂ

PrivateＳubOption2_Ｃlick（）

jx＝2

EndSub

PrivａteSuｂOption３_Clｉck（）

　zy=2

Labｅl４５．Ｖｉｓｉｂｌe＝True

　Texｔ4５.Visiｂle　=Trｕe

ＥndSｕｂ

PrivatｅＳuｂ　Ｏｐtion4＿Ｃｌick（）

zｙ=1

　Labeｌ4５.Visｉbｌｅ　＝　Faｌｓｅ

　　Teｘｔ45.Visiblｅ　＝Faｌse

EndＳub

运行示例如下:

　１、当选择是汽油机且非增压时,如图：

运算结果如下：

此时的示功图如下:

2、当选择是柴油机且非增压时,如图：

运算结果如下:

此时的示功图如下：

１、当选择是柴油机且增压时，如图:

运算结果如下:

此时的示功图如下:

4.实验总结