VB绘制汽车驱动力与行驶阻力平衡图,汽车功率平衡图。
3.2 实验记录（核心代码及调试过程）
1.核心代码：
'.............驱动力与行驶阻力平衡图
Private Sub Command1_Click()
Picture1.DrawWidth = 1
Picture1.Cls
Dim i, j, n As Integer
Dim Ft, Ff, Fw, Umax, Uumax As Single
Picture1.Scale (-20, 18000)-(150, -2000)
Picture1.Line (0, 0)-(125, 0): Picture1.Line (0, 0)-(0, 16000)
Picture1.CurrentX = 125: Picture1.CurrentY = -300: Picture1.Print "Ua(km/h)"
Picture1.CurrentX = 125: Picture1.CurrentY = 192: Picture1.Print "→"
Picture1.CurrentX = 5: Picture1.CurrentY = 16000: Picture1.Print "F/N"
Picture1.CurrentX = -2: Picture1.CurrentY = 16000: Picture1.Print "↑"
Picture1.CurrentX = -5: Picture1.CurrentY = -1: Picture1.Print "0"
For i = 20 To 120 Step 20
Picture1.Line (i, 0)-(i, 200)
Picture1.CurrentX = i - 5: Picture1.CurrentY = -10: Picture1.Print i
Next i
For i = 1000 To 15000 Step 1000
Picture1.Line (0, i)-(1.5, i)
Picture1.CurrentX = -15: Picture1.CurrentY = i + 50: Picture1.Print i
Next i
Picture1.DrawWidth = 2
For n = 600 To 4000
Tq = -19.313 + 295.27 * (n / 1000) - 165.44 * (n / 1000) ^ 2 + 40.874 * (n / 1000) ^ 3 - 3.8445 * (n / 1000) ^ 4
For j = 1 To 5
Ft = Tq * i0 * ηt * x(j) / r
Ua = 0.377 * r * n / (i0 * x(j))
Fw = 0.5 * 1.2258 * ca * (Ua / 3.6) ^ 2
Picture1.Circle (Ua, Ft), 0.1, vbBlue
Picture1.Circle (Ua, (Ff + Fw)), 0.1, vbRed
Umax = Umax
If Ua > Umax Then
Umax = Ua
Else
Umax = Umax
End If
Next j
Next n
Uumax = CInt(Umax * 100) / 100
Picture1.CurrentX = 40: Picture1.CurrentY = 16000: Picture1.Print "最高车速为："; Uumax; "km/h"
End Sub
'...........动力特性图
Private Sub Command3_Click()
Picture1.DrawWidth = 1
Picture1.Cls
Dim i, j, n As Integer
Dim Ft, Fw, D As Single
Picture1.Scale (-20, 90)-(150, -10)
Picture1.Line (0, 0)-(125, 0): Picture1.Line (0, 0)-(0, 85)
Picture1.CurrentX = 125: Picture1.CurrentY = -3: Picture1.Print "Ua(km/h)"
Picture1.CurrentX = 125: Picture1.CurrentY = 0.9: Picture1.Print "→"
Picture1.CurrentX = 5: Picture1.CurrentY = 85: Picture1.Print "D"
Picture1.CurrentX = -2: Picture1.CurrentY = 85: Picture1.Print "↑"
For i = 20 To 120 Step 20
Picture1.Line (i, 0)-(i, 1)
Picture1.CurrentX = i - 5: Picture1.CurrentY = -1: Picture1.Print i
Next i
For i = 20 To 80 Step 20
Picture1.Line (0, i)-(2, i)
Picture1.CurrentX = -10: Picture1.CurrentY = i + 1: Picture1.Print "0"; i / 200
Next i
Picture1.DrawWidth = 2
For n = 600 To 4000
For j = 1 To 5
Tq = -19.313 + 295.27 * (n / 1000) - 165.44 * (n / 1000) ^ 2 + 40.874 * (n / 1000) ^ 3 - 3.8445 * (n / 1000) ^ 4
Ft = Tq * i0 * ηt * x(j) / r
Ua = 0.377 * r * n / (i0 * x(j))
Fw = 0.5 * 1.2258 * ca * (Ua / 3.6) ^ 2
D = 200 * (Ft - Fw) / (m * 9.8)
If D > 0 Then
Picture1.Circle (Ua, D), 0.1, vbBlue
End If
Next j
Next n
End Sub
'.............功率平衡图
Private Sub Command4_Click()
Picture1.DrawWidth = 1
Picture1.Cls
Dim i, j, n As Integer
Picture1.Scale (-20, 80)-(150, -10)
Picture1.Line (0, 0)-(125, 0): Picture1.Line (0, 0)-(0, 75)
Picture1.CurrentX = 125: Picture1.CurrentY = -3: Picture1.Print "Ua(km/h)"
Picture1.CurrentX = 125: Picture1.CurrentY = 0.9: Picture1.Print "→"
Picture1.CurrentX = 5: Picture1.CurrentY = 75: Picture1.Print "Pe/kw"
Picture1.CurrentX = -2: Picture1.CurrentY = 75: Picture1.Print "↑"
Picture1.CurrentX = -5: Picture1.CurrentY = -1: Picture1.Print "0"
For i = 20 To 120 Step 20
Picture1.Line (i, 0)-(i, 1)
Picture1.CurrentX = i - 5: Picture1.CurrentY = -1: Picture1.Print i
Next i
For i = 10 To 70 Step 10
Picture1.Line (0, i)-(2, i)
Picture1.CurrentX = -10: Picture1.CurrentY = i + 1: Picture1.Print i
Next i
Picture1.DrawWidth = 2
For n = 600 To 4000
For j = 1 To 5
Tq = -19.313 + 295.27 * (n / 1000) - 165.44 * (n / 1000) ^ 2 + 40.874 * (n / 1000) ^ 3 - 3.8445 * (n / 1000) ^ 4
Pe = Tq * n / 9550
Ua = 0.377 * r * n / (i0 * x(j))
Pf = (m * f * Ua / 3600) / ηt
Pw = (ca * Ua ^ 3 / 76140) / ηt
Picture1.Circle (Ua, Pe), 0.1, vbBlue
Picture1.Circle (Ua, (Pf + Pw)), 0.1, vbRed
Next j
Next n
End Sub
'...........外特性曲线
Private Sub Command5_Click()
Picture1.DrawWidth = 1
Picture1.Cls
Dim i, n As Integer
Dim Pe As Single
Picture1.Scale (-400, 250)-(5000, -20)
Picture1.Line (0, 0)-(4500, 0): Picture1.Line (0, 0)-(0, 200)
Picture1.CurrentX = 4300: Picture1.CurrentY = -3: Picture1.Print "n(r/min)"
Picture1.CurrentX = 4500: Picture1.CurrentY = 3: Picture1.Print "→"
Picture1.CurrentX = 100: Picture1.CurrentY = 200: Picture1.Print "Ttq(N*m)"
Picture1.CurrentX = 100: Picture1.CurrentY = 80: Picture1.Print "Pe(kw)"
Picture1.CurrentX = -60: Picture1.CurrentY = 205: Picture1.Print "↑"
Picture1.CurrentX = -5: Picture1.CurrentY = -1: Picture1.Print "0"
For i = 500 To 4000 Step 500
Picture1.Line (i, 0)-(i, 3)
Picture1.CurrentX = i - 200: Picture1.CurrentY = -1: Picture1.Print i
Next i
For i = 20 To 180 Step 20
Picture1.Line (0, i)-(50, i)
Picture1.CurrentX = -300: Picture1.CurrentY = i + 3: Picture1.Print i
Next i
Picture1.DrawWidth = 2
For n = 600 To 4000
Tq = -19.313 + 295.27 * (n / 1000) - 165.44 * (n / 1000) ^ 2 + 40.874 * (n / 1000) ^ 3 - 3.8445 * (n / 1000) ^ 4
Pe = Tq * n / 9550
Picture1.Circle (n, Pe), 2, vbRed
Picture1.Circle (n, Tq), 2, vbBlue
End Sub
2.数据图表：
驱动力与行驶阻力平衡图：
动力特性图：
爬坡度图：
发动机外特性曲线：
发动机功率平衡图：
3.2 实验记录（核心代码及调试过程）
1.核心代码：
'.........加速度的倒数曲线
Private Sub Command7_Click()
Picture1.DrawWidth = 1
Picture1.Cls
Dim i, j, n As Integer
Dim Ff, Ft, Fw, a, δ As Single
Picture1.Scale (-20, 160)-(150, -10)
Picture1.Line (0, 0)-(125, 0): Picture1.Line (0, 0)-(0, 150)
Picture1.CurrentX = 125: Picture1.CurrentY = -3: Picture1.Print "Ua(km/h)" Picture1.CurrentX = 125: Picture1.CurrentY = 1.8: Picture1.Print "→"
Picture1.CurrentX = 5: Picture1.CurrentY = 150: Picture1.Print "1/a"
Picture1.CurrentX = -2: Picture1.CurrentY = 150: Picture1.Print "↑"
Picture1.CurrentX = -5: Picture1.CurrentY = -1: Picture1.Print "0"
For i = 20 To 120 Step 20
Picture1.CurrentX = i - 5: Picture1.CurrentY = -1: Picture1.Print i
Next i
For i = 10 To 140 Step 10
Picture1.Line (0, i)-(2, i)
Picture1.CurrentX = -10: Picture1.CurrentY = i + 1: Picture1.Print i / 10
Next i
Picture1.DrawWidth = 2
For n = 600 To 4000
For j = 1 To 5
Tq = -19.313 + 295.27 * (n / 1000) - 165.44 * (n / 1000) ^ 2 + 40.874 * (n / 1000) ^ 3 - 3.8445 * (n / 1000) ^ 4
Ft = Tq * i0 * ηt * x(j) / r
Ua = 0.377 * r * n / (i0 * x(j))
Fw = 0.5 * 1.2258 * ca * (Ua / 3.6) ^ 2
Ff = f * 9.8 * m
δ = 1 + Iw / (m * r ^ 2) + If1 * x(j) ^ 2 * i0 ^ 2 * ηt / (m * r ^ 2)
a = (Ft - Ff - Fw) / (δ * m)
If a > 0 Then
Picture1.Circle (Ua, 10 / a), 0.1, vbBlue
End If
Next j
Next n
End Sub
'..............直接挡加速时间曲线
Private Sub Command8_Click()
Picture1.DrawWidth = 1
Picture1.Cls
Dim i, n As Integer
Dim Ff, Ft, Fw, a, δ, t, tt As Single
Picture1.Line (0, 0)-(55, 0): Picture1.Line (0, 0)-(0, 75)
Picture1.CurrentX = 53: Picture1.CurrentY = -3: Picture1.Print "t(s)"
Picture1.CurrentX = 55: Picture1.CurrentY = 0.9: Picture1.Print "→"
Picture1.CurrentX = 2: Picture1.CurrentY = 75: Picture1.Print "Ua(km/h)"
Picture1.CurrentX = -0.7: Picture1.CurrentY = 75: Picture1.Print "↑"
Picture1.CurrentX = -3: Picture1.CurrentY = -1: Picture1.Print "0"
For i = 10 To 50 Step 10
Picture1.Line (i, 0)-(i, 1)
Picture1.CurrentX = i - 1.5: Picture1.CurrentY = -1: Picture1.Print i
Next i
For i = 10 To 70 Step 10
Picture1.Line (0, i)-(1, i)
Picture1.CurrentX = -5: Picture1.CurrentY = i + 1: Picture1.Print i
Next i
t = 0
For Ua = 0 To 70 Step 0.01
n = Ua * i0 * x(4) / 0.377 / r
If n >= 600 Then
Tq = -19.313 + 295.27 * (n / 1000) - 165.44 * (n / 1000) ^ 2 + 40.874 * (n / 1000) ^ 3 - 3.8445 * (n / 1000) ^ 4
Ft = Tq * i0 * ηt * x(4) / r
Fw = 0.5 * 1.2258 * ca * (Ua / 3.6) ^ 2
Ff = f * 9.8 * m
δ = 1 + Iw / (m * r ^ 2) + If1 * x(4) ^ 2 * i0 ^ 2 * ηt / (m * r ^ 2)
a = (Ft - Ff - Fw) / (δ * m)
t = t + 1 * 0.01 / a / 3.6
Picture1.Circle (t, Ua), 0.1, vbBlue
End If
Next Ua
Picture1.CurrentX = 15: Picture1.CurrentY = 80: Picture1.Print "0-70（km/h）直接挡加速时间为："; tt; "s"
End Sub
2.数据图表：
加速度曲线图：
加速度倒数曲线：
直接挡加速时间图：
3.2 实验记录（核心代码及调试过程）
1.核心代码
'.........最高挡百公里耗油量
Private Sub Command1_Click()
Dim n(1 To 8), B0(1 To 8), B1(1 To 8), B2(1 To 8), B3(1 To 8), B4(1 To 8)
n(1) = 815: n(2) = 1207: n(3) = 1614: n(4) = 2012: n(5) = 2603: n(6) = 3006: n(7) = 3403: n(8) = 3804
B0(1) = 1326.8: B0(2) = 1354.7: B0(3) = 1284.4: B0(4) = 1122.9: B0(5) = 1141#: B0(6) = 1051.2: B0(7) = 1233.9: B0(8) = 1129.7
B1(1) = -416.46: B1(2) = -303.98: B1(3) = -189.75: B1(4) = -121.59: B1(5) = -98.893: B1(6) = -73.714: B1(7) = -84.478: B1(8) = -45.291
B2(1) = 72.379: B2(2) = 36.657: B2(3) = 14.524: B2(4) = 7.0035: B2(5) = 4.4763: B2(6) = 2.8593: B2(7) = 2.9788: B2(8) = 0.71113
B3(1) = -5.8629: B3(2) = -2.0553: B3(3) = -0.51184: B3(4) = -0.18517: B3(5) = -0.091077: B3(6) = -0.05138: B3(7) = -0.047449: B3(8) = -0.00075215
B4(6) = 0.00035032: B4(7) = 0.0002823: B4(8) = -0.000038568
Picture1.Scale (-10, 50)-(150, -10)
Picture1.Line (0, 0)-(120, 0): Picture1.Line (0, 0)-(0, 45)
Picture1.CurrentX = 120: Picture1.CurrentY = -3: Picture1.Print "Ua(km/h)" Picture1.CurrentX = 120: Picture1.CurrentY = 1: Picture1.Print "→"
Picture1.CurrentX = 5: Picture1.CurrentY = 45: Picture1.Print "Qs(L/100km)" Picture1.CurrentX = -2.5: Picture1.CurrentY = 46: Picture1.Print "↑"
Picture1.CurrentX = -3: Picture1.CurrentY = -1: Picture1.Print "0"
For i = 20 To 100 Step 20
Picture1.Line (i, 0)-(i, 1)
Picture1.CurrentX = i - 6: Picture1.CurrentY = -1: Picture1.Print i
Next i
For i = 10 To 40 Step 10
Picture1.Line (0, i)-(3, i)
Picture1.CurrentX = -12: Picture1.CurrentY = i + 1: Picture1.Print i
Next i
g = 3880 * 9.8
f = 0.013
CdA = 2.77
ηt = 0.85
i0 = 5.83
r = 0.367
Uamin = 0.377 * r * n(1) / (0.793 * i0)
Qmin = 100
Picture1.DrawWidth = 2
For ua = Uamin To 100 Step 0.01
Pe = (g * f * ua / 3600 + CdA * ua ^ 3 / 76140) / ηt
nn = 0.793 * i0 * ua / (0.377 * r)
For i = 1 To 7
B0B0 = B0(i) + (B0(i + 1) - B0(i)) * (nn - n(i)) / (n(i + 1) - n(i))
B1B1 = B1(i) + (B1(i + 1) - B1(i)) * (nn - n(i)) / (n(i + 1) - n(i))
B2B2 = B2(i) + (B2(i + 1) - B2(i)) * (nn - n(i)) / (n(i + 1) - n(i))
B3B3 = B3(i) + (B3(i + 1) - B3(i)) * (nn - n(i)) / (n(i + 1) - n(i))
B4B4 = B4(i) + (B4(i + 1) - B4(i)) * (nn - n(i)) / (n(i + 1) - n(i))
End If
Next i
Be = B0B0 + B1B1 * Pe + B2B2 * Pe ^ 2 + B3B3 * Pe ^ 3 + B4B4 * Pe ^ 4
Qs = Pe * Be / (1.02 * 7 * ua)
Picture1.PSet (ua, Qs), vbBlue
If Qs < Qmin Then
Qmin = Qs
Uj = ua
End If
Next ua
Text1.Text = CInt(Uj * 100) / 100
Text2.Text = CInt(Qmin * 100) / 100
End Sub
'..............次高挡百公里耗油量
Private Sub Command2_Click()
Picture1.Cls
Picture1.DrawWidth = 1
Dim n(1 To 8), B0(1 To 8), B1(1 To 8), B2(1 To 8), B3(1 To 8), B4(1 To 8)
n(1) = 815: n(2) = 1207: n(3) = 1614: n(4) = 2012: n(5) = 2603: n(6) = 3006: n(7) = 3403: n(8) = 3804
B0(1) = 1326.8: B0(2) = 1354.7: B0(3) = 1284.4: B0(4) = 1122.9: B0(5) = 1141#: B0(6) = 1051.2: B0(7) = 1233.9: B0(8) = 1129.7
B1(1) = -416.46: B1(2) = -303.98: B1(3) = -189.75: B1(4) = -121.59: B1(5) = -98.893: B1(6) = -73.714: B1(7) = -84.478: B1(8) = -45.291
B2(7) = 2.9788: B2(8) = 0.71113
B3(1) = -5.8629: B3(2) = -2.0553: B3(3) = -0.51184: B3(4) = -0.18517: B3(5) = -0.091077: B3(6) = -0.05138: B3(7) = -0.047449: B3(8) = -0.00075215
B4(1) = 0.17768: B4(2) = 0.043072: B4(3) = 0.0068164: B4(4) = 0.0018555: B4(5) = 0.00068906: B4(6) = 0.00035032: B4(7) = 0.0002823: B4(8) = -0.000038568
Picture1.Scale (-10, 50)-(150, -10)
Picture1.Line (0, 0)-(120, 0): Picture1.Line (0, 0)-(0, 45)
Picture1.CurrentX = 120: Picture1.CurrentY = -3: Picture1.Print "Ua(km/h)"
Picture1.CurrentX = 120: Picture1.CurrentY = 1: Picture1.Print "→"
Picture1.CurrentX = 5: Picture1.CurrentY = 45: Picture1.Print "Qs(L/100km)"
Picture1.CurrentX = -2.5: Picture1.CurrentY = 46: Picture1.Print "↑"
Picture1.CurrentX = -3: Picture1.CurrentY = -1: Picture1.Print "0"
For i = 20 To 100 Step 20
Picture1.Line (i, 0)-(i, 1)
Picture1.CurrentX = i - 6: Picture1.CurrentY = -1: Picture1.Print i
Next i
For i = 10 To 40 Step 10
Picture1.Line (0, i)-(3, i)
Picture1.CurrentX = -12: Picture1.CurrentY = i + 1: Picture1.Print i
Next i
Picture1.DrawWidth = 2
m = 3880
g = 3880 * 9.8
f = 0.013
CdA = 2.77
ηt = 0.85
i0 = 5.83
r = 0.367
Qmin = 100
Pe = (g * f * ua / 3600 + CdA * ua ^ 3 / 76140) / ηt
nn = 0.793 * i0 * ua / (0.377 * r)
For i = 1 To 7
If nn >= n(i) And nn <= n(i + 1) Then
B0B0 = B0(i) + (B0(i + 1) - B0(i)) * (nn - n(i)) / (n(i + 1) - n(i)) B1B1 = B1(i) + (B1(i + 1) - B1(i)) * (nn - n(i)) / (n(i + 1) - n(i)) B2B2 = B2(i) + (B2(i + 1) - B2(i)) * (nn - n(i)) / (n(i + 1) - n(i)) B3B3 = B3(i) + (B3(i + 1) - B3(i)) * (nn - n(i)) / (n(i + 1) - n(i)) B4B4 = B4(i) + (B4(i + 1) - B4(i)) * (nn - n(i)) / (n(i + 1) - n(i)) End If
Next i
Be = B0B0 + B1B1 * Pe + B2B2 * Pe ^ 2 + B3B3 * Pe ^ 3 + B4B4 * Pe ^ 4
Qs = Pe * Be / (1.02 * 7 * ua)
Picture1.PSet (ua, Qs), vbRed
If Qs < Qmin Then
Qmin = Qs
Uj = ua
End If
Next ua
Text1.Text = CInt(Uj * 100) / 100
Text2.Text = CInt(Qmin * 100) / 100
End Sub
2.数据图表：
最高档等速百公里燃油消耗量曲线图：
次高档等速百公里燃油消耗量曲线图：
整车电量平衡计算
谈汽车电平衡的设计计算及验证方法 随着汽车电子电器技术的迅速发展，电器功能日益增多且复杂，对车辆舒适、智能和安全可靠性等要求的提高，整车电平衡的设计及验证尤显重要。整车电平衡是指发电机、蓄电池、整车用电器在一定时间内的电能产生与消耗达到稳定的一种平衡状态，是重要的整车性能指标。它体现了发电机的输出能力与整车用电需求的匹配关系，而不同的整车性能目标定义，对整车电平衡的性能要求也是不同的，所以需要有合适的汽车电平衡设计计算和验证方法。 本文主要结合试验数据，分析改进电平衡的设计计算方法；重点结合整车电平衡试验做出动态特性曲线，对电平衡理论计算结果进行验证。 1 汽车电平衡的设计方法 汽车电平衡的设计需要考虑发动机参数、整车用电器功率和使用频度等，图1为电平衡设计示意图，描述了电平衡关键零部件选型顺序和各关键零部件的影响因素。 2 关键零部件的计算选型 2.1起动机的选型 起动机的作用是起动发动机，一般需要起动机以大电流工作2～5s。发动机的起动特性决定了起动机的性能参数，发动机的起动特性参数包括起动转矩和起动转速。设定试验测定极限低温工况下的起动转矩为M0，起动转速为n0，由M0和n0可得出起动需求功率P0=M0×n0×2π／60。 根据传动比i和齿轮的啮合效率η（η通常为0.9），可计算出发动机起动过程中起动机的输出参数：转矩M1=M0／i，转速n1=n0×i，功率P1=P0／η。 起动机的输出功率会随温度而变化，再根据起动机温度系数修正出常温下起动机输出的转矩和功率，即可完成起动机的参数选择。
蓄电池最主要的作用是起动发动机，故其选型应先分析起动机（或发动机）的特性。蓄电池的低温起动电流应大于起动机输出特性曲线图上功率最大点对应的起动电流，以确保实现起动发动机，同时小于功率曲线与力矩曲线交点处对应的电流，在符合条件的蓄电池中选择容量较大者以增加起动发动机的可靠性。依此原则选择的蓄电池，不会因蓄电池容量选择过大出现浪费及蓄电池体积增大而影响整车 的装配空间及质量。 车辆在长途运输或长时停放后应能起动发动机，所以在蓄电池选型时，需考虑整车静态电流的验证。整车静态电流计算公式为I静=C20×（90%-65%-1‰×T）／（T×24）（1）式中：I静———整车静态电流；90%———下线时，蓄电池的实际容量与额定容量的百分比；65%———确保车辆正常起动的蓄电池最低实际电量与额定电量的百分比；1‰———蓄电池1天的自损耗率；T———储运时间；C20———蓄电池的20h率额定容量，Ah。最后，根据蓄电池的布置位置、车辆销售区域及主要用途等，微调蓄电池的参数。以奇瑞公司某在研车型M为例，根据发动机起动转矩和起动转速选择了1.3kW起动机。该起动机输出特性曲线如图2所示。 根据蓄电池选型方法，结合图2，选择蓄电池放电电流应为260～500A，符合条件的蓄电池容量为45Ah（冷起动电流为425A）和60Ah（冷起动电流为480A），可初选蓄电池的容量为60Ah。根据式（1）可知，若储运时间要求为45天，蓄电池容量为60Ah，得：I静=11.4mA，故整车静态电流须小于11.4mA。
汽车理论图形MATLAB程序
功率平衡图 m=1230;g=9.8; ig=[3.615 2.053 1.393 1.031 0.837]; i0=3.75; r=0.31;yt=0.9;f=0.017;CD=0.31;A=2.2; np=6000;Pemax=83; %绘制汽车驱动力与行驶阻力平衡图 for i=1:56; n=500:100:6000; Pe(i)=Pemax*(n(i)/np+(n(i)/np)^2-(n(i)/np)^3); Tq(i)=9549*Pe(i)/n(i); end for j=1:5 for i=1:56 Ft(i,j)=Tq(i)*ig(j)*i0*yt/r; ua(i,j)=0.377*r*n(i)/(ig(j)*i0); Fz(i,j)=m*g*f+CD*A*(ua(i,j)^2)/21.15; end end figure plot(ua,Ft,ua,Fz); title('汽车驱动力与行驶阻力平衡图'); xlabel('ua(km/h)'); ylabel('Ft(N)'); text(20,6700,'Ft1'); text(40,4000,'Ft2'); text(50,2800,'Ft3'); text(80,2000,'Ft4'); text(100,1600, 'Ft5'); text(100,800,'Ff+Fw'); for k=1:56; n=500:100:6000; Pe(k)=Pemax*(n(k)/np+(n(k)/np)^2-(n(k)/np)^3); Tq(k)=9549*Pe(k)/n(k); Ft(k)=Tq(k)*ig(4)*i0*yt/r; ua(k)=0.377*r*n(k)/(ig(4)*i0); Fz(k)=m*g*f+CD*A*(ua(k)^2)/21.15; E(k)=abs((Ft(k)-Fz(k))); end [Emin,kmin]=min(E); Umax=ua(kmin)
用matlab绘制汽车驱动力 行驶阻力平衡图
汽车驱动力-行驶阻力平衡图m=3880; g=9.8; nmin=600;nmax=4000; G=m*g; ig=[5.56 2.769 1.644 1.00 0.793]; nT=0.85; r=0.367; f=0.013; CDA=2.77; i0=5.83; L=3.2; a=1.947; hg=0.9; If=0.218;Iw1=1.798;Iw2=3.598; n=600:10:4000; Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1 000).^4; Ft1=Tq*ig(1)*i0*nT/r; Ft2=Tq*ig(2)*i0*nT/r; Ft3=Tq*ig(3)*i0*nT/r; Ft4=Tq*ig(4)*i0*nT/r; Ft5=Tq*ig(5)*i0*nT/r; ua1=0.377*r*n/ig(1)/i0; ua2=0.377*r*n/ig(2)/i0; ua3=0.377*r*n/ig(3)/i0; ua4=0.377*r*n/ig(4)/i0; ua5=0.377*r*n/ig(5)/i0; ua=0:5:120; Ff=G*f; Fw=CDA*ua.^2/21.15; Fz=Ff+Fw; plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua,Fz); title('汽车驱动力-行驶阻力平衡图'); xlabel('ua(km*h^-1)'); ylabel('Ft/N'); gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ft5'),gtext('Ff+F w');
汽车功率半导体行业深度报告
汽车功率半导体行业深度报告
报告综述： 汽车功率半导体5年近7倍空间，IGBT最受益 政策支持、节能减排双重驱动，新能源汽车加速渗透，预计 2025 年国内新能源汽车渗透率将达到 20%，2030 年欧盟新能源汽车渗透率将达到 40%。汽车电动化趋势下车用功率半导体单车价值大幅提升。据英飞凌统计，功率半导体 ASP 将从传统燃油车的 71 美元大幅提升至全插混/纯电汽车的 330 美元，是传统燃油车的 4.6 倍。根据我们的测算，预计2025 年全球汽车功率半导体市场规模将达到 80 亿美元，2025 年全球新能源车用功率半导体市场规模将达到 53 亿美元，是 2020 年的 7.3 倍，年复合增速高达 48.8%，未来十年中美欧三地区新能源汽车充电桩用 IGBT 市场将有 94 亿美元增量空间。目前车用功率半导体中主要用到的是 IGBT 和 MOSFET，而 IGBT 在新能源车中是电驱系统主逆变器的核心器件，并可用于辅逆变电路、DC/DC 直流斩波电路、OBC（充电/逆变）等，单车价值达到 273 美元，占车用功率半导体 ASP 的 83%，是绝对大头。我们预计 2025 年全球新能源汽车 IGBT 市场规模将达到 44 亿美元，年复合增速约 48.8%，是电动化趋势下的汽车功率半导体中最受益品种。 产品、工艺、先发优势三大壁垒构筑强护城河 1）产品壁垒：车规级 IGBT 需具备使用寿命长、故障率低、抗震性高等严格要求，能适应“极热”“极冷”的高低温工况、粉尘、盐碱等恶劣的工况环境，承受频繁启停带来的电流频繁变化，对产品要求极高。2）工艺壁垒：车规级 IGBT 设计时需保证开通关断、抗短路和导通压降三者的平衡，参数优化特殊复杂。生产制造时薄片工艺容易碎裂、正面金属熔点限制导致退火温度控制难度大。此外，IGBT 模块封装的焊接和键合环节技术要求同样较高。3）认证周期长、替换成本高、具备经验曲线效应，行业先发优势明显。 a）车规级IGBT 需满足可靠性标准、质量管理标准、功能安全标准，才有资格进入一级汽车厂商的供应链，认证周期一般至少 2 年。b）由于 IGBT 模块是汽车中的关键部件，下游厂商出于安全性、可靠性的考虑，替换时往往呈谨慎态度，只有经过大量验证测试并通过综合评定后，才会做出大批量采购决策，替换成本高。c）IGBT 业务需要长期的经验积累才能达到良好的 know-how 水平。d）IGBT 行业属于资本密集型行业，生产、测试设备基本需要进口。此外，对 IGBT 生产企业的流动资金需求量也较大，新进入者在前期往往面临投入大、产出少的情况，需要较强的资金实力作后盾，才能持续进行产品的研发、生产和销售。综合来看，IGBT 行业中的先行企业具有明显的先发优势。 竞争格局优成为成长行业“优质赛道”，但当前国产化率仍然较低 据 Omdia 2019 年统计数据，全球 IGBT 模块前十大厂商占据了 76%份额，市场份额集
用matlab绘制汽车驱动力-行驶阻力平衡图
汽车驱动力-行驶阻力平衡图 m=3880; g=9.8; nmin=600;nmax=4000; G=m*g; ig=[5.56 2.769 1.644 1.00 0.793]; nT=0.85; r=0.367; f=0.013; CDA=2.77; i0=5.83; L=3.2; a=1.947; hg=0.9; If=0.218;Iw1=1.798;Iw2=3.598; n=600:10:4000; Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445 *(n/1000).^4; Ft1=Tq*ig(1)*i0*nT/r; Ft2=Tq*ig(2)*i0*nT/r; Ft3=Tq*ig(3)*i0*nT/r; Ft4=Tq*ig(4)*i0*nT/r; Ft5=Tq*ig(5)*i0*nT/r; ua1=0.377*r*n/ig(1)/i0; ua2=0.377*r*n/ig(2)/i0; ua3=0.377*r*n/ig(3)/i0; ua4=0.377*r*n/ig(4)/i0; ua5=0.377*r*n/ig(5)/i0; ua=0:5:120; Ff=G*f; Fw=CDA*ua.^2/21.15; Fz=Ff+Fw; plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua,Fz); title('汽车驱动力-行驶阻力平衡图'); xlabel('ua(km*h^-1)'); ylabel('Ft/N'); gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ft5'),gtext( 'Ff+Fw');
汽车理论课后题matlab程序
汽车理论课后题matlab程序
————————————————————————————————作者：————————————————————————————————日期: ?
1.3 n＝600：1：４000; r=0．367; i０=５.83; eff＝0.85; f＝0.０13； ｍ＝3880；g=9.8; G=m*g; CdA=2.77；ａ=１.94７; hg=０.9；L＝3.2； Iw1=1．79８; Iw２＝3.598; Ｉｗ=Iw１+Iw２; Iｆ=0.218； Tｔｑ＝-19.313+295.２7*n/1000－1６5.44*（ｎ／10０0).^2+40.874*(n/1０00).^3－３.８445*（ｎ/１00０)．^4； %驱动力行驶阻力平衡图 for ｉｇ=[5.56,2.７69,１.644,１.０0,0.793] Uａ=0.377*r＊n/ig/ｉ0； Ft=Tｔq*ｉg＊i０*eff/r； ｐlｏt(Ua,Ft)； ｈold on; ｅｎｄ Ff=G*f; ua=0：０.1:max（Uａ); Ｆｗ＝CdＡ*ｕa．^２/２1.15; plot(ua，(Ff+Ｆw)); title('驱动力-行驶阻力平衡图'）; xlａbel(＇Ua／(ｋm/ｈ）'）;ylab ｅl(＇Fｔ/Ｎ'）; gtext('Ｆt1')，gｔext('Fｔ２'),ｇｔext（'Ft3')，gｔexｔ('Ft4＇），ｇteｘt('Fｔ5')，gtext(＇Fｆ+Fw') [x，y]=ginput(1); disp('汽车的最高车速')；disp(x);dｉsp(＇km/h'); %最大爬坡度及最大爬坡度时的附着率 Ua=0.37７*r*ｎ／5.56/ｉ０； Ft=Ｔtq*５.56*i0*eff/r; Ｆw=CｄA*Uａ.^２/21.15； i=tａｎ（asiｎ((Ft-(Ｆf+Fｗ))/Ｇ)); disp（'汽车的最大爬坡度');diｓp（max(ｉ）); Ｃ=maｘ(i)/(a/L＋hg/L＊maｘ(i))； diｓp（＇克服最大爬坡度时的附着率'）;disｐ(C); %加速度倒数曲线 fｉguｒe; ｆoｒ ig=[5.56,２.７６９,1.644，1．00，0．７93] Ua=0．３７７*ｒ*ｎ/ig／i0; q=1+Iｗ/(m*ｒ＾2)+If＊ig^2*i０^２＊eff/(m＊r^2); Ft=Tｔq*ig*i0*ｅfｆ／r; Ｆw=CdA*Ua.^2／２1.15; ａs=(Ｆt-（Fｆ+Ｆw）)/q/ｍ; ｐlｏt(Ua,１．/as)； hold on; ｅnd axiｓ([0 98 0 １0]); titlｅ(＇行驶加速度倒数曲线')；ｘlabel(＇Ua/(kｍ/h）');ylａｂel（'１/a')； gteｘｔ('１/a1'),gｔｅxt('１/a2'),ｇtext('1／ａ3＇),ｇtext('１/a4')，ｇtext('1/
汽车理论习题
汽车理论习题： 1.3 1）绘制汽车驱动力与行驶阻力平衡图。 2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。 3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行驶至70km/h 的车速－时间曲线，或者用计算机求汽车用2档起步加速行驶至70km/h 的加速时间。 分析：本题主要考察知识点为汽车驱动力－行使阻力平衡图的应用和附着率的计算、等效坡度的概念。 应明确道路的概念：坡度的定义tan i α =。求最大爬坡度时可以对行使方程进行适当简化，可以简 化的内容包括两项cos 1α ≈和sin tan αα ≈，简化的前提是道路坡度角不大，当坡度角较大时简化带 来的误差会增大。计算时，要说明做了怎样的简化并对简化的合理性进行评估。 2）求最大爬坡度及克服该坡度时相应的附着率 ②求最大爬坡度： 0=dt du ，()w f t i F F F F +-= 一般汽车最大爬坡度约为30％左右，所以利用行驶方程确定最大爬坡度时应以G sin α作为坡度阻力，即 ()sin t f w G F F F α=-+，tan i α = ()??? ? ?? +-=G F F F i w f t arcsin tan 汽车最大爬坡度max i 为Ⅰ档时的最大爬坡度。利用编程计算可得，352.0max =i 。 ③如是前轮驱动，1?C ＝ q b hg q L L -；相应的附着率1?C 为1.20，不合理，舍去。 如是后轮驱动，2?C ＝ q a hg q L L +；相应的附着率2?C 为0.50。 3）绘制汽车行驶加速度倒数曲线，求加速时间 因各档加速度倒数曲线并无交点，所以从二档开始直接积分即可。 常见错误是未将车速的单位进行换算，时间大3.6倍。 1.7确定上述F.F 轿车在φ＝0.2及0.7路面上的附着力，并求由附着力所决定的极限最高车速与极限最大爬坡度及极限最大加速度（在求最大爬坡度和最大加速度时可设Fw=0）。 （1）求极限最高车速的求解可根据汽车行驶方程得到。 2 W max F 21.15 D a C A u = 具体方法有两个，一是根据地面作用于驱动轮的地面切向反作用力的表达式（1－15），由附着系数得到最大附着力，滚动阻力已知，即可求得最高车速时的空气阻力和最高车速。 X1W f2i j'F =F +F F F ++ （1－15） X11z1F =F F ??=
第2讲 1-3汽车驱动力-行驶阻力平衡图与动力特性图 1-4汽车行驶的附着条件与汽车的附着率
第2讲 2学时 教学目的及要求： 掌握汽车的驱动力－行驶阻力平衡图的绘制，汽车加速度曲线的绘制，汽车的加速度倒数曲线的绘制，汽车的加速时间曲线的绘制，汽车爬坡度曲线的绘制，汽车动力特性图的绘制。 主要内容： §1-3汽车驱动力-行驶阻力平衡图与动力特性图 §1-4汽车行驶的附着条件与汽车的附着率 教学重点： 驱动力－行驶阻力平衡图及利用驱动力－行驶阻力平衡图分析汽车的动力性，汽车的动力因数，动力特性图及利用动力特性图分析汽车的动力性 教学难点： 动力特性图及利用动力特性图分析汽车的动力性 教学过程: §1—3 汽车行驶的驱动与附着条件 一、汽车行驶的驱动与附着条件： 1、驱动条件—首先得有劲 δm du dt = F t – (F f + F W + F i ) ≥ 0 F t ≥F f +F W +F i 2、附着条件—有劲还得使得上 用F φ表示轮胎切向反力的极限，在硬路面上它与驱动轮所受的法向反力成正比：（φ为附着系数） （1）驱动轮的附着力： 前轮驱动汽车： F φ1 = F Z1φ 后轮驱动汽车： F φ2 = F Z2φ 全轮驱动汽车： F φ1 = F Z1φ F φ2 = F Z2φ （2）汽车的附着力： 前轮驱动汽车： F φ = F Z1φ 后轮驱动汽车： F φ = F Z2φ 全轮驱动汽车： F φ = F Z φ = F Z1φ+F Z1φ 对前驱动轮 F x1 ≤ F Z1φ 前驱动轮的附着率： C φ1 = F X1F Z1 则要求 C φ1 ≤φ
对后驱动轮F x2 ≤ F Z2 φ 后驱动轮的附着率： C φ2 =F X2 F Z2 则要求C φ2 ≤φ ∴F t≤F Z2（f+φ）∵f<<φ∴F t≤F Z2φ 一般形式F t ≤F Z φφ 3、驱动与附着条件： F f+F W+F i≤F t≤F Zφφ 二、汽车的附着力：F φ 1、汽车附着力——在车轮与路面没有相对滑动的情况下，路面对车轮提供的切向 反力的极限值。 Fφ=F Zφφ Fφ取决于： ①在硬路面上——可以是最大的静摩擦力， 主要取决于路面与轮胎的性质； ②在软路面上——取决于土壤的剪切强度和车轮与土壤的结合强度 2、Fφ的影响因素： ⑴载重量： 增加驱动轮的法向反力X2，有利于驱动。 例：越野车由货车的F Z2↗（F Z2+F Z1），使Fφ↗ ⑵轮胎结构： 深大花纹——在松软路面上，使土壤与车轮的结合强度提高； 松软路上放气P↘——胎面接地面积大，嵌入土壤的花纹数多，抓地能力强，且沉陷量小，土壤阻力小； ⑶附着系数：φ 取决于路面种类与状况、轮胎结构（花纹、材料等）及u a等因素。 三、驱动轮的法向反作用力 ——汽车行驶时重量再分配 1、根据受力图列方程： 将作用在汽车上的各力对前、后轮接地面中心取矩，则得： F Z1 = G L(bcosα- h g sinα)– 1 L(mh g du dt+∑T j)- F ZW1- 1 L∑T f F Z2 = G L(acosα+ h g sinα)+ 1 L(mh g du dt+∑T j)- F ZW2+ 1 L∑T f 式中，∑T j = T jW1+T jW2 ,∑T f = T f1+T f2 忽略旋转质量的惯性阻力偶矩和滚动阻力偶矩： F Z1 = F ZS1–mh g L du dt- F ZW1 F Z2 = F ZS2 + mh g L du dt- F ZW2 作用在驱动轮上的地面切向反作用力： 前轮驱动：F X1 = F f2 + F W + F i + m du dt 后轮驱动：F X2 = F f1 + F W + F i + m du dt
汽车理论汽车设计课程设计说明书
湖北汽车工业学院 Hubei Automotive Industries Institute 课程设计说明书 课程名称汽车理论 设计题目汽车动力性 班号专业车辆工程学号 学生姓名 指导教师(签字) 起止日期 2011 年 7 月 4 日—— 2011 年 7 月 9 日
目录 1.设计任务及要求．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．．1 2.车辆参数 (2) 3.汽车动力性能计算............. ..... ................... 3.1驱动力-行驶阻力平衡图...................... 3.2最高转速Ｕamax....................... 3.3加速时间t............................... 3.4汽车加速度倒数图............................... 3.5汽车加速时间图............................... 3.6汽车爬坡度图............................... 3.7汽车动力特性图．............................... 3.8汽车功率平衡图.................................. 4.GUI界面设计........................................ 5.归纳与总结........................................ 6.参考文献......................................
1.5 汽车的功率平衡
1.5 汽车的功率平衡 汽车行驶时，不仅存在驱动力与行驶阻力的平衡关系，而且也存在发动机功率和汽车行驶的阻力功率间的平衡关系。即发动机发出的有效功率，始终等于机械传动损失与全部运动阻力所消耗的功率。 1.5.1 功率平衡方程 汽车运动阻力所消耗的功率，有滚动阻力功率、空气阻力功率、坡度阻力功率及加速阻力功率，它们的表达式为 f P W P i P j P 3600cos 10006.3a a f f u Gf u F P α=×= 3600 sin 3600a a i i u G u F P α== 7614036003a D a W W Au C u F P == dt du g Gu u F P a a j j 36003600δ== 功率平衡方程为 )(11 j i W f T T e P P P P P P +++==∑ηη 即 2(cos sin 360021.15a D e T u C Au P Gf G m dt ααδη=+++)a du 当α较小时，i ≈αsin ，1cos ≈α，上式可写成 )15.21(36002dt du g G Au C Gi Gf u P a D T a e δη+++= 1.5.2 功率平衡图及其应用 汽车的功率平衡关系也可以用图解法表示。以纵坐标表示 功率，横坐标表示车速，将发动机功率、汽车经常遇到的 阻力功率e P )(1 W f T P P +η，对应于车速的关系曲线绘在坐标图 上，即得到汽车功率平衡图。图1-20是一紧凑型国产轿车的 功率平衡图。 e P a u ?，可根据发动 发动机功率与行驶车速的关系曲线
机外特性及公式0 0.377a g nr u i i =将发动机转速转换成车速绘得。可见在不同挡位时，功率的大小不变，只是各挡发动机功率曲线所对应的车速位置不同，且低挡时车速低，所占速度变化区域窄；高挡时车速高，所占变化区域宽。 f P 在低速范围内为一直线，在高速时由于滚动阻力系数f 随车速而增大，且比更快的速率加大；而 则是的三次函数。两者叠加后，阻力功率曲线是一条斜率越来越大的曲线。它与挡位无关，只与车速有关，所以高速时，汽车主要克服空气阻力功率。 a u a u W a P max a P u 图1-20中发动机功率曲线(V 挡)与阻力功率曲线相交点处对应的车速便是在良好水平路面上汽车的最高车速。该轿车的V 挡是经济挡位，其发动机最大功率相对应的车速大于u ，所以用该挡行驶时发动机负荷率高，燃油消耗量低。 max a u u 汽车达最高车速时，0j =，0=i ，则 )(1 W f T e P P P +=η 当汽车在良好水平路而上以的速度等速行驶时，汽车的阻力功率为 a u ′1()f w T P P b ηc += 此时，驾驶员给出某一节气门的开度，发动机功率曲线如图中虚线所示，以维持汽车等速行驶。 但是发动机在汽车行驶速度为u a ′时能发出的功率为e P ac = (图1-20)，于是 1()e f w T P P P ac bc ηab ? +=?= 可用来加速或爬坡。 我们称1 (e f T )w P η?+P P 为汽车的后备功率。 就是说，在一般情况下维持汽车等速行驶所需的发动机功率并不 大，发动机节气门开度较小。当需要爬坡或加速时，驾驶 员加大节气门开度，使汽车的全部或部分后备功率发挥作
(完整版)汽车理论习题Matlab程序
1.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选 其中的一种进行整车性能计算）： 1）绘制汽车驱动力与行驶阻力平衡图。 2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。 3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行驶至70km/h 的车速－时间曲线，或者用计算机求汽车用2档起步加速行驶至70km/h 的加速时间。 轻型货车的有关数据： 汽油发动机使用外特性的Tq-n 曲线的拟合公式为 234 19.313295.27()165.44()40.874() 3.8445()1000100010001000 q n n n n T =-+-+- 式中，Tq 为发动机转矩（N?m ）;n 为发动机转速（r/min ）。 发动机的最低转速n min =600r/min,最高转速n max =4000r/min 。 装载质量 2000kg 整车整备质量 1800kg 总质量 3880kg 车轮半径 0.367m 传动系机械效率 ηt =0.85 滚动阻力系数 f =0.013 空气阻力系数×迎风面积 C D A =2.77m 2 主减速器传动比 i 0=5.83 飞轮转动惯量 I f =0.218kg?m 2 二前轮转动惯量 I w1=1.798kg?m 2 四后轮转动惯量 I w2=3.598kg?m 2 质心至前轴距离（满载） a=1.974m 质心高（满载） hg=0.9m 解：Matlab 程序： (1) 求汽车驱动力与行驶阻力平衡图和汽车最高车速程序： n=[600:10:4000]; Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; m=3880;g=9.8;nmin=600;nmax=4000; G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83; L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598; Ft1=Tq*ig(1)*i0*nT/r; Ft2=Tq*ig(2)*i0*nT/r; Ft3=Tq*ig(3)*i0*nT/r;
汽车理论课后习题Matlab程序详解
陕西理工学院 杨飞华 1 / 34 1.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选 其中的一种进行整车性能计算）： 1）绘制汽车驱动力与行驶阻力平衡图。 2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。 3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行 驶至70km/h的车速－时间曲线，或者用计算机求汽车用2档起步加速行驶至 70km/h的加速时间。 轻型货车的有关数据： 汽油发动机使用外特性的Tq-n曲线的拟合公式为23419.313295.27()165.44()40.874()3.8445() 1000100010001000qnnnn T=?+?+? 式中， Tq为发动机转矩（N?m）;n为发动机转速（r/min）。发动机的最低转速nmin=600r/min,最高转速nmax=4000r/min。装载质量 2000kg 整车整备质量1800kg 总质量3880kg 车轮半径0.367m 传动系机械效率εt=0.85 滚动阻力系数 f=0.013 空气阻力系数×迎风面积CDA=2.77m2 主减速器传动比i0=5.83 飞轮转动惯量 If=0.218kg?m2 二前轮转动惯量Iw1=1.798kg?m2 四后轮转动惯量 Iw2=3.598kg?m2 变速器传动比 ig(数据如下表) Ⅰ档Ⅱ档Ⅲ档Ⅳ档Ⅴ档四档变速器 6.09 3.09 1.71 1.00 - 五档变速器 5.56 2.769 1.644 1.00 0.793 轴距 L=3.2m 质心至前轴距
离（满载）a=1.974m 质心高（满载）hg=0.9m 解：Matlab程序：(1) 求汽车驱动力与行驶阻力平衡图和汽车最高车速程序：n=[600:10:4000]; Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10 00).^4; m=3880;g=9.8;nmin=600;nmax=4000; G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83; L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598; Ft1=Tq*ig(1)*i0*nT/r; Ft2=Tq*ig(2)*i0*nT/r; Ft3=Tq*ig(3)*i0*nT/r; 2 / 34 Ft4=Tq*ig(4)*i0*nT/r; Ft5=Tq*ig(5)*i0*nT/r; ua1=0.377*r*n/ig(1)/i0; ua2=0.377*r*n/ig(2)/i0; ua3=0.377*r*n/ig(3)/i0; ua4=0.377*r*n/ig(4)/i0; ua5=0.377*r*n/ig(5)/i0; ua=[0:5:120]; Ff=G*f;
汽车动力性计算matlab程序
%% 汽车动力性计算（自己编的动力性计算程序，供大家计算动力性时参考，具体参数大家根据所 给程序对应输入，并对坐标轴数值按需要进行修改） clc; clear; close all; %%根据所给发动机数据拟合外特性曲线（发动机数据按照你所得到的数据进行输入） n_test=[500 600 700 800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 2000 2100 2200]; T_test=[975 1108 1298 1496 1546 1620 1670 1785 1974 1974 1970 1889 1829 1748 1669 1700 1524 1105]; figure(1) plot(n_test,T_test,'g'); hold on grid on %p=polyfit(n_test,T_test,7); p=polyfit(n_test,T_test,2); n=[450:1:2200]; Ttq=polyval(p,n); plot(n,Ttq,'k'); xlabel('发动机转速n(r/min)'); ylabel('发动机转矩Ttq(N*m)'); title('发动机转矩曲线'); legend('测试曲线','拟合曲线'); %%所给车型动力总成相关参数 ig=[3.07 2.16 1.48 1.0 0.82]; i0=4.0; eta=0.78; r=0.57; M=25000; g=9.8; c=1.5; f0=0.01; f1=0.0002; f4=0.0005; CD=1; A=8;Iw=3.6;If=0.04; %% 发动机外特性曲线图 figure(2) hold on grid on for i=length(n); Pe=Ttq.*n/9550; end [AX,H1,H2]=plotyy(n,Ttq,n,Pe); xlabel('发动机转速n(r/min)'); ylabel('发动机转矩Ttq(N*m)'); ylabel(AX(2),'发动机功率Pe(Kw)'); title('发动机外特性曲线'); %% 各挡位速度曲线 %计算各挡位车速
汽车功率解释
发动机功率 发动机功率是指发动机做功的快慢。 发动机单位时间内所做的功叫做这发动机的功率。 符号为：P 常用单位为：w（瓦特） 发动机的功率并不等于车子的功率，在机械传动中，功率会有中间损失。另外，出于安全考虑而把车身加厚加重的设计，也会导致发动机功率的损失。所以，装有小排量发动机的车子并不一定就比装载大功率发动机的车子慢或者性能差。 比功率和比扭矩当然是越大越好，说明动力储备充足，但是车子如果加了必须的、豪华的配置而使其数值降低，并不能说明更多的问题。 在同等排量的情况下，功率当然是越大越好，有些车厂采用先进的技术，比如进气门行程控制、EGR、5气门等，功率的提升是必然的，但同时，这意味着发动机成本的增加。厂家为某个车型选配的发动机，其实是在动力性、经济性、排放和成本间找一个平衡点，而不仅仅是更大的功率和更高的扭矩，有时为了成本，甚至可以牺牲掉拥有更大功率的发动机。 发动机的功率大了，燃油消耗肯定会高，排放污染也会增加，对于消费者和环境而言，盲目追求马力或者功率是片面的。国内和国外的一些车型之所以能够为消费者提供很大的价格区间，其实是厂家做了各种各样的发动机、变速箱、动力系统和电器系统的匹配才达到的。但装配了好的发动机（功率大，而并不一定是排量大），车的价格就会明显上升。 扭矩和功率一样，是汽车发动机的主要指数之一，它反映在汽车性能上，包括加速度、爬坡能力以及悬挂等。它的准确定义是：活塞在汽缸里的往复运动，往复一次做有一定的功，它的单位是牛顿。在每个单位距离所做的功就是扭矩了。是这样的，扭矩是衡量一个汽车发动机好坏的重要标准，一辆车扭矩的大小与发动机的功率成正比。举个通俗的例子，比如，像人的身体在运动时一样，功率就像是身体的耐久度，而扭矩是身体的爆发力。对于家用轿车而言，扭矩越大加速性越好；对于越野车，扭矩越大其爬坡度越大；对于货车而言，扭矩越大车拉的重量越大。在排量相同的情况下，扭矩越大说明发动机越好。在开车的时候就会感觉车子随心所欲，想加速就可加速，“贴背感”很好。现在评价一款车有一个重要数据，就是该车在0－100公里／小时的加速时间。而这个加速时间就取决于汽车发动机的扭矩。一般来讲，扭矩的最高指数在汽车2000－4000／分的转速下能够达到，就说明这款车的发动机工艺较好，力量也好。有些汽车在5000／分的转速左右才达到该车扭矩的最高指数，这说明“力量”就不是此车所长。 首先谈谈马力。我们知道，马力表示的是功率，也就是单位时间内做功的多少。如果用N表示功率，N=W/t，其中W是功，t是时间。而功可以用力与位移的乘积来表示。就是我们熟悉的物理公式：W=FScosa 其中F表示力，S表示位移，而a是力的方向与位移方向之间的夹角。在研究汽车发动机活塞做功时，主要讨论F和S对W的影响。 在发动机中，F是与活塞面积和气缸内压强相关的。F=PA，P是气缸内的压强，燃料燃烧产生的压强越大，F就越大；A是活塞顶的面积，我们都知道计算圆面积的公式：A=3.14159D2/4（其中D是缸径）。因此，缸径越大，A就越大，从而F就越大。而S是活塞的行程，也就是活塞往复运动的位移，等于从上止点运动到下止点所经过的距离。S越大，W就越大。 功率是单位时间内做功的多少，所以还要考虑做功的快慢。同样缸径和行程的发动机，如果转速不同，马力也是不同的。转速快的马力大。 就是说，衡量发动机马力的主要因素有四个、气缸压强、活塞截面积、活塞行程和发动机转速。其中活塞截面积与行程与我们熟悉的名词"排量"直接相关。活塞截面积与行程的乘积AS是一个缸的排量，如果发动机有n个相同的气缸，它的排量就是nAS。 扭矩表示使物体加速转动的能力。一般来讲，汽车发动机的马力大通常扭矩也大。如果用扭矩来表示
用matlab绘制汽车功率平衡图
用matlab绘制汽车功率平衡图 n=600:10:4000; Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445 *(n/1000).^4; m=3880;g=9.8; nmin=600;nmax=4000; G=m*g; ig=[5.56 2.769 1.644 1.00 0.793]; nT=0.85;r=0.367; f=0.013;CDA=2.77;i0=5.83; L=3.2;a=1.947;hg=0.9; If=0.218;Iw1=1.798;Iw2=3.598; ua1=0.377*r*n/ig(1)/i0; ua2=0.377*r*n/ig(2)/i0; ua3=0.377*r*n/ig(3)/i0; ua4=0.377*r*n/ig(4)/i0; ua5=0.377*r*n/ig(5)/i0; Ft1=Tq*ig(1)*i0*nT/r; Ft2=Tq*ig(2)*i0*nT/r; Ft3=Tq*ig(3)*i0*nT/r; Ft4=Tq*ig(4)*i0*nT/r; Ft5=Tq*ig(5)*i0*nT/r; Pe1=Ft1.*ua1/3600; Pe2=Ft2.*ua2/3600; Pe3=Ft3.*ua3/3600; Pe4=Ft4.*ua4/3600; Pe5=Ft5.*ua5/3600; ua=0:5:120; Pf=G*f.*ua/3600; Pw=CDA.*ua.^3/76140; Pz=(Pf+Pw)/nT; Pe=max(Pe1); plot(ua1,Pe1,ua2,Pe2,ua3,Pe3,ua4,Pe4,ua5,Pe5,ua,Pz,ua,Pe); title('汽车功率平衡图'); xlabel('ua/(km*h^-1)'); ylabel('Pe/kw'); gtext('￠?'),gtext('￠ò'),gtext('￠ó'),gtext('￠?'),gtext('￠?'),gtext('(Pf+ Pw)/nT'),gtext('Pe');
MATLAB画驱动力-行驶阻力图,两个程序
m1=2000 m2=1800; m=3880; r0=0.367 gt=0.85; f=0.013; CDA=2.77; i0=5.83; If=0.218; Iw1=1.798; Iw2=3.598; Ig5=[5.56 2.769 1.644 1.00 0.793]; for i=1:3401; n(i)=i+599; Ttq(i)=-19.313+295.27*(n(i)/1000)-165.44*(n(i)/1000).^2+40.874*(n(i)/1000).^3-3.8445*(n(i)/100 0).^4; end for i=1:3401; for j=1:5; Ft(i,j)=Ttq(i)*i0*Ig5(j)*gt/r0; ua(i,j)=0.377*r0*n(i)./(Ig5(j)*i0); F(i,j)=f*m*9.8+CDA*ua(i,j).^2/21.15; end end plot(ua,Ft,ua,F) xlabel('ua/(km/h)'); ylabel('F/N'); title('汽车驱动力-行驶阻力平衡图'); gtext('Ft1') gtext('Ft2'); gtext('Ft3'); gtext('Ft4'); gtext('Ft5'); gtext('Ff+Fw');
>> m=3880; >> g=9.8; >> G=m*g; >> ig=[5.56 2.769 1.644 1.00 0.793]; >> nT=0.85; >> r=0.367; >> f=0.013; >> CDA=2.77; >> i0=5.83; >> L=3.2; >> a=1.947; >> hg=0.9; >> If=0.218; Iw1=1.798; Iw2=3.598; nmin=600;nmax=4000; >> n=600:10:4000; >> Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; >> Ft1=Tq*ig(1)*i0*nT/r; >> Ft2=Tq*ig(2)*i0*nT/r; >> Ft3=Tq*ig(3)*i0*nT/r; >> Ft4=Tq*ig(4)*i0*nT/r; >> Ft5=Tq*ig(5)*i0*nT/r; >> ua1=0.377*r*n/ig(1)/i0; >> ua2=0.377*r*n/ig(2)/i0; >> ua3=0.377*r*n/ig(3)/i0; >> ua4=0.377*r*n/ig(4)/i0; >> ua5=0.377*r*n/ig(5)/i0; >> ua=0:5:120; >> Ff=G*f; >> Fw=CDA*ua.^2/21.15; >> Fz=Ff+Fw; >> plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua,Fz); >> title('汽车驱动力-行驶阻力平衡图'); >> xlabel('ua(km*h^-1)'); >> ylabel('Ft/N'); >> gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ft5'),gtext('Ff+Fw');