柴油机工作过程计算

1、 发动机热负荷与热强度实验报告Option Explicit Dim n, nn, pr, po As Single Dim tr, pa, s, d, vh, kkk As Single Dim pab, e, paa, pad, pac As Single Dim tt, T, yy, fa, X, Y, ta As Single Dim nv, n1, tc, jx, pc, fx As Single Dim l, vc, Va, gc, gh, g0 As Single Dim vcx, r, lo, A, mt, m1, m2, u0, u As Single Dim khu, cv

2、1, cv2, tz, ccc, hu, qq, aaa, bbb As Single Dim l8, n2, pz, H, tb, pi1, pi, x, , p As Single Dim fi, pm, nm, ni, gi, ne, pe, ge, nne, i, v, pb, pcx, vbx, pbx As Double Dim aa(50), ylsgb As SinglePrivate Sub Command4_Click() Picture1.Cls Picture1.ForeColor = vbGreen Picture1.DrawWidth = 1.5 Picture1.

3、Scale (-0.1, 99)-(1, -10) Picture1.Line (0, 0)-(0, 95) 'y轴 Picture1.Line (0, 0)-(0.95, 0) 'x轴 Picture1.Line (0.93, 1)-(0.95, 0) Picture1.Line (0.93, -1)-(0.95, 0) Picture1.Line (0.015, 94)-(0, 95) Picture1.Line (-0.015, 94)-(0, 95) For i = 1 To 9 Step 1 Picture1.Line (0, i * 10)-(0.015, i *

4、10) Picture1.Line (0.1 * i, 0)-(0.1 * i, -1) Picture1.CurrentX = -0.1: Picture1.CurrentY = 10 * i + 3: Picture1.Print i * 10 Picture1.CurrentX = 0.1 * i - 0.04: Picture1.CurrentY = -3: Picture1.Print "0" & i * 0.1 Next i Picture1.CurrentX = 0: Picture1.CurrentY = 95: Picture1.Print &qu

5、ot; P （atm）" Picture1.CurrentX = 0.8: Picture1.CurrentY = 10: Picture1.Print "V （L）" '压缩过程曲线 If vh = 0 Then MsgBox "请先输入参数,并计算", 0, "提示": GoTo q End If vc = vh / e Dim X As Single: Dim pcx, vcx As Single Picture1.DrawWidth = 1.3 Picture1.Line (vc, 10 * po)-(0.9

6、, 10 * po) Picture1.Line (vc, po)-(vc, 88) Picture1.DrawWidth = 1.5 Picture1.ForeColor = vbYellow For X = v To vc Step -0.0001 vcx = X pcx = pa * (v / vcx) n1 Picture1.PSet (vcx, pcx) Next X '燃烧过程曲线 Picture1.Line (vcx, pcx)-(vc, pz) '膨胀过程曲线 Dim vbx, pbx As Single Picture1.Line (vc, pz)-(vc *

7、 p, pz): vc = vc * p For X = vc To v Step 0.0001 vbx = X pbx = pz * (vc / vbx) n2) Picture1.PSet (vbx, pbx) Next X '排气过程曲线 Picture1.Line (vbx, pbx)-(v, pa) Picture1.DrawWidth = 1 Picture1.ForeColor = vbRed Picture1.CurrentX = 0.1: Picture1.CurrentY = 8: Picture1.Print "po" Picture1.Cur

8、rentX = 0.02: Picture1.CurrentY = 20: Picture1.Print "vc"q:End SubPrivate Sub Command1_Click() s = Val(Text2.Text): d = Val(Text4.Text) n = Val(Text1.Text): po = Val(Text15.Text) / 10: e = Val(Text5.Text): tt = Val(Text17.Text): fa = Val(Text18.Text): r = Val(Text2.Text) l = Val(Text31.Tex

9、t) / 3.14 / (d / 2) 2) gc = Val(Text11.Text): gh = Val(Text12.Text): g0 = Val(Text8.Text) A = Val(Text19.Text): mt = Val(Text20.Text): fi = Val(Text9.Text): i = Val(Text3.Text) n1 = Val(Text14.Text): n2 = Val(Text33.Text) '排气过程 pr = Int(1000 * 10.8 * po) / 1000: tr = 350 / (1.2 / Log(n) * Log(10

10、) + 0.005 * (e - 3) + 0.01 * (fa - 1): hu = 42500 '柴油机 暂时先赋的值 '进气过程 paa = n 2 / 520000000 pac = (e - 0.5) 2) / (e - 1) 2) '残余废气收缩系数 暂时取的0.5 T = 20 't 进气温升 取的 20 pa = 0.9 * po * 10 '进气压力 yy = (tt + T) * pr / tr / (e * pa - pr) ' 残余废气系数 yy ta = (tt + T + yy * tr) / (1 + yy) 

11、9;进气温度 nv = tt * pa * e * (e * pa - pr) / (tt + T) / (10 * po) / (e - 1) / e / pa '充气系数 '压缩过程 '平均压缩多变指数 pc = Int(100 * (pa * e (n1) / 100 '压缩终点压力 tc = Int(10 * ta * e (n1 - 1) / 10 '压缩终点温度 '燃烧过程 lo = 1 / 0.21 * (gc / 12 + gh / 4 - g0 / 32) m1 = A * lo: m2 = A * lo + gh / 4 +

12、g0 / 32 u0 = m2 / m1 u = (u0 + yy) / (1 + yy) '实际分子变更系数 '化学损失 khu = 0 cv1 = (4.815 + 0.000415 * tc) * 4.1868 qq = Val(Text30.Text) '-qq 热量利用系数选取 ccc = (qq * (hu - khu) / m1 / (1 + yy) + cv1 * tc) / u aaa = (3.3 / A + 3.7) * 4.1868 / 10000 '柴油机 bbb = (4.8 + 2.2 / A) * 4.1868 tz = Int(

13、-bbb + (bbb 2 + 4 * aaa * ccc) 0.5) / 2 / aaa) '燃烧终点温度 ylsgb = Val(Text29.Text) p = u * tz / ylsgb / tc '-柴油机预期膨胀比- pz = Int(100 * (ylsgb * pc) / 100 '最高燃烧压力 pb = Int(100 * (pz / (e n2) / 100 '膨胀终点压力 tb = Int(tz / (e (n2 - 1) '膨胀终点温度 pi1 = pc / (e - 1) * (ylsgb * p / (n2 - 1) * (

14、1 - 1 / e (n2 - 1) - 1 / (n1 - 1) * (1 - 1 / e (n1 - 1) pi = fi * pi1 '平均指示压力 pm = 0.1 * (Val(Text32.Text) 0.5 * (1 + n / 1000) '机械损失压力 nm = 1 - pm / pi ' 机械效率 ni = Int(10000 * 8.314 * m1 * tt * pi / hu / 10 / po / nv) / 10000 '指示热效率 gi = 3600000 / hu / ni '指示燃油消耗率 ne = nm * ni &

15、#39;有效热效率 pe = pi * nm '平均有效压力 ge = 3600000 / hu / ne '有效燃油消耗率 vh = Int(10000 * 3.14 * (d / 2) 2 * r / 10 6) / 10000 v = Int(10000 * (vh / e + vh) / 10000 nne = i * v * pe * n / 1224 Text7(1).Text = "0" & Int(10000 * ni) / 10000 Text7(3).Text = Int(100 * gi) / 100 Text7(4).Text

16、 = "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 Text10.Text = pr: Text13.Text = Int(10 * tr) / 10: Text21.Text = "0" & Int(100 * pa) / 100: Text22.Text = Int(ta * 10) / 10 Text

17、23.Text = pz: Text24.Text = tz: Text25.Text = tb: Text26.Text = pb Text27.Text = pc: Text28.Text = tceee:End SubPrivate Sub Command5_Click() Picture1.Cls Picture1.Scale (-10, 120)-(780, -10) Picture1.Line (0, 0)-(725, 0) Picture1.Line (0, 0)-(0, 115) vc = vh / (e - 1) Picture1.Line (0, po)-(720, po)

18、, RGB(0, 220, 0) Picture1.CurrentX = 680: Picture1.CurrentY = 8: Picture1.Print "(°AC)" Picture1.CurrentX = 10: Picture1.CurrentY = 129: Picture1.Print "P (bar)" For Y = 10 To 110 Step 10# Picture1.Line (0, Y)-(7, Y) Picture1.CurrentX = -0.1: Picture1.CurrentY = Y: Picture1.

19、Print Y Next Y Picture1.Line (0, 115)-(-5, 112) Picture1.Line (0, 115)-(5, 112) For X = 0 To 720 Step 60 Picture1.Line (X, 0)-(X, 2) Picture1.CurrentX = X - 20: Picture1.CurrentY = -2: Picture1.Print X Next X Picture1.Line (725, 0)-(716, 2) Picture1.Line (725, 0)-(716, -2) If vh = 0 Then MsgBox &quo

20、t;请先输入参数,并计算", 0, "提示": GoTo q End If Va = vh + vh / (e - 1) vc = vh / (e - 1) For x = 3.14 To 6.28 Step 0.0001 vcx = vh / 2 * (1 - Cos(x) + (1 - Cos(2 * x) * r * l / 4) + vc '压缩过程 pcx = pa * (Va / vcx) n1 = x * 180 / 3.14 Picture1.PSet (, pcx), RGB(0, 400, 200) Next x Va = vh + v

21、h / (e - 1) vc = vh / (e - 1) For x = 9.42 To 6.28 Step -0.0001 vbx = vh / 2 * (1 - Cos(x) + (1 - Cos(2 * x) * r * l / 4) + vc pbx = pb * (Va / vbx) n2 = x * 180 / 3.14 If pbx <= pz Then Picture1.PSet ( + 20, pbx), RGB(400, 0, 220) Else: Picture1.Line (360, pz)-(, pz), RGB(400, 100, 0) '做功过程

22、Picture1.Line (360, pc)-(360, pz), RGB(0, 0, 220) End If Next x Picture1.Line (540, 2.11)-(720, 2.11) Picture1.Line (360, 54.202)-(360, 86.72), RGB(100, 0, 220) Picture1.Line (360, 86.72)-(380, 86.72), RGB(0, 100, 220)q:End SubPrivate Sub Command6_Click() Dim T, ex, A, C, bx, ea, b Picture1.Cls Pict

23、ure1.Scale (-60, 2600)-(780, -150) Picture1.Line (0, 0)-(740, 0) Picture1.Line (0, 0)-(0, 2500) vc = vh / (e - 1) Picture1.CurrentX = 680: Picture1.CurrentY = 120: Picture1.Print "(°AC)" Picture1.CurrentX = 10: Picture1.CurrentY = 2500: Picture1.Print "T(K)" For Y = 200 To 2

24、200 Step 200# Picture1.Line (0, Y)-(7, Y) Picture1.CurrentX = -60: Picture1.CurrentY = Y + 30: Picture1.Print Y Next Y Picture1.Line (0, 2500)-(-5, 2400) Picture1.Line (0, 2500)-(5, 2400) For X = 0 To 720 Step 60 Picture1.Line (X, 0)-(X, 40) Picture1.CurrentX = X - 20: Picture1.CurrentY = -30: Pictu

25、re1.Print X Next X Picture1.Line (740, 0)-(733, 25) Picture1.Line (740, 0)-(733, -25) If vh = 0 Then MsgBox "请先输入参数,并计算", 0, "提示": GoTo q End If b = Val(Text6.Text) For x = 3.14 * 2 To 3.11 * 2.35 Step 0.0001 A = 1 - Cos(x) + 0.5 * 0.25 * Sin(x) 2 C = 1 / (e - 1) ea = e ex = 1 +

26、A / (2 * C) bx = b - (352 - x * 180 / 3.14) / (352 - 400) 3 * (b - 1) T = (ta * ex * 10) / (pa * ea * bx) = x * 180 / 3.14 Picture1.PSet ( - 120, T), RGB(220, 0, 0) Next x For x = 3.14 To 2 * 3.14 Step 0.0001 A = 1 - Cos(x) + 0.5 * 0.25 * Sin(x) 2 C = 1 / (e - 1) ea = e ex = 1 + A / (2 * C) bx = b -

27、 (352 - x * 180 / 3.14) / (352 - 400) 3 * (b - 1) T = (ta * ex * 10) / (pa * ea * bx) = x * 180 / 3.14 Picture1.PSet ( + 120, T), RGB(0, 0, 220) Next xq:End SubPrivate Sub Command7_Click() Dim T, ex, A, Cm, bx, ea, b, Cu, C1, C2, C, h1, h2, H3 Picture1.Cls Picture1.Scale (-60, 2000)-(780, -150) Pict

28、ure1.Line (0, 0)-(740, 0) Picture1.Line (0, 0)-(0, 1900) vc = vh / (e - 1) Picture1.CurrentX = 680: Picture1.CurrentY = 120: Picture1.Print "(°AC)" Picture1.CurrentX = 10: Picture1.CurrentY = 1900: Picture1.Print "h：W/(m2.K)" For Y = 200 To 1800 Step 200# Picture1.Line (0, Y

29、)-(7, Y) Picture1.CurrentX = -60: Picture1.CurrentY = Y + 30: Picture1.Print Y Next Y Picture1.Line (0, 1900)-(-5, 1850) Picture1.Line (0, 1900)-(5, 1850) For X = 0 To 720 Step 60 Picture1.Line (X, 0)-(X, 40) Picture1.CurrentX = X - 20: Picture1.CurrentY = -30: Picture1.Print X Next X Picture1.Line

30、(740, 0)-(733, 25) Picture1.Line (740, 0)-(733, -25) If vh = 0 Then MsgBox "请先输入参数,并计算", 0, "提示": GoTo q End If C = Val(Text34.Text) Cm = s * 2 / 30: C1 = 2.28 + 0.308 * C: C2 = 3.2 * 10 (-3) d = Val(Text4.Text) A = vh * ta / (pa * 0.03719) For x = 3.14 To 6.28 Step 0.0001 Va = v

31、h + vh / (e - 1) vc = vh / (e - 1) vcx = vh / 2 * (1 - Cos(x) + (1 - Cos(2 * x) * r * l / 4) + vc pcx = pa * (Va / vcx) n1 Va = vh + vh / (e - 1) vc = vh / (e - 1) A = 1 - Cos(x) + 0.5 * 0.25 * Sin(x) 2 C = 1 / (e - 1) ea = e ex = 1 + A / (2 * C) bx = 1.04 - (352 - x * 180 / 3.14) / (352 - 400) 3 *

32、(1.04 - 1) T = (ta * ex * 10) / (pa * ea * bx) H = 130 * (d * 0.001) (-0.2) * (p * 10) (0.8) * T (-0.53) * (C1 * Cm + C2 * A * (p - 1) 0.8 = x * 180 / 3.14 Picture1.PSet (, H), RGB(220, 0, 0) H = H + 456.9 Next x h1 = H / 3.14 For x = 6.28 To 450 * 3.14 / 180 Step 0.0001 Va = vh + vh / (e - 1) vc =

33、vh / (e - 1) vcx = vh / 2 * (1 - Cos(x) + (1 - Cos(2 * x) * r * l / 4) + vc pcx = pa * (Va / vcx) n1 Va = vh + vh / (e - 1) vc = vh / (e - 1) A = 1 - Cos(x) + 0.5 * 0.25 * Sin(x) 2 C = 1 / (e - 1) ea = e ex = 1 + A / (2 * C) bx = 1.04 - (352 - x * 180 / 3.14) / (352 - 400) 3 * (1.04 - 1) T = (ta * e

34、x * 10) / (pa * ea * bx) H = 130 * (d * 0.001) (-0.2) * (p * 10) (0.8) * T (-0.53) * (C1 * Cm + C2 * A * (p - 1) 0.8 = x * 180 / 3.14 Picture1.PSet (, H), RGB(0, 0, 220) H = H + 1302.56 Next x h2 = H / 1.07 H3 = (h1 + h2) / 2 Text35.Text = Int(H3) q:End SubPrivate Sub Command8_Click() Dim tmax, C1,

35、C2, f1, t1, tmax1, ti, fx, fy, dyi, yi, xi Picture1.Cls Picture1.Scale (-10, 130)-(130, -10) Picture1.Line (0, 100)-(110, 100): Picture1.Line (52.5, 120)-(52.5, 20) Picture1.CurrentX = 110: Picture1.CurrentY = 105: Picture1.Print "+X" Picture1.CurrentX = 0: Picture1.CurrentY = 105: Picture

36、1.Print "-X" Picture1.CurrentX = 55: Picture1.CurrentY = 120: Picture1.Print "-Y" Picture1.CurrentX = 55: Picture1.CurrentY = 20: Picture1.Print "+Y" Picture1.Line (110, 100)-(108, 98): Picture1.Line (110, 100)-(108, 102) Picture1.Line (52.5, 20)-(50.5, 22): Picture1.Li

37、ne (52.5, 20)-(54.5, 22) 'Picture1.Line (0, 0)-(110, 0) ' Picture1.Line (0, 0)-(0, 110) ' vc = vh / (e - 1) ' Picture1.CurrentX = 680: Picture1.CurrentY = 120: Picture1.Print "(°AC)" ' Picture1.CurrentX = 10: Picture1.CurrentY = 1900: Picture1.Print "h：W/(m2.K)" ' 'For Y = 0 To 105 Step 5 ' Pic