1528
.pdf
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
P |
|
n2 1 |
1 |
||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||
58 |
P |
|
c |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
1 |
|
1 |
|
|
|
|
|
1 |
|
|
|
||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
n 1 |
|
|
|
|
n 1 |
|||||
|
|
|
|
|
|
|
|
1 |
|
|
|
|||||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|||
59 |
P = μ P' |
i |
- φ(P - P ) |
|
|
|
|
|
|
|
||||||
|
i |
|
|
|
|
r |
a |
|
|
|
|
|
|
|
||
60 |
i |
0,008315 |
PiLTk |
|
|
|
|
|
|
|
||||||
Hu V Pk |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
61gi = 3600 / (Hu ηi)
62Pм = 0,089+0,0118 Sp n / 30
63Pe = Pi - Pм
64ηм = Pe / Pi
65ηe= ηi ηм
66ge= gi / ηм
n2 1
PI1:=PC*(LM*(R-1)+LM*R*(1-1/ /exp(NP-1)*Ln(D)))/(NP-1)-(1-1/ /(exp(NC-1)*Ln(S)))/(NC-1))/(-1)
PI:=PI1*MI-FI*(PR-PA)
KI:=0.008315*PI*LT*TK/
/(HU*KV*PK)
GI:=3600/(HU*KI)
PM:=0.089+0.0118*SP*N/30
PE:=PI-PM
KM:=PE/PI
KE:=KI*KM
GE:=GI/KM
Окончание табл. 12
67 |
N'i = N'e / ηм |
||
68 |
D |
|
12Ni1 |
|
ц |
3,14PinNiSp |
|
|
|
|
|
69iVh = 10πDц2Spi / 4
70Ne = P'e iVh n / (30τ)
71Me = 3 104Ne / (π n)
72Ni = Ne/ ηм
73GT = ge Ne/ 1000
74Q = GT Hu
75Qe= 3,6Ne
76qe = Qe 100 / Q
77Qм = 3,6(Ni - Ne)
78qм = Qм 100 / Q
79tr = Tr - 273
80to = To - 273
NI1:=NE1/KM
DC:=sqrt(12*NI1*T/(3.14*PI*N*I*SP)
VHI:=10*3.14*sqr(DC)*SP*IC/4
NE:=PE*VHI*N/(30*T)
ME:=30000Xne/(3.14*N)
NI:=NE/KM
GT:=GE*NE/1000
Q:=GT*HU
QE:=3.6*NE
QE1:=QE*1000/Q
QM:=3.6*(NI-NE)
QM1:=QM*100/Q
TR1:=TR-273
TO1=TO-273
81(mcp)'r = 30,095+0,0031(tr - 500)
82(mcp)''1r = 32,515+0,00418(tr - 500)
83(mcp)''r = 1/α(mcp)''1r+(α-1) / α(mcp)'r
84(mcp)'o = 29,073+0,0008to
85Qог GТ L mcp rtr mcp oto
1000
86qог = Qог1000 / Q
87Qw = c i Dц1+2m nm / α
88qw= Qw 100 / Q
89Qs = Q - (Qe+Qм+Qог+Qw)
90qs = Qs 100 / Q
MCR1:=30.095+0.0031*(TR1-500) MPR1:=32.515+0.00418*(TR1-500) MPR:=MPR1/AF+MCR*(AF-1)/AF
MPO:=29.073+0.0008*TO
QG:=GT*LT*(B*MPR*TR1-
-MPO*TO1)
QG1:=QG*100/Q
QW:=CW*IC*exp((1+2*MW)*Ln(DC))
*
*exp((MW)*Ln(N))/AF
QW1:=QW*100/Q
QS:=Q-(QE+QM+QG+QW)
QS1:=QS*100/Q
Начало
Ввод
исходных
данных
Вычисление
Lo
Задание n(ц)
Вычисление
α, L, M1, MCO,
MH2O, MN2, MO2
|
|
Вычисление |
|
|
|
M2, , TkN, Tk, k, |
|
|
|
Sp, Wвп |
|
|
|
Вычисление |
|
|
|
Pa, PrN, Pr |
|
Рис. 5. |
Схема алгоритма теплового расчёта |
|
|
дизельного двигателя с наддувом (начало) |
|
||
|
|
Задание |
|
|
|
(ц) |
|
|
|
Вычисление |
|
|
|
Tr, ∆T, γr, Ta, ηV |
|
|
|
Вычисление |
|
|
|
(mcv) c, ta, (mcv) a, |
|
|
|
N1, n1, Tc, ∆tc |
|
|
|
< 0 |
> 2 |
|
tc = tc - 5 |
∆tc |
tc = tc + 5 |
|
Вычисление |
|
|
Pr, , U c, U 1c, |
|
|
U c |
|
n ≠ nN |
|
n = nN |
|
Проверка |
ξ = ξN |
|
ξ |
|
|
|
|
|
Вычисление |
|
|
ξ = f(n) |
|
Продолжение рис. 5 |
|
|
|
Вычисление |
|
|
H1, tz, Tz, Pz, |
|
|
, δ |
|
|
Вычисление |
|
|
(mcv) 1z, (mcv) z, |
|
|
(mcv) z |
|
|
Вычисление |
|
|
(mcv) 1b, (mcv) b, |
|
|
(mcv) b |
|
|
Вычисление |
|
|
K2, n2, Tb, ∆Tb |
|
tb = tb + 5 |
∆Tb |
tb = tb - 5 |
<0 >2
<0 >2
Продолжение рис. 5
Вычисление
P i, Pi, ηi, Pм, Pе, ηм, ηе, gе
n ≠ nN |
n = nN |
|
n, IR |
IR≠1 |
IR=1 |
|
Вычисление |
|
NiN, Dц, iVh |
Вычисление
Nе, Mе, Ni, GT, Qе, qе, Qм, qм
Вычисление tr, to, (mcv) r
Вычисление
(mcp) 1r, (mcp) r, (mcp) o
Вычисление
Qог, qог, Qw, qw, Qs, qs
|
|
Продолжение рис. 5 |
|
||
|
|
< 2 |
>5 |
|
|
|
|
|
|
||
|
|
qs |
|
|
|
|
m = m-0,005 |
|
m = m+0,005 |
|
|
|
|
|
|||
|
|
|
|
|
|
Печатать
Pz, gе, ηе, Nе
Печатать
Mе, GT, qе
Печатать
qм, qог, qw,
|
|
>ξmin |
< ξmax |
|
|
|
|
||
|
ξ =ξ-∆ ξ |
ξ |
|
ξ =ξ+∆ ξ |
|
|
|||
|
|
|
|
|
nmin
А |
n = n+∆n |
|
Выход
Окончание рис. 5
5.9. Программа теплового расчета дизельного двигателя с наддувом
C
REAL N1,MT,NK,KCI1,LM,MI,NE1,IC,MW,LO,LT,M1,MC,MH,MN,MO,M2,
*KV,MCC,MCA,KC,NC,KCI,MPZ1,MCZ,MPZ,MPB1,MCB,MPB,KP,NP,KI,KM,
*KE,NI1,NE,ME,NI,MCR,MPR1,MPR,MPO,K,N3,N OPEN(5,FILE='ISX_DAN.DAT')
OPEN(6,FILE='REZYL.TXT')
C OPEN(7,FILE='CON')
READ (5,1)
C,H,O,AF1,N1,MT,TO,PK1,PO,NK,RB,CN,FK,BS,TRN,DT1,S,V,TC1 READ (5,3) L,M,I
READ (5,5) KCI1,HU,LM,TB1,MI,FI,NE1,T,IC,CW,MW,K,N3
WRITE (6,2) C,H,O,AF1,N1,MT,TO,PK1,PO,NK,RB,CN,FK,BS, *TRN,DT1,S,V,TC1
WRITE (6,4) L,M,I
WRITE (6,2) KCI1,HU,LM,TB1,MI,FI,NE1,T,IC,CW,MW,K,N3
1FORMAT (7(/),(63X,F8.0))
2FORMAT (1X,F12.5)
3FORMAT (63X,I8)
4FORMAT (1X,I12)
5FORMAT (63X,F8.0)
C |
ПАРАМЕТРЫ РАБОЧЕГО ТЕЛА |
LO=(C+3*(H-O/8))/(12*0.209) N=N1
10 AF= (AF1-0.6)+(N/N1-0.3)**0.75-0.2357*(N/N1-0.3) LT=AF*LO
M1=LT+1/MT
MC=C/12.0
MH=H/2.0
MN=0.79*AF*LO MO=0.209*(AF-1)*LO M2=MC+MH+MN+MO
B1=M2/M1
WRITE(6,50)M2
50 FORMAT('M2=',E11.4)
C |
ПРОЦЕСС ВПУСКА |
TK1=TO*(PK1/PO)**((NK-1)/NK)
T1K=TK1-273
TK=(T1K-56)+112*(N/N1-0.5)+273
PK=(PK1-0.056)+0.112*(N/N1-0.5)
RK=PK*10**6/(RB*TK)
SP=30*CN/N1
W=0.05433*SP*N*FK
11PA=PK-BS*RK*W**2/2E6 PR1=0.9*PK
PR=(PR1-0.056)+0.07*(N/N1-0.2)
12TR=(TRN+45)-90*(N/N1-0.5)+273
DT=DT1*(110+0.0125*N)/(110+0.0125*N1)
13Y=(TK+DT)*PR/(TR*(PA*S-PR))
14TA=(TK+DT+V*Y*TR)/(1+Y)
15KV=TK*(PA*S-PR)/((TK+DT)*PK*(S-1))
WRITE(6,51)KV
51 FORMAT('KV=',E11.4)
C |
ПРОЦЕСС СЖАТИЯ |
16MCC=22.408+0.00299*(TC1-700)
TA1=TA-273
MCA=20.838+0.00146*(TA1-100) KC=1+8.315*(TC1-TA1)/(MCC*TC1-MCA*TA1)
NC=KC
17TC=TA*S**(NC-1)
DTC=((TC-273)-TC1)*100/(TC-273)
IF(DTC.GT.2.0) GO TO 18
IF(DTC.LT.0.0) GO TO 19
GO TO 20
18TC1=TC1+5
GO TO 16
19TC1=TC1-5
GO TO 16
20CONTINUE PC=PA*S**NC
WRITE(6,52)PC
52 FORMAT ('PC=',E11.4)
C |
ПРОЦЕСС СГОРАНИЯ |
B=(B1+Y)/(1+Y)
IRAB=0
DO 46 JJ=1,M
IRAB= IRAB+1
J=L+(JJ-1)*I
KCI1=J/1000.0
IF(N/N1.GE.1.0) GO TO 21
GO TO 22
21KCI=KCI1 GO TO 60
22KCI=(KCI1-0.2)+(N/N1-0.3)**0.75-0.807*(N/N1-0.3)
60CONTINUE
UCC=15.684+0.02496*(TC1-700)
UPC1=17.585+0.02826*(TC1-700)
UPC=UPC1/AF+UCC*(AF-1)/AF H1=((KCI*HU/LT+UCC+Y*UPC)/(1+Y)+0.008315*LM*TC)/B
TZ1=(H1+(3.59+2.942/AF))/(0.036575+0.00524/AF)
TZ=TZ1+273
PZ=LM*PC
R=B*TZ/(LM*TC)
D=S/R WRITE(6, 53)TZ
53 FORMAT( TZ= , E11.4)
C |
ПРОЦЕСС РАСШИРЕНИЯ |
|
MPZ1=28.638+0.0023*(TZ1-1800) |
|
MCZ=25.003+0.001615*(TZ1-1800) |
|
MPZ= MPZ1/AF+MCZ*(AF-1)/AF |
23 MPB1=25.498+0.00376*(TB1-800)
MCB=22.713+0.00293*(TB1-800)
MPB=MPB1/AF+MCB*(AF-1)/(AF) KP=1+8.315*(TZ1-TB1)/(MPZ*TZ1-MPB*TB1)
NP=KP
TB=TZ/D**(NP-1) DTB=(((TB-273)-TB1)*100)/(TB-273)
WRITE (6, 57) TB
57 FORMAT ('TB=', E11.4)
IF(DTB.GT.2.0) GO TO 24
IF(DTB.LT.0.0) GO TO 25
GO TO 26
24TB1=TB1+5
GO TO 23
25TB1=TB1-5 GO TO 23
26CONTINUE
PB=PZ/D**NP
TRI=TB/(PB/PR)**0.333
DTR=(((TRN+273)-TRI)*100)/(TRN+273)
IF(DTR.GT.2.0) GO TO 27 IF(DTR.LT.0.0) GO TO 28 GO TO 29
27TRN=TRN-5
GO TO 12
28TRN=TRN+5 GO TO 12
29CONTINUE WRITE (6, 54) PB
54 FORMAT ('PB=', E11.4)