打印的vb程序.docx

打印的vb程序.docx

《打印的vb程序.docx》由会员分享，可在线阅读，更多相关《打印的vb程序.docx（17页珍藏版）》请在冰豆网上搜索。

打印的vb程序

3、实验内容

3.2.1开始界面及代码

PrivateSubCommand4_Click（）

IfOption1=TrueThenForm2.Show

IfOption2=TrueThenForm3.Show

EndSub

PrivateSubForm_Load（）

EndSub

PrivateSubCommand3_Click（）

End

EndSub

3.2.2汽油机工作过程计算

DimD,S,i,Ne,n,ε,Vh,R,L,ge,P0,T0,gc,gh,go,mt,ι,α,ξz,Tr,φi,n1,Pz,γ,ΔT,ηm,n2,Vc

DimPr,Pa,Ta,ηv,Pc,Tc,λ,L0,M1,M0,M2,μ0,μ,Hμ,Hμ1,Tz,ρ,Pb,Tb,δ,nh,Nn,ΔHμ,Cv1,Cv2,A,B,C

PrivateSubCommand1_Click（）

k=1.4

D=Val（Text1.Text）

S=Val（Text2.Text）

Ne=Val（Text3.Text）

nh=Val（Text4.Text）

ε=Val（Text5.Text）

L=Val（Text6.Text）

R=Val（Text7.Text）

ge=Val（Text8.Text）

P0=Val（Text13.Text）

T0=Val（Text12.Text）

gc=Val（Text9.Text）

gh=Val（Text10.Text）

go=Val（Text11.Text）

mt=Val（Text14.Text）

i=Val（Text15.Text）

ι=Val（Text16.Text）

α=Val（Text38.Text）

ξz=Val（Text40.Text）

φi=Val（Text41.Text）

ηm=Val（Text42.Text）

n=Val（Text45.Text）

Vh=（S*0.01）*3.14*（D*0.01）^2/4

'排气过程

Pr=1+0.3*n/nh

Tr=850+350*n/nh

Text17.Text=Val（Pr）

Text18.Text=Val（Tr）

'进气过程

δ=0.5

ΔT=ΔTh*（110-0.0125*n）/（110-0.0125*nh）

Pa=Int（10^3*（P0*（1-（（10*n^2/520/10^6）/10.6^2/0.75^2）*（ε-δ）^2/（ε-1）^2）^（k/（k-1））））/10^3

γ=Int（10^3*（T0+ΔT）*Pr/Tr/（ε*Pa-Pr））/10^3

Text44.Text="0"&γ

Text19.Text="0"&Pa

Ta=（T0+ΔT+γ*Tr）/（1+γ）

ηv=ε*Pa*T0/（ε-1）/P0/Ta/（1+γ）

Text20.Text=Val（Int（1000*Ta）/1000）

Text21.Text="0"&Val（Int（1000*ηv）/1000）

'压缩过程

n1=1.46-0.05*Nn/n

Pc=Pa*ε^n1

Tc=Ta*ε^（n1-1）

Text22.Text=Val（Int（1000*Pc）/1000）

Text23.Text=Val（Int（1000*Tc）/1000）

'燃烧过程

L0=（gc/12+gh/4-go/32）/0.21

M1=α*L0+1/mt

Cv1=（4.815+0.415*10^（-3）*Tc）*4.1868

M2=0.79*α*L0+gh/2+gc/12

A=（3.7+3.3*α）*10^（-4）*4.1868

B=（4.4+0.62*α）*4.1868

μ0=M2/M1

μ=（μ0+γ）/（1+γ）

Hμ=44100

Hμ1=58000*（1-α）

C=（ξz*（Hμ-Hμ1）/M1/（1+γ）+Cv1*Tc）/μ

Tz=（-B+（B^2+4*A*C）^（1/2））/2/A

λ=μ*Tz/Tc

Pz=λ*Pc

Text24.Text=Val（Int（1000*λ）/1000）

Text25.Text=Val（Int（1000*Pz）/1000）

Text26.Text="0"&Val（Int（1000*L0）/1000）

Text27.Text=Val（Int（1000*μ）/1000）

Text28.Text=Val（Int（1000*Tz）/1000）

'膨胀过程

ρ=1

n2=1.14+0.035*nh/n

δ=Int（100*ε/ρ）/100

Pb=Int（1000*Pz/δ^n2）/1000

Tb=Int（Tz/δ^（n2-1））

Text30.Text=Val（δ）

Text31.Text=Val（Pb）

Text32.Text=Val（Tb）

'最后结果

Pi1=Pc*（λ*（ρ-1）+λ*ρ*（1-1/δ^（n2-1））/（n2-1）-（1-1/ε^（n1-1））/（n1-1））/（ε-1）

Pi=Int（1000*φi*Pi1）/1000

ηi=Int（100*8.314*α*L0*T0*Pi/Hμ/P0/ηv）/100

gi=Int（100*（3.6*10^6/Hμ/ηi））/100

Pm=0.4*（1+n/1000）

ηe=Int（100*ηm*ηi）/100

ge=Int（3.6*10^6/Hμ/ηe）

Pe=Int（1000*Pi*ηm）/1000

Ne=Int（1000*1*Vh*Pe*n/1200）/1000

Text33.Text=Val（Pi）

Text34.Text="0"&Val（ηi）

Text35.Text="0"&Val（ηe）

Text36.Text=Val（ge）

Text37.Text=Val（Ne）

Command2.Enabled=True

EndSub

PrivateSubCommand2_Click（）

Picture1.Cls

Picture1.Scale（-0.2,100）-（1,-8）

Picture1.ForeColor=vbBlack

Picture1.Line（0,0）-（0.9,0）'画X轴

Picture1.Line（0,0）-（0,90）'画Y轴

Picture1.CurrentX=0.8:

Picture1.CurrentY=-0.02:

Picture1.Print"V（L）"

Picture1.CurrentX=-0.02:

Picture1.CurrentY=95:

Picture1.Print"P（bar）"

ForY=0To85Step5'画刻度

Picture1.Line（0,Y）-（0.02,Y）

Picture1.CurrentX=-0.05:

Picture1.CurrentY=Y+1.5:

Picture1.PrintY

NextY

Picture1.Line（0,90）-（-0.01,88）

Picture1.Line（0,90）-（0.01,88）

ForX=0.1To0.8Step0.1'画刻度

Picture1.Line（X,0）-（X,0.2）

Picture1.CurrentX=X-0.05:

Picture1.CurrentY=-0.02:

Picture1.PrintX

NextX

Picture1.Line（0.9,0）-（0.89,1.5）

Picture1.Line（0.9,0）-（0.89,-1.5）

Vc=Vh/（ε-1）

Picture1.ForeColor=vbBlue

Picture1.CurrentX=Vc+0.01:

Picture1.CurrentY=90:

Picture1.Print"Vc"

Picture1.Line（Vc,90）-（Vc,0）'画Vc线

Picture1.ForeColor=vbRed

Picture1.Line（0,P0）-（0.8,P0）'画P0线

Picture1.CurrentX=0.81:

Picture1.CurrentY=4:

Picture1.Print"P0"

Va=Vh+Vh/（ε-1）

Vc=Vh/（ε-1）

Picture1.ForeColor=vbGreen

Forφx=0To3.14Step0.0001'画压缩曲线

Vcx=Vh/2*（（1-Cos（φx））+（1-Cos（2*φx））*R/4/L）+Vc

Pcx=Pa*（Va/Vcx）^n1

Picture1.PSet（Vcx,Pcx）

Nextφx

Forφx=3.14To6.28Step0.0001'画膨胀曲线

Vbx=Vh/2*（（1-Cos（φx））+（1-Cos（2*φx））*R/4/L）+Vc

pbx=Pb*（Va/Vbx）^n2

Ifpbx

Picture1.PSet（Vbx,pbx）

Else:

Picture1.Line（Vc,Pz）-（Vbx,Pz）

EndIf

Nextφx

Vb1=Vh/2*（（1-Cos（6.4））+（1-Cos（2*6.4））*R/4/L）+Vc

Picture1.Line（Va,Pa）-（Va,Pb）'画放热线

Picture1.Line（Vc,Pc）-（Vc,Pz）'画加热线

Command2.Enabled=False

EndSub

PrivateSubCommand3_Click（）

Picture1.Cls

EndSub

PrivateSubCommand4_Click（）

LoadForm1

Form1.Show

Form2.Hide

EndSub

PrivateSubForm_Load（）

Command2.Enabled=False

EndSub

3.2.3柴油机工作过程计算

DimD,S,i,Ne,n,ε,Vh,R,L,P0,T0,gc,gh,go,α,ξz,Tr,φi,n1,Pz,γ,ΔT,ηm,n2,Vc

DimPr,Pa,Ta,ηv,Pc,Tc,λ,L0,M1,M0,M2,μ0,μ,Hμ,Tz,ρ,Pb,Tb,δ,nh,Nn,ΔHμ,Cv1,Cv2,A,B,C

PrivateSubCommand1_Click（）

k=1.4

D=Val（Text1.Text）

S=Val（Text2.Text）

Ne=Val（Text3.Text）

nh=Val（Text4.Text）

ε=Val（Text5.Text）

L=Val（Text6.Text）

R=Val（Text7.Text）

ge=Val（Text8.Text）

P0=Val（Text13.Text）

T0=Val（Text12.Text）

gc=Val（Text9.Text）

gh=Val（Text10.Text）

go=Val（Text11.Text）

i=Val（Text15.Text）

ι=Val（Text16.Text）

α=Val（Text38.Text）

ξz=Val（Text40.Text）

φi=Val（Text41.Text）

ηm=Val（Text42.Text）

n=Val（Text45.Text）

Vh=（S*0.01）*3.14*（D*0.01）^2/4

'排气过程

Pr=1.08*P0

Tr=850:

mt=0:

ΔTh=20

Text17.Text=Val（Pr）

Text18.Text=Val（Tr）

'进气过程

δ=0.5

ΔT=ΔTh*（110-0.0125*n）/（110-0.0125*nh）

Pa=Int（10^3*（P0*（1-（（10*n^2/520/10^6）/（n^2/520*10^6）^2/0.75^2）*（ε-δ）^2/（ε-1）^2）^（k/（k-1））））/10^3

γ=Int（10^3*（T0+ΔT）*Pr/Tr/（ε*Pa-Pr））/10^3

Text44.Text="0"&γ

Text19.Text=Pa

Ta=（T0+ΔT+γ*Tr）/（1+γ）

ηv=ε*Pa*T0/（ε-1）/P0/Ta/（1+γ）

Text20.Text=Val（Int（1000*Ta）/1000）

Text21.Text="0"&Val（Int（1000*ηv）/1000）

'压缩过程

n1=1.46-0.05*Nn/n

Pc=Pa*ε^n1

Tc=Ta*ε^（n1-1）

Text22.Text=Val（Int（1000*Pc）/1000）

Text23.Text=Val（Int（1000*Tc）/1000）

'燃烧过程

L0=（gc/12+gh/4+gs/32-go/32）/0.21

M1=α*L0

λ=Val（Text14.Text）

Cv1=（4.815+0.415*10^（-3）*Tc）*4.1868

M2=α*L0+gh/4+go/32

A=（3.7+3.3/α）*10^（-4）*4.1868

B=（4.8+2.2/α）*4.1868

μ0=M2/M1

μ=（μ0+γ）/（1+γ）

Hμ=42500

C=（ξz*Hμ/M1/（1+γ）+（Cv1+1.985*λ）*Tc）/μ

Tz=（-B+（B^2+4*A*C）^（1/2））/2/A

ρ=μ*Tz/λ/Tc

Pz=λ*Pc

Text24.Text=Val（Int（1000*λ）/1000）

Text25.Text=Val（Int（1000*Pz）/1000）

Text26.Text="0"&Val（Int（1000*L0）/1000）

Text27.Text=Val（Int（1000*μ）/1000）

Text28.Text=Val（Int（1000*Tz）/1000）

'膨胀过程

n2=1.14+0.035*Nn/n

δ=Int（100*ε/ρ）/100

Pb=Int（1000*Pz/δ^n2）/1000

Tb=Int（Tz/δ^（n2-1））

Text30.Text=Val（δ）

Text31.Text=Val（Pb）

Text32.Text=Val（Tb）

'最后结果

Pi1=Pc*（λ*（ρ-1）+λ*ρ*（1-1/δ^（n2-1））/（n2-1）-（1-1/ε^（n1-1））/（n1-1））/（ε-1）

Pi=Int（1000*φi*Pi1）/1000

ηi=Int（100*8.314*α*L0*T0*Pi/Hμ/P0/ηv）/100

gi=Int（100*（3.6*10^6/Hμ/ηi））/100

Pm=0.04*（1+n/1000）

ηe=Int（100*ηm*ηi）/100

ge=Int（3.6*10^6/Hμ/ηe）

Pe=Int（1000*Pi*ηm）/1000

Ne=Int（1000*1*Vh*Pe*n/1200）/1000

Text33.Text=Val（Pi）

Text34.Text="0"&Val（ηi）

Text35.Text="0"&Val（ηe）

Text36.Text=Val（ge）

Text37.Text=Val（Ne）

Command2.Enabled=True

EndSub

PrivateSubCommand2_Click（）

Picture1.Cls

Picture1.Scale（-0.2,120）-（1.1,-10）

Picture1.ForeColor=vbBlack

Picture1.Line（0,0）-（1,0）'画X轴

Picture1.Line（0,0）-（0,115）'画Y轴

Picture1.CurrentX=0.95:

Picture1.CurrentY=-0.02:

Picture1.Print"V（L）"

Picture1.CurrentX=0.02:

Picture1.CurrentY=115:

Picture1.Print"P（bar）"

ForY=0To110Step5'画刻度

Picture1.Line（0,Y）-（0.02,Y）

Picture1.CurrentX=-0.08:

Picture1.CurrentY=Y+1.5:

Picture1.PrintY

NextY

Picture1.Line（0,115）-（-0.01,113.5）'画箭头

Picture1.Line（0,115）-（0.01,113.5）

ForX=0.1To1Step0.1'画刻度

Picture1.Line（X,0）-（X,0.2）

Picture1.CurrentX=X-0.05:

Picture1.CurrentY=-0.02:

Picture1.PrintX

NextX

Picture1.Line（1,0）-（0.99,2）'画箭头

Picture1.Line（1,0）-（0.99,-2）

Vc=Vh/（ε-1）

Picture1.ForeColor=vbBlue

Picture1.CurrentX=Vc+0.01:

Picture1.CurrentY=10:

Picture1.Print"Vc"

Picture1.Line（Vc,120）-（Vc,0）'画Vc线

Picture1.ForeColor=vbRed

Picture1.Line（0,P0）-（0.95,P0）'画P0线

Picture1.CurrentX=0.95:

Picture1.CurrentY=4:

Picture1.Print"P0"

Va=Vh+Vh/（ε-1）

Vc=Vh/（ε-1）

Picture1.ForeColor=vbGreen

Forφx=0To3.14Step0.0001'画压缩曲线

Vcx=Vh/2*（（1-Cos（φx））+（1-Cos（2*φx））*R/4/L）+Vc

Pcx=Pa*（Va/Vcx）^n1

Picture1.PSet（Vcx,Pcx）

Nextφx

Forφx=3.14To6.28Step0.0001'画膨胀曲线

Vbx=Vh/2*（（1-Cos（φx））+（1-Cos（2*φx））*R/4/L）+Vc

pbx=Pb*（Va/Vbx）^n2

Ifpbx

Picture1.PSet（Vbx,pbx）

Else:

Picture1.Line（Vc,Pz）-（Vbx,Pz）

EndIf

Ne