C/4LL GRBEGK' M0T0R201. C9 in' , 12. XPG, YPG.O, if on)

WRITE<*,«> ' Ждите '.Мгаг. 170 IF(NomVib.eq.l) E=X<170) IF(NomVib.eq.2) N=X<I70> IF<NomVib.eq.4) KDostg=X(170) IF(NomVib.eq.S) Tr=X(170) IF(NomVib.eq.6) P0=X(I70) IF(NomVib.eq.7) TO=X(I70> IF(NomVib.eq.8) ALFA=XH70> IF(NomVib.eq.3) DT=X(I70) IF(NomVib.eq.lO) KDP»=X<I70> ^

IF(NomVib.eq.ll) PSI=X(I70)

Write(4,173) FIO

Write(4.173) PUSTO

write<4.173) PUST01

wri te (4. 175) 1С. DC. SP, LAMK. MP. MS, E, N. NEN

WRITE(4,173) PUST02

WRIТЕ (4. 176) IRADN, UGR/52. JK. TURBO, NK, STEPND, KdOSTG, Тг, Р0, TO

WRITE<4,173) PUST03

WRITE<4.177) ALF/5.DT, KDP/Ч, PSI, L^MD, i f on. Ymaxl, i ri s write(4,74) NomVib, Мгаг, (X < I). 1 = 1, Мгаг). Sh

175 F0RM/1T< I2.2F5.0.1X,F4.3.2F6.3.F6.2,F6.0,F7.2>

176 FORMAT(I2.2X.F4.0.2X. I1,4X, 11.3X.F3.1, IX, F3.1.2X.F5. 1 K,5:.:,f5.0,f7.3,f6.0)

177 FORM/ЧТ<F3.1,2X,F3.0. IX, F3.2.2X, F3. 2. IX. F4.2. i4, f6.1. i4> 74 format <'Nomvib,Мгаг, X,__='_#2i4.8f9.3)

CV ==:—=====" дизель—бензин и р-ыбор исходных данных ■♦■ +++■*•*•++

JD=1 'ДИЗЕАЬ

IF(E.GE.13) GOTO 46 IF(E.LE.ll) GOTO 47

GOTO 200

47 JD=2! КАРБЮРАТОРНЫЙ

C=0.854 H=0.144 0t=0.0 QN=44000 L0m=14.96 46 CONTINUE

IF(JK.EQ.3.AND. JD.EQ.i) GO TO 200 IF<JK.LE.2.AND. JD.EQ.2) GO TO 200

CV «»**«* НАЧАЛО ТЕПЛОВОГО РАСЧЕТА ДВИГАТЕЛЯ **«««» WRITE(4,160) write(4,56)

write (4, 561) L0m,C,H, Ot.QN IF<NomVib.eq.3> Stepnd=X(170) PK=POKSTEPND*«Turbo TK=TO*<PK/P0> «к((NK-1)/NK) PA=PK*<1-KDPA) PR=PK*KdOSTG IF(NoroVib.eq.21) Tt:=X(I70)

JR=(<TK+DT) /TR) K<PR/<E»PA-PR) > IF(NomVib.eq.22) Jr=X<I70)

IF(NomVib.eq.22> Pa= (((TK+DT) /TR) *Pr -e PRitJr >/(E*Jr) TA=(TK+DT+JRkTR) /d-KJR)

ETAV=(E/(E-1))к(РА/РК)к(ТК/(TK+DT+JR*TR)> Nl=1.41-110/N N2=1.22+130/N

IF(JD.EQ. 1) N2=N2-0.01! ДИЗЕЛЬ

IF(NomVib.eq.23> ETAv=X(I70>

IF(NomVib.eq.23> Pa=ETAV*Pk/((E/ (E-l)) к(ТК/ (TK+-DT+JR*TR> >) / IF(NomVib.eq.24> Nl=X(I70)

IF(NoroVib.eq.25) N2=X(I70) WRITE(4,151) ТА,TK,PK,PA,PR,JR,ETAV PC=PA»E»*N1 TC=TA*Ek«<N1-1> CV ========== 1.4.1 Термохимический расчет ============

L0=(C/12+H/4-0t/32)/0.21! Кмоль

L0MR=L0*28.96

ml=ALFAKLO

MR=JR*ml

MC:=ml+MR

M2=ml+H/4+0t/32

IF(JD.EQ.2) M2=M2+0.21kLOk(1-ALFA)! BEN2IN

DM=M2-M1

BO=M2/M1

MZ=M2+MR

B=MZ/MC

Су =====»==== 1.4.2 Термодинамический расчет ============

IF<JD.EQ.2> GO TO 212! БЕНЗИН

AA=B*(13.B/ALFA+15.5> /10000! ДИЗЕЛЬ

BB=B*(8.314+(20.036+O.32/ALFA>) CCC—< (20.16+1.74*TC/1000+8.314*LAMD> *TC+PSI *QN/

* <ALFA*LO*(1+JR> > >

T2=(-BB+SQRT(BB*BB-4*AA*CCC> > /(2*AA> PZ=LAMD«PC RO=B*TZ/(LAMD*TC> GO TO 213

212 AA=B*< 15.431+13.716*ALFA> /10000 III BEN21N BB=B*<18.4+2.6»sALFA)!l!

CCC=-< (20. 1В + 1.74*ТС/1ОО0> *TC+PSI*<QN-120*1000*LO*

* (i-ALFA) > /<ALFA*LO*<i+JR>) > III

T2=(-BB +SQRT(BB *BB-4 «Л» *CCC> > / < 2 кЛЛ) PZ=0.85*B*PC*TZ/TC I I I БЕНЗИН

RO=1 I I I

CV *=—===«.= ПАРАМЕТРЫ ГАЗА В КОНЦЕ РАСШИРЕНИЯ 1.5.=====»•=

213 DEL=E/RO PB=PZ/DEL**N2 TB=TZ/DEL**(N2-i>

WRITE (4,11) Nl,N2,PC,TC,L0,Ml,MR.MC,M2,DM.B0.,M2 WRITE (4, 111) B,RO, T2,.P2.LAMD,DEL, PB,TB,LOMR 555 f or mat (6:>;, i в) CV ******** ПОСТРОЕНИЕ ИНДИКАТОРНОЙ ДИАГРАММЫ P-V 1.6 ******

AB=SP

AO=AB/<E-1> VA=AO*E VB=AO*E V2=A0*R0 R=AB/2 RL2=R*LAMK/2 XlC=Xl+A0+AB/2 Y1C=Y1-15.

CALL PXY(X1C+RL2,Y1C,COK1> CALL CXYR<X1C,Y1C,AB/2,OKR> XA-X1+A0

CALL MOVE(X1C+AB/2,Y1C,O> CALL MOVE<Xl+AO,ViC,l) CALL M0VE<XlC-AB/2,Y1C,O> CALL ARC:IA<R, 180.,35Э.) call setpen(2) call move(xlc,ylc,0> call circd.)

call move<xic.ylc-sp/2-10., 1) call set pen (1) VXD(1)=AO VXD(13)=AB+A0

CV ****** Д/1Я 4 MKR И ИНДИКАТОРНОЙ ДИАГРАММЫ к*к***

CV к» P-V свернутая и P-FI по методу проф. Брикеа 1.6 ****** YMAX=yroaxl CONSTS=PA*VA«*N1 C0NSTR=PB*VBK*N2 WRITE<4,154) AO.AB DO 23 I»1.13 UGOL=<I-1)«15.

IF(I.EQ.1.0R.i.eq.l3> GOTO 124 CALL LPAX<COK1,UGOL,PRM1> CALL PILC(PRMl.OKR.J) CALL SORTP(4,TOC:H1> VXD< I)=TOCH1 (1) -XA+AO CALL MOVE(COKKl),COK1<2),0) CALL MOVE<TOCH1<15 _#TOC:H1<2),1) CALL MOVE<TOC:H1(1>,Y1.1> 124 VXDN1=VXD(I)**N1

VXDN2=VXD <I)**N2 PVPD<I)"PA

PSD(I)=C:ONSTS/VXDN1 I Сжатие PRD(I)=C0NSTR/VXDN2! Расширение IF(VXD<I>.LT.V2) PRD<I)=PZ PVXD <I)=PR XP1=VXD(I)-АО WRITE <4,156) I,UGOL,XP1.PVPD <I>. PSD(I),PRD(I).PVXD <I)

VRl=VXD<13)*1.05 VR2=VR1*1.2 IF(I.EQ.13) GOTO 23 CALL MOVE < VXD П) +X1,Y1,O) YYYY=PRD <I)«VR2/YMAX CALL MOVE<VXD(1)+X1,Yl+YYYY.1) 23 CONTINUE

if<JD.EQ.l) PRD<2>-PRD<1) PRD (13) = (PR+-PB) /2 PVXD<13)=(PR+PB)/2 PVPD<1)=PR YK=YMAX IF<I70.LT. 1) GO TO 560

CALL SYMBOL<XI-10.,Y1+VR2-5.,Э.,'p'.l.O.) CALL SYMBOL<Xl+vrl+6._#Yl-10.,3.,'Xp -.V^,g,O.) call symbol (xl+vi-1+5., yl+100. .Э.. AQ='. 3,0.) call number <xl+vrl+-5. +20, у 1+100.,Э.,АО, 3,0.) call symbol (xl+vrl+5., yl+80,9., ' AB=',3.0.) call number < x 1 +vr 1 +5. +20, yl +60,,9.,AB,1,O.) call symbol (xl+vrl+5.. yl+60.,Э., * SP=' -3.0.) call number<xl+vrl+5.+20.yl+60.,9.,SP,1,0.)

call number<xl+vri+iO.,yl+120.,5.,Nomvib+0.,O,0.>

cal 1 symbol <xl+vrl+5.. yl+120..9.. ' X'. 1.0.) 560 CONTINUE

CALL LGRfXl.Yl,VR1.VR2,0.,VR1,O.,YK,VXD,PVPD,13)!Впуск CALL SETPEN<2>

CALL LGR<X1,Yl,VR1.VR2,0.,VR1,O.,YK, VXD,PSD,13) Гсжатие CALL SETPEN(4)

CALL LGR(X1,Y1,VR1,VR2,O.,VR1,O.,YK,VXD,PRD.13) (расширение CALL ЗЕТРЕЖЗ) ♦

CALL LGR(X1,Yl,VR1,VR2,0.,VR1,O.,YK,VXD,PVXD.13) I выхлоп CALL SETPEN(l) 18=1

IF(JD.EO. 1) 18=2 call maxos(yk,yosd)

CALL LGO<X1, Yl, Vftl, VR2,0.,VR1,O..YK.30.,yo»d) snd«<prd<l)-PSD<3))/3O WRITE (4,504) prd<i),psd(3),'snd

504 FORMAT(' скорость нарастания давления (prd(1)-psd(3))/30 «snd • */,29X,3f9.5)

CV ********** РАСЧЕТ ПОКАЗАТЕЛЕЙ РАБОЧЕГО ЦИКЛА 1.7 ************ IF(JD.EQ.l) Pli = (PC/(E-l>)»(LAMD*(RO-1) +LAMD*RO*

* (i-l/DEL»*(N2-l>)/(N2-l) -(1-1/E**(N1-1) >/(Nl-l))

IF(JD.EQ.2> PI1=(PC:/(E-1)) *

* (LAMD*(1-1/E**(N2-1) >/(N2-l) -(1-1/E**(N1-1))/(Nl-D) CV PIND=PI1*O.94-(PR-PA)

PIND=PIl*0.94 SR=O

ss=o

DO 101 1=2,13

SS=SS+(PSD<1-1)+PSD(I)) * (VXD(I)-VXD(1-1)) /2 101 SR=SR+(PRD(I-1)+PRD(I)) * (VXD(I)-VXD<1-1)) /2 PLOS=SR-SS

WRITE(4,*) ' ПЛОЩАДЬ между кривыми расш. и сжат. мм*мм = ',PLOS

PINDGR=PLOS/AB

S=PI1/PINDGR

WRITE<4,505) PI1,PINDGR,S 505 FORMATC PI1 / PINDGR = s V3f3.5)

PIND=PINDGR*0.34

IF<JK.EQ.3) PM=O.035+5*N/100000 IF<JK.EQ.2) PM=O.067+5.2*N/100000 IF(JK.EQ.1) PM=O.032+1.01»N/10000 PE=PIND-PM ETAM=PE/PIND ROO=PO*1000000/ (287. *TO) ETAI=PIND*L0m*ALFA/(R00»ETAV*QN/1000) ETAE=ETAI*ETAM GI=3600/(ETAI*QN/1000> GE=3SO0/(ETAE*QN/1000> GT=GE*NEN/1000

CV ********** РАСЧЕТ ОСНОВНЫХ РАЗМЕРОВ ДВИГАТЕЛЯ 1.7.3 «к*****-*»-* VH=PI*DC**2*SP/4000000 VL=VH*IC

NE=PE*VL*N/(30*4) ME=3550*NE/N GT=NE*GE/1000. FP=PI»DC*DC/4 VPS=SP*N/(30*1000) NL=NE/VL

NP=NE/(IC*FP)*100*100 SKl=(0.0136*PE*GE*TK/T0)**0.Э8 SK2=(DC/(ETAv*Pk*1000)) **0.38 Qp=l*S4rt(VPS)*SK1*SK2 YG1(170)=Тг YG2(170)=Pz YG3(I70)=Ge YG4 (170) =Ne YG5(170)=Qp

WRI ТЕ < 4, 14) PI 1, PIND, VPS, PM, PE, ETAM, ROO, ETA I, ETAE, GE, GT WRITE(4,142) VL, VH,DC,NE,ME, SP,FP,NP,NL'Qp "t,^,«>

write!*,*) ' Расчет закончен. Нажмите Enter ' Cal 1 grend

j Call Grbe9l('Motor202.C9(n^I,12,200.,200.,0,ifon)

dl —00 ■

CALL SYMBOLdO.,91.,6.,'Мощность Ne, кВт",16,0.) CALL SYMBOH 10., 185.,6., ' Теп.оонапряженность Qn',21.0.) CALL SYMBOL(120.,185.,7.,'Давление Pa',11,О.) CALL SYMB0L(120.,91.,7..'Температура Tz',14,0.)