发动机工作过程计算

第 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 = tav = * 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 If0 = m2 / m1 '理论分子变更系数 = (0 + ) / (1 + ) '实际分子变更系数 Text34.Text = If < 1 Thenhu = 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.Text ne = 10 * i * vh * pe * n / 1200 Text44.Text = neEnd SubPrivate Sub Command3_Click() MsgBox ("确定退出？") EndEnd SubPrivate Sub Command4_Click() Picture1.ClsEnd SubStatic Function Log10(X) Log10 = Log(X) / Log(10#) '求以10为底的对数用的函数End FunctionPrivate Sub Form_Load() Option1.Value = True Option4.Value = TrueEnd SubPrivate Sub Option1_Click() jx = 1End SubPrivate Sub Option2_Click() jx = 2End SubPrivate Sub Option3_Click() zy = 2 Label45.Visible = True Text45.Visible = TrueEnd SubPrivate Sub Option4_Click() zy = 1 Label45.Visible = False Text45.Visible = FalseEnd Sub运行示例如下： 1、 当选择是汽油机且非增压时，如图：运算结果如下：此时的示功图如下：2、 当选择是柴油机且非增压时，如图：运算结果如下：此时的示功图如下：1、 当选择是柴油机且增压时，如图：运算结果如下：此时的示功图如下：4. 实验总结17 / 17