Advisor Document
Advisor Document
ME EN
RT
A
ER
DEP
GY
UNIT
ICA
ER
ED
ST A
AT E S OF
NOTICE
The submitted manuscript has been offered by an employee of the Midwest Research Institute (MRI),
a contractor of the US Government under Contract No. DE-AC36-99GO10337. Accordingly, the US Government and
MRI retain a nonexclusive royalty-free license to publish or reproduce the published form of this contribution, or allow
others to do so, for US Government purposes.
This report was prepared as an account of work sponsored by an agency of the United States
government. Neither the United States government nor any agency thereof, nor any of their employees, makes any
warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or use-
fulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe pri-
vately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trade-
mark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favor-
ing by the United States government or any agency thereof. The views and opinions of authors expressed herein do
not necessarily state or reflect those of the United States government or any agency thereof.
Printed on paper containing at least 50% wastepaper, including 20% postconsumer waste
Page
Opening A ddress
2000 ADVISOR Users Conference: Opening Address. . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
Keith Wipke, National Renewable Energy Laboratory
Post ers
A Logic-Based, Performance-Driven Electric Vehicle
Software Design Tool . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 498
David G. Alexander and Donald M. Blackketter, University of Idaho
Keith Wipke
Senior Engineer, Vehicle Systems Analysis Task Leader,
National Renewable Energy Laboratory
t o ols
3 .0 els sim.
of s 2. 0 I O R od er c e
t I R GU eb I S d M ot h r en
m en gin G U O
IS ed n w ADV rove with nfe
p be . 0 V c o p ng o
lo
ve SOR
1 AD han sed I m k i er s
C
e
D VI OR En lea in
VI
S e L s t Us
AD A D R 1
9 94 9 95 9 96 9 97 9 98 9 99 0 00
1 1 1 1 1 1 2
38 25 10
231
28
Africa
Asia
Australia
Europe
No. America
So. America
459
No Country
922
90
As of 6/7/00
NREL, CENTER FOR TRANSPORTATION TECHNOLOGIES AND SYSTEMS
ADVISOR Downloads
by Type of Organization
64 40
Industry
University
Govt/Labs
Other
597
1012
As of 6/7/00
NREL, CENTER FOR TRANSPORTATION TECHNOLOGIES AND SYSTEMS
ADVISOR Downloads by Industry
5 5
5 Ford Motor Company
5
5 47 DaimlerChrysler Corporation
6 General Motors
7 Visteon
Delphi Automotive Systems
6 Volvo Car Corporation
FIAT Research Center, Italy
6 Hyundai Motor Company
AVL
7 Parametric Technologies Corporation
Siemens Automotive Systems
8 29 AlliedSignal
Honda
Renault
9 Technologies M4 Inc.
Epyx Corporation
Hino Motors, Ltd.
12 Nissan Motor Company
Self-Employed
25 TNO Automotive
13
Legend includes industry with 5 or more users
13 Total number of industries using ADVISOR = 577
22 Total users in all industries = 1012
15
As of 8/18/00
NREL, CENTER FOR TRANSPORTATION TECHNOLOGIES AND SYSTEMS
ADVISOR Downloads by Universities
5
5 5 27
5
5 Ohio State University
6 University of Michigan
University of Maryland
6 University of Bath
George Washington University
6 26 University of California - Davis
Texas Tech University
6 Georgia Institute of Technology
Cornell University
7 University of Tennessee
Texas A&M University
Hanyang University
7 Pennsylvania State University
San Diego State University
22 Seoul National University
7 West Virginia University
Anna University
7 University of Colorado
University of Sheffield
Virginia Tech
7
16 Cranfield university
Institute for Advanced Engineering
8 MIT
University of Kent at Canterbury
University of Leeds
9 15
9
11 15
14 Legend includes universities with 5 or more users
Total number of universities using ADVISOR = 277
Total users in all universities = 597
As of 6/7/00
NREL, CENTER FOR TRANSPORTATION TECHNOLOGIES AND SYSTEMS
ADVISOR Downloads by State
11 11
13 12
15
17
14
23 211 Michigan
California
New York
26 Ohio
Texas
Colorado
27 Maryland
Massachusetts
Virginia
28 Indiana
Pennsylvania
Georgia
Tennessee
33 District of Columbia
Wisconsin
Illinois
Minnesota
36 Florida
110 Washington
36
Legend included states with 11 or more users
38 Total number of states using ADVISOR = 45
41 37 Total users in all states = 841
As of 6/7/00
NREL, CENTER FOR TRANSPORTATION TECHNOLOGIES AND SYSTEMS
ADVSIOR Downloads
by Government Agencies/Labs
2
2
8
2 Argonne National Lab
National Renewable Energy Laboratory
2 Oak Ridge National Lab
Natural Resources Canada
Politecnico di Torino
2 Southwest Research Institute
TASC
USDOE
2 California Air Resources Board
8 EPA
Los Alamos National Laboratory
3
Legend includes government agencies/labs
with 2 or more users
Total number of government agencies/labs
3 7 using ADVISOR = 64
Total users in all agencies/labs = 27
3
4
As of 6/7/00
NREL, CENTER FOR TRANSPORTATION TECHNOLOGIES AND SYSTEMS
Multiple Versions Downloaded as one
Indicator of “Active” Users
• ~20% appear to be “active” with ADVISOR
Number of ADVISOR Versions Downloaded (8/17/00)
2500
2078
2000
1500
# Users
1000
412
500
109
21 5
0
1 or more 2 or more 3 or more 4 or more 5 or more
versions versions versions versions versions
Series1 2078 412 109 21 5
# Versions
450 412
400
350
300
# Users
250
200
150 109
100
50 21 5
0
2 or more 3 or more 4 or more 5 or more
versions versions versions versions
Series1 412 109 21 5
# Versions
First to market . . . .
at low cost.
Quality Product
– Product quality is no longer an option, It is a way of life
– Quality is a measure of satisfaction of the customer’s
requirements & expectations
– Quality is dynamic (changing customer expectations)
Customer Loyalty
– 70 percent of customers bail because of the look/feel/smell/taste
of doing business with a company
Shareholder Value
– mission statement “ … to make a profit…”
The Process
The activities and procedures (and any sub
process) that supports the vehicle
conceptualization, it’s definition, design and
development including :-
– Targets process.
– Quality Operating System.
– Review & decision processes.
– Development process & procedures.
– Sign-off process, etc.
Classic Systems Engineering Process
Vehicle Reqmts
Develop &
System Reqmts Optimize Design
to meet Reqmts
Component Reqmts
Set, cascade
& balance
targets Verify & Confirm
design meets
requirements
A Successful Process requires….
• Customer requirements drive the process
• Well defined process is key to rapid product
development
• One process – no mavericks
• People must understand and use the process
• Compatible tools
• Integrated methods & techniques
• Web based communication
“Knowledge has become the key economic
resource and the dominant, if not the only,
source of competitive advantage”
Peter Drucker
• Virtual prototypes
• Tool selection for Integration
• Directional indications on time is better than
absolute results to late
• Avoid multiple masters - single data base
The People
“ In the digital age, as we move into
the quicker and quicker exchanges of
information . . . and re-inventions of
the world at work, our organizations
and our careers in action will become
more and more closely aligned with
the jazz ensemble . . . we will find
ourselves improvising with greater
and greater confidence and fearing
less and less the imaginative power of
the individual committed to enriching
the whole.”
– Team structure.
– Integrate product attributes & system design & release
responsibilities.
– Minimize the number of teams required.
– Minimize static data with real time communications.
– Provide central program data access.
The People
People Skills provide . . .
WORK TEAMS
Structure/Topology Team Interior Team Exterior Team Cockpit Team Vehicle Dynamics Team Powerpack Team
The People – Team Structure
C o n c e p t Ap p ea ra n c e D ev elo p m en t
C lim a te C o ntrol A rc h ite c tu re
A p p ea ra n c e C o lor H arm on y
W O R K TE AM S
A ttrib u te s C h a s s is N V H P ow e rtrain N V H
S afe ty
S te erin g
S ec u rity/A nti-T h e ft
B ra k ing
M o d u le T e a m s
B o d y S tru c ture In terio r T rim F rt E n d S ys In t. F rt S us p P ow e r U n it
C o m p o n e n t D e s ig n G ro u p s (C D G ) N a tu ra l W o rk G ro u p
C o m po n en t CAD M a nu fac tu ring CAE A ttrib ute A s s em b ly T es t
R e le as e E n gin e er D es ig n ers P ro c e s s E n gin e er E n gin e er D e v. E n gin ee r P ro c es s E n g ine e r E ng in ee r
The People & Vehicle Integration
Structure Topology
Power Propulsion
Vehicle Dynamics
Interior
Exterior
Cockpit
Multi-Functional Attribute Balancing
Manufacturing Shared Structural
Supplier Guidelines Knowledge Performance
t>2, R>1 (Local Buckling, tf/tw) Targets (m,Kt,Kb,ω
ω)
Economic Analysis
Cost (material, Beam Analysis
fabrication, assembly) Behavioral Design of
Modeling Experiments
Styling
Packaging
Optimum Beam Sections Topology
Optimum Connection Shape Optimization
New Approach
CAD Design
FEA / Optimization Meet Final Chassis
Body in White
Body in White Target CAD Design
“ Product quality requires managerial, technological
and statistical concepts throughout all the major
functions of the organization …”
Josheph M. Juran
Mechanism Motion
Dynamic Loading
Space Claim Envelop
John MacBain
August 24, 2000
Co-Simulation of
Page 1 ADVISOR and Saber
Agenda
◆ Overview
◆ Evolution of Automotive Architectures 1
– Traditional ala 2000 - sizing batteries and generators
– Traditional with increased electrical power budget
◆ Co-Simulation Concept
◆ Evolution of Automotive Architectures 2
– Dual Voltage Architecture
– Series Hybrid
– Parallel Hybrid
◆ DOE Contract for Co-Simulation of ADVISOR and
Saber
– Status Report
Co-Simulation of
Page 2 ADVISOR and Saber
Overview
Co-Simulation of
Page 3 ADVISOR and Saber
Traditional 14V Architecture
Drive
IC Engine Transmission Rear Differential Tires Cycle
Other Electrical
14V Battery Loads
Co-Simulation of
Page 5 ADVISOR and Saber
Traditional 14V Architecture
Increased Electrical Power Budget
◆ Three driving factors in future vehicle design
– Projected generator peak power requirements
» 1,800 watts in 1990
» 10,000 watts in 2020
– Tightening CAFÉ standards
– Tightening emission standards
◆ Implications:
– Events in the propulsion system and electrical system have
increasing impact on the other
– Simulations of the total power system of the vehicle, electrical
and propulsion, must be more interconnected as we design to
meet stringent requirements requiring system solutions on a
broader scope
Co-Simulation of
Page 6 ADVISOR and Saber
A Possible Solution
◆ Potentially Ideal solution - model electrical system in
MatLab/Simulink as a part of ADVISOR
Co-Simulation of
Page 7 ADVISOR and Saber
Co-Simulation Concept
Exchange Parameters at Each Time Gate
Saber
Independent Propagation
During Each Time Step
Co-Simulation of
Page 8 ADVISOR and Saber
Co-Simulation Concept
Traditional Vehicle Architecture
Potential Parameters to Pass:
ICE instantaneous rpm
Generator instantaneous required shaft torque
Saber
Propulsion System Propagation
Independent Propagation
During Each Time Step
Co-Simulation of
Page 9 ADVISOR and Saber
Dual Voltage Architectures
IC Engine Transmission Differential Tires Drive Cycle
ADVISOR
Co-Simulation Link
42V Generator Significant 42V Loads 14V Electrical
Related to Propulsion Loads
System
42V Battery
14V Battery
Saber DC to DC
Converter
Co-Simulation of
Page 10 ADVISOR and Saber
Series Hybrid Architecture
Either
ADVISOR
Saber
Battery Generator IC Engine
Co-Simulation of
Page 11 ADVISOR and Saber
Parallel Hybrid Architecture
Either
Significant Electric Loads
Saber
IC Engine Battery
ADVISOR
Co-Simulation of
Page 12 ADVISOR and Saber
Initial DOE Contract Activities
Co-Simulation of
Page 13 ADVISOR and Saber
Restart Demonstration of Saber
1st .5 seconds
Co-Simulation of
Page 14 ADVISOR and Saber
Restart Demonstration of Saber
1st .5 Seconds
Field Voltage at 14V
ADVISOR
MatLab /Simulink Saber
Parameter Exchange
by Direct Communication to
Command Lines
Co-Simulation of
Page 16 ADVISOR and Saber
More Realistic Ideal Co-Simulation Strategy
ADVISOR
MatLab /Simulink
Saber
Parameter Exchange
with AIM Script Shell
Co-Simulation of
Page 17 ADVISOR and Saber
1st Prototype Strategy for Co-Simulation
ADVISOR
MatLab /Simulink Saber
Data Set #2
Prototype AIM Script
Parameter Exchange validated in August by
via Data Sets Joe Conover.
Co-Simulation of
Page 18 ADVISOR and Saber
Conclusions
Co-Simulation of
Page 19 ADVISOR and Saber
Co-Simulation of ADVISOR and Saber - A Solution for Total
Vehicle Energy Management Simulation
ABSTRACT INTRODUCTION
Electrical analysis for vehicles with traditional 14V This paper will explore the largely independent
electrical architectures has often been limited to requirements on electrical system analysis today. Once
engineering functionality simulation for the components this baseline has been established, the paper will move
and systems. This level of analysis ignores the impact through a progression of vehicle architectures from
of the electrical system on the fuel efficiency and traditional architectures with higher electrical power
emissions of the full vehicle, a much larger system. budgets to future hybrid concepts with the perspective of
Even for traditional architecture vehicles that are requirements for electrical system analysis. This
prevalent today, the computational problems of increased electrical power budget comes from the shift
predicting propulsion system performance and electrical of loads from mechanical powered loads to the electrical
system performance are coupled. As the millennium system (EPS, EVA, etc.), the additional of new loads
changes, series and parallel hybrid vehicle concepts are (catalytic converter heaters), and hybrid vehicle
becoming a marketplace reality. These hybrid concepts architectures with electric traction motors. The
accentuate the coupling of the propulsion and electrical propulsion system and the electrical system become
systems making it less viable to simulate either ever more interdependent, moving electrical system
propulsion or electrical systems independently. This analysis from a largely independent activity today to an
paper explores the various goals of electrical system embedded activity with total vehicle propulsion system
analysis and how some electrical analyses require performance. The discussion establishes the argument
acknowledgement of the coupling of the electrical and for co-simulation between traditional electrical analysis
propulsion systems thus creating a requirement for software and traditional propulsion system analysis
simultaneous or co-simulation of the propulsion and software. The paper will end with a discussion
electrical systems. concerning the specific plans of the contract between
Delphi Automotive Systems and the National Renewable
This paper will lay out the plans for developing specific Energy Laboratory (Department of Energy funding) to
co-simulation technologies between NREL’s ADVISOR make co-simulation between ADVISOR and Saber a real
and Saber. This work represents the content of a and available option in future versions of ADVISOR.
contract between Delphi Automotive Systems and the
Department of Energy that establishes a partnership
between the National Renewable Energy Laboratory and
Delphi Automotive Systems. The contract is entitled TRADITIONAL 14V VEHICLES – SIZING
“Modeling and Simulation Development for Advanced ELECTRICAL COMPONENTS
Energy Management and Propulsions Systems.”
A traditional 14V system typical of most passenger car
The final product of the contract will be a co-simulation vehicles on the road today has a single voltage electrical
capability tying Saber (for the electrical architecture system which drives all the electrical loads on the
simulations) and ADVISOR (now for the purely vehicle. A schematic for traditional vehicles is displayed
propulsion system simulations) to produce total vehicle in Figure 1. The energy input to the electrical system is
energy management system simulations for traditional, through the shaft of the generator which is turned by the
series, and parallel hybrid vehicles. The work will internal combustion (IC) engine. The electrical
provide templates to accommodate both single and dual architecture consists of a generator (including voltage
voltage electrical architectures. The co-simulation regulator and rectifier), a storage battery, and a series of
capability will be available from NREL to the global electrical loads. This section will discuss this electrical
automotive engineering community much as ADVISOR architecture with the perspective of performing various
is available today with template interfaces for electrical analyses often faced in vehicle design.
constructing and simulating vehicle energy management
systems. A variety of electrical system problems are regularly
solved without acknowledging the larger vehicle
systems, and often this is well justified. Vehicle
designers always face the challenge of adequately sizing
the generator and the battery to satisfy the electrical
power demands of the vehicle. More detailed models of When a vehicle is sized for a generator and battery, the
the components are required for analyses which would engineer must select several typical driving cycles. The
indicate interference of signals between components engineer also selects several extreme cases such as
and/or the impact of anomalous voltage events triggered being stuck in rush-hour traffic. Typical load scenarios
by the switching of significant electrical loads. are then selected which would cause a strain for the
Challenging problems also exist in the choice of the electrical system. For instance, vehicles must continue
proper cables for the power and signal distribution to perform in the hottest desert environments in Death
system. One must carefully choose the cable gauge, Valley as well as harsh blizzard conditions in Alaska and
insulation class, and circuit fusing methodology (level Canada. Anticipated driver load choices (AC in summer,
and fuse/circuit breaker style) to insure proper and safe etc.) and drive cycle required loads (brake lights, turn
performance of the electrical system with adequate signals, etc.) must be included.
current and voltage being made available to the various
loads in the vehicle. The balance of this section of the The typical analysis methodology proceeds with an
paper will elaborate more on the proper sizing of independent propulsion system analysis showing the
generators and batteries since more awareness of these crankshaft rpm required for the vehicle to meet the
issues will provide the basis for later topics in the paper. requirements of the driving cycle. The crankshaft rpm is
All of these problems are readily addressed as converted to generator shaft rpm reflecting either a direct
independent electrical system analyses. These analysis mechanical drive system or a belt drive. This generator
technologies have already been created with a high level rpm is the sole input to the electrical system analyses.
of sophistication within Delphi Automotive Systems and The electrical system is analyzed for battery and
other companies. generator adequacy for the various load scenarios
coupled with the drive cycles and ignition system
Electrical loads can be categorized as loads which run electrical loading consistent with the drive cycle.
continuously during vehicle operation, loads which run
for prolonged periods, loads which typically run for The engineer adopts a margin of safety for the
briefer periods, and key off loads. Loads which run component specifications and the electrical system
continuously would include the ignition system, the typically works quite well. Analysis has played an
electric fuel pump, electric fuel injection systems, and important role in this engineering methodology for many
engine management computers / controls / sensors. years, and it has proven to be quite successful. For
Loads which run for prolonged periods would include the these purposes, independent electrical system analysis
windshield wipers, headlights, taillights, side marker is quite adequate.
lights, instrument panel lamps, entertainment systems,
engine cooling fans, and the HVAC blower. Loads with As Figure 1 demonstrates, the electrical and propulsion
briefer duty cycles include power windows, power door systems are indeed largely autonomous except for the
locks, power seats, horn, starting motor, brake lights, mechanical connection between the IC engine and the
backup lights, heated rear window, and the power generator where the engine speed is largely determined
antenna. These lists have been presented for by the drive cycle. The luxury of this relative autonomy
clarification purposes and are not exhaustive. is facing extinction. In Figure 1 as with all other figures,
the arrows represent flows of information in ADVISOR
and between ADVISOR and the electrical simulation.
Drive
IC Engine Transmission Rear Differential Tires Cycle
Other Electrical
14V Battery Loads
Now, let us move five years into the future when the Co-simulation is a computational strategy that repeatedly
power budget of a typical car has increased employs a simple algorithm that is very straightforward
substantially. The vehicle architecture may be much the to describe. The schematic of co-simulation is depicted
same for low to mid range vehicles. But, now the in Figure 2. Co-simulation is a process where the two
electrical budget will be much more significant in the solvers move forward independently through a time step.
total power budget of the vehicle. At this point, there is This time step may be user specified or it may be
no alternative to solving the electrical and propulsion theoretically controlled depending upon the dynamics of
problems simultaneously. This would be quite the vehicle system at that point during the drive cycle.
straightforward if all the models and/or differential The time step must be sufficiently small so that
equations resided in a single piece of software with a parameters of mutual interest do not experience
single solver. significant changes. At the end of the time step, key
information must be exchanged between the packages
to update the coupling or linking of the two solutions.
With the electrical analysis updated with refreshed The coupling becomes tighter with the 42V/14V vehicle
propulsion system parameters (IC engine rpm as certain 42V loads are added. Future 42V loads will
determining generator rpm) and the propulsion system certainly include electrical valve actuation (EVA), electric
updated with refreshed electrical system parameters brakes, heated catalytic converter, and electric power
(power drain from the IC engine to drive the generator), steering. These subsystems provide further couplings
a restart capability is exercised for each software between the propulsion system and electrical systems,
package to solve for next time step. This recurring representing a shift from traditional mechanically
sequence of parallel solutions through a time step powered subsystems to electrically powered
followed by an updating of parameters of mutual subsystems. The EVA system presents a speed
significance continues until the drive cycle has been dependent electrical load. Power steering and braking
completed. are considered propulsion for this discussion since they
traditionally draw power directly from the IC engine.
Exchange Parameters at Each Time Gate With this migration from mechanical power sources to
electrical power sources, the electrical power budget of
the vehicle may climb as high as 10 kW. The net result
is that independent propulsions system analyses and
Electrical System Propagation
electrical system analyses will be less representative of
reality and the requirement for co-simulation becomes
more significant.
Dual voltage electrical architectures face potentially The electrical and propulsion systems appear
more complex challenges than traditional 14V significantly more coupled than for the traditional vehicle
architectures. The first difference is the obviously more because the electrical system now embodies a
complex electrical system analysis. Figure 3 displays a significant component of the propulsion system. A
possible schematic for the electrical system. As a matter typical series hybrid will require co-simulation. The
of observation, such architectures will often include two control strategy / algorithm may well derive inputs from
batteries, battery control algorithms to enhance battery both the propulsion system and the electrical system.
life and performance, DC to DC converter(s), and a For instance, the control algorithm for the electric drive
much increased electrical power budget as a proportion motor may well have inputs from an automatic
of the total vehicle power budget. transmission to facilitate smooth shift points. Or, a low
DC to DC
Converter
Figure 3: 42V/14V Vehicle Architecture
charge state might force limitations in torque production PARALLEL HYBRID ARCHITECTURES
that would influence the transmission shift points (limp
home state). In this case, the coupling of the propulsion A parallel hybrid architecture is a propulsion system
and electrical systems requires co-simulation. Other which can draw motive power individually or in
options such as regenerative braking or electric power combination from two distinct powerplants. A possible
steering would also dictate co-simulation. configuration would be an IC engine and an electric drive
motor on the same drive shaft. The electric drive motor
I would be remiss if I did not observe that under certain would provide motive power in the regimes where the IC
limiting assumptions, the particular coupling of systems engine is less efficient. That same drive motor would
in a series hybrid can amount to complete decoupling of then serve as a generator to charge the battery while the
the problems from a computational standpoint. If the IC engine is providing the motive power. This same
control of the IC engine, the generation of electrical drive motor also serves as the flywheel and starter
power, and the battery control algorithm are dependent motor. Figure 5 depicts a typical parallel hybrid
solely upon electrical system parameters, and certain schematic.
electrical loads crossing between the electrical and
propulsion sides are not present (electric brakes, electric Parallel hybrid systems provide the most complete
power steering, etc.), then the problems can be coupling of the propulsion and electrical systems. A
completely decoupled. If one views the propulsion control algorithm shifts the motive power from the IC
system as everything downstream of the electric drive engine to the electrical drive motor to situations where
motor (not including the drive motor), then the necessary both are contributing. The control algorithm governs “IC
speed and torque profile for the electric drive motor can engine off” states, controls the smooth restart, and
be determined for the vehicle to match a drive cycle. A creates a seamless system for motive power generation
subsequent and completely independent analysis of the responsibility. Co-simulation is the only viable approach
electrical system (including the electric drive motor and if true simultaneous solution is not possible.
the IC engine) can determine the viability of the control
algorithm, battery, generator, and IC engine. Please
note that co-simulation is still possible, and the results
should be very comparable. So, the potential for
decoupling would not prevent full integration of the two
simulation techniques.
IC Engine Battery
Similar discussions could be written relative to electric All analyses to determine vehicle energy efficiency over
vehicles, fuel cell / electric series hybrid vehicles, and a drive cycle require careful analyses of the propulsion
the myriad of other nuances of hybrid concepts. The and electrical systems. The optimal situation is the
story will be the same for all. There is a deep coupling simultaneous solution of the equations governing all the
of the propulsion and electrical systems. Any attempt to processes describing the creation, conversion, storage,
analyze the efficiency of the entire vehicle relative to a and application of power. When the opportunity for
drive cycle will require simultaneous solution of all the simultaneous solution in a single solver becomes
relevant equations. When simultaneous solution is not impractical, then co-simulation involving multiple distinct
possible, co-simulation becomes a requirement. solvers becomes the option of choice. With sufficiently
small time steps, the solution resulting from co-
simulation should provide a very good approximation to
the true “simultaneous solution” results.
John Reuyl
Robert Apter
NEVCOR Inc.
NEVCOR, Inc.
Table 1 - Vehicle Fuel Efficiency (mpg)(1,2)
NEVCOR, Inc.
Table 2 - Fuel Efficiency of FCV, SEV, PEV and BEV
Compared to CV(1,2)
Drive Cycle FCV/CV SEV/CV PEV/CV BEV/CV
NEVCOR, Inc.
Figure 1 - Comparison of Range and
Refueling/Recharging Times for CVs, EVs and HEVs
NEVCOR, Inc.
Figure 2 - Distribution of Personal Automobile Use
100%
80%
Cumulative 60%
Percent
of Personal
Automobiles 40%
20%
0%
0 30 60 90 120 150 >155
Average Daily Travel Distance per Vehicle (miles)
Source: 1990 Nationwide Personal Transportation Survey (Ref. 6)
NEVCOR, Inc.
Figure 3 - Personal Automobile Miles Electrified
100%
80%
Vehicle 60%
Miles
Electrified HEV
(%) 40%
BOEV
20%
0%
0 20 40 60 80 100 120 140
Battery Range for BOEV and HEV (miles)
Source:1990 Nationwide Personal Transportation Survey (Ref. 6) and NEVCOR.
NEVCOR, Inc.
Control of NOx emissions is central to the control of
regional air pollution
NEVCOR, Inc.
California NOx standards actually lag the
accomplishments of the auto industry
0.40
Notes:
0.35
1) 1998 V-8 Mercury Grand Marquis is
0.30 a full-sized 6-passenger automobile.
2) California ARB does not require any
0.25 SULEVs until 2004.
NOx, g/mi
0.20
0.15
0.10
0.05
0.00
TLEV LEV/ULEV SULEV 1998 V-8
MERCURY
California ARB NOx Standards Adopted in 1998
NEVCOR, Inc.
ARB still considers the BEV to be the “Gold Standard”
• ARB considers the BEV to be the “Gold Standard” for two reasons:
– 1) inherent emissions durability, and
– 2) extremely low fuel-cycle emissions in California.
• However, the ARB admits that it is quite upbeat regarding
emissions durability of catalysts; ARB observes that: ...with proper
fuel tailoring adjustments, these latest technology palladium-
rhodium designs lose virtually none of their emissions conversion
capability over more than 100,000 miles of aging (1998 LEV II
amendments).
• Regarding low fuel-cycle emissions, the ARB is right; in fact, about
75% of all NOx emissions associated with electricity used by BEVs
in the SoCAB will actually be generated outside California, and
fuel-cycle emissions now dominate total emissions (next slide).
NEVCOR, Inc.
NOx emissions of hypothetical hybrid SUVs that would
match CA NOx of a BEV in the SoCAB in 2010
29 mi/gal Diesel Vehicle, 8.1 l/100 km Total World NOx other than US
Tailpipe NOx 0.01 g/mi.
NEVCOR, Inc.
Key Points
NEVCOR, Inc.
Comparison of Fuel Efficiencies and Fuel Flexibility of
Small Automotive Vehicles
John Reuyl and Robert Apter
NEVCOR, Inc.
In addition, it will be shown that so-called “grid- In addition, BEVs produce zero emissions at the vehicle,
connected” SEVs can reduce fuel use by a surprisingly an attribute that gained political prominence in California
large amount. in 1990 (see Ref. 3). During the past 10 years, since the
mandate of zero-emissions vehicles (ZEV) by the
California Air Resources Board (Ref. 3), BEV
development has been pursued by virtually every major
motor vehicle manufacturer.
Figure 1 - Comparison of Range and Refueling/Recharging Times for CVs, EVs and HEVs (from Ref. 10)
100%
80%
Cumulative 60%
Percent
of Personal
Automobiles 40%
20%
0%
0 30 60 90 120 150 >155
Average Daily Travel Distance per Vehicle (miles)
Source: 1990 Nationwide Personal Transportation Survey (Ref. 6)
The number of vehicles that travel short distances each improving utilization factors and utility efficiency.
day is surprisingly large. Even in the US where average
auto trips are the longest, Figure 2 (from Ref. 4) shows With its own on-board engine, such an HEV has no need
that more than 40% of all personal automobiles on the for daytime recharging or an expensive daytime
road on a typical day (roughly 40 million autos) actually recharging infrastructure. For days and even weeks at a
travel that day 20 miles or less. time, the engines in such HEVs may never be used, but
they would always be available, providing the same “full-
Such short distances are well within the feasibility of tank” feeling of security offered by the conventional
today’s EV and HEV technology using relatively small, vehicle.
low-cost battery packs.
Just as auto manufacturers offer a variety of
However, the majority of these same 40 million vehicles engine/transmission options for today’s buyer, such
will also travel much longer distances on weekends, on HEVs could be offered with a variety of battery and
vacation trips and on occasional side trips or engine options to meet tomorrow’s consumer choices.
emergencies during the work week. For example, the battery size and battery technology
could be selected based on each buyer’s likely average
The dilemma is that short-range BEVs cannot be used daily use. Buyers may choose engine size in much the
for these longer trips. Therefore, BEVs have been fitted same way as in a conventional vehicle; small engines for
with large, heavy (and expensive) battery packs in hopes the economy-minded, larger engines for those expecting
of providing enough range to satisfy at least some to haul trailers or do a lot of mountain driving.
potential users.
According to the research in Ref. 4, even HEVs with very
Yet HEVs can be optimally configured to meet both of small battery packs could electrify a very large fraction of
these very disparate duty cycles; small, low-cost battery total vehicle miles. “For any given battery range, HEVs
packs for the short trips and small engines that can could electrify, on average, more miles than BOEVs
provide the average power for unlimited long-range trips [battery-only electric vehicles] … because HEVs can be
whenever necessary. driven on all trips, and the initial miles every day will be
powered by utility electricity.” (Ref. 4, pg. 12). Figure 3
Figure 1 shows that recharging overnight for the short (next page) is reproduced from Ref. 4; it shows that,
trips would be possible from any 110V receptacle, compared to a fleet of conventional cars, a fleet of HEVs
without the infrastructure expense for a) dedicated 220V that plugged in every night and had only 30-mile battery
circuits, b) high-power chargers or c) elaborate custom packs could electrify roughly 50% of total vehicle miles
interconnect hardware and controls. Such overnight traveled.
charging would help to level utility loads, thereby
100%
80%
Vehicle 60%
Miles
Electrified HEV
(%) 40%
BOEV
20%
0%
0 20 40 60 80 100 120 140
Battery Range for BOEV and HEV (miles)
Source:1990 Nationwide Personal Transportation Survey (Ref. 6) and NEVCOR.
Even if fuel efficiency when using fuels were unchanged remain.... In the presence of anthropogenic NOx and
in these “flexible-fuel” HEVs from today’s fleet average of under favorable meteorological conditions, these
27.5 mi./gal, gasoline use would be reduced by about the background biogenic VOCs can contribute to
same 50% as a fleet of HEVs with doubled fuel efficiency summertime ozone concentrations exceeding the
(to 55 mi./gal) but no ability to recharge from the grid. NAAQS concentration of 120 ppb.” (Ref. 5)
Obviously, were fuel efficiency to be improved in these
HEVs (as is expected), fuel use would be further The California Air Resources Board (CARB) also
reduced. concluded that “...because ozone precursors, such as
NOx also react in the atmosphere to form particulate
The HEV also provides the same benefits of “gasoline matter (PM), reductions in NOx will be crucial to meet
equivalent fuel consumption” as the BEV when using existing state and federal PM10 standards, as well as the
grid electricity. Note from Table 1 that both the SEV and new federal standards for fine particulate matter
the BEV weigh the same 1030 kg. Thus, the SEV using (PM2.5).” (CARB, Ref. 8)
grid electricity would have roughly the same “gasoline
equivalent fuel efficiency” as the BEV when compared to However, “Except in California, NOx emission reductions
a similar CV traveling the same route. have not previously been a major component of most
[State Implementation Plans] SIPs.” (Ref. 5) As a result,
Note that the “grid-connected, flexible-fuel” HEV provides “Despite the major regulatory and pollution-control
the user with mobility even when fuel shortages occur. programs of the past 20 years, efforts to attain the
While longer weekend trips may be curtailed, the National Ambient Air Quality Standard of ozone largely
essential daily work-related trips can be undertaken have failed.” (Ref. 5)
using grid electricity.
Based upon the conclusions of the NRC, the authors
Such HEVs could vastly increase the market share for have focused their emissions work in the remainder of
EVs beyond that for just the BEV. HEVs that can this paper on NOx emissions.
recharge at night from the grid could be the most
practical way in the near-term of displacing transportation VEHICLE NOX EMISSIONS CAN BE REDUCED TO
fuels with utility electricity. NEGLIGIBLE LEVELS
HEVS CAN REDUCE OZONE-PRECURSOR California’s standards for motor vehicle NOx emissions,
EMISSIONS set by CARB, are actually lagging the auto industry’s
capabilities. For example, one of GM’s largest passenger
THE IMPORTANCE OF NOX CONTROL sedans in production, the 1998 V-8 Mercury Grand
Marquis, achieves a NOx level of 0.004 g/mi. as certified
In a major study of urban and regional air pollution, the by the CARB (see Ref. 8). Yet CARB’s tightest NOx
National Research Council (NRC) concluded that “Of the standard is the 0.020 g/mi. SULEV (Super Ultra-low
six major air pollutants for which National Ambient Air Emissions Vehicle) which is five (5) times more lax.
Quality Standards (NAAQS) have been designated under Furthermore, the SULEV standard is not required until
the Clean Air Act, the most pervasive problem continues 2004 and then for only a small fraction of all vehicles.
to be [tropospheric (i.e., near-ground)] ozone, the most
prevalent photochemical oxidant and an important A biennial review of the ZEV Program is scheduled by
component of ‘smog’.” (NRC, Ref. 5) the CARB for September 2000.
Ozone forms when NOx combines with VOCs (volatile Vehicle NOx emissions can now be reduced to such low
organic compounds) in the presence of heat and levels that NOx emissions from both electric- and fueled-
sunlight. VOCs originate from a variety of sources both vehicles would be dominated by the up-stream fuel-cycle
anthropogenic (e.g., paints, vehicle exhaust emissions) emissions. Such emissions are directly proportional to
and biogenic (e.g., plants, trees). NOx, on the other hand, the vehicle’s efficiency in using electricity and fuel and
occurs almost entirely from high-temperature combustion are discussed in the next section.
sources, and these sources are largely anthropogenic
(e.g., vehicle engines, power plants). FUEL-CYCLE EMISSIONS DOMINATE TOTAL NOX
EMISSIONS
The NRC determined that NOx control would be key to
ozone control. “NOx control is necessary for effective When exhaust emissions are at SULEV levels (0.02
reduction of ozone in many areas of the United States. g/mi.) and below, total emissions will be dominated by
...in many urban cores and their environs, even if the emissions associated with the fuel cycle itself. Ref. 9
anthropogenic VOCs are totally eliminated, a high reports that, even for vehicles used in California, total
background concentration of reactive VOCs will NOx emissions (vehicle plus fuel-cycle) for gasoline and
diesel vehicles can be equal to, and even less than, total Energy Laboratory (NREL). Key results from this
fuel-cycle NOx emissions for electric vehicles. The data research are presented in Figure 4 and summarized in
reported in Ref. 9 were developed and reported in Ref. this section.
10 under a research contract for the National Renewable
Figure 4 - Total NOx emissions of a) 1998 Ford Expedition, b) 5 mi./kWh BEV, and c) four hypothetical HEVs that
would have the same 0.022 g/mi. NOx emissions in California as the BEV (all vehicles are operated in
the SoCAB in 2010, Scenario 4) (from Ref. 10)
1)1)Fuel-cycle
Fuel-cycledatadataarearesc aled proportionally
scaled proportionally(based
(basedupon
uponfuel
fueleconomy)
economy)for forSScenario
cenario44from
fromRef.
Ref.7,7,
Evaluation of Fuel-Cycle Emissions...,
Evaluation of Fuel-Cycle Emissions..., 9/19/96. 9/19/96.
2)2)Graph
Graphdepicts
depictsfour
fourhypothetical
hypotheticalvehicles
vehicleswith
withtotal
totalNOx
NOxemissions
emissionsininCalifornia
Californiathat
thatwould
wouldequal
equalthethetotal
total
NOx emissions in California of a 5 mi/kWh EV used in the
NOx emissions in California of a 5 mi/kWh EV used in the S oCAB in 2010. S oCAB in 2010.
3)3)InInthis
thisgraph
graphforforSScenario
cenario44inin2010,
2010,tailpipe
tailpipeNOx
NOxemissions
emissionsforforthe
thehypothetic al fueled
hypothetical fueledvehic
vehic les
lesare
are
assumed
assumed to be either 0.01 g/mi. or 0.005 g/mi. (For comparison purposes, in ARB tests, the full-size1998
to be either 0.01 g/mi. or 0.005 g/mi. (For comparison purposes, in ARB tests, the full-size 1998
Mercury Grand Marquis V8 w ith an advanced catalyst achieved
Mercury Grand Marquis V8 with an advanced catalyst achieved 0.004 g/mi.) 0.004 g/mi.)
4)4)Note
Notethat
thattotal
totalNOx
NOxemissions
emissionsininthe theUS
USfrom
fromthethe55mi/kWh
mi/kWhEV EVarearemore
morethan
thandouble
doublethe thetotal
totalNOx
NOx
emissions in the US from any of the hypothetical fueled
emissions in the US from any of the hypothetical fueled vehicles. vehicles.
In order to determine the magnitude of fuel-cycle primary anthropogenic contributor to ozone (as
emissions, the CARB contracted with Acurex Corporation discussed in the sections above).
for an analysis specifically aimed at emissions associated
with vehicles when they were operated in the Southern The Acurex study concludes that BEVs do offer a unique
California Air Basin (SoCAB). The purpose of the study advantage in the SoCAB, due to regulations capping NOx
is quoted directly from Ref. 7, “This study investigated emissions for SoCAB power plants. “Power plants in the
the emissions associated with the production and South Coast Air Basin are subject to the Regional Clean
distribution of conventional and alternative fuels. Air Initiatives Market (RECLAIM) regulation of the South
Emissions from the production and distribution of fuels Coast Air Quality Management District (SCAQMD),
are known as fuel-cycle emissions and these can be which provides a cap on power plant NOx emissions for
significant in comparison to tailpipe emissions....The each utility. For larger usages of EVs in 2010, power
California Air Resources Board’s (CARB’s) zero-emission generation will result in zero additional NOx in the South
vehicle (ZEV) classification is based on zero emissions Coast Air Basin due to RECLAIM limits.” (Ref. 7)
from electric vehicles; therefore, fuel-cycle emissions
associated with incremental power generation might be However, Ref. 9 concludes that for the rest of California
compared with incremental fuel-cycle emissions from and for the US as a whole, gasoline and diesel-powered
conventional fuels.” (Ref. 7) HEVs with even modest fuel efficiency may offer total
NOx emissions that are even less than the total NOx
Although the Acurex study (Ref. 7) looked at both NOx emissions that are associated with the use of BEVs in the
and VOCs, the research reported in this paper focuses SoCAB.
upon NOx emissions because of their importance as the
Figure 4 (from Ref. 9, 10) shows total NOx emissions for benefit of proper fuel tailoring tends to verify this claim.”
vehicles operated in the Southern California Air Basin (Ref. 8)
(SoCAB) when using diesel, gasoline or electricity. In
addition to the tailpipe NOx levels (which are zero for an The CARB was quite correct that BEVs result in low fuel-
EV), the graph shows the fuel-cycle NOx emissions cycle emissions in California (Ref. 8). In fact, in 2010
associated with each fuel-type and whether those almost three-quarters of all NOx emissions associated
emissions occur a) in the SoCAB, b) elsewhere in with the use of BEVs in the SoCAB (Southern California
California, c) elsewhere in the US or d) outside the US. Air Basin) will actually be produced in other states (see
(In Figure 4, Scenario 4 has been selected for Figure 4).
presentation; Scenario 4 represents additional efficiency
improvements over the Scenario 2 “standards applicable However, fueled vehicles could offer the same low NOx
in 2010.”) emissions in California. For example, the top two bars of
Figure 4 show a 29-mi/gal diesel vehicle and a 46 mi./gal
The bottom bar of Figure 4 depicts the emissions of the gasoline vehicle that have NOx tailpipe emissions of
1998 Ford Expedition, a 3,000-kg sport utility vehicle 0.010 g/mi.; both vehicles achieve the same California
(SUV). It’s 0.03 g/mi. NOx emissions at the tailpipe as NOx emissions (0.022 g/mi.) as the 5 mi./kWh EV.
determined by CARB testing are 0.01 g/mi. above the
ARB’s SULEV standard of 0.020 g/mi. (the tightest Furthermore, total US NOx emissions in the rest of the
standard for fueled vehicles). However, because of its US from these fueled vehicles would be roughly 15% that
relatively low fuel efficiency (13 mi./gal), its total fuel- of the BEVs (see Figure 4).
cycle NOx emissions from Figure 4 are more than 0.23
g/mi., of which 0.073 g/mi. are in California. About half Such fuel efficiency for HEVs in 2010 seems quite
(0.118 g/mi.) occur outside the US (e.g., drilling, tankers feasible given the actual performance of current high-
at sea). The rest (0.044 g/mi.) occur inside the US but mileage diesel- and gasoline-powered prototypes.
outside California.
The assumed NOx tailpipe emissions of 0.010 g/mi. also
The bar immediately above the Ford Expedition seem quite feasible in 2010; according to CARB (Ref. 8),
represents the fuel-cycle NOx emissions of a 5-mi/kWh even the 1998, 6-passenger Mercury Grand Marquis with
BEV. Such efficiency is deemed feasible for a variety of a large V-8 engine has already achieved 0.004 g/mi. NOx
small 2-4 passenger BEV designs in 2010. Note that tailpipe emissions.
while such a BEV would have total NOx fuel-cycle
emissions of almost 0.09 g/mi., only 0.022 g/mi. occur in The middle two bars in Figure 4 depict hypothetical
California. vehicles that have only 0.005 g/mi. NOx tailpipe
emissions (this is still 25% more NOx than the actual NOx
In the next section, HEVs are described that could have emissions of the 1998 Mercury Grand Marquis). If these
the same low 0.022 g/mi. California NOx emissions as the hypothetical vehicles achieved only 20.5 mi/gal (diesel)
5 mi./kWh BEV. or 30.8 mi/gal (gasoline), they, too, would produce the
same California NOx emissions as the 5 mi./kWh BEV.
HEVS USING FUELS OR ELECTRICITY CAN MATCH
THE CALIFORNIA FUEL-CYCLE NOX EMISSIONS These conclusions are especially important given the
FROM BEVS growing popularity of sport utility vehicles (SUV) with
their low fuel economy and correspondingly high fuel-
The CARB states that “ZEVs are the ‘Gold Standard’ cycle emissions. The CARB (1998) LEV II amendments
based upon their extremely low fuel-cycle emissions in to the 1990 LEV Program subject these vehicles in 2004
California and inherent emissions durability” (Ref. 8). and beyond to the same tailpipe standards as passenger
cars, and their is every indication that HEV versions of
Emissions durability of fueled vehicles was an important SUVs could achieve SULEV standards for tailpipe NOx
issue during the early 1990s, and the “inherent durability” emissions. Furthermore, HEV versions of SUVs also
of the BEV ECS was an appealing image. could achieve much better fuel economy and, therefore,
much lower fuel-cycle NOx emissions and CO2 emissions.
More recently, the CARB has been quite upbeat about
the durability of the ECS (emissions control system) of HEVS CAN BE A SOURCE OF ELECTRICITY FOR
fueled vehicles: “Discussions with catalyst suppliers THE GRID
indicate that with proper fuel tailoring adjustment, these
latest technology palladium-rhodium designs lose HEVs may not be just the next generation of mobility.
virtually none of their emissions conversions capability The concept of HEVs as a source of power to the grid
over more than 100,000 miles of aging. The small was first patented in 1980 (Ref. 11). More recently, the
increases in emissions of our two Expeditions without the participants in the Electric Technology Roadmap Initiative
(Ref. 12) have concluded that Distributed Generation HEVS CAN DISPLACE FUEL USE WITH GRID
st
(DG) will be a cornerstone of the 21 Century electricity- ELECTRICITY
production system; as part of a DG system, the grid not
only can deliver power to HEVs for recharging; the grid BEVs displace fuel use by using grid electricity. BEVs
can receive power from the same HEV. offer the greatest “gasoline equivalent fuel economy” as
defined by NREL (see Appendix A), but their limited
Already, stationary “co-generation” systems have been range, lengthy recharge time and high cost have limited
commercialized, and HEVs are often “stationary” for long their marketability to special niches.
periods of time. It is a relatively small step to envision
that valuable assets like HEVs may provide more of a HEVs that can be recharged overnight from the grid also
return-on-investment than just energy-efficient mobility. displace fuel use; when using grid electricity for the local
commutes and around-town trips, such HEVs can match
CONCLUSIONS the BEVs “gasoline equivalent fuel economy.” Such
HEVs also can travel unlimited long distances by using
HEVs can reduce fuel use (and therefore emissions of fuels.
CO2 and NOx) by improving fuel economy and by
displacing fuel use with utility electricity. In the longer term, as renewable energy technologies are
commercialized, HEVs that can recharge from the grid
HEVS CAN REDUCE FUEL USE, NOX AND CO2 BY may get their energy increasingly from renewable energy
IMPROVING FUEL ECONOMY sources, thereby further reducing both carbon emissions
and ozone precursor emissions (e.g., NOx).
“Of the six major air pollutants for which National
Ambient Air Quality Standards (NAAQS) have been HEVS CAN BE A SOURCE OF POWER TO THE GRID
designated under the Clean Air Act, the most pervasive
problem continues to be [tropospheric (i.e., near-ground)] HEVs that can be recharged from the grid may also be a
ozone… (NRC, Ref. 5). The NRC determined that NOx source of power to the grid. Analysis reported in Ref. 13
control was key to ozone control; “… even if indicates that such HEVs represent a sustainable
st
anthropogenic VOCs are totally eliminated, a high technology for the 21 century. These conclusions have
background concentration of reactive VOCs will been corroborated in Ref. 12, a major policy study from a
remain.... In the presence of anthropogenic NOx and broad coalition of participants led by the Electric Power
under favorable meteorological conditions, these Research Institute (EPRI).
background biogenic VOCs can contribute to
summertime ozone concentrations exceeding the MANY HEV DESIGNS ARE LIKELY TO FIND
NAAQS concentration of 120 ppb.” (Ref. 5) PROFITABLE MARKET NICHES
Vehicle NOx exhaust emissions can now be reduced to Toyota’s Prius HEV offers the range, performance and
such low levels that NOx emissions attributed to the fuel comfort of a CV, and it has been a major success in the
cycle dominate total vehicle NOx emissions. By improving Japanese market. Toyota is now offering a version in the
fuel economy, HEVs can reduce total NOx emissions US that is optimized for the American market. Sales of
(exhaust plus fuel cycle) to levels comparable to the fuel- the Prius have exceeded the combined sales of all BEVs
cycle NOx emissions attributable to electric vehicles. from all manufacturers. Honda has entered the HEV
market with the Insight. Although these entries are both
Improvements in fuel economy also reflect directly into a hybrids, they differ significantly in their design and
corresponding reduction in the production of CO2, a operation. Both are competitively priced.
principal greenhouse gas.
In a new industry like HEV technologies and in the
Fuelcell hybrid vehicles (FCV) offer much improved fuel absence of a definitive vision of the future, there can be a
economy when compared to conventional vehicles (CV). tendency for engineering designers to grasp prematurely
If technical and cost barriers can be overcome, FCVs for “the best” HEV design. The authors argue, on the
would be effective in reducing fuel use, thereby reducing other hand, that many different HEV designs can emerge
CO2 and fuel cycle NOx. that could meet regulatory requirements and find
profitable market niches. The future bodes well for those
In the near-term, HEVs using combustion engines offer in the ADVISOR community who can develop these
similar advantages. Such HEVs improve fuel economy by many diverse options for their marketing colleagues.
operating engines only in their most efficient domain, by
capturing and storing braking energy, and by optimizing
each vehicle subsystem and minimizing vehicle weight.
ACKNOWLEDGMENTS 10. Reuyl, J. S., “The impact of HEVs (hybrid-electric vehicles)
upon regional air quality, climate change and sustainable
economic growth,” January 22, 1999, from Reuyl, J. S.,
The authors express their thanks to NREL for the “Why HEVs could be considered ‘Super-ZEVs’ and be
opportunity to present this paper at the ADVISOR Users worthy of multiple credits,” research report for the National
Conference, August 24-25, 2000, Costa Mesa, Renewable Energy Laboratory (NREL) of the US
California. Special thanks are extended to Bob Kost at Department of Energy (DOE). October 28, 1998.
DOE and Terry Penney at NREL for sponsoring the 11. Reuyl, J. S., US Patent No. 4, 182, 960, “Integrated
Residential and Automotive Energy System,” January 8,
research reported in References 4 and 10. 1980.
12. “Electricity Technology Roadmap,” C1-112677-V1, July
CONTACTS 1999, Electric Power Research Institute.
13. Reuyl, John S., “Hybrid-Electric Vehicles, Energy and the
st
Environment - A Sustainable Scenario for the 21
NEVCOR, Inc. th
Century,” Invited Paper for the 17 Congress of the World
Dr. John Reuyl Dr. Robert Apter Energy Council, September 13-17, 1998.
PO Box 18683 28 Ave. 4-Marronniers
Stanford, CA 94309 1400 Yverdon
US Switzerland
Tel/Fax:1 650 701 0464 Tel/Fax: 41 24 426 7160
e-mail: [email protected] e-mail: [email protected]
REFERENCES
The gasoline equivalent fuel consumption for an electric vehicle is calculated as,
mpgge = (distance traveled, mi) / (energy used, J) * ( lower heating value of gasoline, J/g) * ( density of gasoline, g/l) /
(3.785 l/gal)
where “energy used” equals the integral of (the power out of the energy system + discharge losses) / (coulombic efficiency
over all discharge periods. The coulombic efficiency is accounted for here because you would encounter it in recharging
your system back to its original state.
mpgge = (hydrogen fuel economy, mpg) / (density of hydrogen, g/l) / ( lower heating value of hydrogen, J/g) * ( lower
heating value of gasoline, J/g) * (density of gasoline, g/l).
The density of hydrogen in our data file is that of compressed hydrogen stored at 24 Mpa (~3500 psi).
Implementing Optimization in ADVISOR Using VisualDOC
Presentation
Paper
John Garcelon and Valdimir Balabanov,
Vanderplaats Research & Development, Inc.
■ System Responses
– objectives, obj(X): minimize / maximize
– constraints: gk(X) < 0
Optimization Applications
■ Structures
■ Nonlinear Mechanics
■ Computational Fluid Dynamics
■ Simulation of Kinematics
■ Chemical Processing
■ Financial Planning
What is VisualDOC?
■ General-Purpose Optimization System
– Add optimization to almost any analysis
■ Components
– Algorithm Library
– Response Program Interface
– Graphical User Interface
– Application Program Interface (API)
– Design Database
VisualDOC Components
■ Algorithms
– Gradient-Based (MMFD, SLP, SQP)
– Response Surface Approximate Optimization
– Design of Experiments
■ Response Program Interface
– Couple VisualDOC to Analysis Programs
■ Design Database
– Leverage Design Studies
■ API
– Embedding Optimization
Previous Efforts: ADVISOR &
VisualDOC
MATLAB - Workspace 1
Text
ADVISOR File
MATLAB - Workspace 2
VisualDOC
ADVISOR
Text
File
■ Auto-size
■ Responses
– Any combination of…
■ fuel economy
■ gradeability
■ Discrete Optimization
– Leverages design points in the database
■ Faster with fewer analyses
■ API
– Embed VisualDOC into ADVISOR
■ MATLAB calling structure
■ Direct access to the design database
– defining problems
– retrieving results
Embedding VisualDOC in
ADVISOR
Optimization
Problem
Definition
VisualDOC VisualDOC
ADVISOR
Optimization Database
Post-Processing
Design Database Concepts
Database
Interface
Inputs Responses
User
Etc
Control
state = INITIALIZE;
while state != DONE
state = RSA( X, R );
Modify_ADVISOR_Variables( X );
R = Calculate_ADVISOR_Responses();
Post-Processing Results
■ Database access functions
– Allow real-time monitoring of progress
– Complete response and design variable
histories
– Details about convergence, timing, errors
Implementation Schedule
■ New ADVISOR release: August 2000
■ VisualDOC 2 release: September 2000
■ Prototype ADVISOR-VisualDOC API
– fall 2000
■ comparison tests with previous coupling
■ generation and testing of new optimization
MATLAB - Workspace 1
Text
ADVISOR File
MATLAB - Workspace 2
VisualDOC
ADVISOR
Text
File
• Enforced Design Variable Constraints The mixed forward regression model in VisualDOC 2 will
• D-Optimal Designs generate mixed order approximations to minimize
• Mixed Forward Regression Models estimating error for a limited number of analyses. Thus,
mixed second order approximations can be generated
Enforced Design Variable Constraints for a 10 variable problem with significantly less analyses
than the 66 analyses required for a full second order
Designers can easily create irregular design spaces by approximation.
defining design variable constraints and enforcing them
in DOE. All minimum variance designs implemented in
VisualDOC place design points at the limits of design
variables. As a result, most analysis programs including
DISCRETE/INTEGER OPTIMIZATION The foundation of the API like the VisualDOC system is
the database. This means that to use the API,
The branch and bound process [7] for discrete/integer developers must use the VisualDOC database. This is
optimization in VisualDOC has also been improved and not as constricting as it sounds since the database is
takes advantage of the design database to eliminate very efficient. For example, an empty design database
reanalysis of previous design points. This significantly consists of two files, the data and index files, and is less
reduces the computational cost of discrete problems than 15 kilobytes in size. Most average design runs (10
such as auto-size where the number of battery modules design variables and 30 responses) typically require a
must be an integer. total database size of approximately 50 kilobytes.
Furthermore, the database provides a convenient way to
The VisualDOC API is also a new feature in VisualDOC organize and store information that is already required.
2 and it is explained in the next section.
The database organizes design data and provides fast
VISUALDOC API access to VisualDOC parameters and results. The
database is multi-user and platform independent.
The VisualDOC API is designed to allow application Therefore, simultaneous ADVISOR runs may access the
developers to incorporate VisualDOC into their same database, or the VisualDOC GUI could even be
applications. The VisualDOC API is a departure from the used to post-process design runs from ADVISOR.
architecture of previous versions of VisualDOC where
VisualDOC was designed to "wrap" around an analysis The VisualDOC database has three sections, interface
code. Although this wrapping was easy to implement, objects, design task objects, and results objects.
there were two significant drawbacks. First, users Interface objects act like a white board for defining
needed to invest time understanding optimization and design problems. Design task objects are a snapshot of
the VisualDOC GUI. Second, the implementation was the interface objects that define a design optimization
problem dependent; users could not easily develop a problem, and the results objects are the outcome from
general optimization capability. running a design task. Figure 2 depicts the database
sections.
ADVISOR provides an intuitive, well-organized interface
to hybrid vehicle powertrain analysis. Removing that The process of taking a snapshot verifies that a task will
interface and replacing it with a general purpose run successfully. Multiple tasks in a single database
optimization interface made using VisualDOC with collects all runs and their results in a single place,
ADVISOR more difficult. Today, most applications making the design study process more convenient and
include a GUI and the VisualDOC API allows developers less prone to loosing information. The inverse operation
to leverage their GUI and analysis capabilities with the of loading design task data and restarting is also
optimization API from VisualDOC. Users can still easily available.
wrap VisualDOC 2 around an application; however, the
API provides even more flexibility to developers.
Database
Interface
Inputs Responses
User
Etc
Control
• Defining/Modifying Interface Objects. Design modules operate with design task objects to
• Functions That Invoke Design Modules. create results. Design modules encapsulate the
• Retrieving Result Objects. optimization algorithms. There are three primary design
• Utility Functions. modules. These are Direct Gradient-Base Optimization
(DGO), Response Surface Approximations (RSA), and
INTERFACE OBJECTS Design of Experiments (DOE).
ADVISOR primarily works with interface objects in the DGO supplies the MMFD, SLP, and SQP algorithms.
database. Interface objects consist of the following Discrete optimization is integrated into this design
object types: module and is automatically applied when required.
The DOE design module also generates residual objects Figure 4: VisualDOC in the ADVISOR workspace
for all responses when it generates their approximations.
The residual objects help measure the accuracy of the OPTIMIZATION PROBLEM DEFINITION
approximations. There are over 25 attributes for each
residual object. The question that needs to be answered here is how do
we want to present general optimization problems to
UTILITY FUNCTIONS ADVISOR users? The goal is to do this in an intuitive
manner that does not disrupt the continuity of the familiar
The utility functions of the API are used to open and ADVISOR GUI, but still easily allows for the definition of
create databases, check database integrity, and recover optimization problems in ADVISOR. There are a number
data from corrupted databases. Functions also exist to of options here that we are considering since it will have
query the version and build numbers for each module impact on both development and use. Since the
including the database. VisualDOC API allows us to define the problem either
via the ADVISOR GUI (i.e., m-functions), the MATLAB
Because each language including MATLAB has different command window, and even using the VisualDOC GUI,
function calling structures, the VisualDOC API provides we have a great deal of flexibility when prototyping this
language specific calls. The API supports native interface.
MATLAB calls; therefore, ADVISOR can directly call the
VisualDOC API from both m-files and the MATLAB VISUALDOC OPTIMIZATION
command window.
Once the VisualDOC database contains a design task,
COUPLING WITH ADVISOR ADVISOR initiates the optimization run. Because
VisualDOC's design modules are simply another
Coupling VisualDOC and ADVISOR using the API may MATLAB function call, the no_gui_vrd.m file no longer
be schematically described in Figure 3. Each of the three exists. The optimization iteration is directly coded within
arrows connected to ADVISOR define interface points.
Optimization
Problem
Definition
VisualDOC VisualDOC
ADVISOR
Optimization Database
Post-Processing
Figure 3: Coupling VisualDOC and ADVISOR
Before addressing the interface points, it is important to ADVISOR. For example, the RSA design module simply
realize that using the VisualDOC API will eliminate the alters design variable values and requests that
extra MATLAB workspace shown in Figure 1. Since ADVISOR supply response values for these new design
ADVISOR directly calls the VisualDOC design modules, variables. The following pseudo-code shows the basic
VisualDOC shares the ADVISOR workspace (Figure 4). premise.
state = INITIALIZE; ACKNOWLEDGMENTS
while state != DONE
state = RSA( X, R ); The author would like to acknowledge the assistance of
Modify_ADVISOR_Variables( X ); the Systems analysis team in preparing this document
R = Calculate_ADVISOR_Responses(); and in performing the research described here. In
particular, the author would like to acknowledge the
efforts of Tony Markel.
Here, X and R are two MATLAB matrices. X represents
the current design variable values and R represents the CONTACT
corresponding responses that ADVISOR provides.
Modify_ADVISOR_Variables is a function that updates The author may be contacted by email at
the ADVISOR workspace to reflect the new design [email protected].
variable values, and Calculate_ADVISOR_Responses
computes the new responses according to the problem REFERENCES
definition.
1. Vanderplaats, G. N., Numerical Optimization Techniques
rd
for Engineering Design, 3 Edition, Vanderplaats
Because ADVISOR has many "dependent" variables Research and Development, Colorado Springs, CO, 1999.
(i.e., variables whose values dependent on other 2. ADVISOR 2.2 Reference, National Renewable Energy
variables), Modify_ADVISOR_Variables is not a trivial Laboratory, Golden, CO, October, 1999.
function nor is Calculate_ADVISOR_Responses. 3. VisualDOC 1.0 Reference, Vanderplaats Research and
However, we plan to use existing ADVISOR routines Development, Colorado Springs, CO, 1998.
with only minor modifications to generate the necessary 4. DOC 1.3 Reference, Vanderplaats Research and
Development, Colorado Springs, CO, 1995.
responses and to invoke the variable modifications for 5. Myers, R. H. and Montgomery, D. C., Response Surface
optimization. Methodology, Wiley Series in Probability and Statistics,
New York, NY, 1995.
POST-PROCESSING 6. Miller, A. J. and Nguyen, N-K., "A Fedorov Exchange
Algorithm for D-Optimal Design," Journal of the Royal
The Post-Processing connection is relatively easy. The Statistical Society, Volume, 43, No. 4, 1994.
7. Haftka, R. T. and Gürdel, Z., Elements of Structural
VisualDOC API provides functions for querying the Optimization, Kluwer Academic Publishers, Norwell, MA,
database as to the run status and providing intermediate 1993.
results for both inputs and responses. Real-time plotting
capability is available.
IMPLEMENTATION SCHEDULE
CONCLUSION
Exhaust Emissions (grams/mile) Differences in the vehicle thermal models can be seen to
Cold Hot Weighted % Error have an effect on the fuel economy and emissions of the
Average vehicle. The significant difference in the warm-up time
for the engine is the major reason for this. This is
Advisor
understood to be highly dependent on the thermal inertia
CO 12.63 0.728 5.846 244% of the thermal system, especially the volume of coolant
NOx 1.785 0.738 1.188 116% and the mass of the solids directly heated by the cooling
HC 1.981 0.328 1.039 352% system. Design and operation of the thermostat will also
FLOWMASTER2 / ADVISOR Co-Simulation effect the performance of the thermal system.
CO 3.88 0.729 2.084 23% For the hot start simulations, the emissions are seen to
NOx 1.63 0.738 1.122 104% be almost identical between the two cases. This is as
HC 0.807 0.329 0.535 133% expected, because as long as the thermal system can
EPA Test maintain a steady operating temperature, the emissions
of a known powertrain will not be effected by its design.
CO 1.7
Small fluctuations in the coolant temperature, due to the
NOx 0.55 lower system thermal inertia and detailed modeling of
HC 0.23 heat rejection, do not appear to have a noticeable impact
of the modeling of emissions. This detailed modeling is
valuable in that the simulation shows temperatures
above normal design maximums for cooling systems.
ADVISOR Cold Start Co-Simulation Cold Start ADVISOR Hot Start Co-Simulation Hot Start
120.0
110.0
100.0
90.0
Coolant Temperature (C)
80.0
70.0
60.0
50.0
40.0
30.0
20.0
1 101 201 301 401 501 601 701 801 901 1001 1101 1201 1301
Time (s)
ADVISOR Cold Start Co-Simulation Cold Start ADVISOR Hot Start Co-Simulation Hot Start
40000
35000
30000
Radiator Heat Rejection (Watts)
25000
20000
15000
10000
5000
0
1 101 201 301 401 501 601 701 801 901 1001 1101 1201 1301
Time (s)
ACKNOWLEDGMENTS
CONTACT
1
Interfacing of 1D and 3D Fluid Dynamics Programs for the Simulation of a Diesel Engine Cooling System, Keith Austin,
Development Manager, Flowmaster International Ltd
2
www.ctts.nrel.gov/analysis/advisor.html
3
www.mathworks.com/products/simulink
AAT
Vehicle Systems
18
16
14
Millions of Barrels per Day
10 Heavy Trucks
8 Gap
Passenger Vehicles
Light Trucks
6
4
2 Automobiles
0
1970 1980 1990 2000 2010 2020
Source: Transportation Energy Data Book: Edition 18, DOE/ORNL-6941, September 1998, and EIA Annual Energy
Outlook 1999, DOE/EIA-0383(99), December 1998
Dramatic Increases in the use of
AAT
Trucks in the U.S.
Vehicle Systems
160,000
140,000
120,000
100,000
Automobiles
Thousands
80,000 Trucks
60,000
40,000
20,000
0
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
19
19
19
19
19
19
19
19
19
19
19
19
19
19
19
19
19
19
19
19
19
19
19
19
19
19
19
Years
Source: Stacy C. Davis, Transportation Energy Data Book, Edition 18, Sept. 1998, ORNL-6941
Note: FHWA data were used. Starting in 1993, some minivans and sport utility vehicles that were previously
included with automobiles were included with trucks.
Future SUV? AAT
Vehicle Systems
Source: MH Designs, “The Future SUV”, The New SUV Standard, <http://www.poseur.4x4.org/futuresuv.html>,
(May 10, 1999)
Office of Transportation Technology
Objective AAT
Vehicle Systems
A More Energy
Independent
Increased
More Efficient Use of Nation
Vehicles
Alternative
Fuels
OTT/OAAT Is Pursuing Broad Range of
Advanced Technologies AAT
Vehicle Systems
Energy Fuels
Conversion
•Gasoline/ Diesel
•CIDI •Natural Gas
•Fuel Cell •Hydrogen
•SIDI •Dimethyl Ether
•VCR •Ethanol
•Fischer-Tropsch Fuels
Energy
Management Advanced
•Batteries Materials
•Flywheels •Metals
•Ultracapacitors •Composites
•Ceramics
Power Powertrain
Configuration Other
Electronics Attributes
•Inverters •Parallel Hybrid
•Series Hybrid •Accessory Loads
•Motors
•Generators •Electric Vehicle
•Conventional
Government/Industry Working
Together Under PNGV AAT
Vehicle Systems
Suppliers
Universities
Small Business
Daimler Federal Labs
Ford GM
Chrysler DOC DOE NSF DOT
Capabilities Technologies
Goal 1:
Adv. Manufacturing
DOD EPA NASA
Goal 2:
Near-Term Vehicle
Improvements
Goal 3:
Triple Fuel Efficiency
(up to 80 mpg)
OAAT R&D Plan: Strategy
“Systems Driven - Barrier Focused” AAT
Vehicle Systems
Ford P2000
$116M* $151M*
Chrysler ESX
* 50/50 cost
share $52.6M*
Methodology for Managing Vehicle Systems
Consists of Three Integrated Activities AAT
Vehicle Systems
Vehicle Systems
Dodge ESX3
- Body system weighs 46% less* Ford Prodigy Concept Vehicle
- Efficient diesel engine, motor and battery achieve 72 mpg* - Lightweight materials reduce vehicle weight 30% *
- Cost penalty halved to $7500 - Integrated starter/alternator *
- 33% reduction in aerodynamic drag
- Advanced diesel engine with 35% efficiency
2000 Detroit Auto Show improvement *
- High power battery *
PSAT
Ø Control-oriented, transients, many configurations
Ø Forward model
TESTING (APTF):
ØPerform the simulated tests
ØTest and map components
ØTest hybrid vehicles
– 48 parallel hybrids
– 24 series hybrids
Accelerator/Brake pedal
Motor command
TESTING (APTF):
ØPerform the simulated tests
ØTest and map components
ØTest hybrid vehicles
Simulated Actual
Vehicle DOWNLOAD Vehicle
Controller Controller
measures measures
commands
commands
Linked to PSAT
Driver
Dynamometer
Dynamometer
1-Simulation :
We simulate the test procedure with the
Test Procedure
Physical System Model composed by
the components library in order to check
the model.
Test
Procedure
Real World
Simulation
Control
Command dSpace
System
Analysis of Measured
Simulation Prototype
Physical Differences Results
Results Drivetrain
System
Model
Model Modification
clutch
Motor Eng
Dyno Dyno
45 kW 1.7 L
Phase 3 : Conventional Vehicle with CVT Phase 4 : Post-Trans Parallel HEV with CVT
Dyno Dyno
Eng Nissan 4:1 4:1
Eng Motor Nissan
1.7 L CVT 1.7 L CVT
45kW
This approach has been applied to the HEVs by using As PSAT is able to run both performance and
the bond graph methodology. Bond Graphs [Karnopp consumption/emissions tests, the user should also
et al., 1990] are graphical descriptions of dynamics choose what type of test should be done. In the case of
models based on power and informations flows. This an energy consumption test for hybrid configuration, a
technique offers the unique particularity to explicitly State-Of-Charge (SOC) equalization algorithm is also
describe not only the energetic exchanges between the available so that the consumption results of different
base elements of the physical structure of the system, configurations or strategies can be compared with same
but also the structure of its calculation. Moreover, the SOC.
Bond Graph, as it uses a pictorial form of the physical
structure and the calculated structure, appears as a link According to the user’s choices, the software will build
between the perception of the physical studied system, the appropriate model by using the right models,
its models, their exploitation on a microcontroller, and powertrain controller, and initialization files. The model
the interpretation of the results. will then be run automatically, and the results will be
provided in the Matlab prompt command.
In order to validate the models and their commands, we
use the Hardware-In-the-Loop (HIL) method with a CONVENTIONAL CONFIGURATION
Dspace 1103 board. This approach allows us to test and
validate the components and the whole drivetrain on the The main difference among conventional configurations
test bench. is the choice of transmission: manual, automatic, or
CVT. These drivetrains are used to validate PSAT’s
The main objective of this paper is to present the results with existing vehicles and serve as reference.
possibilities and characteristics of PSAT. We begin with
a description of the different drivetrain configurations that HYBRID CONFIGURATIONS
can be modeled. We then describe in more detail the
organization and the capabilities of the software. Finally, HEVs can be classified into two main architectures, as
using the description of the mechatronics approach, we shown in Figure 1 (mechanical or electrical power
explain the HIL methodology used to validate and
HEVs
Series
Parallel (torque addition) Power split (speed addition)
• Three possible transmissions (manual, automatic, The user can choose where to place the electrical motor.
and CVT), Four positions are proposed, as shown in Figure 2:
• Two different thermal sources, and
• Several different axle ratios (e.g., none, fixed gear). • Pos1: between the engine and the link element
(clutch/torque converter),
Parallel Configuration • Pos2: between the link element and the trans-
mission,
When we consider the parallel configuration, we can still • Pos3: between the transmission and the final drive,
divide it into two main architectures [Rimaux, et al., and
1998]: • Pos4: between the final drive and the wheel.
1 2 3 4
PSAT ORGANIZATION AND CAPABILITIES The output ports are used in the power controller (PTC)
for post-treatment, plotting or information in the strategy.
MODEL STRUCTURE The second ports carry the effort (i.e., voltage, torque)
and the last ones the flow (i.e., current, speed).
Organizational Format
To have an optimal reusability of the models, each
In order to exchange easily the models and implement component has its own:
new ones, we use a common format between the
input/output of the power ports by using the Bond Graph • Initialization file,
philosophy, as shown in Figure 5. The first ports are • Scaling file,
used for the information: • Model,
• List of variables (output information), and
• Input: on/off engine, gear number, etc. and • Selector file.
• Output (sensors): torque, rotational speed, current,
voltage, etc.
Command from Info to PTC
PTC
MODEL
Effort Effort
Flow Flow
Figure 5: Global Formalism for the I/O of the Models Using Bond Graph
ORGANIZATIONAL STRUCTURE parameters and variables by using the Tags), the name
of the Tags are defined as follows:
Use of Library
From_’variable name”
To ensure that the models we are using are the last Goto_”variable name’
ones changed or are not modified, we decided to use a
library in which all the models are saved. Libraries Use of Selectors
enable users to copy blocks into their models from
external libraries and automatically update the copied The Selector file allows us to parameterize the location
blocs when the source blocks change. of each variable in an array. Indeed, to carry sufficient
information, three buses are used: one for the
Use of Masks mechanical component, one for the transmission
component, and the last one for the electrical side. Every
We masked the models in order to use the parameters bus consists of the output variables of each component.
as local instead of global variables. This approach allows To find the place of each parameter in these buses, we
us to use the same generic model twice with two use some parameterized value to establish its location.
different data sets. We use the “GOTO-FROM” method For instance, the variable “nb_fc_spd_hist” is located in
th
to carry the information in an array from the model block the 4 position of the array “nb_thermic_variables”
to the “workspace block” where they will be transformed (which consists of the information coming from the
into global variables and thus be used in post-treatment engine, the clutch, and the exhaust).
processing.
Figure 8 shows an example of the use of selectors.
Figure 6 shows an example of a mask with an electrical
motor. NOMENCLATURE OF VARIABLE NAMES
Use of GOTO-FROM Format All the names of the software have been parameterized
and follow some rules.
To simplify the model, we decided to use the GOTO-
FROM format as shown in Figure 7. As far as the At the software level where the computations are made,
models are concerned, all of the GOTO-FROM blocks everything is based on the name of the component (e.g.,
are local and located at the upper level of the model (no ‘compo’ = fc for fuel converter). In fact,
blocks are located in the subsystems). Moreover, to
facilitate the work for HIL (Control Desk access to the • The component model name is defined as
‘compo’_cm (ex: fc_cm),
The same model can be used for
different components
Figure 6: Example of Using a Mask
Variable name
Goto_ptc_fc_trq_hot_mi n_cstr_hi st
Use of to_workspace to be
able to plot the variables
• The initialization file is ccording to the component in which they are used and
‘compo’_init, the type of data they represent.
• The scaling file is
‘compo’_scale, POST-TREATMENT
• The calculation file is
‘compo’_calc, At the end of a simulation, PSAT will display the results
• The parameter used to choose if we scale or not corresponding to user-defined preferences (performance
gui_scale_’compo’, or energy consumption and emissions). The user can
• The parameters used to scale the component also access the information of every model by using an
gui_’compo’, easy and automatic way to plot the variables used for
• The selector file associated with the model that run. For each component, we can plot and add as
selectors_’compo’_cm, and many variables as are located in the menu.
• The name of the main library of a component
lib_’compo’. Figure 9 shows the choices available to the user after
running a simulation with a parallel configuration in
At the component level, all of the variables and position 2 with a reduction block between the motor and
parameters also follow established rules and are named the main shaft.
Mechanical bus
nb_variable
Transmission bus
Electrical bus
Figure 8: Example of Use of Selectors
PSAT VALIDATION USING HIL • The first phase consists of validating one single
component (motor).
MECHATRONIC APPROACH • The second phase allows us to validate the engine
starting and clutching.
The mechatronic approach can be represented as • The third phase will be used to validate the CVT
shown in Figure 11. It consists of defining, according model and its ratio command.
to the customer’s expectations (e.g., acceleration, • The last phase will lead to the validation of the whole
maximum speed), the best HEV architecture and the drivetrain.
During each of the phases, we will use the measured not based on what software it uses but on how the
data (e.g., torque, speed, current, voltage, transient company uses it to accelerate the development of the
next generation of fuel-efficient, environmentally friendly
emission measurements) to come back to the simulated vehicles.
models to validate them.
REFERENCES
Figure 12 shows the different steps of the HIL validation.
Combes, E., and Cottard, C. 1992. Voiture particuliere
CONCLUSION Electrique Hybride (in English), P027, Advanced Electric Drive
Systems for Buses, Vans, and passenger cars to reduce
pollution. EDS, European Association of Electric Road Vehicle
PSAT is a user-friendly simulation toolkit that allows the (A.V.E.R.E).
user to simulate a large number of different HEV
configurations and implement the commands on Rimaux, S.; Delhom M.; Combes E.; and Rault, A. 1998.
prototypes by using an HIL technique. As the Hybrid Vehicle Powertrain: Modeling and Control.
government-financed PSATs develop, our goal is to
provide a nonproprietary version to a wide range of Karnopp, D.; Margolis, D.; and Rosenberg R. 1990. System
nd
people later in 2000. Thus, by incorporating the data, Dynamics: A Unified Approach, 2 ed. John Wiley & Sons,
models, and strategies of users, we can speed up Inc., New York.
PSAT’s development for everyone’s benefit and greatly
DeCharentenay F. Delhom M. and Rault A.. 1996. Seamless
strengthen this powerful simulation tool. In this world of
Mechatronic design of an electric vehicle powertrainm
increasing global competition, the value of a company is Convergence 96 pp413-421
Torque sensor
4 Dyno
4 Dyno :
: 1
Eng Mot C
Eng C 1 V
V T
T
George J. Delagrammatikas
and
Dennis N. Assanis
Automotive Research Center
The University of Michigan
• Application to ADVISOR
– Case Study Implementation
Motivation
• Federal Regulations
– Fuel economy (CAFE)
– Emissions (NOx, smog, and other pollutants)
• Public Awareness
– ‘Green’ movement
– Global warming scare
Objectives
• Design the optimal engine for a given vehicle
Background
– TDES
» Turbocharged Diesel Engine Simulation
(Assanis and Heywood, 1986)
» Zero-dimensional, quasi-static and feed-forward
» Enhanced to calculate an engine map for ADVISOR ‘on-the-fly’
– MATLAB
» SIMULINK - ADVISOR calculations
» Optimization framework - Sequential Quadratic Programming (SQP)
‘Library’ Approach
ADVISOR ENGINE LIBRARY
RUN ADVISOR
CHOOSE A DIFFERENT
OUTPUT METRICS ENGINE TO ACHIEVE
RESULT DESIRED OUTPUTS
DOES ENGINE
PRODUCE DESIRED NO
VEHICLE OUTPUTS?
YES
END
OPERATING
GEOMETRIES CONTROLS CONDITIONS
ENGINE
SIMULATION
CREATE MAPS ‘ON-THE-FLY’ CHOOSE DIFFERENT
ENGINE
CHARACTERISITCS
RUN ADVISOR TO ACHIEVE
DESIRED OUTPUTS
OUTPUT METRICS
RESULT
DOES ENGINE
PRODUCE DESIRED NO
VEHICLE OUTPUTS?
Model Decomposition
• Breaking up the large scale system into more
manageable subsystems
– Ensure conceptual simplicity of overall system
– Decrease computational expense
– Provide a methodical manner to perform trade-off
analyses
3 - Final Execution
1 ADVISOR
3 ADVISOR
2
BSFC TABLE REAL MPG,
MAX TRQ REAL PERFORMANCE
ENGINE
PARAMETERS
• Optimization
– Execute the decomposed, coordinated
subproblems and link their results
• Post-optimality
– Perform trade-off studies to understand the impact
of different parameters on the location of the
optimum design
200
PERFORMANCE
TORQUE
150 RUN
100
FUDS
50
CYCLE
0
0 500 1000 1500 2000 2500 3000 3500 4000 4500
SPEED
250
250
200
200
TORQUE
TORQUE
150
150
100
100
50
50
0
0
0 500 1000 1500 2000 2500 3000 3500 4000 4500
0 500 1000 1500 2000 2500 3000 3500 4000 4500
SPEED SPEED
REVERSE ENGINEERING A DIESEL ENGINE
arc Automotive Research Center NAC
TRQ4 D
TRQ1
CONSTANT
A
CITY HWAY POWER
4
1
PM
PM
PM
PM
R
R
TRQ2
TRQ3
1
TRQ4
4 2
TRQ1
5
1
4
PM
PM
PM
PM
R
R
REVERSE ENGINEERING A DIESEL ENGINE
arc Automotive Research Center NAC
TORQUE
TORQUE
RPM RPM
CORRECT INCORRECT
FINAL RESULT
TORQUE
SPEED
SPEED
REVERSE ENGINEERING A DIESEL ENGINE
arc Automotive Research Center NAC
ADVISOR Subproblem:
Problem Statement
minimize (Rated Torque + Rated Power) =
subject to:
0 – 60 mph 9 12 sec.
40 – 60 mph 9 5.3 sec.
Maximum speed , 100 mph
0 - 85 mph 9 24 sec.
5 sec. Distance , 140 ft.
max acceleration , 0.5 g’s
cruising grade (55 mph) , 6.5%
max launch grade , 30%
max difference between vehicle and driving cycle speed
traces 9 2 mph
number of gear shifts 9 180
TRQ3 9 TRQ2
TRQ4 9 0.7 (TRQ2)
TRQ1 9 0.5 (TRQ2)
TRQ2
TRQ3 VARIABLE LOWER BOUND UPPER BOUND
RPM1 750 850
RPM2 1800 2500
RPM3 2550 3000
RPM4 4350 4450
TRQ1 50 150
TRQ2 70 250
TRQ4 TRQ3 70 250
TRQ4 50 190
TRQ1
2
1
4
PM
PM
PM
PM
R
R
TDES Subproblem:
Problem Statement
subject to:
where i = 1 to 16
Variables:
PIMMAX
2.5
PINLET
2
PIM
1.5
0.5
0
200 5000
150
100
2000
3000
4000
WRPM
50 1000
TORQUE 0 0
RPM
2.5
PEXHAUST 1.5
ENGINE DISPLACEMENT 1
0
WRPM 200 5000
4000
PIM 150
100
3000
2000
PIMMAX 50 1000
RPM
TORQUE 0 0
TDES
200
OPTIMUM
180
TORQUE
160
140
120 ADVISOR
OPTIMUM
100
1000 1500 2000 2500 3000 3500 4000 4500
SPEED
200
231
BASELINE
180
246 277
TORQUE
160 262
140
120
OPTIMUM
100 324
80
356
60 293
434
309 371
402387 465 180
40 340 418 480
512
527 558605 652
449 496 543
574 683
590 621 730 824 715
699 761
777
793 746 808 933
839
637
20 668 855 871949886
902 980918 964 160
1000 1500 2000 2500 3000 3500 4000 262
140
SPEED
TORQUE
120 231
100 246
309
340
identical performance to 40
402
434
324
356
465
496
512 543
449
480
527558 590
621
637
668 699
574
605 746 808 855
777
730 761
683652 886
baseline 20
3500 4000
949
96
9
SPEED
REVERSE ENGINEERING A DIESEL ENGINE
arc Automotive Research Center NAC
180 136 kg
160 200 kg
140 250 kg
120 300 kg
100
80
600 1600 2600 3600 4600
ENGINE SPEED (RPM)
70
POWER (kW) 60
50 136 kg
40 200 kg
30 250 kg
20 300 kg
10
0
600 1600 2600 3600 4600
ENGINE SPEED (RPM)
Research Direction
• Emissions modelling
– Zero-, quasi - dimensional, and multi-zone models
• Transmission parameters
– Gear ratios and upshift/downshift strategies
• Trade-off analyses
– Vehicle mass, driving cycles, and product platform
design
• Methodology refinement
– Different parameterization techniques as well as
coordination methods
ADVISOR
2
An ADVISOR user investigating the effects of engine
INLET PRESS
100
3000 associated torque and speed scale parameters). This
2000
50 1000
procedure, illustrated in Figure 2 (a), is time-consuming
TORQUE
0 0
RPM and does not assure the user that the resulting engine
configuration is optimal, or even realistic. To guarantee
that the engine is a plausible design, we must introduce
(a) TDES as the engine model and link it, either directly or
indirectly, with ADVISOR.
2
subsystem could be performed because the overall
system simulation’s noise level would not adversely
1.5
affect the TDES calculations. The flow of this process,
1 illustrated in Figure 2 (b), is contrasted with the
0.5 traditional approach mentioned earlier in Figure 2 (a).
0
TDES outputs though, may still be masked by underlying
200 5000 vehicle simulation results if it were embedded within
4000
150
3000
ADVISOR as in [Fellini et al 1999][Sasena et al 1999].
100
2000 The following example illustrates this point.
50 1000
0 0
TORQUE RPM Let us consider that the engine designer would like to
understand the effects of varying bore size on the overall
fuel economy of a vehicle driven through the Federal
(b) Urban Driving Schedule (FUDS). Before TDES and
Figure 1: TDES parameterized inputs for inlet manifold ADVISOR were coupled, this type of study was not
pressures (a) and their associated, predicted exhaust possible. Now that the bore can be varied using TDES,
manifold pressures (b). an accurate prediction of engine characteristics can be
produced for an engine by only varying the bore size. A
For this study, TDES was developed to emulate variable new engine map could be produced ‘on-the-fly’ and sent
geometry turbocharging (VNT) with a wastegate. The directly to ADVISOR from TDES for system calculations.
parameters that describe the inlet manifold pressure
(PIM), the maximum boost pressure (PIMMAX), and the
engine speed at which the wastegate becomes activated
models while still keeping the information pertinent to
ADVISOR ENGINE LIBRARY their individual calculations consistent. The techniques
illustrated in the following sections will address these
issues.
YES
Depending on the modularity of the models found in the
large-scale simulations being investigated, the user will
find that one decomposition method may be more
END
attractive to implement than others. In this study, the
natural break-up of the subsystems, or object
(b) decomposition, was most suitable.
Figure 2: Typical engine design method using ADVISOR Figure 3 shows how ADVISOR and TDES were
currently (a). Technique implemented with an integrated separated and what variables linked the two simulations.
engine model like TDES (b). Since the programming styles of the two simulations
(feed-backward and feed-forward, respectively) did not
However, because the approximations present in allow the two simulations to communicate, a modified
ADVISOR’s correlations bring with them inaccuracies to Newton-Raphson techinque was developed to deliver an
the final output, how is the designer to determine the engine map automatically from TDES to ADVISOR in the
true effect of a slight change in bore to the output MPG? required data format [ARC98-1][ARC98-2].
Similarly, if these simulations were placed within a
gradient-based optimization routine, how could the Block 1 and 2 represent the two decomposed
optimizer converge properly when the derivatives found optimization subproblems (ADVISOR and TDES,
at each design point were affected by the variable respectively). In Block 1, ADVISOR interprets the
change and by an unknown noise or inaccuracy level in performance constraints placed on the vehicle as a
the simulation? There is now a need to separate these maximum torque profile required by the engine. Upon
finding the optimum curve to satisfy the vehicle mission, that is achievable by a real engine requires interaction
TDES matches the torque curve in Block 2. After finding with TDES.
the most fuel-efficient engine to do so, TDES then
calculates the entire BSFC map and supplies Block 3 250
with the data. In the final block, ADVISOR simply PERF
calculates the true fuel consumption for the cycle using
200
the data from TDES.
150
1 - ADVISOR Subproblem
VEHICLE
PARAMETERS 2 - TDES Subproblem 100
3 - Final Execution
50
1 ADVISOR
3 ADVISOR
0 FUDS
2 0 500 1000 1500 2000 2500 3000 3500 4000 4500
BSFC TABLE REAL MPG,
MAX TRQ REAL PERFORMANCE
T O R Q U E S (N-m)
PARAMETERS
200
100
0 – 60 mph ≤ 12 sec. 0
0 500 1000 1500 2000 2500 3000 3500 4000 4500
FUDS
40 – 60 mph ≤ 5.3 sec.
Maximum speed ≥ 100 mph
0 - 85 mph ≤ 24 sec.
5 sec. Distance ≥ 140 ft.
max acceleration ≥ 0.5 g’s 250 PERF
cruising grade (55 mph) ≥ 6.5 %
max launch grade ≥ 30%
200
engine does not yet exist for our purposes and, realizing
that the vehicle must have adequate acceleration to be 100
acceptable, we parameterize the maximum torque curve
of that engine. This curve then replaces the data for the
50
engines previously used within ADVISOR and the
associated BSFC and emissions calculations are
0
suppressed. Thus, we can convert the parameters that 0 500 1000 1500 2000 2500 3000 3500 4000 4500 FUDS
define the engine’s torque curve into the variables used
in parametric studies or in an optimization routine. ENGINE SPEED
4
PM
M
by these use points are noted for each cycle.
RP
RP
RP
R
In order to attain the required acceleration demands, line
CD was constrained to be at most tangent to the Figure 6: Representative maximum torque curve and the
constant power curve of approximately 70 kW (the effects of different constraints on the heights and slopes of
vehicle that was selected needed only 66 kW to obtain its three segments.
the performance requirements). This line would be
defined by the minimum rated engine power required by Based on the claim that for the same engine use points,
the performance run. The engine’s required torque a large engine would yield higher fuel consumption than
supply was defined by the height and length of line BC. a similar quality small engine, an optimization framework
Finally, line AB was varied until a minimum required was developed to minimize both these quantities in a
distance was attained between it and the closest use multiobjective function without any weight bias. Further
point(s), that is, the points werer not allowed to exist studies concerning weights are left for future work.
above this line. For AB, the transmission parameters
(denoted by TRANS) at the low speed range were the
driving force behind this line’s slope and height. TRQ2
TRQ3
Note that many maximum torque curves were found that
could satisfy performance mathematically. Figure 7 3
shows how widely these torque curves could vary and
how unrealistic, from an engine perspective, several of
these curves were. Lines 1 and 2 are both realistic 1
curves that can be attained by the engine; we will see TRQ4
shortly which one is more desirable. Lines 3, 4, and 5 4 2
cannot be matched to an engine and are simply TRQ1
parameter combinations at local minima within the
design space. In order to understand the typical shapes 5
of real engine torque curves, TDES was then
implemented in a series of pre-optimality studies.
1
4
PM
PM
PM
PM
R
PM
PM
PM
Torque (N-m)
The TDES subproblem is set up in the following manner,
shown in Table 4. Note that the different sets of curves
reflect the history of the optimization run in Figure 9 (a),
and (b) shows the final output. Let us now perform the
optimization study in the next section.
w h ere i = 1 to 1 6
Torque (N-m)
50 136 kg
Baseline
40 200 kg
30 250 kg 220
20 300 kg TDES
10 200 Optimum
0
600 1600 2600 3600 4600
Torque (N-m)
180
ENGINE SPEED (RPM)
160
(b)
140
Figure 10: Results from post-optimality parametric study
on cargo mass. (a) shows the maximum torque curves of ADVISOR
120
the optimal engines and (b) depicts their corresponding Optimum
power curves.
100
1000 1500 2000 2500 3000 3500 4000 4500
The parameterization method for the maximum torque [Fellini et al 1999] Assanis, D., Delagrammatikas G.,
curve in ADVISOR will also be researched. The Fellini R., Filipi Z., Liedtke J., Michelena, N.,
sensitivity of the objective function in this subproblem Papalambros, P., Reyes, D., Rosenbaum D., Sales A.,
with respect to the ERPM points has raised several and Sasena M., "An Optimization Approach to Hybrid Electric
concerns. Fourth-order quadratic functions, which can Propulsion System Design," Journal of Mechanics of
afford rather flat plateaus and the necessary slopes in Structures and Machines, Vol. 27, No. 4, pp. 393-421,
the low and high speed ranges, are potential candidates 1999.
for this study. In addition, the PIM, PIMmax, and WRPM
parameterization scheme for the TDES subproblem will [Heywood 1988] J. Heywood, “Internal Combustion
also be modified to make the boost scenarios predictive Engine Fundamentals,” McGraw Hill, Inc., New York,
(and more realistic) with turbocharger maps or perhaps NY, 1988.
neural networks.
[Sasena et al 1999] Fellini, R., Michelena, N.,
Next, efforts will be made to incorporate an emissions Papalambros, P., and Sasena, M., "Optimal Design of
modeling and optimization strategy to the overall design Automotive Hybrid Powertrain Systems," Proceedings of the
framework. Though preliminary NOx and dry particulate EcoDesign '99: 1st International Symposium on
matter models have already been implemented within Environmentally Conscious Design and Inverse
TDES, they have not yet been validated and therefore Manufacturing, Tokyo, Japan, February, 1999.
not presented in this paper. Methods similar to those
described earlier will be implemented for this puprose. [Michelena et al 1999] N. Michelena, C. Scheffer, R.
Fellini, and P. Y. Papalambros, ‘A CORBA-based
Finally, this work is the precursor to an optimization Framework for Distributed, Object-Oriented System
strategy for the parallel HEV which will require a greater Design,’ Journal of Mechanics of Structures and
degree of coordination between submodels. Machines, Vol. 27, No. 4,1999.
Transmission and motor design and control quantities
are being screened for their potential incorporation within [Millington and Hartles 1968] B. W. Millington, E. R.
a scheme similar to the one presented above. Hartles, ‘Frictional Losses in Diesel Engines,’ SAE Paper
680590, SAE Transactions, Vol. 77, 1968.
ACKNOWLEDGEMENTS
[Papalambros et al 1997] Wagner, T., Papalambros, P.
The authors would like to thank the generous Selection Families of Optimal Engine Designs Using
contributions of the members of the Automotive Nonlinear Programming and Parametric Sensitivity
Research Center at the University of Michigan. This Analysis. SAE Paper # 971600.
research could not have been performed without the
assistance of Panos Papalambros and his group in the [Papalambros and Wilde 1991] P. Y. Papalambros, D. J.
area of design optimization techniques for large-scale Wilde, Principles of Optimal Design, Cambridge
systems. The works and advice of Nestor Michelena, University Press, New York, 1991.
Zoran Filipi, Michael Sasena, and Ryan Fellini were
influential in our understanding of the issues relating to
[Scheffer 1997] C. Scheffer, “A Software Framework for
Optimal Systems Design,” Master of Science Thesis,
Department of Mechanical Engineering and Applied
Mechanics, University of Michigan, Ann Arbor, MI, 1997.
of ICE
EM
same time
• 2 or 4-wheel-drive layouts
available
• Pure electric operation into BP: battery pack; EM: electric motor; INV: inverter;
MT: internal combustion engine; PGS: planetary gear set;
restricted zones EG: electrical generator; VR: voltage regulator
Planetary Gear Set links
• ICE to the Sun Gear
• Generator to the Planetary Sun gear
Carrier
• Driving Shaft to the Ring to the electrica
Gear Planet gears
generator shaf
with carrier
1 τ −1
ω engine = ⋅ ω driving _ shaft + ⋅ ω generator
τ τ Internal
to the i.c. ring gear
engine shaft
ω driving_sh aft = τ ⋅ ω engine + (1 - τ ) ⋅ ω generator
τ 1
ω generator = ⋅ ω engine − ⋅ ω driving _ shaft
τ −1 τ −1
• The electric motor driving
the rear wheels operates in
the 1st (traction motor) or
4th (regenerative braking)
quadrant
• The electrical generator
connected to the PGS
planetary carrier operates in
the 1st/2nd or 3rd/4th
M: traction torque; MR: resistance torque; quadrants, depending on
N: rotation speed
the sign of the PGS
conversion ratio τ
gal
e x_calc
t total fuel used (gal)
Clock To Workspace AND
fc_emis
emis
HC, CO,
ve hicle controls e xha ust sys NOx, PM (g/s)
sipre hybrid
<vc> control sta te gy <ex>
ICE fuel
<cs>
conve rte r
<fc>
<sdo>
P roduct4
Ground
drive cycle ge a r ge ne ra tor/
<cyc>1 re duction controlle r <gc> powe r e ne rgy
bus <pb> stora ge <e s s>
ve hicle <ve h>1 pla ne ta ry
front whe e l a nd front
ge a r se t <pgs>
a xle <wh> final drive <fd>
veh_front_pwr_a
power available
force ava il. to P roduct2
to front axle
front whe e ls
De mux veh_rear_pwr_a
sum of
availa ble power available
P roduct
Mux tra ctive force s (N) De mux at rear axle
Force /Line a r S pe e d
tau/(tau-1)
1/(tau-1)
Sum3 pgs_trq_out_a
driv ing shaf t to generator speed
To Workspace8
Abs
From <v eh> (rad/s)
Gain
Wheel speed at pgs_spd_out_a
beginning of time
To Workspace9
step
pgs_pwr_out_a
To Workspace11
Product7
Mux 2
22 kW generator/controller scaled
from Mannesmann Sachs 63 kW
PM generator/controller
Electrical-machines
efficiency maps
0.85
0.7
0.6 0.75
0.8
0.65
0.85
0.9
0.75
0.8 0.75
0.85 0.9
0.7
0.65
0.6
0.75 0.85
0.75 0.55
0.8
0.7
0.7
0.8
0.85 0.65 0.85
0.6 0.6
Electrical-machines
0.75
0.65
0.55 0.55
0.65
0.75
0.55
ICE Operation - Fiat - Alfa Romeo 2.4 JTD (100 kW) Turbo Diesel Engine
90
30.1
80 37.1
70 34.1
60
50
Torque (Nm)
34.1
40
36.1
30.1
30
20
22.1
26.1
10 max torque curve
design curve
output shaft
op. pts(includes inertia & accessories)
0
0 500 1000 1500 2000 2500 3000 3500 4000 4500
Speed (rpm)
0.72
82
0.7
80
0.68
MilesPG
0.66
78
NOx
0.64
76
0.62
0.6
74
0.58
72
-5 -4 -3 -2 -1 0.56
-5 -4 -3 -2 -1
tau
tau
73
72 0.1972
71
70 0.197
MPGGE
CO
69
68 0.1968
67
66 0.1966
65
64 0.1964
-5 -4 -3 -2 -1 -5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1
tau tau
HC (grams per mile)
PM (grams per mile)
0.018
0.0196
0.0194
0.0179
0.0192
PM
0.019
HC
0.0179
0.0188
0.0178
0.0186
0.0184
-5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 0.0178
tau -5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1
tau
Accel_0_60 (seconds)
12 Grade_55mph
8
11.5 7.5
11 7
Accel_0_60
Grade_55mph
6.5
10.5
6
10
5.5
9.5
5
9
-5 -4 -3 -2 -1 4.5
-5 -4 -3 -2 -1
tau
tau
Kjara
• Two-seat roadster
• The engine drives
the rear wheels
Leone Martellucci
University of Rome I "La Sapienza" - Italy
RA RB
O RA
is proposed in this paper. RB
O1
B
B R2
B B
O2
R2
O1
R1 R1
O2
τ
Figure 3: PGS links in the SIPRE 3
Tengine = ⋅ Tgenerator
1- τ
From the Willis formula the relationship among the (6)
various speeds can be conveniently derived as follows Tdriving_shaft 1
=-
1 τ −1 Tengine τ
ω engine = ⋅ ω driving _ shaft + ⋅ ω generator
τ τ
The dynamic problem is completely settled when just
one of the three axis torques is fixed. For an automotive
ω driving_shaft = τ ⋅ ω engine + (1 - τ ) ⋅ ω generator (3) power unit the main condition can be the driving torque
Tdriving_shaft needed for the traction, so Tgenerator and Tengine
have to exactly match the Equations 6 and the choice of
the value of the transmission ratio τ affects the kinematic
τ 1
ω generator = ⋅ ω engine − ⋅ ω driving _ shaft and dynamic behaviour of this hybrid traction system.
τ −1 τ −1
It is important to note that the vehicle's dynamic
behaviour is controlled via the brake and accelerator
pedals, acting exclusively on the electrical traction
From these relations is clear that, unlike a traditional equipment arranged on the rear axle. From this point of
generator set used in series hybrid systems, the view this vehicle is managed just like a full electric
electrical generator of SIPRE 3 works at variable speed vehicle. This is possible considering that the hybrid
equipment acting on the front wheels, which supplies an PDS = Tdriving_shaft ⋅ ω driving _ shaft =
additional traction torque to the vehicle and cooperates
with the rear electrical traction, is characterised by a fully Tdriving_shaft ω driving_shaft
automatic operation without the need of an external Tengine ⋅ ω engine ⋅ ⋅ =
management. Furthermore, the engine needs no starter
Tengine ω engine
æ 1 ö æ ω driving_shaft
because the generator can act to crank it. Because ö
there is no reverse gear, vehicle reversing is a motor- P0 ⋅ ç − ÷ ⋅ ç ÷
only drive mode. è τ ø çè ω engine ÷
ø
Looking at the relations among the speeds in the PGS, it (8)
is interesting to observe that the speed of the generator PEG = T generator ⋅ ω generator =
is zero when
Tgenerator ω generator
ω driving _ shaft = τ ⋅ ω engine Tengine ⋅ ω engine ⋅ ⋅ =
ω generator = 0
(7)
Tengine ω engine
æ 1 - τ ö æç ω generator ö
÷
and the two degrees of freedom supported by the P0 ⋅ ç ÷⋅
kinematics of the PGS allow the generator to operate in è τ ø çè ω engine ÷
ø
two different quadrants, depending on the engine’s and
nd th
vehicle’s speeds, in generating (2 and 4 quadrants) or
st rd
motoring (1 and 3 quadrants) mode (Figure 4). The
generator/controller allows thus a bi-directional flow of where PDS is the power that reaches the driving shaft, P0
power from the batteries to the PGS unit. is the power supplied by the ICE and PEG is the power
absorbed by the electrical generator.
From the previous kinetic and dynamic equations it is Assuming the speed ωengine of the ICE fixed, the ICE
easy to calculate the expressions for the power power transmitted directly to the driving shaft increases
distributions over the PGS axes (Eq. 8): linearly with the vehicle speed, related to ωdriving_shaft. At low
speed the traction power is thus mostly supplied by the
electrical drive. The PGS acts like a power-split device
and a CVT (Continuously Variable Transmission) at the
same time.
PEG 1 ω driving_shaft
σ= =1- - ⋅ (10)
P0 τ ω engine
gal
ex_calc
t total fuel used (gal)
Clock To Workspace AND
fc_emis
emis
HC, CO,
vehicle controls NOx, PM (g/s)
sipre hybrid exhaust sys
<vc> control stategy <ex>
ICE fuel
<cs>
converter
<fc>
<sdo>
Product4
Ground
drive cycle gear generator/
<cyc>1 reduction controller <gc> power energy
bus <pb> storage <ess>
vehicle <veh>1 planetary
front wheel and front
gear set <pgs>
axle <wh> final drive <fd>
#
veh_front_pwr_a
power available
force avail. to Product2 to front axle
front wheels
veh_rear_pwr_a
sum of Demux
available power available
Product at rear axle
Mux tractive forces (N) Demux
Force/Linear Speed
When the vehicle is at rest (ωdriving_shaft=0), σ = 1 in Eq. 10 positive power (the electrical machine operates as
and the vehicle is a series hybrid because the whole motor) flows from the batteries to the PGS unit.
power generated by the ICE flows uniquely to the
generator. At higher speeds the powertrain tends to
behave like a parallel hybrid and increases the power SIMULATION RESULTS
flowing to the front wheels.
The vehicle used in the simulations was a PNGV-type
THE SIMULINK MODEL vehicle. Its components and characteristics included:
Figure 5 shows the Simulink model which allows • 28 kW CIDI engine (scaled from Fiat 100 kW 2.4 L
ADVISOR to generate energy flows, emissions and Turbo Diesel Common Rail Direct Injection engine)
consumption of the vehicle. The model describes the
front and rear axles of a vehicle provided with a SIPRE 3 • Auxilec Thomson 32 kW PM traction drive [6]
drivetrain in the 4WD layout and has been built
modifying ADVISOR at the block diagram level, by • 22 kW generator/controller (scaled from
programming in Simulink. Although ADVISOR relies Mannesmann Sachs 63 kW permanent magnet
heavily on a backward-facing approach for its operation generator/controller [7])
[5], it is possible to reconnect ADVISOR block diagrams
to model new vehicle types and add new component • Twenty 18 Ah lead acid batteries (Optima spiral-
block diagrams, like the PGS one, which are related to a wound VRLA, data from NREL tests)
forward-facing approach. The generator/controller block
diagram allows bi-directional power flow, depending on • Mass of 1010 kg.
the sign of power and following the convention in Figure
4: negative power (the electrical machine operates as • Coefficient of drag (Cd) of 0.3
generator) flows from the PGS unit to the batteries;
•
2
Frontal area of 2 m .
MilesPG NOx (grams per mile)
84 0.74
0.72
82
0.7
80 0.68
MilesPG
0.66
NOx
78
0.64
76 0.62
0.6
74
0.58
72 0.56
-5 -4 -3 -2 -1 -5 -4 -3 -2 -1
tau tau
73
72 0.1972
71
70 0.197
MPGGE
CO
69
68 0.1968
67
66 0.1966
65
64 0.1964
-5 -4 -3 -2 -1 -5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1
tau tau
0.0194
0.0179
0.0192
PM
HC
0.019 0.0179
0.0188
0.0178
0.0186
0.0184
-5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 0.0178
tau -5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1
tau
Accel_0_60 (seconds) MPGGE
12 77
76
11.5
75
11
Accel_0_60
74
MPGGE
10.5
73
10
72
9.5 71
9 70
-5 -4 -3 -2 -1 0.2 0.4 0.6 0.8 1
tau tau
0.62
7.5
0.61
7
0.6
Grade_55mph
6.5 0.59
NOx
6 0.58
0.57
5.5
0.56
5
0.55
4.5 0.54
-5 -4 -3 -2 -1 0.2 0.4 0.6 0.8 1
tau tau
0.211
0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
tau
HC (grams per mile) aim of that project was to develop a laboratory vehicle
0.02
useful to experiment various hybrid configurations.
0.0199
The idea is to use the current layout, but adapted to a
SIPRE 3 system, maintaining a big part of the system.
0.0199 Kjara is a roadster and the powertrain characteristics,
the arrangement of mechanical components, the
0.0199 architecture of chassis and suspensions are typical of a
sports-car (see Figg.19-20).
0.0198
HC
0.0198
0.0197
0.0197
0.0196
0.0196
0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
tau
0.0194
PM
0.0193
0.0192
0.0191
0.019
0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
tau Figure 20: Kjara vehicle prototype on Lombardore
(Turin) track.
Figg. 14-18: Parametric study. Performances over The powertrain is composed by a thermal engine and an
electric motor: the turbocharged diesel engine is
the US FTP cycle when τ changes
installed in the rear, mid-mounted, transversally. It drives
between 0.3 and 0.9.
the rear wheels through a conventional 5-speed
gearbox. The engine is a 85 kW Fiat 2.5 L turbocharged
indirect injection diesel engine. The 18.3 kW electric
motor is installed in the front, mid-mounted,
DESIGN APPROACH longitudinally, and drives the front wheels by means of a
transaxle 2-speed gearbox with differential. Its control
In order to give a detailed vision of the hybrid propulsion electronics, manufactured by Celco Profil, is put on the
system here discussed, it is useful to present the design dashboard by the right side of the steering wheel. The
approach for the realization of a working prototype. lead-acid battery pack is composed of 24 cells, each
The Department of Nuclear Engineering and Energy weighting 9.4 kg (total battery weight is 225.6 kg). Being
Conversions of University of Rome “La Sapienza” has the cells sealed, they can be put flatwise on one side
recently realized a hybrid sports-car prototype, called under the two seats.
Kjara, in collaboration with the National Research With such kind of powertrain, it is easy to design a quick
Council (CNR); Kjara [9] is a parallel hybrid car and the conversion of the system to a SIPRE 3 configuration; in
fact, with the simple and quick substitution of the Diesel The design approach here proposed can be realized in
engine currently mounted on Kjara prototype with the short time and at an affordable cost (for a University it is
planetary gear set group (see fig.21: PGS + ICE + EG + very important !) using an existent PGS. The drawing in
final drive, but on rear wheels), it is possible to switch to Fig. 22 represents in fact a feasible solution that uses a
a SIPRE 3 propulsion system. PGS with incorporated automotive differential, all
derived from 1992 4x4 FORD ESCORT. This
component, light and with a good transmission ratio,
matches well the SIPRE 3 demands, and benefits
obviously of a good reliability and quality. The adoption
of this component is a very good opportunity for a
prototype, cutting consistently the design and
construction costs.
In figure 21 it is possible to see the design approach • vehicle reversing is a motor-only drive mode
here proposed.
Maintaining the electric motor, its electronics and the • the generator cranks the engine
battery pack, and mounting the PGS set with a little
diesel engine and the electrical generator, it is possible • the vehicle operates electrically when driving into
in few hours of mechanical work to pass from the restricted zones
original parallel configuration to the SIPRE 3 system.
• 2WD and 4WD layouts are available
DEFINITIONS, ACRONYMS, ABBREVIATIONS
• the generator can spin backwards
ADVISOR:
ACKNOWLEDGMENTS ADvanced VehIcle SimulatOR
Graham T. Reader
University of Windsor, Canada
http://www.ott.doe.gov/hev/gmaccomp.html
http://www.ott.doe.gov/hev/gmaccomp.html
Cold Space
Power Piston
P T
3
QRi Qi Qi
4 3 4
2
QRo QRi QRo
2
1
Qo 1
Qo S
V
ADVISOR Simulation of SE HEV
1 2 3 4
P
T
3
Qi
QRi 4
Qi
2 3 4
QRo QRi
1
QRo
2
1
Qo V
ADVISOR Simulation of SE HEV Qo S
P 2
3
4
1 T
3
Qi
QRi 4
Qi
2 3 4
QRo QRi
1
QRo
2
1
Qo V
ADVISOR Simulation of SE HEV Qo S
• ICE • SE
Energy Input
Energy Input 100%
100%
Preheater
Regenerator 43%
430%
Thermal Loss
3%
Exhaust gases
Exhaust 14%
gases
28%
Mechanical
Radiation Loss
Loss 4%
7%
Cooling
Cooling System
System 45% Work Output
33%
Work 34%
Output
32%
Heat
Transfer
Parameters
SE
Simulation
MARWEISS
Map Generation
(Fuel)
30
25
20
Efficiency [% ]
15
Efficiency Map
10
46.2919
5 42.0836
37.8752
0 GM GPU-3 33.6669
1 1001 2001 3001 4001 5001 6001 7001 35-40
29.4585
6000 Speed [rpm] 30-35
25.2501
Torque (Nm) 25-30
EFFICIENCY 1.1 [Mpa] EFFICIENCY 2.1 [Mpa] 21.0418
5000 EFFICIENCY 3.1 [Mpa] EFFICIENCY 4.1 [Mpa] 20-25
16.8334
15-20
12.6251
4000 10-15
Power [Watts]
8.4167
5-10
3000 4.2084 0-5
0
100
1975
3850
5725
7600
2000
0
1
01
01
01
01
01
01
01
01
01
01
01
01
01
1
50
10
15
20
25
30
35
40
45
50
55
60
65
70
Speed [rpm]
A Ambient Air
E E W A A F Working Fluid
Ta Tcl Ta
Ambient Coolant Ambient Conduction
W W
Convection
Trd Tht Radiation
Radiator heating
Prius
Parallel
Conventional
Stirling
Stirling
Parallel ICE
Series ICE
Series
ADVISOR Simulation of SE HEV
30.00% Overall Efficiency
Overall Efficiency (PE)
25.00%
20.00%
15.00%
10.00%
5.00%
0.00%
Prius
Stirling
Stirling
Fuel Cell
Parallel
Conventional
Parallel ICE
Series ICE
Series
ADVISOR Simulation of SE HEV
a
250% 226%
204% 198%
200%
150% 124% 128%
100%
100%
50%
0%
Conventional
Conventional
Conventional
vehicle ZEV
HEV series
HEV series
with ICE
with SE
Electric
with ICE
with ICE
EZEV
with SE
EZEV
Graham T. Reader
University of Windsor, Canada
Overview of the SE
The inherent advantages of the SE over the ICE are Cooling System
mainly due to the engine’s continuous external heating 33%
[10]. For comparison purposes the typical energy flow -
diagrams of the SE and the ICE are presented below: Work Output
32%
The SE represents a potential alternative to ICE in 4. Third Order Design Methods, or nodal analysis
automobile applications. The main advantages are the: methods simulate the engine by dividing it in a given
§ Ability to use any heating source number of section employing nodes, and solving
§ Highest possible theoretical thermal efficiency numerically the differential equations of
§ Highest possible theoretical mechanical efficiency Conservation of Energy, Mass and Momentum, for
§ Quieter operation than ICE and gas turbines each of the nodes.
§ Lower vibration than ICE
§ Lower emission of pollutants when using fossil fuels, The third order methods are the most accurate and
because of their continuous combustion process appropriate for modeling specific SEs, however, they
§ Cogeneration capability: Heating interior of vehicle require a longer computational and design time than
utilizing the energy kept in the engine’s refrigeration other methods.
system
§ Regenerative braking capability. The engine MARWAISS, a second order method is used in this
operation can be reversed, transforming mechanical project because of its capability to analyze different
power at the output shaft into heat, for deceleration. types of SE within a reasonable computational time. It is
Later on, this heat can be transformed back into validated as a useful computer simulation tool with
shaft power for accelerating the vehicle results close to test bed information, allowing modest
• Higher efficiency than ICEs confidence in the simulation and analysis of energy flow
• Long operating life and efficiency.
• Extended periods between maintenance
• Low lubricating oil consumption MARWEISS or Martini-Weiss is a second order SE
• Reliability (continuous working operation) design aid developed by W. Martini and implemented on
• Higher power to weight ratio when operating at high a PC by M. Weiss at the University of Calgary, Canada.
pressure Based on the Martini isothermal model, performs
decoupled corrections for determining energy losses. It
The main areas of improvement are the: is written in FORTRAN 77, and includes input screens
§ High manufacturing costs with the present low level for the engine parameters. The results of the simulation
of production can be visualized in graphical format and exported to a
§ High materials costs text file. [3]
§ size and weight, somewhat larger than ICE
§ Larger radiator needed compared to ICE, with MARWEISS utilizes the engine general dimensions, heat
closed system cooling exchanger component sizes, drive system dimensions
§ Slower response time compared to Otto engines, and miscellaneous parameters for performing the
similar response as Diesel simulation. It allows the user to redefine a variety of
§ Performance depends on high quality manufacture dimensions, engine configurations, control strategies
§ Failures are catastrophic due to high pressurization and operational parameters.
on engine working fluid
SE Vehicle applications
Simulation of SE
The main thrusts of the modern SE development have
The methods for analysis of SE are classified according been automobile and cryocooling applications. The
to the assumptions and the mathematical theory that history of the modern SE started in the late 1930’s, when
encompass them. They are classified in four main types: Philips researched the development of a 1kW-power
generator. In 1948 Philips designed a 149 kW engine,
1. Zeroth Order Design Methods: simple models based which was seen by Henry Ford II and other major
on an idealized mathematical model (Schmidt cycle automakers in Eindhoven, motivating an interest on the
analysis) and on engine experimental constants engine. Philips used hydrogen and helium for high-
adequate for an initial approximation. speed and high power density SE.
The methods to analyze HEV employ a static or dynamic • PHIC : HEV ICE price, 204% conventional ICE.
approach for vehicle control to optimize several • PICE : ICE price, 10% price of conventional ICE
parameters, such as efficiency, weight or emissions. • PSE : SE price, for a production level of 10,000 units,
would be around 280% the price of an ICE,
The simultaneous optimization of fuel efficiency and according to studies by Jet Propulsion Laboratory
pollution emissions requires of strong mathematical [22].
algorithms in the computer simulations programs. • PEZEV : EZEV ICE equipment price, given the low
ADVISOR performs forward/ backward analysis to emission of pollutants inherent to the SE, there
review the validity of the information. Other methods would be a reduction of 24% basic vehicle price
currently in study, involve computer models with genetic associated to EZEV ICE equipment.
algorithms that combine billions of options on the
different parameters that determine the engine Then:
performance, for finding better results [19]
The ADVISOR program, launched in 1994 by NREL, PHS = PHIC − PICE + PSE − PEZEV
considers dynamic vehicle operation conditions, for Eqn. 1
optimizing both energy efficiency and emissions, in
conventional and advanced vehicles [20]. It is based on PHS = 204 − 10 + 28 − 24 = 198%
Simulink block diagrams and supported by Matlab data In conclusion, the cost of a HEV with SE would be very
files that contain vehicle configuration, control, and similar to the cost of the HEV ICE that complies with
performance data. The use of Simulink provides EZEV legislation. Environmental concerns themselves
accessible documented code in a graphical environment, are not a reason for the general public to pay double for
which makes modification relatively easy, helping an alternative vehicle.
ADVISOR to be well suited to collaboration between
researchers and for distribution to the public. SE and HEV share a common future: In order to be
commercialized successfully, they must first bring
The code incorporates several drive train types including answers to common problems:
electric, conventional ICE (spark ignited or compression • Need to have adequate production levels to reach
ignited) and HEV in series, parallel and split mode with significant cost reductions in tools, materials and
ICE or fuel cell as the APU. One of the objectives of this fixed capital equipment costs, independently of the
project is to include SE to the drive trains available. APU
• Need to internalize environmental costs in
ADVISOR allows the user to interchange a variety of
transportation pricing [23] [24]
components, vehicle configurations, and control
strategies. Modification of data files to represent new or
The implementation of pollution penalty measures (e.g. SE selection
increased fuel taxes, licensing or driving fees) applied to The engine to be simulated can be selected from
vehicle use and ownership, can motivate the public different sources:
interest towards environmentally friendly efficient a) Engine previously studied with ADVISOR,
vehicles. This will increase the research and production selecting the engine in the graphical interface of
of such vehicles, and a subsequent price reduction. ADVISOR
b) Design of a new engine with MARWEISS. Define
SOLUTION DESIGN a new engine based on i) MARWEISS' engine
parameter definition or ii) edition of an existing
GENERAL SIMULATION MODEL engine
c) Engine previously simulated with MARWEISS,
This section describes the steps considered for selecting the text file with the SE simulation
integrating the simulation tools. results (MARWEISS outputs for the desired
engine)
In order to predict the performance of a variety of SE, d) Input experimental results of the SE in a text file,
the simulation model is flexible allowing to model HEV or for selecting this file later in ADVISOR, as
conventional vehicles, in any configuration. described in option a)
Engine
In ICEs, the power output is a function of the mass
Analyze flows of air and fuel. Because of SE nature, the
Analyze MARWEISS Experimental
Advisor model New Engine
model Data
power output can be controlled in different ways. The
most common ways are:[25]
Control Strategy
1) Pressure level of working fluid
2) Stroke variation
Operation 3) Phase angle variation
Parameters
4) Maximum compression / Expansion volume
5) Temperature of Hot / Cold spaces
Heat Transfer
Parameters The control variable determines the type of power
maps that are required for the simulation. Most
commonly, the power maps will be given as a
Stirling Engine function of:
Simulation
i) pressure level and the angular speed, or
ii) phase angle and the angular speed.
Emissions
Simulation Operation Parameters
Iterative
The power maps generated in MARWEISS tabulate the Calculation
output power and the efficiency as functions of the
engine speed, given different pressures of the working
Display Operating
gas, in the case of engines controlled by working gas Conditions
pressure. For phase angle controlled engines, the power
maps tabulate the output power and the efficiency as
functions of the engine speed, given different phase Power Map
angles. Generation
Stirling Engine
File
Reactants Products This module has as function assembly the maps for fuel
consumption, torque, speed and emissions; writing them
in a text file that ADVISOR can utilize as input file. Also
writes the basic identification of the engine. The output
file of this step can be edited for reflecting any other
required change on the engine.
Burner
1
fc = × 3.6 × 10 6 Eqn. 13 Q x 2 y , r = ε * σ * A * (Tx4 − Th4 ) Eqn. 15
η hPR
For convection:
Qin. x − Qout . x
t
Tx = ò Eqn. 18
o m *C
x p. x
Preheater Burner Hot Space Regenerator Cold Space Engine Block Hood
W Cooling Water
E E F F F A
C Combustion reactants
C Tb Trg
Tp Ths Tcs Te Thd E Exhaust mixture
A Ambient Air
E E W A A F Working Fluid
Ta Tcl Ta
Ambient W Coolant Ambient Conduction
W
Convection
Trd Tht Radiation
Radiator heating
RESULTS
The simplicity of ADVISOR for simulating vehicles with The ICE employed for the series and parallel
different configurations and components is an advantage configurations, was a Geo 41 kW engine. For the split
over other HEV simulation programs. Once the model of configuration the ICE is a 1.5L Prius (Atkinson cycle)
a component of ADVISOR has been validated, it is engine.
simple and fast to perform simulations.
The fuel cell is an IFC 50kW net hydrogen fuel cell stack.
As an example for evaluation of the modules developed,
it was performed a cycle of simulation routines with The SE employed is a GPU-3 6kW engine. This is a
different APUs in series and parallel configurations, with smaller engine compared to other engines. It was
a basic automobile arrangement. The following table employed in the Sir-Lec I from GM, the first SE HEV [3].
shows the comparison on energy consumption for It was selected for comparison of simulation results with
performing the same test routine, two cycles of the published data.
Federal Urban Driving Schedule.
Some results from the table, are the percentage of
The types of APU and configurations considered are: energy losses on the engine in descending order:
Conventional vehicle (74%), series ICE (69%), parallel
1. HEV ICE in series configuration ICE (69%), split ICE (66%), fuel cell (63%), series SE
2. HEV ICE in parallel configuration (60%) and parallel SE (54%).
3. Conventional ICE vehicle
4. HEV FC vehicle The overall efficiency, including potential energy:
5. HEV ICE in split configuration Conventional vehicle (17%), series ICE (19%), parallel
6. HEV SE in series configuration ICE (20%), fuel cell (22%), series SE (22%), split ICE
7. HEV SE in parallel configuration (23%), and parallel SE (27%).
Table 3 shows the energy balance after each of the The total energy used for performing the cycle: series
components of the vehicle and the percentage of energy ICE (27816 kJ), fuel cell (25666 kJ), series SE
loss in each component, as a fraction of the total energy. (24721 kJ), parallel ICE (22744 kJ), and parallel SE
(22744 kJ), split ICE (21110 kJ), conventional vehicle
(16825 kJ).
The results recall the importance of having a lighter
vehicle to save energy; Even though the conventional
vehicle had the lower efficiency of all the vehicles, it also
used the least energy of all, because it was also the
lighter vehicle.
Energy Usage [kJ] Series ICE Parallel ICE Conventional Fuel Cell Prius Series Stirling Parallel Stirling
Balance [kJ] Use [%] Balance [kJ] Use [%] Balance [kJ] Use [%] Balance [kJ] Use [%] Balance [kJ] Use [%] Balance [kJ] Use [%] Balance [kJ] Use [%]
Fuel Converter 27816 69.2% 22598 69.2% 16825 74.2% 25666 63.3% 21110 66.1% 24721 59.8% 22744 53.8%
Energy Storage 22967 0.2% 22348 0.2% 16252 24375 1.8% 20190 0.9% 22751 3.2% 21701 2.0%
Motor / Controller (Power) 19692 2.5% 20498 2.5% 16252 20957 4.9% 18874 2.1% 18775 5.2% 18889 5.7%
Gearbox (Power) 18662 1.3% 19561 1.3% 15430 1.3% 19892 1.5% 18874 17771 1.6% 17895 2.0%
Wheel / Axle (Power) 17852 1.1% 18747 1.1% 14759 1.0% 19036 1.2% 18067 1.3% 17094 1.1% 17127 1.6%
Motor / Controller (Regen) 17481 1.6% 17549 1.6% 14270 18651 0.6% 17716 0.6% 16733 0.6% 16361 1.6%
Gearbox (Regen) 17412 0.1% 17502 0.1% 14236 0.1% 18582 0.1% 17716 16664 0.1% 16321 0.1%
Aux Loads 15522 2.6% 15612 2.6% 12368 2.9% 16686 2.7% 15827 3.1% 14778 3.1% 14437 3.9%
Braking 14410 1.5% 14492 1.5% 11133 1.9% 15473 1.7% 14545 2.1% 13739 1.7% 13419 2.1%
Aerodynamic loss 12298 2.9% 12380 2.9% 9023 3.2% 13361 3.0% 12894 2.7% 11627 3.4% 11307 4.3%
Rolling Loss 9459 3.9% 9522 3.9% 6940 3.2% 10276 4.4% 9998 4.7% 8966 4.3% 8697 5.3%
Potential Energy 0 13.0% 0 13.0% -1 10.6% -1 0 10 16.1% 8 14.6% 0 17.7%
Powertrain:
Charge Sustaining Thermostat Series
Motors: 1 90 kW DC electric motors at rear axle
Engine: VW 1.9 liter TDI
Motor/Generator: 1 90 kW motor/generator
Energy Storage System (ESS): NiMH batteries and ultracapacitors
VW TDI specifications:
-Four-cylinder in-line diesel with direct
injection turbocharged (TDI)
-60 kW at 3,300 rpm
-Engine management system with EDC
-Sealed pressurized cooling system with
oil cooler in coolant circuit
-Height: 637 mm, Width: 710 mm,
Length: 545 mm
-Dry weight = 135 kg
Model MDF-375/75-A1-BS
NGM's 90 kW motor characteristic MDF-375/75-A1-BS specifications:
-Design Topology: Axial Flux, Dual Stator Fixed Air gap
800
-Overall size: 445 mm diameter x 189 mm length
700
-Weight: 84 kg
600
-Cooling: Forced Air
Torque (Nm)
500
-Mounting: Face flange per outline
400
300
-Max rated current : 320 amps
200
-Rated DC Bus voltage: 312 V
100 -Mechanical speed limit: 3,000 rpm
0 -Rotor inertia: 0.17 kg-m2
0 500 1000 1500 2000 2500 3000
Motor speed (rpm)
Ni-MH & Ultracapacitors
Pow er_Ultracap
80
Pow er_required
60
40
20
0
0 5 10 15 20
time (s)
Simulation Results
HWFET Cycle
120 1.02
100 1
0.98
80
ESS_SOC
0.96
Speed (kph)
60 0.94
40 0.92
0.9
20
0.88
0 0.86
-20 0 200 400 600 800 1000 0 200 400 600 800 1000
time (s) time(s)
60
50
Po
wer 40
(kW
) 30 Prqd
20 Pregen
10
0
0 200 400 600 800 1000
-10
time (s)
Simulation Results
FUDS Cycle
100 1.02
90
1
80
70 0.98
Speed (kph)
ESS_SOC
60 0.96
50
40 0.94
30 0.92
20
0.9
10
0 0.88
0 500 1000 1500 0 500 1000 1500
time (s) time (s)
70
60 Pregen
Prqd
50
40
Power (kW)
30
20
10
0
0 200 400 600 800 1000 1200 1400 1600
-10
Time (s)
Summary
Pow er_Ultracap
80
Pow er_required
60
40
20
0
Figure 2 - Characteristics of a 1.9 liter VW CIDI 0 5 10 15 20
engine [8] time (s)
ENERGY STORAGE DEVICE SELECTION Figure 3. Power curves vs time during hard
acceleration
The GW team has decided to use a hybrid
configuration of Ovonic NiMH batteries and Pinnacle Table 3 - NiMH battery stack specification
Research Institute, Inc. ultracapacitors to achieve the Battery pack specification Ovonic NiMH
power and energy storage requirements of the vehicle. [ 30 kWh]
During a short period of high power demand,
ultracapacitors will provide the necessary power to the No. of modules 25
motors. We project that the NiMH battery stack alone
will not be able to provide the high power required for Mass (kg) 445
hard acceleration (see Table 2) and would exceed the Volume (liters) 188.3
targeted mass limit. This has dictated the need of an
ultracapacitor pack which has 90 times the power Power@20% SOC (kW) 97.8
density of the batteries [9]. Other advantages of using
Voltage (V) 330
Table 4 - Ultracapacitor stack specification
ELECTRIC DRIVETRAIN
Ultracapacitor specification Pinnacle RIT Inc.
[ 130 kW] The drivetrain will consist of two prototypes
90 kW New Generation Motors Corporation (NGM)
No. of modules 3 DC motors. Each one of the motors will drive the front
drive shaft and the rear drive shaft. During low to
Mass (kg) 1.3 average power demand, only the rear motor is driving
Volume (liters) 1 the vehicle and the front motor is in the generator mode.
Specific power (kW/kg) 100 When the driver demands hard acceleration,
the EM clutch system will disengage the generator from
Voltage (V) 300 the engine and engage it to the frontal drive shaft. The
power circuit is now reversed and has turned the
generator into electric motor mode and will drive the
frontal drive shaft.
CONTROL SYSTEM DESIGN AND ITS
COMPONENTS
NGM's 90 kW motor characteristic
The control system will consiste of power
controllers for the electric motors and the on-board
800
computer system also called a Powertrain Control
Module (PCM). The new NGM-EV-C200 controllers 700
will be used to control and power the 3-phase DC 600
motors. The controllers have programmable logic
Torque (Nm)
500
capability allowing them to be optimized and matched 400
with the motor characteristics. The overall system
300
controller will be implemented on the on-board PCM.
200
This computer will be designed to monitor the state of
charge of the batteries and ultracapacitors, individual 100
motor voltage and current, vehicle speed, temperature 0
of various components, generator output, and will have 0 500 1000 1500 2000 2500 3000
various control functions including activation of the Motor speed (rpm)
CIDI engine and generator, EM clutch engagement or
disengagement, cooling fans, power limited cruise
control, and motor controller. Figure 4 - NGM’s 90 kW DC motors characteristic
The control strategy utilized is the thermostat The prototype NGM 90 kW DC motor is
series, which has the objective of operating mainly on chosen because of its high efficiency (>90%), high
electrical energy. The battery will be discharged until a specific power (each weigh only 55 kg) and capability
low state of charge (SOC) is reached. The engine is to operate efficiently in regenerative mode. Figure 3
then turned on to recharge the battery until a desired shows characteristic of the motor.
SOC is achieved. Advantages of this type of control
strategy are that the engine can be set to operate at one
point of torque and speed at which it is most efficient FUEL SYSTEM DESIGN AND MATERIALS
and least polluting. This prevents the engine from SELECTION
experiencing any transient load, which also should The proposed vehicle will utilize the existing
reduce the emissions. In addition, the operating point fuel system with only minor modifications to
can be set for the optimum performance point of the accommodate the new Volkswagen 1.9 liter CIDI
electric generator by using a simple gear coupling. engine. This has several advantages in that the fuel tank
Another potential reduction to emissions can come from is optimally placed for crashworthiness, the fuel lines
the addition of an electrically heated catalyst, which can and connections are commercially produced and meet
be easily implemented to reduce cold-start emissions all federal standards, and the integration into vehicle is
since there is ample knowledge of when the engine will already achieved. This minimizes the need to make
come on [10]. modifications in the layout and focuses the integration
effort on new components. The only difference in the
fuel system is that the flow rate is no longer regulated
by the gas pedal and fuel pump, but rather will maintain
a constant flow rate to the generator engine and will be
activated by a signal from the PCM. The existing fuel
tank can also be replaced with a smaller tank since the
conversion vehicle will require less fuel to achieve the Assuming a 10% decrease in aerodynamic
same driving range as the original vehicle. drag, an improved rolling resistance of 0.012,
Exhaust Emissions Control Design and Component maintaning vehicle mass of 2,000 kg (model year 1999
Selection Chrysler’s Jeep Cherokee [13]) and using a series
hybrid configuration, the results show that the vehicle is
The catalytic converter, muffler, and the capable of achieving more than 39 mpg fuel economy
tailpipe will have to be resized to accommodate the new (35.5 gasoline equivalent). This fuel economy
CIDI engine. The design will address emission control projection is more than two times that of a typical SUV.
to meet the Federal Tier 2 as required by the With the addition of efficient regenerative braking and
competition rules. Since the engine will operate at a control system design algorithm, and mass reduction,
constant speed for charging purposes only, high this design vehicle has a potential of reaching the 3
emissions typical of today’s vehicle due to acceleration, times current fuel economy. Table 5 shows an example
deceleration, and time at idle are virtually non-existent. of simulation parameters and conditions used to predict
vehicle performance and energy economy projection.
PERFORMANCE AND ENERGY ECONOMY Table 5- Simulation parameters
PROJECTION Simulation Set-up Parameters
Variables
The GWU team has used ADVISOR , a
Drive cycle HWFET (Highway Fuel
MATLAB-based dynamic simulation software to Economy Test Drive)
predict performance and energy economy projections of No. of cycles 25 (256 miles)
the proposed vehicle [11,12]. However, the software Initial SOC of battery pack 95 percent
could not fully simulate the proposed design of the
Maximum SOC 95 percent
hybrid energy storage system of NiMH batteries and
Minimum SOC 30 percent
ultracapacitors. The simulation model does not include Control Strategy Zero delta SOC Correct
a ultracapacitor stack as specified in the vehicle design. (SOCfinal = SOCinitial)
This would not affect the energy economy prediction
since the maximum power demand during the
simulation of FUDS and HWFET cycles is below 50
The use of two 90 kW motors and an
kW and it is achievable by the NiMH battery pack alone
ultracapacitors stack allows the vehicle to meet the high
(see Figures 5 and 6).
power demand during hard acceleration.
Power Requirement (FUDS cycle) Ultracapacitors although having high power density,
70
havw a limited energy storage capacity. Therefore, a
hybrid configuration of NiMH batteries and
60
ultracapacitors is proposed to meet the 300 mile range
50 requirement. The greatest challenges will be in design
40 and implementation of the EM clutch to use the front
Power (kW)
30
motor as the motor/generator and the coupling of the
ultracapacitors and battery stack.
20
10
0 SUMMARY
0 200 400 600 800 1000 1200 1400 1600
-10
Time (s) The GWU innovative series hybrid conversion
approach has a potential to meet and exceed the design
Figure 5- Power requirement vs time for 1 FUDS cycle goals outlined for FutureTruck 2000 competition. The
Power Requirement (HWFET cycle) use of highly efficient and compact DC motors as
generator and driving motors is unique and offers
60 benefits of cost saving and mass reduction. The use of
50 two transaxle motors also allows the vehicle to operate
40 in four wheel drive mode. The hybrid energy storage
system of NiMH batteries and ultracapacitors makes it
Power (kW)
30
possible for the vehicle to achieve high performance
20 and increase its range. A highly efficient CIDI engine
10 will provide power to the generator only in its high
0 efficiency region and low emissions. In addition, non-
0 200 400 600 800 structural body parts will be replaced with lighter
-10
composite materials and the vehicle aerodynamic drag
time (s)
will be reduced through several measures.
Figure 6- Power requirement vs time for 1 HWFET
cycle
CONTACT
12. Senger, R.D., Merkle, M.A. and Nelson, D.J.,
‘Validation of ADVISOR as a Simulation Tool for
Vahid Motevalli, Ph.D., P.E. Series Hybrid Electric Vehicle’. Society of
GW Transportation Research Institute, Automotive Engineers, 1998
The George Washington University,
13. DaimlerChrysler Corporation, ‘Jeep Grand
20101 Academic Way, Ashburn, VA 20147
Cherokee Service Manual’, 1998.
REFERENCES
Simulation of a Heavy-Duty
Series Hybrid Electric Transit Bus
Jean BAVARD / Atef GAYED
INTRODUCTION
Motor Motor
Control 80 kW 80 kW
Alternator Rectifier
supervisor
GOALS
l HYBRID ELECTRIC BUS
– Factory tests 2002
– Prototype validation 2003
– Mass production 2004
l Means :
– Energy Storage System
– Energy Management
– ZEV operating range
– Electric Traction
Time (Sec.)
ADVISOR_MODIFIED
Energy Storage 47157 35184 12027 0.75
Energy Stored -54
Motor/Controller 53934 45939 7995 0.85 24817 23008 1809 0.93
Gearbox 45939 44203 1735 0.96 24951 24817 134 0.99
Final Drive 44203 44203 0 1 24951 24951 0 1
Wheel/Axle 44203 42321 1882 0.96 25618 25279 339 0.99
Braking 328
Atef GAYED
Research & Development Engineer
In general, by pushing the key point "A" to the left Genset is the set of power/speed combinations for a
(minimizing the speed@maxpower), we can have better given vehicle where the maximum efficiency lies, where
bus performances and thus better energy recuperation. power and speed are limited by either the engine or the
generator.
ALSTOM second stage program will be with a
Permanent Magnet Synchronous Motor (PMSM) which The FC is a diesel engine, which have 93 kW maximum
will improve efficiency and regenerative power at high power and 4000 rpm maximum speed, with 38%
speed range. efficiency at 290Nm@2200rpm
The wheel drive line brake fraction (Wh_dl_brake_frac) Figure 4 presents the FUDS, TCL21E_AR urban drive
and the wheel friction brake fraction cycles. The last one represents a French drive cycle in
(Wh_fric_brake_frac) are shown in Figure 3. These Lyon City. These cycles have high dynamic variation.
assume that the maximum brake force will be applied to The maximum speeds are approximately 3 times the
the vehicle by the electrical motor in order to recover the average speed as shown in Table 2.
maximum admissible of the braking power. The
mechanical brakes will do the rest. Speed (Km/h)
FUDS
TCL21E_AR
Friction brake loss energy (kJ) 2242 2242 Motor_regen_energy_in (kJ) 24817 24817
Table 3: Regen simulation results for FUDS cycle for both motor regen_energy_out (kJ) 23008 23399
series control strategies with ADVISOR version
motor_out/veh_regen eff. 89.8% 91.3%
CASE 2 : FUDS WITH ADVISOR MODIFIED FILES
Friction brake energy (kJ) 328 328
The ALSTOM AC induction motor used in this study has
another speed/torque map, Figure 6, in the regenerative Friction brake loss eff . 85.37% 85.37%
operation mode. This allows it to absorbe 70kW
(dot line) more than in traction mode especially in the Fuel consumption (l/100km) 35 36.2
medium and high speed range. Introducing this
difference enabled us to approach the real motor
operation and refine the simulation study and lead to the Table 4: Regen simulation results for FUDS cycle for both
following results. series control strategies with ADVISOR modified version
Table 4 and Figure 7show that while braking : strategy with the electrical motor and the energy storage
system characteristics allow the recuperation of a very
L the mc_trq is not limited
large amount of the braking power.
L the wheel driveline regen power can accept up to
221kW First, the importance of the braking energy for EVs
L brake losses decreased by 85.37% for both control & HEVs was highlighted. Optimizing each individual
strategy. component of the system from one hand and recovering
the kinetic and potential energies on the other hand are
L the regen to motor has increased by 8.18%
the central goals of research for these structures. In the
L the motor energy sent to the battery pack has same time, the electrical motor and the energy storage
increased by 4.6% (ST) and 6.2% (SFO). system characteristics are of major importance to
L Slight improvement in fuel economy. The improve the regenerative recuperation. Consequently,
cs_parameters have to optimise in order to have a introducing the real motor speed/torque map in the
real gain for such point. regenerative operation mode have allowed us to
minimize the mechanical brake losses and thus improve
CASE 3: TCL21E_AR the regenerative energy efficiency as presented for the
urban FUDS and TCL21E_AR cycles.
Table 5 presents the simulation results for a thermostatic
control strategy in the same conditions and parameters ACKNOWLEDGMENTS
for the TCL21E_AR Lyon urban cycle. The first column
presents the original ADVISOR version while the second Special thanks to the ADVISOR team works specially
is with the modified one. We can see the same Valerie JOHNSON for their precious help.
conclusion as found for the FUDS cycles. The difference
is the amount of energy recovered. The gain in the CONTACT
mechanical brake is about 52% instead of 85.3%for the
FUDS cycle. This can be explained by the fact that the
Jean BAVARD was born in Dijon, France in 1958. He
TCL21E_AR cycle has half maximum and average
received the B.Sc. degree in 1980 in general and
speed than the FUDS.
electrical engineering from the Ecole Centrale de Paris,
CYCLES TCL21E_AR TCL21E_AR France. Since 1982 he joined the electrical traction
department of Alstom-transport in Lyon and participated
ADVISOR Adv_Modified in the development of variable speed drives with
alternative machines applied to railways and marines.
max decel (m/ss) 1.415 1.4
Currently, he is in charge of traction drives applied to
max_veh_brake_pwr (kW) 236 236 hybrid electrical bus. The main fields of interest are
system control, energy management and power
max_weh_brake_pwr (kW) 151 222.75 electronics.
veh energy (kJ) 41574 41574 Atef GAYED was born in Ismaïlia, Egypt in 1964. He
received the B.Sc. degree in 1986 in power electrical
motor regen energy_in (kJ) 38989 40018 engineering from Zagazig University, Cairo, Egypt. He
received the Ph. D. degree in 1996 from Nantes
motor_in/veh_regen eff. % 93.78% 96.25% University, France in the application of sliding mode
control to variable speed drives. Since 1999, he
motor energy_out (kJ) 35763 36003
participates as engineer for Alstom-transport in the
motor_out /veh_regen eff. % 86% 86.6% railway and hybrid electrical bus fields. His research
interest includes electrical motor control system, energy
Friction brake loss energy (kJ) 2018 967 management and power electronics.
Table 5: Regen simulation results for TCL21E_AR cycle 1. Yimin Gao, Liping Chen and Mehrdad Ehsani,
'Investigation of the Effectiveness of Regenerative Braking
for series thermostat control strategies for both original & for EV and HEV', SAE 1999-01-2910
modified ADVISOR version
2. Mellor, P.H. 'High Efficiency Drivbe-trains for Electric and
CONCLUSION Hybrid Vehicles' IEE 1999.
Table 1: cs_parameters........................................................................................................................................................................3
Table 2: FUDS & TCL21E_AR cycle data .........................................................................................................................................3
Table 3: Regen simulation results for FUDS cycle for both series control strategies with ADVISOR version ..................................4
Table 4: Regen simulation results for FUDS cycle for both series control strategies with ADVISOR modified version...................4
Table 5: Regen simulation results for TCL21E_AR cycle for series thermostat control strategies for both original & modified
ADVISOR version ................................................................................................................................................................5
T (Nm) Torque
P (kW) Power
ST Series Thermostat
Dart
Vehicle
Passenger Operator
On-board
systems
Vehicle
Infrastructure
30
25
20
m ph
15
10
-5
0 500 1000 1500 2000 2500 3000 3500
S e c s s in c e s t a rt
Payload
Traffic
Driver
Seasonal
…
mph
ØFluid flywheel
ØPayload
ØIntermittent accessory load
ØRWD
ØOptimisation
15
Ø Processing of measurement data was performed in
10
Matlab and was therefore easily linked to ADVISOR.
-5
0 500 1000 1500 2000 2500 3000 3500
Secs since start
Initial vehicle models are derived from the built-in The payload was set by running the model with
ADVISOR sub-system models, suitably modified to increasing load until the best match with the measured
reflect the conventional and Newbus vehicle designs. speed was achieved. In reality the payload changes
As the development progresses it is anticipated that with time and at some time it would be interesting to look
more of the elements will use Newbus-specific models. at changes in this as part of the simulation as it has a
These will be particularly important for the Newbus major influence on fuel consumption and the
custom hybrid controller, which incorporates a number of performance required.
new features, and the battery model, which is to include
a battery damage model to assess in-service battery ACCESSORY LOAD
lifetime. Some particular modelling issues arose during
the work: It would be helpful to be able to model varying accessory
loads as these can vary with engine speed and/or
FLUID FLYWHEEL strategy. Without this the load tends to be overstated on
the average or understated at the peaks. It is planned to
The conventional bus uses a fluid flywheel rather than a capture real operational data in this area.
torque converter and this was modelled as a manual
gearbox with slightly higher numerical ratios in the lower
gears combined with reduced overall efficiency.
It would be helpful to put this in explicitly. Even though Initial results confirm the potential for a lower power
this is almost certainly immaterial here, it would be hybrid design that offers lower purchase and operating
comforting to know that everything was the correct way costs, and reduced emissions in comparison with higher
round. power hybrid designs. Ongoing work will develop this
approach and allow the optimum power level to be
FUEL CONVERTER 'SWEET SPOT' specified taking environmental and commercial
considerations into account.
Some effort was taken to match the generator and fuel
converter to bring the fuel converter operating point to The use of in-service dynamic measurements coupled
maximum efficiency at low operating speed. This may with ADVISOR simulation is essential to allow this
be possible in a more efficient way with more familiarity approach to be carried out with confidence.
with ADVISOR.
By actually capturing measurements from vehicles in
MOTOR/CONTROLLER SPEED OPTIMISATION service it was possible to validate the data on a
conventional vehicle before moving to a hybrid design.
The optimisation tool was found to be very useful in
setting overall reduction ratio. The Newbus reduced envelope hybrid shows a reduction
in fuel consumption of almost 40% over a conventional
INITIAL FINDINGS bus. This is achieved with battery cycling designed for
long service life.
Initial results from the first ADVISOR hybrid models are
shown in Figure 5. Significant fuel savings were
possible in comparison with a 'full-envelope' hybrid
design. As is generally recognised, it is necessary to ACKNOWLEDGMENTS
carefully match the individual elements of the hybrid
driveline to deliver the potential of each hybrid design. The authors would like to acknowledge the assistance
given with access to existing route data by Mike Weston
The initial simulations confirm that a reduced envelope of London Buses and Andy Eastlake of Millbrook Proving
vehicle could service all sections of the route measured Ground. We would also like to thank the ADVISOR
in London. In many cases the acceleration levels team for their advice and response to requests for
currently used offered no savings in route timing, but additions to ADVISOR, in particular the addition of
clearly reduced passenger comfort and increased fuel gradient input to the model.
usage and emissions. There is therefore scope for
further reductions by further reducing performance, CONTACT
though this clearly depends on the results of ongoing
route measurements. Mike Kellaway can be contacted at [email protected],
Alan Ponsford at [email protected]. The Newbus
NEXT STEPS Technology Limited website at www.newbus.com also
contains further material about our work on Hybrid
The data collection technique is to be used in a number Buses.
of UK cities to collect a representative sample of route
data. This will then be characterised and used as
simulation data for a new reduced power hybrid design.
The use of carefully chosen actual route data will allow DEFINITIONS, ACRONYMS, ABBREVIATIONS
detailed assessment of the operation of the new vehicle
and optimisation of its control strategy. GPS:
Mike Ogburn
Mechanical Engineering
Virginia Tech
Hybrid Electric Vehicle Team of Virginia Tech
•Zero emissions
capability
DC/DC Converter
Inverter
Air Compressor
Fuel Cell Controller
Fuel Purge
Battery Traction
Air Pack Motor/
Air Compressor
Exhaust Inverter
Boost Converter
Air
In Air HVAC
DC/DC
Comp Comp HVAC
Fuel Conv
Inv Inv Load
Air
Electrical
Hybrid Electric Vehicle Team of Virginia Tech
Observed
Voltage (VDC)
80
voltages
60
40
20
0
0 50 100 150 200 250 300 350
Current (A)
Hybrid Electric Vehicle Team of Virginia Tech
0.5
(Gross Stack Power)
0.4
0.3
0.2
(Net System Power) Stack Efficiency
0.1
System Efficiency
0
0 2500 5000 7500 10000 12500 15000
Power (Watts)
Hybrid Electric Vehicle Team of Virginia Tech
AD V I S OR
• Proven vehicle modeling platform for:
– Vehicle Energy Use and Emissions
– Powertrain Performance and Control
• Easy to use
– User-Friendly Graphical Interface
– Easily interchange and analyze:
• Components
• System Configurations
• Control Strategies
Hybrid Electric Vehicle Team of Virginia Tech
AD V I S OR
• Uses Matlab and Simulink
– Component building blocks
– Data files load component specifications
• Other Features Include:
– Parametric analysis
– Sensitivity studies
– Combined speed and grade traces
– It’s Free! www.ctts.nrel.gov/analysis
Hybrid Electric Vehicle Team of Virginia Tech
emis
HC, CO,
NOx, PM (g/s)
drive cycle
<cyc>
Version & vehicle <veh> wheel and final drive <fd> gearbox <gb> motor/
Copyright axle <wh> controller <mc> power energy
bus <pb> storage <ess>
Hybrid Electric Vehicle Team of Virginia Tech
6 0
City Driving Cycle
5 0
4 0
S P E E D (m ph)
3 0
2 0
1 0
0
0 2 0 0 4 0 0 6 0 0 8 0 0 1 0 0 0 1 2 0 0 1 4 0 0
T IM E (s e c )
35000
Power (W)
15000
-5000
0 200 400 600 800 1000 1200
Time (sec)
-25000
Motor Controller and Accessory Power (vehicle data and vehicle model)
Hybrid Electric Vehicle Team of Virginia Tech
10000
Data
Model
7500
Power (W)
5000
2500
0 Time (s)
2000
Dyno: 29.1
1500
1000 Battery
Fuel Cell
500
0
Model Dyno **SOC Corrected
Hybrid Electric Vehicle Team of Virginia Tech
60
50
40
D(mp)
h
30
SPEE
20
10
0
0 100 200 300 400 500 600 700 800
TIM E (s e c )
60000
40000
Power (W)
20000
0
0 100 200 300 400 500 600 700 800
-20000
Time (sec)
-40000
Motor Controller and Accessory Power (vehicle data and vehicle model)
Hybrid Electric Vehicle Team of Virginia Tech
8500
Data Model
7000
Power (W)
5500
4000
300 320 340 360 380 400
Time (s)
Hybrid Electric Vehicle Team of Virginia Tech
Time (s)
0
-1 0 200 400 600 800
-2
Amp hours
-3
-4
-5
-6
Data
-7
Model
-8
Hybrid Electric Vehicle Team of Virginia Tech
400
375
350
Voltage (V)
325
300
275
250
Dyno Batt Volts
225
Model Batt Volts
200
0 200 400 600 800
Time (s)
Hybrid Electric Vehicle Team of Virginia Tech
4000
Dyno: 26.7
3000
Battery
2000
Fuel Cell
1000
0
Model Dyno **SOC Corrected
Hybrid Electric Vehicle Team of Virginia Tech
• Accurate for....
– Trends and transients
– Energy generation and energy storage
– Identifying inefficient systems
• Lacking because
– Not easily adaptable
– No forward looking capabilites
– Inter-system impacts not considered
Hybrid Electric Vehicle Team of Virginia Tech
2000
Energy (kW-hour)
500 Battery
Fuel Cell
0
Model
Constant speed case VariableDyno
speed case
-500
**SOC SUSTAINING!
Hybrid Electric Vehicle Team of Virginia Tech
10000
8000
6000
4000
2000
0
0 200 400 Time
600(s) 800 1000 1200 14
Fuel Cell System Gross Power
70
60
50
Temp. (C)
40
30
20
10
0
0 200 400 600Time (sec) 800 1000 1200 1400
Fuel Cell Coolant Output Temperature
Hybrid Electric Vehicle Team of Virginia Tech
100
80
60
40
Air Temp (C)
20 Air Humidity (RH%)
Water Temp (C)
0
0 200 400 600 800 1000 1200 1400
Time (s)
Hybrid Electric Vehicle Team of Virginia Tech
www.futuretruck.org
Hybrid Electric Vehicle Team of Virginia Tech
• Fuel Infrastructure!!
• Design for Manufacturability / Cost
• Weight Reduction = Efficiency + Range
• System Performance
– Current customer expectations worldwide
– Durability over 100,000 mi or 5000 hr
• Education of the public:
If you build it, will they come?
Hybrid Electric Vehicle Team of Virginia Tech
Fuel
Heat Exchanger
Cell
(Recirc.)
Fuel
Cell
Engine Compartment
Hybrid Electric Vehicle Team of Virginia Tech
DC/DC Converter
Contactor Box
HVAC Compressor
(Dashed) Buck Converter Air Compressor
Inverter HX Air (for Air Compressor) Motor Controller
Compressor
Vacuum
Pump
Fuel Cell Cooling Deionized Power Steering
Water Pump Water Tank
Pump (Dashed)
Humidification
Front
Humidification
Recirc. Pump Pump
Note: Components denoted by a dashed line are
located above other components
ADVISOR and Solid Modeler Integration
Mark Henault, ESS West, Inc.
INTRODUCTION
Mark Henault
ESS West, Inc.
OVERVIEW
l Pro/HEV
– First integration of ADVISOR with Pro/Engineer
Optima (Pb-Acid)
60
50
40
30
20
10
0
0 1000 2000 3000 4000
Accessory Load (W)
l ICEM-CFD
– Generates mesh of
volume
– Automatically
launches thermal/fluid
solver
T=0 min
T=10 min
T=20 min
Summary
l Brief Overview of the benefits of ADVISOR &
Solid Modeling Integration
– Packaging studies
– Mass & CG calculations
– Interior volume generation for occupant comfort
analysis & a/c sizing
– Generation of geometry for future work (ADAMS,
Crash analysis, etc.)
◆ Export to ADAMS/Car
■ Method
◆ ADAMS/Car full vehicle model using customized powertrain
template
◆ Use the ADAMS/Controls interface for exporting the plant
◆ Modified ADVISOR model to work with ADAMS/Car model
ADAMS/Car
Full Vehicle Model
Mass, Inertia Handling Results,
& Geometry Animations
Co-simulation
ADVISOR
Hybrid powertrains,
Energy management,
Controls
■ Method
◆ Output mass and inertia properties from ADVISOR to ADAMS/Car
• Optionally, geometry may be specified in web interface
◆ Run standard handling maneuvers in ADAMS/Car
■ Export to ADAMS/Car
◆ Interface for reading in an ADAMS/Car property file into an
ADAMS/Car full vehicle model completed.
◆ Need to develop interface to write out these property files from
within ADVISOR (or any other source)
■ NREL
◆ Aaron Brooker
◆ Terry Penney
◆ Sam Sprik
◆ Keith Wipke
Users Conference for ADVISOR, Costa Mesa, CA, August 2000
Feedback
■ Comments or questions?
■ Contact Information
◆ [email protected]
◆ 949-786-0600 x100
QFD
Target Vehicle
Vehicle Targets Virtual Prototype
Comparison Vehicle
Sub-
System
Targets
Sub-System
Virtual Prototype
Component Targets
Component Virtual Prototype
Full
ADVISOR with ADAMS/Car, a product of Mechanical
Dynamics Inc. (MDI). ADVISOR (Advanced Vehicle
Test Road
SimulatOR) was created by National Renewable Energy Rigs Driver
Lab (NREL) in 1994 to support the DOE (Department of
Energy) Hybrid Program. ADVISOR, created in the Chassi Engine Driveli Body
MATLAB/Simulink environment, simulates conventional,
electric or hybrid electric vehicles during driving cycles to
determine the impact on emissions and fuel economy. Suspe Valvet Trans Body-
ADVISOR is a template driven environment where Steeri
nsion Crank
rain Clutch
missio in-
validated models can be altered to determine the Brake
ng Chain/
train Differ
Tires Acc. n
Axles/ white
Restr
sensitivity of design to changes, a new design analyzed s Belt ential
or a design optimized towards specific goals. It has been Drives CV aints
extensively used by many organizations around the
world, with over 1700 copies downloaded from the Packaging Ride& Loads- Mass Vibrations
world-wide-web. Handling Reduction
Wheel Torque
ADAMS/Controls Vehicle Speed
Co-simulation
ADVISOR
Hybrid Powertrains
Energy
Management
Controls
Figure 12: The ADAMS/CAR Application § 4WD/AWD powertrains (torque split actively
controlled by ADVISOR)
§ Vehicle handling/dynamics with new mass
distribution
Test
Rigs
Full Road
Driver
§
§
Stability issues related to battery placement
Determine energy losses during handling
and durability events
Chassis Power Plant Driveline Body § Determine accessory load energy impact vs.
performance
§ Asses impact of vehicle/component mass
reduction on dynamic performance
Steeri Clutc Fram
ng
Tires Acc. h e EXPORT TO ADAMS/CAR
CO-SIMULATION
CONTACT
REFERENCES
[1] Bidrawn, Les and Kohnle, Thomas, "The New 3," BMW
Magazine, Vol. 2, 1998, pp. 8-14
The power demand from the road is requested of the After a variable has been selected as known or
battery pack and the alternator. The power through the unknown, the Solution Path algorithm is called. This
alternator is adjusted using a constant alternator routine determines whether the known and unknown
efficiency, and the alternator then requests power from variables constitute a solvable equation set. If it is
the APU. determined that an equation set can be solved, the
Solution Path determines the best solution strategy. This
The battery pack was modeled with an open circuit routine is partially based on work completed by Ramirez
battery voltage, internal resistance and load voltage. A and Vestal [4].
voltage loop equation was written for the circuit. The
power through the circuit was determined by multiply the The solution strategy makes decisions based upon how
voltage loop equation by the current. The current was many unknown variables exist in each equation.
then found using the positive root of the quadratic Equations with only one unknown variable are selected
equation determined from the loop equation with internal as being the first equations to solve. Each equation with
resistance and open circuit voltage held constant. only one unknown is sent individually to the numerical
solving routine. The numerical solving routine consists
of a standard Newton-Raphson implicit method for
systems of equations with a Gauss-Jordan
decomposition routine. The results from the equations
with only one unknown variable are then used in Table 1 lists the results from four different design
subsequent calculations. When systems do not have iterations based on the input parameters and
equations containing only one unknown variable, the performance goals listed in Table 2.
Solution Path routine determines the sequence for
solving that minimizes the number of simultaneous Figure 1. SmartHEV GUI
equations sent to the numerical solver. This provides a
high level of stability and efficiency to the numerical
solving routines.
Ricardo is currently creating a link between ADVISOR and WAVE, Ricardo's industry-
standard engine performance and 1-D gas dynamics software package. WAVE models
are based on detailed engine dimensions. The resulting co-simulation will allow
prediction of drive cycle performance for engines without performance data.
Furthermore, the effect of modifications such as turbocharger size and valve timing can
be examined.
Ricardo will present the details of the ADVISOR-WAVE link, including the required
modifications to ADVISOR. Fuel consumption of a Class 7/8 truck modeled with and
without this link will also be presented.
Analysis of PNGV Inverter Power and Current Requirements
Using ADVISOR
Gerald W. Davis
Advanced Technology Labs, Milwaukee, Rockwell Automation
After completion of the Simulation Parameters screen, The power requirement is the fourth plot window in Fig. 3
the user clicks the RUN button at the bottom right. This and can also be viewed by plotting the “ess_power” array
action causes ADVISOR’s Results screen to appear. as shown in
RESULTS FOR VARIOUS SIMULATION SCENARIOS Fig. 4. The plot shows that the peak value required by
the FUDS drive cycle was close to 55 kW, but expansion
This section looks at the simulation results for testing the of the plot shows that the peak value represented only a
PNGV vehicle with the FUDS and US06 drive cycles. single point
The FUDS drive cycle is the likely certification drive cycle
for HEV vehicles. The US06 drive cycle is a higher road A more informative view of the power requirements is
speed drive cycle and may be applied to HEVs at some seen in Fig. 5. This is a histogram of the power samples
point in the future. delivered by the motor controller system to satisfy the
FUDS drive cycle. The bin width in the histogram is Fig. 7 shows a similar power occurrence histogram, but
about 400 W, and the large peak centered at 0 W is due as applied to the US06 drive cycle with 6.5% grade data.
to the zero speed command sections in the drive cycle. The histogram shows a high density of power level
Note that the vehicle accessories file was adjusted for occurrences around 60 kW. The difference between the
this particular simulation to draw zero power. A plot of requested speed command and actual speed command
the histogram of power values without the large peak is was also significant as shown in Fig. 8. Hence, the 55
shown in Fig. 6. kW inverter would not satisfy the US06 drive cycle with
the current PNGV vehicle weight.
RMS Power, Watts
Counts
Seconds
Power, Watts
Fig. 4. RMS Power requirement for results shown in Fig.
3. Fig. 7. Histogram of Motor Controller power occurrence
density vs power level derived from the US06 at 6.5%
performance data.
Counts
Power, Watts
CONCLUSIONS
ACKNOWLEDGMENTS
CONTACT
Individual
Conditions Optimization Evaluation
Suggested
Cost of individual(s)
Demands individual
Indiv.nr.
TABLE 1.
comp Cost of components.
fuel Fuel cost.
pollution Cost of pollution.
grid Cost of grid energy.
wc, wf, wg, wp Cost weights.
performance How well the driving cycle is managed. 1 if the vehicle can follow
the desired speed else 0.1 or 0.00017 etc. .
grams of
pollution per
km
limit