热力学计算VB程序代码

1、附录二：VB程序代码Dim Pa , Va , Pb , n1 , n2 As Single/“计算”命令按钮的单击事件/Private Sub Command1_Click()Dim P0, T0, Ne, n, e, e0, AF, Vr, Pz, Ez, Eb, Q, Tb, C, H, O, Hu, Ta, nv, Tc, Pc, L0, L, B , B0, Xz, Bz, Cvpmz, Cvpmb, Cppmz ,Hu_, i, Tz, p, g, Bb, Pi_, Pi , gi , ni , nm , Pe , ge , ne0 , S , D , S_D , Cm As Si

2、ngle Pa = 0 Va = 0 Pb = 0 n1 = 0 n2 = 0 P0 = Val(Text1(2).Text) T0 = Val(Text1(1).Text) Ne = Val(Text1(3).Text) n = Val(Text1(9).Text) e = Val(Text1(5).Text) e0 = Val(Text1(6).Text) AF = Val(Text1(4).Text) Vr = Val(Text1(13).Text) Pz = Val(Text1(10).Text) Ez = Val(Text1(0).Text) Eb = Val(Text1(7).Te

3、xt) Q = Val(Text1(14).Text) Tb = Val(Text1(8).Text) C = Val(Text1(15).Text) / 100 H = Val(Text1(16).Text) / 100 O = Val(Text1(12).Text) / 100 Hu = Val(Text1(11).Text) /检查输入参数/ Dim k As Integer For k = 0 To 16 If Len(Text1(k).Text) <= 0 Then a = MsgBox("请先输入全部参数", 0 + 0 + 48, "请确定&q

4、uot;) Exit Sub ElseIf Val(Text1(k).Text) <= 0 Then a = MsgBox("参数中有非法数据", 0 + 0 + 48, "请确定") Exit Sub End If Next kCall cqgccs(P0, T0, Tb, Vr, e, e0, Ta, Pa, nv) 调用充气参数计算过程Call ysgccs(n1, Ta, e, Tc, Pc, Pa) 调用压缩参数计算过程Call rsgccs(L0, C, H, O, AF, L, B0, B, Vr, Xz, Ez, Eb, Bz, B

5、b, Cvpmz, Cppmz, Cvpmb, Hu_, Hu, i, Pz, Pc, Tz, Tb, Tc) 调用燃烧参数计算过程Call pzgccs(p, B, Tz, i, Tc, g, e, n2, Bz, Bb, Tb, Hu_, Eb, Ez, L, Cvpmb, Cvpmz, Pb, Pz) 调用膨胀参数计算程Call zscs(Pi_, Pc, e, i, p, n2, n1, Pi, Pb, Pa, e0, gi, nv, P0, T0, ni, Hu, g, L) 调用指示参数计算过程Call yxcs(nm, Pe, Pi, ge, gi, ne0, ni) 调用有效参

6、数计算过程Call qgjbcc(S, D, S_D, Cm, n, Ne, Pe, Va) 调用汽缸基本尺寸参数计算过程/输出计算结果/Text2(0).Text = TaText2(1).Text = PaText2(2).Text = nvText2(3).Text = n1Text2(4).Text = TcText2(5).Text = PcText2(6).Text = L0Text2(7).Text = LText2(8).Text = B0Text2(9).Text = BText2(10).Text = XzText2(11).Text = BzText2(12).Text

7、= CvpmzText2(13).Text = CvpmbText2(14).Text = CppmzText2(15).Text = Hu_Text2(16).Text = iText2(17).Text = Tz Text2(18).Text = pText2(19).Text = gText2(20).Text = n2Text2(21).Text = PbText2(22).Text = Pi_Text2(23).Text = PiText2(24).Text = giText2(25).Text = niText2(26).Text = nmText2(27).Text = PeTe

8、xt2(28).Text = geText2(29).Text = ne0Text2(30).Text = SText2(31).Text = DText2(32).Text = VaText2(33).Text = S_DText2(34).Text = CmEnd Sub/“绘图”命令按钮的单击事件/Private Sub Command2_Click() Dim x0 As Single, y0 As Single, x As Single If (n1 = 0) Then a = MsgBox("请先计算", 0 + 0 + 48, "请确定")

9、 Exit Sub End If Call Load(Form2) Form2.Show x0 = 1920 y0 = 3480 For x = Va To e0 * Va Step 0.01 PSet (x, y0 - Pa * (x - x0) n1) Next x For x = Va To g * Va Step 0.01 PSet (x, y0 - Pa * (x - x0) n2 Next xEnd Sub/定义充气参数计算过程/Sub cqgccs(P0 As Single, T0 As Single, Tb As Single, Vr As Single, e As Singl

10、e, e0 As Single, Ta As Single, Pa As Single, nv As Single) Ta = (T0 + 5 + 1.11 * Vr * Tb) / (1 + Vr) Pa = P0 nv = e * T0 / (e0 - 1) * Ta * (1 + Vr)End Sub /定义压缩参数计算过程/Sub ysgccs(n1 As Single, Ta As Single, e As Single, Tc As Single, Pc As Single, Pa As Single) Dim n10, n11 As Single n10 = 1 Do n11 =

11、 n10 n10 = 8.315 / (19.26 + 0.0025 * Ta * (1 + e (n11 - 1) + 1 Loop While Abs(n11 - n10) > 0.001 n1 = n10Tc = Ta * e (n1 - 1)Pc = Pa * e n1End Sub/定义燃烧参数计算过程/Sub rsgccs(L0 As Single, C As Single, H As Single, O As Single, AF As Single, L As Single, B0 As Single, B As Single, Vr As Single, Xz As S

12、ingle, Ez As Single, Eb As Single, Bz As Single, Bb As Single, Cvpmz As Single, Cppmz As Single, Cvpmb As Single, Hu_ As Single, Hu As Single, i As Single, Pz As Single, Pc As Single, Tz As Single, Tb As Single, Tc As Single) Dim a, Cv_, Cv_, Tz0, Tz1, Cv_0, Cv_1, Cv_0, Cv_1, Cvpmz0, Cvpmz1, Cppmz0,

13、 Cppmz1 As Single a = AF / 14.2 Cv_ = 19.26 + 0.0025 * Tb Cv_ = 20.47 + 0.0036 * Tb L0 = 1 / 0.21 * (C / 12 + H / 4 - O / 32) L = a * L0 B0 = 1 + 0.0639 / a B = (B0 + Vr) / (1 + Vr) Xz = Ez / Eb Bz = 1 + (B0 - 1) / (1 + Vr) * Xz Bb = 1 + (B0 - 1) / (1 + Vr) Cvpmb = (1.064 + Vr) * Cv_ + (a - 1) * (1

14、+ Vr) * Cv_) / (a * (1 + Vr) + 0.064) Hu_ = Hu + L * (1 + Vr) * (B * Cvpmb - 19.26 + 0.0025 * 293) * 293 i = Pz / Pc Cv_0 = 20 Cv_0 = 20 Cvpmz0 = 30 Cppmz0 = 30 Tz0 = 1980Do Cv_1 = Cv_0 Cv_1 = Cv_0 Cvpmz1 = Cvpmz0 Cppmz1 = Cppmz0 Tz1 = Tz0 Cv_0 = 19.26 + 0.0025 * Tz1 Cv_0 = 20.47 + 0.0036 * Tz1 Cvpm

15、z0 = (1.064 * Xz + Vr) * Cv_1 + (a * (1 + Vr) - (Xz + Vz) * Cv_1) / (a * (1 + Vr) + 0.064 * Xz) Cppmz0 = Cvpmz1 + 8.315 Tz0 = (Ez * Hu_ / (a * L0) + (Cv_ + 8.315 * i) * Tc + Vr * (Cvpmb + 8.315 * i) * Tz1) / (Bz * (1 + Vr) * Cppmz1) Loop While Abs(Cv_1 - Cv_0) > 0.001 Or Abs(Cv_1 - Cv_0) > 0.0

16、01 Or Abs(Cppmz1 - Cppmz0) > 0.001 Or Abs(Cvpmz1 - Cvpmz0) > 0.001 Or Abs(Tz1 - Tz0) > 0.001 Cv_ = Cv_0 Cv_ = Cv_0 Cvpmz = Cvpmz0 Cppmz = Cppmz0 Tz = Tz0End Sub /定义膨胀参数计算过程/Sub pzgccs(p As Single, B As Single, Tz As Single, i As Single, Tc As Single, g As Single, e As Single, n2 As Single,

17、Bz As Single, Bb As Single, Tb As Single, Hu_ As Single, Eb As Single, Ez As Single, L As Single, Cvpmb As Single, Cvpmz As Single, Pb As Single, Pz As Single)p = 1 g = e / p n2 = 8.315 * (Bz * Tz / Bb - Tb) / (Hu_ * (Eb - Ez) / (L * (1 + Vr) * B) + Bz / Bb * Cvpmz * Tz - Cvpmb * Tb) + 1 Pb = Pz / g

18、 n2End Sub/定义指示参数计算过程/Sub zscs(Pi_ As Single, Pc As Single, e As Single, i As Single, p As Single, n2 As Single, n1 As Single, Pi As Single, Pb As Single, Pa As Single, e0 As Single, gi As Single, nv As Single, P0 As Single, T0 As Single, ni As Single, Hu As Single, g As Single, L As Single)Pi_ = Pc

19、 / (e - 1) * (i * p * (p - 1) + i * p / (n2 - 1) * (1 - 1 / g (n2 - 1) - 1 / (n1 - 1) * (1 - 1 / e (n1 - 1)Pi = (Pi_ * (e - 1) + (Pb - Pa) * (e0 - e) / (e0 - 1) * 0.92 gi = nv * P0 / (28.96 * L * T0 * Pi) * 12.54 * 10 6 ni = 3600 * 10 3 / (gi * Hu)End Sub/定义有效参数计算过程/Sub yxcs(nm As Single, Pe As Single, Pi As Single, ge As Single, gi As Single, ne0 As Single, ni As Single)nm = 0.8 Pe = Pi * nm ge = gi / nm ne0 = ni * nmEnd Sub/定义汽缸基本尺寸参数计算过程/Sub qgj