柴油机工作过程计算..

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

2、z, ccc, hu, qq, aaa, bbb As SingleDim l8, n2, pz, H, tb, pi1, pi, x, , p As SingleDim fi, pm, nm, ni, gi, ne, pe, ge, nne, i, v, pb, pcx, vbx, pbx As DoubleDim aa(50), ylsgb As SinglePrivate Sub Command4_Click()Picture1.ClsPicture1.ForeColor = vbGreen Picture1.DrawWidth = 1.5 Picture1.Scale (-0.1, 9

3、9)-(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 1Picture1.Line (0, i * 10)-(0.015, i * 10) Picture1.Line

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

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

6、c, po)-(vc, 88)Picture1.DrawWidth = 1.5Picture1.ForeColor = vbYellowFor X = v To vc Step -0.0001vcx = Xpcx = pa * (v / vcx) n1Picture1.PSet (vcx, pcx)Next X'燃烧过程曲线Picture1.Line (vcx, pcx)-(vc, pz) '膨胀过程曲线Dim vbx, pbx As SinglePicture1.Line (vc, pz)-(vc * p, pz): vc = vc * pFor X = vc To v St

7、ep 0.0001发动机热负荷与热强度实验报告vbx = X pbx = pz * (vc / vbx) n2)Picture1.PSet (vbx, pbx)Next X'排气过程曲线Picture1.Line (vbx, pbx)-(v, pa)Picture1.DrawWidth = 1Picture1.ForeColor = vbRedPicture1.CurrentX = 0.1: Picture1.CurrentY = 8: Picture1.Print "po"Picture1.CurrentX = 0.02: Picture1.CurrentY =

8、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.Text) / 3.14 / (d / 2) 2)gc = Val(Text11

9、.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) + 0.005 * (e - 3) + 0.01 * (fa- 1): hu

10、= 42500 ' 柴油机 暂时先赋的值'进气过程发动机热负荷与热强度实验报告paa = n 2 / 520000000 pac = (e - 0.5) 2) / (e - 1) 2)'残余废气收缩系数 暂时取的 0.5T = 20't 进气温升pa = 0.9 * po * 10yy = (tt + T) * pr / tr / (e * pa - pr) ta = (tt + T + yy * tr) / (1 + yy)取的 20'进气压力' 残余废气系数 yy'进气温度'压缩终点压力'压缩终点温度'实际分子

11、变更系数nv = tt * pa * e * (e * pa - pr) / (tt + T) / (10 * po) / (e - 1) / e / pa ' 充气系数'压缩过程'平均压缩多变指数pc = Int(100 * (pa * e (n1) / 100tc = 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 + g0 / 32 u0 = m2 / m1u = (u0 + yy)

12、 / (1 + yy)'化学损失khu = 0 cv1 = (4.815 + 0.000415 * tc) * 4.1868qq = Val(Text30.Text) 'qq 热量利用系数选取ccc = (qq * (hu - khu) / m1 / (1 + yy) + cv1 * tc) / uaaa = (3.3 / A + 3.7) * 4.1868 / 10000'柴油机bbb = (4.8 + 2.2 / A) * 4.1868发动机热负荷与热强度实验报告tz = Int(-bbb + (bbb 2 + 4 * aaa * ccc) 0.5) / 2 / a

13、aa) 燃' 烧终点温度ylsgb = Val(Text29.Text) p = u * tz / ylsgb / tcpz = 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) * (1 - 1 / e (n2 - 1) - 1 / (n1 - 1) * (1 - 1 / e (n1

14、 - 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'有效热效率pe = pi * nm'平均有效压力ge = 3600000 / hu / ne&

15、#39;有效燃油消耗率vh = Int(10000 * 3.14 * (d / 2) 2 * r / 10 6) / 10000v = Int(10000 * (vh / e + vh) / 10000nne = i * v * pe * n / 1224Text7(1).Text = "0" & Int(10000 * ni) / 10000Text7(3).Text = Int(100 * gi) / 100Text7(4).Text = "0" & Int(10000 * ne) / 10000Text7(5).Text = Int

16、(1000 * pe) / 1000Text7(6).Text = Int(100 * ge) / 100Text7(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) / 10Text23.Text = pz: Text24.Text = tz: Text25.Text = tb: Text26.

17、Text = pbText27.Text = pc: Text28.Text = tc eee: End SubPrivate Sub Command5_Click()Picture1.ClsPicture1.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), RGB(0, 220, 0)(°AC)Picture1.CurrentX = 680: Picture1.C

18、urrentY = 8: Picture1.Print "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.Print Y Next YPicture1.Line (0, 115)-(-5, 112)Picture1.Line (0, 115)-(5

19、, 112)For X = 0 To 720 Step 60Picture1.Line (X, 0)-(X, 2)Picture1.CurrentX = X - 20: Picture1.CurrentY = -2: Picture1.Print X Next XPicture1.Line (725, 0)-(716, 2)Picture1.Line (725, 0)-(716, -2)If vh = 0 Then发动机热负荷与热强度实验报告MsgBox "请先输入参数 ,并计算", 0, "提示": GoTo qEnd IfVa = vh + vh /

20、 (e - 1)vc = vh / (e - 1)For x = 3.14 To 6.28 Step 0.0001vcx = vh / 2 * (1 - Cos( x) + (1 - Cos(2 * x) * r * l / 4) + vc ' 压缩过程 pcx = pa * (Va / vcx) n1 = x * 180 / 3.14Picture1.PSet ( , RpcGxB),(0, 400, 200)Next xVa = vh + vh / (e - 1)vc = vh / (e - 1)For x = 9.42 To 6.28 Ste-0p. 0001vbx = vh /

21、 2 * (1 - Cos( x) + (1- Cos(2 * x) * r * l / 4) + vc pbx = pb * (Va / vbx) n2 = x * 1/8 30. 14If pbx <= pz ThenPicture1.PSet ( + 20, pbx), RGB(400, 0, 220)Else: Picture1.Line (360, pz)-(, pz), RGB(400, 100, 0)'做功过程Picture1.Line (360, pc)-(360, pz), RGB(0, 0, 220)End IfNext xPicture1.Line (540

22、, 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, bPicture1.ClsPicture1.Scale (-60, 2600)-(780, -150)Picture1.Line (0, 0)-(740, 0)Picture1.

23、Line (0, 0)-(0, 2500) vc = vh / (e - 1)(°AC)Picture1.CurrentX = 680: Picture1.CurrentY = 120: Picture1.Print "Picture1.CurrentX = 10: Picture1.CurrentY = 2500: Picture1.Print "T(K)"For Y = 200 To 2200 Step 200#Picture1.Line (0, Y)-(7, Y)Picture1.CurrentX = -60: Picture1.CurrentY

24、= Y + 30: Picture1.Print Y Next YPicture1.Line (0, 2500)-(-5, 2400)Picture1.Line (0, 2500)-(5, 2400)For X = 0 To 720 Step 60Picture1.Line (X, 0)-(X, 40)Picture1.CurrentX = X - 20: Picture1.CurrentY = -30: Picture1.Print X Next XPicture1.Line (740, 0)-(733, 25)Picture1.Line (740, 0)-(733, -25)If vh =

25、 0 ThenMsgBox "请先输入参数 ,并计算", 0, "提示": GoTo qEnd Ifb = Val(Text6.Text)For x = 3.14 * 2 To 3.11 * 2.35 Step 0.0001 x) 2A = 1 - Cos(x) + 0.5 * 0.25 * Sin(C = 1 / (e - 1)ea = e发动机热负荷与热强度实验报告ex = 1 + A / (2 * C)bx = b - (352 - x * 180 / 3.14) / (352- 400) 3 * (b - 1)T = (ta * ex * 10)

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

27、0 / 3.14Picture1.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, H3Picture1.ClsPicture1.Scale (-60, 2000)-(780, -150)Picture1.Line (0, 0)-(740, 0)Picture1.Line (0, 0)-(0, 1900) vc = vh / (e - 1)(°AC)Picture1.Current

28、X = 680: Picture1.CurrentY = 120: Picture1.Print "Picture1.CurrentX = 10: Picture1.CurrentY = 1900: Picture1.Print "h： W/(m2.K)"发动机热负荷与热强度实验报告For Y = 200 To 1800 Step 200#Picture1.Line (0, Y)-(7, Y)Picture1.CurrentX = -60: Picture1.CurrentY = Y + 30: Picture1.Print Y Next YPicture1.Li

29、ne (0, 1900)-(-5, 1850)Picture1.Line (0, 1900)-(5, 1850)For X = 0 To 720 Step 60Picture1.Line (X, 0)-(X, 40)Picture1.CurrentX = X - 20: Picture1.CurrentY = -30: Picture1.Print X Next XPicture1.Line (740, 0)-(733, 25)Picture1.Line (740, 0)-(733, -25)If vh = 0 ThenMsgBox "请先输入参数 ,并计算", 0, &q

30、uot;提示": GoTo qEnd IfC = 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.0001Va = vh + vh / (e - 1) vc = vh / (e - 1)vcx = vh / 2 * (1 - Cos( x) + (1- Cos(2 * x) * r * l / 4) + vc pcx = pa * (Va

31、 / vcx) n1发动机热负荷与热强度实验报告Va = vh + vh / (e - 1) vc = vh / (e - 1) x) 2A = 1 - Cos(x) + 0.5 * 0.25 * Sin(C = 1 / (e - 1)ea = eex = 1 + A / (2 * C)bx = 1.04 - (352 - x * 180 / 3.14) / (352- 400) 3 * (1.04 - 1) T = (ta * ex * 10) / (pa * ea * bx)H = 130 * (d * 0.001) (-0.2) * (p * 10) (0.8) * T (-0.53)

32、* (C1 * Cm + C2 * A * (p - 1) 0.8 = x * 180 / 3.14, H), RGB(220, 0, 0)Picture1.PSet (H = H + 456.9Next x h1 = H / 3.14For x = 6.28 To 450 * 3.14 / 180 Step 0.0001 Va = vh + vh / (e - 1)vc = vh / (e - 1)vcx = vh / 2 * (1 - Cos( x) + (1- Cos(2 * x) * r * l / 4) + vcpcx = pa * (Va / vcx) n1Va = vh + vh

33、 / (e - 1)发动机热负荷与热强度实验报告vc = vh / (e - 1) x) 2A = 1 - Cos(x) + 0.5 * 0.25 * Sin(C = 1 / (e - 1)ea = eex = 1 + A / (2 * C)bx = 1.04 - (352 - x * 180 / 3.14) / (352- 400) 3 * (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 * (

34、p - 1) 0.8 = x * 180 / 3.14Picture1.PSet ( , H), RGB(0, 0, 220)H = H + 1302.56Next xh2 = H / 1.07H3 = (h1 + h2) / 2Text35.Text = Int(H3) q:End SubPrivate Sub Command8_Click()Dim tmax, C1, C2, f1, t1, tmax1, ti, fx, fy, dyi, yi, xiPicture1.ClsPicture1.Scale (-10, 130)-(130, -10)发动机热负荷与热强度实验报告Picture1

35、.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: Picture1.Print "-X"Picture1.CurrentX = 55: Picture1.CurrentY = 120: Picture1.Print "-Y"Pic

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