#195柴油机工作过程计算

图1,计算结果;

图2,P-V图

80

70

60

50

40

30

20

0.10.20.30.40.5C.60.70.80.9

图3,P中图

示功图

图4,P-6图

示功图

图5,h-①图

示功图

图6,活塞温度图

示功图

程序代码：

OptionExplicit

Dimn,nn,pr,poAsSingle

Dimtr,pa,s,d,vh,kkkAsSingle

Dimpab,e,paa,pad,pacAsSingle

Dimtt,T,yy,fa,X,Y,taAsSingle

Dimnv,nl,tc,jx,pc,fxAsSingle

Dim1,vc,Va,gc,gh,gOAsSingle

Dimvex,r,lo,A,nit,nil,m2,uO,uAsSingle

Dimkhu,cv1,cv2,tz,ccc,hu,qq,aaa,bbbAsSingle

Dim18,n2,pz,H,tb,pil,pi,(px,tp,pAsSingle

Dimfi,pm,nm,ni,gi,ne,pe,ge,nne,i,v,pb,pcx,vbx,pbxAsDouble

Dimaa(50),ylsgbAsSingle

PrivateSubCommand4_Click()

Picture1.Cis

Picture1.ForeColor=vbGreen

PictureI.DrawWidth=1.5

PictureI.Scale(-0.1,99)-(1,-10)

Picturel.Line(0,0)-(0,95)'y轴

Picturel.Line(0,0)-(0.95,0)'x轴

Picturel.Line(0.93,l)-(0.95,0)

Picturel.Line(0.93,-l)-(0.95,0)

Picturel.Line(0.015.94)-(0.95)

Picturel.Line(-0.015,94)-(0,95)

Fori=1To9Step1

Picturel.Line(0,i*10)-(0.015,i*10)

Picturel.Line(0.1*i,0)-(0.1*i,-1)

Picturel.CurrentX=-0.1:Picture!.CurrenlY=10*i+3:Picture1.Printi*10

Picturei.CuirentX=0.1*i-0.04:Picture1.CurrenlY=-3:Picture1.Print"0"&i*0.1

Nexti

Picturel.CurrentX=0:Picture1.CurrenlY=95:Picture1.Print"P(atm)"

Picturel.CurrentX=0.8:Picture!.CurrentY=10:Picture1.Print"V(L)"

'压缩过程曲线

Ifvh=0Then

MsgBox”请先输入参数，并计算”,0,“提示”:GoToq

EndIf

vc=vh/e

DimXAsSingle:Dimpcx,vexAsSingle

Picturei.DrawWidth=1.3

Picturel.Line(vc,10*po)-(0.9,10*po)

PictureI.Line(vc,po)-(vc,88)

Picturei.DrawWidth=1.5

Picture1.ForeColor=vbYellow

ForX=vTovcStep-0.0001

vex=X

pcx=pa*(v/vex)Anl

Picturel.PSet(vex,pcx)

NextX

'燃烧过程曲线

Picture1.Line(vex,pcx)-(vc,pz)

'膨胀过程曲线

Dimvbx,pbxAsSingle

Picture).Line(vc.pz)-(vc*p.pz):vc=vc*p

ForX=vcTovStep0.0001

vbx=X

pbx=pz*((vc/vbx)An2)

Picturel.PSet(vbx,pbx)

NextX

'排气过程曲线

Picture1.Line(vbx,pbx)-(v,pa)

Picture).DrawWidth=1

Picture1.ForeColor=vbRed

Picturel.CurrentX=0.1:PictureI.CurrentY=8:Picture1.Print"po"

Picturel.CurrentX=0.02:Picturel.CurrentY=20:Picture1.Print"vc"

q-

EndSub

PrivateSubCommandl_Click()

s=Val(Texl2.Texl):d=Val(Text4,Text)

n=Val(Textl.Text):po=Val(Textl5.Text)/10:

e=VaI(Text5.Text):tt=Vai(Text17.Text):fa=Val(Text18.Text):r=Val(Text2.Text)

1=Val(Text3l.Text)/3.14/((d/2)A2)

gc=Val(Textll.Text):gh=Val(Textl2.Text):gO=Val(Text8.Text)

A=Val(Textl9.Texl):mt=Val(Text20.Text):fi=Val(Text9.Text):i=Val(Text3.Text)

nl=Val(Textl4.Texl):n2=Val(Text33.Text)

'排气过程

pr=Int(1000*10.8*po)/1000:tr=350/(1.2/Log(n)*Log(10)+0.005*(e-3)+0.01*(fa

-1)):hu=42500'柴油机暂时先赋的值

'进气过程

paa=nA2/520000000

pac=((e-0.5)A2)/((e-1)A2)'残余废气收缩系数暂时取的0.5

T=20't进气温升取的20

pa=0.9*po*10，进气压力

yy=(tt+T)*pr/lr/(e*pa-pr),残余废气系数yy

ta=(tt+T+yy*tr)/(l+yy)'进气温度

nv=it*pa*e*(e*pa-pr)/(tt+T)/(10*po)/(e-I)/e/pa充气系数

'压缩过程

'平均压缩多变指数

pc=Int(100*(pa*eA(nl)))/100，压缩终点压力

tc=Int(10*ta*eA(nl-1))/10,压缩终点温度

'燃烧过程

10=1/0.21*(gc/12+gh/4-g0/32)

ml=A*io:m2=A*lo+gh/4+g0/32

uO=m2/m1

u=(uO+yy)/(l+yy)'实际分子变更系数

'化学损失

khu=O

cvl=(4.815+0.000415*tc)*4.1868

qq=Val(Text30.Text)'---------------------------qq热量利用系数选取

ccc=(qq*(hu-khu)/ml/(1+yy)+cvl*tc)/u

aaa=(3.3/A+3.7)*4.1868/10000，柴油机

bbb=(4.8+2.2/A)*4.1868

tz=Int((-bbb+(bbbA2+4*aaa*ccc)A0.5)/2/aaa)燃烧终点温度

ylsgb=Val(Text29.Text)

p=u*tz/ylsgb/tc1--柴油机预期膨胀比--

pz=Int(100*(ylsgb*pc))/100'最高燃烧压力

pb=Int(l00*(pz/(eAn2)))/100'膨胀终点压力

lb=Int(tz/(eA(n2-1)))'膨胀终点温度

pil=pc/(e-l)*(ylsgb*p/(n2-1)*(1-1/eA(n2-1))-1/(nl-1)*(1-1/eA(nl-1)))

pi=fi*pil'平均指示压力

pm=0.1*(Val(Text32.Text))A0.5*(1+n/1000)，机械损失压力

nm=1-pm/pi'机械效率

ni=Int(10000*8.314*ml*u*pi/hu/10/po/nv)/10000'指示热效率

gi=3600000/hu/ni'指示燃油消耗率

ne=nm*ni，有效热效率

pe=pi*nm，平均有效压力

ge=3600000/hu/ne'有效燃油消耗率

vh=Int(10000*3.14*(d/2)A2*r/10A6)/10000

v=Int(l0000*(vh/e+vh))/KXXX)

nne=i*v*pe*n/l224

Texl7⑴.Texl="0"&Inl(10000*ni)/10000

Text7(3).Text=Int(100*gi)/100

Text7(4).Text="0"&Int(10000*ne)/10000

Text7(5).Text=Int(1000*pe)/1000

Text7(6).Text=Int(100*ge)/100

Text7(7).Text=Int(100*nne)/100

TextlO.Text=pr:TextB.Text=Int(10*tr)/10:Text21.Text="0"&Int(100*pa)/100:

Text22.Text=Int(ta*10)/10

Text23.Text=pz:Text24.Text=tz:Text25.Text=tb:Text26.Text=pb

Text27.Text=pc:Text28.Text=tc

EndSub

PrivateSubCommand5_CIick()

Picture1.Cis

Picture1.Scale(-10,120)-(780.-10)

Picture1.Line(0,0)-(725,0)

Picture1.Line(0,0)-(0,115)

vc=vh/(e-I)

Picture1.Line(0,po)-(720,po),RGB(0,220,0)

Picturel.CurrentX=680:Picturel.CurrentY=8:Picture1.Print"(p(°AC)"

Picturel.CurrentX=10:Picturel.CurrentY=129:PictureI.Print"P(bar)"

ForY=10To110Step10#

Picture1.Line(0,Y)-(7,Y)

Picturel.CurrentX=-0.1:Picturel.CurrentY=Y:Picturel.PrintY

NextY

Picture!.Line(0,115)-(-5,112)

Picture!.Line(0,115)-(5,112)

ForX=0To720Step60

Picture!.Line(X,0)-(X,2)

Picturel.CurrentX=X-20:Picturel.CurrentY=-2:Picturel.PrintX

NextX

Picture!.Line(725,0)-(716,2)

Picturel.Line(725,0)-(716,-2)

Ifvh=0Then

MsgBox”请先输入参数，并计算”,0,“提示”:GoToq

EndIf

Va=vh+vh/(e-I)

vc=vh/(e-I)

For(px=3.14To6.28Step0.0001

vex=vh/2*((I-Cos(<l>x))+(l-Cos(2*4>x))*r*I/4)+vc'压缩过程

pcx=pa*(Va/vex)An1

(p=(px*180/3.14

Picturel.PSet(<p,pcx),RGB(0,400,200)

Nexttpx

Va=vh+vh/(e-1)

vc=vh/(e-I)

For(px=9.42To6.28Step-0.0001

vbx=vh/2*((1-Cos((px))+(1-Cos(2*(px))*r*1/4)+vc

pbx=pb*(Va/vbx)An2

(p=(px*180/3.14

Ifpbx<=pzThen

Picturel.PSet(<p+20,pbx),RGB(400,0,220)

Else:Picture1.Line(360,pz)-(6,pz),RGB(400,100,0)，做功过程

Picture1.Line(360.pc)-(360.pz).RGB(0.0.220)

EndIf

Next(px

Picturel.Line(540,2.11)-(720,2.11)

Picture1.Line(360,54.202)-(360,86.72),RGB(IOO,0,220)

Picturel.Line(360,86.72)-(380,86.72),RGB(0,100,220)

q：

EndSub

PrivateSubCommand6_Click()

DimT,ex,A,C,bx,ea,b

Picture1.Cis

PictureI.Scale(-60,2600)-(780,-150)

Picturel.Line(0,0)-(740,0)

PictureI.Line(0,0)-(0,2500)

vc=vh/(e-1)

PictureI.CurrentX=680:Picturel.CurrentY=120:Picture1.Print"(p(°AC)"

PictureI.CurrentX=10:Picturel.CurrentY=2500:Picture1.Print"T(K)H

ForY=200To2200Siep200#

PictureI.Line(0,Y)-(7,Y)

Picture!.CurrentX=-60:Picturel.CurrentY=Y+30:Picture1.PrintY

NextY

PictureI.Line(0,25OO)-(-5,2400)

PictureI.Line(0,2500)-(5,2400)

ForX=0To720Step60

Picture1.Line(X,0)-(X,40)

Picturel.CurrentX=X-20:Picture1.CurrentY=-30:Picture1.PrintX

NextX

Picture1.Line(740,0)-(733,25)

Picture1.Line(740,O)-(733,-25)

Ifvh=0Then

MsgBox”请先输入参数.并计算”,0,“提示”:GoToq

EndIf

b=Val(Text6.Text)

For(px=3.14*2To3.11*2.35Step0.0001

A=i-Cos(ipx)+0.5*0.25*Sin((px)A2

C=1/(e-I)

ea=e

ex=1+A/(2*C)

bx=b-((352-(px*180/3.14)/(352-400))A3*(b-1)

T=(ta*ex*10)/(pa*ea*bx)

(p=(px*180/3.14

Picturel.PSet((p-120,T),RGB(220,0,0)

Next(px

For<px=3.14To2*3.14Step0.0001

A=I-Cos((px)+0.5*0.25*Sin((px)A2

C=1/(e-I)

ea=e

ex=1+A/(2*C)

bx=b-((352-<px*180/3.14)/(352-400))A3*(b-1)

T=(ta*ex*10)/(pa*ea*bx)

(p=(px*180/3.14

PictureI.PSet((p+120,T),RGB(0,0,220)

Next(px

q：

EndSub

PrivateSubCommand7_Click()

DimT,ex,A,Cm,bx,ea,b,Cu,C1,C2,C,hl,h2,H3

Picture1.Cis

Picture1.Scale(-60,2000)-(780,-150)

Picture1.Line(0,0)-(740,0)

Picture!.Line(0.0)-(0.1900)

vc=vh/(e-1)

PictureI.CurreniX=680:Picturel.CurrentY=120:Piciurel.Print"(p(°AC)"

PictureI.CurreniX=10:Picturel.CurrentY=1900:Piciurel.Print"h：W/(m2.K)"

ForY=200To1800Step200#

Piciurel.Line(0,Y)-(7,Y)

Picturel.CurrentX=-60:Picturel.CurrentY=Y+30:Piciurel.PrintY

NextY

Picture!.Line(0,1900)-(-5,1850)

PictureI.Line(0,1900)-(5,1850)

ForX=0To720Step60

PictureI.Line(X,0)-(X,40)

Picturel.CurrentX=X-20:Picturel.CurrentY=-30:Piciurel.PrintX

NextX

Picture!.Line(740,0)-(733,25)

Picture1.Line(740,O)-(733,-25)

Ifvh=0Then

MsgBox”请先输入参数，并计算”,0,“提示”:GoToq

EndIf

C=Val(Text34.Text)

Cm=s*2/3():Cl=2.28+().308*C:C2=3.2*10A(-3)

d=Val(Tex(4.Text)

A=vh*ta/(pa*0.03719)

For(px=3.14To6.28Step0.0001

Va=vh+vh/(e-I)

vc=vh/(e-I)

vex=vh/2*((1-Cos((px))+(1-Cos(2*(px))*r*1/4)+vc

pcx=pa*(Va/vex)An1

Va=vh+vh/(e-1)

vc=vh/(e-1)

A=1-Cos((px)+0.5*0.25*Sin((px)A2

C=1/(e-I)

ea=e

ex=1+A/(2*C)

bx=1.04-((352-(px*180/3.14)/(352-400))A3*(1.04-1)

T=(ta*ex*10)/(pa*ea*bx)

AAA

H=130*(d*0.00I)(-0.2)*(p*IO)(O.8)*T(-O.53)*(C1*Cm+C2*A*(p-1))人0.8

o=(px*180/3.14

Picturel.PSetC(p,H),RGB(220,0,0)

H=H+456.9

Next(px

hl=H/3.14

For(px=6.28To450*3.14/180Step0.0001

Va=vh+vh/(e-1)

vc=vh/(e-I)

vex=vh/2*((1-Cos((px))+(1-Cos(2*(px))*r*1/4)+vc

pcx=pa*(Va/vex)An1

Va=vh+vh/(e-1)

vc=vh/(e-1)

A=I-Cos((px)+0.5*0.25*Sin((px)A2

C=1/(e-1)

ea=e

ex=l+A/(2*C)

bx=1.04-((352-(px*180/3.14)/(352-400))A3*(1.04-1)

T=(ta*ex*10)/(pa*ea*bx)

H=130*(d*0.001)A(-0.2)*(p*10)A(0.8)*TA(-0.53)*(Cl*Cm+C2*A*(p-1))A0.8

o=(px*180/3.14

Picturel.PSet(<p,H),RGB(0,0,220)

H=H+1302.56

Next(px

h2=H/1.07

H3=(hl+h2)/2

Text35.Text=Int(H3)

q：

EndSub

PrivateSubCommand8_Click()

Dimtmax,C1,C2,fl,tl,tmaxl,ti,fx,fy,dyi,yi,xi

Picture1.Cis

Picture1.Scale(-10,130)-(130.-10)

Picture1.Line(0,100)-(110,100):Picture).Line(52.5,120)-(52.5,20)

Picturel.CurrentX=110:Picture!.CurrentY=105:Picturel.Print"+X"

Picturel.CurrentX=0:Picturel.CurrentY=105:Picturel.PrintH-X"

Picturel.CurrentX=55:Picturel.CurrentY=120:Picturel.Print""Y"

PictureI.CurrentX=55:Picturel.CurrentY=20:Picturel.Print"+Y"

Picturel.Line(llO,100)-(108,98):Picture1.Line(110,100)-(108,102)

Picturel.Line(52.5,20)-(50.5,22):Picturel.Line(52.5,20)-(54522)

'Picturel.Line(0,0)-(110,0)

'Picturel.Line(0,0)-(0,HO)

'vc=vh/(e-1)

,Picturel.CurrentX=680:Picturel.CurrentY=120:Picturel.Print"(p(°AC)"

,Picturel.CurrentX=1():Picturel.CurrentY=1900:Picturel.Print"h：W/(m2.K)"

''ForY=OTo105Step5