发动机工作过程计算.doc

西华大学实验报告第 2 组西华大学实验报告开课学院及实验室： 实验时间：2012 年 4 月 7 日学 生 姓 名学 号 实 验 成 绩课 程 名 称发动机工作过程的计算课 程 代 码8200250实验项目名称发动机工作过程计算上机实验项 目 代 码指 导 教 师项 目 学 分实验课考勤实验工作表现实验报告1、实验目的2、实验设备、仪器及材料3、实验内容3.1 一般实验（非上机实验）：3.1.1实验方案设计与选择（设计性、创新性实验需写该项内容）3.1.2实验原理及实验步骤（实验工作原理或实验的主要操作过程）3.1.3实验记录（核心代码及调试过程）3.2 上机实验：3.2.1上机实验的内容及要求3.2.2算法设计思想与算法实现步骤3.2.3程序核心代码，程序调试过程中出现的问题及解决方法3.2.4 程序运行的结果注解：理工科实验需记录实验过程中的数据、图表、计算、现象观察等，实验过程中出现的问题；其它如在计算机上进行的编程、仿真性或模拟性实验需记录程序核心代码以及程序在调式过程中出现的问题及解决方法；记录程序执行的结果。4、实验总结4.1实验结果分析及问题讨论4.2实验总结心得体会注解：实验总结的内容根据不同学科和类型实验要求不一样，一般理工科类的实验需要对实验结果进行分析，并且对实验过程中问题进行讨论；在计算机上进行的编程、仿真性或模拟性实验需要对上机实践结果进行分析，上机的心得体会及改进意见。其它实验应总结实验过程写出心得体会及改进意见。说明：各门实验课程实验报告的格式及内容要求，请按照实验指导书的要求手工书写。1. 实验目的2. 实验仪器，设备及材料3. 实验内容 MsgBox 西华大学 2009级 交通与汽车工程学院 热能与动力工程（汽车发动机） 3班 唐俊End Sub运行结果依次如下： 进入界面时：Private Sub Command1_Click() Form1.Show Form2.HideEnd Sub运行出的图像如下：Private Sub Command2_Click() MsgBox 西华大学 2009级 交通与汽车工程学院 热能与动力工程（汽车发动机） 3班 唐俊End Sub当进入主界面之后可对发动机类型，是否需要增压，经验参数以及发动机参数进行选择和修改。选择界面如下：主程序的详细代码如下：Dim As Single, S As Single, n As Integer, nh As Integer, As Single, vc, vh, pa, n1, , pz, vz, n2, jx, zyPrivate Sub Command1_Click() Picture1.Scale (-0.05, -0.5)-(1.05, 11) Picture1.Line (0, 10)-(1, 10) 画X轴 Picture1.Line (0.98, 9.9)-(1, 10) Picture1.Line (0.98, 10.1)-(1, 10) Picture1.CurrentX = 0.98: Picture1.CurrentY = 10.2: Picture1.Print V/L For i = 0 To 0.9 Step 0.1 Picture1.Line (i, 9.9)-(i, 10) Picture1.Print 0 & i Next i Picture1.Line (0, 10)-(0, 0) 画Y轴 Picture1.Line (0, 0)-(-0.01, 0.3) Picture1.Line (0.01, 0.3)-(0, 0) Picture1.CurrentX = -0.04: Picture1.CurrentY = -0.4: Picture1.Print P/Mpa For i = 9 To 1 Step -1 Picture1.Line (0.01, i)-(0, i) Picture1.Print 10 - i Next i Picture1.ForeColor = vbRed 画P0线 Picture1.Line (0, 9.9)-(1, 9.9) Picture1.CurrentX = 1: Picture1.CurrentY = 9.4: Picture1.Print P0 Picture1.ForeColor = vbBlue 画Vc线 压缩曲线 = Val(Text2.Text) d = Text4.Text S = Text3.Text r = S / 2 vh = 3.14 * d 2 * S / 4 * 10 -6 Text16.Text = vh vc = vh / ( - 1) Picture1.CurrentX = vc + 0.01: Picture1.CurrentY = 1: Picture1.Print Vc Picture1.Line (vc, 10)-(vc, 0) Va = vc + vh l = Text15.Text Picture1.ForeColor = vbGreen For x = 3.14 To 0 Step -0.0001 vcx = vh / 2 * (1 - Cos(x) + (1 - Cos(2 * x) * Rl / 4) + vc vcx = (3.14 * d 2 * r * (1 - Cos(x) + r * (1 - Cos(2 * x) / 4 / l) / 4) * 10 -8 + vc pcx = pa * (Va / vcx) n1 Picture1.PSet (vcx, 10 - pcx) Next x 燃烧过程曲线 Picture1.Line (vcx, 10 - pcx)-(vcx, 10 - pcx * ) 膨胀过程曲线 For x = 0 To 3.14 Step 0.0001 vbx = vh / 2 * (1 - Cos(x) + (1 - Cos(2 * x) * Rl / 4) + vc vcx = (3.14 * d 2 * r * (1 - Cos(x) + r * (1 - Cos(2 * x) / 4 / l) / 4 + vc) * 10 -6 pbx = pz * (vc / vbx) n2 Picture1.PSet (vbx, 10 - pbx) Next x Picture1.Line (vbx, 10 - pbx)-(vbx, 10 - pa) Picture1.ForeColor = vbBlackEnd SubPrivate Sub Command2_Click() Picture1.Cls = Val(Text2.Text) d = Text4.Text S = Text3.Text r = S / 2 l = Text15.Text vh = 3.14 * d 2 * S / 4 * 10 -6 Text16.Text = vh 排气过程 n = Val(Text5.Text) nh = Val(Text6.Text) If zy = 1 Then pr = 0.1 + 0.03 * n / nh ElseIf zy = 2 Then pk = Text45.Text pr = 2 * 10 -8 * n 2 pr = pk + pr End If Text18.Text = pr = Val(Text1.Text) = Val(Text2.Text) If = 1 Then MsgBox (压缩比不能为“1”，请重新输入！) Text2.Text = End If Tr = 350 / (1.2 / Log10(n) + 0.005 * ( - 3) + 0.01 * ( - 1) 调用自定义函数就以10为底的对数 Text17.Text = Tr 进气过程 T = 10 * (110 - 0.0125 * n) / (110 - 0.0125 * nh) 10为Tn,选取 = 0.75 = 0.5 k = 1.35 a = vh / f p0 = Text7.Text Pa = P0 * (1 - n 2 / 2*9.8*k*p0 / 2 * ( - ) / ( - 1) 2) k / (k - 1) pa = Int(10 4 * (p0 * (1 - (10 * n 2 / 520 / 10 6) / 10.6 2 / 0.75 2) * ( - ) 2 / ( - 1) 2) (k / (k - 1) / 10 4 Text20.Text = pa t0 = Text8.Text = Int(10 4 * (t0 + T) * pr / Tr / ( * pa - pr) / 10 4 残余废气系数 Text21.Text = ta = (t0 + T + * * Tr * (pa / pr) (n - 1) / n) / (1 + * ) Text19.Text = ta v = * pa * t0 / ( - 1) / p0 / ta / (1 + ) 充气效率 Text22.Text = v 压缩过程 n1 = 1.46 - 0.05 * nh / n Text25.Text = n1 pc = pa * n1 Text29.Text = pc tc = ta * (n1 - 1) Text28.Text = tc 燃烧过程： gc = Val(Text9.Text) gh = Val(Text10.Text) go = Val(Text11.Text) gs = Val(Text12.Text) l0 = 1 / 0.21 * (gc / 12 + gh / 4 - go / 32) 理论空气量 mt = Text13.Text If jx = 1 Then m1 = * l0 + 1 / mt 汽油机新鲜充量 ElseIf jx = 2 Then m1 = * l0 End If If 1 Then m2燃烧产物总量 m2 = * l0 + gh / 4 + go / 32 ElseIf 1 Then m2 = 0.79 * * l0 + gc / 12 + gh / 2 End If 0 = m2 / m1 理论分子变更系数 = (0 + ) / (1 + ) 实际分子变更系数 Text34.Text = If 1 Then hu = 58000 * (1 - ) 化学损失 End If Tz的计算 = Text23.Text If jx = 1 Then hu = 44100 ElseIf jx = 2 Then hu = 42500 End If q1 = * (hu - hu) 加热工质的热量 cv1 = 4.815 + 0.415 * 10 -3 * (tc + t0) / 2 Select Case Case 0.8 To 1 If jx = 1 Then aa = (3.7 + 3.3 * ) * 10 -4 * bb = (4.4 + 0.62 * ) * cc = -(q1 / m1 / (1 + ) + cv1 * tc) tz = (-bb + (bb 2 - 4 * aa * cc) 0.5) / 2 / aa / 2.5 Text31.Text = tz = * tz / tc Text33.Text = pz = * pc Text32.Text = pz ElseIf jx = 2 Then = Text33.Text aa = * (3.7 + 3.3 * ) * 10 -4 bb = (6.385 + 0.62 * ) * cc = -( * hu / m1 / (1 + ) + (cv1 + 1.98 * ) * tc) tz = (-bb + (bb 2 - 4 * aa * cc) 0.5) / 2 / aa / 2.5 Text31.Text = tz pz = * pc Text32.Text = pz = * tz / / tc Text30.Text = End If Case 1 To 2 If jx = 1 Then aa = (3.7 + 3.3 / ) * 10 -4 * bb = (4.8 + 0.22 / ) * cc = -(q1 / m1 / (1 + ) + cv1 * tc) tz = (-bb + (bb 2 - 4 * aa * cc) 0.5) / 2 / aa / 2.5 Text31.Text = tz = * tz / tc Text33.Text = pz = * pc Text32.Text = pz ElseIf jx = 2 Then = Text33.Text aa = (3.7 + 3.3 / ) * 10 -4 * bb = (6.785 + 0.22 / ) * cc = -( * hu / m1 / (1 + ) + (cv1 + 1.98 * ) * tc) tz = (-bb + (bb 2 - 4 * aa * cc) 0.5) / 2 / aa / 2.5 Text31.Text = tz pz = * pc Text32.Text = pz = * tz / / tc Text30.Text = End If End Select 膨胀过程： n2 = 1.2 + 0.039 * nh / n Text26.Text = n2 = / Text37.Text = pb = pz / n2 Text36.Text = pb tb = tz / (n2 - 1) Text35.Text = tb 发动机性能分析 = 1 Text30.Text = = 0 Text37.Text = Pii = pc * ( * (1 - 1 / (n2 - 1) / (n2 - 1) - (1 - 1 / (n1 - 1) / (n1 - 1) / ( - 1) Pii = pa * n1 * ( * ( - 1) + * / (n2 - 1) * (1 - 1 / (n2 - 1) - (1 - 1 / (n1 - 1) / (n1 - 1) / ( - 1) i = Text24.Text pi = i * Pii Text39.Text = pi i = 8.314 * m1 * t0 * pi / hu / p0 / v Text38.Text = i gi = 3.6 * 10 6 / hu / i Text40.Text = gi pm = 0.04 * (1 + n / 1000) m = 1 - pm / pi Text27.Text = m e = i * m Text42.Text = e pe = pi * m Text41.Text = pe ge = 3.6 * 10 6 / hu / e Text43.Text = ge i = Text14.