`
`
`
`
`
`Computer modelling of the automotive
`
`
`
`energy requirements for internal
`
`
`
`combustion engine and battery
`
`electric-powered vehicles
`J.R. Bumby, B.Sc., Ph.D., C.Eng., M.I.E.E., P.H. Clarke, B.Sc., Ph.D., C.Eng.,
`
`
`
`
`
`
`
`
`
`
`F.inst.E., and |. Forster, B.Sc.
`
`
`
`Indexing terms: Measurement and measuring, Instrumentation and measuring science, Computer simulation
`
`
`
`
`
`
`
`
`
`
`Abstract: In the paper the road vehicle simulation package Janus, developed in the Engineering Departmentat
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`Durham University, is described. Janus is a flexible simulation package that allows internal combustion engine
`
`
`
`
`
`
`
`
`
`
`
`
`
`vehicles, electric vehicles and hybrid vehicles to be simulated, and their performance and energy consumption
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`evaluated over standard driving cycles. The simulation techniques used in these programs are described and the
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`simulation program shownto produce results comparable with experimental data.
`©
`
`
`
`
`
`
`
`
`
`
`List of symbols
`
`
`a
`=
`vehicle acceleration, m/s?
`
`
`
`A
`= vehicle projected frontal area, m?
`
`
`
`
`
`BMEP
`= brake meaneffective pressure, kPa
`
`
`
`
`
`cc
`= engine cubic capacity, cm?
`
`
`
`
`Cy
`= aerodynamic drag coefficient
`
`
`
`Cp
`= discharged ampere-hours. A h
`
`
`
`
`CG
`= coefficient of rolling resistance
`
`
`
`
`Cr
`= engine compression ratio
`
`
`
`C,, C,, C; = motorloss constants
`
`
`
`
`= distance travelled, km
`d
`
`
`
`= EMF, V
`E
`
`
`= battery energy density, kJ/kg
`Ey
`
`
`
`
`= battery open-circuit
`terminal voltage when
`E,
`
`
`
`
`
`fully charged, V
`;
`
`
`= battery terminal voltage, V
`E,
`
`
`
`
`= cold engine fuel flow factor
`Sr
`
`
`
`
`
`= acceleration due to gravity, 9.81 m/s”
`g
`
`
`
`
`
`
`= armature current, A
`I,
`
`
`
`= battery current, A
`I,
`
`
`
`= field current
`I;
`
`
`= battery discharge current
`1,
`
`
`
`= wheelinertia, kg m?
`ly
`
`
`
`
`= DC motor EMFconstant
`ky
`
`
`
`
`= DC motor torque constant
`ky
`
`
`
`
`= polarisation resistance, Q
`K
`
`
`
`= rotational speed, rev/min
`n
`
`
`
`= Peukert index
`n
`
`
`= DC machine armature conduction loss, W
`Por
`
`
`
`
`
`
`= DC machine brush loss, W
`Pur
`
`
`
`
`
`= DC machinecore loss, W
`Poy
`
`
`
`
`
`= DC machinefield resistance loss, W
`Pr
`
`
`
`
`
`
`= input power, W
`P;
`
`
`
`= DC machine mechanical loss, W
`Pat
`
`
`
`
`
`= DC machinestray loss, W
`Poy
`
`
`
`
`
`= battery power density, W/kg
`Py:
`
`
`
`
`= gearbox losses, W
`Pear
`
`
`
`= battery ampere-hour capacity, Ah
`
`
`
`
`
`rp
`= wheel rolling radius, m
`
`
`
`
`R
`= ohmicresistance, 2
`
`
`
`Paper 4048A (SI), first received 1st November 1984 and in revised form Ist May
`1985
`
`
`
`
`
`
`
`
`
`
`
`
`
`The authors are with the Department of Engineering, University of Durham,
`
`
`
`
`
`
`
`
`
`
`Science Laboratorics, South Road, Durham DH] 3LE, United Kingdom
`
`
`
`
`
`
`
`
`1EE PROCEEDINGS,Vol. 132, Pt. A, No.5, SEPTEMBER 1985
`
`
`
`
`
`
`
`
`
`
`= battery state of charge
`
`
`
`
`= time
`
`= ambient temperature, °C
`
`
`
`= airgap torque, Nm
`
`
`
`
`= engine compression braking torque, N m
`
`
`
`
`
`
`= tractive effort at the road wheels, N
`
`
`
`
`
`
`
`= vehicle velocity, m/s
`
`
`
`= armature voltage, V
`
`
`
`= battery voltage, V
`
`
`
`= head(ortail wind) velocity m/s
`
`
`
`
`
`
`= vehicle mass, kg
`
`
`
`= modified vehicle accelerative mass, kg
`
`
`
`
`
`= time interval
`
`
`= chopperefficiency
`
`
`= battery charge efficiency
`
`
`
`= gearbox partialefficiency
`
`
`
`= battery discharge time at constant current
`
`
`
`
`
`
`I,,h
`
`= battery discharge time at constant power
`
`
`
`
`
`
`density Pp;,h
`
`
`= air density, kg/m*
`
`
`
`= hill severity, degrees
`
`
`
`= field flux, Wb
`
`
`
`
`Soc
`t
`T
`T,
`Tes
`Ty
`V
`Vi
`v,
`Vv,
`w
`Ww’
`At
`He
`Nee
`gb
`Tei
`
`Tpi
`
`p
`d
`br
`Suffixes
`max
`a
`f
`i
`5
`k
`pu
`
`= maximum
`
`= armature
`
`= field
`
`= discharge rate
`
`
`= 5 hourdischarge rate
`
`
`
`
`= step number
`
`
`= per unit
`
`
`NBcc denotes engine cubic capacity, 1 cc = 1 cm?
`
`
`
`
`
`
`
`
`
`
`
`1
`
`Introduction
`
`
`In 1976 the Department of Industry commissioned The
`
`
`
`
`
`
`
`International Research and Development Co. Lid., New-
`
`
`
`
`
`
`castle upon Tyne, to produce a worldwide survey of hybrid
`
`
`
`
`
`
`
`
`
`electric vehicles
`[1-3]. This report, prepared by A.J.
`
`
`
`
`
`
`
`Mitcham and J.R. Bumby, described numerous types of
`
`
`
`
`
`
`
`hybrid vehicles and their operating philosophies. It was
`
`
`
`
`
`
`
`evident that hybrid vehicles, and hybrid electric vehicles in
`
`
`
`
`
`
`
`
`particular, could be designed to meet a numberofdifferent
`
`
`
`
`
`
`
`
`
`265
`
`BMW1014
`Page 1 of 15
`
`BMW1014
`Page 1 of 15
`
`
`
`operating objectives. However, as to which was the best
`
`
`
`
`
`
`
`
`design to use to meet a particular operating objective was
`
`
`
`
`
`
`
`
`
`muchless obvious. Consequently, the report recommended
`
`
`
`
`
`
`that initially the performance potential of different hybrid
`
`
`
`
`
`
`
`vehicle designs should be examined using computer simu-
`
`
`
`
`
`
`
`lation.
`this
`this
`by
`paper
`recommendation,
`Stimulated
`
`
`
`
`
`describes a general road vehicle simulation program for
`
`
`
`
`
`
`
`evaluating the performance and energy efficiency of inter-
`
`
`
`
`
`
`
`nal combustion engined vehicles, battery electric vehicles
`
`
`
`
`
`
`or hybrid electric vehicles. This simulation package was
`
`
`
`
`
`
`
`conceived as a user-friendly interactive program, capable
`
`
`
`
`
`
`of evaluating the performance of a vehicle that is charac-
`
`
`
`
`
`
`
`
`
`teristic of a particular vehicle group, e.g. small car, medium
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`car, light commercial van ete. The program was given the
`name Janus. However, during the development of Janus,it
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`became evident that additional software options must be
`
`
`
`
`
`
`
`
`
`
`included to allow the user to define precisely the vehicle to
`
`
`
`
`
`
`
`be simulated, thereby extending the program to specific
`
`
`
`
`
`
`
`
`
`
`vehicle studies. Indeed,
`is this option that
`is used in
`it
`
`program verification.
`Since the publication of the Mitcham and Bumby
`
`
`
`
`
`
`
`the use of computer simulation has played an
`report,
`
`
`
`
`
`
`
`
`increasingly important role in evaluating the potential of
`
`
`
`
`
`
`
`new vehicle technologies, and in the design of specific vehi-
`
`
`
`
`
`
`
`
`
`cles [4-6]. A number of simulations have been written spe-
`
`
`
`
`
`
`
`
`
`cifically for one type of vehicle [7-10], while others, such
`
`
`
`
`
`
`
`
`
`as ELVEC [11-12] and HEAVY[13-14], are general road
`
`
`
`
`
`
`
`
`this
`vehicle
`simulation packages. The majority of
`
`
`
`
`
`
`published work has been in, and from, the USA,andthisis
`
`
`
`
`
`
`
`
`
`
`
`particularly true of the two general packages ELVEC and
`
`
`
`
`
`
`
`
`HEAVY,both of which are powerful design tools. One of
`
`
`
`
`
`
`
`
`
`these, HEAVY, has been mounted on the SERC comput-
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`ing network, but, at the present time, cannot be run inter-
`
`
`
`
`
`
`
`
`
`actively. In its conception, Janus was designed to be an
`interactive simulation package directly applicable to Euro-
`
`
`
`
`
`
`
`pean vehicles,
`In this paper, the characteristics of Janus are described
`
`
`
`
`
`
`
`
`in detail and used to demonstrate how a particular vehicle
`
`
`
`
`
`
`
`
`
`can be quickly and easily assembled from a standard sub-
`
`
`
`
`
`
`
`
`
`routine
`library.
`Initially the
`fundamental
`equations
`
`
`
`
`
`describing the vehicle dynamics are reviewed. In the final
`
`
`
`
`
`
`
`
`part of the paper, the use of Janus in simulating the per-
`
`
`
`
`
`
`
`
`
`
`
`formance of both electric and internal combustion engined
`
`
`
`
`
`
`
`road vehicles is demonstrated. In the case of the electric
`
`
`
`
`
`
`
`
`
`vehicle, this is achieved by simulating the General Electric
`
`
`
`
`
`
`
`
`ETV-1 electric car [15]; while when considering the inter-
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`nal combustion engined vehicle its usefulness in examining
`new vehicle concepts is demonstrated by investigating the
`
`
`
`
`
`
`
`influence of‘fuel off at idle’,
`
`
`
`
`
`
`2
`
`Vehicle dynamics
`
`
`To provide the necessary propulsion power, any vehicle
`
`
`
`
`
`
`
`drive train must be able to provide sufficient tractive effort
`
`
`
`
`
`
`
`
`
`at the road wheels to overcome aerodynamic drag,rolling
`
`
`
`
`
`
`
`
`resistance and hill gradient effects, while still providing the
`
`
`
`
`
`
`
`
`necessary vehicle acceleration. Consequently, at any parti-
`
`
`
`
`
`
`cular velocity and acceleration,
`tractive effort
`the net
`
`
`
`
`
`
`
`required at the road wheels can be expressed as the alge-
`
`
`
`
`
`
`
`
`
`
`braic sum of these components, i.e.
`
`
`
`
`
`
`
`
`
`
`E=]%+T+%Z+N
`where the respective componentsoftractive effort are
`
`
`
`
`
`
`
`Ty = 1/2pC, AV + VP N
`
`
`
`
`
`
`T,=C,WgN
`
`
`266
`
`(1)
`
`(2)
`(3)
`
`T, = Wg sin 6 N
`
`
`
`
`
`T.=WaN -
`
`
`and
`
`
`
`(4)
`(5)
`
`
`
`C, = drag coefficient
`
`
`
`A = vehicle frontal area, m?
`
`
`
`
`
`C, = coefficient of rolling resistance
`
`
`
`
`
`W = vehicle mass, kg
`
`
`
`
`@ =hill severity and percentage grade = 100 tan }, %
`
`
`
`
`
`
`
`
`
`
`V = vehicle velocity, m/s
`
`
`
`
`V,, = head wind velocity, m/s
`
`
`
`
`
`a = vehicle linear acceleration, m/s”
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`p = density of air = 1.226 kg/m? (at 15°C
`
`
`
`
`
`
`and 10° Pa (1 bar) ambient conditions)
`In eqn. 3 the coefficient of rolling resistance is dependent
`
`
`
`
`
`
`
`
`
`on the type of tyre used, the tyre pressure and the vehicle
`
`
`
`
`
`
`
`
`
`
`
`speed. However,
`it
`is difficult
`to quantify reliably for
`
`
`
`
`
`
`
`
`vehicle simulation studies, and a number of authors have
`
`
`
`
`
`
`
`
`used a quadratic, velocity-dependent equation for rolling
`
`
`
`
`
`
`resistance [7, 9, 11, 16], while others assumethe coefficient
`
`
`
`
`
`
`
`
`
`of rolling resistance to be a constant [17]. In general, the
`
`
`
`
`
`
`
`
`
`
`variation in the coefficient of rolling resistance with veloc-
`
`
`
`
`
`
`
`
`ity is small up to about 90 km/h, with the major changes
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`occurring at speeds significantly greater than this. Indeed,
`in tests on the General Electric ETV-1 electric car the
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`energy required to overcome
`rolling resistance was
`
`
`
`
`
`
`
`
`observed to decrease with speed [15], possibly due to
`
`
`
`
`
`
`
`
`
`increased tyre heating at higher speeds. To allow for simu-
`
`
`
`
`
`
`
`
`lation flexibility the coefficient of rolling resistance used in
`
`
`
`
`
`
`
`Janus contains both a constant and:a velocity-dependent
`term, both of which can be defined by the user.
`
`
`
`
`
`
`
`
`
`The accelerative tractive effort in eqn. 5 relates solely to
`
`
`
`
`
`
`
`
`
`the linear acceleration of the vehicle and takes no account
`
`
`
`
`
`
`
`
`
`of the rotational inertia of the road wheels and engine. The
`
`
`
`
`
`
`
`
`
`
`inertia of the road wheels can increase the effective vehicle
`
`
`
`
`
`
`
`
`
`weight by about 2% during acceleration, and is included
`
`
`
`
`
`
`
`
`in eqn. 5 by using an effective accelerative weight W’given
`
`
`
`
`
`
`
`
`
`
`by
`
`
`
`6)
`
`I
`
`
`
`
`
`W=W+ Ske
`YD
`where
`I, = inertia of the road wheels, kg m?
`
`
`
`
`
`
`
`
`rp = rolling radius of the wheel, m
`
`
`
`
`
`
`
`The engine inertia is not included directly in the tractive
`
`
`
`
`
`
`
`
`
`effort but as a local energy demand within the engine algo-
`
`
`
`
`
`
`
`
`
`
`rithm itself. It is interesting to note that, if referred to the
`
`
`
`
`
`
`
`
`
`
`
`wheels through the overall gear ratio, the engine inertia
`
`
`
`
`
`
`
`
`can increase the effective accelerative weight by about 20%
`
`
`
`
`
`
`
`
`in first gear and 2%in top gear.
`
`
`
`
`
`
`
`
`3
`
`
`
`Janus—its structure and operation
`
`
`
`
`the use of separate
`A fundamental feature of Janus is
`
`
`
`
`
`
`
`
`
`Fortran subroutines to represent individual vehicle com-
`
`
`
`
`
`
`ponents. The necessary subroutines can then be easily
`
`
`
`
`
`
`
`assembled into a master program to simulate a particular
`
`
`
`
`
`
`
`
`type of vehicle. This approach is adopted as many of the
`
`
`
`
`
`
`
`
`
`
`components required by the conventional
`internal com-
`
`
`
`
`
`
`bustion engined vehicle are encountered in an electric
`
`
`
`
`
`
`
`vehicle while components used in both these appear in the
`
`
`
`
`
`
`
`
`
`hybrid electric vehicle. Each component block is written to
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`ANSI Fortran VII standard. The master program can also
`
`
`
`
`
`
`include additional Fortran statements if required. Such
`
`
`
`
`
`
`
`
`
`
`versatility is of prime importance, as it allows the user to
`
`
`
`
`
`
`
`
`
`increase the degree of output information to suit his own
`TEE PROCEEDINGS,Vol. 132, Pt. A, No. 5, SEPTEMBER 1985
`
`
`
`
`
`
`
`
`
`
`BMW1014
`Page 2 of 15
`
`BMW1014
`Page 2 of 15
`
`
`
`particular requirements, or to restructure the control com-
`
`
`
`
`
`
`
`mands when simulating a hybrid vehicle.
`
`
`
`
`
`The advantages of this block structure approach can be
`
`
`
`
`
`
`
`
`demonstrated by considering the layout of the internal
`
`
`
`
`
`
`
`combustion engined vehicle shown schematically in Fig. 1.
`
`
`
`
`
`
`
`WHEELS
`
`COUPL
`
`TRANS
`
` ICENG
`
`1c engine [OI
`
`
`
`Fig.1
`
`
`
`
`Conventionalinternal combustion engine vehicle drive train
`
`
`
`
`
`
`
`The software blocks representing this vehicle are shown in
`
`
`
`
`
`
`
`
`the flow diagram of Fig. 2, and it is only the components
`
`
`
`
`
`
`
`
`
`
`
`within the dotted frame that are assembled by the user. To
`
`
`
`
`
`
`
`
`
`
`simplify program writing, the names of the software blocks
`
`
`
`
`
`
`
`
`have been selected to be similar to those of the com-
`
`
`
`
`
`
`
`
`
`
`ponents they represent. For operational
`flexibility each
`
`
`
`
`
`
`subroutine is divided into three sections termed, respec-
`
`
`
`
`
`
`
`the initial section, dynamic section and output
`tively,
`
`
`
`
`
`
`
`
`section. The initial and dynamic sections are further
`
`
`
`
`
`
`
`divided into two subsections. During the initialisation
`
`
`
`
`
`
`phase, each of the component subroutines is entered in
`
`
`
`
`
`
`
`
`turn, and the vehicle power train and driving cycle par-
`
`
`
`
`
`
`
`
`
`ameters
`specified. Unfortunately,
`in some subroutines
`
`
`
`
`
`information may be required from components upstream
`
`
`
`
`
`
`of the one currently being accessed, and as such will have
`
`
`
`
`
`
`
`
`
`
`yet to be defined. To overcome this problem a two pass
`
`
`
`
`
`
`
`
`
`
`system is used. On thefirst pass full details of the individ-
`
`
`
`
`
`
`
`
`
`
`
`ual componentsare specified, while on the second pass any
`
`
`
`
`
`
`
`
`
`calculations requiring information from an upstream com-
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`ponent are completed. Such a system is necessary,
`for
`
`
`
`
`
`
`
`example, when automatic weight generation is used. In
`
`
`
`
`
`
`
`
`such circumstances,
`the weight of the gearbox and final
`
`
`
`
`
`
`
`
`drive depends on the specified engine torque and power,
`
`
`
`
`
`
`
`one of the last componentsto be defined.
`Having established the initial operating conditions the
`
`
`
`
`
`
`simulation enters the dynamic section, the main computa-
`
`
`
`
`
`
`
`tional part of the program.In this part of the program, the
`
`
`
`
`
`
`
`
`
`
`
`fuel efficiency of the vehicle and all the component effi-
`
`
`
`
`
`
`
`
`
`ciencies are calculated. For the majority of driving cycles,
`
`
`
`
`
`
`
`
`such as the ECE-15 urban cycle shown in Fig. 3,
`the
`
`
`
`
`
`
`
`
`
`
`vehicle velocity is explicitly defined as a function of time,
`
`
`
`
`
`
`
`
`
`
`
`set t=0
`
`
`set flag for
`
`
`
`initial pass
`
`
`
`
`VEHICLE
`
`
`specifies all
`vehicle
`
`
`
`parameters.
`
`DCYCLE
`
`
`acceleration and
`
`velocity at time
`
`
`
`t for specific
`
`
`
`driving cycle
`
`
`
`WHEELS
`
`
`torque and speed
`reset
`initial
`at wheels
`
`
`
`
`conditions
`
`
`
`AXLE
`
`
`torque and speed
`at drive shaft
`
`
`
`
`account of axle
`
`
`
`
`
`efficiency
`
`
`TRANS
`
`
`transmission type,
`
`
`
`gear,torque and
`
`speed at engine
`
`
`
`
`COUPL
`
`
`friction clutch
`
`or torque con.
`
`
`
`
`ICENG
`
`
`fuel consumption
`over interval
`
`
`
`
`at
`
`yes
`
`
`
`set flog for
`set flag for
`
`
`dynamic initial
`
`
`dynamic
`
`
` Fig. 2
`
`IEE PROCEEDINGS,Vol. 132, Pt. A, No. 5, SEPTEMBER 1985
`
`
`
`
`
`
`
`
`
`
`BMW1014
`Page 3 of 15
`
`Flowchart for a conventional internal
`
`
`
`
`
`
`combustion engine
`
`
`267
`
`BMW1014
`Page 3 of 15
`
`
`
`enabling the program to step through the driving cycle at
`
`
`
`
`
`
`
`
`
`one second intervals (default value) calculating vehicle
`
`
`
`
`
`
`40
`
`60
`
`50
`
`30
`
`40r .
`g 20
`:
`= 30 E
`
`S
`G
`
`8778 4
`¥
`$
`10
`
`0
`
`o
`
`Fig.3
`
`
`
`
`
`
` ECE-15 urban driving cycle
`
`
`
`
`20
`
`40
`
`60
`
`8&0
`
`100
`
`time,s
`
`
`
`120
`
`
`
`140
`
`
`
`160
`
`
`
`180 200
`
`
`fuel used during the last time step and the energy loss in
`
`
`
`
`
`
`
`
`
`
`
`each of the drive-train components is recorded. The timeis
`
`
`
`
`
`
`
`
`
`then updated and the process repeated until the driving
`
`
`
`
`
`
`
`
`cycle is completed.
`
`
`the simulation
`Once the driving cycle is completed,
`
`
`
`
`
`
`
`enters the output section, where full details of the vehicle,
`
`
`
`
`
`
`
`
`
`driving cycle and the individual drive-train componentsis
`
`
`
`
`
`
`
`given with graphical presentation of time-varying quan-
`
`
`
`
`
`
`tities, engine fuel maps etc. Besides displaying component
`
`
`
`
`
`
`
`efficiencies, losses and the overall vehicle fuel economy, the
`
`
`
`
`
`
`
`
`percentage of the total cycle time spent in each area of the
`
`
`
`
`
`
`
`
`
`
`
`engine fuel mapis also given. Such fuel map informationis
`
`
`
`
`
`
`
`
`
`
`invaluable, particularly when detailed studies on the effect
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`of the vehicle component sizing and control] on fuel effi-
`
`
`
`
`
`
`
`ciencies are being undertaken. Alternatively, if a detailed
`
`
`
`
`
`
`
`component breakdownis not required, a reduced output
`
`
`
`
`
`
`
`
`facility may be selected which simply details the vehicle
`
`
`
`
`
`
`range and/or fuel consumption when operated contin-
`
`
`
`
`
`uously over a particular driving cycle.
`Once a simulation run has been completed, the control-
`
`
`
`
`
`
`
`
`ling software allows further runs to be conducted over a
`
`
`
`
`
`
`
`
`
`different driving cycle. Alternatively, modifications may be
`
`
`
`
`
`
`madeto the individual power-train components and/or the
`
`
`
`
`
`
`
`vehicle parameters.
`
`The different component subroutines available within
`
`
`
`
`
`Janus are listed in Table 1 and, in many cases, contain
`
`
`
`
`
`
`
`
`
`
`Table 1: Component subroutine names
`
`
`
`
`Component
`Simulation name
`
`_ Vehicle definition
`VEHICLE
`
`Oriving cycle
`DCYCLE
`
`Wheels
`WHEELS
`Final drive
`AXLE
`
`Transmission
`TRANS
`Clutch or torque convertor
`COUPLE
`
`
`
`internal combustion engine
`ICENG
`
`
`Series DC motor
`DCSER
`
`
`Separately excited OC motor
`DCSHUNT
`
`
`
`DC switched reluctance motor
`OCREL
`
`
`
`AC inductance motor
`ACINDUC
`
`
`DC generator
`OCGEN
`
`Field chopper
`FCHOPR
`
`Armature chopper
`ACHOPR
`
`Traction battery
`BATTERY
`
`‘Gearing’ far connecting two
`DRIVE
`
`
`
`prime movers
`
`Battery switching
`BATSWCH
`
`Summing block
`SUM
`
`Torque splitting module (hybrids)
`TOROQSPLT
`
`
`
`Powersplitting module (hybrids)
`POWSPLT
`
`
`
`
`Vehicle controllers VEHCONT
`
`
`performance details and simulation algorithms on a
`
`
`
`
`
`
`numberof related items, For example, the tree structure of
`
`
`
`
`
`
`
`
`
`Fig. 4 describes the specific components available within
`
`
`
`
`
`
`
`the transmission module, any of which can be selected
`
`
`
`
`
`
`
`
`from within the program.
`
`
`
`
`velocity and acceleration directly from the driving-cycle
`
`
`
`
`
`
`data. The tractive effort at the road wheels is calculated at
`
`
`
`
`
`
`
`
`
`
`each time instant using eqn. 1 and converted into a torque
`
`
`
`
`
`
`
`
`
`
`and rotational speed demand in the subroutine WHEELS.
`
`
`
`
`
`
`
`This torque requirement
`is then reflected back through
`
`
`
`
`
`
`
`each of the drive train components to the engine output
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`shaft. At each stage, account is taken of the gear ratio and
`instantaneous loss within each of the transmission com-
`
`
`
`
`
`
`
`ponents. Fuel usage is now obtained from the engine fuel
`
`
`
`
`
`
`
`
`
`map by assuming the engine load to be constant over the
`
`
`
`
`
`
`
`
`
`
`step interval. By sequentially repeating this process,
`the
`
`
`
`
`
`
`
`total fuel used over the driving cycle can be found.
`
`
`
`
`
`
`
`
`
`Such a direct calculation method is adequate for the
`
`
`
`
`
`
`
`
`majority of situations, as the performance of each drive
`
`
`
`
`
`
`
`
`train component can be considered to be independent of
`
`
`
`
`
`
`
`
`the other components. Normally the componentefficiency
`
`
`
`
`
`
`will depend only on the torque and speed it is required to
`
`
`
`
`
`
`
`
`
`
`
`transmit at that time instant. However, in some cases, the
`
`
`
`
`
`
`
`
`
`performance and behaviour of two or more of the drive
`
`
`
`
`
`
`
`
`
`train components may be
`strongly dependent. For
`
`
`
`
`
`
`example, in an electric vehicle, the behaviour of both the
`
`
`
`
`
`
`
`
`
`electronic controller and the traction motor depends on
`
`
`
`
`
`
`
`the battery terminal voltage. In turn, the battery terminal
`
`
`
`
`
`
`
`
`voltage depends on load current. In such instances, it may
`
`
`
`
`
`
`
`
`
`be necessary to iterate around the vehicle drive train until
`
`
`
`
`
`
`
`
`
`a stable operating condition is achieved, when component
`
`
`
`
`
`
`
`efficiencies can then be calculated. To accommodatethis, a
`
`
`
`
`
`
`
`
`two-pass system is used at each time step. On the first
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`pass, all the vehicle variables such as torques, speeds etc.
`
`
`
`
`
`
`
`
`
`are calculated from the road wheels to the energy source,
`
`
`
`
`
`
`
`
`using an iterative process to establish the stable operating
`
`
`
`
`
`
`
`condition. Once the stable operating condition has been
`
`
`
`
`
`
`
`
`
`reached the second dynamic pass is started, in which the
`
`TRANS
`
`gearbox
`
`CVT
`
`automatic
`(as manual)
`
`
`268
`
`manual
`(numberof gears: 2-6,
`
`
`
`gearratios, upchange
`
`
`schedule, downchange
`
`schedule)
`
`hydromechanical
`
`
`
`Perbury
`
`
`
`'V’ belt
`
`
`
`cone
`
`rofler
`
`
`constant
`efficiency
`
`Fig.4
`
`
`
`
`Subdivision oftransmission module
`
`
`
`IEE PROCEEDINGS,Val. 132, Pt. A, No.5, SEPTEMBER 1985
`
`
`
`
`
`
`
`
`
`
`BMW1014
`Page 4 of 15
`
`BMW1014
`Page 4 of 15
`
`
`
`and speed increments and stored as a 20 x 20 two-
`
`
`
`
`
`
`
`
`
`dimensional array. At each time interval through the cycle,
`
`
`
`
`
`
`
`
`
`140,
`
`1000
`
`%000
`
`5000
`
`4
`Componentsimulation
`
`
`417° Transmission
`
`the
`simulates
`The transmission subroutine (TRANS)
`
`
`
`
`
`behaviour of the vehicle transmission system and includes
`
`
`
`
`
`
`
`models of manual and automatic, fixed ratio, gearboxes
`
`
`
`
`
`
`
`and a variety of continuously variable transmissions
`
`
`
`
`
`
`(CVTs). The fixed ratio gearbox model is based on aneffi-
`
`
`
`
`
`
`
`
`
`
`ciency algorithm that predicts
`the gearbox efficiency,
`
`
`
`
`
`
`depending on the transmitted torque and speed. The algo-
`
`
`
`
`
`
`
`
`rithm also allows the number of gears, gear ratios and the
`
`
`
`
`
`
`
`
`
`
`gear change schedule to be specified. Both the up-change
`
`
`
`
`
`
`
`
`and the down-change schedule can be specified as a func-
`
`
`
`
`
`
`
`
`tion of vehicle speed or, alternatively,
`the ECE-15 gear
`
`
`
`
`
`
`
`
`change schedule maybe used. As the time taken during the
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`gear change is small compared with the total cycle time,
`
`
`
`
`
`the simulation assumes instantaneous gear changes.
`Theefficiency of a gearbox depends on the powertrans-
`
`
`
`
`
`
`
`
`mitted, the operational speed, the gear ratio and the gear
`
`
`
`
`
`
`
`
`
`3000
`2000
`
`speed rev/min
`profile and will vary with different gearbox designs.
`
`
`
`
`
`
`
`
`Although efficiency is high at normal full-load operating
`Fig.5=Enginefuel mapfor a typical 1850 cc engine
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`See Reference 18
`conditions, typically in excess of 96%, at part load, parti-
`
`
`
`
`
`
`
`
`
`
`
`@ Maximum torquecurve
`cularly high-speed part-load operation, efficiency drops
`
`
`
`
`
`b Constant SFC 400g/kWh
`
`
`
`and cansignificantly effect the vehicle fuel economy.
`
`
`
`
`
`
`
`
`
`
`c Power for constant 115 km/h
`
`
`
`
`
`ad Power far constant 80 km/h
`Gearbox lasses can be broadly divided into two com-
`
`
`
`
`
`
`
`
`
`
`
`
`
`e Power for constant $0 km/h
`ponents, a no-load loss due to gears churning in the
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`gearbox oil and a load friction loss due to the transfer of
`
`
`
`
`
`
`
`
`
`
`load between the gears. This last componentof lass is also
`
`
`
`
`
`
`
`
`
`
`dependent on the gear profile. Such losses are approx-
`
`
`
`
`
`
`
`
`imated in Janus by the power loss expression
`
`
`
`
`
`
`
`
`Poor = (1 — nealP; + 1.14 x 10782” W
`
`
`
`
`
`
`
`
`
`
`(7)
`
`where P, is the gearbox input power and n¢, is the partial
`
`
`
`
`
`
`
`
`
`
`
`efficiency relating to thefriction loss. The second velocity-
`
`
`
`
`
`
`
`
`dependent term is associated particularly with the churn-
`
`
`
`
`
`
`
`ing loss. The friction loss component is dependent on gear
`
`
`
`
`
`
`
`
`
`ratio and in the simulation is changed such that non-
`
`
`
`
`
`
`
`
`
`meshing gears,
`i.e. gear ratio of 1:1, have a higher effi-
`
`
`
`
`
`
`
`
`
`
`ciency than the meshing gears. Typically, 7g, is taken as
`
`
`
`
`
`
`
`
`
`0.99 for nonmeshing gears and 0.98 for meshing gears. For
`
`
`
`
`
`
`
`
`
`an automatic gearbox,a single value of 0.97 is assumed.
`
`
`
`
`
`
`
`
`At high rotational speeds, particularly at low load, the
`
`
`
`
`
`
`
`
`value of the churning loss speed index v becomescritical as
`
`
`
`
`
`
`
`
`
`
`gearbox efficiency is dominated by this loss term. To
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`obtain gearboxefficiency values typical of those quoted in
`the literature [17, 22], a speed index of 2.1 is assumed.
`
`
`
`
`
`
`
`
`
`The transmission is coupled to the engine by the
`
`
`
`
`
`
`
`
`COUPLEroutine, which contains models for both a
`
`
`
`
`
`
`
`torque convertor anda friction clutch.
`
`
`
`
`Efficiency curves similar to those described by eqn. 7 are
`
`
`
`
`
`
`
`
`
`also used in the final drive model. However, as the final
`
`
`
`
`
`
`
`
`
`
`drive may be constructed from either bevel or spur gears, a
`
`
`
`
`
`
`
`
`
`
`different friction loss componentis required for each gear
`
`
`
`
`
`
`
`
`type, bevel gears having a higherfrictional loss. The final
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`drive simulation block also contains efficiency values rep-
`resentative of both chain and belt drives, as these have
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`been used previously in electric vehicles (15, 20].
`
`4.2 Internal combustion engine
`
`
`
`Accurate engine maps giving the specific fuel consumption
`
`
`
`
`
`
`
`at different loads and speeds are essential if fuel consump-
`
`
`
`
`
`
`
`
`
`tion over urban driving cycles is to be accurately com-
`
`
`
`
`
`
`
`
`
`puted. Such fuel maps are commonly presented,either in a
`
`
`
`
`
`
`
`
`
`power/speed,
`brake-mean-effective-pressure
`(BMEP)/
`
`
`speed, or torque/speed form. An example ofthe latter is
`
`
`
`
`
`
`
`
`
`shownin Fig. 5 for a typical 1850 cc engine. For computa-
`
`
`
`
`
`
`
`
`
`
`tional purposes the fuel map is divided into twenty BMEP
`
`
`
`
`
`
`
`
`
`TEE PROCEEDINGS,Vol. 132, Pt. A, No. 5, SEPTEMBER 1985
`
`
`
`
`
`
`
`
`
`
`the fuel consumption is calculated for the particular load
`
`
`
`
`
`
`
`
`and engine speed by linearly interpolating between the
`
`
`
`
`
`
`
`four specific fuel values adjacent to the operating point.
`
`
`
`
`
`
`
`
`The position within the fuel map array is then recorded. At
`
`
`
`
`
`
`
`
`
`
`the end ofthe cycle simulation, the total fuel used is calcu-
`
`
`
`
`
`
`
`
`
`
`
`lated and the fuel usage information presented graphically
`
`
`
`
`
`
`
`in the form of a fuel usage map. This map displays that
`
`
`
`
`
`
`
`
`
`
`percentage of the cycle time the engine spent in different
`
`
`
`
`
`
`
`
`
`parts of the operating region. As engine efficiency maps are
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`obtained from steady-state load tests to obtain net engine
`
`
`
`
`
`
`
`torque during an acceleration interval, the engine inertia
`
`
`
`
`
`
`
`
`torque must be added algebraically to the output torque.
`
`
`
`
`
`
`
`
`
`As the engine inertia torque is proportional to the engine
`
`
`
`
`
`
`
`angular acceleration, this is readily achieved, except when
`
`
`
`
`
`
`
`
`
`
`a gear change-takes place. As a gear change is assumed to
`
`
`
`
`
`
`
`
`occur instantaneously, a step change occurs in the engine
`
`
`
`
`
`
`
`
`
`
`speed which can lead to very high, false, inertia torques. In
`
`
`
`
`
`
`
`
`
`general, engine inertia has only a small
`influence on the
`
`
`
`
`
`
`vehicle fuel economy;
`thereby, allowing engine inertia
`
`
`
`
`
`
`
`
`effects to be neglected during the computational timestep,
`
`
`
`
`the gear change takes place.
`the
`use of BMEP/
`that
`Experience has
`shown
`
`
`
`
`
`
`
`normalised speed fuel maps are to be preferred to power
`
`
`
`
`
`
`
`
`
`maps, as, with the latter, significant errors can be intro-
`
`
`
`
`
`
`
`
`
`duced at low-torque low-speed conditions; conditions that
`
`
`
`
`
`
`are commonin urban driving cycles. This arises because
`
`
`
`
`
`
`
`
`under such conditions the power demand is low, possibly
`
`
`
`
`
`
`
`
`less than 10% of maximum power, so that if equal power
`
`
`
`
`
`
`
`
`
`
`increments are used in the storage process, accurate spe-
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`cific fuel values corresponding particularly to the lower
`
`
`
`
`
`
`
`
`
`
`left-hand corner of Fig. 5 can be lost. By normalising the
`
`
`
`
`
`
`
`
`
`
`engine speed, the fuel maps can be ‘stretched’ to either a
`
`
`
`
`
`
`
`
`different engine cubic capacity, and hence power output, or
`
`
`
`
`
`
`
`
`to a different maximum engine speed other than their
`
`
`
`
`
`
`
`
`design value, Although stretching of the engine maps over
`
`
`
`
`
`
`
`
`too large an engine range, and particularly over a different
`
`
`
`
`
`
`
`compression ratio range,
`is not
`recommended,
`it does
`
`
`
`
`
`provide substantial additional
`flexibility with negligible
`
`
`additional programming complexity.
`At some stage during the driving cycle, it is possible for
`
`
`
`
`
`
`
`
`
`
`the demand torque to be greater than the fuil
`throttle
`
`
`
`
`
`
`
`
`
`269
`
`BMW1014
`Page 5 of 15
`
`BMW1014
`Page 5 of 15
`
`
`
`Fourtypes ofelectric traction motor are currently avail-
`
`
`
`
`
`
`
`
`able within the simulation;
`the DC separately excited
`
`
`
`
`
`
`
`motor, the DC series motor, the DC switched reluctance
`
`
`
`
`
`
`
`
`motor and the AC squirrel-cage induction motor. The DC
`
`
`
`
`
`
`
`
`series motor hastraditionally been used for electric vehicle
`
`
`
`
`
`
`
`
`traction as its
`torque characteristic, varying
`inherent
`
`
`
`
`
`
`inversely as a function of speed, best suits the requirements
`
`
`
`
`
`
`
`
`of the road vehicle. However, the series motor is not so
`
`
`
`
`
`
`
`
`
`
`well
`suited to regenerative braking, as
`the separately
`
`
`
`
`
`
`
`excited motor and also tends to have a lowerefficiency at
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`low speed. Consequently, with the advent of power elec-
`
`
`
`
`
`
`
`
`
`tronics there is an increasing use of the separately excited
`
`
`
`
`
`
`
`DC motor, particularly with transistor chopper control of
`
`
`
`
`
`
`
`
`field current and SCR armature chopper controlas, for
`
`
`
`
`
`
`
`
`
`example, used in the Lucas hybrid electric car [24], The
`
`
`
`
`
`
`
`
`
`advances in powerelectronics have also paved the way for
`
`
`
`
`
`
`
`more advanced traction motor types such as electronically
`
`
`
`
`
`
`commutated DC motors,
`the DC switched reluctance
`motor and the AC induction motor. In the Janus simula-
`
`
`
`
`
`
`
`
`
`tion the DC switched reluctance and the AC induction
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`motors are both represented as being favourable advanced
`traction motor drives.
`
`
`The torque/speed characteristic and efficiency map for
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`the separately excited motor used in the General Electric
`ETV-1 electric car is shown in Fig. 6 [25, 26]. Such a char-
`
`
`
`
`
`
`
`
`
`
`
`
`acteristic can be split into four operating regions; regions |
`
`
`
`
`
`
`
`
`
`and 2 associated with motoring action and regions 3 and 4
`
`
`
`
`
`
`
`
`
`
`with regenerative braking. When operating in regions |
`
`
`
`
`
`
`
`and 4 the field current, and hence flux, is maintained con-
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`stant and speed control of the motor is achieved by chop-
`
`
`
`
`
`
`
`
`
`
`ping the armature supply voltage. In regions 2 and 3 the
`battery is connected directly across the motor armature
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`
`terminals and speed control is ach