C |
IRAB3=O |
29 MP1B=MPB1+0.002826*(TB1-1400) |
|
C-------------------------- |
|
C |
IRAB3=IRAB3+1 |
C |
WRITE(6,102)TB1,DTB |
C 102 FORMAT(1X,' TB1=',E11.4,' DTB=',E11.4) |
|
C |
IF(IRAB3.GT.500) STOP |
C-------------------------- |
|
|
AB=MP1B/MPB1 |
|
MPB=MP1B+AB**1.7*(6*(AF-0.7)-3.51*(AF-0.7)**1.75) |
|
KP=1+8.315*(TZ1-TB1)/(MPZ*TZ1-MPB*TB1) |
|
NP=KP |
|
TB=TZ/S**(NP-1) |
C----------------------------------- |
|
C |
WRITE (6,755) TZ1,MPZ,MPB |
C 755 FORMAT(' TZ1=',E11.4,' MPZ=',E11.4,' MPB=',E11.4) |
|
C |
WRITE (6,655) TZ,S,KP,TB |
C |
655 FORMAT(' TZ=',E11.4,' S=',E11.4,' KP=',E11.4,' TB=',E11.4) |
C----------------------------------- |
|
|
DTB=((TB-273)-TB1)*100/(TB-273) |
|
IF(DTB.GT.2.0) GO TO 30 |
|
IF(DTB.LT.0.0) GO TO 31 |
GO TO 32
30 TB1=TB1+5
GO TO 29
31TB1=TB1-5 GO TO 29
32CONTINUE
PB=PZ/S**NP
TRI=TB/(PB/PR)**0.333
DTR=(TRN-TRI)*100/TRN
IF(DTR.GT.2.0) GO TO 33
IF(DTR.LT.0.0) GO TO 34
GO TO 35
33TRN=TRN-5 GO TO 24
34TRN=TRN+5
GO TO 24
35 CONTINUE
C |
ПОКАЗАТЕЛИ РАБОТЫ ЦИКЛА |
|
PI1=PC*(LM*(1-1/S**(NP-1))/(NP-1)-(1-1/S**(NC-1))/(NC-1))/(S-1) |
|
PI=MI*PI1-FI*(PR-PA) |
|
GI=3600*KV/(24.4*LT*PI) |
|
KI=3600/(GI*HU) |
|
PM=0.039+0.0132*SP*N/30 |
|
PE=PI-PM |
|
KM=PE/PI |
|
GE=GI/KM |
|
KE=KI*KM |
|
IF(N.EQ.N1.AND.IRAB.EQ.1) GO TO 36 |
|
GO TO 37 |
36 |
NI1=NE1/KM |
|
DC=SQRT(12*NI1*T/(3.14*PI*N1*IC*SP)) |
|
VHI=3.14*DC**2*SP*IC*2.5 |
37 |
CONTINUE |
|
NE=PE*VHI*N/(30*T) |
|
ME=3*10**4*NE/(3.14*N) |
|
GT=GE*NE/1000.0 |
|
NI=NE/KM |
C |
ТЕПЛОВОЙ БАЛАНС ДВИГАТЕЛЯ |
|
Q=GT*HU |
|
QE=3.6*NE |
|
QE1=QE*100/Q |
|
QM=3.6*(NI-NE) |
|
QM1=QM*100/Q |
|
TR1=TR-273 |
|
TO1=TO-273 |
|
MP1R=MPR1+0.00367*(TR1-700) |
|
AR=MP1R/MPR1 |
|
MPR=MP1R+AR**2.9*(4.956*(AF-0.7)-2.906*(AF-0.7)**1.75) |
|
MRP=MPR+8.315 |
|
MCO=20.759+0.0008*TO1 |
|
MOP=MCO+8.315 |
|
QG=GT*LT*(B*MRP*TR1-MOP*TO1)/1000 |
|
QG1=QG*100/Q |
|
QH=DHU*GT |
|
QH1=QH*100/Q |
C |
IRAB4=O |
38 |
QW=CW*IC*DC**(1+2*MW)*N**MW*(HU-DHU)/(AF*HU) |
C-------------------------- |
|
C |
IRAB4=IRAB4+1 |
C |
WRITE(6,104)MW |
C 104 FORMAT(1X,'MW=',E11.4) |
|
C |
IF(IRAB4.GT.500) STOP |
C--------------------------
QW1=QW*100/Q
QS=Q-(QE+QM+QG+QH+QW)
QS1=QS*100/Q
IF(QS1.GT.3.0) GO TO 39
IF(QS1.LT.1.0) GO TO 40
GO TO 41
39 MW=MW+0.005
GO TO 38
40MW=MW-0.005 GO TO 38
41CONTINUE WRITE(6,45)
45 FORMAT (1X,30('*'))
WRITE(6,55)N,J
55 FORMAT ('N=',E11.4,' J=',I5) WRITE (6,46) PZ, GE, KE
WRITE (6,47) NE, ME, GT
WRITE (6,48) QE1, QM1, QG1
WRITE (6,49) QH1, QW1, QS1
46FORMAT (2X,'PZ=' ,F5.3,2X,'GE=' ,F5.1,2X,'KE=' ,F6.4)
47FORMAT (2X,'NE=' ,F5.1,2X,'ME=' ,F5.1,2X,'GT=' ,F5.2)
48FORMAT (2X,'QE1=',F5.2,2X,'QM1=',F5.3,2X,'QG1=',F6.3)
49FORMAT (2X,'QH1=',F5.2,2X,'QW1=',F5.2,2X,'QS1=',F4.2)
50CONTINUE
N=N-K
IF(N. GE. N3) GO TO 18
STOP
END
Контрольные вопросы
1. Что называется коэффициентом наполнения, коэффициентом остаточных газов и как они влияют на мощностные и экономические показатели двигателя?
2. Чем ограничивается степень сжатия для двигателей легкого топлива?
3. Что называется коэффициентом избытка воздуха
и в каких пределах он лежит для бензиновых двигателей?
4.Чем объяснить более низкую экономичность бензинового двигателя
всравнении с дизельным, несмотря на большую теплотворную спо-
собность бензина и более низкие механические потери?
5.Влияние какого параметра на мощностные и экономические показатели двигателя определяется в составленной программе?
6.Схема введения необходимого параметра в программу, чтобы определить его влияние на показатели работы двигателя.
Глава 3. ТЕПЛОВОЙ РАСЧЕТ ДВИГАТЕЛЯ,
РАБОТАЮЩЕГО НА ПРИРОДНОМ ГАЗЕ
Произвести тепловой расчет четырехтактного двигателя, предназ-наченного для грузового автомобиля, работающего на природном газе и имеющего следующие исходные данные:
эффективная мощность на номинальном режиме работы двигателя
Ne = 110 кВт;
частота вращения коленчатого вала nN = 3200 мин-1;
число цилиндров i =8;
расположение цилиндров V-образное;
степень сжатия =8,5;
применяемое топливо – природный газ.
Состав топлива: метан 94 %, этан – 4 %, пропан – 2 %. Теплотворная способность топлива по формуле Менделеева
Hu = 35,8МСН4 + 63,4 МС 2Н6+91МС 3Н8 = = 35,8 0,94+63,4 0,04+91 0,02 = 38 МДж/м3.
Теплотворная способность одного моля
Hu = Hu 22,4 = 38 22,4 = 851,2 МДж/кмоль.
3.1. Определение параметров рабочего тела
Теоретический расход воздуха, необходимый для сгорания одного кмоля газообразного топлива
L = |
1 |
(n |
m |
)М |
CnHm |
кмоль. |
0,21 |
|
|||||
o |
4 |
|
|
|||
Обозначим:
n число атомов углерода у метана; n число атомов углерода у этана; n число атомов углерода у пропана; m число атомов водорода у метана; m число атомов водорода у этана; m число атомов водорода у пропана.
Теоретическое количество воздуха, необходимое для сгорания метана
|
1 |
|
|
|
|
1 |
|
|
4 |
|
|
||
|
|
|
m |
|
= |
|
|
|
= 8,952 кмоль. |
||||
0,21 |
4 )MCH4 |
0,21 |
(1 |
4) 0,94 |
|||||||||
Lo |
(n |
||||||||||||