1524

4.9. Программа теплового расчета двигателя, работающего на сжиженном газе

REAL MT1,N1T,M1T,MC3,MT2,N2T,M2T,MC4,KT1,KT2,N1,KCI1,MI,NE1,IC, *MW,K,LO1,LO2,LO,LT,MT,M1,MCO,MC,MH2,MH,MN,V2,KV,MCC,MPC,MCM, *MCA,MPA,MAC,KC,NC,KCI,LM,MPZ,MPB,KP,NP,KI,KM,KE,NI1,NE,ME, *MPR,MOC,MOP,N3,N,M2

OPEN(5,FILE='ISX_DAN.DAT')

OPEN(6,FILE='REZYL.TXT') C OPEN(7,FILE='CON')

READ (5, 1)MT1,N1T,M1T,MC3,MT2,N2T,M2T,MC4,AF1,KT1,KT2,PO,TO,S READ (5,3) L,M,I

READ (5,5)RB,CN,N1,FK,BS,TRN,FN,DT1,V,TC1,KCI1,HU,TB1,MI,FI *,NE1,T,IC,CW,MW,N3,K

WRITE (6,2)MT1,N1T,M1T,MC3,MT2,N2T,M2T,MC4,AF1,KT1,KT2,PO,TO,S WRITE (6,4)L,M,I

WRITE (6,2)RB,CN,N1,FK,BS,TRN,FN,DT1,V,TC1,KCI1,HU,TB1,MI,FI *,NE1,T,IC,CW,MW,N3,K

SP=30*CN/N1

LO1=(N1T+M1T/4)*MC3/(0.209*MT1)

LO2=(N2T+M2T/4)*MC4/(0.209*MT2)

LO=LO1+LO2

N=N1

9 A=0.7-N/N1

IF (A.GE.0.0) GO TO 10 AF=AF1

GO TO 11

10AF=AF1-0.78*A**2.36

11CONTINUE LT=AF*LO MT=MC3*MT1+MC4*MT2 M1=LT+1/MT

MCO=0.42*(1-AF)*(LO1/(1+KT1)+LO2/(1+KT2)) MC=N1T*MC3/MT1+N2T*MC4/MT2-MCO MH2=MCO*(KT1+KT2)/2 MH=M1T*MC3/(2*MT1)+M2T*MC4/(2*MT2)-MH2 MN=0.79*AF*LO

M2=MCO+MC+MH2+MH+MN

B1=M2/M1 WRITE (6,60) M2

60FORMAT (1X,'M2=',E11.4)

RO=PO*10**6/(RB*TO)

W=0.05433*SP*N1*FK

12PA=PO-BS*W**2*RO/2E6 PR1=0.098+0.539*N1/10**5

13PR=1.035*PO+(PR1-1.035*PO)*(N/N1)**2 IRAB=0

20TC1=TC1+5

GO TO 18

21TC1=TC1-5 GO TO 18

22CONTINUE

24TZ=TZ1+273

25PZ=PC*B*TZ/TC

LM=PZ/PC WRITE (6,63) TZ

63FORMAT(1X,'TZ=',E11.4)

IF (DTB.GT.2.0) GO TO 31

IF (DTB.LT.0.0) GO TO 32

GO TO 33

31TB1=TB1+5 GO TO 30

32TB1=TB1-5 GO TO 30

33CONTINUE PB=PZ/S**NP WRITE (6, 64) PB

64FORMAT (' PB = ', E11.4) TRI=TB/(PB/PR)**0.333 DTR=(TRN-TRI)*100/TRN IF (DTR.GT.2.0) GO TO 34 IF (DTR.LT.0.0) GO TO 35 GO TO 36

34TRN=TRN-5

GO TO 14

35TRN=TRN+5 GO TO 14

36CONTINUE

GO TO 46

48MW=MW-0.005 GO TO 46

49CONTINUE

WRITE(6,50)

50 FORMAT (1X,30('*')) WRITE (6,70) S,N

70 FORMAT (1X,'S=',E11.4,' N=',E11.4) WRITE (6,51) PZ, GE, KE

WRITE (6,52) NE, ME, GT WRITE (6,53) QE1, QM1, QG1 WRITE (6,54) QH1, QW1, QS1

51 FORMAT (3X,'PZ=' ,F5.3,2X,'GE=' ,F5.1,2X,'KE=' ,F6.4)

52 FORMAT (3X,'NE=' ,F6.1,2X,'ME=' ,F5.1,2X,'GT=' ,F5.2)

53 FORMAT (3X,'QE1=',F6.3,2X,'QМ1=',F5.3,2X,'QG1=',F6.3)

54 FORMAT (3X,'QH1=',F5.2,2X,'QW1=',F5.2,2X,'QS1=',F4.2)

55 CONTINUE N=N-K

IF (N. GE. N3) GO TO 9 STOP

END

Контрольные вопросы

1.Экономическая эффективность использования сжиженного газа в двигателях с воспламенением рабочей смеси от электрической искры.

2.Как изменится токсичность двигателя при переводе его на питание

сбензина на сжиженный газ?

3.Достоинства и недостатки двигателя, работающего на сжиженном

газе.

4.Чем ограничивается степень сжатия при использовании легких топлив (бензина, сжиженного газа, природного газа)?

5.Что называется литровой мощностью двигателя, и какие факторы влияют на её величину?

6.Влияние какого параметра на мощностные и экономические показатели работы двигателя определяется в программе?

7.Какие показатели работы двигателя выведены на печать в составленной программе?

Глава 5. ТЕПЛОВОЙ РАСЧЕТ ДИЗЕЛЬНОГО ДВИГАТЕЛЯ С НАДДУВОМ

Выполнить тепловой расчет дизельного двигателя, предназначенного для грузового автомобиля, и имеющего следующие исходные данные:

эффективную мощность Nе = 190 кВт;

частоту вращения коленчатого вала nN = 2600 об /мин;

число цилиндров i = 8;

применяемое топливо – дизельное;

цетановое число топлива – 45;

состав топлива C = 0,87 кг/кг топл., H = 0,126 кг/кг топл., О = 0,004 кг/кг топл.;

молекулярная масса топлива mТ = 190;

низшая теплота сгорания Hu =42,5 МДж/кг;

степень сжатия ε = 16;

коэффициент избытка воздуха α =1,8;

давление наддува PkN = 0,18.

5.1. Определение параметров рабочего тела

Теоретический расход воздуха, необходимый для сгорания 1 кг топ-

Так как при работе двигателя по внешней скоростной характеристике изменяется интенсивность трубонаддува, то коэффициент избытка воздуха для различных скоростных режимов можно определять по формуле

Действительный расход воздуха для номинального режима работы двигателя

L = αLо = 1,8 0,497 = 0,8946 кмоль.

Количество свежего заряда

M1 = L + 1/mТ = 0,8946 + 1/190 = 0,8998 кмоль.

Количество отдельных компонентов, входящих в состав продуктов сгорания:

Mco2 = C / 12 = 0,87/12 = 0,0725 кмоль; MH2O =H / 2 = 0,126/2 = 0,063 кмоль;

MN = 0,79α Lо = 0,79 1,8 0,497 = 0,7067 кмоль;

2

MO = 0,209(α-1) L = 0,209 (1,8-1) 0,497 = 0,0831 кмоль.

2

Количество продуктов сгорания

M2 = MCO +MH O+MN +MO = 2 2 2 2

= 0,0725+0,063+0,7067+0,0831= 0,9253 кмоль.

Химический коэффициент молекулярного изменения

β = М2 / М1= 0,9253 / 0,8998 = 1,0283.

5.2. Процесс впуска

Параметры окружающей среды: Ро = 0,1 МПа, То = 293 K.

Степень подогрева свежего заряда на номинальном режиме двигателя ТN= 10 . Среднюю скорость поршня на этом же режиме принимаем CN= =10,4 м/с.

Температура воздуха после компрессора

гдеnk показательполитропы сжатия компрессора,принимаем равным 1,7. Плотность воздуха, поступающего в двигатель на различных скорост-

ных режимах будет изменяться в зависимости от изменений значений Тk и

Pk, которые можно определять по формулам:

n

Тk = (tkN - 56)+112 ( nN - 0,5)+273; tkN = TkN - 273;

n

Pk = (PkN - 0,056)+0,112 ( nN - 0,5).

Плотность воздуха, поступающего в двигатель, на номинальном режиме будет равна

Ход поршня при номинальной частоте вращения коленчатого вала

Sp = 30 CN / nN = 30 10,4/2600 = 0,12 м.

Отношение площади поршня к площади клапана Fn / fкл = 6.

Максимальная скорость свежего заряда в проходном сечении клапана

Wвп = 0,05433 0,12 2600 6 = 101,7 м/с.

Сопротивление выпускаемой системы с учетом коэффициента затухания скорости

(β2+ ) = 2,7.

Давление в цилиндре двигателя в конце впуска

Pa =PkN (β2+ ) Wвп2 ρk / (2 106) = 0,18 2,7101,72 1,6814/(2 106) = 0,1565 МПа.

Давление остаточных газов на номинальном режиме работы двигателя

PrN = 0,9 PkN = 0,9 0,18=0,162.

Давление остаточных газов для различных скоростных режимов в первом приближении будет изменяться по тому же закону как и давление наддува.

n

Pr =(PrN 0,056)+0,07 ( nN 0,2).

Температура остаточных газов для различных скоростных режимов

Tr = (trN + 45) 90 ( n 0,5) + 273. nN

Принимаем TrN = 780 К, тогда trN = 780 273 = 507 С.

Изменение степени подогрева свежего заряда в зависимости от скоростного режима двигателя можно принять таким же, как и для карбюраторных двигателей, т.е.

110 0,0125n

Т ТN 110 0,0125nN .

Коэффициент остаточных газов

здесь ν = Сpr / Cpо = 1,14 (Срr теплоемкость продуктов сгорания при

постоянном давлении, Cpо – теплоемкость воздуха при постоянном давлении).

Коэффициент наполнения

При определении параметров рабочего тела в конце впуска не учитывались дозарядка и продувка.

5.3. Процесс сжатия

Для определения показателя адиабаты сжатия принимаем нижний предел температуры рабочего тела в конце сжатия tc1 = 700 C, расчетную температуру tс = 810 C.

Теплоемкость свежего заряда в точке С (см. рис.1) согласно таблицы 10 будет равна

(mcν)'c= 22,408+0,00299 (tc 700) = 22,408+ +0,00299 (810 700) = 22,7369 кДж/кмоль.

Теплоемкость свежего заряда в точке А

(mcν)'а = 20,838+0,00146 (tа 100) = 20,8380,00146 (127 100) = 20,878 кДж/кмоль.

Показатель адиабаты сжатия

Принимаем n1 = К1.

Температура рабочего тела в конце сжатия

Тс = Та ε n1 1 = 400 161,36 1 = 1085 К.

Ошибка в выборе расчетной температуры рабочего тела в конце сжа-тия составит

Если ошибка будет больше 2 %, необходимо задать новое значение tс и расчет повторить.

Давление в конце сжатия

Рс = Ра εn1= 0,1565 161,36= 6,794 МПа.

В расчетах не учитывалась разница между теплоемкостью остаточных газов и теплоемкостью воздуха.