Программа для выполнения расчетов и построения графиков на ПЭВМ

CV Motor27.FOR —> 1,2,3.4,5,6,7.8,9,10 6.03.98

IMPLICIT REAUK—X)

RKAL JR.JP

DIMENSION СОК I (2), OKR<3), PRM1 (3),TOCH1(2),VXD(13),PVI'D(13), *PSD<I3),PRD<13),PVXD(13),PRP(49),FK49),PJ(49),PSUM(49), T(49) *,KS(49),MK(49),RSS(49),RKS(49),T1(49),XP(49),VP(49),JP(49)

*,NX(20),NEX(20).MEX(20),GEX(20),GTX(20),SI7.N(25) DATA yl, y2, уЗ, у4, у5, уб, у7, xl, x2, VR2, VR3

* / 100.,100,100.,100., 20., 20., 20., 65.,20., 120., 240./ data ic, SP, DC, LAMK, MP, MS, E, N, NEN

* / 6,140.,130.,0.264,4.150,4.623,16.7,1700,133.2/

DATA TR, C, H, Ot, QN, LOm, ETAV, PO, TO, PR, NK

* /700..0.864.0.126,0.01,41700,14.35, 0.8, 0.1, 298.,0.12,1.6/ DATA ALFA, DT, KDPA, PSI, LAMD, RO

* /1.5, 20., 0.05, 0.95, 1.1, 1.5/

PI-3.14159 RAD-PI/180

FORMAT(1X,7F10.3)

3 ' FORMATC/' КАКОЙ ДВИГАТЕЛЬ? [1—РЯДНЫЙ/2—V-образный]'?: ',$) ■ 4 FORMATC Введите угол развала в градусах: ',$)

5 FORMATC ТУРБОНАДДУВ ЕСТЬ? [1—ДА/0—НЕТ]?: ',$)

И FORMAT(/16X,""***** НАЧАЛО РАСЧЕТОВ ****** ' ♦/вх.'Тешювой РАСЧЕТ '

*/6X,'Nl ¹.7X,^>N2',8X,^>PC\8X,TC^I,8X,'LO\8X,^>M1\8X,'MR' */7X,'MC',8X,'M2',8X,'DM',8X,'B0',8X,'MZ',8X,'B ',8X,'RO' VSX.TZ'.eX.'PZ'JX.'LAMD'JX.'DEL'.SX.'PB'JX.'TB'.SX,' L0MR') 14 FORMAT(8X,'Индикаторные и эффективные показатели, размеры' */6X,'Pir,8X,'PIND',7X,'VPS',7X,'PM',7X,

* 'PE',8X,'ETAM',6X,'ROO76X,'ETAI',7X,'ETAE',7X,

* 'GE ',7X,'GT',7X,'VL',8X,'VH',8X,'DC' */6X,'NE ME SP R FP NP NL')

FORMAT (1X,I2.2(1X,F4.O), F5.3,2(1X,F6.3), lx,f6.2,lX,F6.0,f8.2)

FORMATC 1С SP DC LAMk MP MS E N ',

* ' NEN — ДИЗЕЛЬНЫЙ ДВИГАТЕЛЬ')

FORMATC 1С SP DC LAMk MP MS E N ',

* • neN — КАРБЮРАТОРНЫЙ ДВИГАТЕЛЬ')

FORMAT (1X,I2,1X,2(F4.O,1X), 3(F5-3,1X),2(F4.1,1X), ЗЦ2.1Х))

FORMATU0A4)

FORMAT(2(F3.0,1X),9(F4.2,1X))

FORMAT(5(F5.3,1X),F4.0,1X,F6.0)

Продолжение прил. 15

KORMAT(1X,5(F5.3,1X),F4.0,1X,F6.0)

FORMAT(1X,2(F3.0,1X),9(F4.2,1X))

FORMAT(

*6X,'ALFA',7X,DT',7X,KDPA',7X,'PSI',7X, 'LAMD'.8X, RO',8X, 4.0m' */8X,'C \8X,'H \8X,'Ot',7X,'QN',8X,'ETAV\8X,PO',8X,'TO' */8X,'PR',8X,'TR\8x,Nk')

57 fonnatC Продолжить 1 —да/О—выход/2—возврат на ввод данных: ',$)

58 FORMAT(/' Введите значения 7 *' ALFA — коэффициент избытка воздуха '

*' карб — 0.85 —0.95, диз — 1.45 — 27

•' DT — температура подогрева свежего заряда '

*' карб — -5 — +25, диз — +5 — +40 7

*' KDPA — коэффициент понижения атмосферного давления на впуске7

*' сум.марные гидравлические потери (0.05—0.20)7

*' PSI — коэффициент использования теплоты на участке видимого7

*' сгорания; карб — 0,85 — 0,95; дизель — 0,7 — 0,97

*' LAMD — степень повышения давления карб — 3 — 4; диз — 1.2 — 2.57

*' RO — степень предварительного расширения карб — 1, диз — 1.2 — 2.47

*' Введите значения по образцу '

*/' ALFA DT KDPA PSI LAMD RO '.

*/' 1.5 20. 0.05 0.85 2.0 1.5")

59 FORMATUX.ll,' — ТУРБОНАДДУВ'

*,' NK-1.6 PK-PO*2 PA-0.9*PK PR-0.8*PK') 591 FORMATU2,' — ТУРБОНАДДУВА нет')

60 FORMAT(12,' — V-образный двигатель\5Х,Р5.1,' — угол развала')

FORMAT(12,' — рядный двигатель',5X,F4.1,' — угол развала') 102 FORMAT (1X.9F8.2)

105 FORMATC MEmax-',F8.1,' MKTE-',F8.1,' Кприспос- \F7.4) 111 FORMATC Какой у Вас двигатель: '

* /' 1 — дизель вихрекамерный СМД-14, Д-50'

* /' 2 — дизель с неразделенными камерами'

* /' 3 — карбюраторный —: ',$)

151 FORMATC ТА ТК РК PA PR JR ETAV '

* /1X,2F6.1,5F7.4)

152 FORMAT (8X, 'СКОРОСТНАЯ ХАРАКТЕРИСТИКА'/ *' N NEX МЕХ GEX GTX ')

FORMAT(1X,F6.0,F6.1,F7.1,F8.1,F5.1)

FORMAT(8X,'Индикаторная диаграмма 7

*' I UGOL Ход_П PVPD(I) PSD(I) PRD(I) PVXD(I) — МПа', *' AO-',F4.2 /' ГРАД мм. ВПУСК СЖАТИЕ РАСШИР ВЫПУСК' *,' AB-SP-',F7.2)

FORMAT(///36X,F3.1)

FORMAT(I3,F7.0,f7.2,4F7.3)

Продолжение прил. 15

157 FORMAT(8X,"CH.ibi инерции и нагрузки на шейки к/в 7 *' mj Мнгш Микш СС РСшат.ш РСкш ' */3F8.3,3F9.1)

FORMATC Ruimcp Runu48O RKincp Rkui480 74F9.1)

159 FORMAT(8X,' КРУТЯЩИЙ МОМЕНТ'/' MKSR=',F8.2,' MKTE=',F8.2)

FORMAT(8X,'Исходные данные для расчета ДВС) 162 FORMAT(8X,' 0 Турбонаддува нет ')

171 FORMAT(/' 1С SP DC LAMK MP MS E N '

*,' NEN — введите данные')

181 FORMAT(///' Курсовая работа по ТЕОРИИ ДВИГАТЕЛЯ VER 3.4 1998 г.' */' 1С — число цилиндров, '

*/' SP.DC — ход поршня, диаметр цилиндра, мм; ' */' LAMk — отношение радиуса кривошипа к длине шатуна; ' */' MP.MS — масса поршня, масса шатуна, кг; ' */' E.N.NEN — степень сжатия, об/мин, мощность, кВт ') CV - Ввод исходных данных ——————

OPEN(4,FILE-'REZULT') write(*,181) wrile(M71)

wrile(M7) IC,SP,DC,LAMK,MP,MS,E,N,NEN READ(*,*) IC,SP,DC,LAMK,MP,MS,E,N,NEN SPDC-SP/DC Cv - дизель-бензин и выбор исходных данных 'I I I I I I I I

JD-1!ДИЗЕЛЬ

IF(E.GE13) GOTO 46 lF(E.LE.ll) GOTO 47 write(*,*) ' ДАЛЬНЕЙШИЙ РАСЧЕТ НЕВОЗМОЖЕН'

GOTO 200

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

C-0.854 H=0.144 Ot-0.0 QN-44000 TR-1000. L0m-14.96 46 CONTINUE

WRITE (4,160) IF(E.GE.13) WRITE(4,18) IF(E.LE.ll) WRITE(4,19)

WRITE(4,17) IC,SP,DC,LAMK,MP,MS,E,N,NEN CV ****** ВВОД ИСХОДНЫХ ДАННЫХ ****** 402 IRADN-2

Продолжение прил. 15

Ugraz-60.

WRITE (*,3)

READ (*,*) IRADN

if(iradn.eq.2) WRITE(*,4)

if(iradn.eq.2) READ (*,*) ugraz

If(iradn.eq.l) WRITE(4,61) IRADN.UGRAZ

If(iradn.eq.2) WRITE(4,60) IRADN,UGRAZ

iturbo-1 WRITE (*,5) READ (*,*) ITURBO if(iturbo.eq.l) write(4,59) ITURBO if(iturbo.eq.O) write(4,591) ITURBO JK-2

WRITE(*,1U) READ(*,*) JK

WRITE(V) ' Двигатель — Jk- ',JK WRITE(4,*) ' Двигатель — Jk- ',JK IFUK.EQ.3.AND. JD.EQ.l) GO TO 200 IF(JK.LE.2..AND. JD.EQ.2) GO TO 200 write(*,158) write (*,58)

read<*,*) ALFA,DT,KDPA,PSI,LAMD,RO write(*,56) write (4,56)

WRITE(4,2) ALFA,DT,KDI>A,PSI,LAMD,RO,L0m WRITE(4,2) C,H,Ot,QN,ETAV,PO,TO WRITE(*,2) ALFA,DT,KDPA,PSI,LAMD,RO,L0m WRITE(*,2) C,H,Ot,QN,ETAV,PO,TO WRITE(4,2) PR.TR.Nk WRITE(?,2) PR.TR.Nk

_cv ****** НАЧАЛО ТЕПЛОВОГО РАСЧЕТА ДВИГАТЕЛЯ ****** ТК-ТО РК-РО

PA-PK-KDPA*PO PR-1.1*PO

IFUTURBO.NE.l) GO TO 203! *** TURBO ДО 203 " PK-PO*2

TK-TO*(PK/PO)"((NK-1)/NK) PA-0.9*PK PR-0.8*PK 203 JR-((TK+DT)/TR)*(PR/(E*PA-PR))

> Продолжение прил. 15

TA-(TK+DT+JR*TR) / (1+JR)

ETAV-(E/(E-1))*(PA/PK)*(TK/TA)*(1/(1+JR))

WRITE(4,15I) TA,TK,PK,PA,PR,JR,ETAV

WRITEC.151) TA,TK,PK,PA,PR,JR,ETAV

NI-1.41-110/N

N2-1.22+130/N

IFUD.EQ.l) N2-N2-O.O1! ДИЗЕЛЬ

PC-PA*E**N1

TC-TA*E**(NI-1) CV - Термохимический расчет 1.4.1 —=———■,

LO-(C/12+H/4-Ot/32)/0.21! кмоль

L0MR-L0*28.96

ml-ALFA*L0

MR-JR*ml

MC-ml+MR

M2-ml+H/4+Ot/32

IFOD.EQ.2) M2-M2+O.21*W*(1-ALFA)! BENZIN

DM-M2-M1

B0-M2/M1

MZrM2+MR

B-MZ/MC CV — Термодинамический расчет 1.4.2.

IFOD.EQ.2) GO TO 212! БЕНЗИН

AA-B*(13.8/ALFA+15.5)/10000! ДИЗЕЛЬ

BB-B* (8.314+(20.096+0.92/ALFA))

CCC—((20.16+1.74*TC/1000+8.314*LAMD)*TC+PSI*QN/

* (ALFA*LO«(1+JR)))

TZ-(-BB+SQRT(BB*BB-4*AA*CCC))/(2*AA) PZ-LAMD*PC

RO-B*TZ/ (LAMD*TC)

GO TO 213

212 AA-B*(15.491+13.716*ALFA)/10000!!! BENZIN

BB-B*(18.4+2.6*ALFA)!!!

CCC=-((20.16+1.74*TC/1000)*TC+PSI*(QN-120*1000*LO*

* (1-ALFA))/(ALFA*LO*(1+JR)))!!!

TZ-(-BB+SQRT(BB*BB-4*AA*CCC))/(2*AA)

PZ=0.85*B*PC*TZ/TC!!! БЕНЗИН

RO-1!!!

Cv - — ПАРАМЕТРЫ ГАЗА В КОНЦЕ РАСШИРЕНИЯ 1.5--------____

DEL-E/RO

PB-PZ/DEL**N2

Продолжение прил. 15

. TB=TZ/DEL**(N2-1) WRITE (4,11) WRITE(Vl)

WRITE(4,2) M,N2,PC,TC,LO,M1,MR WRITE<*,2) N1,N2,PC,TC,I.O,M1,MR WRITE(4,2) MC,M2,DM,B0,MZ,B,RO WRITE(*,2) MC,M2,DM,B0,MZ,B,RO WRITE(4,2) TZ,PZ,LAMD,DEL,PB,TB,LOMR WRITE(*,2) TZ,PZ,LAMD,DEL,PB,TB,LOMR write (4,*) ' ' write(*,*) ' ' iprod-1

wite(',57)

read(*,*) iprod if(iprod.eq.O) go to 200 if(iprod.eq.2) go to 402

CV ********** ПОСТРОЕНИЕ ИНДИКАТОРНОЙ ДИАГРАММЫ р—V 1.6> AB-SP

AO-AB/(E-1) VA-AO*E VB-AO*E VZ-AO*RO R-AB/2 RL2-R*LAMK/2

CALL GRBEG('MOT01.cgm',9,300.,250.,2) X1C-X1+AO+AB/2 Y1C-Y1-15.

CALL PXY(X1C+RL2,Y1C,COK1) CALL CXYR(X1C,Y1C,AB/2,OKR) XA-X1+AO

CALL MOVE(X1C+AB/2,Y1C,0) CALL MOVE(X1+AO,Y1C,1) CALL MOVE(X1C-AB/2,Y1C,0) CALL ARCIA(R,180.,359.)

YMAX=PZ*1.05

CONSTS=PA*VA**N1

CONSTR-PB*VB**N2 VXD(1)-AO VXD(13)=AB+AO

CV ****** ДДЯ 4 MKR И ИНДИКАТОРНОЙ ДИАГРАММЫ CV ** p—V свернутая и p—FI развернутая по методу проф. Крикса 4.1

WRITE(4,154) A0.AB

' Продолжение прил. 15

DO 23 1-1,13

UGOL-(I-1)*15.

IFU.EQ.l.OR.i.eq.l3) GOTO 124

CALL LPAX(COK1,UGOI.,PRM1)

CALL PILC(PRM1,OKR,J>

CALL SORTP(4,TOCH1)

VXD(I)-TOCH1 (П-ХЛ+АО

CALL MOVE(COK1(1).COK1(2),0)

CALL MOVE(TOCH1(1),TOCH1(2),1)

CALL MOVE(TOCHI(1),Y1,1) 124 VXDN1-VXD(I)**N1

VXDN2=VXD(1)**N2

PVPD(I)-PA

PSD(I)-CONSTS/VXDN^I1! Сжатие

PRD(I)-CONSTR/VXD\2! Расширение

IF(VXDd).LT.VZ) PRD(I)-PZ

PVXD(I)-PR

XP1-VXD(I)-AO

WRITE(4,156) I,L'GOL,XP1,PVPD(I),PSD(I),PRD(I),PVXD(I) IF(I.EQ.13) GOTO 23

CALL MOVE(VXD(I)+X1,Y1,0)

YYYY=PRD (I) «VR2/YMAX

CALL MOVE(VXD(I)+X1,Y1+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

CALL SYMBOL(X1,Y1+VR2+5.,10., * 'Индикаторная диаграмма р—V*,28,0.) VR1-VXD(13)*1.O5

CALL LGR(X1,Y1,VR1,VR2,O.,VR1,O.,YK,VXD,PVPD,13)!впуск

CALL SETPEN(2)

CALL LGR(X1,Y1,VR1,VR2,O.,VR1,O.,YK,VXD,PSD,13)!сжатие

CALL SETPEN(4)

CALL LGR(X1,Y1,VR1,VR2,O.,VR1,O.,YK,VXD,PRD,13) расширение

CALL SETPENO)

CALL LGR(X1,Y1,VR1,VR2,O.,VR1,O.,YK,VXD,PVXD,13) 'выхлоп

CALL SETPEN(l)

CALL LGO(X1,Y1,VR1.VR2,0.,VR1,0,YK,40.,1)

CALL GREND

Продолжение прил. 15

snd-(prd(l)-PSD(3))/30 WRITE(4,504) prd(l),psd(3),snd 504 FORMAT С Скорость нарастания давления (prd(l)-psd(3))/30 - snd

*/,3f9.5)

CV ******«**. РАСЧЕТ ПОКАЗАТЕЛЕЙ РАБОЧЕГО ЦИКЛА 1.7 ************ IF(JD.EQ. 1) РИ-<РС/ (Е-1))*(LAMD*(RO-1)+LAMD*RO*

* (1-1/DEL*«(N2-1))/(N2-1)-(1-1/E**(N1-I))/(N1-1))

IF(JD.EQ2) PI1-(PC/(E-1))*

* (LAMD*(1-1/E**(N2-1))/(N2-1)-(I-1/E**(N1-1))/(N1-1)) CV PIND-PI 1*0.94-(PR-PA)

PIND-PI1*O94 SR-0 SS-0

DO 101 1-2,13

SS-SS+(PSD(I-1)+PSD(I)) * (VXD(l)-VXD(I-D) /2 101 SR-SR+(PRD(I-1)+PRD(I)) * (VXD(I)-VXD(I-D) /2 PLOS-SR-SS PINDGR-PLOS/SP

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

IFOK.EQ.3) PM-0.035+5*N/100000

IFUK.EQ.2) PM=0.067+5.2*N/100000

IFUK.EQ.l) PM-0.092+1.01*N/10000

PE-PIND-PM

ETAM-PE/PIND

ROO-PO*1000000/(287.*TO)

ETAI-PIND*L0m*ALFA/(ROO*ETAV*QN/1000)

ETAE=ETAI*ETAM

GI-36OO/(ETAI*QN/1OOO)

GE-3600/(ETAE*QN/1000)

GT-GE*NEN/1000 CV ********** РАСЧЕТ ОСНОВНЫХ РАЗМЕРОВ ДВИГАТЕЛЯ 1.7.3 **********

VL-30*4*NEN/<PE*N)

VH-VL/IC

DCH00*(4*VH/(PI*SPDC))**(l 73.)

SP-DC*SPDC

R-SP/2

VL-PI*DC*DC*SP*IC/4000000

NE-PE*VL*N/(30*4)

ME-9550*NE/N

Продолжение прил. 15

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

NP-NE/(IC*FP)*100*100 WR1TE(4,14)

WRITE(4,2) РП,PIND,VPS,PM,PE,ETAM,ROO WRITE(4,2) ETAI,ETAE,GE,GT,VL,VH,DC WRITE(4,2) NE,ME,SP,R,FP,NP,NL

CV ********** КИНЕМАТИЧЕСКИЙ РАСЧЕТ З ************ W-PI*N/30 DO 21 1-1,49 FI(I)-(I-1)*15. F1R-FI(I)*RAD

XP(I)-R*((l-COS(FIR))+(l-COS(2*FIR))*LAMK/4) Vl'(l)-W*R*(SIN(FIR)+LAMK*SIN(2*FIR)/2)/1000 JP(I)-W*W*R* (COS (FIR)+LAMK*COS (2*FIR)) /1000 21 CONTINUE

YK-W*R/1000

YY-W*W*R/1000

CALL GRBEG('MOT09.cgm',9,300.,250.,2) CALL SYMBOL(Xl,Y2+VR2+5.,9.,Трафики Xn,Vn,Jn',16,0.) CALL LGR(X1,Y2,VR3,VR2,O.,72O.,O.;2*R*1.1, FI.XP.49) ' CALL VAX IS (-60.,'Xn ',3,40.,4,0,1)

CALL SETPEN(2)

CALL LGR(Xl,Y2,VR3,VR2,0.,720,-YK,YK, FI.VP.49) CALL LGO(X1,Y2,VR3,VR2,0.,720.,-YK,YK, 90..5.)

call setpen(4)

CALL LGR(X1,Y2,VR3,VR2,O.,72O.,-2*YY,2*YY,F1,JP,49) CALL YAXlS(-120.,'Jn ',3,1000.,5,0,D

CALL GREND

CV ********** ДИНАМИЧЕСКИЙ РАСЧЕТ 4 ************ CV ********** СИЛЫ ИНЕРЦИИ (КОНСТАНТЫ) 4.1 ************ W-PI*N/30 MJ-MP+O.275*MS WRl-W*W*R/1000 CC—MJ*WR1

MNGST-0.725*MS

PCS—MNGST*WR1

MIKS-423.3*FP/(1000*1000)+MNGST

PCK—MIKS*WR1 CV *** СИЛЫ ИНЕРЦИИ, СУММАРНЫЕ, КАСАТЕЛЬНЫЕ, ТАНГЕНЦИАЛЬНЫЕ ****

Продолжение прия- 15

CV ***** РЕАКЦИИ МЛ ШАТУННОЙ И КОРЕННОЙ ШЕЙКАХ 4.3 и 4.4 *** CV ****** СИЛА PRP ИЗ ИНДИКАТОРНОЙ ДИАГРАММЫ *******

DO 24 1-1.12 PRP(I+1)=PVPD(I+1)*FP PRP(I+13)=PSD(13-I)*FP PRP(r+25)=PRD(I+l)*FP PRP (1+37)-P VXD (13-D *FP 24 CONTINUE

PRP(1)=PR*FP PRP(25)-PZ*(l-0.05)*FP PRP(26)-PZ*FP PRP(37)-(PA+PB)*FP/2 S1S-0 S1K-0

DO 25 1-1,49 FI(I)-(I-1)*15. FIR=FI(I)*RAD

PJ(I)-CO(COS(FIR)+LAMK*COS(2*FIR)) PSUM(I)-PRP(I)+PJ(I) BETfASIN (LAMK*SIN (FIR)) COS (BET)

T(I)PSUM(I)«SIN<FIR+BET)/ KS (I)-PSUM(l) 'COS (FIR+BET) /COS (BET) RSSO)-SQRT((KS(I)+PCS)**2+T(I)**2/ RKS(I)-SQRT((KS(I)+PCK)**2+T(I)**2) S1S-S1S+RSS(I) S1K-S1K+RKSU) 25 CONTINUE

RSSR-S1S/49

RKSR-S1K/49

WRITE(4,157) MJ,MNGST,MIKS,cc,pcs,PCK

WRITE(4,158) RSSR,RSS(32),RKSR,RKS(32)

CV ----— ПОСТРОЕНИЕ ГРАФИКОВ СИЛ И РЕАКЦИЙ 4.6 —-----------

call minmax(Pj,49,YN2,YKl) call minmax(PRP,49,YNl,YKl) YK-YK1*1X)5 YN-0.-YK/3

CALL GRBEG('MOT02.cgm',9,300._>250.,2) CALL SYMBOL(X2,Y2+VR2+5.,10., *'Диафамма сил Prp, Pj, Psum',27,0.)

CALL LGO(X2,Y2-VR2/3,VR3,4*VR2/3,0.,720.,YN,YK,90.,20000.) CALL LGR(X2,Y2-VR2/3.VR3,4*VR2/3,0.,720.,YN,YK,FI,PRP,49) CALL SETPEN(2)

Продолжение прия. 15 CALL LGR(X2,Y2-VR2/3,VR3,4*VR2/3,0.,720.,YN,YK,FI,PJ,49)

CALL SETPE.N44)

CALL LGR(X2,Y2-VR2/3,VR3,4*VR2/3,0.,720.,YN,YK,FI,PSUM,49) CALL GREXD

CALL GRBEG('MOT03.cgm',9,300.,250.,2) CALL SYMBOL(x2,Y3+VR2+5.,10., *Трафик сил на шатунной шейке KS, T ',39,0.)

CALL SETPEN42) CALL LGR(X2,Y3-VR2/3,YR3,4*VR2/3,O.,72O.,YN,YK,FI,KS,49)

CALL SETPEN(l)

CALL LGR(X2,Y3-VR2/3,VR3,4*VR2/3,0.;720.,YN,YK,FI,T,49) CALL LGO(X2,Y3~VR2/3,VR3,4*VR2/3i0.,720.,YN,Yk,90.,20000.) CALL GREND

CALL GRBEG('MOT04.cgm',9,300.,250.,2) CVV TMP-VR2*4/3

TMP-200.

CALL SYMBOL(Xl,Y6+TMP+5.,6.,

"Полярная диаграмма нагрузки на шатунную шейку',48,0.) CALL LGR(X1,Y6,TMP,TMP,YN,YK,YK,YN,T,KS,49) CALL LGO(X1,Y6,TMP,TMP,YN,YK,YK,YX,20000.,20000.) CALL TMF(O.,-PCS,X11,Y11) CALL TMF(0.,-PCK,X12,Y12) CALL TMF(T(32),KS(32),X13,Y13) CALL SYMBOL(Xll+5.,Yll,7.,'Or,2,0.) CALL SYMBOL(X12+5.,Y12,7.,^>O2-,2,0.) CALL MOVE(X11,Y11,0) CALL MOVE(X13,Y13,1) CALL MOVE(X12,Y12,1) CALL MOVE(X12-20.,Y12,0) CALL MOVE(X12+20.,Y12,I) CALL MOVE(X11-20.,YI 1,0) CALL MOVE(X11+20.,Y11,1) CALL GREND

CALL GRBEG('MOT05.cgm',9,300.,250.,2) YN-0.

CALL SYMBOL(x2,Y4+VR2+5.,9., *'Рсакции на шатунной и коренной шейках RSS и RKS',47,0.)

CALL SETPEN(2)

CALL LGR(X2,Y4,VR3,VR2,0.,720.,YN,YK,FI,RSS,49) CALL TMF(0., RSSR,X11,Y11) CALL TMF(720.,RSSR,X12,Y12) CALL MOVE(XI 1,Y11,0)

Продолжение прил. 15

CALL MOVE(X12,Y12,1)

CALL SYMBOL(X12+l.,Y12+2.,7.,'Riuuicp',5,0.)

CALL SETPEN44)