Modelling and Control of Diesel Engine Modelling and Control of Diesel Engine
Modelling and Control of Diesel Engine Modelling and Control of Diesel Engine
Modelling and Control of Diesel Engine Modelling and Control of Diesel Engine
Items in Figshare are protected by copyright, with all rights reserved, unless otherwise indicated.
PUBLISHER
© M.A. Eissa
LICENCE
CC BY-NC-ND 4.0
REPOSITORY RECORD
Eissa, M.A.. 2019. “Modelling and Control of Diesel Engine”. figshare. https://hdl.handle.net/2134/10318.
This item was submitted to Loughborough University as a PhD thesis by the
author and is made available in the Institutional Repository
(https://dspace.lboro.ac.uk/) under the following Creative Commons Licence
conditions.
- ?Ht~R 1996
'-~5-
17 NOV 1992
1 f DEC'1992
,/
I
000 7421 02
~I~I I I I I I ~I I I I ~I I I I~I I I ~I I ~
MODELLING AND CONTROL OF DIESEL ENGINE
by
M. A. EISSA B.SC.
A Doctoral Thesis
MAY 1984
PAGE
~cRnowledgements 1
Synopsis 2
Chapte,Jr 1 3
:tntroduction 3
Chapter 2 6
List ef symliols 6
2.1 Introduction 14
conditions
multiple plot
each cycle
diesel engine
theary.
Appendtx (3 [ 256
microprocessor.
Engineering for the experimental work which was carried out in the mechanical
Also I wish to thank Mr. S.K. Sinha of the G,E.C. for providing
Loughborough University who have made the four years' work so enjoyable.
1
SYNOPSIS
which takes into account in cycle calculations the input temperature and
ture and pressure, and gives as outputs the average torque and speed. A
typical heat release curve from a practical engine test is used with
ignition delay considered constant in time. Results from the program are
matrix relating input manifold pressure, mass of fuel injected, and angle
and zero placement techniques, the system inputs and outputs being demanded
and achieved torque and speed .•
2
Chapter 1 Introduction
is relatively recent (14), and many models are based on the emptying
3
based on the relation between the injection rate and heat release
(41,69,82,84).
engine, the dual combustion engine cycle is linearized, with the help
of the non. linear digital computer model, and cast in the form of
and parameters must be considered constant which are not truly so,
4
of control, so the contribution in this thesis is the utilization of
given its original transfer function matrix (G) and unknown feedback
~trb(H).
This approach does not use matrix inversion or devision.
lacking that, a program has been written to help produce the algebraic
basic rules laid down to ease the design procedure (17,52,68). The
5
CHAPTER 2
LIST OF SYMBOLS
feedback controller. (0 or 1)
AREAI elements of area under the heat release curve (kg fuel.deg)
2
AREAMO moment of area of heat release curve (kg fuel.deg )
6
CVAIR specific heat of air as function of the (J/kg.k )
o
ambient temp.
combustion zone
per step
DWORK work done in one step (bar.lt.)
ECC cooling coefficient of engine block (kg.J/kg.k o )
EETAV volumetric efficiency of the engine, defined as
7
EETAVO volumetric efficiency at idling speed (used
compressor
8
IHP indicated horse power of the engine (hp)
nearest integer
changes from (0 + I)
each step
of the cylinder
of the engine
9
MWGAS mean molecular weight of combustion gases (gm/mol. )
MWFUEL. molecular weight of diesel fuel (gm/mol.)
MXCMLT max. value of combustion length of max. engine (deg)
10
P5EXH average value of exhaust pressure over (bar)
exhaust stroke
the T.O.C.
input temp.
11
TS gases temp. at the beginning of exhaust stroke (ko)
stroke
in compression stroke
in expansion stroke
function
in each step
12
WORKC amount of work done in compression stroke (bar.lt.)
WORKEX amount of work done in expansion stroke (bar.lt. )
WORKI amount of work done each step (bar.lt.)
X current position in the heat release (deg. )
13
2.1 Introduction
suitable for control studies are few; every diesel engine manufacturer
has a number of them but for reasons of industrial secrecy they are
rarely published.
" . temperature
. injection advance angle
.. mass of fuel
the fuel as CSH IS ' which is petrol, while the constant pressure combustion
14
period and is obtained from an empirical formula which is rather
of 1 Bar is· quoted though in fact any value will give the same
results:
inadequate.
15
the cycle, a realistic heat release curve is used
at each step.
the next page and the digital model is explained line by line after
16
( nUT
"1
2.2.2 DOJl! sP!CIruno';I
DATA AIIO Jlrnu. COlfDrfIOIIS
, -I CtCLr; CALCIIUTIOIfS
I 1:. 2•4• 2
I CALC'OLAtIon
COMnUSOII rovt
. J
I
I
1:. 2 .,.) MODtLLIIC or:1
IIJICTIO. ADVA'Ct rll~r
I
I e·2 .,., CALC'OUTto. or .sl
I!lUU SUOK! III 110 snrs
I j
I 2.2.'. ,
co~rnsstOI
CALCOUTIO. or
snOK! 11 110 snP'S.
I
J!IIIJ!:CTtOJf' COIDttlO.S.
CO~JUSTtO' UJUI!~ CO.STAllt YOLII'll.
VO!tK DOl! U CO!'lrlllSSIOII snOK!.
2.2.'.6 CIoLC'Ou,nOIl or
1
I
VOIII DOll III mlollStOIl snon.
lOll 1'0 DtTlCT Mu-conUStIO!l '~!s.s'OIII.
O"Io'OST "Ato' oronc.
~
.
2.2.'.1 CIoLCUUTIOI 0' IIIIIoOS't1
S'tltOIl u 110 snps.
~
1,2.2., •• CloLCU~Io'ttO' or T!'III
110111 DOl! rE!!: ClCLI.
1
2.2.'.9 101f TO CLlU IoUAlS TO
OSlO IoCUI.
III
l2.2.,.lD !lIlloUSt ploIIAlfnns JUDItII
Till tvlllIJI •
2.2.4.U
~
CALcvuno. OP PCII!
OOrPUT PlllrDn.UCt ,
I
J
.1 • .1.,.12 CALCIILAfIO. 01 TUIIII'~
oo'tPlfT ram •
Gb
(2.1) Block diagram of the digital computer model.
17
2.2 The digital computer program.
The main purpose for which the model has been used is to investigate
the engine such as specific heat, and to check any assumptions made when
(1.1) •
explanation.
18
REAL FUEL (300) ,FUELl (300) ,FUEL2 (300) ,FUEL3 (300)
REAL·BEATRL(100),IBP(3oo),IMEP,INJLTH,IPTMAS,Ml,MlMAX
REAL MB20,MININJ,MN2,M02,M02BRN,MS,MSMAX,~IR,MWGAS,MWFUEL
REAL MXIJLT,Nl,N2,NE,NEMAX,NEMIN
REAL TCY (800) ,TORQ (800) ,TORQl (800) ,TORQ2 (800) ,TORQ3 (800)
Explanation :
The real var:ta.I>les and arrays are def:tned :tn the beg:tnning of
EPSLUM=14.4
NCLNDR-6
BORE-<). 1 27 M
STROK-D.l 49 M
STROKS-4
CMCHLT-STROK*«EPSLUM+l)/EPSLUM) M
PA-l.03 BAR
TA=288 K
19
C THE MAIN INPUT VARIABLES
Pl=1.25 BAR
RACK=l.
MFMAX=0.OOO7 KG/CYCL
C ANGLIJ=-30 DEG
T1IPUT=300. K
R=8.2057E-2
C NE=ENGINE SPEED
NEMAX=2000.0 RPM
NEMIN=500. RPM
NE=500 RPM
._ PREVNE=500 RPM
NPT=4
CWRITE=lO
NSTP =180
EVOBDC=-50 DEG
MXIJLT=24 DEG
ADVMCH=l.
TURPWR=4 H.P.
ETATUR=O.25
20
STROK length of stroke (m)
STROKS number of strokes in cycle
CMCHLT combustion chamber length (m)
PA ambient pressure (Bar)
TA ambient temp. (Ko )
PI compressor output pressure (Bar)
RACK fuel control rack (0 ? 1)
MFMAX max. amount of fuel per cycle (kg)
ANGLIJ injection advance angle (deg.)
TlIPUT after cooler output temp. or engine input temp. (K)
R
universal gas constant. (LT.Bar/deg.mol.)
NE Engine speed (r.p.m)
NEMAX max. engine speed (r.p.m.)
NEMIN min. engine speed (r.p.m.)
PREVNE engine speed at previous cycle (r.p.m.)
NPT
number of·points (or cycles) need to be calculated.
CWRITE=lO
parameters used to write out the engine variables
C THE DIESEL FUEL CONTAINS (13 OR 14) CARBON ATOMS IN THE MOLECULE
MWFUEL=198. GM/MOL
EXES=l. 7
AFR=EXES*(21.S(16.*2.)+4.*21.S*(14.*2»/198.0 KG/KG
V2=VlMAX/EPSLUM LT
VSMAX=VlMAX-V2
C .H + 21.S O 414 CO + IS H 0
14 30 2 2 2
198 + 21.S(16x2) 414(12+16x2) + lS(2+16)
22
The theoretical amount of oxygen needed for one mol of fuel is
In fact the real amount of air needed for complete combustion is much more
than the theoretical amount, So the excess air factor was chosen to be
EXES=l. 7
Now the real amount of air needed for one mol can be calculated as
In the program
ADVCON=O.O DEG
DO 7 J=l,NPT
SPOIT=pl*MNAIR*lE-3/ CR*TlIPUT)
EETAVo=0.90-Ce.E-8)*CNEMIN**2» KG/CYCL
M1MAX=SPCWT*VIMAX*EETAVO KG/<-'YCL
23
MSMAX=SPCWT*VSMAX*EETAVO KG/CYCL
C MFMAXcMSMAX/AFR KG
TRAMAS.=SPCWT*V2*EETAVO
C MF=MUST BE GIVEN
WRITE(l,5)
5 FORMAT (lHl,/////
WRITE (1,5)
WRITE (1,200)
WRITE (1, 5)
MF:RACK*MFMAX KG/CYCL
Explanation.
(ADVCON), input fuel (RACK) and air pressure (Pl). (See the end of the
loop 7-continuel.
has an optimum value depending mainly upon input fuel and output speed.
So the optimum advance is obtained for different values of speed and fuel by
plotting the engine output power for different values of speed and certain
24
TRAMAS the trapped mass of gas in the cylinder Ckg)
MF
is the variable mass of fuel per cycle . Ckg)
as RACK changed from Co ~ 1)
Ml=SPCWT*VlMAX*EETAV
KG/CYCL
MS=SPCWT*VSMAX*EETAV
KG/CYCL
C DISASSOCIATION OF INPUT GASES
C INPUT MASS OF 02
MC02=o. 0
MH20=0. 0
M02BRN=O. 0
TI=T1IPUT
K
PI=Pl*EETAV
BAR
CVAIR=ClOOO/29.0)*C20.75S+0.00224S*(TA_273.)) J/KG. K
CVGAS=CVAIR
po 10 1=1, 3
NWRITE=CWRITE
WHITEU, 5)
25
210 FORMAT(lHl,14BENGlNE SPEED= ,F7. 2,2X,6HR. P. H , 10HAND RACK= ,F4. 2)
AIRFUL=HS/MF
250 FORMAT(lHO, 26HMASS OF AIR PER CYCLE Hl= ,F10. 8,3HKG ,10X
Explanation:
to the speed
26
MC02 initial mass of carbon dioxide it is
considered = 0 (KG)
equal to CVAIR
certain conditions
27
WRITE (l,225) the current· value of angle of injection
C COMPRESSOR POWER
ETACSR=ETATUR
GAMAC=1.4
CSRPWR=TURPWR
Pl=PA*(l+(736.*CSRPWR*ETACSR/(AIRMAS(I)*CVAIR*TA»)**
* (GAMAC/(GAMAC-l» BAR
WRITE(l,265).PL.
WRITE (l,285) TC
CSRPWR=«AIRMAS(I)*CVAIR *TA)/(ETACSR*736.»*
«(Pl/PA)**«GAMAC-l)/GAMAC1-l) H.P
WRITE (1,2901 CSRPWR
28
C TURBINE DYNAMIC EQUATION IS
ETACOL--Q.5
TlIPUT=TC-ETACOL*(TC-TA) K
WRITE{1,295) TlIPUT
, lX, 6HTEMP. , IX, 51lMWGAS , lX,8HSP. HEAT, 3X, 5HliORK , 5X, 5HFUEL )
WRITE (1,305)
- - - - - - - - - - - - - - - - - - ______ 1
Explanation:
conditions.
29
CSRPWR compressor power as function to air mass
of the compressor
30
C ALSO THE END OF COMBUSTION COMES A BIT LATE SO
CURVl=O. 0
DO 150 L= 1 , IRF
X=L*TIMSTP
UNIT
CURV2=(X**0.25) *EXP (-0.13*ABS (X-6. 5»
AREAI(L)= (CURV2+CURVl)/2. 0
CURVl=CURV2
AREAFN=AREAFN+AREAI (Ll
150 CONTINUE
Explanation:
31
INJLTH fuel injection period length as function
T.O.C. (OEG)
CMBLTH the total period of. combustion (see fig. (2.1) (DEG)
I
CMBOLY PRETDC PS'ITDC
_I
combustion
. begins ..... !
~ ,
'///// /////////
~,," .'-""'.'
/ '.LL
injection
/ INJL--'l'H I
CMBOLY CMBEND
begins I
- ANGLIJ
T. .C.
Fig. (2.l)
TIMSTP The variable TIMSTP describes the increase
release curve
32
:--:it
....-.~t
. --:-
.;._. -::
,--
33
AREAI the average area of each fuel sample
WRITE (1,315)
DEGRE=O
DEG
WORKI=O.O
VI=V2
LT
PRVVCY=V2
LT
SPCWT =Pl*MWAIR*lE-3/(R*TlIPUT} . KG/LT
AREAMO=O.O
KG.DEG.2
DDEGRE=O.O
DEGRE=DEGRE+DDEGRE
DEG
DDEGRE= (180. O/FI.oAT (NSTP JJ DEG
PHI=(DEGREl*3.14/180.0) RADIAN
VCY(K}=V2+«VlMAX-V2l/2.0l*
PRVVCY=VCY (Kl
IPTMAS=DVCY*SPOwr*EETAV
KG/STP
34
TRAMAS=TRAMAS+ABS (IPTMAS) KG
DM02=IPTMAS*16./(14.*4.+l6.) KG/STP
M02=M02+DM02 KG
DMN2=IPTMAS*(14.*4.)/(14.*4.+l6.) KG/STP
MN2=MN2+DMN2 KG
TCY(K)=(T1IPUT*IPTMAS*CVAIR+TI*TRAMAS*CVGAS)/(IPTMAS*CVAIR+TRAMAS*
*CVGAS) K
MWGAS=(M02*32.+MN2*28.+MC02*44.+MH20*18.)/(M02+MN2+MC02+MH20) . G/MOL
PCY(K)=IPTMAS*R*TCY(K)/(VCY(K)*MWAIR*lE-3 ) +
+TRAMAS*R*TCY(K)/(VCY(K)*MWGAS*lE-3 BAR
TI=TCY(K)
PI=PCY(K) BAR
VI=VCY(K) LT
CV02=(lOOO/32.0)*(20.959+o.oo295*(TCY(K)-273.» J/KG.K
CVN2=(lOOO/28.0)*(20.704+o.OO211*(TCY(K)-273.» J/KG.K
CVC02=(lOOO/44.0) * (27.545+0.00922* (TCY (K)-273.» J/KG.K
CVH20=(lOOO/18.0) * (25.184+o.ooS2S*(TCY(K)-273.» J/KG.K
CVGAS=(M02*CV02+MN2*CVN2+MH20*CVH2O+MC02*CVC02)/
/(M02+MN2+MH20+MC02) J/KG.K
CVGAS=«704.9+o.0708*TCY(K»*Ml+(1833.4+o.4931*TCY(K»*
*MAXMF)/(Ml+MAXMF) J/KG.K
C CVGAS=838. J/KG.K
IF (NWRITE. LT. CWRlTE) GO TO 300
WRITE (1, 3201 DEGRE, VCY (Kl,TRAMAS ,PCY (K) ,TCY (K) ,MWGAS,
320 FORMAT(lX,F4.0,2X,F8.S,2X,F8.7,lX,F7.2,lX,F8.1,2X,FS.2,lX,
,F7.2,lX,F7.2,2X,E12.6 J.
NWRlTE=O
35
300 CONTINUE
NWRlTE=NWRlTE+l
500 CONTINUE
Explanation:
To find the relation between volume of combustion chamber and crank angle.
B.D.C. V1MAl{
..•.. VCY(~l
R
..
A area
xo = R+r
T.D.C.
Fig. (Z.3)
36
PHI - crank anqle
- ABS (DEGRE*3. l4/l80.0) radian
From the drawing Fiq. (2. 3)
-----
sin a
sin ~
R
~
r
r
R
~l
----~-
Xo - X - R + r - R cos a - r cos ,
37
DVCY=VCY(K)-PRVVCY the change in the volume of the combustion
38
CV02, CVNN2, CVC02, CVH20 the linearized values of specific
the pure air and the mass of burned fuel, see Appendix _ (2 )
WRITE (1, 320) to write the current values of crank angle, cylinder
previous variables.
write (1,320)
MAXMF=o.o
liORKC=o.O
liORKEX=O.O
DEGCMB=O.O
39
WRITE (1,3231
WRITE (1,310)
DEGRE=DEGRE+DDEGRE DEG
DDEGRE=(180.0/FLOAT(NSTP » DEG
Explanation:
to the variables, max. mass of fuel, work done in compression stroke, work
PHI=ABS(DEGRE*3.l4/l80.0) RADIAN
• (4.333-COS(PHI)-3,3333*SQRT(1.-0.09*(SIN(PHI)**21» LT
C INTAKE VALVE CLOSURE
SPCWT=PI*MWGAS*lE-3/(R*TI) KG/LT
C DVCY=VCY(K)-PRVVCY
C PRVVC¥=VCY(K)
C IPTMAS=DVCY*SPQiT KG/STP
C TRAMAS=TRAMAS+IPTMAS KG
PCY (Xl =PI* (VI/VCY (K) I **Nl BAR
TC¥(KI=TI*(VI/VCY(KII**(Nl-ll K
'40
CV02=(1000/32.0)*(20.959+0.00295*(TCY(K)-273.)) J/KG.K
CVGAS=(M02*CV02+MN2*CVN2+MH20*CVH20+MC02*CVC02)/
/M02+MN2+MH20+MC02) J/KG.K
CVGAS=«704.9+0.0708*TCY(K))*M1+(1833.4+0.4931*TCY(K))*MAXMF)/
/(M1+MAXMF) J/KG.K
C CVGAS=838. J/KG.K
MWGAS=(M02*32.+MN2*28.+MC02*44.+MH20*l8.)/(M02+MN2+MC02+MH20) G/MOL
Explanation:
average molecular weight, for more detail see intake stroke. para (2.2.4.4)
C INJECTION CONDITIONS
41
Explanation:
NWRlTE=CWRlTE
DEGCMB=DEGCMB+DDEGRE
DEG
NCOMB=IFIX(DEGCMB+0.5)
APPROXMT
XCOMB=NCOMB*TIMSTP
BETFN2=(XCOMB**0.25)*EXP(-0.13*ABS(XCOMB_6.5»
BEATRL(NCOMB)=(BETFN2+BETFN1)*MF/(2.0*AREAFN) KG/STEP
C BEATRL(NCOMB)=«2./CMBLTB)-(2./CMBLTB**2.)*(OEGCMB_0.5»*MF KG/STEP
HETFNl=HETFN2
MFI(K)=BEATRL(NCOMB)
KG/STEP
AREAMO=AREAMO+BEATRL (NCOMB) *DEGCMB KG.OEG 2
TRAMAS=TRAMAS+MFI (Kl KG
MFCUR(KI=MFI(K)*1.E6/00EGRE MG/OEG
MAXMF=MAXMF+MFI(K) KG
CNTROO= (AREAMO/MAXMFI +PRETOC OEG
C SP. BEAT AT CONSTANT VOLUM
J/MOL. K
C OUTPUT MASS OF C02 IN TBIS STEP IS (SEE COMBUSTION EQUATION)
C02I=(14*44/l98.0) *MFI (K)
KG
MC02=MC02+C021
KG
C OUTPUT MASS OF H20 IN TBIS STEP IS (SEE COMBUSTION EQUATION)
H20I= (15*18/198.0) *MFI (K)
KG
MH20=MH20+H201
KG
CVGAS=(M02*CV02+MN2*CVN2+MH20*CVH20+MC02*CVC02)/
/(M02+MN2+MH20+MC021
J/KG.K
42
MWGAS=(M02*32.+MN2*28.+MC02*44.+MH20*18.)/(M02+MN2+MC02+MH20) G/MOL
TI2=TCY (K) K
PI2=PCY (K)
TCY(K)=TI2+«HEATRL(NCOMB)*CALV*ETAUSE)/(TRAMAS*CVGAS» K
PCY(K)=PI2*(TRAMAS+HEATRL(NCOMB»*TCY(K)/(TRAMAS*TI2) BAR
710 CONTINUE
Explanation:
combustion.
of injection pump)
release function.
43
MFI mass of fuel burned each step (KG/STP)
release curve
combustion (KG)
cycle (KG)
CVGAS the mean specific heat of the gas at current step (J/KG.K)
the lctw
. y-l
= constant , and TV = constant.
44
But during the combustion an external heat added to the system.
is the adiabatic changes, and second stage is the heat addition to the
Fig. (2.4)
and
V
P
2
= P (
1 V
-1:: )y where y is adiabatic exponent
2
V
1 y-l
= T ( -)
1 V
2
In stage (rI) the combustion heat added to the system increases the
T T + Jlheat
3 = 2 m .CV (the temp. increase due to heat release)
3
V2 T3m3
- --
V3 T2m2
=
45
710 The end of combustion calculation under constant volume.
Explanation:
IF (NP3MAX.EQ.l) GO TO 705
NP3MAX=NP3MAX+l
P3 (I) =PI
705 CONTINUE
NWRITE=O
327 CONTINUE
NWRITE=NWRITE+1
VI=VCY(K)
LT
TI=TCY(K)
K
PI=PCY(K)
BAR
700 CONTINUE
46
Explanation:
between two successive pressures and when this sign is changed then the
so it is passed over.
IF(PCY(K)-PI).GT.OI GO TO 705
pass over.
detection.
IF statement to control writing the variables
47
temp. of cycle, molecular weight of gases, specific heat, work done so far
TI initial temp
PI initial pressure
WRITEC1,329}
DEGT=O.O
00 900 K= 361,540
DDEGRE=(180.0/FLOAT(NSTP )} DEG
DEGRE=DEGRE+DDEGRE DEGRE
Explanation:
and C2.2.4.5}
48
2.2.4.6.1 Conditions of Next Step Without Combustion
VCY(K)=V2+«VlMAX-V2)/2.0)*
/(M02+MN2+MH2O+MC02) J/KG.K
CVGAS= «704. 9+0. 0708*TCY (K) f*Ml+ (1833.4+0. 4931*TCY (K)) *MAXMFJ /
/ (Ml+MAXMF) J/KG.K
C CVGAS=838 J/KG.K
Explanation:
..
IF (DEGCMB.GT. (CMBLTH-DDEGRE)) GO TO 910
NWRITE=CWRlTE
DEGCMB=DEGCMB+DDEGRE DEG
NCOMB=IFIX(DEGCMB+O.5) APROXMT.
XCOMB=NCOMB*TIMSTP
HETFN2=(XCOMB**0.25)*EXP(-0.13*ABS(XCOMB-6.5))
49
HEATRL (NCOMB) = (HETFN2+HETFNl) *MF/(2.0*AREAFN)
KG/STEP
C HEATRL(NCOMB)=«2./CMBLTH)-(2./CMBLTH**2.)*(DEGCMB_0.S»*MF KG/STEP
HETFNl=HETFN2
MFI(K)=HEATRL(NCOMB)
KG/STEP
AREAMO=AREAMO+HEATRL (NCOMB) *DEGCMB
KG.DEG 2
TRAMAS=TRAMAS+MFI (K)
KG
MFCUR(K)=MFI(K)*1.E6/DDEGRE
MG/DEG
MAXMF=MAXMF+MFI (K)
KG
CNTROD= (AREAMO/MAXMF) +PRETDC DEG
C SP. HEAT .AT CONSTANT VOLUM
J/MOL. K
C OUTPUT MASS OF C02 IN THIS STEP IS (SEE COMBUSTION EQUATION )
C02I=(14*44/198.0)*MFI(K)
KG
MC02=MC02+C02I
KG
C MASS OF 02 WHICH BURNED IN THIS STEP IS
M02BRN=(21.S*(2*16)/198.0}*HFI(K) KG
C THE RESIDUAL MASS OF 02 IS
M02=M02-M02BRN
KG
C OUTPUT MASS' OF H2o IN THIS STEP IS (SEE COMBUSTION EQUATION )
H20I=(15*18/198.0} *MFI (K)
KG
MB20=MH20+H20I
KG
CVGAS=(M02*CV02+MN2*CVN2+MH20*CVH2O+MC02*CVC021/
/(M02+MN2+MH20+MC021
J/KG.K
CVGAS=«704.9+O.0708*TCY(KI)*Ml+(1833.4+O.4931+TCY(K»*MAXMF)/
/ (Ml+MAXMF)
J/KG.K
C CVGAS=838.
J/KG.K
MWGAS=(M02*32.+MN2*28.+MC02*44.+MH20*18.1/(M02+MN2+MC02+MH20) G/MOL
C COMBUSTION TEMP. CALCULATIONS
TI2=TCY(KI
K
50
PI2=PCY(K)
PCY(K)=PI2*(TRAMAS+IlEATRL(NCOHB»*TCY(K)/(TRAMAS*TI2) BAR
910 CONTINUE
Explanation:
release function
release curve 2
(KGFUEL. DEG )
TRAMAS trapped mass of gases along the cycle (KG)
MFCUR similar to MFI, used for plotting.
51
For obtaining the specific heat, disas.sociation of gases, molecular
weight, temp. and pressure see the Combustion Under Constant Volume para.2.2.4.s.3.
IF (NP3MAX.HQ.l) GO TO 905
NP3MAX=NP3MAX+l
P3(I)=PI
905 CONTINUE
DDEGT=(DDEGRE*60.0)/(NE*360.0) SEC
TI3=TCY(K)
PI3=PCY(K)
RVALV=O.018
DLTAP=(PI3-Pl)*RVALV*EXP(-DDEGT) BAR
PCY(K)=PI3-DLTAP BAR
TCY(K)=T13* «PCY(K)/PI3)** «N2-l)/N2» K
TRAMAS=PCY(K)*VCY(K)*MWGAS*lE-3/(R*TCY(K» . KG
903 CONTINUE
IF(K.NE.(3*NSTP» GO TO 324
52
NWRlTE=CWRlTE
324 CONTINUE
NWRlTE=O
328 CONTINUE
NWRlTE=NWRlTE+1
VI=VCY(K) LT
TI=TCY(K) K
PI=PCY(K) BAR
900 CONTINUE
Explanation:
function. (BAR)
adiabatically (K)
53
TRAMAS The residual mass of gases each step (KG)
903 CONTINUE the end of exhaust valve model.
expansion stroke
PRVVCY=VI
WRITE {1,3l01
DDEGRE={l80,0/FLOAT{NSTP 11
DEG
DEGRE=DEGRE+DDEGRE
DEGRE
OLDMAS=TRAMAS
PHI=ABS{DEGRE )*3.14/180.01
RADIAN
VCY{KI=V2+{{V1MAX-V21/2.01*
DVCY=VCY(KI-PRVVCY
OUTMAS=TRAMAS*DVCY/PRVVCY
KG
PRVVCY=VCY(K)
TRAMAS=TRAMAS-ABS (OUTMASI
KG
54
C MASS ESCAPED DUE TO HIGH PRESSURE
DLTAP=(PI-P1)*RVALV*EXP(-DDEGT) BAR
PCY(K)=PI-DLTAP BAR
TCY(K)=TI*«PCY(KI/PI)**«N2-1)/N2» K
TRAMAS=PCY(K)*VCY(K)*MWGAS*lE-3/(R*TCY(K» KG
VI=VCY(K)
PI=PCY{K)
TI=TCY(K)
RESDUL--TRAMAS/OLDMAS 0/0
M02=M02*RESDUL KG
MN2=MN2*RESDUL KG
MC02=MC02*RESDUL KG
MH20=MH20*RESDUL KG
CV02=(lOOO/32.01*(20.959+O.00295*(TCY(K)-273.» J/KG.K
CVN2=(lOOO/28.0)*(20.704+0.00211*(TCY(K)-273.» J/KG.K
CVGAS=(M02*CV02+MN2*CVN2+MH20*CVH2O+MC02*CVC02)/
/(M02+MN2+MH20+MC02) J/KG.K
CVGAS=«704.9+O.0708*TCY(K»*M1+(1833.4+o.4931*TCY(K»*MAXMF)/
/ (Ml+MAXMFI J/KG.K
CVGAS=838 J/KG.K
IF (NWRlTE.LT. CWRlTE) GO TO 335
NWRlTE=O
335 CONTINUE
NWRlTE=NWRITE+1
1100 CONTINUE
55
NP3MAX=O
(2.2.4.6.5)
56
IF (NWRITE ••• ) To control the writing of the main variable in
See (2.2.4.5.5)
WORK=WORKEX+WORKC BAR. LT
IMEP=WORK/(VlMAX-V2) BAR
WRITE (l,36B) IMEP
C CALL Cl05lN
C NCY=NSTPC+NSTPEX
C CALL GRAPOL(VCY,PCY,NCY)
C CALL GRAPOL(VCY,MFCUR,NCY)
C CALL OEVEND
WORK=WORKEX+WORKC
57
The indicated mean effective pressure CIMEP) in this cycle is the
nett work over the swept volume, i.e.
WORK
IMEP = VIMAX-V2
writing the indicated mean effective pressure in the current cycle. Calling
the line printer to plot the pressure of the cycle against the volume and
CALL GRAPOL C'-) to plot the volume of the cycle against the pressure,
CALL GRAPOL (--) to plot the mass of fuel injected Cor burned) through
VCY(K}=O.O
PCY(K)=O.O
1200 CONTINUE
next cycle, otherwise spurious values will be written by means of WRITE (1,320) ,
58
2.2.4.10. Exhaust"Parameters (Before Turbine)
P5EXH ~ (PI+P5(I»/2.0
T5EXH~ (TI+T5(I»/2.0
WRITE(2,363) T5EXH
BMEP(I)~IMEP-FMEP
BAR
WRITE(2,375) BMEP(I)
Writing out the value of friction mean effective pressure in this cycle
Writing out the brake mean effective pressure of the current cycle.
RPM(I)~NE
R. P .M.
BHP(I)~BMEP(I)*V1MAX*NE/225.0*STROKS) H.P.
IHP(I)~IMEP*V1MAX*NE/(225.0*STROKS) H. P.
'BHP(I) is the brake horse power of the engine at the current cycle
59
It depends upon the swept volume engine, the speed and the number
effective pressure.
of the engine.
C TURBINE POWER
GAMAEX=l. 30
C TURBINE BLADESPEED=6.28*NTUR*RADUS/60.0
BLSPED=6.28*(30.*NE)*0.05/60.0
EXPVEL=SQRT(2.*CVGAS*T5E>H*(1.-(PA/P5EXH )**«GAMAEX_l)/GAMAEX»)
ETATUR=1.2*(BLSPED/EXPVEL) **(0.81)
TURPWR=ETATUR*GJ\SMAS*CVGAS *TSEXH *
WRITE(1,435) TMPOUT
Explanation:
60
BLSPED The turbine average blade speed
M/SEC
EXPVEL The isentropic expansion velocity of the gases
ECC=300
T2W(I)=T1W+ALLHEAT/(RCC+ECC)
K
WRITE(2,4501 T2W(I)
IF(T2W(I).GT.350) GO TO 1000
T2W(II=353
1000 CONTINUE
,
61 • i
ALLHEAT is the total amount of heat released by the engine; it depends
upon mass of fuel consumed per second, the calorif1c value of the fuel, the
turbine output temp., atmospheric air temp.
TlW=TA+lO is the radiator output temp. (or the engine cold water
temp. ) It is considered to be 10 degrees higher than the atmospheric temp.
RCC, ECC is the cooling coefficient of the radiator and engine block.
RCC depends upon engine speed (or the water pump speed).
i.e. when this temp. increases above 350 Ko then this thermostat
.
opens further to increase the rate o'f flow of hot water (but this instrument
is not accurate in practice}.
engine brake horse power, engine torque and specific fuel consumption.
This is the end of the calculation of the cycle variables. (NPT) cycles
are calculated at constant engine speed and variable injection advance; the
output variables (such as the power) can be plotted against the injection
62
To obtain the optimum advance for another value of engine speed, simply
change the speed, calculate a few cycles and deduce the optimum again.
IF (J.NE.l) GO TO 61
DO 15 1=1, NPT
BHP1(I)=BHP(I)
TORQ1(I)=TORQ(I)
15 CONTINUE
61 CONTINUE
IF (J.NE.2) GO TO 66
DO 4i\I=l,NPT
BMEP2(I)=BMEP(I)
TORQ2(I)=TORQ(I)
FUEL2(I)=FUEL(I)
41 CONTINUE
66 CONTINUE
IF (J.NE.3) GO TO 69
63
DO 27 I=l,NPT
P33(I)=P3(I)
BMEP3(I)=BMEP(I)
TORQ3(I)=TORQ(I)
FUEL3(I)=FUEL(I)
T53(I)=T5(I)
P53(I)=P5(I)
27 CONTINUE
69 CONTINUE
NE=NE+5oo
7 CONTINUE
necessary to store all the output variables related to each engine speed,
IF(J.NE.l) GO TO 61
means that if it is the first value of engine speed (500 r.p.m.) then
NPT -- is the number of times which the advance has been tested at
constant engine speed. So this DO statement stores those results in more
arrays similar to the original arrays.
AIRMSl(I)=AIRMAS(I)
P31(I)=P3(I)
BMEPl(I)=BMEP(I)
BMEP(I) is an array to store the values of brake mean effective
TORQl(I)=TORQ(I)
FUELl(I)=FUEL(I)
PSI (I)=PS(I)
15 - continue
61 - continue
65
00 41 I=l, NPT
NPT -- is the number of times which the advance has been tested at
BMEP2(I)=BMEP(I)
BBP2(I)=BBP(I)
TORQ2(I)=TORQ(I)
FUEL2(I)=FUEL(I)
T52(I)=T5CI)
are the variables related to the second value of engine speed (1000 r.p.m.)
Those variables are: max. combustion pressure, brake mean effective pressure,
brake horse power, output torque, specific fuel consumption, exhaust manifold
41 - continue
66- continue
DO 27 I=l, NPT
P33(I)=P3(I)
BMEP3(I)=BMEP(I)
BBP3(I)=BBP(I)
66
TORQ3(I)=TORQ(I)
P53(I)=P5(I)
are the variables related to the third value of engine speed (1500
r.p.m.) Those variables are, max. combustion pressure, brake mean "
27 - continue
69 - continue
(2000-500=1500 r.p.m.) and 4 speeds were tested therefore the step increase
1500
is --3- = 500 r.p.m.
7 - continue
CALL C1051N
CALL PENSEL(l,O.4,3)
CALL PENSEL(2,0.4,3)
CALL GRAPOL(ANGLIJ,BMEP1,NPT)
67
CALL MOVT02 (25. ,110.)
CALL CHAHOL(9HGRAPH-l*.)
CALL PENSEL(3,0.4,3)
CALL GRAPOL(ANGLIJ,BMEP2,NPT)
CALL CHAHOL(9HGRAPH-2*.)
CALL PENSEL(4,0.4,3)
CALL GRAPOL(ANGLIJ,BMEP3,NPT)
CALL CHAHOL(9HGRAPH-3*.)
CALL PENSEL(5,0.4,3)
CALL CHAHOL(9HGRAPH-4*.)
CALL PICCLE
CALL PENSEL(1,0.4,3)
CALL AXIPLO(0,150.,100.,1,1,10,10,MAXINJ,MININJ,0.,260.,
CALL PENSEL(2,0.4,3)
CALL CHAHOL(9HGRAPH-l*.)
CALL CHAHOL(9HGRAPH-2*.}
CALL PENSEL(4,0.4,3)
CALL CHAHOL(9HGRAPH-3*)
CALL PENSEL(S,0.4,3)
CALL CHAHOL(9HGRAPH-4*.)
CALL PICCLE
CALL PENSEL(1,0.4,3)
68
CALL AXIPLO(0,150.,100.,1,1,10,10,MAXINJ,MININJ,,1,,170,'INJECTION
CALL PENSEL(2,0.4,3)
CALL GRAPOL(ANGLIJ.FUEL1,NPT)
CALL CHABOL99HGRAPH-l*.)
CALL PENSEL(3,0.4,3l
CALL GRAPOL(ANGLIJ,FUEL2,NPTl
CALL CHABOL(9HGRAPH-2*.)
CALL PENSEL(4,0.4,3)
CALL CHABOL(9HGRAPH-3*.1.
CALL PENSEL(5,0.4,31.
CALL CHAHOL(9HGRAPH-4*.)
CALL PICCLE
CALL PENSEL(1,0.4,3)
CALL AXIPLO(0,150.,100.,1,1,10,10,MAXINJ,MININJ,50.,99.,
CALL PENSEL(2,0.4,3)
CALL CHABOL(9HGRAPH-l*.)
CALL PENSEL(3,0.4,3)
CALL GRAPOL(ANGLIJ,TORQ2,NPT}
69
CALL CHAHOL(9HGRAPH-2*.)
CALL PENSEL(4,O.4,3)
CALL GRAPOL(ANGLIJ,TORQ3,NPT)
CALL CHAHOL(9HGRAPH-3*.)
CALL PENSEL(5,O.4,3)
CALL CHAHOL(9HGRAPH-4*.)
CALL PICCLE
CALL PENSEL(1,O.4,3)
CALL PENSEL(2,O.4,3)
CALL CHAHOL(9HGRAPH-l*.)
CALL PENSEL(3,O.4,3)
CALL CHAHOL(9HGRAPH-2*.)
CALL PENSEL(4,O.4,3)
CALL CHAHOL(9HGRAPH-3*.)
CALL PENSEL(5,0.4,3)
CALL CHAHOL(9HGRAPH-4*.)
CALL DEVEND
STOP END
70
Plotting of the variables produced from each speed, but each family
of curves is plotted on one common axis; these families are: brake mean
CALL MOVT02 ( ---I is to move the monitor* to a proper position for writing
the number of the curve
CALL CHAHOL ( ---I is to write the number of the curve
CALL PENSEL ( ---) calling another colour to plot second curve
CALL GRAPOL ( ---I plotting relation between angle of injection and brake
CALL PENSEL ( ---) to call another colour to plot the fourth curve
CALL GRAPOL ( ---I to plot the relation between angle of injection and
speed
CALL PICCLE to clear all the picture because another family will be
plotted.
are four curves, each one related to one value of engine speed; they are
plotted on one common axis but they are different in colour, and each colour
pressure.
Now the family of fuel consumption curves is plotted, there are four
curves, each one related to one value of engine speed. They are plotted on
one common axis but they are different in colour, and each colour is numbered.
For more detail, see plotting of brake mean effective pressure. Also the
family of output torque curves is plotted; they are four curves plotted on
one common axis, each one is related to one value of engine speed and has a
they are four curves plotted on one common axis, each one is related to
one value of engine speed and has a distinctive colour and number.
STOP
END
72
2.3 Discussion
No. of strokes = 4
No. of cylinders = 6
- Combustion pressure
- Combustion temp.
- work done.
Turbine power
73
Specific fuel consumption
3 to 5 degrees.
The mass of output gases depends upon the amount of fuel injected
per cycle.
the T.D.C.
2. The second column represents the volume of the cycle at each step,
the max. volume at the B.O.C. is 12.04 (Lt), while the min. volume at the
3. The 3rd column is the trapped mass of gases through the cycle.
4. a)The 4th column represents the pressure at each step of the cycle;
c)The max. combustion pressure never occurs before T.O.C. tas some
books describe) but always at the T.O.C. or after depending upon the
injection advance.
d)In the optimum case the max.· combustion pressure occurs a few
74
:
degrees after T.O.C.
5. a)The 5th column is the temp. which behaves like the pressure.
b)The initial temp. depends upon the residual gase"stemp. and the"
d)The max, temp. like max. pressure never occurs before T.O.C. but
8. a)The seventh column is the work done so far at the current step
75
work in the expansion stroke starts to count. Its max. value is
at the B.O.C. and the nett work in each cycle is the algebraic sum
curve.
1. On Graph 2.1 the first set of curves is the BMEP against the injection
advance angle.
It is clear that the BMEP varies with the injection advance, and
(optimum advance).
Also the max. engine speed can be increaseR, i.e. the power to
76
2. Graph 2.2 shows the output power against injection advance. It
proper choice of injection angle, and that the value of the optimum
4. Graph 2.4 gives the variation of output torque with injection advance.
present.
5. Graph 2.5 gives the max. combustion pressure with the injection advance.
Also the value. of max. pressure decreases with increasing engine speed.
affected.
The computer model never reaches 100% accuracy and I believe that
the fuel and power saving could be greater than the model shows.
77
TABLE (2.1) 'IHE DIGITAL mMIUTER ClJTPUT FOR CYCLE CALQJIATION
78
347. 1.03286 .0176545 34.71 715.9 28.94 754.79 -60.69 0.106551El-04
348. 1.00498 .0176798 38.75 776.5 28.97 758.89 -61.72 0.252913El-04
349. 0.97924 .0177134 43.85 854.6 29.01 766.00 -62.78 o.336343El-04
350. 0.95566 .0177561 50.12 950.9 29.06 775.22 -63.89 0.426649El-04
351. 0.93425 .0178091 57.70 1067.1 29.12 786.81 -65.04 o.529475El-04
352. 0.91503 .0178691 66.24 1195.7 29.19 800.75 -66.23 o.600511El-04
353. 0.89800 .0179274 74.66 1318.3 29.25 815.77 -67.43 0.583167El-04
354. 0.88316 .0179793 82.34 1425.9 29.31 830.02 -68.60 0.518865El-04
355. 0.87053 .0180253 89.28 1520.0 29.36 842.77 -69.68 0.459544&-04
356. 0.86012 .0180658 95.44 1601.8 29.41 854.08 -70.64 o.40 5566El-04
357. 0.85192 .0181015 100.82 1672.7 29.45 864.07 -71.45· 0.356918El-04
358. 0.84595 .0181328 105.42 1733.8 29.48 872.83 -72.06 0.313382El-04
359. 0.84220 .0181603 109.24 1786.0 29.51 880.47 -72.46 0.274629El-04
360. 0.84068 .0181843 112.30 1830.3 29.54 887.09 -72.63 0.240278El-04
CALCllLATIOO OF EXPANSION SrRa<E
80
,rr (Si',,::"
[l
-
. \ 1
,
) .
) .
'.,
. :f
\1
'.;
:
!
GRftJ'l+- .2 •1
.1.. .. __ ....... _._. ______ .._...
._-,--------''''''::-,
-15 -10 -3~ -38 -25 -20 -1:1 -12 -5
I NJfcr 1O~, ANRfC ces )
81
spriT' SOO r.p.m.
GRAPH - 2.3
.- .. _-.-._----_._
. .•--------_._-_.
. . _ - - - - - - '--'
-15 -10 -2~ -20 _ 1 c:
. ..J -18 -5 G
C·
f,
. ,--- ---------- 4
.--" ---,
1000 r.p.m •
----
-- .----------l---------.~
----·--_SP[[ r.:
500 r.p.m.
GRAPH - 2.2
--l'; "l,'
." . ..J
_ 1 (
~
--
. ," .
~',
._. ~)
J:-, ( ,
,
1 I ~ •
!
, t·,
r" . r'
I ,:
-'-----
-- -.
'----- ---
----
:;)UL (KGr. ~~ )
\
\
GRAPH - 2.4 ,,
/
\
\
srC[ L' - 1 500 r.p.m.
.......
1000 r.p.m.
"
------ ... -~,
1500 r.p.m. ..... ,' .....
2000 r.p.m.
_. ") r
.... _--_ ... _._._.- --,
.' , _. :·~8 .- i
1 C·
.....)
1 r",
') G
"
83 . T'! T!' ( '1 T:
/"
• 'r
, ,
!"..J l L '.
Chapter 3 Verification of the Diesel Engine Model
3.1 Introduction
compression ignition engine, the aim being to find the optimum rMEP and
of Chapter (2) and the experimental results obtained from the laboratory,
as shown later.
Another input variable tested was the input air temp., the experiment
being carried out in different weather conditions, and with good agreement
engine.
the value of the load is adjustable to keep the engine at steady state
speed for variable angle of injection; the value of the air/fuel ratio,
load torque, B.H.P., max. combustion pressure, and carbon monoxide (CO)
84
INJEcrION ADVANCE TEST
Trial No Advance
- angle Air/Fuel Torque
B.H.P B.M.E.P. SP. Fuel Thermal
0
1 (deg. ) ratio Cons tunption Max. / Co
KgF.m BAR 0/07 0
gm/B.HP.HR PresSurE
1__ - - - Bar
Table (3.1a)
INJEcrION ADVANCE TEST
Table (3.1b)
INJEcrION ADVANCE TEST
Trial No Advance
-
angle Air/Fuel Torque B.H.P B.M.E.P.
SP. Fuel Thermal
0/
0
1 (deg.)
ratio
KgF.m BAR
Consumption Max.
Pressure 07
/ Co
0
gm/B.HP.HR __!l_a];"_
--- ---- ----- --- ----,. - . - --t------ 1----
30 19.8 2.52 5.28 6.25 240 64 25.15 0.08
-- r-------- - .. ~.
Table (3.1e)
INJEcrION ADVANCE TEST
Trial No Advance
- angle
Air/Fuel Torque
B.H.P B.M.E.P. SP. Fuel Thermal
0/
o /o~
1 (deg. ) ratio Consumption Max. Co
KgF.m BAR 0
Pressure
grn/B.HP.HR
_._. --- Bar
Table (3.ldf.
INJEcrION ADVANCE TEST
,
40 " 2.39 8.0 5.92 253.4 75 23.82 0.11
Table (3. le )
INJECI'ION ADVANCE TEST
Trial No Advance
Air/Fuel Torque B.H.P SP. Fuel Thermal
angle
ratio
B.M.E.P. 0
Consumption Max. / Co
o'" 1 (deg. ) KgF.m BAR
Pressure 0/07 0
gm/B.HP.HR -1L<u:.. __
-. -- ..---. - r---
30 16.2 2.08 7.97 5.15 291 50 20.75 0.11
- ._.
35 " 2.24 8.59 5.55 270 62 22.36 0.11
.
Table (3.1fl
.-;
.--.~.
'. P (bar)
!.
max. combustion
':(
,
. .:j. 'cc ,•• \
I. .i :1 .
I:'
+~. -1,.1
...
~.:: 'I>. "
• I.~
---:i
:1 ,
..... ---:.
~.,
f -i""
;:---
i:~
'i-,'+
~--
'.
.. ,
,
. I ..
"!:":
.. '-',
I"
:FfY!:'lit~''-+-:+''+:;+'L*::++:.\l;;4B4:-+:-':'·.;.j'''~·'·J\mti,''Ij (deg.)
Injection advance
6e:~rl±~T:±~~T=~=r~~~~~b±~±±~~~F±~~~"~":
-1--- • - • -_.- - .
·B.M.E.P. (Bar)
._0 ___ .... _. _._. ____ .....
..... _--.,.- --
Fig. (3-2b)
1250 r.p.m .
.---r -.-----.-.---~-j::::::::.:-~--~.~-~-~--=--:·--:--=::~:==~::...11ooo r. p. m.
1500 r.p.m.,
5-+--______---.;..___-.;...___...:-___+-__---!_ _ _ANGLIJ ~
(deg)
INJECTION
-45 -40 -35 -3b ADVANCE
;
, i
._-----------1.------- t· .. I
- . -, ._. -- "T-'
---r ---f-~--- ----'.-_.._. ----.-
I
:,. ·"·~l
Fig. (3 - 2 Cl
/
4
7::11
5J j
GRftPH • 2.1
G0.l. ___ ....____ .._ _ _ _ _:--_ _:--_---:_ _-...,._ _-...,._ _--,-,_::::.....",
-15 -10 -3~ -38 -25 -20 -15 -13 -5
INJrCTIOrl ANGLECOEG)
..
'1AX 1 rr:[·c('·
. • i'\. .• ) ;,)
'<13'
18·
16
1 .,!
12 1
13.
Cr'[·[·[J
..:;to.. . '))
1500 r.p.m.
93
3.2.2 Comments and Conclusion
From graph No. (3.2.a) it can be seen that, when increasing the injection
advance the max. pressure increased, and when decreasing the injection advance,
the max. pressure has no chance to build up, so it drops sharply, and that
.,
is typically showed by the digital computer model (end of Chapter 2), or Fig(3.2C).
angle, see graph No. (3.2.b) and the digital computer graphs at the end
of the pressure curve graph (3.2.a) and the position of the optimum advance
air which helps to increase the engine output power, provided this decrease
Tables (3.3.1) and (3.3.2) give the performance of the laboratory engine
The digital computer was tested and the input air temp. was reduced
o 0
from 298 K to 288 K so the output IMEP calculated was 12.414 and 12.458
Bar respectively. i.e. decreasing the input air temp. bYIOoC causes
increasing the IMEP by less than 0.5%, thus agreeing with the practical
observations above.
The effect of input air temp. upon engine performance is very small,
provided that the change is around room temp. The effect of the input
95
INJECTION ADVANCE TEST
Trial No Advance -
-
angle Air/Fuel Torque B.H.P B.M.E.P. SP. Fuel Thermal
1 0/
ratio Consumption Max.
%~
Co
(deg. ) KgF.m BAR
gm/B.HP.HR PressurE °
--- ---_..•. .-- --.---.-- .- -'La.L +-
List of Symbols
A. element of area
~
DNE" !!
dt
NE angular acceleration
98
P4 or P pressure at the end of constant pressure combustion
4 (bar)
P4V4 or P4V4 - P4 * V4 pressure and temp. at the end of constant
99
work done in compression stroke without
combustion (bar.Lt.)
W2 or W extra work done due to combustion (bar.Lt)
2
W3 work done at constant pressure combustion (bar. Lt)
W4 work done in expansion stroke and after constant
100
4.1 Introduction
matrix.
(See below)
b - Inlet air temp: the effect of variable input air temp. can
constant.
101
e - Instants of opening and closing of valves: they are assumed
are carried out to linearize the real heat release function and
study its effect upon combustion pressure rise and engine performance.
is used.
102
4.3 The Error Generated by Using a Constant Specific Heat
The non linear dependance of specific heat upon the temp. and upon
simplified further.
required, see eqn. (4.5) be low, and this value was initially chosen by
the output performance, Table (4.2) compared with that in Table (4.1).
represents the average value of specific heat over the combustion period.
This value can be obtained from the digital computer model which gives
almost the same theoretical engine performance after a few trial runs
of the model.
103
have its greatest effect on the cylinder pressure and temperature but
104
TABLE (4.1) ':mE DIGITAL OOMRJTER aJTRJT FOR CYCLE C'ALaJIATION.
CYCLE WI'lli VARIABLE SPCIFIC HEAT.
105
320. 2.52340 .0176439 10.56 532.2 28.93 742.58 -33.44 O.OOOOOOE 00
330. 1.81805 .0176439 16.14 586.2 28.93 746.40 -42.61 O.OOOOOOE 00
340. 1.28687 .0176439 25.25 649.1 28.93 750.86 -53.29 O.OOOOOOE 00
347. 1.03286 .0176545 34.71 715.9 28.94 754.79 -60.69 0.106551E-04
348. 1.00498 .0176798 38.75 776.5 28.97 758.89 -61.72 0.252913E-04
349. 0.97924 .0177134 43.85 854.6 29.01 766.00 -62.78 0.336343E-04
350. 0.95566 .0177561 50.12 950.9 29.06 775.22 -63.89 0.426649E-04
351. 0.93425 .0178091 57.70 1067.1 29.12 786.81 -65.04 0.529475E-04
352. 0.91503 .0178691 66.24 1195.7 29.19 800.75 -66.23 0.600511E-04
353. 0.89800 .0179274 74.66 1318.3 29.25 815.77 -67.43 0.583167E-04
354. 0.88316 .0179793 82.34 1425.9 29.31 830.02 -68.60 0.518865E-04
355. 0.87053 .0180253 89.28 1520.0 29.36 842.77 -69.68 o.459544E-04
356. 0.86012 .0180658 95.44 1601.8 29.41 854.08 -70.64 0.405566E-04
357. 0.85192 .0181015 100.82 1672.7 29.45 864.07 -71.45 o.356918E-04
358. 0.84595 .0181328 105.42 1733.8 29.48 872.83 -72.06 0.313382E-D4
359. 0.84220 .0181603) 109.24 1786 .0 29.51 880 .47 -72.46 0.274629E-04
360. 0.84068 .0181843 112.30 1830.3 29.54 887 .09 -72.63 0.240278E-04
CALOJLATIOO- OF EXPANSION SIROKE
106
- - - - - ----"
107
TABLE (4.2) _ ENGlliE CYQ.E WI'IH OONSTANT SPCIFIC HEAT
108
348. 1.00498 .0176767 38.88 779.1 28.97 830.00 -62.53 0.252913&-04
349. 0.97924 .0177104 43.70 851.7 29.01 830.00 -63.59 o.336343E-04
350. 0.95566 .0177530 49.64 942.0 29.06 830.00 -64.69 0.426649E-04
351. 0.93425 .0178060 56.90 1052.5 29.12 830.00 -65.83 0.529475&-04
352. 0.91503 .0178660 65.18 1176.8 29.19 830.00 -67.01 0.600511&-04
353. 0.89800 .0179243 73.45 1297 .3 29.25 830.00 -68.19 0.583167E-04
354. 0.88316 .0179762 81.11 1404.8 29.31 830.00 -69.34 0.518865&-04
355. 0.87053 .0180222 88.10 1500.1 29.36 830.00 -70.40 0.459544E-04
356. 0.86012 .01806Z7 94.37 1584.1 29.41 830.00 -71.35 0.405566E-04
357 • 0.85192 .0180984 99.90 1657.7 29.45 830.00 -72.15 0.356918&-04
358. 0.84595 .0181297 104.67 1721.7 29.48 830.00 -72.76 0.313382&-04
359. 0.84220 .0181572 108.67 1777.0 29.51 830.00 -73.16 0.274629E-04
360. 0.84068 .0181812 111.91 1824.2 29.54 830.00 -73.33 0.240Z78E-04
CALaJIATIOO OF EXPANSION STROKE
109
427. 4.98027 .0183407 '14.44 1382.7 29.71 830.00 158.80 O.OOOOOOE 00
437. 6.00319 .0183407 11.45 1320.9 29.71 830.00 171.94 O.OOOOOOE 00
447. 7.01784 .0183407 9.42 1271.3 29.71 830.00 182.47 O.OOOOOOE 00
457. 7.98916 .0183407 8.02 1231.5 29.71 830.00 190.91 O.OOOOOOE 00
467. 8.88859 .0183407 7.02 1199.8 29.71 830.00 197.65 O.OOOOOOE 00
477. 9.69484 .0183407 6.30 1174.5 29.71 830.00 203.01 O.OOOOOOE 00
487. 10.39351 .0183407 5.78 . 1154.6 29.71 830.00 207.23 O.OOOOOOE 00
497. 10.97612 .0173686 4.88 1116.8 29.71 830.00 210.38 O.OOOOOOE 00
507. 11.43856 .0158125 4.12 1080.5 29.71 830.00 212.47 O.OOOOOOE 00
517. 11.77959 .0145213 3.58 1050.6 29.71 830.00 213.79 O.OOOOOOE 00
527. 11.99944 .0134480 3.18 1026.3 29.71 830.00 214.54 O.OOOOOOE 00
537. 12.09878 .0125510 2.88 1007.1 29.71 830.00 214.84 O.OOOOOOE 00
540. 12.10517 .0123107 2.81 1002.2 29.71 830.00 214.86 O.OOOOOOE 00
CALalIATION OF EXHAUST S'lROKE
110
4. 4 The Linearization of the Heat Release Function.
In the digital computer program Chapter (2) the heat release curve
used during combustion of the fuel is a realistic one such as that shown
in Fig. (4.3).
of two straight lines as in Fig. (4.5), and the third of one straight line
Fig. (4.6). The three shapes develop reasonable accuracy if their centroids
are coincident with that of the real heat release curve, (where the centroid
In this case the simplest shape is chosen, Fig. (4.61 to replace the
The centre of gravity of the linearized heat release curve (i.e. the
111
EL BURNED/DEG
mg/deg.
inear
330 335 340 345 350 355 360 365 370 375 380 385 390
CRANK ANGLE
(Deg. )
Fig. (4.3)
112
Heat release/deg Heat release/deg
deg
deg
Fig (4.4) Fig (4.5)
Heat release/deg.
Heat release/deg.
deg
deg 2L :-
I
centroid
Fig (4.6) Fig (4.7)
113
but for the real heat release curve, the centioid can be calculated using
the formula:
i=CMBLTH
E Ai • Xi
Centroid = =i_=~l~~~=-_____
i=CMBLTH
E
i=l
From the results in table (4.1) and table (4.3) it can be seen that
the two cycles are different in injection advance angle only. This is
necessary to make the centroids of the heat release curve coincident. Note
that the largest error appears in the max. combustion pressure which is about
3%.
in the shape of a triangle which" is very easy to handle. The only effect
position of the centroid of each heat release model with respect to the
T.D.C.
The effect of linearization upon the max. combustion temp. is very small
the heat release curve upon the dynamic performance of the engine.
114
TABLE (4.3) f.I ENGINE CYCLE WITH TRIANGULAR HEAT RREASE CURVE
115
343. 1.16558 .0174126 38.33 904.6 29.05 773.83 -56.38 0.255671E-04
344. 1.12925 .0174376 42.34 966.5 29.08 780.60 -57.84 0.250288E-04
345. 1.09502 .0174621 46.47 1027 .2 29.11 787 .30 -59.36 0.244906E-04
346. 1.06288 .0174860 50.70 1086 .4 29.14 793.91 -60.92 0.239523E-04
347. 1.03286 .0175094 55.02 1144.2 29.16 800.42 -62.51 0.234141E-04
348. 1.00498 .0175323 59.41 1200.5 29.19 806.82 -64.11 0.228758E-04
349. 0.97924 .0175546 63.83 1255.2 29.22 813.11 -65.69 0.223376E-04
350. 0.95566 .0175764 68.26 1308.3 29.24 819.28 -67.25 0.217993E-04
351. 0.93425 .0175977 72.65 1359.7 29.27 825.31 -68.76 0.212611E-04
352. 0.91503 .0176184 76.97 1409.2 29.29 831.19 -70.20 0.207228E-04
353. 0.89800 .0176386 81.17 1456.8 29.31 836.91 -71.54 0.201845E-04
354. 0.88316 .0176582 85.21 1502.4 29.34 842.46 -72.78 0.196463E-04
355. 0.87053 .0176774 89.04 1545.8 29.36 847.82 -73.88 0.191080E-04
356. 0.86012 .0176959 92.63 1587.1 29.38 852.99 -74.82 0.185698E-04
357. 0.85192 .0177139 95.91 1626.1 29.40 857.96 -75.60 0.180315E-04
358. 0.84595 .0177314 98.86 1662.7 29.42 862.71 -76.18 0.174933E-04
359. 0.84220 .0177484 101.44 1696.9 29.44 867.24 -76.55 o.169550E-04
360. 0.84068 .0177648 103.62 1728.7 29.46 871.53 -76.71 o.164168E-04
CALalLATION OF EXPANSION STRa<E
116
440. 6.31016 .0180070 10.80 1334.0 29.73 863.57 173.76 O.OOOOOOE 00
450. 7.31522 .0180070 8.98 1286.6 29.73 859.45 183.65 O.OOOOOOE 00
460. 8.26764 .0180070 7.71 1248.6 29.73 856.15 191.57 O.OOOOOOE 00
470. 9.14103 .0180070 6.81 1218.2 29.73 853.52 197.90 O.OOOOOOE 00
480. 9.91621 .0180070 6.15 1194.2 29.73 851.43 202.91 O.OOOOOOE 00
490. 10.58074 .0183146 5.60 1172.4 29.73 849.79 206.84 O.OOOOOOE 00
SOO. 11.12757 .0165841 4.65 1130.1 29.73 846.10 209.65 O.OOOOOOE 00
510. 11.55362 .0151487 3.96 1095.2 29.73 843.04 211.50 O.OOOOOOE 00
520. 11.85824 .0139583 3.47 1066.6 29.73 840.53 212.64 O.OOOOOOE 00
530. 12.04186 .0129682 3.10 1043.6 29.73 838.51 213.24 O.OOOOOOE 00
540. 12.10517 .0121392 2.84 1025.6 29.73 836.93 213.43 O.OOOOOOE 00
CALClJLATION OF EXHAUST S'lROKE
117
4.5 The process of Linearization of the model
introduced to help calculate the values of the expressions. using the same
T2 P 2
L
p ..
l.n) f
P I.
exh f
volume
118
4.5.1 Calculation of the Analytic Expression Describing the P-V Diagram
Fig. (4.8)
P
inj = [
V l{l
• P
l (4.1)
inj
T
inj = t~ V
inj
f l • Tl (4.2)
P
2 = [:~rl • P
l (4.3)
T2 =
rr
V~
l
• Tl
~r
(4.5)
so T3
P3 P2 (4.6)
T2
119
SIMULATION ~ W. (4.8) 'ID IESCRIBE 'mE DIESEL EN>INE CYCLE OF FIG. (4.8)
C FIXED mTA
Vl=12.105
V2=O.84068 LT
EPSLUM= (Vl/V2)
Nl.=1.28
N2=1.22
S'JROI<S=4
INERTE=0.5368 z 5 H:;M.S2/R
INER'l'IP2.209
MFtWC:=O .0007
~. ,". S !,...d. llGM.S2/R
RG/CYClo
W0lnl=-50
C INRJT VPRI1\BLES
Pl=I.3 BI\R
Tl=300 K
Y=100 EUEr. 0/0
ProLLJ=29 ra;
C INrrIAL a:NmTIOOS
NE=2000. R.P.M
C OF ADJUSTJBLE PARAl£~
CAI.<lJIJ\TJ:CE
IX)99 1=1,5 .
EXLOAD=O.005*(Vl-V2)*NE
ProLLJ=29
WRlTE(I,10) NE , EXLOIID
10 FORM!\T(lHO,'ENGINE SEEED=' ,F5.0,'lU'M' ,5X, 'EXLOAD=' ,E7.1,'lGM')
MF=Y*MElWV'100.0
IX) 101 J=1 ,11
S'lRTCM= (miLLJ-O.006*NE) om
EHI=ProLLJ*3.1~180.0 RAD.
VINJ=V2+«Vl-v2)/2.0)*(4.3333-aDS(EHI)-3.3333*SQRr(
(1-o.09*(SIN(EHI)**2»)
MroJ=+(SlRTOJ,/24.) *MFtWC: RG/CYCL
MFCP=MF-MEW
WO=(3.1~180.0)*(I80.+WOPEN) .
~=V2+«V1-V2)/2.0)*(4.33333-<DS(EVO)-3.3333*SQRr(
. (1.-0.09*(SIN(EVO)**2»)
C HIRJ!Il.E'l'Eru OF mE CYCLE
P2=(EPSLUM**Nl)*Pl
T2=(EPSLUM** (Nl. -1) ) *Tl
T3=T2+(2.10E6)*MroJ
P3=P2+O.1857E6*MFCV
C P3=P2+( 0.07272) * ('1'3-'1'2) .
. P4=P3
T4=T3+(1.72E6)*MFCP
V4=V2+( 820.0) *MFCP
P5=«VZ/VEXH)**N2) *P3+(6500.) *MFCP
P4V4=P2*V2+O.1527E6*V2*MFCV+128.58E3*MFCP
120·
P4V4=P.2*V2+O.1527E6*V2*MFCV+128.58E3*MFCP
Wl=(l/ (Nl-l) )* (P2*V2-Pl*Vl.)
1'12=( (P3-P.2)/2.0) * (VlNJ-V2)
W3=(P4V4-P4*V2)
W4=(l/(N2-1»* (P4*V4-15*VEXII)
*
DElb(l/ (Vl-V2) ) (-+ll-W2iW3ofW4)
FMEP=O .0025*NE
B!ElbIMEP-FMEP
'rolQ=716 .2* (VI.-V2) *B!EP/ (225*S'lRORS)
WE-(60/(6.28* (INERmI-INER'lL»)* ('ro~EXJ:.<W)
WRlTE(1,20) PN:;LLJ , Y , BmP,MFCV,MEcrl
20 FORM1\.T(lHO,'/ROLLJ=' ,F5.1,3X, 'MF=' ,F4.0,' 0/0' ,3X, 'BmP=' ,F7.4,
,3X,'MFCV=',F9.7,3X,'MfCP=',F9.7)
l\N3LLJ=J\N3LLJ-1
101 cmTINUE
NE=NFr400
99 <INTINUE
O\LL EXIT
END
121
The effect of the temp. T2 and 'pressure P2 in the second term'is not
so big. and the equation Can be linearized around the mean value of
P
3
= P
2
+C
2
MFCV ........ (4.7)
6
where C2 = 0.15271 10 in the model considered
P = P + C (T - T ) (4.8)
3 2 3 3 2
C
2
C = - = 0.07272
3 Cl
(4.9)
V2 V3
W = pdv + f pdv + f ~P. dv
c
Vinj V
inj
where ~P = the rise in the pressure due to the heat released in the system.
V
2
= ! pdv
V.1n j'
=
(4. lOa)
122
V3
and W2 = ! 6P.dv
V
inj
and volume during combustion, but the amount of work done can be approximated
P _.
3
V• •
l.nJ
Fig. (4.9)
1 I
W
2 = 3(P3-P 2 )· (V
inj
-V )
2 = 3(C2MFCV) (V -V )
inj 2 (4. lOb)
I 1
andW
c = n 1 - 1 (P 2V2 - PlVl ) + 3(C2MFCV) • (Vinj - V2 ) (4.11)
5
W = 48.006 PI - 0.4279 10 *MFCV + 0.509 105 V MFCV (4.12)
c • inj·
4.5.3 Computation of V
4
Ml+MF
=
Ml
123
Analogue Co.puter Si_ut.tion Chart
p1
Fig. ,(4.10)
C8
AVERAGE WORK
Cll -. 0.79578
crcuo
C9 ". + 716.2
•q
225 • STROKES
Hr
CH - 0.028161
+
...
-0.0025 (VI-V2)
Cl3 - 0.056322
HE
NE
--------{OoOO6\--------i.-------------------------
where Ml, MF are masses of air and fuel burned in one cycle.
Ml+MF
The ratio Ml can be considered equal = 1
but T4 = T3 + Cs MFCP so
3 + Cs
r
MFCP
V3
V = = V3 + T' Cs • MFCP
4 T3 ]-V3
3
V3 V2
but -T3 =
T2+C MFCV = constant because when MFCP has a value more than
I
zero it means that MFCV has the max. value, and the effect of T2 isnegiigible •
= V3 + C6 ·MFCP (4.13)
V
2
= C - = 1000 in the model considered.
S T3
V4 as is shown later.
Ps = [V f2
Ye: 'P4
but P P
3 and
=
4
T4
V = so
4 T"3
3
= [V3 T4(2
Ps Vexh T3 'P 3
= [v:~Jn2. [::(2, P3
12S
T4 = T3 + CS.MFCP
r
We have thus
+ Cs·MFCp'
2
Ps = [2f2 [T3
Vexh 'P 3 T3
n
2
=
[v:~r2'P3+ + Cs.MFCP]
T3
=
[v:~r2'P3' [1 + n 2·cS·MFCP J
T3 approximately
. [v:~r2'P3 + [[2-]"'
V exh
.p)"
",.c,.~
3T 1
This equation must be linearized, and only the most important variable
(notice that the model is not valid for small value of MF)
..
126
4.5.5 The Work done in the Expansion Stroke.
= =
= (4.15)
= =
From the simulation program (4.8) see Equ~tion (4.7) and (4.13)
6
= 0.1527.10
= 1000
6
= (P 2 + 0.1527 10 MFCV) (V
2
+ 1000 MFCP) =
6
= P2V2 + 1000 P2 MFCP + 0.1527 10 .V .MFCV +
2
It can be seen that the 1st and 3rd terms are linear.
The fourth term is non linear, but the value of MFCV may be considered
9 30
0.1527(10 ) (}6 0.0007) MFCP
3
= 89.08 (10 ) MFCP
127
3 3 3
39.5 10 MFCP + 89.08 10 MFCP = 128.58 i0 MFCP.
= = =
= C . (MFCP) (4.17)
7
3
where C = 128.58 10
7
= - PS· V
- exh)
- [v::~r2.P3.Vexh
=
0.""" • "'",]
n
= [n,:,] [,",. ,,,,., '0\ • [ j V
--
V
exh
2
2
.V h.0.1s27E6)MFCV +
ex
+ (128.58 10
3
- 7000 V h)MFCP- -V V: n .V ] J2
ex exh ·P 2
e h
128
=
= 3 3
[n2=1] [10.916 P1 + 191.52 10 MFCV + 54.436 10 MF(;..,]
3
W4 = 38.986 P1 + 684 10 MFCV + 194.41 103 MFCP
or
(4.17a)
where
C
8 = 38.986
3
C
9 = 684. 10
C 3
10 = 194.41 10
IMEP =
Nett work
Swept volume
n~-l
1
=
[ -
(4.19)
129
4.7 The Friction Mean Effective Pressure F:M.E.P.
Tq = (4.22)
716.2 (V -V )
Tq = l 2
225 STROKES (IMEP -- 0:0025 NE) (4.23)
-2. (NE) =
60
dt 2IT(INERTE + INERTL) (Tq - 0.005 NE) (4.24)
where INERTE, INERTL are the inertia of the rotating parts of the engine
From equations (4.12) (4.18), the summation of work done in one cycle
(two rev.) is
5
EW = - (48.006 PI - 0.4279-l0 MFCV + 0.509 105 V. . MFCV) +
- 1n]
3
+ (323 10 MFCP + 38.986 PI + 684 103 MFCV)
l30
or
3
LW = 9.02 PI + (726.79 10 - 50.9 10 3 V ) MFCV + 323 10 3 MFCP (4.25)
inj
This can be represented as in Fig. (4.10).
Notice that Cll, C12 have numerical values derived from the values
From the analogue simulation diagram Fig. (4.10) the output speed
NE is given by
C12
NE = S (Tq - C13.NE)
NE ~ (S + C12.Cl3) = C12 Tq
NE
=
C12
Tq S + C12.C13 (4.26)
NE =
Cll.C12
(S + Cll C13) (LW - C14.NE) (4.27)
and
C12
NE = (S + Cll.Cl3) Tq (4.28)
or Tq = S + Cll Cl3 NE
C12
NE =
2.7675
S + 0.1227 LW (4.31)
Substituting for the value of the constants Cll C12 in Equation (4.30)
Tq = 0.79578 (S + 0.04482) LW
(S + 0.1227) (4.32 )
131
Substituting for EW from equation (4.25) in equation (4.32) gives
0.79578 (5 + 0.04482) [ 3
Tq = (5 + 0.1227) - 9.02 P1 + (726.79 - 50.9 V ) 10 MFCV +
inj
3
-7.1779 (5 + 0.04482) (578.4 - 40.5 V ) 10 (5+0.04482)
Tq = inj
(5 + 0.1227) P1 + (5 + 0.1227) MFCV +
3
+ 257.04 10 (5 + 0.04482)MFCP
(5 + 0.1227) (4.33)
NE =
2.7675
(5 + 0.1227)
rL-9.02 P 3
1 + (726.79 - 50.9 Vinj ) 10 MFCV +
3
+ 323 10 MFCpJ
3
893.9 10
+ (5+0.1227) MFCP (4.34)
The transfer function matrix between outputs (Tq, NE) and inputs
(P1, MFCV, MFCP) can be formed from equations (4.33), (4.34) as follows:
where
.
g 1.)
.. in this (G) matrix is
24.963
(5+0.1227)
1·1119 (5+0.04482)
=- (5+0.1227) r-111'J ( )
3
893.9 10
= - '(";;5-+:--::0:"""."12:::-:2"7")
3
257.04 10 (5 + 0.04482)
(5 + 0.1227)
132
This transfer function matrix is valid for great ·v&lue of MF.
of (ANGLIJ), so the variables MFCV, MFCP are not the fuel input variables,
combustion starts PREXDC ---- (deg), which depends upon many factors
such as the engine speed and temp. It was modelled in the digital
The effect of temp. is not taken into account because the injection
of the fuel.
shown in Fig.(4.10)
the optimum BMEP (BAR) and find the corresponding· values of MFCV --- (kg),
and MFCP --- (kg) at various engine speeds and injection advances.
133
COMPUTER OUTPUT TABLE -(4.8a) OF PROGRAM NO. (4.8)
134
COMPUTER OUTPUT TABLE (4. 8b) OF PROGRAM NO. (4.8)
'1lIE INPUT AIR PRESSURE pI = 1.1 BAR AND RACK = 1.0.
135
· COMPUTER OUTPUTTABLE (4.8c) OF PROGRAM (4.8)
THE INPUT AIR PRESSURE pI = 1.3 BAR AND RACK = 0.8
OENGlNE SmED=2000.Rl'M EXLOIID= 112.6KGM
OlmLIJ= 29.0 MF= 80. % Bl£P= 5.3575 Mrol=0.0004958 MFCP=O .0000642
OlmLIJ= 28.0 MF= 80. % Bl£P= 5.7153 Mrol=0.0004567 MOCP=0.0000933
OlmLIJ= 27.0 MF= 80. % BffiP= 5.9334 Mrol=0.0004375 MFCP=0.0001225
OlmLIJ= 26.0 MF= 80. % Bl£P= 6.0132 Mrol=O.0004083 MOCP=0.0001517
OlmLIJ= 25.0 MF= 80. % BffiP= 5.9558 Mrol=0.0003792 MFCP=0.0001808
OlmLIJ= 24.0 MF= 80. % BffiP= 5.7625 Mrol=0.0003500 MFCP=0.0002100
OlmLIJ= 23.0 MF= 80. % BffiP= 5.4346 Mrol=0.0003208 MFCP=0.0002392
OlmLIJ= 22.0 MF= 80. % BffiP= 4.9734 Mrol=0.0002917 MOCP=0.0002683
OlmLIJ= 21.0 MF= 80. % BffiP= 4.3803 Mrol=0.0002625 MFCP=0.0002975
OlmLIJ= 20.0 MF= 80. % BffiP= 3.6567 Mrol=0.0002333 MOCP=0.0003267
OlmLIJ= 19.0 MF= 80. % BffiP= 2.8040 Mrol=0.0002042 MFCP=0.0003558
OENGlNE SmED=1600.Rl'M EXLQ.lID= 90.lKGM
OlmLIJ= 29.0 MF= 80. % BffiP= 5.6202 Mrol=0.0005658 MFCP-.0000058
OlmLIJ= 28.0 MF= 80. % BffiP= 6.2631 Mrol=0.0005367 MFCP=O .0000233
OlmLIJ= 27.0 MF= 80. % BffiP= 6.7654 Mrol=0.0005075 MFCP=0.0000525
OlmLIJ= 26.0 MF= 80. % BffiP= 7.1282 Mrol=O.0004783 MOCP=0.0000817
OlmLIJ= 25.0 MF= 80. % BffiP= 7.3527 Mrol=0.0004492 MFCP=0.0001108
OlmLIJ= 24.0 MF= 80. % BffiP= 7.4403 Mrol=0.0004200 MFCP=0.0001400
OlmLIJ= 23.0 MF= 80. % BffiP= 7.3921 Mrol=0.0003908 MFCP=0.0001692 .
OlmLIJ= 22.0 MF= 80. % BffiP= 7.2096 Mrol=0.0003617 MFCP=O.0001983
OlmLIJ= 21.0 MF= 80. % BffiP= 6.8941 Mrol=0.0003325 MFCP=0.0002275
OlmLIJ= 20.0 MF= 80. % BffiP= 6.4468 Mrol=0.0003033 MOCP=0.0002567
OlmLIJ= 19.0 MF= 80. % BffiP= 5.8693 Mrol=O.0002742 MFCP=0.0002858
OENGlNE SmED=1200.Rl'M EXLQ.lID= 67.6KGM
OlmLIJ= 29.0 MF= 80. % BffiP= 5.2807 Mrol=0.0006358 MFCP-.0000758
OlmLIJ= 28.0 MF= 80. % BffiP= 6.2088 Mrol=0.0006067 MOCP- .0000467
OlmLIJ= 27.0 MF= 80. % BffiP= 6.9951 Mrol=0.0005775 MFCP-.0000175
OlmLIJ= 26.0 MF= 80. O/OBffiP= 7.6410 Mrol=0.0005483 MFCP=0.00001l7
OlmLIJ= 25.0 MF= 80. % BffiP= 8.1474 Mrol=0.0005192 MFCP=0.0000408
OlmLIJ= 24.0 MF= 80. % BffiP= 8.5159 Mrol=0.0004900 MFCP=0.0000700
OlmLIJ= 23.0 MF= 80. % BffiP= 8.7475 Mrol=O.0004508 MFCP=0.0000992
OlmLIJ= 22.0 MF= 80. % BffiP= 8.8437 Mrol=0.0004317· MOCP=O .0001283
OllOOLIJ= 21.0 MF= 80. % BffiP= 8.8056 Mrol=O.0004025 MFCP=0.0001575
OlmLIJ= 20.0 MF= 80. % BffiP= 8.6347 Mrol=0.0003733 MOCP=0.0001867
OllOOLIJ= 19.0 MF= 80. % BffiP= 8.3323 Mrol=0.0003442 MFCP=0.0002158
OENGlNE SmED= 800.Rl'M EXLQ.lID= 45 .1KGM
OllOOLIJ= 29.0 MF= 80. % BffiP= 4.3391 Mrol=0.0007058 MFCP-.0001458
OllOOLIJ= 28.0 MF= 80. % BffiP= 5.5523 Mrol=0.0006767 MOCP-.0001167
OlmLIJ= 27.0 MF= 80. % Bl'EP= 6.6227 Mrol=0.0006475 MFCP-.0000875
OlmLIJ= 26.0 MF= 80. % Bl'EP= 7.5516 Mrol=O .0006183 MOCP- .0000583
OlmLIJ= 25.0 MF= 80. % Bl'EP= 8.3400 Mrol=0.0005892 MFCP-.0000292
OlmLIJ= 24.0 MF= 80. % BffiP= 8.9893 Mrol=O.0005600 MFCP=O .0000000
OllOOLIJ= 23.0 MF= 80. % BffiP= 9.5007 Mrol=0.0005308 MFCP=0.0000292
OllOOLIJ= 22.0 MF= 80. % BffiP= 9.8755 Mrol=0.0005017 MFCP=0.0000583
OlmLIJ= 21.0' MF= 80. % BffiP=10.1150 Mrol=0.0004725 MFCP=0.0000875
OlmLIJ= 20.0 MF= 80. % Bl£P=10.2204 Mrol=0.0004433 MOCP=0.0001167
OllOOLIJ= 19.0 MF= 80. % BffiP=10.1933 Mrol=O.0004142 MFCP=0.0001458
OENGlNE SIEED= 400.Rl'M EXLQ.lID= 22.5KGM
OllOOLIJ= 29.0 MF= 80. % BffiP= 2.7952 Mrol=0.0007758 MFCP-.0002158
OlmLIJ= 28.0 MF= 80. % Bl'EP= 4.2936 Mrol=0.0007467 MFCP-.000HI67
OllOOLIJ= 27.0 MF= 80. % Bl'EP= 5.64l!l Mrol=0.0007175 MFCP-.0001575
OlmLIJ= 26.0 MF=80~ % Bl'EP= 6.8600 Mrol=0.0006883 MFCP-.0001283
OlmLIJ= 25.0 MF= 80; % BffiP= 7.9304 Mrol=0.0006592 MFCP-.0000992
OlmLIJ= 24.0 MF= 80. % Bl'EP= 8.8606 Mrol=0.0006300 MFCP-.0000700
OlmLIJ= 23.0 MF= 80. 0/0' Bl£P= 9.6518 Mrol=0.0006008 MFCP-.0000408
OlmLIJ= 22.0 MF= 80. % Bl'EP=10.3052 Mrol=0.0005717 MFCP-.00001l7
OlmLIJ= 21.0 MF= 80. % Bl'EP=10.8221 Mrol=0.0005425 MFCP=0.0000175
OlmLIJ= 20.0 MF= 80. % Bl-EP=11.2040 Mrol=0.0005133 MFCP=0.0000457
OlmLIJ= 19.0 MF= SO. % Bl£P=11.4520 Mrol=0.0004842 MFCP=0.0000758
136
<X>MPUTER OUTPUT TABLE (4.8d) OF PROGRAM (4.8)
WE INPUT AIR PRESSURE Pl = 1.3 BAR AND RACK = 0.6.
.137
COMPUTER OUTPUT TABLE (4.8e). OF PROGRAM NO. (4.8)
THE INPUT AIR PRESSURE pl =
1.3 BAR AND RACK 0.4. =
OFNGlNE SfEED=2000.RPM EXLCW>= 112.6KGM
OllOOLLJ= 29.0 MFa 40. % Bl£lb-2.2504 MFC\T=0.0004958 MFClb-.0002158
OllOOLLJ= 28.0 MFa 40. % BMElb-1.3909 MFC\T=0.0004567 MOClb- .0001867
OllOOLLJ= 27.0 MF= 40. % Bl£lb-0.6709 MFC\T=0.0004375 MFClb-.0001575
OllOOLLJ= 26.0 MF= 40. % BMElb-0.0893 MroT=0.0004083 MOClb- .0001283
OllOOLLJ= 25.0 MFa 40. % BMEP= 0.3551 MFC\T=0.0003792 MFClb-.0000992
OllOOLLJ= 24.0 MF= 40. % BMEP= 0.6636 MroT=O.0003500 MFClb-.0000700
OllOOLLJ= 23.0 MF= 40. % BMEP= 0.8375 MFOl=0.0003208 MFClb- .0000408
OllOOLLJ= 22.0 MFa 40. % BMEP= 0.8781. MFC\T=0.0002917 MOClb- .0000117
OllOOLLJ= 21.0 MF= 40. % BMEP= 0.7869 MFOl=0.0002625 MFCP=O .0000175
OllOOLLJ= 20.0 MF= 40. % BMEP= 0.5651 MFOl=0.0002333 MOCP=0.0000457
OllOOLLJ= 19.0 MF= 40. % BMEP= 0.2142 MFC\T=0.0002042 MFCP=0.0000758
OFNGlNE SfEED=16 OO.RPM EXLCW>= 90.IKGM
OllOOLLJ= 29.0 MF= 40. % BMElb-3.1921 MFC\T=0.0005658 MFClb- .000285 8
OllOOLLJ= 28.0 MF= 40. % BMElb-2.0474 MFOl=0.0005367 MOClb-.0002567
OllOOLLJ= 27.0 MF= 40. % BMElb-1.0433 MroT=0.0005075 MFClb-.0002275
OIR;LLJ= 26.0 MFa 40. % BMElb-O.1787 MroT=O.0004783 MOClb- .0001983
OllOOLLJ= 25.0 MFa 40. % BMEP= 0.5476 MroT=0.0004492 MFClb-.0001692
OllOOLLJ= 24.0 MF= 40. % BMEP= 1.1370 MFC\T=0.0004200 MFClb- .0001400
OllOOLLJ= 23.0 MF= 40. % BMEP= 1.5907 MroT=O.0003908 MFClb-.0001108
OllOOLLJ= 22.0 MF= 40. % BMEP= 1.9100 MFC\T=0.0003617 MOClb-.0000817
OIR;LLJ= 21.0 MF= 40. % BMEP= 2.0963 MroT=0.0003325 MFClb-.0000525
OllOOLLJ= 20.0 MF= 40. % Bl£P= 2.1508 MFOl=0.0003033 MOClb-.0000233
OllOOLLJ= 19.0 MF= 40. % BMEP= 2.0751 MFC\T=0.0002742 MFCP=0.0000058
OFNGlNE SfEED=1200.RPM EXLCW>= 67.6KGM
OllOOLLJ= 29.0 MFa 40. % BMElb-4.7359 MFOl=0.0006358 MFCP=-.0003558
OIR;LLJ= 28.0 MF= 40. % BMElb-3.3061 MFOl=0.0006067 MOCP=-.0003267
OIR;LLJ= 27.0 MF= 40. % BMElb-2.0179 MFOl=0.0005775 MFCP=-.0002975
OllOOLLJ= 26.0 MF= 40. % BMElb-O.8703 MroT=0.0005483 MOCP=- .0002683
OllOOLLJ= 25.0 MF= 40. % BMEP= 0.1380 MFOl=0.0005192 MFCP=-.0002392
OIR;LLJ= 24.0 MF= 40. % BMEP= 1.0083 MroT=0.0004900 MFCP=-.0002100
OllOOLLJ= 23.0 MF= 40. % BMEP= 1.7417 MroT=O.00045 08 MFClb-.0001808
OllOOLLJ= 22.0 MF= 40. % BMEP= 2.3397 MroT=O.0004317 MOCP=-.0001517
OllOOLLJ= 21.0 MF= 40. % BMEP= 2.8034 MFOl=O.0004025 MFCP=- .0001225
OllOOLLJ= 20.0 MF= 40. % BMEP= 3.1344 MFOl=0.0003733 MOCP=- .0000933
OllOOLLJ= 19.0 MF= 40. % BMEP= 3.3338 MFC\T=0.0003442 MFClb-.OO00642
OFNGlNE SmED= 800.RPM EXLCW>= 45 .IKGM
OllOOLLJ= 29.0 MF= 40. % BMElb-6.8820 MFC\T=O .0007058 MFClb-.0004258
OllOOLLJ= 28.0 MF= 40. % Bl£lb-5.1670 MFOl=0.0006767 MOClb-.0003967
OllOOLLJ= 27.0 MF= 40. % Bl£lb-3.5947 MroT=0.0006475 MFClb-.0003675
OllOOLLJ= 26.0 MF= 40. % Bl£lb-2.1640 MFC\T=0.0006183 MOClb- .0003383
OllOOLLJ= 25.0 MF= 40. % BMElb-0.8737 MroT=0.0005892 MFClb- .0003 092
OllOOLLJ= 24.0 MF= 40. % BMEP= 0.2774 MFOl=0.00056 00 MFClb- .0002800
OllOOLLJ= 23.0 MF= 40. % Bl£P= 1.2906 MroT=O.0005308 MFClb-.0002508
OllOOLLJ'" 22.0 MF= 40. % BMEP= 2.1672 MFOl=0.0005017 MOClb-.0002217
OllOOL1J= 21.0 MF= 40. % BMEP= 2.9085 MFU/=0.0004725 MFClb-.0001925
llOOLLJ'" 20.0 MF= 40. % BMEP= 3.5157 MroJ=0.0004433 MOClb- .0001633
OllOOLLJ= 19.0 MF= 40. % BMEP= 3.9904 MFU/=0.0004142 MFClb-.0001342
FNGlNE SfEED= 400.RPM EXLCW>= 22.5&>M
llOOLLJ'" 29.0 MF= 40. % BMElb-9.6302 MFC\T=0.0007758 MFClb-.0004958
llOOLLJ= 28.0 MFa 40. % B/-Elb-7.6300 MFU/=0.0007467 MOClb-.0004567
llOOLLJ= 27.0 MF= 40. % 8I£lb-5.7736 MroJ"'0.0007175 MFClb-.0004375
llOOLLJ= 26.0 MF= 40. % BMElb-4.0599 MEO/=0.0006883 MOClb- .0004083
llOOLLJ= 25.0 MF= 40. % BMElb-2.4877 MroJ=0.0006592 MFClb-.0003792
llOOLLJ= 24.0 MF= 40. % 8I£lb-1.0557 MKV=0.0006300 MFClb-.0003500
llOOLLJ= 23.0 MFa 40. % BMEP= 0.2373 MroT=0.0006008 MFClb-.0003208
llOOLLJ= 22.0 MF=4O •. % Bl£P= 1.3925 MFOl=0.0005717 MOClb-.0002917
llOOLLJ= 21.0 MF= 40. % BMEP= 2.4113 MroT"'0.0005425 MFClb-.0002625
llOOLLJ= 20.0 MF= 40. % Bl£P= 3.2949 MEO/=0.0005133 MEClb-.0002333
llOOLLJ= 19.0 MF= 40. % Bl£P= 4.0448 MEO/"'0.0004842 MFClb- .0002042
138
The output is illustrated in the simulation program output tables
(4.8a) to (4. Ba). The initial conditions and the load are taken to
the optimum BMEP occurs at different injection advance when the speed
The input pressure depends upon the method of charging of the engine
variables are:
-28~Anglij~-14 deg. 1
400~NE~2000 rpm (4.37)
O. B~Pl~1. 5 BAR
volume MFCV is constant at the optimum conditions and constant load. Its
Within the ranges indicated above the volume of the combustion chamber
V
inj
V2 + (V1MA~ + V2) (4.333 - cos ~ - 3.3331'1- 0.09 sin2~ )=
= 1. 7455BLT (4.39)
4.13 Obtaining the Relation Between Angle of Injection and MFCV, MFCP
The idea is to eliminate MFCV, MFCP from equation (4.25) and to use the
139
(4.36) (4.39) giving:
3 3
Ew = - 9.02 PI + [726.79 10 - 50.9 10 (1.74558)] (PRECMB.~~) +
3
+ 323 10 (MP - MPCV)
3 0.0007 3
EW = -9.02 PI + 637.94 10 (-ANGLIJ + 0.006 NE) 2(24) + 323 10 •
with
MP = .0007 kg/cycle
3
EW = -9.02 PI - 4.593 ANGLIJ + 323 10 MP + 0.02766 NE (4.40)
2.7675 3
NE = (5 + 0.1227) (- 9.02 PI 4.593 ANGLIJ + 323 10 MF + 0.02756NE)
2.7675(323 10 3 )
+ (S + 0.1227) MF
3
NE = _ 24.963 1 12.711 893.9 10
(S + 0.04644) P - (S + 0.04644) ANGLIJ + (S + 0.04644) MF (4.41)
Notice that due to the inherent feedback of engine speed (0.006 NE,
see the computer simulation chart Fig. (4. 10) the term (S+O.1227) is
changed to (S+0.04644).
0.79578(S+0.04482)
Tq = (S + 0.1227) ( -9.02 PI - 4.593 ANGLIJ +
3
+ 323 10 MP + 0.02756 NE)
140
7.l779(S+0.04402) PI _ 3.655(S+0.04402)
1"'1 = IINGLIJ +
(S+0.04644) (StO. 04(44)
3
257.04 10 (S+0.04402)
+ MF (4.42)
(S+0.04644)
It is clear now that the pole has been shifted from (S+0.1227)
12.711
(S+0.04644) - 3.655 (S+0.04482)
(S + 0.04644)
It can be seen that the transfer function matrix has two dependent
columns, because the output speed~f the engine depends u?On t~e output
torque for a given load.
The first row is negative, because the increase in the input pressure
The second row has very high gain to transform a small amount of fuel
negative.
many difficulties arise, and the conclusion must be that the linear model
can not be realized with any great accuracy. However, the overall error
is only about 5 to 10\ as can be seen from the results of the nonlinear
141
Although the linear model is not accurate, it still has acceptable
it is very small.
test.
follow.
142
CHAPTER 5
MATRIX FOR MULTI INPUT MULTI OUTPUT (MIMO) SYSTEMS USING LOOPS THEORY
5.1 INTRODUCTION
Signal flow graphs and matrix methods nave been used for some
71
considerable time in the fi'eld of network analysis ,40, referring
9-28
chiefly to electrical networks, and in the area of control In
the latter field, much thought has been given in recent years to the
and feedback paths. The instructions thus spelled out would in fact
very true, but it is a step forward if we can state clearly just what
this complex relationship i'9, and that is what this chapter sets out
following chapter will allow the reader to see clearly why this should
be so,
143
-------------------------- .........
5.2 PATHS, LOOPS, AND MATRICES
transfer function matrix, which describes the way in which all the
Yl gll g21
Y3 g13 g23
L
To resolve the conflict and still maintain the notation, which
[:1 Y2 Y3] = h r
x2J gll
(5.2)
1921
Yet another form of description for the flow-graph of Fig. (5.1)
leaving a particular node are found in the appropriate row while the
the connection matrix appears to offer no more than, say, the matrix
of eqn. (5.2), but when we add to the plant of Fig. (5.1) a compensating
144
X·~------~rJ~~ ______~~________~
Xl x Yl Y2 Y3
2
Yl 0 0 0 0 0
Y2 0 0 0 0 0
Y3 0 0 0 0 0
145
feedback network as in Fig. (5.3) we obtain the connection matrix of
branches of Fig. (5.3) that have a gain of unity we may observe that
Yl = ul ' Y2 =u2 ' etc., and reduce the unwieldy matrix of Fig. (5.4)
with one another, and likewise for the output quantities, so the shape
which are null, along the diagonal, and the other two which are identical
with the transfer function matrices for the plant, G(s), and the feed-
expressed in terms of G(sl and HCsI and their minors which does not
146
Fig.(S.3) A compensated feedback network
147
"'1 z2 xl x Yl Y2 Y3 u u u v2
2 1 2 3 vI
zl 0 0 1 0 0 0 0 0 0 0 0 0
z2 0 0 0 1 0 0 0 0 0 0 0 0
xl 0 0 0 0 0
gll g12 g13 0 0 0 0
x 0 0 0 0
2 g21 g22 g23 0 0 0 0 0
Y1 0 0 0 0 0 0 0 1
.... 0 0 0 0
""
CD
Y2 · 0 0 0 0 0 0 0 0 1 0 0 0
Y3 0 0 0 0 0 0 0 0 0 1 0 0
u 0 0 0 0 0 0 0
1 0 0 0 hll h12
u 0 0 0 0 0 0 0
2 0 0 0 h21 h22
u 0 0 0 0 0 0 0 0
3 0 0 h31 h32
vI 0 0 1 0 0 0 0 0 0 0 0 0
v 0 0 0 1 0 0
2 0 0 0 0 0 0
Fig. (5.41 The connection matrix for the diagram of Fig. (5.31
Xl x YI Y2 Y3
2
x 0 0 gn
2 g22 g23
YI h11 hl2 0 0 0
Y2 h21 h22 0 o. 0
Y3 h31 h32 0 0 0
H
11
22
0:.=-----.. . _____ <<<'-___-=="""
G
11
•
149
Describing an m-input n-output system henceforth as an m x n system,
=
(5.3)
so that
g11 (1-g22 h 22) + g12 h 22 g 21
tll = ~
g11 - h 22 (det G}
1 + detGH - trGH (5.5)
150
In a similar manner one may derive
g12 + h 12 (det G)
t12 =1 + detGH trGH (5.6)
g22 - hU (detG)
t22 = 1 + detGH - trGH (5.8)
case, the two ports of (H) that are not directly connected to the input
W=
etc.
(Note that W in this context need not have the shape of Fig. (5.5)
are derived from the n rows and n columns in such a way that no
151
(iil The sum of the terms wd that constitutes the determinant covers
For example w13 w21 w32 is the gain of a single loop Fig. (5.7),'
while wll w32 w23 is the gain of two disjoint loops, Fig. (5.8).
051 Det W gives the surn of the gains of all the loops containing
n nodes, together with the surn of the gain products of those
0 w w13 w
12 14
0 0 w w
W =
23 24
0 w 0 w
32 34
1<41 0 0 0
whence
= + +
152
I / I \
I
/ \
\
/
I
I
w
\
,,
I
3 2
W
32
Fig. (5.7) A 3-node system w
wi th one loop 32
w13 w w Fig. (5.8) A 3-node system
21 32 incorporating two
loops
wll ' w32 w23
153
The terms wl and w2 are 4-node loops, while w3 is the product
(cl' If one node is removed from the flow graph, which means
one row and column from W, the det. of the reduced matrix,
give the gains of all the loops and disjoint groups with
n - 1 nodes.
wi'll giVe tlie sum of all the loop gains together with the
from the graph. all inputs to node i and all outputs from
node j. Tlie resulting minor (the det of the reduced matrix)
0 0 g11 g12
0 0 g2l g22
W=
hll h12 0 0
h2l h22 0 0
154
Removing all inputs to node 1 and outputs from node 3 means
,0 gll g12
W* = g21 g22
l:22
° °
and the minor by
multiplied by the gain of the loop that does 'not touch the
path, and the second term is the gain of the other path
(e[ To reveal every path from node i to node j with its associated
155
1 If
__- r 1 .
156
row j, and the nodes in all possible groups of ~ always
0 0 w 0
13
0 0 0 w
W= 24
w w 0 0
3l 32
w w w 0
41 42 43
and let us discover the paths and disjoint loops from node 4
to node 2.
o o
W*= (eliminating col.4, row 2)
~=
1
[.W32
w
42 :J (eliminating col.4,
row 2, plus col. 1, row 1)
W* =
3
[:41 :J (eliminating col. 4,
row 2, plus col.3, row 3.)
det WJ: = w w
32 43
det W* = 0
3
det Wi,3 = w
42
157
It is easily verified that these terms give all the necessary
where required.
transfer function.
(which is the shape of Fig. (5.111 with the addition of node 5) and
output respectively.
are four loops in the diagram which may be described in terms of the
nodes comprising them as: L13, L24, L243, L1324. Since the first two
do not touch,
5 4 1
Fig. (5.12)
3 2
158
There are three paths from node 5 to node 2: p542, P5432, P54132,
one of which (P542! has the disjoint loop L13. Hence the numerator
matrix
:;
0 0 w13 0
0 0 0 w
24
W = w w 0 0 (5.11 )
31 32
w
0
41
w
0
42
w
0
43
0
w
54
:j
we first note that because column S is zero det Wand all minors that
include node S are sero. In fact the only submatrices with nonzero
determinants are
-0 0 W13 0
0 0 0 w
24
Ws = deleting node S
w w 0 0
31 32
w w w
41 42 43 0
J
W
5,3,1 deleting nodes 5, 3 and 1
159
W deleting nodes 5, 4 and 2
5,4,2
where
q5 q5 1 q
(-11 . det W5 + (-11 ' ·det W5 1 + (-11 5,3'~det W +
, 5,3,1
q
+ (-11 5,4'~det W
5,4,2
(5.12)
So
more intricate because we are now concerned with loops and paths, and
p = r + c + q - 1 (5.14)
160
0 0 W13 0
w w 0 0
31 32
W* =
w w w 0
41 42 43
0 0 0 w
54
set W* =
det Wl: =
det Wl:,3 =
are the only nonzero minors. For each of these minors the value
Po = 4 + 2 + 4 - 1 = 9 (from W*)
P1 3 + 1 + 3 - 1 =6 (from Wi)
P1 ,3 2 + 1 + 2 - 1 =4 (from Wl:,3)
=
and evaluate the numerator of the transfer function from the expression
Po P1 P
S* (-11 . det W* + (-1) .det W* + (-1) l,3.det W*
n 1 1,3
161
5.3.4 Summary .
given by
S*
n
= = (5.16)
matrix W*.
having the shape of Fig. (5.51: all the forward gains appear in one
and CS.1511 in terms of (G1 and CHI, removing the need to form the
W-matrix.
real cost by omitting the input nodeCsl and assuming that all inputs
162
5.5 Calculation of an Element t .. of the Total Transfer Function Matrix
1)
transfer function matrix (T), given the plant matrix (G) and the
S*
n
t .. =
1) 1 + Sd
resulting minor is' zero, and in fact all submatrices will be singular
unless they are formed By removing an equal number of input and output
nodes'.
th
i
Remov'ng the a th.rnput an dtheb output node f rom (W)has t h e
th th
effect of removing the a row from (G) and the a column from (H)
and the b th column from (Cl: and the b th row from CHI. If we delete
row a, column I:i from (C) the resulting suhmatrix is denoted (G) a
b
163
and it follows from the observation above that when we form (G)~ we
must also form (H)b the H-matrix less row b, column a. It is found
a'
furthermore that
and
det (Wl =
ab
th
where CWl ab represents (w) less the a input node and the b th output
a
node, and q (or qab) is the dimension of (G) (or (G)b).
qab a
S = E(-l) .det (G)b.det (5.17)
d
where the summation is taken over all submatrices of (G) and (H)
exceptions.
a
Omit row a column b from (G) to form (G)b. Omit column a row b
b
from (H) to form (H) a where a, b will take every possible combination
164
(a) (b) order of det. (q)
Group 0 0 0 3
Group I 1 1 2
1 2 2
1 3 2
2 1 2
2 2 2
2 3 2
3 1 2
3 2 2
3 3 2
Group IX 12 12 1
12 13 1
12 23 1
13 12 1
13 13 1
13 23 1
23 12 1
23 13 1
23 23 1
a Hb
Rows and columns to be omitted in the formation of Gb'
a
TABLE (5.13)
165
total number of entries NCO over all the groups is given by:
r=(n-l) [ n! j2
NCO = l: ' ( _ ) ' when forming the denominator. (5.18)
r-
=0 .r. n r .
3'
NCO = [ O! (3!)]
12 r 3'
+ li!(2n]
12
+
[3' 12
2!(i!)] = 19 entries
(5.171, Le.
the minor.
In the given example when calculating t , the minors of (G) and
23
corresponding (H) are:
I
11, ~.
2
166
2
G =
2 l'n 'Bj
g31 g33
H2 =
2 rn OBj
h31 h33
r"
h22
3
G =
1
['n 'B] HI =
3
g22 g23 h31 h32.
on]
3
G =
2
r"
g21
'B]
g23
2
H =
3
r" h31 h32
on]
3
G =
3 ['" g21
'n]
g22
a33 =
r
u
h21 h22
12 12
G
12
= g33 a12 = h33
12 H13
G = g32 = h23
13 12
12 23
G = g31 H12 =
23 h13
13 H12
G = g23 =
12 13 h32
G13 = H13 =
13 g22 . 13 h22
13 H23
G = g21 = h12
23 13
23 H12
G = g13 =
12 23 h31
23 H13
G = g12
13 23 h21
23 23
G = g11 H h11
23 23
167
Now the denominator term Sd is
9 11
-
9 12 913 h11 h12 h13
3
Sd = (-1) 9 21 9 22 9 23 h21 h22 h23
9 31 932 9 33 h31 h32 h33.
+(_1)2
r" '
932
'j ~~' ~'j [~" ' 'j
9 33 h32
l'' "D]
h33
(-1) 2
9 31 9 33 h32 h33
+ (-1) 2
r ' '
9 31
j ~"u ""j
932 r,'" '
h22
uj r"
h23
",,](-1) 2
932 9 33 h31 h33
+(_1)2
r" r" r" '''j. r" ""j
931
' 'j "Uj
9 33 h31 h33
(-1) 2
9 31 932 h21 h23
+(_1)2
[~" ' 'j ["" ""j r" ' 'j r "')
9 22 9 23 h31 h32
(-1) 2
9 21 9 23 h31 h32
+(_1)2
1'"
9
,
21
,]. r" :,,]
922 21 22
1
+(-1) (933.h33 + 9 2h23 + 931h13 + 923 h 32 +
3 922h22 +
168
S.S.2 The Numerator Term Sn of the Element t
ij
it will be seen that this process leaves (G) untouched but deletes
a row and column from (HI. The process of deleting from W* all the
nodes, except those designated input and output, in all possible ways
has its counterpart in terms of (G) and (H) in the following process:
Take the list given above of rows and columns to be deleted from
-
(Gl and (HI and modify it to form (G) a by removing all lines that
b
involve either ~ i or column j. If, for instance, in the 3 x 3
example quoted above, the input were to node 2 and the output from
a = 2 or b:= 3.
Group 0 0 0
Group I I I
1 2
3 1
3 2
Ib9
The submatrices of (G) are then formed from this list by the same
of (HI which is formed as in the case of (G) but row i and column j
0,2 0,3
1,2 1,3
1,2 2,3
3,2 1,3
3,2 2,3
1,3,2 1,2,3
Table (5.151
The submatrices of (H). are then formed from this list by the
p--
. a • det (H) b
S*
n
= gij +"~ (-1'L ab. det (G) (5.20)
ii a
or if
p--
S !;(-l)
ab • det (G) a ii
• det(H)_
n (5.21)
b a
then
S* = gOj + Sn (5.22)
n 1·
170
In the case of group 0 in table (5.14), there is only one
=
entry which corresponds to no deletions.
In general (when (G) and (H) are of dimension n x n), the total
r= (~-l) 1 (n-l) : 12
NCN = r:o r:(n-r-l): when forming the numerator (5.23)
NCN = 2: 12 1 2: 12
+ I!.l!
1 2: 12 = 6 entries
1O!.2! + 2!.O!
Sn = 1: (-l)Pdet(G)a.det(H)a
all minors .b b
The value of parameter p is deduced from (G)a according to
b
p = r+C+q-l (5.24)
-
where q is the order of (G) a
b
r, C depend upon the position of the pivot, which is the element
a above by the element
gij of the matrix (G)_ and 1, j are defined as
b
t of the total transfer function matrix (T).
ij
a
r, C are the row and column of (G) in which the pivot gij lies.
b
In the given example when calculating t 23 the basic matrices are:
gll
(G)a = CJ21
b
g31
b
The degree of matrix (H) is always one less than the degree of
a a
(G) matrix.
b
Thus (G) ~ and the corresponding (H) ~ are: (notice that the pivot is
b a
underlined)
and H~ = [:::
P = 3+2+3 -1 = 7
P = 1+2+2-1 = 4.
P = 1+2+2-1 = 4.
P = 2+2+2-1 = 5.
P = 2+2+2-1 = 5.
12
and H13= . o.
172
Now the element t of the total transfer function matrix can be
23
calculated according to equation (5.16) and (5.22)
gij + Sn
t'1. j .
= (5.27)
1 + Sd
. Le.
g23 + Sn
t =
23 1 + Sd
5.6 Digital program for calculation of total transfer function matrix (T)
173
A computer program was therefore written to produce the necessary submatrices
so that the reader has only to sum .the necessary matrix products to form
Sn and Sd· See appendix (4.·1) .'ll1e requirerl input data for the program is simply
e.g. NN =3
- The position I, J of the total transfer function element t
ij
e.g. I = 2
J =3
The output of the program will be sets of corresponding minors of
(G) and (H), and it is stated that each minor product must be multiplied
by either'. (+1) or (-1) when performing the summation. See comt>uter output table (5.'. .1)
number (94 ) which helps performing the hand work and saves enormous time
and mistakes.
Fig. (5.16)
174
described above and the program used previously modified accordingly
(see Appendix 4.2 also see computer output table 5.2 at the end of
this chapter.
and feedback matrix (H) which depends on obtaining the total transfer
may be obtained.
th
* In general when the denominator of t is of n order, n-
ij
conditions must be satisfied (i.e. according to the coefficients of
n-l
S ,and the specification of the numerator constant adds
h where
ji
NH = n+l.
m m
NF = 1: 1: Rij
i=l j=l
the rank of the final T-matrix, and in this case no F-matrix can be
obtained.
• It is not always possible when fi' and h" are constant to satisfy
. ) )1
• In general when h
ji
is a constant it can not appear in the numerator,
This method opens the way to the design of feedback controllers for
Appendix (5) lay down the rules, and illustrate its capability and
The reader need not know much about the theory or the procedures
matrix when dealing with poles, because the denominator of each term
contains all the poles of the system, see Appendix (51 ,but it is
When calculating SN, SO the hand work needed is· not so big as most
176
TO shift n poles it is necessary to consider at least n feedback
Though often the (F) matrix is unable to place the zeroes exactly
where they are required, but a satisfactory value for (F) can be
stability.
in the system.
177
COMPUTER OUTPUT TABLE (5.1) TO CALCULATE THE ELEMENT t
23
KlLTIH.IED BY -1
DE'lfmIImNT OF
FGl.l ROl.2 ROl.3
R02.1 R02.2 R02.3
ro3.1 ro3.2 R03.3
KlLTIH.IED BY ':mE IE'IERMlNllNT OF
Hl.l m.3
m.l m.3
MDLTIH.IED BY +l
DE'lmMINllNT OF
002.2 R02.3
ro3.2 ro3.3
KlLTIH.IED BY ':mE lE'l'mMINllNT OF
H2.3
K1LTIR.IED BY +l
DE'lfmIImNT OF
FG2.l R02.3
ro3.l ro3.3
KlL'l'IH.IED BY '!HE lE'l'mMINllNT OF
Hl.3
KJLTIR.IED BY -1
DE'lmMINllNT OF
FGl.2 ROl.3
FG2.2 FG2.3
MDLTIH.IED BY ':mE lE'l'mMINllNT OF
H2.l
/oUL'l'IH.IED BY -1
DE'lmMINllNT OF
FGl.l R01.3
R02.l FG2.3
MUL'l'IH. IED BY 'llIE DETERMlNllNT OF
l7ti
lULTIH.mD BY 'ruE IE'reRMlNlINT OF
Hl.l
IULTIR.omD BY -1
DE'llRMINJ\NT OF
FGl.l FGl.2 FGl.3
FG2.1 FG2.2 FG2.3
FG).1 FG3.2 FG3.3
KJLTIR.mD BY 'ruE rE'IERMlNlINT OF
Hl.l HI.2 HI.3
H2.1 H2.2 H2.3
ID.l ID.2 ID.3
KJLTIH.mD BY +1
DE'lmMINlINT OF
FG2.2 FG2.3
FG).2 FG).3
KJLTIR.mD BY 'ruE IETERMINlINT OF
H2.2 H2.3
ID.2 ID.3
KJLTIH.IED BY +l
DE'llRMINJ\NT OF
EG2.1 FG2.3
FG).1 FG).3
KJLTIH.IED BY 'ruE IETERMINlINT OF
ffi.2 HI.3
ID.2 ID.3
KJLTIR.IED BY +l
DE'lmMINlINT OF
FG2.1 FG2.2
FG).1 FG3.2
IULTIR.IED BY 'ruE IETERMINlINT OF
Hl.2 HI.3
H2.2 H2.3
KJLTIR.mD BY +1
DE'lERMINlINT OF
FGl.2 FGl.3
FG3.2 FG3.3
MULTIR.IED BY '!HE IETERMINlINT OF
H2.1 H2.3
179
H2.1 H2.3
H3.1 H3.3
KJLTIH.IED BY +1
DE'lmMINJ\NT OF
R;l.l R;1.3
R;3.1 R;3.3
KJLTIH.IED BY 'mE IETERMlNllNTOF
Hl.I m.3
H3.1 H3.3
KJLTIH.IED BY +1
DE'lnU«NPNT OF
R;l.l R;1.2
R;3.1 R;3.2
!{)LTIH.IED BY 'mE I:ETERMlNllNT OF
Hl.1 m.3
H2.1 H2.3
KJLTIH.IED BY +1
DE'lmMINJ\NT OF
R;l.2 R;1.3
R;2.2 R;2.3
!{)LTIH.IED BY +1
DE'lmMINJ\NT OF
R;l.l R;1.3
R;2.1 R;2.3
!{)LTIH.IED BY +1
DE'lfRMINPNT OF
R;l.l R;1.2
R>2.1 R;2.2
!{)LTIH..IED BY -1
DE'lffiMINllNT OF
180
lE'ImMINNlT OF
FG3.3
KJLTIH.IED BY WE IE'lERMlNNlT OF
m.3
KJLTIH.IED BY -1
DE'lmMINllNT OF
FG3.2
MJLTIH.IED BY WE IE'lERMlNllNT OF
H2.3
KJLTIH.IED BY -1
DE'nRMINllNT OF
Rn.I
MJLTIR.IED BY WE IE'lERMlNllNT OF
HI.3
KJLTIR.IED BY -1
DE'nRMINllNT OF
FG2.3
KJLTIR.IED BY WE IETERMINllNT OF
m.2
MJLTIR.IED BY -1
DE'lmMINllNT OF.
FG2.2
MJLTIR.IED BY WE IETERMINllNT OF
H2.2
KJLTIH.IED BY -1
DE'lmMINllNT OF
FG2.1
MJLTIR.IED BY WE IE'lmMlNllNT OF
H1.2
MJLTIH.IED BY -1
DE'illRMINllNT OF
roI.3
KlLTIH.IED BY WE IE'lmMlNNlT OF
H3.1
MJLTIH.IED BY -1
181
KJLTIR.mo BY -1
OE'lmMINllNT OF
R;1.2
MJLTIR.mO BY 'mE ~l1NT OF
H2.l
KJLTIR.mo BY -1·
OE'lfRMINllNT OF
FGl.l
KJLTIR.mo BY 'mE ~l1NT OF
HI.l
COMPUTER OUTPUT TABLE (5.2) TO CALCULATE '!HE ELEMENT t WHEN FORWARD
23
COMPENSATION MATRIX (F) IS USED.
KJLTIH.IED BY -1
DE'.lmMINlNr OF
Gl.l Gl.2 Gl.3
G2.l G2.2 G2.3
G3.l G3.2 G3.3
MDLTIH.IED BY +1
DE'lmMINmT OF
G2.2 G2.3
G3.2 G3.3
KJLTIH.IED BY mE lE'IERMINlINT OF
H2.3
KJLTIH.IED BY +1
DE'.lmMINlNr OF
G2.l G2.3
G3.l G3.3
KJLTIH.IED BY mE IE'IERMINlINT OF
m.3
MDLTIH.IED BY -1
DE'nRMINmT OF
Gl.2 Gl.3
G2.2 G2.3
MJLTIH.IED BY -1
DE'lmMINl'lNT OF
Gl.l Gl.3
G2.l G2.3
183
KJLTIILIED BY 'mE IE'ffiRMlNllNT OF
H1.1
KlLTIILIED BY -1
DE'lnUITNllNT OF
Gl.1 Gl.2 Gl.3
G2.1 G2.2 G2.3
G3.1 G3.2 G3.3
MULTIILIED BY +1
DE'lnUITNllNT OF
G2.2 G2.3
G3.2 G3.3
MULTIILIED BY +1
DE'lnUITNllNT OF
G2.1 G2.3
G3.1 G3.3
MULTIILIED BY +1
DE'lfRMINllNT OF
G2.1 G2.2
G3.1 G3.2
MULTIILIED BY +1
DE'lERMINilNT OF
Gl.2 Gl.3
G3.2 G3.3
184
H2.l H2.3
1D.1 1D.3
KJLTIH.mo BY +1
OE'mOONllNT OF
Gl.l Gl.3
G3.l G3.3
KJLTIH.mo BY +1
OE'mOONllNT OF
Gl.l Gl.2
G3.l G3.2
KlLTIH.mo BY +1
OE'mOONllNT OF
Gl.2 Gl.3
G2.2 G2.3
KJLTIH.mo BY +1
OE'mRMINllNT OF
Gl.I Gl.3
G2.1 G2.3
IULTIR.mO BY +1
OE'mRMINllNT OF
Gl.I Gl.2
G2.1 G2.2
IULTIR.mO BY -1
DETI':RMINl'NT OF
185
IE'lERM!NllNT OF
Gl.3
IDLTIR.IED BY '!HE IE'lERMlNllNT OF
H3.3
MULTIR.IED BY -1
DE'lmMINllNT OF
Gl.2
MJLTIR.IED BY '!HE lE'lERMlNllNT OF
H2.3
MULTIR.IED BY -1
DE'lERMINllNT OF
Gl.1
MJLTIR.IED BY '!HE lE'lERMlNllNT OF
H1.3
MULTIR.IED BY -1
DE'lERMINllNT OF
G2.3
IDLTIR.IED BY -1
DE'lEOONllNT OF
G2.2
MJLTIR.IED BY '!HE lE'.IERMlNllNT OF
H2.2
KJLTIR.IED BY -1
DE'IERMINllNT OF
G2.1
MJLTIR.IED BY '!HE lE'lERMlNllNT OF
H1.2
KJLTIR.IED BY -1
DE'IERMINllNT OF
Gl.3
. l.uLTIR.IED BY '!HE IETERMllillNT OF
ID.l
MJLTIR.IED BY -1
186
l-llLTIPLIED BY -1
DE'lERMlNliNT OF
Gl.2
l-llLTIPLIED BY THE IEiERMlNliNT OF
H2.1
l-llLTIPLIED BY -1
DE'lERMlNliNT OF
Gl.1
l-llLTIPLIED BY THE IEiERMlNliNT OF
Hl.1
18.1
CHAPTER 6.
6.1 Introduction
The diesel plant is described by the matrix (G), see Chapter (~),
equation ( 4.43).
And the signal flow graph of the system is shown in Fig. (6.1)
PI
NR.
MF
g22
T
q
ANGLIJ
.Fig. (6.1)
188
-VI
-001
PI
Fi. (6.lA)
AVERAGE WORK
CYC:U: ell a, 0.79578
716.2
•q
LW
225 .. STROKES
-0.0025 (VI-Vl)
MGLIJ NE
-0.005 (VI-V2)
o
e13. 0.056322
- W2 NE
NE
--------~o.-~--------------------------------------------~
Let the feedback link connect the output sp~ed with the input fuel
NE
MF
ANGLIJ
Fig. (6.2)
In this case the feedback matrix (H) is
(m= [ :
h12
0 :] (6.2)
TO achieve this constant steady state gain in the engine, the total
using the theory illustrated in Chapter 5, see pro\1ram output table (6.1).
g21 + SN
=
1 + SO (6.3)
In this case
r" r
g12 g13
h21 23
SN = h ]
g2l g22 g23
h3l h33
g3l g32 g33
[ g21 g23] •
h + [ g2l g22] • h23
+ - g3l g33 33
g3l g32
190
PROGRAM OUTPUT TABLE (6.1) TO CALCULATE t21
HJLTmIED BY -1
DE'lEOONllNT OF
Gl.1 Gl.2 Gl.3
G2.1 G2.2 G2.3
G3.1 G3.2 G3.3
MDLTmIED BY -1
DE'lfmIINllNT OF
G2.1 G2.3
G3.1 G3.3
mLTIR.mD BY -1
DE'lfmIINliNT OF
. G2.1 G2.2
G3.1 G3.2
MDLTIILmD BY +1
DE'lERMINllNT OF
Gl.1 Gl.3
G2.1 G2.3
MDLTIILmD BY mE IETERMlNliNT OF
83.1
191
KJLTIR.IED BY +1
DEifRMINJINT OF
Gl.1 Gl..2
G2.1 G2.2
KJLTIPLIED BY '!HE IE'lERMlNJINT OF
H2.1
KJLTIR.IED BY -1
DEifRMINJINT OF
Gl.1 Gl..2 Gl..3
G2.1 G2.2 G2.3
G3.1 G3.2 G3.3
mLTIPLIED BY +1
DE'lERMINllNT OF
. G2.2 G2.3
G3.2 G3.3
mLTIPLIED BY +1
DE'lEUIlNJINT OF
G2.1 G2.3
G3.1 G3.3
KJLTIPLlED BY +1
DE'lEUIlNJINT OF
G2.1 G2.2
G3.1 G3.2
KJLTIRolED BY +1
192
MJLTIH.IED BY +1
DEmRMINANT OF
Gl.2 Gl.3
G3.2 G3.3
MJLTIH.IED BY +1
DE'lmMINANT OF
Gl.l Gl.3
G3.l G3.3
MULTIH.IED BY +1
,
DE'lmMINANT OF
. Gl.l G1..2
G3.l G3.2
MJLTIH.IED BY +1
DE'lmMIWNT OF
. Gl.2 Gl.3
G2.2 G2.3
MJLTIH.IED BY +1
DE'lmMINANT OF
Gl.l G1..3
G2.l G2:.J
KlLTIPLIED BY +1
DE'IERMINllNT OF
Gl.l G1..2
G2.l G2.2
193
IDLTIILIED BY 'mE lE'lERMlNlINT OF
HI.1 m.2
H2.1 H2.2
IDLTIR.IED BY -1
DE'lERMINANT OF
. G3.3
MlJLTIR.IED BY -1
DE'lERMINANT OF
. G3.2
MlJLTIR.IED BY -1
DE'lmMlNlINT OF
G3.1
MlJLTIR.IED BY -1
DE'lmMlNl\NT OF
G2.3
KlLTIR.IED BY 'mE lE'mRMlNlINT OF
H3.2
IDLTIR.IED BY -1
DE'lERMINANT OF
G2.2
IDLTIR.IED BY 'mE IETERMllillNT OF
H2.2
IDLTIR.IED BY -1
DE'lmMINlINT OF
G2.1
IDLTIILIED BY 'mE lE'mRMlNlINT OF
HI.2
KlLTIR.IED BY -1
DE'lmMlNlINT OF
Gl.3
194
G1..3
KlLTIH.IED BY mE lE~llNT OF
83.1 -
MJLTIH.IED BY -1
DE'nmIINllNT OF
. Gl.2
KlLTIR.IED BY mE·IE~llNT OF
H2.1
KlLTIH.IED BY -1 .
DE'lfRMINllNT OF
Gl.1
KlLTIR.IED BY mE IE~T OF
Hl.1
195
AS only h12 is non zero, the numerator term SN is zero.
table (6.1) •
i
+ det'l
r g22
g32 ['" "" ] ['n 'n 1["" ~31
g23
g33
1
J h32 h33
+ det.
g31 g33 _ h32 h33
+ det. [ 'n '" 1['H ,,,] ['" ,,,] [On + det. "n ]
g31 g32 h22 h23 g32 g33 h31 h33
+ det.
[,,, ,,,] [On "" j ['H '" ]["H 'n ] + det.
g22 g23 h31 h32 . g21 g23 h31 h32
+ det • .['Hg21
'n ]
g22
['H
h21
"n ]
h22
-g12 h 21 - gllh ll =
SO = -g21,h12 (6.4)
196
Therefore the total transfer function element t2l is
3
893.9 10
=
S + 0.04644
=
893.9 10 3
(6.5)
893.9 103 S + 0.04644 - 893.9 10 h12
3
1 - S + 0.04644,h
12
2000
at steady state t2l = (NE) = Le.
MF max 0.00035
6
h12 = -0.123 10- (kg fuel/ ) (6.6)
r.p.m
1
-24.963
3
-7.1779 (S+0.04491) 1
(T) = (S + 0.1564) 893.9 10 256.04 103 (S+O.04482) (6.8)
[
-12.711 -3.655 (S+0.04482)
It can be compared with equation (6.1) >bich shows that the single pole
197
PROGRAM OUTPUT TABLE (6.2) TO CALCULATE THE NUMERATOR OF t11
!fJLTIRoIED BY +1
DE'IlruIlNlINT OF
Gl.1 Gl.2 Gl.3
G2.1 G2.2 G2.3
G3 .1 G3 .2 G3 .3
WLTIRoIED BY -1
DE'IlruIlNlINT OF
. Gl.1Gl.3
G3.1 G3.3
MULTI~D BY -1
DETERMINANT OF
G1.1 Gl.2
G3.1 G3.2
;. ~. MJLTIPLIED BY '!HE;IETERMlNlINT Of
H2.3
WLTIRoIED BY -1
DE'JERMllillNT OF
Gl.1 Gl.3
G2.1 G2.3
MJLTIRoIED BY -1
DE'IlruIlNlINT OF
. Gl.1 Gl.2
G2.1 G2.2
MJLTIRoIED BY '!HE IETERMlNlINT OF
H2.2
198
PROGRAM OUTPUT TABLE (6.3) TO CALCULATE THE NUMERATOR OF t
12
KlLTJ:R.IED BY -1
DE'IERMINlINT OF
Gl.l Gl.2 Gl.3
G2.1 G2.2 G2.3
Gl.l Gl.2 Gl.3
KJLTIR.IED BY -1
DE'IERMINlINT OF
Gl.2 Gl.3
Gl.2 Gl.3
MULTIIt.IED BY +1
DE'IERMINlINT OF
Gl.l Gl.2
Gl.l Gl.2
MULTIIt.IED BY -1
DE'IERMINlINT OF
Gl.2 Gl.3
G2.2 G2.3
KlLTIR.IED BY +l
DE'lmMINlINT OF
Gl.l Gl.2
G2.1 G2.2
199
PROGRAM OUTPUT TABLE (6.4) TO CALCULATE '!liE NUMERATOR OF t
22
WLTIPLIED BY +l
DE'lERMINllNT OF
Gl.1 Gl.2 Gl.3
G2.1 G2.2 G2.3
'G3.1 G3.2 G3.3
MULTIPLIED BY WE DETERMINllNT OF
HI.1 m.3
1D.1 1D.3
MJLTIPLIED BY -1
DE'lERMINllNT OF
"<;2.2 G2.3
G3.2 G3.3
MJLTIPLIED BY WE DETERMINllNT OF
ID.3
MULT:[pIjED BY +l
DE'lERMINllNT OF
G2.1 G2.2
G3.1 G3.2
MJLTIPLIED BY WE DETERMINllNT OF
HI.3
MULTIPLIED BY +l
DE'reRMINllNT OF
G1.2 Gl.3
G2.2 G2.3
MULTIPLIED BY WE DETERMlNllNT OF
ID.I
MlLTIPLIED BY -1
DE'reRMINllNT OF
Gl.1 Gl.2
G2.1 G2.2
MlLTIPLIED BY WE DETERMlNllNT OF
HI.I
200
.... PROGRAM OUTPUT TABLE (6.5) TO CALCULATE THE NUMERATOR OF t31
MULTIPLIED BY +1
DE~lINTOF
Gl.l Gl.2 Gl.3
G2.1 G2.2 G2.3
G3.1 G3.2 G3.3
I-IJLTIPLIED BY mE IE'lERMINlINT OF
H2.1 1l2.2
113.1 113.2
MJLTIPLIED BY +l
DE'lfIDIlNlINT OF
G2.1 G2.3
G3.1 G3.3
MULTIPLIED BY mE DETERMINlINT OF
113.2
MJLTIPLIED BY +l
DETERMINllNT OF
G2.1 G2.2
G3.1 G3.2
MJLTIPLIED BY +l
DE~llNTOF
Gl.l Gl.3
G3.1 G3.3
MULTIR.IED BY mE DETERMINllNT OF
113.1
MlLTIPLIED BY +l
DE~lINTOF
Gl.l Gl.2
G3.1 G3.2
MJLTIPLIED BY mE DETERMINlINT OF
H2.1
201
PROGRAM OUTPUT TABLE (6.6) TO CALCULATE THE NUMERATOR OF t32
MJLTIPLIED BY -1
DE'lERMINANT OF
. G1.1 Gl.2 Gl.3
G2.1 G2.2 G2.3
G3.1 G3.2 G3.3
MJLTIPLIED BY +l
DE'lERMINllNT OF
G2.2 G2.3
G3.2 G3.3
MJLTIPLIEIJ BY -1
DE'lERMINllNT OF
G2.1 G2.2
G3.1 G3.2
MJLTIPLIED BY +l
DE'lERMINllNT OF
Gl.2 Gl.3
G3.2 G3.3
MJLTIPLIED BY -1
DE'lERMINANT OF
. Gl..1 Gl..2
G3.1 G3.2
MJLTIPLIED BY THE IETERMlNllNT OF
Hl.I
202
the feedback loop, the calculated value is zero, because of the dependency
between the output torque and speed at the given constant value of external
Fig. (6.3).
PI
I NE
NED
MP'
Tq
Fig. (6.3)
pure oxygen. Therefore the action of the link fl must settle down to zero
at steady state after a step in the demanded speed NED, so the simplest
transfer function of fl is
fl = kl·S (k l is gain).
fl = S (6.10)
203
In the simplest arrangement the independent fuel input is controlled
k2 can be obtained.
NE
From the matrix (6.8) and fig. (6.3) the term NED' is given by
3
893.9 10 k2
NE
= = 1 (at 5=0)
NED (S + 0.1564)
0.1564 6
= = 0.17496 10- (kg fuel/RPM) (6.11)
3
893.9 10
fuel, the effect of demanded load (or speed at constant external load)
needs to be considered.
As the load increases, the optimum injection angle moves far from
T.D.C., e.g. at 40% demanded load the optimum injection angle is -22 deg.
at 2000 r.p.m. (See Chapter 4, computer output table (4.8e) , while the
optimum injection is -28 deg. when the load is 100% (See output table (4.8a)
-28 - (-22)
= - 0.1
100 -40
The controlling link f3 must achieve this shift due to the variation
of load, i.e.
6
0.17496 10-
%LOAD = *loo*NED
MFMAX
204
and MI"MAX = 0.00035 (DATA) (kg./R,P,M. )
Substitute in (6.12)
NED , ~~L1J
I f3
Fig. (6.4)
(F) = [ s 6
0.17.5 10- -0.005] (6.14)
for obtaining more unknowns like the condition of optimum injection advance
* Use of the compensation controller (F) in~;ide the loop is not possible,
as the output speed and torque of the engine are dependent on each other
• For a single pole transfer function matrix with common pole in all
gij' the pole is shifted in all terms due to the effect of any feedback
link h • No extra poles are generated when using one or more feedback
ji
links h" •
J~
• When obtaining a feedback link for constant steady state gain the pole
• In the matrix No. (6.2), only one feedback element is non zero;
205
in fact the other h might take any constant real value, not necessarily
ji
zero, so in this case there are many possible values for the CH) matrix,
control of the pressure can be achieved through the use of air or pure
oxygen injected into the engine, this is not discussed here because it is
206
CHAPTER 7
Specifically:
* The residual gases are about 3% of the max. mass of gas per
cycle; the main observable effect is to raise the temp. of the fresh
air charge.
* The max. temp. usually occurs (2 0 to 50) degrees after the max.
pressure.
* The exhaust pressure and temp. reach min. value at the optimum
injection advance.
* The max. value of specific heat appears just after the T.D.C. or
of the pressure curve (graph 2.5) at the end of Chapter 2, and the
position of optimum IMEP (graph 2.1) for each engine speed; the
207
experimental results show similar behaviour (fig. (3.2.a) and fig.
very small, (provided that the change is around room temp.) and thus
can be neglected.
and temp., although the max. pressure and temp. occur in the same
effect of the heat release curve upon the dynamic performance of the
engine.
208
* In Chapter 5 an easy method of obtaining the forward compensation
matrix (F) and feedback matrix (H) is discussed. The method depends
on obtaining the total transfer function matrix (T), then getting the
knowing their basic transfer function matrix (G) without using matrix
inversion or division.
* The reader need not know much about the theory or the procedures
* It is found that the use of the (F) matrix inside the feedback
loop is impossible, as the output torque and speed at the engine are
209
1, Transient response modelling using the digital computer model
needed).
aspects:
210
3. computer aided feedback design on the basis of pole-zero
paragraph 5.5.
211
R E FER E N C E S
1. Anonymous.
Latest EMD engines aim at cutting fuel consumption.
Railway gazette inter. Jan.1980, Vol.136, No.l, pp. 58-59.
2. Anonymous.
Diesel electronic fuel injection investigated.
Automot.Eng. (USA) Vol.89, No.2, pp.62-67 (Feb. 1981)
3. Anonymous.
Factors affecting the fuel consumption of road vehicles.
Report - British Rail Library - Derby, MRTL 14100, 621436-01
R.C.C. Publication No.171961
London British Transport Commission.
4. Aizerman, M.A.
Theory of Automatic Control.
Pergamon 1963.
9. Barnett, s.
Introduction to mathematical control theory.
Clarendon Press, Oxford, 1975.
10. Barnett, S.
Matrices in control theory.
Van Nostrand Reinhold Company, London, 1971.
212
12. Belletrutti, J.J. and MacFarlane, A.G.J.
Characteristic Loci techniques in multi-variable-control-
system design.
Proc. lEE Vol.118 No.9, Sept. 1971, pp. 1291-1297.
18. Bowns, D.E., Cave, P.R., Hargreaves M.R.D. and Wallace F.J.
(Bath Univ.)
Transient characteristics of turbo-charged diesel engine.
Conference on engine performance modelli~g 22-23, May 1973,
London I.M.E. pp.13-23.
213
24. El-Imam Y.A.
Simulation of a high-gain adaptive speed regulating system for
an internal combustion engine.
Simulation (USA) Vol.39, No.4, pp.133-139 Oct. (19B2).
27. Emre, E.
Pole-zero cancellations in dynamic feedback system.
Proc. of the 21st IEEE conference on decision and control
Orlando FL, USA B-10 Dec.19B2, pp. 422-423, Vol.l.
2B. Fallside, F.
Control system design by pole-zero assignment.
Academic Press 1977.
29. Furukawa, H.
Digital simulation of turbocharged diesel engine with air
injection into the compressor.
M.S.C. Thesis, 1972 UMIST
30. Glassman, I.
Combustion.
Academic Press 1977.
214
36. Ingham, D.R.
Heat release in a medium speed diesel engine.
M.Sc. Thesis UMIST, Jan.1969.
37. Judge, w.
High speed diesel engines.
Chapman & Hall Ltd., 1948.
39. Khovakh, M.
Motor vehicle engines.
MlR Publishers, Moscow, 1971.
215
48. MacFarlane, A.G.J.
Return-difference and return-ratio matrices and their use in
analysis, and design of multivariable feedback control systems.
Proc.IEE, Vol.117, ..No.lO,
.- - --- - _.
Oct. 1970, pp.2037-2049 •
51. Malone, F.
Diesel fuel : the squeeze is on.
Railway Age Aug.13th 1979. Vol. 180, No.15.
56. Meystel, A.
Decoupling and decentralized control of intelligent multiactuator
system.
Proc. 1982 IEEE International, Large scale systems symposium,
Virginia Beach, USA, 11-13, Oct.1982, pp.495-498, Vol.l.
57. Mikhail, S.
A filling and emptying model of a turbocharged diesel engine
for dynamic calculations.
M.Sc. 1975 UMIST
216
60. Powell, B.K. and Powers, W.F.
Linear quadratic control design for non linear I.C. Engine
systems.
ISATA 81 Proceedings, lOth anniversary international symposium
on automotive technology and automation, Stockholm, Sweden,
7-11, Sept. 1981.
(Croydon, Surrey, England: Automotive Automation 1981) P.29/1-18
Vol.l.
64. Richards, J.
An introduction to dynamics and control.
Longman 1979.
217
71. Seshu, S. and Reed M.B.
Linear graphs and electrical network.
Addison-Wesley Publishing Company Inc. 1961.
B1. Whatton M.
Existing design C.I. Engines increase specific outputs.
Report: 13th CIMAC congress in Vienna
I.M.E. 1979 Railway Engineering International, pp.46-5l
July/Aug. 1979.
21B
83. Whitehouse N.D. and Stotter, A., Goudie G.O. and Prentice B.W.
Method of predicting some aspects of performance of a diesel
engine with a digital computer.
Proc. I.M.E. Vol.176, No.9 1962.
90. Winterbone D.E., Bendon E.S., Mortimer A.G. and Kenyon (UMIST)
and Stotter A. (TIlT) Israel
Transient response of turboeharged diesel engine.
S.A.E 770122 1978
91. Winterbone D.E., Benso~ R.S., Closs G.D. and Mortimer A.G.
A comparison between experimental and analytical transient
test results for a turboeharged diesel engine.
I:"M.E Proe. vol.190 22/76, 1976.
93. Wosehni, G.
A universally applicable equation for the instantaneous heat
transfer coefficient in the internal combustion engine.
S.A.E. trans. Vol.76 Paper 676931.
219
APPENDIX n.ll.
C mm ItmJ~ VI"UABLI::S
'lllE
C PI= FRON Q)HffiESSOR OlARECl'ARISTICS BAR
Pl=1.5 BAR
RACK=l.
HF11AX=O.0007 KG/CYCL
C ANGLLJ=-30 Dffi
TlIPUT=300. K
220
NSWY=4
00UTE=10
NSTP =180
E'lOOOC=-50
11XIJLT=24 DEX>
MXOILT=60. 0 DEX>
JIDVl.Xl!=1.
WRfWR=8
ETATUR=0.35 H.P.
221
MF=RACK*MFMl\X
KG/CYCL
C BooINING OF 'nIE CYCLE CALQJIATIONS
EETAV=0.90 -CC3.E-8)*CNE**2»
Ml=SPa'lT*VlMI\X*EETAV
f.1S=Sl'OIT*VSMAX*EETAV KG!CYCL
C DISSOCIATIOO OF INPUT GASES KG/CYCL
C INPUT MASS OF 02
MD2=C16/C16.+4.*14.»*CMS)
C INPUT MASS OF N2 KG
MN2=C4*141C16.+4.*14.»* CMS)
C INPUT MASS OF 002 KG
la)2=O.0
C INPUT' MASS OF H20
NH20=0.0
C INPUT ~lI\SS OF 02 WHIClI BURNED
ID2Bm=0.0
TI=TlIPUT·
PI=P1 *EETAV K
CVAIR=C1000/29.0)*C20.758+O.002248*CTA_273.» BAR
CVGAS=CVAIR J/KG.K
222
GAMAO=l.4
CSRmR='lURlWR
Pl=PA*(1+(736.*CSRlWR*E~CSR/(AIRMAS(I)*CVAIR*TAl»**
*(GAMAC/(GAMAC-l» BAR
C Pl=1.2
WRITE (1,265) PI
265 FORMATClHO,28H<DMPRESSOR aJTroT PRESSURE= ,F4.2,lX,4HBAR )
TO=TA*(PlIPA)**«GAMAC-l)/GAMAC) K
WRITE Cl,285) TC
285 FORH1\T(1X,25H<DlolPRESSOR aJTroT mlP.= ,F6.2,2HK )
CSRlWR= ( (AIRMAS (I) *CVAIR *TAl 1 (ETACSR*736.) ) *
«(PlIPA)**«GAMAC-l)/GAMAC)-l) H.P
WRITE Cl, 290) CSRlWR
290 FORl-1l1.T (lX, 1 8H<DlolPRESSOR IDlER= ,F6.3 ,4HH. P )
C 'lURBINE DYNAMIC mUATION IS
C «'lUlUWR-CSRIWRl /6 .28) * (6 Q/ N'lUR) =ITCH* (0 N'lUR/D T)
C AFTER COOLER IDIP.
ETA<DL=0.5
TlIEUT=TC-~<XlL*(TC-TA) K
WRITE Cl,295) TlIEUT
295 FORH.l\T(1X,24HAFTER CXlOLER AIR TEl-IP.= ,F6.2,3H K )
WRITE Cl,302)
302 FOIDlAT(1X,5IlS'lEP ,2X,6H1TOLUH ,1X,1lH'lRAPED GAS ,lX,9HPRESSURE ,
,lX,611'lEMP. ,lX,5IDWGAS ,lX,8HSP. HEAT,3X,5!1'lORK ,5X,5HFUEL )
WRITECl,305)
305 FORMAT(1X,6HDEGRE ,lX,4HCLT),5X,5HCKG) ,6X,5HCBAR),5X,3HCK),2X,
,7HCGlOOL) ,8HCJ/KG.K) ,lX,9H(BAR.LT) ,IX,9HCKGlSTP) )
WRITE Cl,310)
310 FORMATClX, , - - - - - - - - - - - - - - . - - - - - - - -
----, )
223
C CALalLATION OF INTAKE STROKE
WRITE Cl,315)
315 FORMI'.T ClOX, 'CALalLATION OF INTAKE STROKE' )
OEGRF?-O.
OEG
I'ORKI=O.O
VI=V2
LT
PRWC:l=V2
LT
AREAMJ=O.O
KG.OEG 2
DDEGRF?O.O
00 500 K= 1 , 181
OEGRF?OEGRE+ODEX>RE
OEG
OOEGRE=ClSO.O/F'I.CWrCNSTP » om
C OONDINTIONS OF NEXT STEP
PHI=C COmRE) *3.141180.0) RADIAN
VC:lCK) =V2+C CVlM/\X-V2) 12.0) *
.C4.3333-QOSCPHI)-3.3333*SQRTCl.-0.09*CSINCPHI)**2 ») LT
OVC:l=VCYCK)-PRVVCY
PFNVC'i=VCY CK)
IPl'MAS=OVC:l*SroIT*EETAV KG! STP
'l'RAMAS='lRlIMAS+lIBSCIPl'MAS) KG
DM02=IPl'MAS*16.ICI4.*4.+16.) KGlSTP
M:l2=M:l2+DM02 KG
DMN2=IPl'MAS*CI4.*4.)ICI4.*4.+16.) KGlSTP
NN2=r.tl2+1Nl2 KG
TCYCK)=CTlIPUT*IPTMAS*CVAIR+TI~*CVGAS)ICIPl'MAS*CVAIR+TRAMAS*
*CVGAS) K
~~ CM:l2*32. +~1N2*2S. +f.KXl2* 44. +MH20*IS.) 1 CM:l2+f.1N2+MOJ2+MH20) GlM:lL
PCYCK)=IPl'MAS*R*TCYCK)ICVCYCK)*l1'lAIR*IE-3 ) +
+'lRlIMAS*R*TCYCK)aVCYCK)*~*IE-3 ) BAR
TI=TCY(K)
PI=PCY (K) BAR
VI=VC:l(K) LT
CV02=Cl000I32.0) * (20.959+0 .00295* (TCY ([0 -273.) ) J/KG. [(
CVN2= ClOOOI 28.0) * (20.704+0. 00211 * ('I'CY (K) -273.) ) J/KG. K
cvm2=ClOOOI 44.0) * (27 .545+0 .00922* (TCY (K) -273.) ) J/KG.I(
CVH20=(1000/18.0) * (25.184+0 .00525* (TCY(K) -273.) ) J/I<G. K
CVGAS= (OO2*CV02~2*CVN2#1H20*CVH20HKXl2*CVQ)2) 1
1(M:l2#1N2+MH20+~2) J/KG.K
CVGAS= (C704.9+O.0708*TCY(K» *H1 +Cl833. 4+0 .4931 *TCY(K) ) *f.lAXHF) I
I (Ml+MAXMF) J/KG.K
CVGAS=CVOONS J/KG.K
IF (NVRITE.LT. GlRITE) GO 'ID 300
WRITECl,320) OEGRE,VCY(K),'lRlIMAS ,PCY(K) ,TCY(K) ,M'lGAS,CVGAS,
,WORKIP ,MFI(K)
320 FORMI'.TClX,F4.0,2X,FS.S,2X,FS.7,lX,F7.2,lX,FS.l,2X,FS.2,IX,F7.2,IX
"F7.2,2X,El2.6 )
NVRITE=O
300 COOTINUE
U'lRITE=NVRITE+1
500 CONTINUE
224
C CALClJLATION OF COMPRESSION STROKE
MAXMF=O.O
1'mKC=0.0
IDRKEX=O.O
0000lB=0.0
WRITE <1,323)
323 FORMAT (10X, 'CALalLATION OF OmffiESSION STROKE' )
WRITE <1,310)
00 700 K= 181 , 360
OFX>RE=OFX>RE+DOEX;RE
DDrx;RE= <180. O/FLOll.T CNSTP »
C INJECl'ION cnIDITIONS
IF «OBGHE-360.l .LT. PRI:inx::l GO 'ID 710
IF CDffiQ1B. GE. OlBL'lHl GO 'ID 71 0
225
C SP. HEAT AT Q)NSTJINT VCLUM J/OOL. K
C OUTFUT IWlS OF CO2 ill 'lHIS STEP IS (SEE CDMBUSTION mUATION )
C02I=(14*441198.0) *MFI (K) KG
oo:>2=MCD2-+<Xl2I KG
C MASS OF 02 WHIOi BURNED IN 'lHIS STEP IS
002BBN=(22.5*(2*16)/198.0)*MFI(K) KG
C THE RESEOOAL IWlS OF 02 IS
1{)2=/>02-~02Bm KG
C OUTFUT ~IASS OF H20 ill 'lHIS STEP IS (SEE CDMBUSTION mUATION )
H20I=(15*11V198.0) *MFI (K) KG
~1l2CP~H20!-H20I KG
OVGAS= (OO2*0V02#1N2*0VN2#H20*CVH2O!-HC02*0VCD2) /
/ (1{)2-Hfi20020!-MC02) J/KG. K
OVGAS=«704.9+O.0708~(K»*Ml+(1833.4+O.4931*TCY(K»*MAXMF)/
/(Ml+MAXMF) J/KG.K
OVGAS=CYCDNS
SUI-1SPC=SJMSl'C+0VGAS
AVACY=SUMSPClDEG01B J/KG. K
~(OO2*32.#1N2*28.+MOD2*44.+~1l20*18.)/(OO2#1N2+OOJ20020) G/OOL
ETAUSE=0.65* (EXES*MFMAXIMF) **0.1
226
I
L -________________________ I
DEXiT=O.O
DO 900 K= 361 , 540
c FOR EXPANSION CP*V) ** CN2l =<XNSTANT
DDEXiRE=ClSO.0/FLOATCNSTP II D&i
DEXiRE=DEXiRE+DD&iRE DEXiRE
227
CVGAS=CV~S
SUI1SOC=SUnSpc+cvGlIS J/KG.K
AVACV=SUl1SPClDOOO!B
J/KG.K
r-w:;AS= COO2*32. +MN2*28. HKX>2* 44. +MH20* 18.) / COO2+MN2+1-KX>2+MH20) G/MJL
C <XlBBUSTION TEMP. CAL<lJLATIONS
i
TI2=TCYCK)
PI2=PCYCK) K
TCYCK) =TI2+C CHFA'lRLCN<XlI13) *CALV*ETAUSE)/ C'lRAWIS*CVGAS» K
PCY CK) =PI2* C'lRAMAS+HFA'lRL CN<XlI13) ) *TCY CK) / C'lRAWIS*TI2) BAR
910 <XlNTINUE
228
WRITE (1 ,310)
DO 1100 K= 541 , 720
DDEX;RE=Cl80.0/FLOAT(NSTP » DEX;
DEX;RE=DEX;RE+DDEX;RE . DEX;RE
OLOOAS='lID\MAS
VI=VCY(K)
PI=PCY(K)
TI=TCY(K)
C RESOOAL ~IASS OF EArn GAS IS
RESIlJL='lID\MASI 0Lr.W!S 0/0
ro2~02*RESOOL KG
MN2=MN2*RESOOL KG
l>KXl2=l>KXl2*RESOOL KG
NH20=~iI20*RESOOL KG
229
C 350 FORl-lAT(1H ,'AVAREX>E SPC.HEAT AT Q)~IBUSTION=' ,P6.2,' J/KG.K')
C I'K>RK IXXm reR CYCLE
WORK:mRKEXiOORRC BAR. LT
IME~mRK/ CVUll\X-V2) BAR
WRITE <l,368) IMEP
368 FOmlAT(lOx,6HIME~ ,P6.3,41lBAR )
C CALL C10S1N
C CALL AXIPLOCO,lS0. ,200.,1,1,10,10,0. ,10.00,0. ,120., 'CYCLE VCLUl1E (L
C .T)' ,16, 'CYCLE PRESS. (BAR)' ,17)
C NCY=NSTPC+~X
C CALL GRAroLCVCY,PCY,NCY)
C CALL GRAroLCVCY,MFCUR,NCY)
C CALL DEVEND
2Ju
.,
,
I
231
AIN<DN=AIN<DN-DAINNC
123 OONTlNUE
PRE.VHP=BHP (I)
PRE.VNE=NE
RPH
PREITRK=RACK
12 OONTlNUE
AIN<DN=AIN<DN+2
10 <DNTlNUE
232
& llOOLE (OEXi) , ,20,' BMEP (BI\R)' ,10)
CALL PENSEL(2,0.4,3)
CALL GRAroLCllOOLLJ,BMEP1,NPl')
CALL MDV!02 (25. , 110. )
CALL OlAHOLC9HSPEED-1*.)
CALL PENSELC3,0.4,3)
CALL GRAroLCllOOLLJ, BMEP2 ,NPl')
CALL MDV!02 C25 ~ , 100. )
CALL OlAHOLC9HSPEED-2*.)
CALL PENSELC4,0.4,3)
CALL GRAroL CllOOLLJ , Bl-lEP3 , NPl')
CALL ~0IJ'102 C25 • , 90. )
CALL OIAHOLC9HSPEED-3*.)
CALL PENSELC5,0.4,3)
CALL GRAroL (llOOLLJ , Bl-lEP , NPl')
CALL MDV!02 (25 • , 80. )
CALL OlAHOL C9HSPEEO-4*.)
CALL Plca.E
CALL PENSELC1,0.4,3)
CALL AXIPLOCO,lSO.,100.,1,1,10,10,MAXINJ,MININJ,0.,260.,'INJECTION
& llOOLECDEXi) ',20, 'B.H.P.' ,6)
CALL PENSELC2,0.4,3)
CALL GRAroLCllOOLLJ ,BHPl ,NPl')
CALL OIAHOLC9HSPEED-1*.)
CALL PENSELC3,0.4,3)
CALL GRAroLCANGLLJ ,BHP2 ,NPl')
CALL OlAHOLC9HSPEED-2*.)
CALL PENSELC4,0.4,3)
CALL GRAroLCANGLLJ,BHP3 ,NPl')
CALL OlllHOL(9HSPEED-3*.)
CALL PENSELC5,0.4,3)
CALL GRAroL (ANGLLJ , Bffi!. , NPl')
CALL Oll\HOL (9HSPEEO-4*.)
CALL PlCCLE
CALL PENSEL(1,0.4,3)
CALL AXlPLOCO ,150.,100. ,1,1,10,10 , !\'\XI1U , t!TInm, O. , .17 0, , I! U1':C'1'1(} 1
& l\NGLECDEG) ',20,'fUu.. CUJ.O(G/IHP.lut) ',20)
CALL PENSEL(2,O.4,3)
CALL GRAroL(ANGLtJ,FUELl,NPr)
CALL 11J\I'102 (25. , llO. )
CALL OlllHOL(9HSPEEO-1*.)
CALL PENSELO,0.4,3)
CALL GRAroL(ANGLLJ,FUEL2,NPl')
CALL ~0IJ'102 (25. , 100. )
CALL CHlIHOL(9HSPEEO-2*.)
CALL PENSELC4,0.4,3)
CALL GRAroL (ANGLLJ , EUEL3 ,NPl')
CALL ~0IJ'102 C25. , 90 • )
CALL OIAHOL C9HSPEED-3*. )
CALL PENSELC5,0.4,3)
CALL GRAroL (ANGLTJ , EUEL , NPl')
CALL MDV!02 C25 • , 80 • )
CALL OIAHOL C9HSPEED-4*.)
CALL Plca.E
CALL PENSEL(1,0.4,3)
233
CALL PENSEL(1,0.4,3)
CALL AXIPLO(0,150.,100.,1,1,10,10,MAXINJ,MININJ,0.,100.,' INJECTION
& AOOLE(Om) ',20, ''IDRQUE (KGF.M) , ,14)
CALL PENSEL(2,0.4,3)
CALL GRAroL (AOOLLJ , 'IDaJl ,NPl')
CALL ~0\1'ID2 ( 25 • , 55. )
CALL OIMOL(9HSPEED-l*.)
CALL PENSEL(3,0.4,3)
CALL GRAroL(AOOLLJ, 'IDRQ2 ,NPl')
CALL ~0\1'ID2 (25. , 45. )
CALL QllIHOL(9HSPEEO-2*.)
CALL PENSEL(4,0.4,3)
CALL GRAroL (AOOLLJ, 'IDRQ3 ,NPl')
CALL ID.rro2 (25. , 35. )
CALL OIMOL(9HSPEED-3*.)
CALL PENSEL(5,0.4,3)
CALL GRAroL(AOOLLJ,'IDRQ ,NPl')
CALL MJ<J'102 (25. , 25. )
CALL OIMOL(9HSPEED-4*.)
CALL PICCLE
CALL PENSEL(1,0.4,3)
CALL AXIPLO(0,150.,100 .,1 ,1,10 ,10,MAXINJ ,11ININJ, O. ,160.,' INJECTION
& AOOLE (om) , ,20,' MI\X PRESS (BAR)' ,15)
CALL PENSEL(2,0.4,3)
CALL GRAroL(AOOLLJ,P31 ,NPl')
CALL MJ<J'102 (25. , 55. )
CALL OIMOL(9HSPEED-l*.)
CALL PENSEL(3,0.4,3)
CALL GRAroL(AOOLLJ ,P32 ,NPl')
CALL ~0\1'ID2 (25. , 45. )
CALL OIMOL(9HSPEED-2*.)
CALL PENSEL(4,0.4,3)
CALL GRAroL (ANGLLJ, P33 ,NPl')
CALL ~0\1'ID2 (25. , 35. )
CALL QllIHOL (9HSPEED-3 *.)
CALL PENSEL(5,0.4,3)
Cl\LL GRAroL (I\NGLLJ, 1'3 , NPl')
CALL 10\1'ID2 (25. , 25. )
CALL QlAHOL (9HSPEED-4*.)
CALL DEVEND
CALL EXIT
rno
arro~1
234
l\PPENDI X 1. 2
THE SIMPLIFIED DIGITAL COMPUTER MODEL OF DIESEL ENGINE CONTAINING TRANSIENTS
RESPONSE OF TURBOCHARGER. (THE INTAKE AND EXHAUST STROKE IS NOT CONSIDERED).
236
TIM=O.O
C BIDlliIN3 OF mE CYCLE CALaJLATIONS
00 10I=1;NPr
NVRITE=<WRl'l'E
E~ =0.90-«4.2E-S) * «NE -500.)**2.»
V1=V1M1\X*~ LIT
VS=VSMllX*~ LIT/CYCL
MI. =SP<WT*Vl. KG/CYCL
M3=SPCWT*VS KG/CYCL
C DISSOCIAT1".OO OF INIUT GASES
C INIUT MllSS OF 02
MD2=(16/(16.+4.*14.»*MS KG
C INruT MllSS OF N2
MN2=(4*14/(16.+4.*14.»*MS KG
C INIUT MllSS OF . 0)2
M:D2=O.O
C INIUT MllSS OF H2O
Ml20=0.0
C . INIUT MllSS OF 02 WHIQI BURNED
ID2BR-f=O.O
WRlTE(1,210) NE , RACK
210 FORMAT(Im,14HENGINE SFEED= ,F7.2,2X,6HR.P.M ,10HllND RACK= ,F4.2)
C ArJl<XN= 5.
l\N3LLJ(I)=ArJlCDN-(20.+10. * (MF/MFMAX) )
WRlTE(1,225) ANGLLJ(I)
225 FORMlI.T(lIIJ ,20Hm:;LE OF INJEm'ION= ,F6.2 ,2X,4HIE3 )
E~1.25+(1.0E-5)*NE
EXFCtl'l=1.26
IF (NE. GT. NEMm) GO 'ID 303
~=O.O
303 OJNTrnUE
AIruUlF I£/MF
WRITE(1,235) AIRFUL
235 FORMlI.T(lHO,21HREAL AIR/FUEL RATIO= ,F6.3 )
AIRMI\S (I) =1£*NE*2./ ( srnOKS*60 .0) KG/SEC
FUELI1:FMF*NE*2.0/ (SUl.OKS*60.) KG/SEC
GASWIS=AIRMAS (I) +FU ELI£ n:;/SEC
WRITE(1,250) Ml ,Cu\SMI\S
250 FORMlI.T(1IIl,26HM1\SS OF Am FER cyo:.E Ml= ,F10.8,3!OO ,1OX
• ,SmASWIS= ,F7.5,7!OO/SEC )
C O)MffiESSffi RWER
C EETA <DMBUSSOR MJST BE TAKEN FRCl1 ITS OIARACI'ERISTICS
ETJlCSR=O .5
TC=TA* (Pl/PA) ** ( (GAMI\-l) /GAMI\) K
WRITE (1,306) TC
306 FORMlI.T(lOX,25H<DMIruSSOR ClJTRJT TEMP.= ,F6.2,2HK )
C mE SFECIFIC HEAT IEPENOO CN TENP.
~AIR=(1000/29.0)*(20.75S+O.002248*(TC-273.» J/KG.K
<.vGAS=~Am J/KG.K
CSRlWR =«AmMI\S(I)*<.vGAS *T1)/(ETACSR*736.»*
« (Pl/PA) ** «GAMI\-l) /GAMI\) -1) H.P
WRITE (1,3 07) CSRlWR
307 FORMlI.T(lOX,lSH<DMIruSSOR RWER= ,F6.3,4HH.P)
WRITE(1,30S)
237
WRITE (1,308)
308 R)RM!IT(lHO,12HSTEP(DEl3RE) ,lX,16I!1lCL.OF CYCr..(LT) ,lX,17HmFS.OF C'l
•a.. (BAR) ,lX,1611TEMP OF CYCr..(K) ,lX,6HmGAS ,1X,llHOl(J/KG.K) ,IX,
,13lWOOK(BAR.LT) ,21HEUEL BtJmED/Sl'EP(KG) )
238
C OTH2o. (1000/18.0)* (25.184-1{) .00525* (TCY(K) -2:13.) ) J/KG.K
C OTGl>S= (ID2*OlOZ+m2*<.VN2+MI2O*CVH20HlD2*0l0)2) /
C / (ID2+m2+MI20HlD2) J/KG.K
C OOl£Usrr.oo CALO.JLATIOOS
IF( (DEXiREl-ERIDB) .LT.360.0) GO 'ID 910
IF(DEXi<M3.GT. (OBL'lH-IDEXiRE» GO 'ID 910
N'/'RITE=(}lRITE
DEXiCM3=I:'00<M3+IDEX:iRE 003
NODM3= IFIX (DEXi<M3+0 .5) APROX
moo (NODl£) =DEXiRE
XOOM3=NODl£*TDSTP
HETFN2=(XO)l£**0.25)*EXP(-D.13*ABS(XODMB-6.5»
HEA±RL(NODl£)=«2./CMBL~)-(2./CMBL~**2.)*(DEGCMB-0.5»*MF le/SW
H~(NODl£)=(HETFN2+HETFNl)*MF/(2.0*~) IC/SI'EP
HE:rFm =HET1iN2
MFI(K)=~(NODl£) IC/SI'EP
AREAM:)=l\REAIDt~(NODl£) *DEXiCMB le. 003
MIOOIF=M\XMF+-MFI (K) I<G
<N'rnOD= (AREAM)/MAXMF) +mID13 om
MfCUR(NODl£)=MFI(K)*l.E6/DDEXiRE M>/DEXi
MWGAS=(ID2*32.+m2*28.+M002*44.+MI20*18.)/(ID2+m2+M002+MI20) QlIDL
esp. HEAT NI: aNsrJNT Va.UM J/IDL. K
C OO'mJT MISS OF 0)2 m ~IS STEP IS (SEE OOM3UsrJON B;lUATJON )
C02I=(14*4~198.0)*MFI(K) I<G
Ml)2=M:D2+cn2 I le
C MASS OF 02 WHICH BURNED ill ~IS STEP IS
M02BBN= (22.5* (2*16)/198.0) *MFI (K) I<G
C 'lHE RESEIlJAL MASS OF OZ IS
M02=/02-/02Bm le
C OO'mJT MISS OF H2O ill ~IS STEP IS (SEE OOM3UsrION B;lUATION )
H20l= (15*18/198.0) *MFI (K) ~ I<G
MI20=MI201- H2O! le
0lGl>S=«704.9-1{).0708*TCY(K»*Ml+(1833.4-1{).4931*TCY(K»*MAXMF)/
/ (Ml +MI'\XMF) J/KG.K
C 0lGl>S=838. J/KG.K
C ml>BUSl'ION 'lUll'. CALQJIJl'l'lONS
ET!'IJSE=O .85
'IRAMI\S= Ml. +MI'\XMF KG
TIO='lcr (K) K
PIO=OCY(K)
TCY (K) =TIO+ ( (HEA'IRL (NCOl£) *CALV*ET!'IJSE) / (TRAMIIS *OlGAS) ) K
C mIS B;lUATIOO IS P1*Vl/Tl=I'2*V2/T2
PCY(K)=PIO*(TRAMl\StH~(NODl£) ) *TCY(K)/('lRAMI\S*TIO) BAR
C PCY(K)=PIO+(0.0036)*(2500.-'lcr(K» BAR
C PCY(K)=PI*«VI/VCY(K»**2.8) . FAILED BAR
C PCY(K)=PIO+( (1.085ES)*H~(NODM3» BAR
C PCY(K)=PIO+O.0025 *(TCY(K)-554.25) BAR
910 CDNTlNUE
23')
905 CDNTlNUE
C OOliJST VMNE OIENIN3
IF(OmRE. LT. (540 .+E.VCBTC)) GO 'ID 903
oOm'J.!:lDmRE*60 .0/ (NE*360 .0) UNIT T
F£X=OCY(K)
C l'CY (K) =F£X* (EXP (-100. *DDmT)) BAR
l'CY(K)=F£X*(O.97S+(NEVIOOOOO.0)) BAR
903 CDNTlNUE
240
MFI(K)~.O
VCY(K)~.O
PCY(K)~.O
MroJR(K)~.O
1100 OONTINUE
C 'lUffiINE IUiffi
GNW:X=1.30
C ETA'lUR= FRCM 'lUffiINE OlARACl'ERI5'rICS (tUl' CDNSl'lNT)
ETA'lUR=O .4
'lURlWR=E1'A'lUR*GASMIIS*ClJGAS *'IS (I) *
«l-(PA/PS(I) )**( (Gruw:X-l)/Gruw:X» / (736.0) H.P
WRITE (1,425) 1URH'iR
425 FORMAT(lOX,15H'lURBINE RWER: ,FS.3,3HHP)
TMroU'l'='15 (1) * ( (PA/PS (1) ) ** ( (GAM1IEX-l) /GAM1IEX) ) K
WRITE (1,435) TMFOUT
435 FORMAT(10X,22H'lURBINE CXJTRJT TDIP.= ,F6.l,2X,2HK )
C COCLING CALaJIATIONS
ALLHE'l'=FUELI£*CALV-BHP(I) *736-{;ASMIIS*ClJGAS * (TMFOUT-TA) J/SEC
TlW=TAl-lO K
C RCC , ECC IS OXLING <DEFFICIENTS OF RllDIA'IOR JlND EmINE BIDeK ro.J/I<GK
RCC=0.8*NE
ECO=300
T2W( I) ='llWf-ALLHET/ (RCCI- OCC) K
WRITE (1,450) T2W(I)
450 FORMAT(lOX,13IWATER 'lEMP.= ,F6.2)
IF(T2W(I) .GT.350) GO 10 1000
C CLQ3E 'mERMJSTATE M (l)LD CDNDITIONS
T2W(I) =353 K
1000 OJNTINUE
241
C AFl'ER a:xx.ER TEl-IP.
E."l'P£OIFO .5
TlIlU'J.l='ro-El'JI(l)L* (TC-TlW) K
WRITE(I,465) TlIlUT
465 FORMM'(IOX,24HAFl'ER a:xx.ER AlR TEl-IP.= ,F6.2,21lK )
WRITE (1,468)
468 FORMAT(IHO,5HTIME ,4X,4HRPM ,4X,71lB.H.P. ,2X,14I1'lUR;lUE(KGF.M) ,
.IX,20llEXT&lNlIL LOIID(KGF.M) ,2X,18H EUEL(KG/mp.HR) ,
,2X,10H'lUmO RIM )
TRPM(I) =IRITRPM
WRITE(I,472) TIM.RPM(I) ,BHP(I) ,~(I) ,EXLOAD(I) ,roa(I) ,'lRPM(I)
472 FORMM'(F5.2,FI0.2,lX,FS.2,IX,FI0.2,10X,FS.2,10X,FI0.5,10x,F7.1)
C lNTEGRlITIO OF ENGlNE R. P. M-
C UNrI'S OF IDMENr OF lNERTIA IS IG.M.SEC**2/RADIlN
lNER'lI.=2.209
lNER'l'E=O .5368
Z(I)=NE *(2.0*3.14*(INERTEtlNERTL»/60.0
N=1
'lOL=10. **(-4)
TEND=TIMH71'IM
IFAlD=O .
CALL D02BIIF(TIM, TEND, N, Z, 'lOL, DERlV, WSPC, !FAlL)
NE =Z(I)*60.0/(2.0*3.14*(INERTEl+-lNERlL»
WRITE (1,478) NE,TIM
478 roRMM'(IHO,4HNE= ,F6.1,4HRPM ,2X,6H'l'n£= ,F5.2,4HSEC )
IF (IFAIL.HJ.O) GO ID 15
WRITE(I,479) IFAlL
479 FORMM'(IHO, ' ERROR WlSSID IFAIIF ',I2)
15 <DNTlNUE
c wmOOIlIffiER Sffim
WRLOD(I)=3.25El-6 * TRPM(I) KG.M
'lURLD='lURLOD (I)
'1'Ul<rnQ(I)=736.*(WlUWR 1'l'RPN(I)) KG.l:\
WRlQ='lURrnQ (I)
COMPlU(I)=736. * (CSRlWR ITRPM(I)) KG.M
OIPIQ=CDMPlU ( I)
INERT'J.l=3.2El-4
ZT(l) =TRPM(I) * (2.0*3.14*INERTT)/60.0
N'J.l=1
TIM=TIM-DTIM
IFAlL'J.l=O
C CALL 002l'BF (TIM. ZT, GI', !TT, Nl', IFAlLT, H'lO, Iilll', rERIVT,
C , ZOT, ET, PI',QT, RRr, SSl')
CALL 002BAF(TIM. TEND, Nl', ZT, 'lOL, DERlVT, WSPCl', IFAlLT)
TRPM(I)=ZT(1)*60.0/(2.0*3.14*INERTT)
PR'IRPM='lRPM (I)
TIM='ImD
WRITE(I,550) TRPM(I) , TIM
550 FORMM'(IOX,13H'lUmINE RFlT= ,FI0.l,4HRPM ,6HTIMEl= ,F6.2,2X,4HSEC )
IF(IFAILT. HJ.O) GO 'ID 31
. WRITE (1,552) IFAIL
552 FORMM'(lHO, 'ERROR WlSSID IFAn. = ' , 12)
242
552 FQRMM' (llJJ, 'ElUlCR MASSOO lFAIL = ' ,12)
31 <mTINUE
243
CALL PICa.E
CALL PENSEL(1,0.4,3)
CALL AXIl'LO(0,150 .,100.,1,1,10,10,0. , TIf.RUN,'lURMlN, 'lURMAX, 'TIME(S)
&',7, ''lUmlNE Rm' ,11)
CALL PENSEL(2,0.4,3)
CALL GRAroL(TIME, 'lRPM, NPr)
CALL ClI1'BCL(900RAm-l*.)
CALL PICa.E
CALL PENSEL(1,0.4,3)
CALL AXIl'LO(0,150 .,100.,1,1,10,10,NEMlN,NEMAX,0 .,0.25,' EN3lNE RPM'
&,10,'AIR ELOi (I<G/SEC) ',17)
CALL PENSEL(2,0.4,3)
CALL GRAroL(RPM,AIRMIIS,NIT)
CALt OIl'BCL(900RAm-1*.)
CALL PENSEL(3,0.4,3) .
CALL GRAroL(RPM,AIruSl ,NPr)
CALL OIl'BOL(900RAm-2*.)
CALL PICa.E
CALL PENSEL(1,0.4,3)
CALL AXIl'LO(0,150 .,100.,1,1,10,10,NEMlN,NEMAX, .0,170.,' EN3lNE RPM'
&,10,' (mBusr~ mES (BIlR' ,20)
CALL PENSEL(2,0.4,3)
CALL GRAroL(RPM,PJ ,NPr)
CALL OIl'BCL(900RAm-1*.)
CALL PENSEL(3,0.4,3)
CALL GRAroL(RPM,PJ1 ,NIT)
CALL OIl'BCL( 900RAm-2*.)
CALL PICa.E
CALL PENSEL(1,0.4,3) .
CALL AXIl'LO(0,150 .,100 .,1,1,10,10,NEMlN,NEMAX, 7. ,12.0,' EN3lNE RPM'
&,10,' BMEP (BIlR)' ,10)
CALL PENSEL(2,0.4,3) .
CALL GRAroL(RPM,BMEP,NPr)
CALL CllNlOL (9a;RAm-1 * • )
CALL PENSEL(3,0.4,3)
CALL GRAroL(RPM,BMEP1 ,NPl')
Cl\.LL CllNlOL(9IGHAPH-2*.)
CALL PICCLE
CALL PENSEL(1,0.4,3)
CALL AXIl'LO(0,150 .,100.,1,1,10 ,10,NEMlN,NEMAX, 90. ,260.,' EN3lNE RP
&M' ,10,'B.H.P.' ,6)
CALL PENSEL(2,0.4,3)
CALL GRAroL(RPM,BHP,NPr)
CALL CllNlCL(9a;RAm-1*.)
CALL PENSEL(3,0.4,3)
CALL GRAroL(RPM,BHP1 ,NPr)
CALL ClI1'BCL(9a;RAm-2*.)
CALL PICCLE
CALL PENSEL(1,0.4,3)
CALL AXIl'LO(0,150 .,100.,1 ,1,10,10,NEMlN,NEMAX, 70.0,120.,' EN3lNE RP
&M' ,10,' 'lO~UE (KGF. M) , ,14)
CALL PENSEL(2,0.4,3)
CALL GRAroL(RPM,~,NPr)
CALL OIl'BCL( 91J3RAm-1*.)
CALL PENSEL(3,0.4,3)
CALL GRAroL(RPM,'lO~l ,NIT)
244
CALL GRAPOL( RPM, 'lD~1 ,NPr)
CALL OIl\IICL( 91(;RAm-2* • )
CALL PIca.E
CALL PENSEL(1,0.4,3)
CALL AXIH.O(0,150 .,100.,1,1,10,10,NEMlN,NEMAX, .08, .16,' ENGINE RPM'
&,10,'FUEL CON. (RGVIHP.HR) ',20)
CALL PENSEL(2,0.4,3)
CALL GRAPOL( RPM, FUEL, NPr)
CALL OIllHCL( 91(;RAm-1*.)
CALL PENSEL(3,0.4,3)
CALL GRAPOL( RPM, FUELl ,NPr)
CALL OIl\IICL( 91(;RAm-2*.)
CALL PICCLE
CALL PENSEL(1,0.4,3)
CALL AXIH.O(0,150.,100.,1,1,10,10,NEMlN,NEMAX,900.,1400.,'ENGINE R
&PM' ,10,' EXHUsr 'lEMP. (K) , ,15) .
CALL PENSEL(2,0.4,3)
CALL GRAPOL (RPM, '15 ,NPr)
CALL OIllHCL(91(;RAm-1*.)
CALL PENSEL(3,0.4,3)
CALL GRAroL(RPM,'151 ,NPr)
CALL OIl\IICL( 91(;RAm-2*.)
CALL PICCLE
CALL PENSEL(1,0.4,3)
CALL AXIH.O(0,150 .,100.,1,1,10,10,NEMlN,NEMAX,2., 5.0, 'ENGINE RPM'
&,10,'EXHUsr ERESS. (BAR) ',18)
CALL PENSEL(2,0.4,3)
CALL GRAPOL(RPM,PS ,NPr)
CALL OIl\IICL ( 91(;RAm-1* • )
CALL PENSEL(3,0.4,3)
CALL GRAOOL(RPM,PS1 ,NPl')
CALL OIllHCL (91(;RAm-2*.)
CALL PICCLE
CALL PENSEL(1,0.4,3)
CALL AXIH.O(0,150 .,100 .,1,1,10,10,NEMlN,NEMAX, 0 .,500.,' ENGINE RPM'
&,10,'WATER 'IEMP. (K)' ,14)
CALL PF1JSEL(2,O.4,3)
CALL GRAOOL (RPM. 'I2W , NPl')
CALL GlllHCL(9ffiRAm-1 *.)
CALL D1NEND
CALL EXIT
lliD
SUBRalTINE IERIV(ATIM, Z,A)
RFAL*8 A,Z,ATIM
DlMENSlON A(l),Z(l)
mMlm /TE/ ro~UE, EX'lLOD
A(l) =(ro~UEJ- EXTLOD) KGF.M
RE'lUm
END
SUBRaJTINE IERIVT(ATIMl', ZT, AT)
RFAL*8 AT,ZT,ATIMl'
DlMENSlON AT(l) , ZT(l)
mr-mv'lUIV 'lURLD, 'lURlQ, CMPlQ
AT(l) = ('lURIQ-'lURLD-<MPlQ) KGF.M
RE'lUm
END
245
APPENDIX (2)
realistic assumptions.
at various temperatures.
CnH2n+2 + 02 + CO 2 + H20
The number of carbon atoms in each molecule is (13 and 14), i.e.
the mean molecular weight of the DERV is (180 + 200) gram/mol,while the
mainly. C H •
8 18
It can be noticed that the fuel does not consist of only one hydro-carbon
chains.
heat under constant volume is function not only to the type of gas but to
the temp. of the gas also, so the relation has been linearized in a suitable
246
(TABLE 2.1)
...
COl HIO xlde CO
..'". .. ..
C
•
.:
~
..
too
t;;
"
~
~
.'" .
t;; ~
t;;
.'".
~
~
..'" '"" .
t;; ~
- '""
t;;
too ~
c
'" '" '".!l '" '""
'0
'""
...: "
"0
"0
.!l
.. .!!
. " " "0
-s "
"0
e ...
-6
.!!
0
...s
e
0
--
.!l
0
...
0
...s "0
- ....-...-
0
...
S
-..
"0
"0
...e
.!!
s
.!!
0
::::
..3
'4
S
>-
~
::::
u
39
37
•
35
•
•
•
33
31
~-=op= 0.005125
/
29 '
/
.
/
•
27
•
•
25
_. __.----_.......
0 200 400 800 1200 1600 2000 2400
Graph - a
248
Mean molar heat capacity of O at constant volume
2
CV02
(J/mol. deg)
28
27
26
•
25
Slop = 0.00288
24
23
/
.'
22
/
/.
/
21
20 >-
o 200 400 800 1200 1600 2000 2400
Graph - b
249
Mean molar heat capacity of N2 at constant volume
CIIN2
(J/mol. deg)
26
25
24
23 Slop 0.00198
22
21 /..
20
o 200 400 800 1200 1600 2000 2400
Graph - c
250
Mean molar heat capacity of CO at constant volume
2 •
CVC02
•
(J!mol. deg)
45
•
43
•
41 •
Slop ~ 0.00515
39
37
35 •
33
31
29
27
o 200 400 800 1200 1600 2000 2400
Graph - d
251
Mean molar heat capacity of AIR at constant volume
CVAIR
(J!mol. deg)
28
27
26
•
•
2$ •
•
•
24
23
Slop 0.6022
22
21
20
o 200 400 800 1200 1600 2000 2400
Graph - e
CVN2 = (1~0 )*120.704 + 0.00198 (T-273) 1 J/KG.K
- (5)
(HEATRL) is
(,88
MO~ ~ 198 IIE/\TH!.
M02 O. 2 * M1 - 688
198 HEATRL (7 )
253
The average specific heat of the gases is (use eq. (7), (8), (9), (10»
Emass*CV
CVGAS = Emass
= (0.2Ml-?~8HEATRL)*CV02+(0.8Ml)CVN2+(~EATRL)*CVH20+(~EATRL)CVC02
688 270 616
(0.2Ml- 198IIEATRL) + 0.8Ml + f98IIEATRL~EATRL
= 0.2Ml+CV02+0.8Ml*CVN2+ ( -72OCV02+27OCVH20+616CVC02)IIEATRL
198
Ml + HEATRL
I
--- call it equation ( 11 )
0.2Ml*CV02+0.8MlCVN2 = (0.2*CV02+0.8CVN2)Ml
[0.2(664+0.07125(T-273»+0.8(739.4+0.7071(T-273» [Ml =
-720CV02+27OCVH20+6l6CVC02
=
198
0.117 (T-273) [ =
l833.38+0.4g3l*T ( 13 )
For the uses of this equation in the digital computer program, the
amount of fuel (HEATRL) described by the heat release curve can be obtained
practically. i.e.
KMAX
HEATRL = f HEATRL(K)
K=O
254
In the Fortran prog:-am this will be as follows:
FUEL=FUEL+HEATRL(K)
use the final value of combustion temp., i.e. at the end of each step. But
be made.
Conclusion
combustion gases and to the temp. However it is clear that the dependance
upon the temp. can be 1inearised and average specific heat for gas mixture
of intake stroke.
especially in the case of petrol fuel, but the composition C H1a must be
a
used.
255
APPENDIX ( 3 )
mechanism.
idling speeds.
mechanism.
256
ADVCON= the new value of advance angle equal to the old
all jumps.
is functioning properly.
Advance -----..
0"1 diesel
eng~ne
digital
~.---- controller
Fig. ( 1
Block diagram of the digital feedback system for the diesel engine.
2')7
APPENDIX (4.1)
.tuLL.
C 'lHIS m£X;RAM CALOlLATE 'lHE 'IDTlL 'rnlNSFER FUocrION ELEMEm' TIJ
INTEGER L(lOO) , MAXL(lOO) , GRCPA(200,25) , GRCEB(200,25)
INTEGER HA(200,25) , HB(200,25) ,GA(200,25) , GB(200,25)
INTmER HA!, IlBI , GAI , GBI
WRrrE(l,lO) I , J
10 FORMII.T(lH , 'CALQJLATION OF 'lHE 'lOTlL mlNSFER ruOCl.'ION ELEMEm'
.T(',n,',' ,n,') , )
WRrrE (1,12)
12 FORMII.T(lH ,10X, 'WHERE T(I,J)=(G(I,J)+EN)/(l+SO) ')
WRrrE(l,lS)
15 FORMII.T(lH ,10X,' EN , SO ARE PS reFINED BELQol ')
WRITE (1,22)
WRrrE(1,22)
WRrrE(1,18)
18 EORMI\.T(lH " EN IS 'lHE ALGEBRAIC SUM OF 'lHE IRa:oCI'S OF ALL 'lHE
• EOLLCWOO MINORS, , )
WRrrE(1,20)
20 FORMII.T(lH ,'EAQI MJLTIILIEDBY+! OR-l PS INDICATED' )
WRrrE(1,22)
22 FORMII.T (lH " ')
25U
25 K=K+l
DO 50 Il=l, NNl.
50 IF(K. EO. L(Il) ) GO 'lO 25
HA(n:tA, HAI) =K
HB (IHB, IlBI) =K
GA(IGA, GAl) =K
GB(IGB,GBI)=K
KI=KI+1
HAI=HAI+l
IlBI=IlBI+1
GAl=GAl+l
GBI=GBI+1
IK=IK+1
HAI=l
IlBI=l
GAl=l
GBI=l
KI=l
K=O
IS'lOR=O
JSlOR=O
IF(NP.IQ.O) GO 'lO 90
GO 'lO 95
25 C)
GO 'ID 95
90 OJNTINUE
NIbl
95 OJNTINUE
1OOL= 1OJItl-l
100 OJNTINUE
L(1COL)=L(1COL)+1
IF(L(ICOL) .NE.l) GO 'ID 140
LMI\X=NLIWH
nIw.x=nIh-l
DO 110 00=1, nIMlIX
DO 120 1X2=l, nIMlIX
RClV=0
OJLffi=O
128 ISIGN=RCW+<DLmtNLMI'.X-l
RClV=0
OJLffi=O
WRl'l'E(1,22)
WRl'l'E(1,22)
IF(MOD(ISIGN,2).NE.0) GO 'ID 133
WRl'l'E (1,13 0)
130 FORMI\T(lH " IVLTIILIED BY +1 ')
GO 'ID 131
133 WRl'l'E(1,132)
132 FORMI\T(lH " IVLTIILIED BY -1 ')
131 OJNTINUE
120 OJNTINUE
110 OJNTINUE
NLw.x=NLMI'.X-l
IF(R.MI'.X. B;2.1) GO 'ID 200
IK=l
nIlFl
nIB=1
260
lllB=1
IGh=1
IGB=1
140 <XlNTINUE
IF(L(I<DL) .GT.MAXL(I<DL» GO ID 95
ISO <XlNTINUE
200 <XlNTINUE
WRITE(1,22)
WRITE(I,22)
WRITE (1,205)
205 FORMI\T(lH " SO IS 'lHE ALGEBRAIC SUM OF ALL 'lHE FtlLr.cwOO MlNORS
•mowCl'S' )
NLWIX=NN
K=O
NP=O
KI=1
IK=1
NNl.=NN-l
IQ)L=O
00 210 10=1 , NN
L(IO) =0
2J.O <XlNTINUE
225 K=K+l
IK=IK+l
KI=1
K=O
IF(NP.EQ.O) GO ID 290
GO ID 295
290 <XlNTINUE
NP=1
295 <XlNTINUE
261
295 <DNTINUE
1Q)I;= IQ)Itt1
300 <DNTINUE
L(IODL) =L(IQ)L) +1
IF(L(IQ)L) .NE.l) GO 'ID 340
IKMAX=IK-l
00 310 IXl=I, 1KMAX
00 320 IX2=I,IKMAX
WRITE(I,22)
WRITE(I,22)
IF (IDD(N:.MAX,2) .NE.O) GO 'ID 333
WRITE (1,330)
330 FORMI\T(lH " MJLTIEt.1FD BY +1 ')
GO 'ID 331
333 WRITE(I,332)
332 FORMI\T(lH ,'MlLTIEt.1ED BY -1 ')
331 <X>NTlNUE
CALL !l;'lMAT(IXl,IX2,NLMAX,NLMAX,GROPA,GROm,GROPA,GRaB)
320 <X>NTlNUE
310 <X>NTlNUE
NLf.1lIX=NLf.1lIX-l
1K=1 ~
340 <X>NTlliUE
350 CONTlNUE
400 SIDP
END
WRITE(I,510)
262
WRlTE(1,510)
510 FORMM' (lH ,5X, 'IE'l.ERMlNlflT OF ')
WRITE (1,528)
WRITE (1,529)
529 FORMM' (lH ,5X, 'MJLTIPLIED BY lliE IE'lliRMlNl\NT OF ')
528 FORMAT (lH " ')
00 540 IYl=l , LMAX
WRITE (1,535) (TITLE(2), (HB(IX2,IYl) ,HA(IXl,IY2» ,IY2=1,LMAX)
540 (l)NTINUE
26)
APPENDIX (4.2)
WRXTE(l,lO) I , J
10 FORMI\T(lH " oo.aJLlITICN OF 'llIE 'IOTN. 'lRlNSFER FUN::l'ION ELEMENT
.T(',n,',',ll,') , )
WRXTE(1,12)
12 FORMI\T(lH ,1OX, 'WHERE T(I,J)=(FG(I,J)+EN)/(l+SO) ')
WRXTE(1,15)
15 FORMI\T(lH ,10X,' EN , SO ARE PS IEFINFD BELQV' ')
WRXTE(1,22)
WRXTE(1,22)
WRITE (1,18)
18 FORMI\T(lH , 'o!:N IS 'lllE ALGEBRAIC SUM OF 'llIE PRWJcrs OF ALL 'lllE
• FOLL(Wm:; MINORS, , )
WRITE (1,20)
20 FORMI\T(lH " FAOI MJLTIlLIED BY +l OR -1 PS mDIOI.TFD ' )
WRITE (1,22)
22 FORMI\T(lH " ')
264
22 FORWl.T(lH " ')
25 K=K+l
00 50 Il=l, NNl.
50 IF(K. EO. L(Il» GO 'ID 25
HA(IHA, HAI) =K
HB(IHB,lBI)=K
GA(IGA,GAl) =K
GB (1GB, Gn) =K
IF(K.NE.I) GO 'ID 60
IS'IDR=l
HAI=HAI-l
60 ro:JTlNUE
IF(K.NE.J) GO 'ID 70
JS'IOR=l
HB I=HB I-I
70 ro:JTlNUE
KI=KI+l
HAI=HAI+l
HBI=IBI+l
GAl=GAl+l
GBI=GBI+l
IF(KI.LE.NLMAX) GO 'ID 25
IF(IS'IDR.NE.l) GO 'ID 85
IHA=IHMl
IGA=IGMl
85 al'ITlNUE
IK=IK+l
HAI=l
HBI=l
GAl=l
GBI=l
KI=l
K=O
IS'IDR=O
JS'IOR=O
GO 'ID 95
265
GO 'ID 95
90 <nmNUE
. NIGl
95 <mTINUE
IOOL= I<DIJI-l
100 <XX'ITINUE
L(IOOL)=L(IOOL)+1
IF(L(IOOL) .NE.l) GO ID 140
UI!\X=NLM1\X-l
mM1\X=IHl'r-l
00 110 !Xl=1, mMAX
00 120 IX2=1, mMAX
RCW=O.
())IJ.N=O
00 122 IR=1, NLMI\X
RCW=0Cl#1
IF(GA(IXI,m) .IQ.I) GO 'ID 124
122 cmTINUE
128 ISIGN=rooft<DLmtNLM1\X-l
RCW=O
<DLIN=O
WRITE(1,22)
WRITE (1,22)
IF(IDD(ISIGN,2) • NE. 0) GO 'ID 133
WRITE(1,130)
130 FORMII.T(1H , I f.1JLTIILIED BY +l ')
GO 'ID 131
133 WRITE(1,132)
132 FORMII.T(1H " MlLTIILIED BY -1 ')
131 <DNTINUE
120 <mTINUE
110 <XX'ITINUE
IF(IOOL.IQ.NN) GO ID 200
NLM1\X=NLM/\lH
IF(R.M1\X.IQ.l) GO ID 200
!K=1
mA=l
IHB=rl
IGA=l
2GG
lGlf=1
lGB=1
140 OJNTINUE·
IF(L(ICOL).GT.MAXL(ICOL» GO ~ 95
150 OJNTINUE
IF(ICOL.EQ.l) GO ~ 25
IOLD=ICOL
I<DL=I<DIt-l
L(ICOL) =L(IOLD) +1
GO,~ 150
-
200 OJNTINUE
WRlTE(I,22)
WRlTE(I,22)
WRlTE(I,205)
205. FORMAT(lH " SO IS mE ALGIDRAIC SUM OF ALL mE FOLILWOO MlNORS
.1ROWCl'S')
NLM!\X=NN
K=O
NP=O
KI=1
IK=1
NNI.=NN-l
I<DL=O
00 210 10=1 , NN
L(IO)=O
210 OJNTINUE
225 K=K+l
GROPA(1K, KI) =K
GROIB(1K, KI) =K
K1=KI+1
1F(K1. LE. NLM!\X) GO 'ID 225
IK=IK+l
KI=1
K=O
IF(NP.EQ.O) GO ~ 290
IF(L(1COL).NE.~(IOOL» GO ~ 300
GO ~ 295
290 OJNTINUE
NP-l
295 OJNTINUE
1OOL= room-I
300 OJNTINUE
267
300 CWTlNUE
L(ICOL)=L(ICOL)+1
IF(L(ICOL).NE.1) GO m 340
II<MAX=IK-1
DO 310 !Xl=1, II<MAX
DO 320· IX2=1, IKMAX
WRlTE(1,22)
WRITE (1,22)
IF(/oOD(lUlAX.2) .NE.O) GO 'ID 333
WRITE (1,330)
330 FORM1\T(1H·,' KJLTIRoIED BY +1 ')
GO m 331 -
333 WRITE(1,332)
332 FORM1\T(lH.,' KJLTIRoIED BY -1 ')
331 CWTlNUE
320 CWTlNUE
310 CWTlNUE
IF(ICOL.EQ.NN) GO m 400
. NLMI\X=Nr.MI\X-1
IK=l
340 <XlNTINUE
IF(L(ICOL).GT.MAXL(IooL» GO m 295
350 OONTlNUE
400 SJDP
END
SUBROJTINE OO'lliAT(IXl.,IX2,NLWIX,I&X,GA,GB,HA,IB)
WRITE (1,510)
510 FORM1\T(lH ,5X, 'lE'lERMINllNT OF ')
268
IX) 525 m=l, NLMAX
WRlTE(1,530) (TITLE(l) ,GA(IXl.,m) ,GB(IX2,IY2) ,IY2=l,tLMAX)
525 <mTINUE - .
WRlTE(1,528)
WRlTE(1,529)
529 FORW.T(lH ,5X, 'MJLTIH.IED BY WE IE'IERMlNllNT OF ')
528 FORMAT(lH " ')
DO 540 m=l , uwt
WRlTE(1,535) (TITLE(2), (1B(IX2,m) ,HA(IXl.,IY2)) , IY2=1,UWC)
540 <mTINUE
269
APPENDIX (5.)
I - A typical Example for using loops theory to demonstrate the pole and
function matrix.
1
2
s+l
(G)= _1_
(I-I)
[ 3
s+2
s + 4
where det. (G) = (s+l) (s+2) (I-2)
~ Xl ______~~____~~l
____~~______~ Y1
x
2
~~----~~----~--~~------~Y2
g22
Fig. (I. 1)
It is desired to move the poles to (5+3) & (5+4) by applying the forward
and (I-3)
27u
FIG. (1.2)
and then calculate the forward and feed back matrices (F) and (H) needed to
From the loop theory, it is observed that every element in the total
transfer function matrix has all the poles of the system (because they have
( T)=
2 (s+3)
(1-4)
3 (s+3) (s+4)
(s+3) (s+4) (s+3) (s+4)
271
when relocatin9 two poles in 9
11 and 9 21 , only two feedback elements
are needed h11 and h , Le.
12
h12
(H)=
[:11 (1-5)
0 J
Also we need at least two forward compensation elements, i.e.
fll
s+l 2f11 1
(FG) =
(1-6)
f22
S+2 3f22 J
fll f22
so = 0- - h
s+l 11
--h
s+2 12
272
COMPUTER PROGRAM OUTPUT TABLE (I-I)
MULTIPLmD BY· -1
DE'.lERMIw\NT OF
FGl.l FGl~2
FG2.l FG2.2
-
mLTIPLmD BY +1
D~·OF
-FGl.l FGl.2
FG2.l FG2.2
KlLTIPLmD BY -1
DETERMINANT OF
FG2.2
MULTIPLmD BY -1
DE'reRMINANT OF
FG2.l
MULTIPLIED BY -1
DE'reRMINANT OF
FGl.2
273
112.1
MULTIPLmO BY -1
DE'lERMIW\NT OF
FG1.1
274
=
f22
--h
s+2 12
fU (s+2)
= (1-9)
Also
f (s+1)
22
= (I-U)
1+S0
2
3f
22
S + (9f
22
- hll)S+ (6f
n - 4hU)
1 + SO
(1-12) .,."
7 (1-13)
12 (1-14)
= 12 (1-15)
= 3 (1-16)
(1-17)
275
6fU + h12 = 14 (1-18)
f22 = 2 (1-20)
f22 = 6
(1-21)
3f
22 = 3 (1-22)
9f - hU = 21
22 (1-23)
6f 22 - 4hU = 36 (1-24)
achieve the required poles and keeping a suitable value for the zeros in
From tll the equations (1-13) and (1-14) cannot be altered, but it
can be deduced from (I-IS) (1-16) and (1-17) that the value for fll
3 (1-25)
Also from equations (1-20), (1-21), (1-22) the range for the values of
(1-26)
Substitute in (1-13) and 1-14) and solve together to get the value
of h
ji
3 - 3h - 2h12 = 7
U
2 - 6h
U - 2h12 = 12
i.e. hU = 3" 4
and h12 = -1 (I-27)
27G
Substitute in (1-9) to obtain tll
3 (5+2)
= (1-28)
(5+3) (5+4)
Substitute in (1-12) to obtain t12
6(5+0.5967) (5+2.2367)
= (1-29)
(5+3) (5+4)
2(5+1)
= (1-30)
(5+3) (5+4)
From (1-28), (1-29), (1-30), (1-31) the total transfer function matrix
is
(T)= (1-12)
It is clear that six new poles originated due to forward and feedback
--=--':::'=::."::=-
compensation.
277
II. Use of forward compensation matrix to eliminate cross coupling effects.
~
(G)=
S!ll (II-l)
s+2 :1
Let us use the (F) matrix to obtain diagonal dominance in the (G) matrix
(II-2)
r fll
1 2
l,,, 0.5 3
=
and
Also
= 1.5 (rr-7)
fn
[,.,
=2 (rr-8)
f12
i:e. (F) =
-1
-:]
Then the forward compensated (G) matrix becomes:
0.5 (S+4)
0
(s+1) (s+2)
(FG) = (rr-9)
S
4
(s+1) (s+2)
279
III Example on feedback design.
(G)=
S!l S!3j
[--
2
s+l --
3
s+3
1
where det.(G) =(s+1) (s+3)
where
( H)=
[0
hll
:]
First step it is necessary to specify the matrix (T)
3 (s+3) s+3
(s+3) 2 (s+3) 2
(T)= (Ill.l)
§ (s+3) 3 (s+3)
(s+3) 2 (s+3) 2
Then the (T) matrix must be calculated according to the loop theory
1
tll =
(s+1) S + a S + a
l 2
3 (s+3) 3
tll = =
(s+3) 2 S+3
1 3
(l+a ) S + (l+a ) - s+3
l 2
280
COMPUTER PROGRAM OUTPUT TABLE (III.2)
MULTIPLIED BY -1
DE'lERMINl\NT OF
G1.1 Gl.2
G2.1 G2.2
mLTIPLIED BY TIlE DETERMINANT OF
H2.2
MULTIPLIED BY +1
DE'mRMINlINT OF
G1.1 Gl.2
G2.1 G2.2
MULTIPLIED BY TIlE DETERMINANT OF
H1.1 Hl.2
H2.1 H2.2
mLTIPLIED BY -1
DE'mRMINlINT OF
G2.2
mLTIPLIED BY '!HE IETERMINlINT OF
H2.2
MULTIPLIED BY -1
DE'mRMINlINT OF
G2.1
mLTIPLIED BY TIlE IETERMINANT OF
H1.2
MULTIPLIED BY -1
DE'mRMINlINT OF
G1.2
mLTIPLIED BY THE DETERMINANT OF
H2.1
281
H2.1
MULTIPLmD BY -1
DE'JERMINlINT OF
Gl.l
**** S'lOP
282
2
"3 and "2 = 0
2 (III-3)
then \ 1 - - S
3
check tll
1 3 (III-4)
=
1
= s+3 3 (s+1) (III-5)
= =
1 2
1 + - ( - - S) (s+3) 2
s+l 3
Similarly t21 is
2
s+l 6
= = -(s+3)
- - (III-6)
1 +
Finally t22 is
1
3
g22 - hll det. IGI ~ s( )
sf3 + 3 (s+1) (s+3)
=
1 + gllh 2
ll "3 s
1 ---
s+l
11 s + 9 (III-7)
2
(s+3)
Equations (II!. 4) ,(UI ,5), (III ,6) ,(In, 7) give the total transfer
function matrix.
3 s (s+l)
s+3 2
(s+3)
(T) =
6 11s + 9
s+3 2
(s+3)
2uJ