Electric and PlugIn Hybrid Vehicles Advanced Simulation - Methodologies
Electric and PlugIn Hybrid Vehicles Advanced Simulation - Methodologies
Electric and
Plug-In Hybrid
Vehicles
Advanced Simulation Methodologies
Green Energy and Technology
More information about this series at http://www.springer.com/series/8059
Bogdan Ovidiu Varga Florin Mariasiu
•
123
Bogdan Ovidiu Varga Dan Moldovanu
Automotive and Transport Department Automotive and Transport Department
Technical University of Cluj-Napoca Technical University of Cluj-Napoca
Cluj-Napoca, Cluj Cluj-Napoca, Cluj
Romania Romania
When Bogdan Varga asked me to write the Foreword to his book on “Electric and
Plug-In Hybrid Vehicles Advanced Simulation Methodologies” I was indeed
honored: Generally, people who are asked to write Forewords are those who matter
to the author.
I met Bogdan Varga in 2007 at a conference for automotive engineers organized
by the Technical University of Cluj. I had just taken over the job of Sales Manager
for Central and Eastern Europe for Instrumentation and Test Systems. I was
impressed by the determination and ardor Bogdan showed in his scientific work as
well as his professionalism and effectiveness in preparing big research projects.
In 2010, AVL together with Technical University of Cluj built a Laboratory for
Testing, Research, and Certification of Internal Combustion Engines, AVL deliv-
ering not only testing equipment but also research tools as well as software for data
processing and simulation.
Since then the Technical University of Cluj and AVL have been in close contact,
exchanging knowledge and experience in research, focusing mainly on methods to
deal with the increasing complexity of powertrains.
The complexity and diversity of future powertrain systems elevates dramatically.
One driving factor is given by enhanced legislative requirements, such as reduction
in CO2 fleet emissions and pollutant emissions under real-world driving conditions
(RDE).
Another factor is increased customer demands as positive driving experience,
safety, agility, comfort, and confidence in driver assistance systems. On top of this,
shorter model cycles and moderate cost of ownership are demanded.
One key to managing this complexity is to use and integrate virtual development
environments into the development process of powertrain systems.
In the past, for example, the functional integration of subsystems into a complete
vehicle system was done within a system integration using only real hardware.
Today this is combined with virtual environments, too. In future this integration
will be done in purely virtual environments in early stages of development.
v
vi Foreword
With the book’s approach Bogdan provides an excellent tool for researchers
working on these challenges: he not only describes the simulation method of
complex powertrains, but also explains how to build a working virtual development
environment for classical, hybrid, and electrical powertrains.
I am convinced that Bogdan’s book will be a big help for any researcher to better
understand the process of modeling and simulation of complex powertrains and to
start out in creating a virtual laboratory.
The writing process of this book consumed a huge amount of work, research, and
documentation, but in the end, the effort is dedicated to those who are interested in
applicative use of advanced modeling and simulation technologies for the optimi-
zation of energy efficiency in automotive electric and hybrid powertrains.
First of all, we are deeply grateful to AVL LIST GmbH, A-8020 Graz,
Hans-List-Platz 1 (www.avl.com) for their essential logistic support, availability,
and also for provision of their high-level professional expertise and technical
support in the writing process of the book. Their advice definitely increased the
scientific and applicative quality of book.
Still, all of this would not have been possible if we did not have a strong and
permanent support of leadership (Prof. Nicolae Burnete, Ph.D.) and our colleagues
from Automotive Engineering and Transport Department of Technical University
of Cluj-Napoca, Romania. Therefore, we would like to extend our sincere thanks to
all of them.
Nevertheless, we express our gratitude toward our families for their support,
understanding, and encouragement, which helped us permanently throughout the
preparation of the book.
vii
Contents
ix
x Contents
xv
xvi Abbreviations
Once the observations based on the source system can provide a model to
characterize its evolution, according to the variation of the variables considered, it
can proceed to the next phase, namely the simulation process.
The simulation process can be defined, in general terms, as all functional and/or
operational activities on a model that representing the image of a source system
necessary to be analyzed.
Because the model attached to a source system consists of one or more math-
ematical expressions (in the form of equations that describe the relationships
between model components), simulation processes are performed using computer
simulation.
There are a total of three main activities in the realization of a simulation
process: development, testing, and operation (Fig. 1.2).
The development work refers to the integration of mathematical models that
define the system state and variables into informatics codes (lines of code, pro-
grams) that can be processed further using computers.
In the testing phase of the model is required a check of its behavior at specified
conditions (if implementation of the model takes into account the initial
assumptions).
The results should be analyzed carefully to achieve assessing the fairness of the
model, by analogy with the data obtained direct from source system (under certain
known conditions).
If the model is valid and the simulation process is functional then the third phase
is run, the operating activities of the model.
The operation of model may possibly impose some conditions on the source
system studied, conditions that would not have been possible through direct
experiments on the physical source system.
The results should be analyzed and used both to improve the model and to
optimize source system behavior depending on the initial objectives of the research.
4 1 Principles of Modeling and Simulation Processes
into account multiple hypotheses and wide range of determinants factors) using
ultrafast informatics programs on electronic processing in computerized
environments.
For this, numerical methods are applied to discretize equations that define the
condition and/or evolution of the considered technical system. Using high potential
and processing speed of information available from today’s computers (for analysis
of complex data sets), equations that offer multiple connects between different
factors of different sizes can be determined in real time.
Using this method, the solving procedures of complex problems in processing
are substantially reduced, to temporal sequences of discretized sizes, with imme-
diate benefit to achieve economies of time and labor. Mathematical developed code
is implemented further as one of the appropriate programming languages.
Some of these programming languages have the advantage of object-oriented
programming capacity and accessible syntax for user. Applications developed using
programming codes allow simulation of functional systems and sub-systems that
compose the functional structure of a vehicle.
Using software tools existing in this area, simulations can be developed, con-
trolled, and coordinated (at the same time) to achieve results and solutions to
optimize and improve the performance of vehicle’s technical systems and are
analyzed in terms of performance and operating efficiency.
These applications allow or build linear or nonlinear models (depending on the
complexity and processing potential of software program), computer processing
lines and matrix, characterized by the boundary conditions at various points (areas),
and by the conditions of the initial state of the considered technical system.
Modeling and detailed analysis of the powertrain system has been, is, and will be
considered permanently in research worldwide. Due to complexity of the technical
structure of a powertrain, research focus is on the study of a single component, or of
the entire powertrain system, or by studying the interactions of this with the general
constructive-functional system of a vehicle [5, 6].
Walker and Zhang [7] studied the influence of engine torque harmonics, model
degrees of freedom, and dual mass flywheels on the transient response of a vehicle
powertrain equipped with a dual clutch transmission for two powertrain models.
Their results by modeling and simulation process indicate that degrees of freedom,
engine model, and flywheel model all contribute significantly to variance in pow-
ertrain response under each configuration.
Kouroussis et al. [8] present an effective formulation to model an automatic
transmission (based on Del’s method), to evaluate the vehicle performance in the
early development stage of a powertrain design. The efficiency and accuracy of
modeling process was validated taking into consideration the application of sim-
ulation results on two realistic cases (Chrysler 45RFE and Aisin Warner 55–50 SN
powertrains).
To analyze the advantages of the hybrid powertrain based on the component
characteristics and to define an optimum powertrain configuration with respect
to the applied drive cycle, a combined simulation and analytical analysis of the
6 1 Principles of Modeling and Simulation Processes
be investigated and then used further to optimize the functional and technical
parameters of the studied system.
Modeling and simulation approach of technical and functional characteristics
that influence the powertrain is widely presented in the subsequent chapters of this
book. For modeling and simulation processes the software packages CRUISE
(AVL List GmbH) and CarMaker (IPG GmbH) were used. These software pack-
ages were developed for precision applications in analysis and dynamic perfor-
mance of a vehicle’s technical systems.
Detailed methods, methodologies, and approaches to the immediate application
of modeling and simulation processes are presented, to analyze the vehicle’s
technical systems and/or sub-systems of which it is composed, in order to provide a
powerful tool to any researcher that wants to use immediately and apply these
modeling and simulation software packages in research activities.
The structure of each chapter allows any researcher to make easy the process of
modeling and simulation of a vehicle’s technical system, according to its specific
needs and requirements. This is possible due to detailed presentation of the algo-
rithms for construction, operation, and running of the presented models.
References
1. Mariasiu F (2013) Consumers’ Attitudes Related to Biofuel Use in Transportation. Int Rev
Manage Mark 3(1):1–9
2. Mariasiu F (2012) Energy Sources Management and Future Automotive Technologies:
Environmental Impact. Int J Energy econ Policy 2(4):342–347
3. Iclodean CD, Mariasiu F (2014) Possibility to increase biofuels energy efficiency used for
compression ignition engines fueling, TEM J 3(1):36–41
4. Moldovanu D, Burnete N (2013) Computational fluid dynamics simulation of a single cylinder
research engine working with biodiesel. Therm Sci J 17(1):195–203
5. Varga BO (2012) Energy management of electric and hybrid vehicles dependent on
powertrain configuration. Central Eur J Eng 2(2):253–263
6. Varga BO, Iclodean C (2014) Numerical investigation of hybrid electric vehicle’s
performance. Proc Int Cong Sci Manage Automot Transp Eng SMAT 1:27–24
7. Walker PD, Zhang N (2013) Modelling of dual clutch transmission equipped powertrains for
shift transient simulations. Mech Mach Theory 60:47–59
8. Kouroussis G, Dehombreux P, Verlinden O (2015) Vehicle and powertrain dynamics analysis
with an automatic gearbox. Mech Mach Theory 83:109–124
9. Katrasnik T (2007) Hybridization of powertrain and downsizing of IC engine—A way
to reduce fuel consumption and pollutant emissions—Part 1. Energy Convers Manag
48:1411–1423
10. Damiani L, Repetto M, Prato AP (2014) Improvement of powertrain efficiency through energy
breakdown analysis. Appl Energy 121:252–263
11. Klein CE, Bittencourt M, dos Santos Coelho L (2015) Wavenet using artificial bee colony
applied to modeling of truck engine powertrain components. Eng Appl Artif Intell 41:41–55
12. Li L, Zhang Y, Yang C, Jiao X, Zhang L, Song J (2015) Hybrid genetic algorithm-based
optimization of powertrain and control parameters of plug-in hybrid electric bus. J Franklin
Inst 352:776–801
8 1 Principles of Modeling and Simulation Processes
13. Wang B, Xu M, Yang L (2014) Study on the economic and environmental benefits of different
EV powertrain topologies. Energy Convers Manag 86:916–926
14. Ryu J, Park Y, Sunwoo M (2010) Electric powertrain modeling of a fuel cell hybrid electric
vehicle and development of a power distribution algorithm based on driving mode recognition.
J Power Sources 195:5735–5748
15. Xie C, Xu X, Bujlo P, Shen D, Zhao H, Quan S (2015) Fuel cell and lithium iron phosphate
battery hybrid powertrain with an ultra capacitor bank using direct parallel structure. J Power
Sources 279:487–494
Chapter 2
Mathematics Behind the Models
The Vehicle is one of the main objects in a model. This component contains general
data of the vehicle, such as nominal dimensions and weights. One and only one
vehicle component is needed in the model. Road resistances and dynamic wheel
loads are calculated for road and dynamometer runs based on the dimensions and
the load state. The wheel loads are calculated, taking motion into account (e.g.,
from acceleration, aerodynamic drag, rolling resistance). The aerodynamic, rolling,
climbing, acceleration, and total resistance are calculated [1].
2.1.1 Properties
Switch Variation
With this switch the vehicle can be given free for variation. For the vehicle the drag
coefficient and the vehicle mass can be varied. The setup of the variation parameters
is done in the folder.
Switch Output
If this switch is activated a result output for this component is made.
Driving Resistance Switches
Activation of one or more of the following five driving resistance models. The
associated resistance model can be calculated in one or more AVL CRUISE
Calculation Task when selected. Also the required data input for the activated
resistance models has to be done. In every task the selection button ‘Driving
Resistance’ determines which of the activated models will be really used in the task
calculation. Every model which is chosen in a task, has also to be activated in the
vehicle properties window, otherwise the calculation stops with an error message.
Switch Driving Resistance: Physical
The driving resistance is defined by the physical values (drag coefficient, rolling
resistance).
Switch Driving Resistance: Function with Reference Vehicle
The driving resistance is defined by three parameters. The defined values are
converted from the reference vehicle to the actual one.
Switch Driving Resistance: Characteristic with Reference Vehicle
The driving resistance is defined in a characteristic as function of the vehicle
velocity. The defined values are converted from the reference vehicle to the actual
one.
Switch Driving Resistance: Function without Reference Vehicle
The driving resistance is defined by three parameters. The defined values are used
without a reference vehicle.
Switch Driving Resistance: Characteristic without Reference Vehicle
The driving resistance is defined in a characteristic as function of the vehicle
velocity. The defined values are used without a reference vehicle.
Selection Button Driving Resistance Input Options
Standard
From Coasting Characteristic
From Deceleration Interval Characteristic
‘Normal’ input of the required values for Resistance Function or Resistance Table.
Selection Button Driving Resistance Input Options
Standard
From Coasting Characteristic
From Deceleration Interval Characteristic
When the switches ‘Driving Resistance: Function with Reference Vehicle’ or
‘Driving Resistance: Characteristic with Reference Vehicle’ are activated, the
option ‘Driving Resistance from Coasting Characteristic’ can be chosen. By
selecting a conversion button (located at the Coasting Characteristic), the
Resistance Function/Resistance Table is calculated using the values of the Coasting
Characteristic. The Coasting Characteristic is defined by vehicle velocity,
depending on the time since the coasting measurement started.
2.1 Vehicle (V) 11
Remark: for stationary iterations the upper hysteresis values are used.
Selection Button Aerodynamic Coefficients
Drag Coefficient Constant
Additional Aerodynamic Coefficients—Data Bus Dependent
Drag Characteristic
In this case the drag coefficient is interpolated from the drag characteristic,
depending on the value of the Data Bus input x-axis for drag modification.
12 2 Mathematics Behind the Models
The Tank Volume is the gas (fuel) tank volume of the vehicle. It is used to
calculate the distance range the vehicle can reach with the actual fuel consumption.
The temperature and the pressure difference serve for the correction of the intake
conditions. The differences in the conditions of the environment and the engine
compartment have to be defined. These values are used in the power correction on
environment conditions (see Chapter ‘Power Correction on Environment
Conditions’ in engine component).
Nominal Dimensions
The distance between hitch and front axle is the horizontal distance between the
coupling point for the trailer and the front axle. In the case of a Semi-trailer Truck,
this value can be smaller than the wheel base.
Wheel Base is the horizontal distance between the front and the rear axle
(Fig. 2.1).
h V,cog
h V,vcp
l V,cog
z
l V,fr
+M
l V,vcp
x
+M
x
h V,cd
Every calculation can be run with three different load states (zV,load): empty, half,
and full
• Empty: means that the weight of the vehicle equals the Curb weight
• Full: means that the weight of the vehicle equals the Gross weight
• Half: means that the weight of the vehicle equals the middle of the Curb and the
Gross weight
Load-Dependent Characteristics
lV,cog (zV,load) Distance of Gravity Center of the Vehicle Dependent on the Load mm
State
hV,cog (zV,load) Height of Gravity Center of the vehicle dependent on the Load State mm
hV,vcp (zV,load) Height of Hitch dependent on the Load State mm
PV,front (zV,load) Tire Inflation Pressure Front Axle bar
PV,rear (zV,load) Tire Inflation Pressure Rear Axle bar
If an additional mass is used the dimensions depending on the load state (dis-
tance and height of gravity center, height of hitch) will be interpolated between the
given values.
Nominal Weight
The Curb weight is the weight of the empty vehicle including a full fuel tank.
Av Frontal Area m2
The Frontal Area is the cross-sectional area of the vehicle that is used for the
calculation of the air resistance. With the frontal area the influence of the size of the
vehicle is considered.
The Drag Coefficient is the factor of the air resistance that depends on the shape
of the vehicle.
The lift coefficients of the front and the rear axle consider the influence of the
vehicle velocity on the wheel loads. Usually, the wheel loads are getting smaller
with higher velocity. So a positive lift coefficient means that the wheel loads are
decreasing. If the lift coefficients are negative then the wheel loads are increasing
with increasing velocity.
The air resistance acts in the track surface. Therewith the horizontal air resistance
has no influence on the wheel loads. To consider the effects of the different heights
of the real acting point the lift coefficients of the axles have to be used.
Additional Aerodynamic Coefficients—Data Bus Dependent
Further option for modification of aerodynamic parameters (drag coefficient, lift
coefficient front axle, and lift coefficient rear axle) which can be modified according
to the Data Bus input value ‘x-axis for drag modification’:
that means opening and closing of aerodynamic spoiler as function of any Data Bus set
value (with hysteresis).
16 2 Mathematics Behind the Models
Threshold Values
In these table columns the values are defined which the aerodynamic coefficients
will have at the end of the modification time interval.
Lateral Dynamics Vehicle
Chassis Data
Dependent on curvature radius the vehicle must generate lateral forces, which
are impacted by the lateral acceleration. For calculation of the sideslip angle and
thus the steering angle the stiffness of the front and rear axle is needed. For the
calculation of the downforces of each wheel the track width is also needed.
Crosswind Force Data
The curvature radius is not only responsible for the lateral forces. If cornering is
chosen, the influence of sidewind can be calculated. Dependent on the angle of
approach the above subscripted values are recommended.
2.1 Vehicle (V) 17
Drag Characteristic
Here the time-dependent temperature characteristics for the 1st cycle are defined.
Temperature Curves Cycle 2
Here the time-dependent temperature characteristics for the 2nd cycle are
defined.
18 2 Mathematics Behind the Models
Here the time-dependent temperature characteristics for the 3rd cycle are
defined.
Temperature Curves Cycle 4
Here the time-dependent temperature characteristics for the 4th cycle are defined.
Temperature Curves Cycle 5
Here the time-dependent temperature characteristics for the 5th cycle are defined.
Temperature Values—Cycle Independent
The Resistance Function can also be converted from the Coasting Characteristic,
see “Conversion from Coasting Characteristic to Resistance Function/Resistance
Table.”
2.1 Vehicle (V) 19
Resistance Table
mV,ref Weight kg
The Resistance Table can also be converted from the coasting characteristic, see
“Conversion from Coasting Characteristic to Resistance Function/Resistance Table.”
Conversion from Coasting Characteristic to Resistance Function/Resistance
Table
This feature converts the required values from the so-called ‘Deceleration
Interval Characteristic’; every input row of this characteristic defines the time for a
20 2 Mathematics Behind the Models
aV Acceleration m/s2
sV Distance m
vV Velocity m/s
TV,temp,cycle,i Temperature Curve 1,…,5 for Cycle i K
PV Power W
PV,loss Power Loss W
For evaluating the forces, which are acting on the vehicle, the current load conditions
are important. They are specified for each calculated task in two different ways.
(a) Load conditions fixed as load stage [2
The actual vehicle mass
mV;act ¼ mV zV;load ð2:1:1Þ
The actual values for the height and the horizontal distance are determined by
interpolation between the three load stages.
(b) For the chassis dynamometer [2]:
Front Axis
The summation of all moments around the contact point between rear wheels and
road has to be zero [2]:
X
Mr;ax ¼ 0 ð2:1:8Þ
Out of this equation we can calculate the load of the front axis [2].
lV;cog;act hV;cog;act
Fw;x;f;ax ¼ mV;act 1 g cos aU ðaV þ g sin au Þ
lV;fr lV;fr
FV;lift;f
ð2:1:9Þ
The load on the front axis is allocated proportionate onto the single front wheels
for the right and the left side.
Right front wheel load [2]:
FW;x;f;ax
FW;x;f;r ¼ ð2:1:10Þ
2 NW;f;r
FW;x;f;ax
FW;x;f;l ¼ ð2:1:11Þ
2 NW;f;l
Rear Axis
The summation of all moments around the contact point between front wheels and
road has to be zero [2]:
X
Mf;ax ¼ 0 ð2:1:12Þ
From this equation we can calculate the load of the rear axis [2].
lV;cog;act g cos aU hC;cog;act
Fw;x;r;ax ¼ mV;act þ ðaV þ g sin aU Þ FV;lift;r
lV;fr lV;fr
ð2:1:13Þ
The load on the rear axis is allocated proportionate onto the single rear wheels for
the right and the left side.
24 2 Mathematics Behind the Models
FW;x;r;ax
FW;x;r;r ¼ ð2:1:14Þ
NW;r;r
FW;x;r;ax
FW;x;r;l ¼ ð2:1:15Þ
NW;r;l
Additional traction or pushing forces are taken into account with the relative
forces kV,add,trac and kV,add,push. They are related to the weight of the vehicle and
their size is fixed by the equation solver for special Calculation Tasks (e.g., Max
Traction Force).
Summation [2]
FV;res ¼ FV;air þ FV;incl þ kv;add;trac þ kv;add;push mV;act g ð2:1:20Þ
mact mact cW A
FV;res ¼ cA þ cB vV þ cC v2V ð2:1:21Þ
mref mref cW;ref Aref
ð vV Þ
F ð vÞ ¼ A þ B v V þ F ð2:1:22Þ
B ¼ F 0 ðvV ¼ 0Þ ð2:1:24Þ
These parts are related to the actual vehicle mass and they can be handled
separately.
The function part of higher order is calculated [2]:
ðvV Þ ¼ F ðvV Þ A B vV
F ð2:1:25Þ
Now the actual resistance force can be evaluated for each time step [2]:
mact ð vV Þ
FV;res ¼ ðA þ B vV Þ þ F ð2:1:26Þ
mref
FV;res;ges ¼ FV;res þ FV;incl þ km;add;trac þ km;add;push mV;act g ð2:1:27Þ
ð vV Þ
F ð v Þ ¼ A þ B vV þ F ð2:1:29Þ
B ¼ F 0 ðvV ¼ 0Þ ð2:1:31Þ
26 2 Mathematics Behind the Models
These parts are related to the actual vehicle mass and they can be handled
separately.
The function part of higher order is calculated [2]:
ð v V Þ ¼ F ð vV Þ A B vV
F ð2:1:32Þ
Now the actual resistance force can be evaluated for each time step [2]:
ðv V Þ
FV;res ¼ ðA þ B vV Þ þ F ð2:1:33Þ
FV;res;ges ¼ FV;res þ FV;incl þ km;add;trac þ km;add;push mV;act g ð2:1:34Þ
The single parts of the lateral forces can be calculated by the following
formulas:
Front (rear) axle force by lateral acceleration [2]:
The angle αCompass is the angle of the vehicle since starting (which could
change). The start value can be chosen under Wind Velocity.
Lateral force influence by sideslip angle [2].
2
vVehicle
Fs;b ¼ cs q AS sin bslip ð2:1:38Þ
2
The pressure point of the crosswind is important for the splitting on the axles
of crosswind force. It is measured from the front of the vehicle [2].
cMZ
dlateral ¼ l ð2:1:39Þ
cs
dlateral lfront;frontaxle
FS;front ¼ Fs;Wind þ Fs;b ð2:1:40Þ
bwheelbase
hGC
FD;f;rðlÞ ¼ FD;f;rðlÞ ðþÞ FS;front ð2:1:43Þ
tf
hGC
FD;r;rðlÞ ¼ FD;r;rðlÞ ðþÞ FS;rear ð2:1:44Þ
tr
(3) Sideslip angel front (rear) under the influence of the radial force [2]:
FS;fðrÞ
ca
afðrÞ ¼ Fradial;fðrÞ ð2:1:45Þ
1 þ ca
bWheelbase
aackermann ¼ ð2:1:46Þ
q
(5) Steering angle [2]:
aackermann
d¼ Fradial;f þ af ar ð2:1:47Þ
1þ ca
(6) Slip angle [2]:
bwheelbase lGC
b¼ ar ð2:1:48Þ
d
Stationary idle, transition to motion and interruption of the power flow are all made
possible by the clutch. The clutch slips to compensate for the difference in the
rotational speeds of engine and drivetrain when the vehicle is set in motion. When a
change in operation conditions makes it necessary to change gears, the clutch
disengages the engine from the drivetrain for the duration of the procedure.
The clutch contains a model for a friction clutch as it is used in cars with manual
gearboxes. In this case the clutch is controlled by the driver via the cockpit [1].
The last possibility of using the clutch is together with a CVT gearbox. There it
is only necessary for starting. In this case the clutch is controlled by the control
module CVT control [1].
2.2 Clutch (C) 29
2.2.1 Properties
Switch Output
If this switch is activated a result output for this component is made.
Selection Button Model
simple
detailed
For the simple model of the clutch the maximum transferable torque has to be
defined.
Selection Button Model
simple
detailed
At the detailed model of the clutch the transferable torque is derived from the
geometric data of the clutch as well as from the friction characteristic.
Selection Button Control Variable
Desired Clutch Release
Desired Torque
This is the ‘usual’ selection where the clutch gets the value of the desired clutch
release. Therefore, the Data Bus input Desired Clutch Release has to be connected.
Selection Button Control Variable
Desired Clutch Release
Desired Torque
This mode considers the desired torque to be transmitted. Therefore, the Data Bus
input.
Torque Demand has to be connected.
Switch Dynamic Mode
If this switch is activated the clutch calculation always stays dynamic and never
switches to a kinematic connection.
Switch Free Definable Friction Characteristic
If this switch is activated and depending on whether ‘simple’ or ‘detailed’ model
has been chosen, the characteristic ‘Maximum Transferable Torque’ overslip, or the
characteristic ‘Friction Coefficient’ overslip can be defined.
30 2 Mathematics Behind the Models
The number of sets of frictional surfaces is equal to the number of the clutch
disks. In the calculation every disk is set with two frictional surfaces.
From these values the friction characteristic of the clutch is determined. The
input of the two friction coefficients is only necessary if the switch ‘Free definable
Friction Characteristic’ is deactivated (Fig. 2.3).
The maximum transferable torque is the maximum torque the clutch can transfer
between the drive and the power take-off side. Its input is only necessary if the
switch ‘Free definable Friction Characteristic’ is deactivated.
Pressure Force
The pressure force is the axial force the frictional surfaces are pressed together.
This pressure force depends on the clutch release. With the pressure force it is
possible to calculate the transmitted torque by considering the dimensions.
θ C,in θ C,out
MC,in MC,out
MC
MC;max u_ c;rel Maximum Transferable Torque overslip –
Friction Coefficient
Friction coefficient has to be evaluated for stick and slid friction. Which value is
taken in consideration depends on the circumstances that the clutch is locked or
unlocked.
They are fixed in the code with the following values:
MC;max
rC;m ¼ ð2:2:1Þ
NC lC;st FC
The friction gradient is also fixed in the code for typical conditions:
CC ¼ 0:01
The actual friction coefficient is calculated by using the following formula [2]:
ju_ c;rel jCC
l
c;st lc;sl
lC;act ¼ lC;sl þ lC;st lC;sl e ð2:2:3Þ
2.2 Clutch (C) 33
For the adhering clutch the clutch torque is determinate by the value which is
applied by the system [2]:
Torque converters employ the force represented by a moving fluid to transmit engine
torque. Because these devices compensate for differences in the rotating speeds of
engine and drivetrain, they are ideal for effecting the transition from stationary to
mobile operation. The torque converter also multiplies torque. First an impeller
(pump) converts the mechanical energy emanating from the power unit into fluid
energy (hydraulic fluid—ATF—is the preferred medium); a second transformation,
back into mechanical energy, occurs at the blades within the turbine [1].
Among the benefits offered by torque converters are the following attributes:
infinitely variable, stepless variations in torque and rpm, vibration insulation,
absorption of torque peaks, and virtually wear-free power transfer [1].
These devices offer economical operation when used together with mechanical
variable-ratio transmissions; operation with nonpositive engagement (slip) is lim-
ited in the interests of efficiency [1].
The converter contains a model for a torque converter as it is, e.g., used in
vehicles with automatic gearboxes. For computing the output variables, the various
characteristic curves of the converter are used. The characteristic curves are usually
34 2 Mathematics Behind the Models
Lock-up Clutch
MT,lc
MT,in MT,out Power Take-off
Drive Side Side
MT,p MT,h
Torque Converter
2.3.1 Properties
Switch Output
If this switch is activated a result output for this component is made.
Switch Use for Reference
At the task Climbing Performance the maximum gradients are only calculated for
the coupling part (torque multiplication equals zero). If more than one torque
converter is in the model, the reference torque converter has to be determined with
this switch. This means that from this torque converter the coupling part is deter-
mined, for which the gradients are calculated.
Selection Button Lock-up Friction Model
Linear
Simple
In this case the ‘default’ model is used: the maximum transferable torque of the
lock-up clutch is calculated linearly through the clutch release value.
Selection Button Lock-up Friction Model
Linear
Simple
In this case the friction model of the clutch component is used.
2.3 Torque Converter (T) 35
The mass moments of inertia of the pump and the turbine side are given with the
oil share.
The reference angular velocity is the speed where the table for the torque con-
version is measured. Reference speed means that throughout the whole measuring
process the pump speed was fixed at this special reference angular velocity. For
other pump speeds the corresponding values are extrapolated.
The maximum lock-up torque is the maximum torque that can be transferred via
the lock-up clutch.
This torque can be transferred if the clutch is completely closed. If the clutch is
closed partly (controlled operation; refer to component Clutch Program (CP) for
more information), a correspondingly lower torque can be transferred.
Torque Conversion
Torque
Ratio
0
0 1 Speed Ratio
The torque ratio is the division between the turbine torque and the pump torque.
It depends on the speed ratio which is the division between the turbine and the
pump speed. The torque ratio is nearly one throughout the whole speed ratio. Only
when the speed ratio is nearly zero the torque ratio is increasing (refer to Fig. 2.5).
Torque Curve
These values are required when the selection button Input options is set to
torque curve. The pump torque is the transferable torque between the pump and the
turbine. It is a function of the speed ratio. At the speed ratio one, i.e., there is no
difference in speed between the pump and the turbine, the pump torque is zero. If
the speed ratio is between zero and one (the pump drives the turbine), the pump
torque is positive. If the speed ratio is above one (the turbine drives the pump), the
pump torque is negative.
In addition, there is a dependence on the ratio between the pump speed and the
reference speed. This is calculated in the following way [2]:
!2
u_ T;Pump
MT;Pump;act ¼ MT;Pump ð2:3:1Þ
u_ T;P;n
Cf-Curve
These values are required when the selection button Input options is set to
Cf-Curve. In the calculation core, the Cf-values are converted to the pump torque
assuming the reference speed as being 2000 rpm.
The relation between the Cf-values and the pump torque at reference speed is
described by the following equation [2]:
M T;Pump;act 106
Cf T ¼ ð2:3:2Þ
u_ 2T;Pump
K-Curve
These values are required when the selection button Input options is set to
K-Curve. In the calculation core, the K-factors are converted to the pump torque
assuming the reference speed is 2000 rpm.
The relation between the K-factors and the pump torque at reference speed is
described by the following equation [2]:
u_ T;Pump
KT ¼ pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ð2:3:3Þ
MT;Pump;act
The K-values can only be positive (due to their definition). AVL CRUISE calcu-
lates pump torque for speed ratios above 1 by mirroring the pump torque around the
point (1/0) in the speed ratio—pump torque—plane.
The torque on the converter input and output depends on the speed ratio between
pump and turbine as well on the absolute angular velocity of the pump.
Speed Ratio
This is the difference between the angular velocities of the clutch input and output
side.
The actual speed ratio is evaluated:
(1) For pump speed is zero u_ T;in ¼ 0 [2]:
(a) (u_ T;out \0)
mT;act ¼ 1 ð2:3:6Þ
(2) With an input speed not equal zero u_ T;in 6¼ 0 [2]:
u_ T;out
mT;act ¼ ð2:3:7Þ
u_ T;in
Pump Torque
Now the transferable torque can be evaluated with a liner interpolation for the pump
torque map [2]:
MT;pump;nom ¼ MT;pump vT;act ð2:3:8Þ
The scale factor KT,scale is for transformation from the nominal pump to the actual
pump speed [2]:
u_ T;in
KT;scale ¼ ð2:3:9Þ
u_ T;p;n
Turbine Torque
With a linear interpolation in the torque ratio map we evaluate the actual torque
ratio between pump and turbine [2]:
lT;act ¼ CT iT;Tor vT;act 1 þ 1 ð2:3:11Þ
Clutch Torque
The clutch torque is the transmitted torque of the lock-up clutch. If the lock-up
clutch is acting the whole transmitted torque is split into the pump-turbine part and
the clutch part. The size of the transmittable clutch torque depends linear on the
clutch release [2].
Transmittable Torque
The transmittable torque of the converter is the summation of the actual pump and
the clutch torque [2].
available gears, the spacing of the individual ratios within the required conversion
range and the engine’s full load torque curve [1].
The component Gearbox contains a model for a gearbox with different gear steps.
You can define as many gears as you want. For every gear it is possible to define the
transmission ratio, the mass moments of inertia, and the moment of loss [1].
In the component for manual gearboxes, the engine torque will be turned into a
power take-off torque by considering the transmission, the mass moments of inertia,
and the moment of loss [1].
The gearbox can be used for a manual or automatic gearbox. When used as an
automatic gearbox, the gear shifting process will be controlled by the control
module gearbox control or gearbox program. The driver will do this task when used
as a manual gearbox [1].
2.4.1 Properties
Switch Variation
The gearbox can be given free for variation with this switch. For the gearbox, the
transmission ratio can be changed by multiplication with an additional factor. The
setup of the variation parameters is done in the folder.
Switch Output
If this switch is activated a result output for this component is made.
Switch Gear Shift Time
If this switch is activated, the map for the gear shift times is used during calculation.
Selection Button Losses
Deactivated
Efficiency
Efficiency and Torque Loss Correction
Torque Loss Map
Torque Loss Map with Torque Division
Torque Loss Map with Torque and Speed Division
Torque Loss Map temperature dependent
In this case no losses are calculated in the gearbox.
Selection Button Losses
Deactivated
Efficiency
Efficiency and Torque Loss Correction
Torque Loss Map
Torque Loss Map with Torque Division
Torque Loss Map with Torque and Speed Division
Torque Loss Map temperature dependent
42 2 Mathematics Behind the Models
In this case, the losses are calculated by a constant efficiency which only depends
on the gear position.
Selection Button Losses
Deactivated
Efficiency
Efficiency and Torque Loss Correction
Torque Loss Map
Torque Loss Map with Torque Division
Torque Loss Map with Torque and Speed Division
Torque Loss Map temperature dependent
In this case, the losses are calculated by the efficiency which only depends on the
gear position and by a torque loss map which depends on the gear position and the
drive speed.
Selection Button Losses
Deactivated
Efficiency
Efficiency and Torque Loss Correction
Torque Loss Map
Torque Loss Map with Torque Division
Torque Loss Map with Torque and Speed Division
Torque Loss Map temperature dependent
In this case the losses are interpolated out of a torque loss map, which depends on
the gear position, the drive speed, and the drive torque. At the torque and speed
division a special approximation of the maps is made by cutting through the map at
the defined torque and speed values.
Selection Button Losses
Deactivated
Efficiency
Efficiency and Torque Loss Correction
Torque Loss Map
Torque Loss Map with Torque Division
Torque Loss Map with Torque and Speed Division
Torque Loss Map temperature dependent
In this case, the losses are calculated with the usage of up to 5 torque loss maps
Gear Losses Temperature Dependent where each map is associated with a
user-defined constant temperature value. An interpolation between these maps is
done so that for every time step a torque loss can be generated.
The temperature used in the calculation comes through the Data Bus input
Temperature of the component Gearbox. For example, it could be connected with
one of the up to 5 cycle-dependent temperature characteristics in the component
Vehicle.
2.4 GearBox (G) 43
This selection is only available if the selection button Losses is set to Torque Loss
Map temperature dependent.
Table Numerical Values of the Torque Loss Partition
For the torque and the speed range the borders for activation and deactivation can
be defined here. The values are defined in percentage of the maximum torque and
speed values.
NG Gear Position –
The input stands for the gear number (1 for 1st gear, 2 for 2nd gear,…). Zero
means the neutral gear. It has always the transmission one. This neutral gear is
needed to consider the mass moments of inertia when the vehicle is only rolling.
The program determines the number of gear steps automatically.
The transmissions can be defined for the different gear positions. However the
transmission for the neutral gear (gear position 0) has to be 1. If the number of teeth
is defined, the transmission ratio of each gear step is automatically determined.
The mass moments of inertia on the drive and the power take-off sides can be
defined different for each single gear step. But it is also possible to define the same
mass moments of inertia for every single gear step.
The number of teeth can be defined here. If the user puts in the transmission
ratio, AVL CRUISE detects automatically the best fitting number of teeth to get the
defined transmission ratio. If it is not possible to get the transmission ratio exactly a
suggestion is made for number of teeth which gives the closest possible trans-
mission ratio.
2.4 GearBox (G) 45
If the switch Gear Shift Time is activated, the map of the times for upshifting
and downshifting dependent on the actual gear number is used. For gears where
times are not defined, a default value of 0.1 s is assumed.
Gear Shifting Efficiency for AT Models
When the selection button Gear Shifting Losses is set to Efficiency, the losses
during the shifting processes are calculated.
Efficiency and Torque Loss Correction
The efficiency can be defined separately for the single gear steps. This infor-
mation is used if the selection button Losses is set to Efficiency or Efficiency and
Torque Loss.
M G;V;3 u_ G;in ; N G Torque Loss in the single gear steps Nm
The additional torque loss is defined dependent on the gear step and the drive
speed of the gearbox. These values are used if the selection button Losses is set to
Efficiency and Torque Loss.
Torque Loss Correction Temperature Dependent
M G;V;4 u_ G;in ; N G;TG Torque Loss temperature dependent Nm
The torque losses are defined for every single gear step also considering the
drive speed and the drive torque of the gearbox.
This map is used if the selection button Losses is set to Torque Loss Map,
Torque Loss Map with Torque Division, or Torque Loss Map with Torque and
Speed Division.
Torque Loss Map Temperature Dependent
TG,i Temperature i = 1, …, 5 °C
For each activated, temperature dependent torque loss map, the temperature level
has to be specified.
M G;V;5;i u_ G;in;i ; M G;in ; N G; i Torque loss map temperature dependent i = 1,…,5 Nm
Here for each activated temperature-dependent torque loss map and for the
specified temperature level, the losses are defined depending on speed, input torque
and gear. This is only required if the selection button Losses is set to Torque Loss
Map temperature dependent.
When the switch Efficiency Map (for Torque Loss Maps) is activated, the
efficiency for gear losses and temperature-dependent gear losses can be entered in a
separate table instead of the torque loss. After pushing the button Conversion
Efficiency in Torque Loss, the torque loss map is created automatically.
u_ in
u_ out ¼ ð2:4:2Þ
iG;act
€ in
u
€ out ¼
u ð2:4:3Þ
iG;act
Inertia Moments
The mass moments of inertia are also selected out of the table [2].
Torque Loss
There are four kinds to take the torque loss of the gearbox in consideration.
(1) Calculation without losses (ZG,i = 1) [2]
MG;loss ¼ 0 ð2:4:6Þ
gG ¼ 1 ð2:4:7Þ
MG;loss;nom
gG ¼ 1 ð2:4:10Þ
MG;in
2.4 GearBox (G) 49
(3) Calculation with use of the efficiency value and the torque loss (ZG,i = 3)
The rotational speed contribution to the torque loss is linear interpolated in
the loss curve [2]:
MG;loss;n ¼ MG;V;3 u_ G;in ; NG;act ð2:4:11Þ
Is the simulation done with temperature-dependent losses then the torque loss
is interpolated out of the map MG,V,4 (u_ G;in ; NG ; TG ) [2]:
MG;loss;n ¼ MG;V;4 u_ G;in ; NG;act ; TG ð2:4:12Þ
Also the efficiency value is selected for the actual gear [2]:
calculation time is decreasing. The disadvantage is that the transitions might be not
real harmonic and that there is the need to change the equation system for every
change in transmission [1].
Because there is no possibility to define an unlimited transmission there is a
clutch needed for starting. This clutch will be controlled by the control module
CVT control which is also controlling the CVT gearbox [1].
2.5.1 Properties
Switch Output
If this switch is activated a result output for this component is made.
Selection Button Power Transmission
without slip
with slip
The calculation of the CVT does not include chain slip. A purely kinematic cou-
pling is performed.
Selection Button Power Transmission
without slip
with slip
The calculation of the CVT includes also chain slip. A dynamic coupling between
input and output is made.
Selection Button Losses at “Without Slip”
Deactivated
Efficiency
Efficiency and Torque Loss Correction
Efficiency and Torque Loss Correction Temperature Dependent
Torque Loss Map
Torque Loss Map and Torque Loss Correction Temperature Dependent
Torque Loss Map Temperature Dependent
At this selection no losses are considered in the CVT when calculating without slip.
Selection Button Losses at “Without Slip”
Deactivated
Efficiency
Efficiency and Torque Loss Correction
Efficiency and Torque Loss Correction Temperature Dependent
Torque Loss Map
Torque Loss Map and Torque Loss Correction Temperature Dependent
Torque Loss Map Temperature Dependent
At this selection only the efficiency map is considered when calculating without
slip.
2.5 CVT—Continuously Variable Transmission (H) 51
through the Data Bus input Temperature. For example, it could be connected with
one of the 5 cycle-dependent temperature characteristics in the component Vehicle.
Switch Efficiency Map (for Torque Loss Maps)
When this switch is activated, the efficiency can be entered in a separate table
instead of the torque loss.
After pushing the button Conversion Efficiency in Torque Loss, the torque loss
map is created automatically.
Switches Torque Loss Map 1, …, 5
With these switches, the associated loss maps can be activated for up to 5 tem-
perature levels. This selection is only available if the selection button Losses is set
to Torque Loss Map Temperature Dependent.
Selection Button Adjustment Time
Fixed Value
Speed Dependent
Speed Dependent from Data Bus
The adjustment time is the time the CVT needs to adjust from minimum to max-
imum transmission or vice versa. At this selection the adjustment time is a fixed
value.
Selection Button Adjustment Time
Fixed Value
Speed Dependent
Speed Dependent from Data Bus
The adjustment time depends on the input speed of the CVT.
Selection Button Adjustment Time
Fixed Value
Speed Dependent
Speed Dependent from Data Bus
The adjustment time depends on the speed connectivity from the Data Bus.
TH Adjustment Time s
The adjustment time is the time the CVT needs to adjust from minimum to
maximum transmission or vice versa.
Here the minimum and the maximum transmission have to be defined. Between
these two threshold values the transmission can change infinitely, i.e., the trans-
mission can reach every possible value between the two thresholds.
The mass moments of inertia have to be defined separately for the drive and the
power take-off side (Fig. 2.6).
Slip Description
Is active when the CVT is calculated with slip.
The Slip Characteristic describes the slip behavior of the CVT as relative map
which depends on its nominal values.
MH,pd.in
A iH,act
MH,in
B
Slip s
MH,out
θ H,in
.
iH,act θ H,out
MH,pd.out
ƞH,V,2(iH) Efficiency –
The additional torque loss is defined as function of the transmission ratio and the
drive speed.
These data can only be used for a calculation of a CVT without slip.
Torque Loss Correction Temperature Dependent
MH;V;3 u_ H;in ; iH ; TH Torque loss temperature dependent Nm
The torque lost in the transmission is defined in a table depending on the drive
speed, the drive torque and the actual transmission ratio. These data can only be
used for a calculation of a CVT with slip.
Torque Loss Map Temperature Dependent
TH,i Temperature i = 1, …, 5 °C
For each activated temperature-dependent torque loss map, the temperature level
has to be specified.
M H;V;5;1 u_ H;in;i ; M H;in;i ; N H;i Torque loss map temperature dependent i = 1, …, 5 Nm
Here for each activated temperature-dependent torque loss map and for the
specified temperature level, the losses are defined depending on speed, input torque
and gear. This is only required if the selection button Losses is set to Torque Loss
Map Temperature Dependent.
When the switch Efficiency Map (for Torque Loss Maps) is activated, the
efficiency for gear losses and temperature-dependent gear losses can be entered in a
separate table instead of the torque loss. After pushing the button Conversion
Efficiency in Torque Loss, the torque loss map is created automatically.
56 2 Mathematics Behind the Models
TH Adjustment Time s
TH u_ Adjustment time Speed Dependent s
The adjustment time is the time the CVT needs to adjust from minimum to
maximum transmission or vice versa.
If this gear shifting step is great enough and the gear ratio limits are not reached
the gear shifting process will be done.
(a) Calculation without slip
The calculation without slip is done with kinematic equations only. The input
and output speed relation is fixed. Also the input and output torque are defined
by the transmission ratio.
(b) Calculation with slip
In this case also the force equations are used (dynamic simulation):
For the actual ratio the reference slip is calculated with the following formula
out of the Slip Characteristic [2]:
sH;diff
sH;act ¼ ð2:5:2Þ
sH;ref
sH;ref ¼ sH;nom sH;diff iH;act iH;nom þ 1 ð2:5:3Þ
Using the relation between actual and nominal CVT-torque the slip relation
can be evaluated.
Now the actual slip is determined for the reference slip and the slip relation [2].
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
2ffi
u_ ref ¼ u_ 2break þ ju_ Belt j þ Du_
Torque Loss
There are three kinds to take the torque loss of the gearbox in consideration.
(1) Calculation without losses (ZH,i = 1) [2]
MH;loss ¼ 0 ð2:5:5Þ
gH ¼ 1 ð2:5:6Þ
MH;loss;nom
gH ¼ 1 ð2:5:9Þ
MH;in
(3) Calculation with use of the efficiency value and the torque loss (ZG,i = 3)
The nominal loss is linear interpolated in the loss curve [2]:
MH;loss;nom ¼ MH;V;3 u_ H;in ; iH;act ð2:5:11Þ
MH /MH,nom
0
0 1 SH /sH,nom
Also the efficiency value is selected for the actual gear [2]:
gH ¼ gH;V;3 iH;act ð2:5:13Þ
The Single Ratio Transmission is a gear step with fixed ratio. It can be used, e.g., as
transmission step of the differential (final drive unit).
A drive torque will be transferred to a power take-off torque of the transmission
step by considering the transmission, the mass moments of inertia and the moment
of loss [1].
2.6.1 Properties
Switch Variation
With this switch the single ratio transmission can be given free for variation. For the
single ratio transmission the transmission ratio can be changed by multiplication
with an additional factor. The setup of the variation parameters is done in the folder.
Switch Output
If this switch is activated a result output for this component is made.
60 2 Mathematics Behind the Models
In this case a torque partition is done. A special approximation of the maps is made by
cutting through the map at the defined torque values. This selection is only available if
the selection button Losses is set to Torque Loss Map temperature dependent.
Selection Button Torque Loss Partition
Deactivated
Torque Partition
Torque and Speed Partition
In this case a torque and speed partition is done. A special approximation of the
maps is made by cutting through the map at the defined torque and speed values.
This selection is only available if the selection button Losses is set to Torque Loss
Map temperature dependent.
Table Division Values for the Torque Loss
For the torque and the speed range the borders for activation and deactivation can
be defined here. The values are defined in percentage of the maximum speed and
torque values.
iD Transmission Ratio –
The overall ratio can be defined if this case of input is chosen in the properties.
The user can put in the transmission ratio or the number of teeth depending on
the selection in the properties (Fig. 2.8).
MD,z
θ D,out
MH,out
2.6 Single Ratio Transmission (D) 63
ƞD,V,3 Efficiency –
Torque Loss on the drive side. The torque loss is defined considering the drive
speed and the drive torque of the single ratio transmission.
Torque Loss Map Temperature Dependent
TD,i Temperature i = 1, …, 5 °C
Here for each activated, temperature dependent torque loss map and for the
specified temperature level, the losses are defined depending on speed and input
torque. This is only required if the selection button Losses is set to Torque Loss
Map temperature dependent.
When the switch Efficiency Map (for Torque Loss Maps) is activated, the
efficiency can be entered in a separate table instead of the torque loss.
64 2 Mathematics Behind the Models
TD Temperature °C
If the numbers of teeth are defined for the single transmission step, then the gear
Z
ratio can be calculated: iD ¼ ZD;2
D;1
.
u_ D;in
u_ D;out ¼ ð2:6:1Þ
iD
€ D;in
u
€ D;out ¼
u ð2:6:2Þ
iD
There are four kinds to take the torque loss of the gearbox in consideration.
(1) Calculation without losses (ZD,i = 1) [2]
MD;loss ¼ 0 ð2:6:3Þ
gD ¼ 1 ð2:6:4Þ
MD;loss;nom
gD ¼ 1 ð2:6:7Þ
MD;in
66 2 Mathematics Behind the Models
(3) Calculation with use of the efficiency value and the torque loss (ZD,i = 3)
The rotational speed contribution to the torque loss is linear interpolated in the
loss curve [2]:
MD;loss;nom ¼ MD;V;3 u_ D;in ð2:6:8Þ
The differential unit compensates for discrepancies in the respective rotation rates of
the drive wheels: between inside and outside wheels during cornering and between
different drive axles on 4WD vehicles [1].
With rare exceptions for special applications, the differential is a bevel gear drive
unit. When the output bevel gears on the left and right sides (most common
arrangement) are of equal dimensions, the differential gears act as a balance arm to
equalize the distribution of torque to the left and right wheels [1].
When unilateral variations in road surfaces result in different coefficients of
friction at the respective wheels, this balance effect limits the effective drive torque
to a level defined as twice the traction force available at the wheel (tire) with the
lower coefficient of friction. This wheel then responds to the application of
excessive force by spinning. To avoid such effects a positive lock is available at the
component [1].
In the differential it is possible to define a torque split factor (Fig. 2.9).
This is required if it is used as a central differential for a four-wheel drive, as
sometimes you want to have either more torque transferred to the front or the rear
axle. This is used for preserving the driving quality of the corresponding two wheel
driven car which is either front wheel driven or rear wheel driven [1].
2.7 Differential (N) 67
Drive Side
MN,i
MN,o.1 MN,o.2
Power Power
Take-off Take-off
Side 1 Side 2
2.7.1 Properties
Switch Output
If this switch is activated a result output for this component is made.
Selection Button Losses
Deactivated
Efficiency
In this case, the calculation of losses is deactivated.
Selection Button Losses
Deactivated
Efficiency
In this case, the losses are calculated through a user-defined input of a constant
‘Stationary Efficiency’ value.
This button allows the user to run the differential in a locked or an unlocked
mode or to get the torque split factor from Data Bus.
Unlocked mode means that the two power take-off torques are the same and the
speeds can be different. In the locked mode there is a rigid connection between the
68 2 Mathematics Behind the Models
unlocked locked
two power take-off sides. That means that the speeds are the same and the torques
can be different (Fig. 2.10).
To have a different torque output on the power take-off side, choose the cor-
responding splitting ratio. To have the same torques on both outgoing sides, the
splitting ratio equals 1.
Stationary Efficiency
gN;stat Efficiency –
The stationary efficiency η is the efficiency of the differential with fixed input
shaft.
If the properties switch ‘Losses’ is set to ‘Efficiency’, AVL CRUISE calculates
in every time step the actual losses depending on the fixed efficiency η.
2.7 Differential (N) 69
Gear Ratio
At first the split ratio for the inversion is determined [2]:
1
iN;inv ¼ ð2:7:1Þ
iN
The second ratio which is calculated is the ratio for the summation [2]:
iN;sum ¼ 1 þ iN ð2:7:2Þ
Split Equations
Unlocked differential
For the unlocked differential the following equation for the torque is used [2]:
MN;lock ¼ 0 ð2:7:4Þ
Locked differential
For the locked differential both takeoffs have the same speed [2]:
The Planetary Gearbox belongs to the power split devices. It consists of three main
components: the sun gear, the planet carrier (with planets), and the ring gear. This
gearbox can be used in hybrid systems, where vehicles be driven using different
power suppliers. With the help of a planetary gear it is possible to add torque and
speed with different directions at the transmission input shafts [1].
In a planetary gear system, the speed ratio and the direction of rotation can be
changed according to which member is fixed. There are three types of planetary
gear mechanisms, depending upon which member is locked (locked ring gear, or
locked planet carrier, or locked sun gear) [1].
2.8 Planetary GearBox (PG) 71
2.8.1 Properties
Switch Output
If this switch is activated a result output for this component is made.
Selection Button Input Mode
Base Ratio
Number of Teeth
In this case, the gear ratio between ring gear and sun gear, the so-called base ratio of
a planetary gearbox, are required.
Selection Button Input Mode
Base Ratio
Number of Teeth
In this case, the number of sun gear’s teeth and the number of ring gear’s teeth are
required.
Selection Button Losses
Deactivated
Efficiency
In this case, the calculation of losses is deactivated.
Selection Button Losses
Deactivated
Efficiency
In this case, the losses are calculated through a user-defined input of a constant
‘Stationary Efficiency’ value.
This button allows the user to run the planetary gearbox in a locked or an
unlocked mode.
Unlocked mode means that on the three shafts (planet carrier, sun, and ring)
torques can be different and the speeds can also be different. In the locked mode
there is a rigid connection between two of them.
That means that the speeds on all shafts are the same (because of the speed’s
equation at the planetary gear) and the torques can be different.
Base ratio is quotient of ring gear’s teeth to the sun gear’s teeth.
The number of teeth for the ring gear is negative because of its internal toothing.
Stationary Efficiency
ƞPG,stat Efficiency –
Under the balance conditions, the summation of all torques acting on planetary
gear is equal to zero by the stationary state of motion [2]:
Two of the three torques always possess the same sign. The summation of their
absolute values is equal to the third torque. The shaft, which this summation leads
to, is defined as summation shaft. The others are difference shafts.
There are two characteristic rules, which are very useful to analyze the complex
interaction on double pinion gearbox:
The summation shaft’s torque and difference shafts torques has opposite signs.
The torques of both Difference Shafts have the same sign.
The component engine contains a model for an internal combustion engine. As the
characteristic curves for the full load, the fuel consumption and others can be freely
defined by the user. It is possible to define a gasoline engine as well as a diesel
engine [1].
74 2 Mathematics Behind the Models
2.9.1 Properties
Switch Variation
With this switch the internal combustion engine can be given free for variation.
For the internal combustion engine the engine displacement can be varied. The
setup of the variation parameters is done in the folder.
Switch Output
A result output for this component is made if this switch is activated.
Switch Use for Reference
This switch is used to select one reference engine if more engines are used within a
model. Out from this engine the following data are calculated, which can then be
used in the Calculation Tasks:
• Minimum speed
• Maximum speed
• Maximum torque
• Speed at maximum torque
• Maximum power
• Speed at maximum power
Selection Button Intended for Calculation of
Performance
Performance and Consumption
Performance, Consumption, and Emission
Performance, Consumption, and CO2 Emission
Performance, Consumption, Emission, and CO2 Emission
By selecting this switch, different calculations are executed for which the corre-
sponding input data have to be defined. For example, if “Performance” is selected
just the general data for the engine like displacement, the Full Load Characteristic
and the Motoring Curve have to be defined.
Selection Button Temperature Model
Deactivated
Cold Start Correction
Others…
In this case no temperature calculation and no correction of the fuel consumption or
the emissions will take place.
2.9 Internal Combustion Engine (E) 75
In this case a power correction on environment conditions is made. For this cor-
rection the standard 97/21/EG (April 1997) is used, which was originally developed
for the standardization of power measurements.
Switch Specific Consumption Map
When this switch is activated the specific Fuel Consumption Map is active.
Switch Full Load Gear Dependent 1–5
For defined gear positions, extra Full Load Characteristics can be defined. They are
used instead of the original ones.
Selection Button Full Load Reduction Gear Dependent
Off
Factor
From Characteristic
In this case the ‘Full Load Reduction Characteristics’ have no influence.
Selection Button Full Load Reduction Gear Dependent
Off
Factor
From Characteristic
In this case for every activated characteristic ‘Full Load Gear Dependent 1,…,5’ a
reduction factor has to be defined. The values from the gear dependent character-
istic are multiplied with this factor.
Selection Button Full Load Reduction Gear Dependent
Off
Factor
From Characteristic
In this case for every activated characteristic ‘Full Load Gear Dependent 1,…,5’ an
additional ‘Full Load Reduction Characteristic’ has to be defined. The values from
the reduction characteristic are subtracted from the gear dependent characteristic.
Switch Engine Brake Curve
This switch activates the engine brake curve. With this brake curve it is possible to
simulate a jake brake system. This works with changing of the valve opening
period, when the vehicle brakes. With the Data Bus input Jake Brake Activation the
consideration of the brake curve can be controlled during calculation (1..activated,
0..deactivated).
Switch Starter Current
When this switch is activated, up to 5 temperature-dependent starter current char-
acteristics can be used when the engine is started. Each characteristic is associated
with a user-defined constant temperature value. An interpolation between these
maps is done so that for every time step a current depending on the actual tem-
perature can be determined.
82 2 Mathematics Behind the Models
The engine working temperature is needed for a calculation with hot start. If hot
start is chosen the engine temperature at the beginning of the calculation is equal to
the engine working temperature. If cold start is chosen the engine temperature at the
beginning of the calculation is equal to the ambient temperature. This engine
working temperature means always the temperature of the equivalent mass which is
used for the calculation of the engine temperature.
The mass moment of inertia contains all parts of the engine like the crankshaft,
the flywheel, a possible fan, the camshaft, and others.
The response time is the time the engine needs to build up the full power.
It is possible to define a maximum speed for the engine. If the engine speed is
above this speed the fuel injection will be stopped.
Fuel Type
Here the characteristic data of the fuel like the density or the heating value are
input. They are used for temperature and emission calculations of the motor.
Idle
Here the fixed value for the idle speed can be defined. This value is used when
the selection button Idle Speed is set to Fixed Value.
Here the fixed values for the idle consumption and the idle emissions can be
defined.
These values are used when the selection button Consumption/Emissions at
Idle is set to fixed value.
Fuel Shut-Off
The fuel shut-off is used to switch off the fuel injection while the engine is in thrust
operation. With this the fuel consumption can be decreased.
Selection Button Fuel Shut-Off
Absolute Speed Limits
This activates the definition of absolute lower and upper speed values for fuel shut-off.
Relative Speed Differences
Selection of this option activates the input fields ‘Speed Difference above Idle
Speed for Lower Speed Limit’ and ‘Speed Difference above Idle Speed for Upper
Speed Limit.’
2.9 Internal Combustion Engine (E) 85
With these speed difference definitions, a definition of the speed limits relative to
the idle speed can be done.
This is useful if the idle speed is not fixed to a constant value during the
calculation, e.g., when a temperature-dependent idle speed characteristic is used.
Cut-Off Signal from Data Bus
The user can create his own cut-off control and connect the cut-off signal via
Data Bus to the engine component via the new Data Bus input ‘Fuel Shut-Off
Activation.’
When the Engine is in thrust mode and the Data Bus signal is not equal to zero,
the engine is set to cut-off mode and the ‘Residual Fuel Consumption’ is taken for
the consumption calculation. In all other cases, the engine is not in cut-off mode.
Comments
Engine-specific information can be entered, stored, and viewed in a ‘comments’
table. The table consists of the columns ‘Description | Value | Unit.’ When adding a
row to the table, the window of a global, ‘pre-defined’ comments table appears.
In the first step, definitions for the pre-defined comments table have to be done in
the form ‘Attribute | Description | Data Type | Unit.’
In the second step, rows from this global table can be selected to be added as
rows in the engine-specific comments table for the actual engine.
Friction Mean Effective Pressure (FMEP)
CE,bore Bore mm
CE,hub Hub mm
rE,compress Compression Ratio –
SAE 10W/30/10W/40/10W/60
SAE 15W/30
Selection Button Fuel Injection Pump installed
no
yes
Pre-defined Temperature Curve
Here the temperature characteristic of the engine can be defined by the user.
Cooling Characteristic
The cooling characteristic is used for the determination of the amount of heat
which is lost to the environment.
The heat capacity has to be given as a function of the actual temperature of the
engine.
Warm-up Enrichment
ΔbE,start,enrich(t) Consumption enrichment during the first time after starting l/h
The start enrichment is defined as function of time. The start enrichment can be
switched on and off with the switch for the start enrichment.
Acceleration Enrichment
Here the idle speed can be defined as function of the engine temperature. This
characteristic is only active when the selection button Idle Speed is set to From Idle
Speed Characteristic.
Boost Pressure Build-up Time
or
This is the time the charger needs to build up the full boost pressure.
Boost Pressure
This value is required if the option ‘‘Boost Pressure at Idle” in the Properties
window is activated.
pE;ch u_ E;out Boost Pressure as absolute pressure bar
This is the engines BMEP at the end of BMEP build-up time. It has to be defined
as percentages of maximum BMEP at a certain speed. This speed is user-defined (in
case that build-up time as function of speed is chosen), otherwise (when constant
build-up time is chosen) a reference speed of 1000 rpm is considered.
This value has to be defined if the selection button Response Behavior Turbo
Charger is set to Constant BMEP Build-up Time.
tE;ch:buildupBMEP u_ E;out BMEP Build-up Time as function of speed s
This value has to be defined if the selection button Response Behavior Turbo
Charger is set to BMEP Build-up Time From Characteristic.
Charger Outlet Temperature
tE;ch u_ E;out Temperature behind charger °C
The idle consumption map is defined as function of the brake mean pressure with
the engine speed as parameter. It is only used when the selection button
Consumption/Emissions at Idle is set to From Idle Map.
Idle Emission Maps
eE;idle;NOX u_ E;out ; peff Idle Emission Map NOx kg/h
eE;idle;CO u_ E;out ; peff Idle Emission Map CO kg/h
eE;idle;HC u_ E;out ; peff Idle Emission Map HC kg/h
eE;idle;Soot u_ E;out ; peff Idle Emission Map Soot kg/h
The idle emission maps are defined as function of the brake mean pressure with
the engine speed as parameter. They are only used when the selection button
Consumption/Emissions at Idle is set to From Idle Map.
90 2 Mathematics Behind the Models
The idle emission maps are defined as a function of the brake mean pressure with
the engine speed as parameter.
eE;idle;exh u_ E;out ; peff Idle Exhaust Mass Flow Characteristic kg/h
TE;idle;exh u_ E;out ; peff Idle Exhaust Temperature Characteristic °C
At the Fuel Consumption Map the actual fuel consumption for stationary
operation is listed depending on the brake mean pressure with the engine speed as
parameter. This kind of representation is used to make a unique relation between the
brake mean pressure and the fuel consumption for a given engine speed.
For the idle consumption there are different possibilities. By choose it can be
interpolated out of the Fuel Consumption Map or the defined value is used. If the
mean pressure is out of the defined range the corresponding values are extrapolated.
Specific Consumption Map
bE;spec u_ E;out ; peff Specific Fuel Consumption Map l/kWh or kg/kWh
Emission Maps
eE;NOX u_ E;out ; peff NOx Emission Map kg/h
bE;CO u_ E;out ; peff CO Emission Map kg/h
bE;HC u_ E;out ; peff HC Emission Map kg/h
bE;Soot u_ E;out ; peff Soot Emission Map kg/h
In the emission maps the emissions for the different exhaust gases are listed
dependent on the mean pressure with the engine speed as parameter.
Extended Engine Maps
The following maps are only used when the selection button Exhaust System
Model is set to Advanced AVL Exhaust System Model.
eE;O2 u_ E;out ; peff O2 Emission Map kg/h
eE;H2 O u_ E;out ; peff H2O Emission Map kg/h
eE;CO2 u_ E;out ; peff CO2 Emission Map kg/h
eE;H2 u_ E;out ; peff H2 Emission Ma kg/h
In the emission maps the emissions for the different exhaust gases are listed
dependent on the mean pressure with the engine speed as parameter.
eE;exhaust u_ E;out ; peff Exhaust Mass Flow Characteristic kg/h
TE;exhaust u_ E;out ; peff Exhaust Temperature Characteristic °C
The environment conditions at the test bench run are used for the power cor-
rection on environment conditions.
Full Load Characteristic optionally:
PE;vk u_ E;out Full Load Characteristic as power W
or
ME;vk u_ E;out Full Load Characteristic as torque Nm
92 2 Mathematics Behind the Models
or
pE;eff;vk u_ E;out Full Load Characteristic as BMEP bar
The full load characteristic can be defined either as power, as torque or as brake
mean pressure dependent on the engine speed.
Full Load Reduction
Constant Reduction Factor
The Full Load Characteristic is multiplied with this factor over the whole speed
range.
Full Load Reduction Characteristic optional:
PE;vk;red u_ E;out Full Load Reduction Characteristic as power W
or
ME;vk;red u_ E;out Full Load Reduction Characteristic as torque Nm
or
pE;eff;vk;red u_ E;out Full Load Reduction Characteristic as BMEP bar
The full load reduction characteristic is subtracted from the original Full Load
Characteristic.
Full Load Gear Dependent 1–5
NE,vk,geardep,i, Lowest Gear Position for which the gear dependent Full Load –
low Characteristic is active
NE,vk,geardep,i, Highest Gear Position for which the gear dependent Full Load –
high Characteristic is active
Lowest and highest gear position has to be defined for every activated gear
dependent Full Load Characteristic.
Optionally:
PE;vk;red;i u_ E;out Gear dependent Full load 1–5 as power W
or
ME;vk;red;i u_ E;out Gear dependent Full Load 1–5 as torque Nm
2.9 Internal Combustion Engine (E) 93
or
pE;eff;vk;red;i u_ E;out Gear dependent Full Load 1–5 as BMEP bar
The gear dependent Full Load Characteristics can be defined separately for the
single gears either as power, as torque or as brake mean pressure dependent on the
engine speed.
If one gear is associated with more than one of the five maps, then the first
associated map is used for calculation.
Gear Dependent Full Load Reduction Characteristic 1–5
Constant gear dependent Reduction Factor 1–5:
The associated gear dependent Full Load Characteristic is multiplied with this
factor over the whole speed range.
PE;vk;red;geardep u_ E;out Gear dependent Full Load Reduction Characteristic 1–5 W
The gear dependent full load reduction characteristics are subtracted from the
associated gear dependent Full Load Characteristics.
Motoring Curve
Optionally:
PE;sk u_ E;out Motoring Curve as power W
or
ME;sk u_ E;out Motoring Curve as torque Nm
or
pE;eff;sk u_ E;out Motoring Curve as BMEP bar
In the Motoring Curve, the friction of the engine is described considering the idle
torque the engine can deliver at the special engine speed. The Motoring Curve can
be defined similar to the Full Load Characteristic either as power, as torque or as
brake mean pressure dependent on the engine speed.
94 2 Mathematics Behind the Models
or
or
In the engine brake curve, the back pressure of the engine is described con-
sidering the brake torque the engine can deliver in the case of a braking. The power
loss curve can be defined similar to the load characteristic either as power, as torque
or as brake mean pressure dependent on the engine speed.
With the Data Bus input Jake Brake Activation the consideration of the brake
curve can be controlled during calculation (1…activated, 0…deactivated).
Starter Current
I E;start;i ðtÞ Sequence of the Starter Current at the starting process, defined for up to 5 A
temperature levels
The starter current characteristics are used for the start-stop automatic (if this
special component is used).
For each characteristic a temperature level has to specified; the current is defined
depending on time.
Load Signal Map
aE;th u_ E;out ; peff Load Signal Map %
The throttle valve positions can be defined dependent on the engine speed and
the brake mean pressure.
2.9 Internal Combustion Engine (E) 95
PE;vk
ME;vk ¼ ð2:9:1Þ
u_ E;out
PE;vk NE;stroke p
peff;vk ¼ ð2:9:2Þ
VE;h u_ E;out
The whole full load curve is multiplied with this factor [2].
Reduction curve
The reduction curve is subtracted from the full load curve [2].
If there exists a special full load curve for the selected gear then this curve is used
instead of the original one [2].
kE;operate ¼ 1 ð2:9:8Þ
Dp
peff;sk ¼ fE;sk u_ E;out þ pE;sk;help þ Dp fE;sk ð2:9:10Þ
u_ E;idle
Dp Dp
peff;vk ¼ peff;sk þ fE;sk 1 u_ 3E;out 2 fE;sk 2 u_ 2E;out
u_ E;idle
3 u_ E;idle
ð2:9:11Þ
2.9 Internal Combustion Engine (E) 99
kE;operate ¼ 1 ð2:9:13Þ
NE;stroke p
peff;vk ¼ pE;vk;help u_ E;max u_ E;out ð2:9:14Þ
VE;h
NE;stroke p
peff;sk ¼ pE;sk;help u_ E;max u_ E;out fE;sk ð2:9:15Þ
VE;h
kE;operate ¼ 0 ð2:9:17Þ
The power correction on environment conditions serves for the correction of the
Full Load Characteristic when the environment conditions the user puts in are
different from the conditions on the test bench. In this case a correction is made at
which the standard 97/21/EG (April 1997) is used, which was originally developed
for the standardization of power measurements.
For the power correction two selection buttons are used:
Selection Button Engine Type
Gasoline
Diesel
Selection Button Charger Type
Without
Turbo Charger
TC with Intercooler
With these selection buttons the factors for the power correction are fixed and the
needed input maps and characteristics are activated.
To make a power correction, the standard conditions must first be defined.
100 2 Mathematics Behind the Models
These are called environment conditions during test bench run in the page of the
Full Load Characteristic. If a charger or an intercooler is used the corresponding
characteristics for the charger pressure the temperature behind charger and the
temperature behind the intercooler have to be defined too. With these values the
conditions at the inlet for the test bench run and the efficiencies for the charger and
the intercooler can be calculated.
With the characteristics for the charger and the intercooler efficiency, the actual
inlet conditions can be calculated out of the actual environment conditions by
considering the temperature and pressure differences between the environment and
the power house.
With the actual inlet conditions and the ones on the test bench the correction
factors can be calculated and the Full Load Characteristic can be corrected.
(a) Inlet Conditions on the test bench
First the dry pressure of the environment at the test bench is calculated [2]:
with [2]:
58:73895 6852:493 5:262: ln TU;air
p0W ¼ 9:80665:e ð2:9:24Þ
T U;air
where the factors A, B, X, and Y depend on the engine and charger type.
• Engine Factor:
– For gasoline engines [2]:
fm ¼ 1 ð2:9:28Þ
apow;corr ¼ fA fm ð2:9:31Þ
Time [2]
Response
Dt
If peff;theo [ peff;theo;old þ peff;max tE;add
Dt
peff;theo ¼ peff;theo;old þ peff;max ð2:9:34Þ
tE;add
The charger response behavior works very similar to the power correction. Also a
comparison between different inlet conditions is made. These conditions are the
condition with delayed response of the turbo charger and an infinite response. The
infinite response means that always the full boost pressure for this speed is available.
The correction is made with the same basis as the power correction, i.e., also the
correction parameters form the standard 97/21/EG (April 1997) are used.
The calculation of the charger response behavior is done after the power
correction.
For the calculation without power correction the boundary conditions are defined
by the chassis dynamometer [2]:
For calculation with power correction the conditions from the engine trunk has
to be used [2]:
The boost pressure increases after a response time, which is calculated as fol-
lowed [2]:
pE;ch;over;dyno ¼ pE;inlet;dyno pE;env;dyno RWet TE;inlet;dyno qE;W;env;dyno
ð2:9:43Þ
Dt
pE;ch;over;act ¼ pE;ch;over;act þ pE;ch;over;dyno ð2:9:44Þ
tE;ch;buildup
For the response time a constant value or the value out of the map is used.
If the requested engine torque comes from the Data Bus as ‘Desired Torque,’ this
value is directly converted into peff,theo and then limited by peff,vk.
Dt
peff;theo ¼ peff;theo;old þ peff;max ð2:9:46Þ
tE;add
In AVL CRUISE, there are three options to determine the actual temperature of the
engine:
104 2 Mathematics Behind the Models
The heat produced by the combustion can be computed by considering the fuel
consumption. This heat is distributed among the engine block and exhaust gas by a
factor. The heating of the engine can be computed by considering the temperature
loss in the radiator. Thus, the engine temperature can be computed.
This technique has two advantages: First, the actual heat flow in the engine is
represented. Second, the required data can be measured on the testbed with but little
expenditure. No transient maps to be determined in a complex way are required.
For the computation of the engine temperature, the following ideas are used:
• The heat capacity of the injected fuel (fuel consumption of the warm engine plus
additional fuel consumption caused by the increasing friction mean pressure) is
the input to the engine system (engine + cooling system). To estimate the heat
capacity of this system, an equivalent mass of water (mE,eq) with the same heat
capacity has to be given.
• One part of the fuel energy is given off as work. The rest is divided into two
parts: The first part goes to the exhaust gas system (ZE,EE), the second part is the
input to the equivalent mass.
2.9 Internal Combustion Engine (E) 105
The heat produced by the combustion will be computed from this formula [2]:
PE;cool aA TE TW;air Dti
TE ¼ TE þ ð2:9:52Þ
mE;eq cp;H2 O
In AVL CRUISE, the following models and methods can be used for the calculation
of the higher fuel consumption while the engine is cold:
• Warm-up Enrichment
• Cold Start Correction
• FC External and FC Coefficient External
• Increasing Friction Mean Pressure
• Increase by Mean Pressure Factor
The actual fuel consumption is multiplied with cold start factor which is defined the
in the component vehicle. The factor can be defined for different cycles as a
function of time. In the task, the switch ‘Cold Start Correction’ has to be activated.
FC External
An additional fuel consumption from the Data Bus input channel ‘Additional FC
(Mass Flow)’ is added to the originally determined fuel consumption.
FC Coefficient External
The determined fuel consumption (inclusive additional external FC) is multiplied
with the FC coefficient delivered by the Data Bus input channel ‘FC Coefficient
External.’
The characteristic for the friction mean pressure is given by the friction mean
pressure at minimum (pE,min) and maximum speed (pE,max) and a curvature factor
(cE,fric,p) between them (cE,fric,p = 1 means a linear characteristic; refer to Fig. 2.11).
At first, a reference value for the mean friction pressure is calculated [2]:
2.9 Internal Combustion Engine (E) 107
pE;fric;max pE;fric;min
fB ¼ c c ð2:9:53Þ
u_ E;max E;fric;p u_ E;idle E;fric;p
c
fA ¼ pE;fric;min fB u_ E;idle E;fric;p ð2:9:54Þ
cE;fric;p
pE;fric;ref ¼ fA þ fB u_ E;out ð2:9:55Þ
The dependence on the temperature is done via the dynamic oil viscosity. Here,
only the oil viscosity at 40 °C (ηE,40) is needed. The oil viscosity at other tem-
peratures is calculated by means of an empirical equation [2]:
gE;40
fhelp ¼ 159:56 ln ð2:9:56Þ
1:8 104
fhelp fhelp
gðT Þ ¼ 1:8 104 e T273:15þ95 e 887 ð2:9:57Þ
gE;act ¼ gðTE ÞgE;ref ¼ g TE;N ð2:9:58Þ
For the conversion of the oil viscosity into the friction mean pressure, the
exponent of friction mean pressure (cE,vis,exp) is needed. The equation for is as
follows:
!
cE;fric;exp
gE;act
pE;fric;add ¼ pE;fric;ref 1 ð2:9:59Þ
gE;ref
The friction mean pressure, the Fuel Consumption Maps, and the emission maps
are given at the nominal temperature (TE,N). For the correction of the fuel con-
sumption the working point in the engine is changed. For this the change in friction,
mean pressure between the actual temperature and the nominal temperature is
108 2 Mathematics Behind the Models
calculated and added in the Fuel Consumption Map. So the working point is
changed to a higher (actual temperature is lower than the nominal temperature) or to
a lower (actual temperature is higher than the nominal temperature) mean pressure.
With this other working point, the fuel consumption is also increasing or
decreasing (Refer to Fig. 2.12).
Introduction
The principal factors influencing a vehicle’s fuel consumption can be determined
on the basis of this definition of consumption per distance.
They can be categorized as:
1. Contribution factors related to motor (Engine Losses),
2. Drivetrain Losses, and
3. Exterior resistance
There is continuing interest in better understanding the friction losses in internal
combustion engines. These adversely affect the maximum work output and fuel
2.9 Internal Combustion Engine (E) 109
B Bore m
S Stroke m
Db Bearing diameter m
Lb Bearing length m
nb No. of bearings –
nc No. of cylinders –
N Engine rotational Speed rpm
Ccb Coefficient of the hydrodynamic losses in main bearings Pa/rpm m
Ccs Coefficient of friction losses in main bearing seals Pa m2
Ctd Coefficient of friction losses due to viscous dissipation Pa s2/m2 rpm
110 2 Mathematics Behind the Models
3. Valve Train
The friction means effective pressure in the valve train is calculated using
the following equation [2]:
1:5 0:5
Nnb L N nv
Fmepcv ¼cvb 2 þ cvo þ Cvh v
B Snc BSnc
3
ð2:9:62Þ
10 Lv nv
þ Cvm 1 þ þ Fmepcam follower
N Snc
where:
nv No. of valves –
Cvf Flat cam follower constant Pa m
Cvr Roller cam follower constant Pa m
Both of the cam followers’ constants depend on the valve train configuration.
4. Auxiliary Losses
The following equation can be used to calculate the FMEP due to the
auxiliaries [2]:
Fmepaux ¼ a þ bN þ cN 2 ð2:9:65Þ
The constants α, β and γ are used to be 6.23 × 103 Pa, 5.22 Pa/rpm,
−1.79 × 10−4 Pa/rpm2, respectively.
B. Pumping Losses
The correlation used by Bohic has been used. That is [2]:
Disp
Pmep ¼ 105
0:0025
0:0785 þ 4:02105 N þ 1:06 108 N2 þ 4:64 108 IMEP þ 2:171010 NIMEP
ð2:9:66Þ
where: IMEP, Disp are the indicated mean effective pressure and engine
displacement, respectively.
III. Engine Losses: Modified PNH Model
The rate at which the engine structure and lubricating oil warm-up influences
friction losses over the cycle, because these depend on oil viscosity, which in turn
depends upon temperature. The friction power losses drop sharply over the first 50 s
of engine warm-up [4], after which the rate of fall is greatly reduced. Such behavior
indicates that friction level changes are likely to be influenced by the oil viscosity.
The effect of viscosity changes is taken into account by scaling the instantaneous
friction loss in proportion to values for fully warmed up [2].
0:24
FmepToil vToil
¼ ð2:9:67Þ
FmepToil ¼ 90 C vToil ¼ 90 C
3. Valve Train
Means effective pressure due to the friction in the piston group is expressed
as follows [2]:
l n
N 0:6 nb
Fmepvalvetrain ¼ cvb
B2 Snc
lref
1:5 0:5 n
Lv N nv l 10 Lv nv
þ Cvs þ Cvh þ Cvm 2 þ
BSnc lref 5 þ lN Snc
þ fmepcam=follower
ð2:9:70Þ
2.9 Internal Combustion Engine (E) 113
For oil pump, the constants α, β, γ & viscosity index (n) are: 1.28 kPa,
7.9 × 10−3 kPa mm3/rpm, −8.4 × 10−7 kPa mm3/rpm2, 0.3 (viscosity index),
respectively. For water pump α, β, γ & viscosity index (n) 0.13 kPa,
2 × 10−3 kPa mm3/rpm, 3 × 10−7 kPa mm3/rpm2, 0.7, respectively. They are:
1.72 kPa, 6.9 × 10−4 kPa mm3/rpm, 1.2 × 10−7 kPa mm3/rpm2, 0.5,
respectively for the oil pump.
5. Pumping Losses
The model was developed for motored engines. Therefore, the pumping
losses are zero.
The fuel consumption is linearly interpolated for the engine speed u_ E;out and the
mean effective pressure peff,act out of the fuel consumption map bE,h u_ E;out ; pE .
For u_ E;out \u_ E;idle and αth < 0, the fuel consumption can optionally be taken
out of the overall map, out of the detailed map, or can be the fixed idle
consumption.
114 2 Mathematics Behind the Models
The fuel shut-off is used to cut off the fuel injection while thrust operation. This is
used in modern cars to decrease the fuel consumption
For the fuel cut-off, two speeds are necessary to create a hysteresis so that the
fuel shut-off is not always activated and deactivated in too short times.
This speed is the lowest speed the shut-off can work on. If the engine speed falls
below this lower speed, the fuel shut-off is deactivated.
This speed has to be reached once that the fuel shut-off can be activated.
u_ E;SA;low;rel Speed Difference above Idle Speed for Lower Speed Limit rpm
u_ E;SA;high;rel Speed Difference above Idle Speed for Upper Speed Limit rpm
If ‘Relative Speed Differences’ are chosen, the upper and lower speed borders
are calculated with these values relative to the actual idle speed.
The residual fuel consumption is the fuel consumption while the fuel shut-off is
active.
This is, for example, needed if you have an additional heater.
The button for Consumption Increase after Deactivation is to change the way
the fuel consumption is increasing after the fuel shut-off is deactivated (refer to
Figs. 2.13 and 2.14).
(a) for u_ E;act u_ E;SA;low and ZE,SA,fc = activated [2]:
bE;idle bE;SA
bE;act ¼ bE;idle u_ E;act u_ E;idle ð2:9:74Þ
u_ E;SA;low u_ E;idle
2.9 Internal Combustion Engine (E) 115
(b) for u_ E;act [ u_ E;SA;low [2]:
2.9.10.1 Emissions
At first the factor of fuel consumption is evaluated. It is the relation between the
actual and the map consumption [2]:
bE;act
ffuel ¼ ð2:9:76Þ
bE;h
The emissions are linear interpolated out of the maps for u_ E;act and PE,act.
116 2 Mathematics Behind the Models
This values out of the map have to be transformed with ffuel [2]:
Here, the additional fuel consumption after starting the engine is interpolated out of
a user-defined characteristic.
Here, the additional fuel consumption while stepping on the acceleration pedal is
interpolated out of a user-defined characteristic.
The charge of battery charge caused by the engine starter is calculated [2]:
0:5 IE;start ðti1 Þ þ IE;start ðtiþ1 Þ
QE;Q ¼ QE;Q ð2:9:78Þ
IE;start ðti Þ IE;start ðti2 Þ
The generator (alternator) must provide the vehicle electrical system with a sufficient
supply of current under all operating conditions in order to ensure that the State of
Charge (SOC) in the engine storage device (battery) is consistently maintained at an
adequate level. The object is to achieve balanced charging, i.e., the curves for
performance and speed–frequency response must be selected to ensure that the
amount of current generated by the alternator under actual operating conditions is at
least equal to the consumption of all electrical equipment within the same period [1].
Automotive alternators are designed to supply charge voltages of 14 V (with
28 V for heavy utility vehicles) in order to maintain an adequate charge in 12 V (or
24 V) batteries.
The alternator produces alternating current. The vehicle’s electrical system, on
the other hand, requires direct current to recharge the battery and operate the
electrical equipment; it is thus direct current that must ultimately be supplied to the
electrical system. For this case a rectifier must be provided to convert the alter-
nator’s three-phase alternating current into DC.
This arrangement also provides the battery from discharging when the vehicle is
stationary. In the input data for the component generator the already into direct
2.10 Generator (L) 117
current converted characteristics have to be given. The characteristics used here are
measured behind the rectifier [1].
In the component generator a simple regulator is integrated. The regulator works
in that way that a threshold voltage has to be defined. If the voltage of the onboard
network is below this threshold the generator is activated.
The generator will be separated from the network if the voltage of the onboard
network reaches the threshold [1].
The current delivered by the generator results from a map dependent on the
speed and the mains voltage. By considering the internal resistance, the instanta-
neous current consumption helps to acquire the torque absorbed by the generator
with the corresponding moments of loss [1].
2.10.1 Properties
Switch Output
If this switch is activated a result output for this component is made.
Selection button Definition
Efficiency Map
Torque and Efficiency Characteristic
The description of the losses in the generator is done by the efficiency map.
Selection button Definition
Efficiency Map
Torque and Efficiency Characteristic
The description of the losses in the generator is done by a speed dependent torque
loss (mechanical losses) and a current dependent efficiency (electrical losses).
The nominal voltage is the threshold voltage for the regulator of the generator. If
the voltage in the onboard network is below this threshold the generator is acti-
vated. If the voltage is above the threshold the generator is deactivated.
The equivalent resistance regulator defines the gradient of the regulator voltage
(see Fig. 2.15).
UG,reg
UG,Nom
IG
Idle Voltage
The idle voltage is the voltage the generator is producing without an electrical
consumer, means that there is no flow of current. This idle voltage is a function of
the speed.
Efficiency Map
The efficiency map describes all losses of the generator (mechanical and elec-
trical). It depends on speed and current.
Torque Loss
M L;pd;gen u_ L;gen Torque Loss of the generator as a function of the generator speed Nm
The moment of loss is the mechanical loss in the generator due to friction at the
roller bearings and at the collector ring, aerodynamic friction at the fan, and, above
all, the power required to run the fan itself, which increases dramatically at higher
speeds.
Together with the torque that is transferred into electrical energy the moment of
loss builds the absorbed torque of the whole generator. The moment of loss depends
on the generator speed.
Efficiency
The efficiency is used to consider the copper and the iron losses. The copper
losses are produced by resistance in the rotor and stator windings. Their extent is
proportional to the power-to-weight ratio, i.e., the ratio of generated electrical
power to the mass of the effective components.
The iron losses result from the hysteresis and eddy currents produced by the
alternating magnetic fields in the iron in the stator and rotor. This efficiency is a
function of the delivered current.
Maximum Current
I L;max u_ L;gen Maximum Generator Current as a function of the generator speed A
The maximum generator current is the maximum current the generator can
deliver. The characteristic depends on the kind of generator (AC, DC).
ZL Switch –
1
CL;cont ¼ RG;equ ð2:10:2Þ
IL;max;cont
The current flow depends on the net and the controller minimum and maximum
voltages as well as the controller factor [2]:
UL;net UL;min;cont ½1 UL;net UL;max;cont
IL ¼ ð2:10:4Þ
CL;cont UL;net UL;max;cont
In the Generator Characteristic, the current flow is fixed by the net and the con-
troller minimum and maximum voltages as well as the maximum available gen-
erator current [2]:
sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi!2
IL;max;act
IL ¼ IL;max;act RG;equ UL;net UL;max;act þ UL;net UL;min;act
RG;equ
1 UL;net UL;max;act
ð2:10:5Þ
sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi!
IL;max;act 1
b ¼ 2 UL;max;act þ RG;equ ð2:10:8Þ
RG;equ CL;cont
sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi!2
IL;max;act UL;max;cont
c ¼ IC;max;act RG;equ UL;max;act þ þ ð2:10:9Þ
RG;equ CL;cont
Efficiency
The generator efficiency ƞL,act for the actual loading conditions is evaluated out of
the efficiency curve for the actual generator current IL.
Torque
IL UL;net
ML;in ¼ ML;pt;act ð2:10:10Þ
ZL u_ L;in
1 IL UL;net
ML;in ¼ ð2:10:11Þ
gL;act gL;el u_ L;in
Electric consumers are represented as ohmic resistors in the onboard network. They
represent an electric current loss. The number of the resistors that can be defined is
user-dependent. The resistors can be fixed by a constant value or by means of
2.11 Electrical Consumer (X) 123
2.11.1 Properties
Switch Output
If this switch is activated a result output for this component is made.
Selection Button Definition
Resistance
Current (Power)
The Resistance is defined depending on the value of the Data Bus input ‘Set Value X’.
Selection Button Definition
Resistance
Current (Power)
In this case, the Current or Power has to be defined depending on the values of the
Data Bus inputs ‘Set Value X’ and ‘Set Value Y’.
The nominal voltage is the definition of the voltage in the onboard network.
Resistance Curve
The Resistance Curve has to be defined if the selection button ‘Definition’ is set
to Resistance.
Switch: This is used to switch the resistor on and off depending on exceeding a
definable value of the switch. For this connected switch you can define a threshold
(Ux, border) which is the border for the on and off switching.
For on and off switching there are different possibilities:
Direction (positive/negative): Positive means that the switch is turned on when
the input of the switch is above the threshold value and it is turned off when the
input of the switch is below the threshold value. If negative is chosen it is vice
versa.
Reference (absolute/relative): The resistance is defined dependent on an input
value. If absolute is chosen the resistance is switched on and off when the switch is
turned on and off. If relative is chosen the resistance starts every time the switch is
turned on at zero.
Exceeding of value range (admissible/inadmissible): Sometimes it can occur that
the switch is longer open than the resistance is defined or that the set value is
outside the defined range (e.g. the resistance is defined as a function of the engine
speed with a maximum of 6000 rpm and the engine speed reaches 6500 rpm, i.e.,
that the resistance is outside the defined value range. If the exceeding of value range
is inadmissible the switch will be turned off when the defined border is reached. If
exceeding of value range is admissible the resistance will be extrapolated to get data
for the resistance.
Current (Power) Map
The Current (Power) map has to be defined if the selection button ‘Definition’ is
set to Current (Power).
The Current or Power is defined depending on the values of the Data Bus Inputs
‘Set Value X’ and ‘Set Value Y’.
Additionally, the map interpolation mode can be selected between ‘Continuous’
and ‘Steplike’. For further description of these modes, please refer to the description
of component ‘Map’.
SX,act Switch –
ZX,act Set Value X –
ZY,act Set Value Y –
IX Current A
UX,net Net voltage V
Px Power W
UX;net
IX ¼ ð2:11:1Þ
RX;act
The electric motor can be used for electrically driven vehicles, cars with hybrid
drive or to operate auxiliaries such as a fan or an oil pump.
The electric motor is defined by means of characteristic curves. Therefore, dif-
ferent motor type models can be constructed [1].
126 2 Mathematics Behind the Models
2.12.1 Properties
Switch Output
If this switch is activated a result output for this component is made.
Selection Button Definition
Efficiency
Efficiency Map
Load-dependent Efficiency Maps
Current Map
The efficiency is defined as a constant value. In addition to the torque map, the
efficiency-current curves (idle current map, starting current map) have to be given.
Selection Button Definition
Efficiency
Efficiency Map
Load-dependent Efficiency Maps
Current Map
Here the efficiency map has to be defined. Then the current–voltage characteristic is
calculated with the efficiency map and the torque-speed map.
Selection Button Definition
Efficiency
Efficiency Map
Load-dependent Efficiency Maps
Current Map
In addition to the torque map and the efficiency–current curves, there can be up to 5
load-dependent efficiency maps defined.
The efficiency value of the actual voltage is interpolated from this data.
Selection Button Definition
Efficiency
Efficiency Map
Load-dependent Efficiency Maps
Current Map
The Motor Characteristic is defined by the torque map and the current–voltage map.
5. Switches: Load-Dependent Efficiency Map U1, …, U5
If the selection button Definition is set to Load-dependent Efficiency Maps, up to
5 efficiency maps can be defined for up to 5 voltage levels.
The minimum speed is the minimum angular velocity the electric motor can run at.
The maximum speed is the maximum angular velocity the electric motor can run at.
The nominal voltage is defined for every electric motor. Normally it should be
the same as in the whole onboard network.
ηJ Efficiency %
This input can only be provided if the selection button Definition is set to
Efficiency.
M J u_ J;out ; UJ Torque–Voltage Map of Motor Nm
Map of the output torque dependent on the speed and the voltage.
I J u_ J;out ; U J Current–Voltage Map of Motor A
Map of the absorbed current dependent on the speed and the voltage.
gJ u_ J;out ; U J Efficiency Map –
Efficiency map is dependent on the speed and the voltage. It has to be defined
when selected in the property window. In addition the starting current value has to
be given.
gJ;i u_ J;out;i ; PJ;i Load-dependent Efficiency Maps i = 1, …,5 –
For each voltage level selected by the Load-dependent Efficiency Maps switches,
efficiency map depending on the speed and the electrical power can be defined.
I J u_ J;out Idle Current Curve A
I J;A hJ;out Starting Current Curve A
The last 2 curves must be defined if the selection button Definition is set to
Efficiency or Load-dependent Efficiency Maps.
αJ Load Signal –
ZJ Switch –
For the actual angular velocity nact and the actual voltage the actual values of MJ are
calculated:
n, nact... (actual) speed of the motor [2]:
For the actual angular velocity nact and the actual net voltage, the actual value of
the mechanical power is determined from the torque map M (nact, UJ,act) [2]:
Current–Voltage Map
Calculation by efficiency value, idle current map, and starting current map. Actual
electrical power [2]:
130 2 Mathematics Behind the Models
PJ;mech;act
PJ;el;act ¼ n 6¼ 0; M 6¼ 0 ð2:12:4Þ
gJ
PJ;el;act
IJ;act ¼ : ð2:12:5Þ
UJ;act
Special cases:
n = 0, M = 0 → Iact = 0.
n ≠ 0, M = 0 → by the idle current curve: IJ,act = I(nact).
n = 0, M ≠ 0 → by the starting current curve: IJ,act = I(MJ,act).
Efficiency Map
Calculation of the actual power and the actual current by the efficiency map and the
starting current [2]:
PJ;mech;act
PJ;el;act ¼ ð2:12:6Þ
gJ;act
PJ;el;ct
IJ;act ¼ ð2:12:7Þ
UJ;act
PJ;mech;ct
PJ;el;act ¼ ð2:12:9Þ
gJ;act
PJ;el;act
IJ;act ¼ ð2:12:10Þ
UJ;act
Special cases:
n = 0, M = 0 → IJ,act = 0.
n = 0, M ≠ 0 → by the starting current curve: IJ,act = I(MJ,act).
n ≠ 0, M = 0 → by the idle current curve: IJ,act = I(nact).
Minimal angular velocity u_ min maximal angular velocity u_ max ... are given by the
user; Minimal
Mmin, maximal torque Mmax … are determined by the torque
torque
map MJ u_ J;out ; UJ
If the actual values of angular velocity/torque are beyond their 1 %-limit of their
maximum values, the engine is in starting/idle status.
Starting limit u_ 0 ¼ u_ max 0:01
When u_ act ¼ u_ 0 we consider the motor to be in standstill, and the calculation
uses the starting current curve.
Idle limit M0 = Mmax · 0.01
When Mact < M0, we consider the motor to be in idle status, and the calculation
uses the idle current curve.
dependent on the actual temperature of the motor and the maximal moment of
inertia is determined according to these values [1].
2.13.1 Properties
Switch Output
If this switch is activated a result output for this component is made.
Selection Button Losses
Efficiency
Power Loss
The efficiency has to be defined depending on speed and power (or torque).
Selection Button Losses
Efficiency
Power Loss
In this case, the power loss has to be defined depending on speed and power (or
torque).
Selection Button Temperature
Calculated
From Data Bus
The machine temperature is calculated by the component using a special thermal
model. Selection button Temperature
Calculated
From Data Bus
The actual value of the machine temperature is taken from the Data Bus and used
for further calculations in the component.
Selection Button Current Limit
On
Off
If the button is switched on, the actual current is set to the user-defined limit every
time it exceeds this limit.
Selection Button Current Limit
On
Off
When the button is switched off, the actual current has no user-defined limit.
Selection Button Control Variable
Load Signal
Desired Torque
In this case the Electric Machine is controlled by Load Signal. The Data Bus Input
‘Load Signal’ has to be connected.
2.13 Electric Machine (EM) 133
The user can select between an ASM (Asynchronous Motor) and a PSM
(Permanent magnetic Synchronous Motor).
The maximum speed is the maximum angular velocity the electric motor can run at.
Thermal Model
This is the temperature which the electric machine has at calculation start.
The specific heat transition sums up all influences (such as material, surface state
etc.), which influence the transmission of heat.
This is the maximum time that the machine can run with maximum power.
The thermal time constant of maximum power is the time duration where the
maximum power can be delivered. During this time, the ambient temperature must
not increase.
The specific heat capacity defines the different heating behavior of the material.
It is the energy, which is needed in order to warm up material of 1 kg mass by 1 K
(temperature).
The layout temperature is the temperature where the electric machine has its
highest efficiency.
With this coefficient, the influence of the temperature on the remanence induc-
tion is defined.
2.13 Electric Machine (EM) 135
UEM Voltage V
PEM;UEM ;max u_ EM;out Maximum Power (Torque) mechanical kW
The Efficiency of the electric machine is defined dependent on speed and power
(or torque). This has to be defined if the selection button ‘Losses’ is set to
Efficiency.
PEM;loss PEM;UEM;act ;_uEM;out . Power Loss W
The Power Loss of the electric machine is defined dependent on speed and
power (or torque). This has to be defined if the selection button ‘Losses’ is set to
Power Loss.
€ EM;out
MEM;dt ¼ MEM HEM;nom u ð2:13:1Þ
2
MEM ¼ MEM;drag u_ EM =u_ EM;max ð2:13:3Þ
For permanent field machines the losses of iron have to be taken into
consideration.
If the drivetrain is turned on the following calculation is used [2].
The loss of power PEM,loss consists of losses of iron, copper, and losses
according to friction.
It is completely transformed into heat. The mechanical power of the drivetrain is
defined by [2]:
There are two possible cases for the loss of power and the efficiency in
motor-related mode. There MEM,1 and u_ EM;1 are the smallest data points in the
characteristic map of the machine.
Motor-related losses [2]
If M EM;out [ M EM;1 ; u_ EM;out [ u_ EM;1
PEM;loss;act ¼ PEM;mech;act 1=EM MEM;out ; u_ EM;out 1 ð2:13:7Þ
If jM EM j 0; u_ EM ffi 0
ð2:13:11Þ
If MEM;out _ EM;out ffi 0
0; ug
h i
PEM;loss ¼ u_ EM MEM;1 aðMEM Þ 1=g MEM;1 ; u_ EM;1 1 þ ð1 aðMEM ÞÞ 1 gEM MEM;1 ; u_ EM;1
EM
ð2:13:12Þ
where
aðMEM Þ ¼ MEM MEM;drag u_ EM 2 = MEM;1 MEM;drag u_ EM 2 ð2:13:13Þ
The difference between the losses in motor-related mode and the losses in
generator-related mode is given by the lowest torque in the Characteristic Map and
the drag torque. In the domain above the smallest positive torque the motor-related
mode is used and the efficiency is unique defined. In the domain below the negative
torque and below the drag torque the generator-related mode is used. In the domain
2.13 Electric Machine (EM) 139
between the described domains again the motor-related mode is being used. There a
smooth transition in the loss of power.
Furthermore the terminal current of the drivetrain is given by [2]:
For the maximal torque the following is defined using the loss of power and [2]:
Here the complete Characteristic Map is being searched beginning at the max-
imal permitted torque, until the loss of power is smaller than PEM,loss,max. This value
is the maximal permitted torque. The same holds for the generator-related mode.
There the maximal torque is calculated using the maximum of all angular
velocities using [2]:
The basic equations of the transition of heat in the electric machine are [2]:
PEM,W is the complete power, which is transformed into heat. It consists of the heat
power due to electric losses and due to the heat transition to the environment [2]:
and [2]:
PEM;env ¼ aEM;th TEM TEM;env ð2:13:23Þ
140 2 Mathematics Behind the Models
Together with the components supercapacitor and electric machine, the component
battery H is used for modeling hybrid vehicles.
The basic model consists of a voltage source and a resistance. The resistance is
constructed that way that complex processes within the battery can be considered.
2 RC elements can be added optionally. They describe the concentration over-
voltage and the transition overvoltage.
Also, the resistance’s dependence on the temperature can be activated optionally.
There can be modeled single cells as well as a combination of them. Therefore,
any modules can be constructed [1].
The thermal behavior of the battery is described by a thermal model. Here, the
warming caused by the losses and the cooling caused by convection is taken under
consideration [1].
2.14.1 Properties
Switch Output
If this switch is activated a result output for this component is made.
Selection Button Temperature
Calculated
From Data Bus
The battery’s temperature is calculated by the component using a thermal model.
Selection Button Temperature
Calculated
From Data Bus
The actual value of the battery’s temperature which has been calculated externally
is taken from the Data Bus and used for further calculations in the component.
Selection Button Ohmic Resistances
Constant
Temperature Dependent
Temperature and SOC Dependent
The constant values are used for the calculation.
Selection Button Ohmic Resistances
Constant
Temperature Dependent
Temperature and SOC Dependent
In this case the temperature-dependent maps for the internal resistances are used.
2.14 Battery H (QH) 141
The initial charge is the charge the cell has at the beginning of the calculation. It
is defined as percentage of the maximum charge.
The actual voltage of the cell should not exceed the defined value of the
maximum voltage.
The actual voltage of the cell should not be lower than the defined value of the
minimum voltage.
Number of Cells
The battery’s cell structure is kind of a matrix with a constant number of rows
and a constant number of cells in each row.
Thermal Model
This is the temperature where the battery has its highest efficiency.
The specific heat transition summarizes all influences (such as material, surface
state, conductivity) which influence the transmission of heat.
The specific heat capacity defines the different heating behaviors of the
material. It is the energy which is needed in order to warm up material of 1 kg mass
by 1 K.
Characteristic Maps of Machine
The idle voltage—charger is the idle voltage of one cell in charge mode. There is
no electrical consumer. It is a function of the SOC.
2.14 Battery H (QH) 143
TG,i Temperature i = 1, …, 5 °C
Through an internal calculation it is possible to use one single cell as well as any
combination of cells. Therefore the user can construct any desired module.
Here the current influences the idle voltage and the resistance only by its sign.
The maximal currents of the battery cell are calculated by [2]:
IQH;max;charge ¼ ðUQH;idle;charge IQH;terminal ; SOC
UQH;max
QQH;conc =CQH;conc
QQH;trans =CQH;trans Þ=RQH IQH;terminal :
ð2:14:4Þ
and [2]:
IQH;max;discharge ¼ ðUQH;idle;discharge IQH;terminal ; SOC UQH;min;discharge
QQH;conc =CQH;conc QQH;trans =CQH;trans Þ=RQH IQH;terminal
ð2:14:5Þ
with [2]:
The maximum power in and maximum power out are calculated in a numerical
way by searching for the maximum within the functions [2]:
PQH;out;max UQH;min;discharge ; IQH;terminal and PQH;in;max UQH;max;charge ; I QH;terminal
148 2 Mathematics Behind the Models
Here, PQH,th is the total power that is converted into heat. It consists of the heat
power inside the cell caused by electrical losses and of the heat transfer to the
environment [2]:
and it is [2]:
PQH;th;ambient ¼ aQH;th;trans TQH;bat TQH;ambient ð2:14:10Þ
and [2]:
PQH;th;el ¼ IQH;ohmic
2
RQH þ IQH;ohmic;trans
2
RQH;trans þ IQH;conc
2
TQH;act RQH;conc
þ 0:5IQH;terminal UQH;idle IQH;terminal SOCQH UQH;idle IQH;terminal ; SOCQH
ð2:14:11Þ
with the last term describing the losses caused by the polarization voltage.
2.15.1 Properties
Switch Output
If this switch is activated a result output for this component is made.
2.15 Gearbox Control (O) 149
The velocities for up- and downshifting can be read e.g. in the following way:
The upshifting velocity of the 2nd gear means that at this velocity the gearbox
control is upshifting from the 2nd into the 3rd gear. The downshifting velocity for
the 2nd gear means that at this velocity the gearbox control is downshifting from
the 3rd into the 2nd gear. Therefore it is necessary to define the up- and down-
shifting velocities always only for one gear less than are available in the gearbox
(i.e., for a five step gearbox, only for four gears the up- and downshifting velocities
have to be defined).
Gear Shifting according to Speed
u_ O;incr ½i Upshifting Speed u_ O;incr ½i for the nth gear rpm
u_ O;decr ½i Downshifting Speed u_ O;decr ½i for the nth gear rpm
Similar to the shifting according to velocity the speeds for up- and downshifting
can be read e.g. in that way:
The upshifting speed of the 2nd gear means that at this speed the gearbox control
is upshifting from the 2nd into the 3rd gear. The downshifting speed for the 2nd
gear means that at this speed the gearbox control is downshifting from the 3rd into
the 2nd gear.
Therefore it is necessary to define the up- and downshifting speeds always only
for one gear less than are available in the gearbox (i.e., for a five step gearbox, only
for four gears the up- and downshifting speeds have to be defined).
Gear Shifting according to Speed of next Gear
u_ O;incr;next ½i Upshifting Speed of next Gear for the nth gear rpm
u_ O;decr;next ½i Downshifting Speed of next Gear for the nth gear rpm
Here the upshifting and downshifting speeds of the next gear are defined sep-
arate for each gear.
Delay Time Gear Dependent
Here the upshifting and downshifting time delays are defined separately for each
gear. Shifting procedures which are not described in the table get the time delay 0.
2.15 Gearbox Control (O) 151
The gear shifting strategy up (Zcont,up) and down (Zcont,dn) has the information from
where the gearbox control should take the target gear:
• profile shifting
• velocity shifting
• speed shifting
• program shifting
152 2 Mathematics Behind the Models
The target gear position comes from the driving profile [2]:
If jG;target jG;act ! jG;new;up ¼ jG;target ð2:15:1Þ
If jG;target jG;act ! jG;new;dn ¼ jG;target ð2:15:2Þ
Each angular velocity could be selected via Data Bus for the reference speed [2].
If jG;act \ NG ^ u_ ref;act u_ o;incr bjG;act b ! jG;new;up ¼ jG;act þ 1 ð2:15:5Þ
If jG;act [ 1 ^ u_ ref;act u_ o;deu bjG;act c ! jG;new;dn ¼ jG;act 1 ð2:15:6Þ
2.16.1 Properties
Switch Output
A result output for this component is made if this switch is activated.
Switch Optimized Shifting Program
By activating this switch the optimized shifting table is used in the gearbox
program. (An optimized shifting table must be established using “Start
Optimization of Shifting Program”.)
Selection Button Shifting Program Selection
Default
Optimized
Advanced
In this case up to five shifting programs can be defined.
Selection Button Shifting Program Selection
Default
Optimized
Advanced
In this case a shifting program is activated which can be optimized with iSIGHT.
Selection Button Shifting Program Selection
Default
Optimized
Advanced
In this case up to five shifting programs and up to five kick down tables can be
defined.
Remark: the advanced shifting programs import gear shifting data generated by the
task ‘GSP Generation’.
154 2 Mathematics Behind the Models
The gear shifting process is done similar to the gearbox control. The only
difference is the dependence on the load signal of the engine.
If the shifting program should be defined not by speed but by velocity, the Data
Bus input ‘Velocity’ must be connected. The units of the shifting program and the
optimized shifting program influence each other. That means, if one of these pro-
grams is loaded with a different unit, also the other program gets this new unit.
The speeds for up- and downshifting can be read e.g. in the following way:
The upshifting speed of the 2nd gear with the special load signal means that at
this speed with this load signal the gearbox control is upshifting from the 2nd into
the 3rd gear.
The downshifting speed for the 2nd gear with the special load signal means that
at this speed with this load signal the gearbox control is downshifting from the 3rd
into the 2nd gear. Because it is not possible to define the shifting speeds for every
possible load signal the characteristic will be interpolated between the given points.
It is only necessary to define the up- and downshifting speeds always for one
gear less than are available in the gearbox (i.e., for a five step gearbox, only for four
gears the up- and downshifting speeds have to be defined).
2.16 Gearbox Program (P) 155
u_ P;incr n; aE;dk or Upshifting Speed 2–5 (Velocity) dependent on the gear rpm or
vP;incr n; aE;dk and the load signal km/h
u_ P;decr n; aE;dk or Downshifting Speed 2–5 (Velocity) dependent on the rpm or
vP;decr n; aE;dk gear and the load signal km/h
Regarding the minimization of the fuel consumption and the minimization of the
Full Load Acceleration time in standard driving cycles, it is possible to optimize the
Gear Shifting Program. Furthermore the driving performance should be improved
and the shifting process which should be optimized is equal to the normal gear
shifting process (see Fig. 2.16).
Figure 2.16 shows a typical gear shifting map. The green line is for downshif-
ting, the red one for upshifting of the gear. The first lines are the change between
the gears one and two, second pair between two and three and so on. This means, if
the actual speed is lower than the most left green line, then the gear is changed from
the second to the first gear, and if it is larger than the first red one, then the gear is
changed from the first to the second gear.
2 to 1 1 to 2
100
90
80 1 Up
70
Load signal [%]
2 Up
60
3 Up
50
4 Up
40
30 1 Down
20 2 Down
10 3 Down
0
4 Down
0 1000 2000 3000 4000 5000 6000
Gear output speed [1/min]
Fig. 2.16 Gear shifting map example. (dashed red down shifting, black up shifting) [2]
156 2 Mathematics Behind the Models
100
90
80
Load signal [%] 70
60
50
40
30
20
10
0
0 1000 2000 3000 4000 5000 6000
Gear output speed [1/min]
The discretized gear shifting map has between 80 and 100 unknowns, depending
on the refinement level of the gear shifting map in the input data.
The following figure denotes the possible movements of the unknowns (dis-
cretization node). The horizontal lines can move to the right and left and the inner
points may move up and down.
The limits for the movements are for the inner points the neighbor points of the
same line and for the horizontal lines the corresponding other line (the up for the
down line and vice versa) such that there is still a hysteresis in the curve.
This is needed for the decision of the Gear Shifting Program about up or
downshifting for an actual speed of the vehicle (Fig. 2.17).
The fuel consumption is given by the calculation of the driving cycle in liters per
100 km as calculated in the task for the cycle. To the optimization objective
different penalty terms are added to introduce some regularization of the perfor-
mance and the result.
The objective function for the fuel optimization is [2]:
Z Z Z
J ¼ afc þ b1 ðvactual vdesired Þ2 dt þ b2 jgear0 jdt þ b3 jacceleration0 jdt:
time time time
ð2:16:1Þ
file system
data check adaptation
data
AVL
data adaptation
CRUISE iSIGHT
GUI preparation optimization
kernel calls
iSIGHT call
results
x 2 IRn : gj ð xÞ ¼ 0; j ¼ 1; . . .; me ð2:16:2Þ
gj ðxÞ 0; j ¼ me þ 1; . . .; m
x l x xu
X
m
Lð x; uÞ :¼ f ð xÞ uj gj ðx): ð2:16:3Þ
j¼1
quadratically, where x ϵ IRn, and where u = (u1,…, um)T ϵ IRm is the multiplier
vector.
To formulate the quadratic subproblem, it will be proceeded from given iterates
xk ϵ IRn, an approximation of the solution, vk ϵ IRm and approximation of the
multipliers, and Bk ϵ IRnxn an approximation of the Hessian of the Lagrangian
function.
Then the following quadratic programming problem has to be solved [2]:
2.16 Gearbox Program (P) 159
1
min d T Bk d þ rf ðxk ÞT d ð2:16:4Þ
2
xkþ1 xk dk
:= þ ak ð2:16:5Þ
vkþ1 vk uk vk
wk :¼ xkþ1 xk ð2:16:2Þ
Y
Bkþ1 :¼ ðBk ; qk ; wk Þ ð2:16:7Þ
Y qqT BwwT B
ðB; q; wÞ :¼ B þ T T ð2:16:8Þ
q w w Bw
The above formula yields a positive definite matrix Bk+1 provided that Bk is
positive definite and qTk wk [ 0.
As already stated, the given optimization problem is not smooth. To avoid the
direct calculation of the gradients the results calculated in the line search are
approximated with a quadratic function. Since this function is differentiable, the
gradients used in the SQP method could be calculated. The optimization of the gear
shifting map does not include any constraints, but some boundaries for the
parameter vector (discretization nodes) which are chosen according to the possible
movement of the nodes.
The optimized shifting program can be defined and optimized in the unit for
speed as well as in the unit for velocity. The units of the shifting program and the
optimized shifting program influence each other. That means, if one of these pro-
grams is loaded with a different unit, also the other program gets this new unit.
160 2 Mathematics Behind the Models
If the optimized shifting program is defined in the unit for velocity, the Data Bus
input ‘Velocity’ must be connected.
Kickdown Table
The Kickdown Table is used when the driver steps full on the acceleration pedal
(load signal = 100 %). Instead of this criterion, the optional Data Bus input
‘Kickdown Signal’ can be used to control the kickdown event.
Velocities can be defined instead of speeds. In this case the Data Bus input
‘Velocity’ has to be connected.
Shifting Programs Advanced
Shifting Programs 1–5
Up to five shifting programs can be defined. With the Data Bus input ‘Shifting
Program Selector’ (which provides a double value), interpolations between the
shifting programs can be done.
Kickdown Tables 1–5
When in the properties window the switch ‘Shift Remaining Time’ is activated,
the remaining time for upshifting and downshifting for every gear can be defined.
The remaining time is the minimal time after the shift decision, in which no further
shifting is allowed.
Exception: a kickdown situation. In this case there is no remaining time.
These channels with the actual speeds for up and downshifting are available for
user-defined controls (in function component, Black Box, MATLAB® API,
MATLAB® DLL).
For advanced control algorithms, for current load signal and gear, the gearbox
program calculates upper and lower speed limits for shifting.
In order to ensure the existence of upper and lower speed limits, in the lowest gear,
lower speed limit is calculated as (current speed—INFINITESPEED), and in the
highest gear, upper speed limit is calculated as (current speed + INFINITESPEED).
INFINITESPEED is defined in the Cruise.ini file with the default value of
104.72 rad/s which equals 1000 rpm.
162 2 Mathematics Behind the Models
NG Number of gears –
Gear shifting to the next or next lowest gear will be done according to the gear
shifting criteria by using the defined maps.
The CVT Control is an easy way to pretend a target transmission considering the
velocity and the engine load. In addition there is a possibility of controlling a
clutch. This clutch is simply added to the model but it is not necessary [1].
2.17.1 Properties
Switch Output
If this switch is activated a result output for this component is made.
Selection Button Control Value
Transmission Ratio
Input Speed
Input Speed (ECO)
In this case the CVT is controlled by definition of Desired Transmission Ratio
depending on velocity and load signal.
Selection Button Control Value
Transmission Ratio
Input Speed
Input Speed (ECO)
In this case the CVT is controlled by definition of Desired Input Speed depending
on velocity and load signal. With this information the actual transmission ratio is
then calculated. When using this option, the Data Bus input channel ‘Input Speed’
has to be connected.
2.17 CVT Control (H) 163
Here the time constant of the PT1-retarding element for the load signal can be
defined.
Time Constant Velocity
Here the time constant of the PT1-retarding element for the velocity can be
defined.
Desired Transmission Ratio
iH,soll (vV, αE,dk) Desired Transmission Ratio as a function of velocity and load signal –
Here the demand transmission of the CVT gearbox can be defined as a function
of the velocity and the load signal.
Desired Input Speed
u_ H;soll vV ; aE;dk Desired Input Speed as a function of velocity and load signal rad/s
Here the demand input speed of the CVT gearbox can be defined as a function of
the velocity and the load signal.
To use this option, the Selection button ‘Control Value’ has to be set to ‘Input
Speed’ and the Data Bus input channel ‘Output Speed’ has to be connected.
Desired Input Speed (ECO)
u_ H;soll vV ; aE;dk Desired Input Speed as a function of velocity and load signal rad/s
164 2 Mathematics Behind the Models
This map works just like the ‘Desired Input Speed’ map. To use it, the selection
button ‘Control Value’ has to be set to ‘Input Speed (ECO)’ and the Data Bus input
channel ‘Output Speed’ has to be connected.
Desired Clutch Release
_ aE;dk or
LH;soll u; Desired Clutch Release as a function of speed (or %
LH;soll vV ; aE;dk velocity) and load signal
The demand clutch release of the starting clutch can be defined as a function of
the speed (or velocity) and the load signal.
Fast Linear Interpolation for CVT Control Maps
The interpolation in the maps of CVT control (e.g. interpolation of transmission
ratio) is done by an alternative method which does not use the map preprocessing
approximation to equidistant points.
This method is the default setting and can be changed in the ‘Interpolation and
Approximation’ section in the ‘Expert Mode’ settings of the Project Settings
(switch ‘Fast Linear Interpolation for CVT-Control-Maps’).
The demanded transmission step is selected from the map for the current vehicle
velocity and the throttle position.
The same thing is done for the new clutch release. Here, the reference speed and
the throttle position are taken into consideration.
The Anti-Slip Control checks the Force Transmission Factor (ratio between force
that has to be transmitted and maximum transmittable force) of all connected
wheels (driven wheels). If the force that has to be transmitted exceeds the adhesion
limit, the load signal or the clutch release changes [1].
The Anti-Slip Control is only active when a model is calculated quasi-stationary
in combination with special drive train configurations [1].
2.18.1 Properties
Switch Output
A result output for this component is made if this switch is activated.
Switch Use for Simulation
If the ASC should be used for calculations in simulation mode this switch can be
activated. In this case the parameters for the PID controller can be set.
PID Controller
Here the parameters for the PID controller can be defined. This controller is used
for the controlling of the ASC.
166 2 Mathematics Behind the Models
At first, the wheel with the highest value for the load transmitting factor is selected
(Fig. 2.19). If this value greater than 1, the Anti-Slip Control is activated (ZASC = −1)
and the load position will be reduced as long as the wheels have slip conditions.
For the task Full Load Acceleration from rest, the value for the clutch release
will be increased.
For the modification of the signals a PID controller is used. The parameter values
have to be entered in the input dialog.
2.19.1 Properties
Switch Output
If this switch is activated a result output for this component is made.
Selection Button PID Mode
Standard PID Control
Advanced Signal Control With Limitation
The functionality of a standard PID control is supported.
Selection Button PID Mode
Standard PID Control
Advanced Signal Control With Limitation
The second operation mode of the PID control is calculated.
Selection Button Actual Value Limitation
Lower Limit
Upper Limit
The actual value is limited from the lower side. This option is only available in the
second operation mode of the PID control.
Selection Button Actual Value Limitation
Lower Limit
Upper Limit
The actual value is limited from the upper side. This option is only available in the
second operation mode of the PID control.
168 2 Mathematics Behind the Models
CPID,act,min Minimum –
If the selection button Actual Value Limitation is set to Lower Limit, the
possible minimum of the actual value has to be defined.
CPID,act,max Maximum –
If the selection button Actual Value Limitation is set to Upper Limit, the
possible minimum of the actual value has to be defined.
Fixed Desired Value
If the switch Fixed Desired Value is chosen, the desired constant has to be
defined.
Output Value Limitation
CPID,out,min Minimum –
CPID,out,max Maximum –
If the switch Output Value Limitation is chosen, the possible minimum and
maximum of the output value are defined here.
2.19 PID Control (PID) 169
The output value is sent to the Data Bus and can be used by other components.
The difference between target value and actual value is calculated by [2]:
If there is a limitation of the actual value by definition of a lower limit, the sign
of ΔC is changed.
The controller output is defined by [2]:
Zt dðDC Þ
CPID;out ¼ CPID;P DC þ CPID;I ðDCÞðsÞds þ CPID;D ð2:19:2Þ
0
dt
In the 2nd mode (‘Advanced Signal Control With Limitation’), the entire con-
troller output is calculated by [2]:
If the PID Control is in 2nd mode (‘Advanced Signal Control With Limitation’),
ΔC < 0 and CPID,out = 0, the control is set to ‘neutral’ position, that means ΔC,
ΔCold and the
Integral are set to 0.
2.20.1 Properties
Switch Output
A result output for this component is made if this switch is activated.
Switch Dynamic Mode
This switch is activated by default. The brake calculation always stays dynamic and
never switches to a kinematic connection. Deactivation might be done for modeling
a sticking brake in kinematic mode.
Selection Button Control Variable
Brake Pressure
Braking Torque
2.20 Brake (B) 171
In this case the brake is controlled by the brake pressure which has to be supplied
via Data Bus input ‘Brake Pressure’.
Selection Button Control Variable
Brake Pressure
Braking Torque
In this case the engine is controlled by the desired braking torque value which has
to come from the connected Data Bus input ‘Braking Torque.’
The brake piston surface is the area of the hydraulic cylinder. Multiplied with the
brake pressure and the efficiency, it gives the axial brake force.
µB Friction Coefficient –
The friction coefficient is between the brake drum, respectively the friction disc
and the brake shoes.
The specific brake factor is a factor that depends on the design of the brake.
For disc brakes it is always one, for drum brakes it is usually larger than 1.
The effective friction radius is the radius where the braking force applies.
ƞB Efficiency –
The efficiency considers the effects of the conversion of the hydraulic into the
mechanical part of the brake (Fig. 2.20).
ΘB Inertia Moment kg m2
172 2 Mathematics Behind the Models
cB 2.7 cB 3.7
Drum Brake
(Duplex)
cB 1.7
Fig. 2.20 Specific brake factors for different brake types [2]
PB Brake Pressure Pa
MB Braking Torque Nm
2.20 Brake (B) 173
u_ B Speed rad/s
MB Braking Torque Nm
PB Power W
M B ¼ 2 p B A B g B lB r B c B ð2:20:1Þ
The cockpit links the driver and the vehicle. In this component, there are only
connections made via the Data Bus. On one hand, the driver receives information
such as the vehicle velocity and the vehicle acceleration. On the other hand,
information about of the driver, such as the pedal positions, are delivered to other
components. The pedal positions (e.g. clutch pedal position) are transferred into
corresponding indicators (e.g.. clutch release) via the pedal characteristics (e.g.
clutch pedal characteristic) [1].
2.21.1 Properties
Switch Output
A result output for this component is made if this switch is activated.
Switch Actual Distance and Velocity From Data Bus
The actual distance and velocity are taken from the Data Bus channels ‘Actual
Distance extern’ and ‘Actual Velocity extern’.
Switch Desired Velocity From Data Bus
When using AVL CRUISE as a dynamic linked library (CruiseNT.dll) with the
component ‘Interface’, the internal driver can control the velocity which is preset
over the Data Bus channel ‘Desired Velocity extern’.
It is important to connect the channel on the Data Bus system with an external
time depending signal. The driver will follow this desired velocity in an appropriate
range. Note that there must be an adapted profile. e.g.. gear shifting in the task
Cycle Run.
174 2 Mathematics Behind the Models
The switch is used for the Calculation Tasks, so that the calculation knows if the
gear shifting process should be controlled by the driver (manual) or the gearbox
control.
This is used to indicate to the driver how many gears the gearbox contains, i.e., if
there is still another gear for upshifting or not.
The maximum brake pedal force is needed for the calculation of the brake pedal
position. Output of the driver component is the brake pedal force. This brake pedal
force divided by the maximum brake pedal force is the brake pedal position.
It is also possible to define a maximum brake force in the component driver.
With this maximum brake force, the behavior of a driver can be considered. You
can have a driver who is stepping with a maximum of e.g. 100 N on the brake pedal
(defined in the driver component) and in the cockpit the maximum brake pedal
force is defined as 200 N. That means that the driver can have a maximum of 50 %
of the brake pedal position (100 N divided by 200 N). Thus, he has a lower brake
pressure and a longer brake distance.
This is used to give the driver the information how many retarder steps the
retarder is containing, i.e., if there is still another step with a higher braking torque
or not.
Acceleration Pedal Characteristic
pCO,brake (LBr) Brake Pressure as a function of the Specific Brake Pedal Force bar
2.21 Cockpit (CO) 175
The characteristics of the brake pedal, the clutch pedal, and the acceleration
pedal can be defined. It is possible to define an idle path or a special shape of the
characteristic.
u_ CO Speed rad/s
ZCO,cont,in (i) Operation Control 0–4 –
jCO,act Gear Indicator –
jCO,target Target Gear –
u_ CO;measure 1 Measuring Speed 1 rad/s
u_ CO;measure 2 Measuring Speed 2 rad/s
sextern Actual Distance extern m
vextern Actual Velocity extern m/s
vextern Desired Velocity extern m/s
In the module cockpit, switch, and selection functions were done. Input data were
transformed depending to the calculation settings and maps.
After that the changed values were put onto the output channels.
Exhaust systems consider the effects of the catalytic converter and soot filter on the
raw emissions of the engine. Starting with the temperature of the catalytic con-
verter, factors for the conversion will be computed for the single emissions com-
ponents on the basis of the temperature-dependent maps [1].
2.22.1 Properties
Switch Output
A result output for this component is made if this switch is activated.
Selection Button Temperature
Temperature calculated
Temperature from Data Bus
Calculation of temperature by using the thermal model
Selection Button Temperature
Temperature calculated
Temperature from Data Bus
The temperature is calculated by an external simulation tool and is transferred into
the catalyst component via Data Bus.
2.22 Exhaust System (EX) 177
The heat loss coefficient is for determining which amount of the exhaust energy
the engine has delivered reaches the exhaust system and which amount is lost
between the engine and the exhaust system to the environment.
The heat loss coefficient is the percentage of the exhaust energy delivered by the
engine which reaches the exhaust system.
This is the weight of the catalytic converter. This weight is needed for the
calculation of the temperature of the catalytic converter and as consequence the
efficiency of the emission conversion (Fig. 2.21).
This is the specific heat capacity of the catalytic converter which is needed for
the calculation of the temperature of the catalytic converter.
The reference area of the exhaust system is the outside surface area of the
catalytic converter. It is needed for the calculation of the heat transfer between the
catalytic converter and the ambient air.
This is the coefficient for the heat transfer between the catalytic converter and the
ambient air. With this heat transfer, the energy lost to the environment can be
determined.
The operating temperature is needed for a calculation with hot start. If the hot
start is chosen, the temperature of the catalytic converter at the beginning of the
calculation is equal to the operating temperature. If cold start is chosen, the tem-
perature of the catalytic converter at the beginning of the calculation is equal to the
ambient temperature.
Conversion Efficiency
This is the proportion of the heat of the exhaust gas which is transferred to the
exhaust system.
In relation to the entered exhaust energy PEX, the heat transmitting conditions, the
heat capacity of the exhaust system, and the temperature of the exhaust system TEX
are calculated.
In the next calculation step, the efficiency values ƞEX,act concerning TEX are
evaluated out of the maps.
The rest emissions are calculated with the following formula [2]
eEX;out ¼ eEX;...;in 1 gEX;...;act ð2:22:1Þ
2.23.1 Properties
Switch Output
A result output for this component is made if this switch is activated.
180 2 Mathematics Behind the Models
The input values are transferred to the MATLAB® function. From the module, the
results are received and input to the output channels.
Refer to the interfaces manual for more information.
The component Function can be used for calculating with user-defined functions.
There are two possibilities to define the function:
(a) Function definition in RPN-notation style (Reverse Polish Notation, UPN,
Postfix Notation, Stack Logic). No need for programming, the component has
a user-friendly selection table. Only one function output y can be defined.
(b) Programming in C code. Up to 99 function outputs can be defined.
Up to 99 input values can be read through the Data Bus. The values are the
function arguments, together with up to 10 constants defined by the user.
Those arguments must be connected to form the function term by using
operators which can be selected from a selection menu. Also the most
important physical constants (e.g., n) are in this menu, which can be selected
as additional function arguments. After the function value has been calculated,
it can be used on the Data Bus as an output value. When defining the function
in C code, it is possible to deliver values up to 99 function outputs [1].
2.24 Function (FU) 181
2.24.1 Properties
Switch Output
If this switch is activated a result output for this component is made.
Selection Button Formula Mode
RPN-Formula
C-Function
In this case the formula definition has to be done by the so-called Reverse Polish
Notation (UPN, postfix notation, stack logic).
Selection Button Formula Mode
RPN-Formula
C-Function
The function has to be defined in C code.
RPN-Formula Mode
This mode is active if the selection button ‘Formula Mode’ is set to RPN-formula.
The most important pre-defined constants can be selected from a selection menu
and used as function arguments.
In every time step this variable has the same value as the Data Bus output channel
‘Real Time’ of component cockpit. It can be used in System Analysis Mode
(SAM) in models without a cockpit.
182 2 Mathematics Behind the Models
The operators are selected from a selection menu. They form a function term
with the constants and input values from the Data Bus as arguments. The function
definition format is the so-called Reverse Polish Notation (UPN, postfix notation,
stack logic), i.e., the term “A + a[0] =” is written as “A a[0] + ENTER”. The
definition is done in a table and should be completed by typing Enter.
Click on the Function Preview button to open a window showing the function
term converted into the “normal” prefix notation. Additionally, there is a copy
button so that the user can copy this prefix term and paste it for example into the
comment line.
C-Function Mode
This mode is active if the selection button ‘Formula Mode’ is set to C-Function.
This is a pre-defined variable of type ‘double’ which can be used in the C code
without a declaration. In every time step this variable has the same value as the Data
Bus output channel ‘Real Time’ of component cockpit. It can be used in System
Analysis Mode (SAM) in models without a cockpit.
Additionally, most of the ANSI C operators and statements can be used, e.g., the
if-then statement.
Up to 99 Data Bus output variables y[0], y[1], …, y[98] can be delivered with
output values.
Up to 99 values can be read from the Data Bus and used as function arguments
(for RPN- mode only the first 20).
The function values are sent to the Data Bus and can be used by other com-
ponents. When defining the function in RPN-Formula mode, only the output y[0]
can be used.
The user-defined constants and the input values are transferred to the function
which calculates the output values by using the formula that has been defined by the
user in a table, or by using the user-defined C-function code. The output value(s)
can be used by other components through the Data Bus.
Example
The input value a[0] should be multiplied by 9.81. The result should be multiplied
by the factor 0.92 which has been defined as a constant (B).
In the parameter field the letter B is assigned the value 0.92. The function is
defined in the formula table as follows:
B:= 0.92
184 2 Mathematics Behind the Models
Position Y=
1 a[0]
2 9.81
3 *
4 B
5 *
6 ENTER
The expressions of the table are selected in a certain selection menu. At the end
of the table ENTER should always be written.
Derivation
This operation has 2 arguments. For example, the following input table means that
the signal of Data Bus input a[0] should be differentiated by the signal of Data Bus
input a[1].
Position y=
1 a[0]
2 a[1]
3 Derivation
4 ENTER
Integration
This operation has 3 arguments:
Position y=
1 a[0]
2 a[1]
3 a[2]
4 Integration
5 ENTER
The above input table means that the signal of Data Bus input a[0] should be
integrated by the signal of Data Bus input a[1], whereas Data Bus signal a[4]
represents the reset signal.
The meaning of the reset signal is as follows:
Value Action
>0 “normal” summation
0 Reset; integration output is set to 0 and stays 0 as long as the reset value is <=0
<0 No further summation
In the Constants component the user can define up to 99 constant values which can
be used by other components through the Data Bus. The values can be of type
integer, double or string [1].
2.25.1 Properties
Switch Output
If this switch is activated a result output for this component is made.
The constant output values are sent to the Data Bus and can be used by other
components.
186 2 Mathematics Behind the Models
2.26 Monitor
The monitor can be introduced if the calculation run must be detected. It is possible
to show some results of the calculation while the calculation is running. This is
done via the online-monitor [1].
It is possible to connect up to ten input channels of the monitor to output
channels of different components that are available at the AVL CRUISE Data Bus.
Currently, the output in the online-monitor is only available in the form of
tables [1].
2.26.1 Properties
Switch Output
A result output for this component is made if this switch is activated.
These input values can be read in via the Data Bus. During the calculation, they
were written on the online-monitor for each time step.
The wheels and tires link the vehicle to the road. The component wheel allows to
consider many influencing variables and their effect on the rolling state.
The moment of rolling drag can be computed on the basis of the wheel load, the
corrected dynamic rolling radius, and the coefficient of rolling drag [1].
The longitudinal tire force (circumferential force) results from the friction
coefficient, the wheel load as well as from the wheel load factor and the slip factor.
It is possible to define variable friction coefficients along the driving profile when
different road conditions are considered [1].
The wheel component also includes a detailed rolling resistance furnished by
Michelin. This model describes the rolling resistance of the tire depending on the
tire inflation pressure, load, translation speed, the ambient temperature and the time
[1].
2.27 Wheel/Tire (W) 187
2.27.1 Properties
Switch Variation
The Wheel can be given free for variation with this switch. The setup of the
variation parameters is done in the folder.
Switch Output
If this switch is activated a result output for this component is made.
Selection Button Dynamic Rolling Radius
Constant Value
From Characteristic
In this case the input constant value is used for the dynamic rolling radius.
Selection Button Dynamic Rolling Radius
Constant Value
From Characteristic
In this case the dynamic rolling radius is interpolated out of the characteristic.
Selection Button Wheel Slip
deactivated
Function
Function with Limit
Map
Map with Limit
In this case no consideration of the wheel slip in the calculation is possible.
Selection Button Wheel Slip
deactivated
Function
Function with Limit
Map
Map with Limit
Here the definition of the wheel slip is done by the characteristics “Maximum in the
Slip Characteristic” and “Asymptote in the slip characteristic.” If Function with
Limit is selected the maximum slip is limited with a slip which is slightly higher
than the maximum value from the slip characteristic.
This limitation is also used for Function without limitation for all stationary
calculations (Climbing Performance, Full Load Acceleration/Acceleration in all
Gears, Maximum Traction Force). This is done to avoid instabilities in the calcu-
lation and therefore results which can be misunderstood.
Selection Button Wheel Slip
deactivated
Function
Function with Limit
Map
Map with Limit
188 2 Mathematics Behind the Models
Here the specific traction (as percentage of the wheel load) is defined as function
of the wheel slip and the friction coefficient of the road. If Map with Limit is
selected the maximum slip is limited with a slip which is slightly higher than the
maximum value from the traction map.
This limitation is also used for Map without limitation for all stationary cal-
culations (Climbing Performance, Full Load Acceleration/Acceleration in all Gears,
Maximum Traction Force). This is done to avoid instabilities in the calculation and
therefore results which can be misunderstood.
Switch Rolling Resistance Model
If this switch is activated, the rolling resistance is calculated by a detailed resistance
model.
Switch Rolling Resistance
wheel load dependent
wheel pressure dependent
velocity dependent
temperature dependent
With these switches the input of the rolling resistance can be defined. It is possible
to activate more than one influence. All activated rolling resistances are added to the
overall rolling resistance. These switches can only be activated if the Switch
Rolling Resistance Model is deactivated.
Selection Button Dimensions
calculated
defined
In this case the tire diameter and the surface area are calculated with the input data
of tire section width, tire aspect ratio and tire seat diameter.
Selection Button Dimensions
calculated
defined
In this case the tire diameter and the surface area are defined by the user.
Selection Button Wheel Location
Vehicle: Front Left
Vehicle: Rear Left
Vehicle: Rear Right
Vehicle: Front Right
Trailer: Left
Trailer: Right
Two Axle Trailer: Front Left
Two Axle Trailer: Rear Left
Two Axle Trailer: Rear Right
Two Axle Trailer: Front Right
With this selection button the location can be defined where this wheel is located on
the vehicle or the trailer.
2.27 Wheel/Tire (W) 189
Wheel Slip
The friction coefficient of the tire depends on the material used in the tire. With
this friction coefficient it is possible to consider the effects of different wheels on the
behavior of the vehicle.
The curves for the maximum slip in the slip curve and the transmission ratio at
infinite slip are measured at the reference wheel load.
The calculation of the longitudinal force is done like follows [2]:
FL Longitudinal force
µRoad Friction Coefficient of the road
cS Slip Correction Factor
cF Wheel Load Correction Factor
FW,s Wheel Load
cW,s,F Wheel Load Correction Coefficient –
For the real wheel load is never the reference wheel load the longitudinal force
has to be corrected with the wheel load correction factor (cF). For this case the
correction coefficient is used. It determines the gradient of the curve (refer to
cW;S;F
Fig. 2.22; the gradient is FW;S;norm ÞAlthough the gradient is normally negative, a
positive correction coefficient has to be entered in the input field [2].
FW;act FW;S;norm
cF ¼ 1 cW;S;F ð2:27:2Þ
FW;S;norm
cF
FW,s,norm FW,s
Rolling Radius
Static Rolling Radius
The static rolling radius is the distance between the center of the wheel and the
road surface for the loaded vehicle but without moving.
The dynamic rolling radius is the distance between the center of the wheel and
the road surface for the loaded vehicle with moving. The dynamic rolling radius is
usually a little bit bigger than the static rolling radius.
Dynamic Rolling Radius (Constant)
The constant value is used when the selection button Dynamic rolling radius is
set to constant value.
Dynamic Rolling Radius (Characteristic)
Here the dynamic rolling radius can be defined as function of the velocity. This
characteristic is active when the selection button Dynamic rolling radius is set to
“from characteristic.”
Slip Curve
Maximum in the Slip Characteristic
cs
1
µ W,U
sW,grenz sw
The slip at maximum traction is together with the specific traction at infinite slip
the Characteristic Values for the slip correction factor (refer to Fig. 2.23). The slip
correction factor depends on the friction coefficient of the road.
Traction Map
This part of the rolling resistance is converted form the nominal wheel pressure
to the actual wheel pressure.
192 2 Mathematics Behind the Models
cW,r
vV
The Tire Section Width is the linear distance between the outside sidewalls of an
inflated tire without any load (exclusive of protruding side ribs and decorations).
The tire seat diameter is the diameter of the rim seat supporting the tire bead.
2.27 Wheel/Tire (W) 193
The tire design diameter is the diameter of the inflated tire without any load.
The surface area is the total external surface area of the tire.
Depending on the selection button dimensions, the tire section width, tire aspect
ratio and tire seat diameter, or the tire design diameter and tire surface area have to
be defined.
Measured data at ISO-test conditions:
Tire-specific data:
Initial conditions:
sW Slip Signal %
u_ W Speed rad/s
MW,trans Wheel Torque Nm
TW,act Tire Temperature K
FW,act Wheel Load Nm
PW Power W
PW,loss Power Loss W
FW,roll Rolling Resistance Force N
SW,peak,act Actual slip of peak friction point for current road condition –
µW,asym,act Actual Asymptote of slip correction factor for current road condition –
CW,coeff1 Wheel slip coefficient 1 –
CW,coeff2 Wheel slip coefficient 2 –
vW,ref Reference Velocity m/s
FW,a Longitudinal force of tires N
FW,roll Rolling drag N
rW,dyn Corrected dynamic rolling radius m
(continued)
2.27 Wheel/Tire (W) 195
There are two possibilities for determining the dynamic wheel radius.
1. Dynamic wheel base radius (ZW=1)
The dynamic wheel radius is constant for the whole speed range [2]:
The rolling resistance coefficient is evaluated out of the map Cw,f(Vv) for the
actual vehicle velocity.
The rolling resistance is defined as followed [2]:
c ¼ cv þ cq cRL þ cT þ cp ð2:27:7Þ
cv = cv(v)
cρ Rolling resistance influenced by side slip angle [2]:
mvehicle v2vehicle
cq ¼ abs sin af ðrÞ ð2:27:8Þ
qg
cT ¼ cT ðTU Þ ð2:27:10Þ
Slip Moment
Wheel Slip function
For the current road condition
) the Slip Characteristics can be evaluated
(µu,road
from the SW;peak lu;road and lw;asym lu;road
With this, we calculate the slip coefficients [2]:
2 arcsin lw;asxm;act
Cw;coeff1 ¼2 ð2:27:12Þ
p
1 p2
Cw;coeff2 ¼ tan ð2:27:13Þ
SW;peak;act 4 p arcsin lw;asym;act
Slip difference
The slip difference is defined as difference between the actual slip and the infinite
slip with the peak friction value [2].
DSW \0
ð2:27:18Þ
Mw;slip ¼ 0
DSW 0:03
ð2:27:19Þ
MW;slip ¼ SIGNðSW Þ 5 105 DS2W
DSW [ 0:03
ð2:27:20Þ
MW;slip ¼ SIGNðSW Þ 5 105 0:03 ð2 DSw 0:03Þ
For the actual wheel slip, the maximum transmittable force is [2]:
The moment of the vehicle (acting at the virtual wheel with r =1 m) in defined as
followed [2]:
FW;a
CW;act ¼ ð2:27:25Þ
FW;pot
FW;a
Cslip ¼ ð2:27:26Þ
FW;pot Cw;coeff1 Cw;coeff2
Wheel power
The wheel power is calculated out of the transmitted wheel torque and the wheel
speed [2]:
2RW WW
DW;design ¼ DW;seat þ ð2:27:29Þ
100
!
D2W;design D2W;seat
AW ¼ DW;design p WW þ 2p ð2:27:30Þ
4 4
The fundamental assumption of this detailed tire model is that the transient tire
temperature response can be modeled in the same way as the steady state response
to the ambient temperature [2]:
FW;roll ¼ FW;roll;stab 1 þ kW;emp TW;act TW;stab ð2:27:31Þ
In order to calculate the instantaneous rolling resistance at each time step, the
subsequent algorithm has to be followed:
2.27 Wheel/Tire (W) 199
(1) The steady state rolling resistance value FW,roll,stab at current conditions as a
function of tire load FW,act, tire inflation pressure zV,load, tire translation speed
vW,act and ambient temperature TU [2]:
!aW
pact FW;act bW
FW;roll;stab ¼ cW;roll;ISO FW;act
pW;ISO FW;ISO
" #
vW;act vW;ISO vW;act vW;ISO 2
1þb þc
vW;ISO vW;ISO
1 þ cW;ambient T W;ISO T U
ð2:27:32Þ
in which the coefficients cW,roll,ISO, pW,ISO, vW,ISO, TW,ISO, bW,lin, cW,quad, αW,
βW and cW,ambient are unique to each tire. The subscript ISO simply refers to
the parameters value under ISO (or actual) test conditions. Naturally these
values must be furnished along with the model coefficients.
In the above equation, we have not yet added the aerodynamic drag. Since the
aerodynamic drag is assumed independent of temperature, it is added after the
instantaneous rolling resistance is calculated.
(2) The stabilized tire temperature TW,stab is equal to [2]:
vW;act FW;roll;stab
TW;stab ¼ TU þ ð2:27:33Þ
CW;conv AW
where [2]:
c
vW;act
cW;conv ¼ cW;conv;ISO ð2:27:34Þ
vW;ISO
and cW,conv,ISO, γ and AW are unique to each tire. For the case of vW,act = 0 it is
necessary
To apply a free convection exchange coefficient. In this case, we first calculate
the Rayleigh number [2]:
g bW TW;act TU D3W;design
RaD ¼ ð2:27:35Þ
va
with [2]:
g gravitational constant (9.81 m/s2)
βW volumetric thermal expansion coefficient of air (which can be
approximated as TW;act2þTU ½K 1
v [m/s2] kinematic viscosity of air (evaluated at TW;act þTU through interpola-
2
tion in a table)
200 2 Mathematics Behind the Models
Then the Nusselt number and subsequently the heat transfer coefficient are
calculated [2]:
NuD ¼ 0:15RaD
0:333
ð2:27:36Þ
NuD kW;therm
cW;heat ¼ ð2:27:37Þ
DW;design
with
kW,therm [W/mK] thermal conductivity of air (evaluated at TW;act þTU
2 through
interpolation in a table)
(3) The instantaneous average tire temperature at each time increment is then
obtained by integrating [2]:
1 Zt
TW;act ðtÞ TW;act;last ¼ FW;roll vW;act cW;heat AW TW;act;last TU d~t
mW CW;heat;equiv 0
ð2:27:38Þ
with
cW,heat,equiv equivalent specific heat of the tire [J/kg K]
(4) Once the instantaneous tire temperature is calculated, the instantaneous rolling
resistance is calculated by [2]:
FW;roll ¼ FW;roll;stab 1 þ kW;emp TW;act TW;stab ð2:27:39Þ
References
1. AVL Cruise (2011) Users guide, AVL list GmbH, Graz, Austria, Document no. 04.0104.2011,
Edition 06.2011
2. AVL Cruise (2009) Theory manual, AVL list GmbH, Graz, Austria, Document no.
04.0108.2009, Edition 06.2009
References 201
3. Patton KJ, Nitschke RG, Heywood JB, Development and evalu-ation of a friction model for
spark-ignition engines. SAE Paper No. 89-0836
4. Shayler PJ, Christian SJ, Ma T, A model for the investigation of temperature heat flow and
friction characteristics during engine warm-up. SAE Paper 93-1153
5. Shayler PJ, Leong DKW, Murphy M (2003) Friction teardown data from motored engine tests
on light duty automotive diesel engines at low temperatures and speeds. In: ASME 2003, Fall
Technical Conference
Chapter 3
Virtual Powertrain Design
The drivetrain configuration is freely definable and the user can build in as many
clutches, torque converters, or elastic shafts as required.
The drivetrain configuration is like a standard vehicle with a manual gearbox. In
this configuration, only one clutch is allowed. In front and after the clutch, there can
only be one degree of freedom. In addition, just one gearbox is allowed in this
configuration. The standard manual gearbox was developed for shifting in the next
gear according to the engine speeds. In addition, check that the friction clutch is
directly connected to the gearbox.
The drivetrain configuration is like a standard vehicle with an automatic gear-
box. In this configuration, only one torque converter is allowed. In front and after
the torque, converter there can only be one degree of freedom. In addition, just one
gearbox or one Continuously Variable Transmission (CVT) is allowed in this
configuration. For the different configurations, different calculation algorithms are
used. Thus the calculation times can be decreased for standard drivetrains [1].
Standard automatic was developed for shifting according to speeds in the next
gear. In addition, the torque converter is directly connected to the gearbox or the
CVT.
In the User Area (Fig. 3.1) click on the Projects folder with the right mouse
button and select new project. Enter a name and select OK. The new name will
appear in the projects tree. Select this new project with the right mouse button and
select new version. Enter the new version name or keep the default version name.
Select OK. Select the version with the right mouse button and select save. All files
and subdirectories required for the project management are created [2].
In the navigation frame there is current default project directory with tree
structure, showing all existing subfolders, projects, and versions in that directory.
On the right side of the window, there are detailed data regarding objects in the tree
structure.
To start a new project, in User tab with right mouse click popup menu, select
new and chose project or version and create project or version with new name. For
already existing project, load project or version by selecting the desired version of
the project from project tree structure in User view [1].
After starting AVL CRUISE platform, select Desk tab and click to Vehicle
Model to open Desk Area to begin building the new model.
© Springer International Publishing Switzerland 2015 203
B.O. Varga et al., Electric and Plug-In Hybrid Vehicles,
Green Energy and Technology, DOI 10.1007/978-3-319-18639-9_3
204 3 Virtual Powertrain Design
To insert component into model click on the Desk Area (Fig. 3.2) and then click
on the Vehicle Model.
The components used to build the models are located under the Desk Area—
Modules Sidebar, and are described in Table 3.1.
In addition to the components of the drivetrain, whose single modules are
defined by a Configuration File at the program start, there are general modules on
each view which corresponds to a vehicle.
3 Virtual Powertrain Design 205
In the Vehicle Modeler, the components are organized in component groups and
they can simply be dragged and dropped onto the working area, where they can be
linked together to represent the desired drivetrain configuration.
To build a new model in AVL CRUISE, the recommended work process flow to
create a project/version is described in Fig. 3.3.
To connect the components, the user must click with the right mouse button on
connector of the first component to be connected and chose connect (or break the
established connections, or cloning the pin if are more of them needed). Drag the
connection with the mouse toward to desired component. Click with the left mouse
button on connector of the second component to be connected.
There are two types of connections: energetic connections (mechanical con-
nection—blue, exhaust connection—gray, electrical connection—red), and infor-
mational connection (red–green–blue).
Energetic connections which transfer energy are divided by the type of energy
which is transferred into the following groups (mechanical, exhaust, electrical) and
are established by polylines between the pins.
Information to the Data Bus is the only informational connection. These con-
nections are established by an arrow which shows the information flow. The
energetic connections are represented in form of lines on the screen [3].
The components with informational connections have only a small colored
arrow on the left bottom corner of the symbol. For connecting the different modules
with information flow, open the Data Bus dialog by double clicking with the left
mouse button on the colored lines (red—green—blue) at the bottom of the vehicle
modeler (Fig. 3.4).
The connections in the data bus are described as follows: in the first column
(component … requires) all components are listed which need information from
other components (e.g., IC Engine). In the second column (input information …
from), all values required by this component will be displayed (e.g., Throttle
Position). Values that are marked in blue in this column have to be connected;
values marked in black are optional. After selecting a value in the second column,
214 3 Virtual Powertrain Design
in the third column (component … delivering) all components appear which can
deliver such information (e.g., Cockpit; in this case the only component that can
deliver a throttle position). In the fourth column (output information …), the values
appear that can deliver this information (e.g., Throttle Position).
In AVL CRUISE a preselection is done so that only the components that can
provide the input information as requested in the second column are presented in
the third and fourth columns (in this case, only the Cockpit can deliver an infor-
mation like the throttle position). For this, the Cockpit delivers the information
about the actual throttle position to the engine, where it is used to determine the
output torque.
Click on a components Data Bus symbol (input or output) with the left mouse
visualizes the Data Bus connections between this component and other components
(Fig. 3.5). The channels appear as colored lines and also their names are written.
Freezing visualization is possible by clicking on the connection icon until other
connection is selected.
Components can be encoded using the button in the upper left corner of the
components window. The user can define a password with a maximum length of
eight characters. Encrypted components cannot be opened until the correct pass-
word is entered. If a model which contains encrypted components is calculated, all
results can be viewed except the results of the encrypted components.
Macrocomponents can also be encrypted using the right mouse button [1].
Components and their connections can be grouped into subsystems. The sub-
systems can be activated or deactivated so that several configurations of an active
subsystem can be defined. With these three subsystems it is possible to define two
3 Virtual Powertrain Design 215
vehicle configurations, which are from now on called systems. The structure of the
systems/subsystems and their activation status is shown i
• Conventional System, containing the subsystems:
– Main (activated);
– Manual Transmission (activated);
– Hybrid (deactivated).
• A Hybrid System, containing the active subsystems:
– Main (activated);
– Manual Transmission (deactivated);
– Hybrid (activated).
n the navigation tree and also in the “Systems” section of the sidebar of the vehicle
modeler.
Each component of the model has a popup menu accessible by right mouse
button click (Fig. 3.6).
Special components which have the sheet description of Data Bus (Function,
Map, Interface Components, etc.) have user-defined channel names.
The Kinematic Chain Browser can be launched in the popup menu of a selected
component (Fig. 3.7). All mechanical connections starting from the selected
component are shown in a tree structure in the left side of the window. When a
second component is selected in the tree structure, all speeds between the selected
components are shown in the table in the right side of the window, each row
representing one gear of the gearbox [1].
216 3 Virtual Powertrain Design
Direction is shown with arrows above component icons. The direction can be
changed by sorting it in the tree structure. This can be done by double clicking on a
module in the tree or by clicking on the button. The selected module will then
become the new root module. In the head of the table, symbols indicate in which
direction the modules are calculated: >> X >> forward, << X << backward.
If the component data has been changed, the Kinematic Chain Browser can
calculate once more and update its results when pressing the update button.
The Kinematic Chain Browser also takes into account the ratio changes of the
Component Variation or Matrix Calculation. When opening the Kinematic Chain
Browser, the ratios of the basic model are taken for the calculation. When a certain
variation case is selected in the variations DoE plan, the chain browser calculates
with the ratios of that case.
Activation and deactivation of Task Folders and tasks within Task Folders is
always related to the actual system. Definition can be done in the sidebar of the
vehicle modeler: after the selection of a system, its Task Folders and tasks can be
activated/deactivated in the ‘Tasks’ sheet of the sidebar (Fig. 3.8).
Cycle Run (Fig. 3.9) is a calculation task for the fuel consumption and emissions
calculation. This is done for a fixed driving profile. The determination of this profile
has to be entered as function of the time. This task should be used especially for
The Driver should imitate the behavior of a real driver as well as possible
(Fig. 3.12). For this purpose, the behavior can be defined by different input data.
The used data can be divided into the following groups:
• Data defining the shifting behavior;
• Data defining the starting behavior;
• Data defining the driving behavior.
The data for the starting and driving behavior can be defined by the user (expert
mode) or default values can be used (standard) [1].
Simulation allows modeling of the real behavior if the overall model is har-
monized correspondingly. This will automatically lead to a control loop in which
the driver (modeled as Driver Module in the program) acts as a controller when
simulating a real driving task. For this the driver has to change the acceleration,
brake, or clutch pedal position.
In the elastic connections, the torsional stiffness of the connection (shafts) can be
taken into account.
To solve the differential equation system, a single-step integration method has
been implemented. This method calculates the actual conditions and the conditions
of the next time step. If the accuracy limits cannot be reached with the chosen time
step width (check with a second integration method), the time step width is
decreased [1].
If the temporary results cause a change of the lock conditions, the equation
system is switched to a new one.
3 Virtual Powertrain Design 221
General Configurations are models without any restrictions in the build-up of the
model. The calculation is done without special knowledge about the drivetrain.
In general, all component specific calculations are done in the respective mod-
ules. The values which are determined in the modules are combined together to an
equation system with the smallest possible number of degrees of freedom. The
integration for this reduced system can be done independent of the external modules
in the internal equation solver.
References
1. AVL Cruise (2011) Users Guide, AVL List GmbH, Graz, Austria, Document no. 04.0104.2011,
Edition 06.2011
2. AVL Cruise (2011) Release Note, AVL List GmbH, Graz, Austria, Document no.
04.0102.2011, Edition 06.2011
3. AVL Cruise (2011) Database, AVL List GmbH, Graz, Austria, Document no. 04.0112.2011,
Edition 06.2011
Chapter 4
Classical Powertrain Configuration Model
and Simulation
This chapter describes how to create and run the model of a manual front wheel
drive. To place a component in the working area, click on the component group
with the left mouse button, move the mouse over the required component, and then
click and drag the component into the vehicle model area. Arrange the components
as shown in Fig. 4.1 and connect the components as specified in the diagram.
Place the following components in the working area: Vehicle (1), IC Engine (2),
Torque Converter (3), CVT Continuously Variable Transmission (4), Final Drive
(5), Vehicle Front Right (6), Vehicle Rear Right (7), Vehicle Front Left (8), Vehicle
Rear Left (9), Front Disk Brake (10), Rear Drum Brake (11), Front Disk Brake (12),
Rear Drum Brake (13), Differential (14), Cockpit (15), CVT Control (16), ASC
Anti-Slip Control (17), Catalyst (18), and Monitor (19).
There are three types of energy connections: mechanical, exhaust, and electrical
and are made by connecting polylines between the pins.
All the connections in the Data Bus for this model are described in Table 4.1 [1].
Validation of the connections are marked using bold and are required for the
functioning of the model.
Vehicle (1) is one of the main objects in a model. This component contains
general data of the vehicle, such as nominal dimensions and weights. Only one
vehicle component is required in a model. Road resistances and dynamic wheel
loads are calculated for road and dynamometer runs based on the dimensions and
the load state. The wheel loads are calculated considering motion. The aerody-
namic, rolling, climbing, acceleration, and total resistance are calculated [2].
To define component input data (Table 4.2) double click on the icon or click on
it the right mouse button and select edit to open the following window (Fig. 4.2).
The properties window is accessed through properties button in the input data
window, or from the pop-up menu by selecting “Properties”. Selected parameters
have an influence on the needed component data input for calculation (Fig. 4.3).
With these settings it is also possible to run a calculation with minimal input data,
as the input fields of the unavailable data can be switched off.
IC Engine (2) component contains a model for an internal combustion engine.
The characteristic curves for the full load, fuel consumption, and others can be
freely defined by the user. It is possible to define a gasoline engine as well as a
diesel engine. In this component a temperature model is included to consider the
influence of the temperature on the fuel consumption and emissions while the
engine is cold. The engine will be modeled by a structure of characteristic curves
and maps [2]. To define component input data (Table 4.3) double click on the icon
or click on it with the right mouse button and select edit to open the following
window (Fig. 4.4).
To load IC Engine characteristic click on the Full Load Characteristic tab
(Fig. 4.5). Select (kW) unit from the b column using pull-down menu. Add rows by
clicking on and enter the following data (Table 4.4).
The Motoring Curve (Fig. 4.6) will be defined similar to the Full Load
Characteristic as power versus engine speed (Table 4.5). Engine Maps Basic data
can be loaded from the folder projects\CVT_FWD\CVT_FWD.dbf [1].
In the Fuel Consumption Map (Fig. 4.7, Table 4.6) the actual fuel consumption for
stationary operation is listed depending on the brake mean pressure with the engine
speed as parameter (a—speed (1/min), b—BMEP (bar), c—volume flow (l/h)).
4.1 Classical Vehicle Model 225
This representation is used to make a unique relation between the Brake Mean
Effective Pressure (BMEP) and the fuel consumption for a given engine speed. The
Fuel Consumption Map can also be defined as a Specific Consumption Map (Fig. 4.8,
Table 4.7).
4.1 Classical Vehicle Model 227
50
45
40
35
Power (kW)
30
25
20
15
10
5
0
800 1800 2800 3800 4800 5800
Speed (1/min)
0
-0.5
-1
Power (kW)
-1.5
-2
-2.5
-3
-3.5
-4
800 1800 2800 3800 4800 5800
Speed (1/min)
35 Engine
Speed
30 (1/min)
850
Volume Flow (l/h)
25
1250
20 1700
2300
15 2900
10 3500
4200
5 5100
0 6000
-5 0 5 10 15 20
BMEP (bar)
stationary to mobile operation. The torque converter also multiplies torque. First, an
impeller (pump) converts the mechanical energy from the power unit into fluid
energy (hydraulic fluid); a second transformation, back into mechanical energy,
occurs using the turbine blades [2].
These devices offer economical exploitation when used together with mechan-
ical variable ratio transmissions; operation with nonpositive engagement (slip) is
limited for better efficiency. The torque on the converter input and output depends
on the speed ratio between pump and turbine as well on the absolute angular
velocity of the pump (a—speed ratio (-), b—torque ratio (-), c—pump torque (-) as
shown in Fig. 4.13 and Table 4.11) [1].
To define component input data (Table 4.12) double click on the icon or click on
it with the right mouse button and select edit to open the following window.
CVT (4) can convert every point on the engine operating curve to an operating
curve of its own, and every engine operating curve into an operating range within
the field of potential driving conditions (Fig. 4.14). Its advantage over conventional
230 4 Classical Powertrain Configuration Model and Simulation
4.00 Engine
Specific Volume Flow (l/kWh)
3.50 Speed
(1/min)
3.00 850
2.50 1250
1700
2.00
2300
1.50 2900
1.00 3500
4200
0.50
5100
0.00 6000
0.00 5.00 10.00 15.00
BMEP (bar)
fixed-ratio transmissions lies in the potential for enhancing performance and fuel
economy while reducing exhaust emissions [2].
To define component input data (Table 4.13) double click on the icon or click on
it with the right mouse button and select edit to open the following window [1].
4.1 Classical Vehicle Model 231
1.2 Engine
Speed
(1/min)
1
850
0.8 1250
Mass Flow (kg/h)
1700
0.6 2300
2900
0.4
3500
0.2 4200
5100
0 6000
-4 1 6 11 16
BMEP (bar)
14 Engine
Speed
12 (1/min)
850
10
1250
Mass Flow (kg/h)
8 1700
2300
6
2900
4 3500
4200
2
5100
0 6000
-4 1 6 11 16
BMEP (bar)
0.4 Engine
Speed
0.35 (1/min)
850
0.3
1250
Mass Flow (kg/h)
0.25
1700
0.2 2300
0.15 2900
3500
0.1
4200
0.05
5100
0 6000
-4 1 6 11 16
BMEP (bar)
The Efficiency Fixed curve can be defined as function of the transmission ratio
(Fig. 4.15). The additional torque loss correction fixed curve is defined as function
of the transmission ratio and the drive speed (Fig. 4.16).
Final Drive (5) is a gear step with fixed ratio. It can be used as a transmission
step of the differential (final drive unit). A drive torque will be transferred to a
power take-off torque of the transmission step by considering the transmission, the
mass moments of inertia, and the moment of loss (Fig. 4.17) [2]. To define
4.1 Classical Vehicle Model 235
2.5
2
800
1.5
1
0.5 300
0
0 0.5 1 1.5 2 2.5 -200
-0.5
-1
-1.5 -700
-2
-2.5 -1200
Speed Ratio (-)
component input data double click on the Single Ratio Transmission icon and enter
the following data (Table 4.14).
The Wheels and Tires link the vehicle to the road. The wheel components take
into consideration many influencing variables and their effect on the rolling state
(Fig. 4.18).
236 4 Classical Powertrain Configuration Model and Simulation
The moment of rolling drag can be computed on the basis of the wheel load, the
corrected dynamic rolling radius, and the coefficient of rolling drag.
The longitudinal tire force results from the friction coefficient, the wheel load as
well as from the wheel load factor and the slip factor. It is possible to define
variable friction coefficients along the driving profile when different road conditions
are considered [2].
Information data for the four wheels are the same (Vehicle Front Right (6),
Vehicle Rear Right (7), Vehicle Front Left (8), and Vehicle Rear Left (9)).
Double click on the wheel and enter the following data (Table 4.15). Select
Properties menu and define the appropriate wheel location for each wheel (Vehicle:
Front Right, Rear Right, Front Left, and Rear Left) [1].
4.1 Classical Vehicle Model 237
0.92
0.91
0.9
Efficiency (-)
0.89
0.88
0.87
0.86
0.85
0.85 1.35 1.85 2.35 2.85 3.35 3.85
Transmission Ratio (-)
4.5
4.3
4.1
Torque Loss (Nm)
3.9
3.7
3.5
3.3
3.1
2.9
2.7
2.5
1000 2000 3000 4000 5000 6000
Speed (1/min)
The static rolling radius is the distance between the center of the wheel and the
road surface for the loaded, not moving vehicle. The static rolling radius is not used
for AVL CRUISE calculation, but only for the Kinematic Chain Browser.
4.1 Classical Vehicle Model 239
The dynamic rolling radius is the distance between the center of the wheel and
the road surface for the loaded, moving vehicle. The dynamic rolling radius is
usually a little bit bigger than the static rolling radius.
The dynamic rolling radius can be defined as function of the velocity (Fig. 4.19
and Table 4.16). This characteristic is active when dynamic rolling radius is set
from characteristic (a—velocity (km/h), b—rolling resistance factors (%)).
The Brake component is described by brake data and dimensions. It is possible
to define drum brakes as well as disk brakes. The retarder is used for heavy vehicles
and is described below (Fig. 4.20).
The braking torque is computed considering the braking dimensions and the
input brake pressure. This brake pressure can come from the cockpit component or
brake control. If the vehicle is standing still, the degrees of freedom will be reduced
as this also reduces the calculation time. This reduction is done in a way that the
equation system is switched if a small velocity threshold is reached. In this case
movement is suppressed [2].
240 4 Classical Powertrain Configuration Model and Simulation
At the same time, the instantaneous compulsive force is checked (if it is smaller
than the braking torque). If this condition is no longer true the brake is released.
Input data for the two front brakes are the same and also the data for the two rear
brakes are the same (Front Disk Brake Right (10), Rear Drum Brake Right (11),
Front Disk Brake Left (12) and Rear Drum Brake Left (13)). Double click on the
brake icon and enter the following data (Table 4.17).
Differential (14) unit compensates for discrepancies in the respective rotation
rates of the drive wheels: between inside and outside wheels during cornering and
between different drive axles on 4WD vehicles (Fig. 4.21).
When there is a unilateral variation in road surfaces it results in different coef-
ficients of friction at the respective wheels, this balance effect limits the effective
drive torque to a level defined as twice the traction force available at the wheel (tire)
with the lower friction coefficient [2].
This wheel then responds to the excessive force by spinning. To avoid such
effects a positive lock is available at the component.
In the differential it is possible to define the torque split factor. This is needed if
it is used as a central differential for a four wheel drive.
4.1 Classical Vehicle Model 241
To define component input data double click on the differential icon and enter
the following data (Table 4.18).
The Kinematic Chain Browser can be launched in the pop-up menu of a selected
component (Fig. 4.22). All mechanical connections starting from the selected
component are shown in a tree structure on the left side of the window. When a
second component is selected in the tree structure, all speeds between the selected
components are shown in the table on the right side of the window, each row
representing one gear of the gearbox [2].
Cockpit (15) links the driver and the vehicle (Fig. 4.23). In this component
connections are only made via Data Bus.
The driver gets information such as vehicle velocity and vehicle acceleration, but
also generates information such as the pedal positions for other components. The
pedal positions are transferred to corresponding indicators via the pedal charac-
teristics (Table 4.19) [2].
To define the input data click on the acceleration pedal characteristic (a—
acceleration pedal travel (%), b—load signal (%)) and brake pedal characteristic (a
—specific brake pedal force (%), b—brake pressure (bar)) tabs and enter the fol-
lowing data (Fig. 4.24).
242 4 Classical Powertrain Configuration Model and Simulation
100
90
The CVT Control (16) is an easy way to provide target transmission consid-
ering the velocity and the engine load. In addition, it is possible to control the
clutch. This clutch is simply added to the model (Fig. 4.25).
To define component input data double click on the CVT control icon and enter
the following data (Table 4.20).
244 4 Classical Powertrain Configuration Model and Simulation
To define the input data for desired transmission ratio click on the tabs and enter the
following data (Table 4.21 and Fig. 4.26). The demand transmission of the CVT
gearbox can be defined as a function of a—velocity (m/s), b—load signal (%) and c—
transmission ratio (–). The interpolation of transmission ratio in the maps of CVT
control are done by an alternative method and can be changed in the interpolation and
approximation section in the Expert Mode settings of the Project Settings [3].
To define the input data for desired clutch release click on the tabs and enter the
following data (Fig. 4.27 and Table 4.22). The demand clutch release of the
starting clutch can be defined as a function of a—speed (1/min), b—load signal (%)
and c—clutch release (%) [3].
The ASC Anti-Slip Control (17) checks the force transmission factor (ratio
between force that should be transmitted and maximum transmittable force) of all
connected wheels (driven wheels). If the force that has to be transmitted exceeds the
adhesion limit, the load signal or the clutch releases changes [2].
4.1 Classical Vehicle Model 245
4.00 Velocity
(m/s)
3.50
Transmission Ratio (-) 0
5
3.00
10
2.50 15
20
2.00 25
30
1.50 35
40
1.00 45
50
0.50
0 20 40 60 80 100
Load Signal (%)
The Anti-Slip Control (Fig. 4.28) is only active when a model is calculated
quasi-stationary in combination with special drivetrain configurations and does not
require data entry.
If the Anti-Slip Control is used for calculations in simulation mode this switch
can be activated. In this case the parameters for the PID controller can be set
(proportional (–), differential (s) and integral parameter (1/s)).
246 4 Classical Powertrain Configuration Model and Simulation
Catalyst (18) or exhaust systems consider the effects of the catalytic converter
and soot filler on the raw emissions of the engine.
Starting with the temperature of the catalytic converter, factors for the conver-
sion will be computed for the single emission components using temperature-
dependent maps (Fig. 4.29).
4.1 Classical Vehicle Model 247
The heat loss coefficient is the percentage of the exhaust energy delivered by the
engine which reaches the exhaust system [2].
To define component input data double click on the catalyst icon and enter the
following data (Table 4.23).
Conversion efficiencies of the different emissions (NOx, CO, HC and Soot) can
be defined dependent on the actual temperature of the catalytic converter.
Conversion efficiencies curve can be defined as function of a—temperature (°C)
and efficiency (–) (Table 4.24 and Fig. 4.30—NOx conversion efficiencies,
Table 4.25 and Fig. 4.31—CO conversion efficiencies, Table 4.26 and Fig. 4.32—
HC conversion efficiencies).
The heat transfer Gas/Cat (Table 4.27 and Fig. 4.33) is the proportion of the
exhaust gas which is transferred to the exhaust system and can be defined as a
function of a—temperature (°C) and b—heat-up proportion of exhaust gas
energy (%) [1].
The Monitor (19) can be introduced, if the calculation run must be detected. It is
possible to show some results of the calculation while the calculation is running.
This is done via the online monitor. It is possible to connect up to ten input channels
of the monitor to output channels of different components that are available at the
AVL CRUISE—Data Bus [2].
248 4 Classical Powertrain Configuration Model and Simulation
1.000
0.900
0.800
0.700
Efficiency (-)
0.600
0.500
0.400
0.300
0.200
0.100
0.000
0 100 200 300 400 500
Temperature (°C)
1.000
0.900
0.800
Efficiency (-) 0.700
0.600
0.500
0.400
0.300
0.200
0.100
0.000
0 100 200 300 400 500
Temperature (°C)
1.000
0.900
0.800
0.700
Efficiency (-)
0.600
0.500
0.400
0.300
0.200
0.100
0.000
0 100 200 300 400 500
Temperature (°C)
70
Heat-Up Proportion of Exhaust
65
Gas Energy (%)
60
55
50
45
40
35
-50 150 350 550 750 950
Temperature (°C)
Currently, the output in the online monitor is only available in the form of tables.
More than one monitor component can be used in a AVL CRUISE model. In this
case, more than one instance of the online monitor is activated during calculation.
Click on the description of Data Bus tab and open the following dialog. To select
a Data Bus channel, click to access the available list. Enter the text for the
description and then click on to select the unit (Fig. 4.34).
This chapter describes how to create and run the model of an automatic front wheel
drive vehicle. To place a component in the working area, click on the component
group with the left mouse button, move the mouse over the required component and
then click and drag the component into the vehicle model area. Put the component
with release of the mouse button on the desired place.
Place the following components in the working area: Vehicle (1), IC Engine (2),
Torque Converter (3), Gearbox (4), Final Drive (5), Vehicle Front Right (6),
4.1 Classical Vehicle Model 251
Vehicle Rear Right (7), Vehicle Front Left (8), Vehicle Rear Left (9), Front Disk
Brake (10), Rear Drum Brake (11), Front Disk Brake (12), Rear Drum Brake (13),
Differential (14), Cockpit (15), Gearbox control (16), Gearbox program (17), ASC
Anti-Slip Control (18), Catalyst (19) and Monitor (20).
Arrange the components as shown in Fig. 4.35 and connect the components as
specified in the diagram.
The electric connections for energy transfer are splinted as a function of energy
that is transferred into mechanical, for exhaust gases and electrical connections. The
links between elements are specified in red, green, or blue, depending on their type.
All the connections in the Data Bus for this model are described in Table 4.28
[1]. Validation of the connections marked in bold are required for functioning of the
model.
The Vehicle (1) is one of the main objects in a model. This component contains
general data of the vehicle, such as nominal dimensions and weights. Only one
vehicle component is required in a model.
Road resistances and dynamic wheel loads are calculated for road and dyna-
mometer runs based on the dimensions and the load state. The wheel loads are
252 4 Classical Powertrain Configuration Model and Simulation
angular velocity of the pump (a—speed ratio (–), b—torque ratio (–), c—pump
torque (–), Fig. 4.45 and Table 4.38).
Gearbox (4) transmissions featuring several fixed ratios can maintain a corre-
spondence between the respective performance curves for engine and vehicle
256 4 Classical Powertrain Configuration Model and Simulation
129
109
Power [kW]
89
69
49
29
9
800 1800 2800 3800 4800 5800
Speed [1/min]
(Fig. 4.46). The correspondence with the hyperbola for maximum engine output
will be acceptable or indeed quite good, depending on a multiplicity of factors
including the number of available gears, the spacing of the individual ratios within
the required conversion range, and the engines’ full-load torque curve [2].
4.1 Classical Vehicle Model 257
-2
Power [kW] -4
-6
-8
-10
-12
-14
800 1800 2800 3800 4800 5800
Speed [1/min]
The component Gearbox contains a model for a gear box with different gear
steps. The user can define as many gears as needed. When it is used as an automatic
gearbox, the gear shifting process will be controlled by the control module gearbox
control and program. To define component input data for Gear Ratio Table double
click on the icon and select edit to open the following window [1].
258 4 Classical Powertrain Configuration Model and Simulation
50 Engine
45 Speed
(1/min)
40
Volume Flow (l/h) 700
35
1200
30 1700
25 2550
20 3000
15 3600
10 4150
5 5050
6000
0
-4 1 6 11 16
BMEP (bar)
1.80
Engine
1.60 Speed
(1/min)
1.40
700
Mass Flow (kg/h)
1.20 1200
1.00 1700
0.80 2550
3000
0.60
3600
0.40
4150
0.20 5050
0.00 6000
-5.00 0.00 5.00 10.00 15.00 20.00
BMEP (bar)
The efficiency and torque loss correction can be defined as function of the gear
efficiency (Fig. 4.47).
Final Drive (5) is a gear step with fixed ratio (Fig. 4.48). It can be used as a
transmission step of the differential (final drive unit). A drive torque will be
transferred to a power take-off torque of the transmission step by considering the
transmission, the mass moments of inertia and the moment of loss [2].
To define component input data double click on the final drive icon and enter the
following data (Table 4.39).
The Wheels and Tires link the vehicle to the road. The wheel component
considers many influencing variables and their effect on the rolling state.
4.1 Classical Vehicle Model 259
20
Engine
18 Speed
(1/min)
16
700
14
1200
Mass Flow (kg/h)
12 1700
10 2550
8 3000
3600
6
4150
4
5050
2 6000
0
-5 0 5 10 15 20
BMEP (bar)
0.6 Engine
Speed
(1/min)
0.5
700
0.4 1200
Mass Flow (kg/h)
1700
0.3 2550
3000
0.2
3600
4150
0.1
5050
6000
0
-5 0 5 10 15 20
BMEP (bar)
Information data for the four wheels are the same (Vehicle Rear Right (6),
Vehicle Front Right (7), Vehicle Rear Left (8), and Vehicle Front Left (9)).
Double click on the wheel (Fig. 4.49) and enter the following data (Table 4.40).
Select Properties menu and define the appropriate wheel location for each wheel
(Vehicle: Front Right, Rear Right, Front Left, and Rear Left).
The static rolling radius is the distance between the center of the wheel and the
road surface for the loaded, not moving vehicle.
262 4 Classical Powertrain Configuration Model and Simulation
2.5 1200
2
800
1.5
0 0
0 0.5 1 1.5 2 2.5
-0.5
-400
-1
-1.5
-800
-2
-2.5 -1200
Speed Ratio (-)
0.96
0.955
Efficiency (-)
0.95
0.945
0.94
0 1 2 3 4 5
Gear (-)
The dynamic rolling radius can be defined as function of the velocity (Fig. 4.50
and Table 4.41). This characteristic is active when dynamic rolling radius is set
from characteristic (a—velocity (km/h), b—rolling resistance factors (%)).
The slip at maximum traction (Fig. 4.51) is, together with the specific traction at
infinite slip (Fig. 4.52), the characteristic values for the slip correction factor. The
slip correction factor depends on the friction coefficient of the road [1].
The Brake component is described by brake data and dimensions (Fig. 4.53). It
is possible to define drum brakes as well as disk brakes.
Input data for the two front brakes is the same and also the data for the two rear
brakes is the same (Rear Drum Brake Right (10), Front Disk Brake Right (11), Rear
Drum Brake Left (12) and Front Disk Brake Left (13)). Double click on the brake
icon and enter the following data (Table 4.42).
264 4 Classical Powertrain Configuration Model and Simulation
2.3
2.2
Rolling Resistence Factor (%)
2.1
2
1.9
1.8
1.7
1.6
1.5
1.4
1.3
0 40 80 120 160 200 240 280
Velocity (km/h)
The Kinematic Chain Browser can be launched in the pop-up menu of a selected
component (Fig. 4.55).
The Cockpit (15) links the driver and the vehicle (Fig. 4.56). In this component
connections are only made via the Data Bus. The pedal positions are transferred to
corresponding indicators via the pedal characteristics (Table 4.44) [2].
To define the input data click on the acceleration pedal characteristic (a—
acceleration pedal travel (%), b—load signal (%)) and brake pedal characteristic (a
—specific brake pedal force (%), b—brake pressure (bar)) tabs and enter the fol-
lowing data (Fig. 4.57).
The Gearbox Control (16) is required to define an automatic gearbox
(Fig. 4.58). In the gearbox control, the gear shifting process can be defined auto-
matically without any influence of the driver.
The gearbox control shifts the gears dependent on a speed or velocity. In the
calculation tasks it is decided which of the two shifting strategies is used. The speed
or velocity can be defined via the Data Bus. It is possible to define the engine speed
and the wheel speed as Ref. [2].
266 4 Classical Powertrain Configuration Model and Simulation
15
Slip at Maximum Traction (%)
14
13
12
11
10
9
8
7
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Friction Coefficient of Road (-)
90
Specific Traction at Infinite Slip
88
86
84
82
(%)
80
78
76
74
72
70
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Friction Coefficient of Road (-)
100
90
80
Load Signal (%)
70
60
50
40
30
20
10
0
0 20 40 60 80 100
Acceleration Pedal Travel (%)
curves. The curves are given as a function of the load signal and the engine speed.
The target gear determined by the program is transmitted to the gearbox control that
transmits it to the gearbox.
270 4 Classical Powertrain Configuration Model and Simulation
With the properties switches shifting program 2–5 up to four additional shifting
programs can be defined (a—gear (–), b—load signal (%), c—speed upshifting
(1/min) and d—speed downshifting (1/min)).
With the Data Bus input shifting program selector (which provides a double
value), interpolations between the defined shifting programs can be done
(Fig. 4.60a, b and Tables 4.46, 4.47) [3].
The speeds for upshifting and downshifting can be checked in the following
way: the upshifting speed of the 2nd gear with the special load signal means that at
this speed with this load signal the gearbox control is upshifting from the 2nd into
the 3rd gear. The downshifting speed for the 2nd gear with the special load signal
means that at this speed with this load signal the gearbox control is downshifting
from the 3rd into the 2nd gear.
The ASC Anti-Slip Control (18) checks the force transmission factor of all
connected wheels (Fig. 4.61).
The ASC is only active when a model is calculated quasi-stationary in combi-
nation with special drivetrain configurations and does not require data entry.
If the ASC is used for calculations in simulation mode this switch can be
activated. In this case the parameters for the PID controller can be set proportional
(–), differential (s), and integral parameter (1/s).
Catalyst (19) or exhaust systems consider the effects of the catalytic converter
and soot filler on the raw emissions of the engine (Fig. 4.62).
To define component input data double click on the catalyst icon and enter the
following data (Table 4.48).
Conversion efficiencies of the different emissions (NOx, CO, HC, and Soot) can
be defined dependent on the actual temperature of the catalytic converter.
Conversion efficiencies curve can be defined as function of a—temperature (°C)
and efficiency (–) (Fig. 4.63 and Table 4.49—NOx conversion efficiencies,
Fig. 4.64 and Table 4.50—CO conversion efficiencies, Fig. 4.65 and Table 4.51—
HC conversion efficiencies).
The heat transfer Gas/Cat (Table 4.52 and Fig. 4.66) are the proportion of the
exhaust gas which is transferred to the exhaust system and can be defined as a
function of a—temperature (°C) and b—heat-up proportion of exhaust gas energy
(%) [1].
272 4 Classical Powertrain Configuration Model and Simulation
(a) 100
Gear
90 (-)
1 2 3 4
80 1
70 2
Load Signal (%)
60 3
50 4
40
30
20
10
0
500 1500 2500 3500 4500 5500
Speed Upshifting (1/min)
(b) 100
Gear
90 (-)
1 2 3 4
80 1
Load Signal (%)
70 2
60 3
50 4
40
30
20
10
0
0 1000 2000 3000 4000
Speed Downshifting (1/min)
The Monitor component (20) can be introduced, if the calculation run must be
detected. It is possible to show some results of the calculation while the calculation
is running.
274 4 Classical Powertrain Configuration Model and Simulation
1.000
0.900
0.800
0.700
Efficiency (-)
0.600
0.500
0.400
0.300
0.200
0.100
0.000
0 100 200 300 400 500
Temperature (°C)
Click on the description of Data Bus tab and open the following dialog. To select
a Data Bus channel, click to access at the available list. Enter the text for the
description and then click on to select the unit (Fig. 4.67).
4.1 Classical Vehicle Model 275
1.000
0.900
0.800
Efficiency (-) 0.700
0.600
0.500
0.400
0.300
0.200
0.100
0.000
0 100 200 300 400 500
Temperature (°C)
1.000
0.900
0.800
0.700
Efficiency (-)
0.600
0.500
0.400
0.300
0.200
0.100
0.000
0 100 200 300 400 500
Temperature (°C)
70
Heat-Up Proportion of Exhaust Gas Energy (%)
65
60
55
50
45
40
35
-50 150 350 550 750 950
Temperature (°C)
To run a simulation for each model, there are five different options for calculation
[2]:
• Single Calculation is a regular calculation with the input data supplied in the
components. The output of the result defines which components result will be
written and available in the Result Manager.
• Matrix Calculation is used for multiple calculations and to evaluate the effect of
varying the values of input parameters on vehicle performance, consumption
and emissions. The results can be viewed for each varied calculation.
• Component Variation is used for multiple calculations and the calculation
consists of sub-calculations in which components are replaced by components
of the same type. The results can be viewed in the Result Manager.
• System Variation is used for multiple calculations from a list of all the model
systems. The results can be viewed similar to Component Variation.
• Batch Calculation is used for calculations series of different projects, either
Single Calculation, Matrix Calculation with DoE option and Component
Variation with DoE option.
278 4 Classical Powertrain Configuration Model and Simulation
To run the calculation in the Desk Area click on Calculation Center or click on
AVL blue button in the upper right corner (Fig. 4.68). Click on to “Start
Calculation” button to run the simulation.
After the completion of the simulation click on the Result Manager in the
Navigation Frame (Fig. 4.69) and expand the task of which the results should be
viewed. The selected objects are shown in the right section of the window.
Project information can be viewed in the main window by selecting info in the
tree structure or by double clicking info in the right window. The window shows
general information: project name, author, vehicle model, drivetrain model, date,
and state (regarding errors), then the tree structure of the folder with defined
Calculation Tasks (different status for run or not run is shown using blue/gray
color). The picture of the Vehicle Model is also displayed.
One output file is produced for each component in the Vehicle Model and each
Calculation Task. Data for diagrams are extracted from these files and are shown as
Energy Flow Graphics, presentation diagrams for current tasks and standard dia-
grams for used components.
In the tree structure of the Result Manager all defined Calculation Task folders
are shown separately for each task. The folders with results are shown in color;
tasks which are not simulated are gray. Files for each Calculation Task are located
inside the respective folder. In Result Manager, all components are listed where
results have been written in pages mode or data sources mode (Fig. 4.70).
The result output for components is selected with output in the properties
window of the respective component [4].
4.2 Run Simulation 279
Select as defined in module for output of results in the folder of the corre-
sponding calculation task. There are two types of result outputs available with the
Result Manager: standard diagrams and presentation diagrams.
Standard diagrams (Fig. 4.71) are created separately for each component and
only contain results of the individual component and represent time histories of the
single components.
Standard diagrams have a predefined representation page, which can only be
changed in some parts. Standard diagrams are opened by selecting the respective
calculation task in the tree structure of Result Manager and then selecting the
corresponding component.
The pages mode sidebar is divided into a list of predefined standard diagram
pages for the respective component and TeimView Assistant with a tree structure of
components with all available curves and the current parameter for current diagram
of respective component with the list of all used curves. TeimView has two tabs:
pages and data sources. Pages are used to view results as curves in standard diagram
and data sources in tables as numerical data for the respective component for
analyzing and extracting [2].
The data sources mode sidebar (Fig. 4.72) has a tree structure where the data
source, data file, and table are shown. In data sources view there is a list of all curve
data that can be shown.
280 4 Classical Powertrain Configuration Model and Simulation
The Energy Flow Graphics feature can be selected in the Result Manager
(Fig. 4.73). An animation shows the flow between the drivetrain components over
time of the following parameters: Current (A), Energy Stored (kJ), Energy (kJ),
Power (kW), Speed (1/min), Velocity (km/h), Torque (Nm), Force (N), Voltage
(V).
4.2 Run Simulation 281
The animation can be controlled by using the play, pause, stepwise, stop, repeat,
velocity switches. A time-bar, showing the actual simulation time, can be moved to
a certain time-step with the mouse.
Below the flow graphics a monitor is shown which displays the result curves of
the selected components. With the zoom-option, the horizontal can be scaled so that
long cycles can also be viewed properly.
In addition to the graphical representation, the associated values can also be
viewed in the monitor as well as in the drivetrain flow graphics.
The values of the Energy Flow Graphics can also be displayed in a pie chart.
Clicking the button EFG pie chart opens an area where the segments can be added
or removed with the buttons “+” and “−”.
The association between a segment and a component parameter can be done by
adding the segment in the segment list and then a left mouse-click on a connector of
the component (Fig. 4.74). Several connectors can be chosen for one segment. In
that case, the segment represents the sum of all associated values. The segment list
consists of the absolute segment values as well as of the values in terms of per-
centages [2].
The task of the mode trimap is to construct a 3D surface and describe it by
isolines. The input data (3D points) are projected into the x-y-surface and are
triangulated (in 2D). The isolines are created by an intersection list (Fig. 4.75) [2].
282 4 Classical Powertrain Configuration Model and Simulation
To evaluate the fuel economy characteristics of the IC Engine motored bus a typical
NEDC cycle was simulated (Fig. 4.76).
The total distance of this cycle is 1.0931 m, the total time is 1180 s, idle standing
time is 267 s, average speed including stops is 33.35 km/h, average driving speed
excluding stops is 43.10 km/h, average acceleration is 0.506 m/s2 and maximum
acceleration is 1.042 m/s2.
The statistical percentages of IC Engine operating region are found in conden-
sation map and it is show in Fig. 4.77.
Driving time distribution in consumption map is show in Fig. 4.78, the engine
speed for IC Engine is show in Fig. 4.79, the engine torque for IC Engine is show in
Fig. 4.80, engine power for IC Engine is show in Fig. 4.81 and the mean effective
pressure for IC Engine is show in Fig. 4.82.
4.2 Run Simulation 283
References
1. AVL Cruise (2011) Primer. AVL List GmbH, Graz, Austria, Document no. 04.0103.2011,
Edition 06.2011
2. AVL Cruise (2011) Users Guide. AVL List GmbH, Graz, Austria, Document no. 04.0104.2011,
Edition 06.2011
3. AVL Cruise (2011) Gear Shifting Program (GSP). AVL List GmbH, Graz, Austria, Document
no. 04.0114.2011, Edition 06.2011
4. AVL Cruise (2011) Notes. AVL List GmbH, Graz, Austria, Document no. 04.0102.2011,
Edition 06.2011
Chapter 5
Hybrid Powertrain Configuration Model
and Simulation
This chapter describes how to create and run the model of a hybrid CVT
front-wheel drive. To place a component in the working area, click on the com-
ponent group with the left mouse button, move the mouse over the required
component, and then click and drag the component into the vehicle model area.
Arrange the components as shown in the following figure (Fig. 5.1) and connect the
components as specified in the diagram. All the connections in the Data Bus for this
model are described in Table 5.1 [1].
Place the following components in the working area: Vehicle (1), IC Engine (2),
Clutch (3), CVT Continuously Variable Transmission (4), Final Drive (5), Vehicle
Front Right (6), Vehicle Front Left (7), Vehicle Rear Right (8), Vehicle Rear Left
(9), Rear Drum Brake (10), Front Disk Brake (11), Front Disk Brake (12), Rear
Drum Brake (13), E-Machine (14), Differential (15), Cockpit (16), CVT control
(17), Battery H (18), MATLAB® DLL (19), Delta Torque CE (20), Delta Torque
Brake (21), Total Torque (22), and Monitor (23).
The Vehicle (1) is one of the main objects in a model. This component contains
general data of the vehicle, such as nominal dimensions and weights. Only one
vehicle component is required in a model. Road resistances and dynamic wheel
loads are calculated for road and dynamometer runs based on the dimensions and
the load state. The wheel loads are calculated considering motion. The aerody-
namic, rolling, climbing, acceleration, and total resistance are calculated [2].
To define component input data (Table 5.2) double click on the icon or click on
it with the right mouse button and select edit to open the following window
(Fig. 5.2).
The properties window is accessed through properties button in the input data
window, or from the pop-up menu by selecting “properties.” Selected parameters
have influence on needed component data input for calculation (Fig. 5.3). With
these settings, it is also possible to run a calculation with minimal input data, as the
input fields of the unavailable data can be switched off.
IC Engine (2) component contains a model for an internal combustion engine.
The characteristic curves for the full load, the fuel consumption, and others can be
freely defined by the user [2]. It is possible to define a gasoline engine as well as a
diesel engine. In this component a temperature model is included to consider the
influence of the temperature on the fuel consumption and emissions while the
engine is cold. The engine will be modelled by a structure of characteristic maps.
To define component input data (Table 5.3) double click on the icon or click on
it the right mouse button and select edit to open the following window (Fig. 5.4).
The properties window is accessed through properties button in input data
window, or from pop-up menu (Fig. 5.5).
To load IC Engine characteristic click on the Full Load Characteristic tab
(Fig. 5.6). Add rows by clicking on and enter the following data (Table 5.4) [1].
The Motoring Curve (Fig. 5.7) will be defined similar to the Full Load
Characteristic either as power dependent on the engine speeds (Table 5.5).
At the Fuel Consumption Map (Table 5.6 and Fig. 5.8) the actual fuel con-
sumption for stationary operation is listed depending on the brake mean pressure
with the engine speed (a—speed (1/min), b—BMEP (bar), c—volume flow (l/h)).
In the emission maps the emissions for the different exhaust gases are listed
dependent on the mean pressure with the engine speed as parameter (a—speed
(1/min), b—BMEP (bar), c—mass flow (kg/h), namely NOx emission map in
Table 5.1 Data Bus connections
Component Requires Input Information Component Delivering Output Information
Battery H Ambient Temperature Cockpit Course Ambient
Brake Rear Drum Brake Pressure Delta Torque Brake Brake Pressure
Brake Front Disk Brake Pressure Delta Torque Brake Brake Pressure
Brake Rear Drum Brake Pressure Delta Torque Brake Brake Pressure
Brake Front Disk Brake Pressure Delta Torque Brake Brake Pressure
Clutch Desire Clutch Release Cockpit Desire Clutch Release
Cockpit Gear Indicator CVT Current Gear
Operation Control 0 IC Engine Operation Control
5.1 Hybrid Vehicle Model Creation
Fig. 5.9 and Table 5.7, CO emission map in Fig. 5.10 and Table 5.8, HC emission
map in Fig. 5.11 and Table 5.9) [1].
The Clutch (3) makes stationary idle possible, transition to motion, and inter-
ruption of the power flow (Fig. 5.12). The clutch slips to compensate for the
294 5 Hybrid Powertrain Configuration Model and Simulation
difference in the rotational speeds of engine and drive train when the vehicle is set
in motion. When a change in operation conditions makes it necessary to change
gears, the clutch disengages the engine from the drivetrain for the duration of the
procedure [2].
The clutch can be used together with a CVT gearbox and it is only necessary for
starting. In this case, the clutch is controlled by the control module CVT control.
To define component input data (Table 5.10) double click on the icon and select
edit to open the following window (Fig. 5.13).
The pressure force is the axial force when the frictional surfaces are pressed
together. This pressure force depends on the clutch release. With the pressure force
the transmitted torque can be calculated by considering the dimensions (a—clutch
release (%), b—pressure force (N)).
CVT (4) can convert every point on the engine operating curve to an operating
curve of its own, and every engine operating curve into an operating range within
5.1 Hybrid Vehicle Model Creation 295
135
130
125
Torque (Nm)
120
115
110
105
100
95
90
750 1750 2750 3750
Speed (1/min)
Table 5.4 Full Load Speed (1/min) Torque (Nm) Speed (1/min) Torque (Nm)
Characteristic data
800 90 2500 120
1000 100 3000 125
1500 108 3500 130
2000 115 5000 120
-1
-2
Torque (Nm)
-3
-4
-5
-6
-7
0 500 1000 1500 2000 2500 3000 3500
Speed (1/min)
the field of potential driving conditions (Fig. 5.14). With the model of the CVT
gearbox it is possible to change the transmission between two user-defined
threshold values. The adjusting speed between different transmissions is internal
fixed at a constant value [2].
To define component input data (Table 5.11) double click on the icon and select
edit to open the following window. The efficiency fixed curve can be defined as
function of the transmission ratio (Fig. 5.15).
Final Drive (5) or Single Ratio Transmission is a gear step with fixed ratio. It
can be used as a transmission step of the differential (Fig. 5.16).
298 5 Hybrid Powertrain Configuration Model and Simulation
35
Engine
Speed
30 (1/min)
850
Volume Flow (l/h)
25
1250
20 1700
2300
15
2900
10 3500
4200
5
5100
0 6000
-4 0 4 8 12 16
BMEP (bar)
1.2
Engine
Speed
1 (1/min)
850
Mass Flow (kg/h)
0.8 1250
1700
0.6 2300
2900
0.4 3500
4200
0.2
5100
6000
0
-4 1 6 11 16
BMEP (bar)
To define component input data double click on the final drive icon and enter
the following data (Table 5.12).
The Wheels and Tires link the vehicle to the road. The wheel component
considers many influencing variables and their effect on the rolling state (Fig. 5.17).
Information data for the four wheels is the same (Vehicle Front Right (6),
Vehicle Front Left (7), Vehicle Rear Right (8), and Vehicle Rear Left (9)). Double
click on the wheel and enter the following data (Table 5.13) [1].
5.1 Hybrid Vehicle Model Creation 299
14 Engine
Speed
(1/min)
12
850
10
1250
Mass Flow (kg/h)
8 1700
2300
6 2900
3500
4
4200
2 5100
6000
0
-4 1 6 11 16
0.4
Engine
Speed
0.35 (1/min)
0.3 850
Mass Flow (kg/h)
1250
0.25
1700
0.2 2300
0.15 2900
3500
0.1
4200
0.05 5100
6000
0
-4 1 6 11 16
The Brake component is described by brake data and dimensions (Fig. 5.18). It
is possible to define drum brakes as well as disk brakes.
Input data for the two front brakes is the same and also the data for the two rear
brakes is the same (Rear Drum Brake Right (10), Front Disk Brake Right (11), Rear
302 5 Hybrid Powertrain Configuration Model and Simulation
6000
5000
Pressure Force (N)
4000
3000
2000
1000
0
0 20 40 60 80 100
Clutch Release (%)
0.895
Efficiency (-)
0.893
0.891
0.889
0.887
0.885
0.883
0.881
0.879
0.877
0.875
0.45 0.95 1.45 1.95 2.45
Transmission Ratio (-)
Drum Brake Left (12), and Front Disk Brake Left (13)). Double click on the brake
icon and enter the following data (Table 5.14).
Using E-Machine component (14) together with the battery H, the user can
simulate hybrid systems. The model of the electric machine contains two compo-
nents, the inverter and the electric motor. For this type of model a characteristic map
for the efficiency is used to calculate the loss of power. The thermic model takes the
warm up of the electric machine into account regarding the occurring losses. The
5.1 Hybrid Vehicle Model Creation 305
warm up of the environment with respect to the cooling system due to the electric
machine is not considered in the electric machine component.
The maximal power should be restricted to avoid exceeding the given limit for
the temperature due to the occurring losses. Therefore, the permissible losses are
dependent on the actual temperature of the motor and the maximal moment of
inertia is determined according to these values [2].
Double click on the e-machine icon and enter the following data (Fig. 5.19).
The value of the voltage level should be entered at which the following char-
acteristic maps of the machine have been measured. In addition to the U1 level
(Fig. 5.20 and Table 5.15) the U2 level (Fig. 5.21 and Table 5.16) will be activated
in the property window.
For each level the corresponding voltage value must be defined and the corre-
sponding characteristic maps (voltage for U1 = 80 V, voltage for U1 = 280 V). The
maximum mechanical (power) torque is defined as a function of the machine speed
(a—speed (1/min), b—torque (Nm)) [1].
306 5 Hybrid Powertrain Configuration Model and Simulation
70
60
Torque (Nm)
50
40
30
20
10
0 1000 2000 3000 4000 5000 6000 7000 8000 9000
Speed (1/min)
Table 5.15 Maximum Speed (1/min) Torque (Nm) Speed (1/min) Torque (Nm)
mechanical torque
variation U1 0 65 4500 28
2000 65 5500 23
2500 50 6500 19
3500 36 8500 15
5.1 Hybrid Vehicle Model Creation 307
70
60
Torque (Nm)
50
40
30
20
10
0 1000 2000 3000 4000 5000 6000 7000 8000 9000
Speed (1/min)
Table 5.16 Maximum Speed (1/min) Torque (Nm) Speed (1/min) Torque (Nm)
mechanical torque variation
U2 0 65 4500 28
2000 65 5500 23
2500 50 6500 19
3500 36 8500 15
98
Engine
93 Speed
(1/min)
88
500
83
Efficiency (%)
1000
78 1500
2000
73 2500
68 3500
4500
63 5500
58 6500
8500
53
48
0 10 20 30 40 50 60 70
Torque (Nm)
(a—clutch pedal travel (%), b—clutch release (%)), and brake pedal characteristic
(a—specific brake pedal force (%), b—brake pressure (bar)) tabs and enter the
following data for all (Fig. 5.27).
The CVT Control (17) is an easy way to provide target transmission consid-
ering the velocity and the engine load (Fig. 5.28). In addition it is possible to
control a clutch. To define component input data double click on the CVT control
icon and enter the following data (Table 5.20).
310 5 Hybrid Powertrain Configuration Model and Simulation
100
The desired transmission ratio (Fig. 5.29 and Table 5.21) can be defined as a
function of a—velocity (m/s), b—load signal (%), and c—transmission ratio (−).
The desired clutch release (Fig. 5.30 and Table 5.22) of the starting clutch can be
defined as a function of a—speed (1/min), b—load signal (%), and c—clutch
release (%) [3].
The Battery H (18) component is used for simulating hybrid and electric
vehicles. The basic model consists of a voltage source and an ohmic resistance
(Fig. 5.31) [2].
312 5 Hybrid Powertrain Configuration Model and Simulation
2.9 Velocity
(m/s)
0
2.4 3
Transmission Ratio (-)
6
8
14
1.9 17
19
22
1.4 25
28
31
0.9 33
36
39
42
0.4
44
0 20 40 60 80 100
47
Load Signal (%)
100
Velocity
90 (km/h)
80
Clutch Release (%)
70 0
10
60
20
50
30
40 40
30
20
10
0
0 10 20 30 40 50 60 70 80 90 100
Load Signal (%)
Double click on the battery H icon and enter the following data (Table 5.23).
The resistance is constructed in a way that a big part of the complex processes
within the battery are taken into consideration. 2 RC elements which describe the
concentration overvoltage and the transfer overvoltage can be added optionally.
The thermal behavior of the battery H is described by a thermal substitution
model. The warming caused by the losses inside the battery H and the cooling
caused by convection is taken into consideration.
Single cells can be modelled as well as any combination of cells. Therefore, the
user can construct any desired model.
The idle voltage—charge (b) is the idle voltage of one cell in charge mode.
There is no electrical consumer; it is a function of the State of Charge SOC (a)
(Fig. 5.32 and Table 5.24).
314 5 Hybrid Powertrain Configuration Model and Simulation
The Idle voltage—discharge (b) of one cell is the voltage without an electrical
consumer, there is no flow of current, and it is a function of the State of Charge
SOC (a) (Fig. 5.33 and Table 5.25).
MATLAB® DLL (19) Engineering Measurement System (EMS) interface
allows to integrate controllers and/or mechanical elements, such as special drives,
which have been created using MATLAB®/Simulink™, into the computational
5.1 Hybrid Vehicle Model Creation 315
165
160
Voltage (V)
155
150
145
0 20 40 60 80 100
State of Charge (%)
160
155
150
145
0 20 40 60 80 100
State of Charge (%)
procedure is that AVL CRUISE itself does not need to be compiled again if
MATLAB®/Simulink™ modules are to be integrated (Fig. 5.34) [2].
The input and output variables will be transferred by means of in-ports and out-
ports. These ports create the interface to the MATLAB® module and are available at
the AVL CRUISE Data Bus.
The compiled DLL file will be loaded from folder: project/version/userdata/
_CRUISE_R121_ctrl_3_DLL.dll.
The input values are transferred to the MATLAB® function. From the module,
the results are received and input to the output channels (Table 5.26).
The Delta Torque Combustion Engine (20) component is a function used to
calculate with user-defined functions, in this case the adaption of the acceleration
pedal for the IC Engine with Integrated Starter Generator (ISG).
The Integrated Starter Generator (ISG) merges both the conventional starter and
generator of an automobile in a single electric device. The ISG produces electric
power when the vehicle is running, which is used to supply electric devices and to
charge the battery.
The ISG can help to decelerate the vehicle by generating electric power, which is
called regenerative braking. The electric power generated charges the battery,
reducing fuel consumption. If a clutch disconnects the ISG and the compressor
from the engine during the idle stop, the ISG can drive the air-conditioning com-
pressor via a belt [2].
To define component input data double click on the delta torque combustion
engine icon and enter the following data for function (Fig. 5.35), and for description
of Data Bus (Table 5.27).
The Delta Torque Brake (21) component is a function used to calculate (with
user-defined functions), the adaption brake pedal for mechanical brake with ISG.
Fig. 5.35 Delta torque combustion engine input data specifications window
318 5 Hybrid Powertrain Configuration Model and Simulation
To define component input data double click on the delta torque brake icon and
enter the following data for function (Fig. 5.36), and for description of Data Bus
(Table 5.28).
The Total Torque (22) component is a function used to calculate with
user-defined functions, the summation of torque from IC Engine and ISG. To define
component input data double click on the total torque icon and enter the following
data for function (Fig. 5.37), and for description of Data Bus (Table 5.29).
The Monitor component (23) can be introduced, if the calculation run must be
detected.
It is possible to show some results of the calculation while the calculation is
running. This is done via the online monitor. It is possible to connect up to ten input
channels of the monitor to output channels of different components that are
available at the AVL CRUISE—Data Bus [2].
Click on the description of Data Bus tab and open the following dialog. To select
a Data Bus channel, click to access at the available list. Enter the text for the
description and then click on to select the unit (Table 5.30).
This chapter describes how to create and run the simulation for a hybrid FWD
model with planetary gearbox.
To place a component in the working area, click on the component group with
the left mouse button, move the mouse over the required component, and then click
and drag the component into the vehicle model area. Put the component with
release of the mouse button on the desired place.
Place the following components in the working area: Vehicle (1), IC Engine (2),
Final Drive (3), Vehicle Front Left (4), Vehicle Front Right (5), Vehicle Rear Left
(6), Vehicle Rear Right (7), Front Disk Brake (8), Rear Drum Brake (9), Rear Drum
Brake (10), Front Disk Brake (11), E-Machine (12), Generator (13), Differential
(14), Planetary Gearbox (15), Cockpit (16), ASC (17), Battery H (18), MATLAB®
DLL (19), and Monitor (20).
Arrange the components as shown in the following figure (Fig. 5.38) and con-
nect the components as specified in the diagram [1].
All the connections in the data bus for this model are described in Table 5.31.
The Vehicle (1) is one of the main objects in a model. This component contains
general data of the vehicle, such as nominal dimensions and weights.
To define component input data (Table 5.32) double click on the icon and select
edit to open the following window (Fig. 5.39).
IC Engine (2) component contains a model for an internal combustion engine.
To define component input data (Table 5.33) double click on the icon and select
edit to open the following window (Fig. 5.40).
The Motoring Curve (Fig. 5.42) will be defined similar to the Full Load
Characteristic either as power dependent on the engine speeds (Table 5.35) [1].
At the Fuel Consumption Map (Fig. 5.43 and Table 5.36) the actual fuel con-
sumption for stationary operation is listed depending on the brake mean pressure
with the engine speed as parameter (a—speed (1/min), b—BMEP (bar), c—volume
flow (l/h)) [1].
5.1 Hybrid Vehicle Model Creation 323
Final Drive (3) or Single Ratio Transmission is a gear step with fixed ratio
(Fig. 5.44). It can be used as a transmission step of the differential [2].
To define component input data double click on the final drive icon and enter the
following data (Table 5.37).
The Wheels are: Vehicle Front Left (4), Vehicle Front Right (5), Vehicle Rear
Left (6), and Vehicle Rear Right (7) with the data defined in Table 5.38 and
Fig. 5.45.
324 5 Hybrid Powertrain Configuration Model and Simulation
105
Torque (Nm)
100
95
90
85
80
1000 1500 2000 2500 3000 3500 4000 4500 5000
Speed (1/min)
Table 5.34 Full Load Speed (1/min) Torque (Nm) Speed (1/min) Torque (Nm)
Characteristic data
1000 80 3500 113
1500 90 4000 115
2000 98 4200 115
2500 104 4500 112
3000 109 4700 110
-10
-15
Torque (Nm)
-20
-25
-30
-35
-40
1000 1500 2000 2500 3000 3500 4000 4500
Speed (1/min)
Table 5.35 Motoring Curve Speed (1/min) Torque (Nm) Speed (1/min) Torque (Nm)
data
1000 −8.00 4500 −40.00
The Brake component is described by brake data and dimensions (Fig. 5.46 and
Table 5.39). It is possible to define drum brakes as well as disk brakes (Front Disk
Brake Right (8), Rear Drum Brake Left (9), Rear Drum Brake Right (10), and Front
Disk Brake Left (11)).
5.1 Hybrid Vehicle Model Creation 325
10.00 Engine
Speed
9.00 (1/min)
8.00
1000
Volume Flow (l/h)
7.00 1250
1500
6.00
1750
5.00 2000
2250
4.00
2500
3.00 2750
3000
2.00
3250
1.00 3500
4000
0.00
0.00 1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 9.00
BMEP (bar)
E-Machine (12) characteristic map for the efficiency is used to calculate the loss
of power. Double click on the e-machine icon and enter the following data
(Fig. 5.47 and Table 5.40).
The value of the voltage level should be entered at which the following char-
acteristic maps of the machine have been measured (U1 = 280 V).
The maximum mechanical (power) torque (Table 5.41 and Fig. 5.48) is defined
as a function of the machine speed (a—speed (1/min), b—torque (Nm)).
5.1 Hybrid Vehicle Model Creation 327
240
Torque (Nm)
200
160
120
80
40
0 1000 2000 3000 4000 5000 6000
Speed (1/min)
The efficiency map (Fig. 5.49 and Table 5.42) is defined dependent on a—speed
(1/min), b—power (kW), and c—efficiency (%) [1].
The Generator (13) must supply the vehicle electrical system with a sufficient
current under all operating conditions in order to ensure that the state of charge in
the engine storage device is consistently maintained at an adequate level [2].
Double click on the generator icon and enter the following data (Fig. 5.50 and
Table 5.43).
The maximum mechanical (power) torque (Fig. 5.51 and Table 5.44) is defined
as a function of the machine speed (a—speed (1/min), b—torque (Nm)).
The alternator systems are designed to supply charge voltages of 14 V (with
28 V for heavy utility vehicles) in order to maintain an adequate charge in 12 V (or
24 V) batteries.
The current delivered by the generator results from a map dependent on the
speed and the main voltage. By considering the internal resistance, the
330 5 Hybrid Powertrain Configuration Model and Simulation
100
Engine
95 Speed
(1/min)
90
500
Efficiency (Nm)
85 1000
80 1500
2000
75
2500
70 3000
65 3500
4000
60
4500
55
50
0.00 50.00 100.00 150.00 200.00
Power (kW)
60
40
20
Torque (Nm)
0
-7000 -5250 -3500 -1750 0 1750 3500 5250 7000
-20
-40
-60
Speed (1/min)
velocities of its sun and annular gears. To define component input data double click
on the differential icon and enter the following data (Fig. 5.53 and Table 5.46).
The Kinematic Chain Browser can be launched in the pop-up menu of a selected
component (Fig. 5.54—E-Machine, Fig. 5.55—IC Engine).
5.1 Hybrid Vehicle Model Creation 333
90
Speed
(1/min)
85 -4000
-3500
-3000
80
-2500
Efficiency (%)
-2000
75 -1500
-1000
-500
70 500
1000
1500
65
2000
2500
60 3000
3500
4000
55 5500
-32 -24 -16 -8 0 8 16 24 32
Power (kW)
The Cockpit (16) links the driver and the vehicle. In this component connections
are only made via the Data Bus. The pedal positions are transferred into corre-
sponding indicators via the pedal characteristics (Table 5.47).
To define the input data click on the acceleration pedal characteristic (a—
acceleration pedal travel (%), b—load signal (%)) and brake pedal characteristic (a
—specific brake pedal force (%), b—brake pressure (bar)) tabs and enter the fol-
lowing data for all (Fig. 5.56).
The ASC Anti-Slip Control (17) checks the force transmission factor of all
connected wheels (Fig. 5.57). The Anti-Slip Control is only active when a model is
calculated quasi-stationary in combination with special drive train configurations
and does not require data entry [2].
If the Anti-Slip Control should be used for calculations in simulation mode, this
switch can be activated. In this case the parameters for the PID controller can be set
(proportional—1 (−), differential—0 (s), and integral parameter—1 (1/s)).
The Battery H (18) component is used to simulate hybrid vehicles (Fig. 5.58).
Double click on the battery H icon and enter the following data (Table 5.48).
The idle voltage—charge (b) is a function of the State of Charge SOC
(a) (Fig. 5.59 and Table 5.49), and the idle voltage—discharge (b) is a function of
the SOC (a) (Fig. 5.60 and Table 5.50).
334 5 Hybrid Powertrain Configuration Model and Simulation
100
90
80
Load Signal (%)
70
60
50
40
30
20
10
0
0 20 40 60 80 100
Acceleration Pedal Travel (%)
ten input channels of the monitor to output channels of different components that
are available at the AVL CRUISE—Data Bus [2].
Click on the description of Data Bus tab and open the following dialog. To select
a Data Bus channel, click to access at the available list. Enter the text for the
description and then click on to select the unit (Table 5.52).
338 5 Hybrid Powertrain Configuration Model and Simulation
8.20
Voltage (V)
8.00
7.80
7.60
7.40
7.20
0 20 40 60 80 100
State of Charge (%)
8.40
8.20
7.80
7.60
7.40
7.20
0 20 40 60 80 100
State of Discharge (%)
To run a simulation for each model, there are five different options for calculation
[2]:
• Single Calculation is a regular calculation with the input data supplied in the
components. Output of result defines which components result will be written
and available in the Result Manager.
• Matrix Calculation is used for multiple calculations and to evaluate the effect of
varying the values of input parameters on vehicle performance, consumption,
and emissions. The results can be viewed for each varied calculation.
• Component Variation is used for multiple calculations and the calculation
consists of sub-calculations in which components are replaced by components
of the same type. The results can be viewed in the Result Manager.
• System Variation is used for multiple calculations from a list of all the model’s
systems. The results can be viewed similar to Component Variation.
• Batch Calculation is used for calculations of series of different projects, either
Single Calculation, Matrix Calculation with DoE option and Component
Variation with DoE option.
To run the calculation in the Desk Area click on Calculation Center or click on
AVL blue button in the upper right corner (Fig. 5.63). Click on to “Start
Calculation” button to run the simulation.
After the simulations are complete click on the Result Manager in the Navigation
Frame (Fig. 5.64) and expand the task of which the results should be viewed. The
selected object is shown in the right section of the window.
Project information can be viewed in the main window by selecting info in the
tree structure or by double clicking info in the right window. The window shows
general information: project name, author, vehicle model, drivetrain model, date
and state (regarding errors), then the tree structure of the folder with defined
Calculation Tasks (different status for run or not run is shown using blue/gray
color). The picture of the Vehicle Model is also displayed.
One output file is produced for each component in the Vehicle Model and each
Calculation Task. Data for diagrams are extracted from these files and are shown as
Energy Flow Graphics, presentation diagrams for belonging task and standard
diagrams for used components.
In the tree structure of the Result Manager all defined Calculation Task folders
are shown separately for each task. Those folders with results are shown in color;
tasks which are not run are gray. Files of the each Calculation Task are located
inside the respective folder. In Result Manager, all components are listed where
results have been written in pages mode or data sources mode (Fig. 5.65).
The result output for components is selected with output in the properties
window of the respective component. Select as defined in module for output of
342 5 Hybrid Powertrain Configuration Model and Simulation
results in the folder of the corresponding calculation task. There are two types of
result outputs available with the Result Manager: standard diagrams and
presentation diagrams.
Standard diagrams (Fig. 5.66) are created separately for each component and
only contain results of the individual component and represent time histories of the
single components.
They have a predefined representation page, which can only be changed in some
parts. Standard diagrams are opened by selecting the respective calculation task in
the tree structure of Result Manager and then selecting the corresponding
component.
The pages mode sidebar is divided into a list of predefined standard diagram
pages for the respective component and TeimView Assistant with a tree structure of
components with all available curves and the current parameter for current diagram
of respective component with the list of all used curves. TeimView has two mode
tabs: pages and data sources. Pages are used to view results as curves in standard
diagram and data sources in tables as numerical data for the respective component
for analyzing and extracting [2].
The data sources mode sidebar (Fig. 5.67) has a tree structure where the data
source, data file, and table are shown. In data source view is a list of all curve data
that can be shown.
Energy Flow Graphics feature can be selected in the Result Manager (Fig. 5.68).
An animation shows the flow between the drivetrain components over time of the
344 5 Hybrid Powertrain Configuration Model and Simulation
following parameters: Current (A), Energy Stored (kJ), Energy (kJ), Power (kW),
Speed (1/min), Velocity (km/h), Torque (Nm), Force (N), Voltage (V).
The animation can be controlled by appropriate switches (play, pause, stepwise,
stop, repeat, velocity). A time bar, showing the actual simulation time, can be
moved to a certain time step with the mouse.
Below the Energy Flow Graphics a monitor is shown which displays the result
curves of selected components. With a zoom option, the horizontal direction can be
scaled so that long cycles can also be viewed properly.
In addition to the graphical representation, the associated values can also be
displayed in the monitor as well as in the drivetrain flow graphics.
The values of the Energy Flow Graphics can also be displayed in a pie chart.
Clicking the button EFG pie chart opens an area where the segments can be
added or removed with the buttons “+” and “−”.
The association between a segment and a component parameter can be done by
marking the segment in the segment list and then doing a left mouse click on a
346 5 Hybrid Powertrain Configuration Model and Simulation
so-called connector of the component (Fig. 5.69). Several connectors can be chosen
for one segment. In that case, the segment represents the sum of all associated
values. The segment list consists of the absolute segment values as well as of the
values in terms of percentages [2].
The task of the mode trimap is to construct a 3D surface and describe it by
iso-lines. The input data which is data of 3D points is projected into the x-y-surface
and are triangulated (in 2D). The iso-lines are created by an intersection list
(Fig. 5.70).
To evaluate the fuel economy characteristics of the IC Engine motored bus was
running a typical NEDC cycle (Fig. 5.71). The total distance of this cycle is
10.931 m, the total time is 1180 s, idle standing time is 267 s, average speed
including stops is 33.35 km/h, average driving speed excluding stops is 43.10 km/h,
average acceleration is 0.506 m/s2, and maximum acceleration is 1.042 m/s2.
The statistical percentages of IC Engine operating region are found in conden-
sation map and are show in Fig. 5.72. Driving time distribution in consumption map
is show in Fig. 5.73.
The fuel consumption of all engines can be evaluated overall, and at idle speed,
acceleration, constant ad deceleration as shown in Fig. 5.74, the electrical power,
maximum power out and maximum power input are shown in Fig. 5.75 and for the
5.2 Run Simulation 347
hybrid version also the charge, voltage, current, and State of Charge (SOC) are
important to monitor (Fig. 5.76).
During the test cycle, the electrical power and power loss are also monitored
(Fig. 5.77), and in Fig. 5.78 the total output and total input energy variation can be
seen. The clutch and brake travel are shown in Fig. 5.79, and in Fig. 5.80 the brake
pedal travel is correlated to the power of all brakes, also for the braking, in
348 5 Hybrid Powertrain Configuration Model and Simulation
Fig. 5.81, the braking torque and braking power are shown with the vehicle speed.
The raw emissions are shown in Fig. 5.82, while the speed input and output are
shown in Fig. 5.83.
5.2 Run Simulation 349
Rotational acceleration input and output (Fig. 5.84), clutch torque and maximum
torque (Fig. 5.85); acceleration, velocity, and distance (Fig. 5.86), load signal, brake
pressure, and clutch release (Fig. 5.87); operation control and speed (Fig. 5.88);
desired gear and current gear (Fig. 5.89) are all automatically generated graphs.
350 5 Hybrid Powertrain Configuration Model and Simulation
Fig. 5.75 Electrical power, maximum power output, and maximum power input evaluation
Fig. 5.76 Charge, voltage, current, and State of Charge (SOC) monitoring
The real time, current gear, and velocity are shown in Fig. 5.90; speed input and
output in Fig. 5.91, torque input, output, and loss in Fig. 5.92, torque loss and
efficiency in Fig. 5.93, desired transmission ratio and transmission ratio in Fig. 5.94,
input and output speeds in Fig. 5.95, torque inputs and outputs in Fig. 5.96, torque,
5.2 Run Simulation 351
speed, and mechanical power in Fig. 5.97, load signal, desired torque, and start
switch in Fig. 5.98.
From the electric motor, the voltage and current (Fig. 5.99) can be extracted, but
also: status overvoltage, operating mode, and operation control (Fig. 5.100), speed
352 5 Hybrid Powertrain Configuration Model and Simulation
and rotational acceleration (Fig. 5.101), torque, maximum torque generator, and
maximum torque motor (Fig. 5.102), efficiency, electrical power, and power loss
(Fig. 5.103), total input and output energy (Fig. 5.104), engine torque, power, and
speed (Fig. 5.105), engine torque and mean effective pressure (Fig. 5.106), engine
5.2 Run Simulation 353
activation, mean effective pressure, and fuel consumption (Fig. 5.107), engine
temperature and load signal (Fig. 5.108), HC, NOx, CO, and Soot emissions
(Fig. 5.109), engine rotational acceleration (Fig. 5.110), cold start factor and fuel
consumption (Fig. 5.111), power correction factor and engine torque (Fig. 5.112),
load signal, operator control and start switch (Fig. 5.113).
354 5 Hybrid Powertrain Configuration Model and Simulation
For all inputs and outputs, a specific graph can be made: brake pressure—
cockpit (Fig. 5.114), brake pressure—EMS (Fig. 5.115), brake pressure—output
(Fig. 5.116), load signal—cockpit (Fig. 5.117), load signal—EMS (Fig. 5.118),
load signal—output (Fig. 5.119), torque—input (Fig. 5.120), torque behind fly-
wheel (Fig. 5.121), velocity (Fig. 5.122), battery charge (Fig. 5.123).
5.2 Run Simulation 355
Further inputs and outputs can generate the following graphs: current trans-
mission ratio (Fig. 5.124), actual clutch release (Fig. 5.125), current—input
(Fig. 5.126), load signal—EMS—output (Fig. 5.127), load signal—output
(Fig. 5.128), vehicle distance—input (Fig. 5.129), engine load signal—input
(Fig. 5.130), speed input and output (Fig. 5.131), efficiency and torque loss
356 5 Hybrid Powertrain Configuration Model and Simulation
(Fig. 5.132), air resistance force, resistance force, and traction force (Fig. 5.133),
wheel loads front (Fig. 5.134), wheel loads rear (Fig. 5.135), axle loads
(Fig. 5.136), speed and rotational acceleration (Fig. 5.137), friction coefficient and
friction coefficient exploit (Fig. 5.138), longitudinal force and slip (Fig. 5.139),
5.2 Run Simulation 357
clutch pedal travel and acceleration pedal travel (Fig. 5.140), brake pedal travel and
acceleration pedal travel (Fig. 5.141), distance, velocity, and acceleration of the
vehicle (Fig. 5.142), electrical consumption and fuel consumption (Fig. 5.143), load
signal and vehicle acceleration (Fig. 5.144).
358 5 Hybrid Powertrain Configuration Model and Simulation
Fig. 5.102 Torque, maximum torque generator, and maximum torque motor
364 5 Hybrid Powertrain Configuration Model and Simulation
Fig. 5.107 Engine activation, mean effective pressure, and fuel consumption
Fig. 5.133 Air resistance force, resistance force, and traction force
References
1. AVL Cruise (2011) Primer. AVL List GmbH, Graz, Austria, Document no. 04.0103.2011,
Edition 06.2011
2. AVL Cruise (2011) Users guide. AVL List GmbH, Graz, Austria, Document no. 04.0104.2011,
Edition 06.2011
3. AVL Cruise (2011) Gear shifting program (GSP). AVL List GmbH, Graz, Austria, Document
no. 04.0114.2011, Edition 06.2011
Chapter 6
Electric Powertrain Configuration Model
and Simulation
This chapter describes how to create and run the model of an electric front-wheel
drive. To place a component in the working area, click on the component group
with the left mouse button, move the mouse over the required component, and then
click and drag the component into the Vehicle Model area. The components must
be arranged as shown in the following figure (Fig. 6.1) and connected as specified
in the diagram. All the connections in the Data Bus for this model are described in
Table 6.1 [2].
Place the following components in the working area: Vehicle (1), Final Drive
(2), Vehicle Rear Left (3), Vehicle Front Left (4), Vehicle Rear Right (5), Vehicle
Front Right (6), Rear Disk Brake (7), Front Disk Brake (8), Rear Disk Brake (9),
Front Disk Brake (10), E-Machine (11), Differential (12), Cockpit (13), ASC (14),
Electrical Consumer (15), Battery H (16), E-Drive Function (17), E-Brake &
M-Brake Function (18), Monitor (19), and Constants (20).
The Vehicle (1) is the main object in a model. This component contains general
data of the vehicle, such as nominal dimensions and weights. The aerodynamic,
rolling, climbing, acceleration, and total resistance are calculated [1].
To define component input data (Table 6.2), double click on the icon or click on
it by using the right mouse button and select edit (Fig. 6.2).
In the properties window the parameters which have influence on the needed
component data input for calculation are selected (Fig. 6.3).
Final Drive or Single Ratio Transmission (2) is a gear step with fixed ratio
(Fig. 6.4). It can be used as a transmission step of the differential. To define
component input data, double click on the final drive icon and enter the following
data (Table 6.3).
The Wheel component considers many influencing variables and their effect on
the rolling state (Fig. 6.5). The wheels and tires link the vehicle to the road.
Information data for the four wheels is the same (Vehicle Rear Left (3), Vehicle
Front Left (4), Vehicle Rear Right (5), and Vehicle Front Left (6)). Double click on
the wheel and enter the following data (Table 6.4) [2].
The Brake component is described by brake data and dimensions (Fig. 6.6).
Input data for the two front brakes are the same and also the data for the two rear
brakes are the same (Rear Disk Brake Left (7), Front Disk Brake Left (8), Rear Disk
Brake Right (9), and Front Disk Brake Right (10)). Double click on the brake icon
and enter the following data (Table 6.5).
The E-Machine (11) component can be used either as an electric motor or as a
generator. There are separate characteristic maps for each mode. Together with the
battery H, the user can simulate an electric system. For this model, a characteristic
map for the efficiency is used to calculate the loss of power [1]. Double click on the
e-machine icon and enter the following data (Table 6.6 and Fig. 6.7).
The value of the voltage level has been measured and should be entered in the
following characteristic maps of the machine. For each level, the corresponding
voltage value must be defined and the corresponding characteristic maps (voltage
for U1 = 100 V, voltage for U1 = 320 V). The maximum mechanical (power)
torque (Table 6.7 and Fig. 6.8) is defined as a function of the machines speed (a—
speed (1/min), b—torque (Nm)) [2]. Efficiency map is dependent on speed and
power (or torque) and is shown in Table 6.8 and Fig. 6.9.
The Differential (12) is a particular type of simple planetary gear train that has
the property that the angular velocity of its carrier is the average of the angular
6.1 Electric Vehicle Model Creation 389
velocities of its sun and annular gears (Fig. 6.10). To define component input data,
double click on the differential icon and enter the following data (Table 6.9).
The Kinematic Chain Browser can be launched in the pop-up menu of a selected
component. All mechanical connections starting from the selected component are
shown in a tree-structure in the left side of the window (Fig. 6.11).
390 6 Electric Powertrain Configuration Model and Simulation
The Cockpit (13) links the driver and the vehicle. The pedal positions are
transferred into corresponding indicators via the pedal characteristics (Table 6.10).
In this component, connections are only made via the Data Bus (Fig. 6.12).
To define the input data, click on the acceleration pedal characteristic
(a—acceleration pedal travel (%), b—load signal (%)), clutch pedal characteristic
6.1 Electric Vehicle Model Creation 391
(a—clutch pedal travel (%), b—clutch release (%)) and brake pedal characteristic
(a—specific brake pedal force (%), b—brake pressure (bar)) tabs and enter the
following data for all (Fig. 6.13) [2].
The ASC Anti-Slip Control (14) checks the force transmission factor of all
connected wheels (Fig. 6.14). The Anti-Slip Control is only active when a model is
392 6 Electric Powertrain Configuration Model and Simulation
Table 6.7 Maximum Speed (1/min) Torque (Nm) Speed (1/min) Torque (Nm)
mechanical torque data
0 240 7000 105
500 240 7500 95
1000 240 0 −229
1500 240 750 −229
2000 240 1500 −236
2500 240 2250 −238
3000 240 3000 −239
3500 205 3750 −191
4000 180 4500 −159
4500 160 5250 −136
5000 140 6000 −119
5500 130 6750 −106
6000 120 7500 −95
6500 110 – –
250
200
150
100
Torque (Nm)
50
0
-50
-100
-150
-200
-250
0 1500 3000 4500 6000 7500
Speed (1/min)
95 Engine
Speed
90 (1/min)
500
1000
85 1500
2000
Efficiency (Nm)
2250
80 2500
3000
75 3500
3750
4000
70 4500
5000
5250
65 5500
6000
6500
60 6750
-240 -200 -160 -120 -80 -40 0 40 80 120 160 200 240 7000
Torque (Nm) 7500
To define component input data, double click on the C code icon and enter the
following data for function (Fig. 6.22), and for description of Data Bus (Table 6.17)
[1] (Fig. 6.23).
The Monitor component (19) can be introduced, if the calculation run must be
detected (Fig. 6.24). Click on the description of Data Bus tab and open the
6.1 Electric Vehicle Model Creation 397
following dialog. To select a Data Bus channel, click to access at the available list.
Enter the text for the description and then click on to select the unit (Table 6.18).
Constants (20) component enables the user to define up to 99 constant values
which can be used by other components through the Data Bus. The values can be of
type integer, double or string (Table 6.19) [1].
This chapter describes how to create and run the model of an electric front wheel
drive with range extender. To place a component in the working area, click on the
398 6 Electric Powertrain Configuration Model and Simulation
100
90
80
Load Signal (%)
70
60
50
40
30
20
10
0
0 20 40 60 80 100
Acceleration Pedal Travel (%)
component group with the left mouse button, move the mouse over the required
component and then click and drag the component into the vehicle model area. Put
the component with release of the mouse button on the desired place.
Place the following components in the working area: Vehicle (1), IC Engine (2),
Final Drive (3), Vehicle Rear Left (4), Vehicle Front Left (5), Vehicle Rear Right
(6), Vehicle Front Right (7), Rear Disk Brake (8), Front Disk Brake (9), Rear Disk
Brake (10), Front Disk Brake (11), E-Machine (12), Generator (13), Differential
(14), Cockpit (15), ASC (16), Electrical Consumer (17), Battery H (18), E-Drive
Function (19), E-Brake & M-Brake Function (20), Range Extender Function (21),
Monitor (22), PID Control (23), and Constants (24).
6.1 Electric Vehicle Model Creation 399
Arrange the components as shown in the following figure (Fig. 6.25) and con-
nect the components as specified in the diagram [2]. All the connections in the data
bus for this model are described in Table 6.20.
The Vehicle (1) contains general data of the vehicle (Fig. 6.26), such as nominal
dimensions and weights (Table 6.21).
IC Engine (2) component contains a model for an internal combustion engine.
To define component input data (Table 6.22), double click on the icon or click on it
with the right mouse button and select edit to open the following window
(Fig. 6.27) [1].
400 6 Electric Powertrain Configuration Model and Simulation
600
550
500
Resistance ( )
450
400
350
300
250
200
150
220 270 320 370 420
Voltage (V)
Table 6.12 Resistance Table Voltage (V) Resistance (Ω) Speed (1/min) Torque (Nm)
data
220 160 340 385
240 190 360 430
260 225 380 480
280 260 400 530
300 300 420 590
320 340 – –
350
345
340
Voltage (V)
335
330
325
320
315
310
20 40 60 80 100
State of Charge (%)
The properties window is accessed through properties button in the input data
window, or from the pop up menu (Fig. 6.28).
To load engine characteristic click on the Full Load Characteristic tab (Fig. 6.29
and Table 6.23).
The Motoring Curve (Fig. 6.30) will be defined similar to the Full Load
Characteristic either as power dependent on the engine speeds [2].
402 6 Electric Powertrain Configuration Model and Simulation
350
345
340
Voltage (V)
335
330
325
320
315
310
20 40 60 80 100
State of Charge (%)
In the Fuel Consumption Map (Fig. 6.31 and Table 6.24), the actual fuel con-
sumption for stationary operation is listed depending on the brake mean pressure
with the engine speed as parameter (a—speed (1/min), b—BMEP (bar), c—volume
6.1 Electric Vehicle Model Creation 403
flow (l/h)). In the emission maps, the emissions for the different exhaust gases are
listed dependent on the mean pressure with the engine speed as parameter (NOx
emission map—Fig. 6.32 and Table 6.25, CO emission map—Fig. 6.33 and
Table 6.26, HC emission map—Fig. 6.34 and Table 6.27) [2].
Final Drive (3) or Single Ratio Transmission is a gear step with fixed ratio
(Fig. 6.35). To define component input data, double click on the final drive icon
(Table 6.28).
The Wheel component considers many influencing variables and their effect on
the rolling state (Fig. 6.36). Information data for the four wheels is the same
(Vehicle Rear Left (4), Vehicle Front Left (5), Vehicle Rear Right (6), and Vehicle
Front Left (7)) [2].
Double click on the wheel and enter the following data (Table 6.29).
404 6 Electric Powertrain Configuration Model and Simulation
The Brake component is described by brake data and dimensions (Fig. 6.37).
Input data for the two front brakes is the same and also the data for the two rear
brakes is the same (Rear Disk Brake Left (8), Front Disk Brake Left (9), Rear Disk
Brake Right (10), and Front Disk Brake Right (11)).
Double click on the brake icon and enter the following data (Table 6.30).
Electric Machine (12) together with the battery H the user can simulate an
electric system. For this type of model, a characteristic map for the efficiency is
used to calculate the loss of power [1]. Double click on the electric machine icon
and enter the following data (Fig. 6.38 and Table 6.31).
The value of the voltage level should be entered at which the following char-
acteristic maps of the machine have been measured. For each level, the corre-
sponding voltage value must be defined and the corresponding characteristic maps
(voltage for U1 = 100 V, voltage for U1 = 320 V).
The maximum mechanical (power) torque (Fig. 6.39 and Table 6.32) is defined
as a function of the machines speed (a—speed (1/min), b—torque (Nm)) [2].
Efficiency map is defined dependent on speed and power (or torque) and is shown
in Fig. 6.40 and Table 6.33.
6.1 Electric Vehicle Model Creation 405
The Generator (13) must supply the vehicles electrical system with a sufficient
current under all operating conditions in order to ensure that the State of Charge in
the engine storage device is consistently maintained at an adequate level [1].
Double click on the generator icon and enter the following data (Fig. 6.41).
The maximum mechanical (power) torque (Fig. 6.42 and Table 6.34) is defined
as a function of the machines speed (a—speed (1/min), b—torque (Nm)) [2]. The
efficiency map (Fig. 6.43 and Table 6.35) is defined dependent on a—speed
(1/min), b—power (kW) and efficiency (%).
Differential (14) unit compensates for discrepancies in the respective rotation
rates of the drive wheels (Fig. 6.44). To define component input data, double click
on the differential icon and enter the following data (Table 6.36).
The Kinematic Chain Browser can be launched in the pop-up menu of a selected
component. All mechanical connections starting from the selected component are
shown in a tree-structure in the left side of the window (Fig. 6.45) [1].
6.1 Electric Vehicle Model Creation 407
The Cockpit (15) links the driver and the vehicle. In this component, connec-
tions are only made via the Data Bus (Fig. 6.46). The pedal positions are transferred
into corresponding indicators via the pedal characteristics (Table 6.37).
To define the input data, click on the acceleration pedal characteristic (a—
acceleration pedal travel (%), b—load signal (%)), clutch pedal characteristic (a—
clutch pedal travel (%), b—clutch release (%)), and brake pedal characteristic (a—
specific brake pedal force (%), b—brake pressure (bar)) tabs and enter the following
data for all (Fig. 6.47) [2].
The ASC Anti-Slip Control (16) checks the force transmission factor of all
connected wheels (Fig. 6.48). The Anti-Slip Control is only active when a model is
calculated quasi-stationary in combination with special drivetrain configurations
and does not require data entry.
410 6 Electric Powertrain Configuration Model and Simulation
105
95
85
75
Torque (Nm)
65
55
45
35
25
15
5
1000 2000 3000 4000 5000 6000
Speed (1/min)
Table 6.23 Full Load Speed (1/min) Torque (Nm) Speed (1/min) Torque (Nm)
Characteristic data
800 7 3900 77
1350 20 4250 85
1800 30 4750 94
2100 37 5200 100
2500 45 5500 103
3000 57 5780 103
3300 63 6000 102
3550 69 – –
-2
Power (kW)
-4
-6
-8
-10
-12
1000 2000 3000 4000 5000 6000
Speed (1/min)
45
Engine
40 Speed
(1/min)
35 700
Volume Flow (l/h)
30 1200
25 1700
2550
20
3000
15
3600
10 4150
5 5050
0 6000
-4 0 4 8 12 16
BMEP (bar)
1.400 Engine
Speed
1.200 (1/min)
700
1.000
1200
Mass Flow (kg/h)
0.800 1700
2550
0.600 3000
3600
0.400
4150
0.200 5050
6000
0.000
-5 0 5 10 15 20
BMEP (bar)
18.000
Engine
16.000 Speed
(1/min)
14.000 700
12.000 1200
Mass Flow (kg/h)
1700
10.000
2550
8.000
3000
6.000 3600
4.000 4150
2.000 5050
6000
0.000
-5 0 5 10 15 20
BMEP (bar)
0.500 Engine
Speed
0.450 (1/min)
0.400 700
0.350 1200
Mass Flow (kg/h)
0.300 1700
0.250 2550
0.200 3000
0.150 3600
4150
0.100
5050
0.050
6000
0.000
-5 0 5 10 15 20
BMEP (bar)
The idle voltage—charge (b) is the idle voltage of one cell in charge mode.
There is no electrical consumer; it is a function of the State of Charge SOC
(a) (Fig. 6.52 and Table 6.40). The idle voltage—discharge (b) of one cell is the
voltage without an electrical consumer, there is no flow of current, and it is a
function of the SOC (a) (Fig. 6.53 and Table 6.41) [1].
6.1 Electric Vehicle Model Creation 417
E-Drive (19) component is a function that can be used for calculations (with
user-defined functions) (Fig. 6.54). To define component input data, double click on
the C code icon and enter the following data for function (Fig. 6.55), and for
description of Data Bus (Table 6.42) [1].
418 6 Electric Powertrain Configuration Model and Simulation
250
200
150
100
Torque (Nm)
50
0
0 1000 2000 3000 4000 5000 6000 7000 8000
-50
-100
-150
-200
-250
Speed (1/min)
95 Engine
Speed
90 (1/min)
500
1000
85 1500
Efficiency (%)
2000
2250
80 2500
3000
75 3500
3750
4000
70 4500
5000
5250
65 5500
6000
6500
60 6750
-240 -180 -120 -60 0 60 120 180 240 7000
Torque (Nm) 7500
E-Brake & M-Brake (20) component is a function can be used for calculating
with user-defined functions used to conversion of e-drive torque to brake pressure
(Fig. 6.56). To define component input data, double click on the C code icon and
enter the data for function (Fig. 6.57), and for Data Bus (Table 6.43) [1].
Range Extender (21) component is a function can be used for calculating with
user-defined functions used to range extender control (Fig. 6.58). To define com-
ponent input data, double click on the C code icon and enter the following data for
function (Fig. 6.59), and for Data Bus (Table 6.44) [1].
The Monitor element (22) can be introduced, if the calculation run must be
detected (Fig. 6.60). Click on the description of Data Bus tab and open the fol-
lowing dialog. To select a Data Bus channel, click to access at the available list [1].
Enter the text for the description and then click on to select the unit (Table 6.45).
420 6 Electric Powertrain Configuration Model and Simulation
250
200
150
100
Torque (Nm)
50
0
0 1000 2000 3000 4000 5000 6000 7000 8000
-50
-100
-150
-200
-250
Speed (1/min)
Table 6.34 Maximum Speed (1/min) Torque (Nm) Speed (1/min) Torque (Nm)
mechanical torque data
1 240 1 −230
500 240 750 −230
1000 240 1500 −235
1500 240 2250 −237
2000 240 3000 −239
2500 240 3750 −190
3000 240 4500 −160
3500 205 5250 −135
4000 180 6000 −120
4500 160 6750 −105
5000 140 7500 −95
5500 130 – –
6000 120 – –
6500 110 – –
7000 105 – –
7500 95 – –
PID Control (23) may be used to build up more powerful control subsystems in
AVL CRUISE, combined with other signal processing components (Fig. 6.61).
Double click on the PID control icon and enter the data (Table 6.46).
422 6 Electric Powertrain Configuration Model and Simulation
95 Engine
Speed
90 (1/min)
85 500
1000
Efficiency (%)
80 1500
2000
75 2250
2500
70 3000
3500
65 3750
4000
60 4500
55 5000
5250
-250 -150 -50 50 150 250 5500
Torque (Nm) 6000
100
80
60
b)
40
20
0
0 20 40 60 80 100
a)
To run a simulation for each model, there are five different options for
calculation [1]:
6.2 Run Simulation 425
• Single Calculation is a regular calculation with the input data supplied in the
components. Output of result defines which components result will be written
and available in the Result Manager.
• Matrix Calculation is used for multiple calculations and to evaluate the effect of
varying the values of input parameters on vehicle performance, consumption,
and emissions. The results can be viewed for each varied calculation.
• Component Variation is used for multiple calculations and the calculation
consists of subcalculations in which components are replaced by components of
the same type. The results can be viewed in the Result Manager.
426 6 Electric Powertrain Configuration Model and Simulation
600
550
500
Resistance ( )
450
400
350
300
250
200
150
220 270 320 370 420
Voltage (V)
• System Variation is used for multiple calculations from a list of all the model’s
systems. The results can be viewed similar to Component Variation.
• Batch Calculation is used for calculations series of different projects, either
Single Calculation, Matrix Calculation with DoE option, and Component
Variation with DoE option.
6.2 Run Simulation 427
8.40
8.20
8.00
Voltage (V)
7.80
7.60
7.40
7.20
0 20 40 60 80 100
State of Charge (%)
To run the calculation in the Desk Area, click on Calculation Center or click on
AVL blue button in the upper right corner (Fig. 6.63). Click on to “Start
Calculation” button to running the simulation.
After the simulations were complete, click on the Result Manager in the
Navigation Frame (Fig. 6.64) and expand the task of which the results should be
viewed. The selected object is shown in the right section of the window.
428 6 Electric Powertrain Configuration Model and Simulation
8.40
8.20
8.00
Voltage (V)
7.80
7.60
7.40
7.20
0 20 40 60 80 100
State of Discharge (%)
Project information can be viewed in the main window by selecting info in the
tree structure or by double clicking info in the right window. The window shows
general information: project name, author, vehicle model, drivetrain model, date
and state (regarding errors), then the tree structure of the Folder with defined
Calculation Tasks (different status for run or not run is shown using blue/grey
color). The picture of the Vehicle Model is also displayed.
6.2 Run Simulation 429
One output file is produced for each component in the Vehicle Model and each
Calculation Task. Data for diagrams are extracted from these files and are shown as
Energy Flow Graphics, presentation diagrams for belonging task and standard
diagrams for used components.
430 6 Electric Powertrain Configuration Model and Simulation
In the tree structure of the Result Manager, all defined Calculation Task folders
are shown separately for each task. Those folders with results are shown in color;
tasks which are not run are grey. Files of the each Calculation Task are located
inside the respective folder. In Result Manager, all components are listed where
results have been written in pages mode or data sources mode (Fig. 6.65).
The result output for components is selected with output in the properties
window of the respective component. Select as defined in module for output of
results in the folder of the corresponding calculation task. There are two types of
result outputs available with the Result Manager: standard diagrams and presen-
tation diagrams.
Standard diagrams (Fig. 6.66) are created separately for each component and only
contain results of the individual component and represent time histories of the
single components.
They have a predefined representation page, which can only be changed in some
parts. Standard diagrams are opened by selecting the respective Calculation Task in
the tree structure of Result Manager and then selecting the corresponding
component.
The pages mode sidebar is divided into a list of predefined standard diagram
pages for the respective component and TeimView Assistant with a tree structure of
components with all available curves and the current parameter for current diagram
of respective component with the list of all used curves. TeimView has two modes
tabs: pages and data sources. Pages are used to view results as curves in standard
6.2 Run Simulation 431
/*Conversion routine*/
/*Converting brake Torque into Pressure for all Brakes (Front & Rear)*/
double E-Brake;
/* Check whether Torque of E-Drive is negative that means the load signal to the
E-Motor is negative.
In this case the corresponding brake Pressure will be calculated*/
if (a[0] < 0 )
{
/* The corresponding Brake pressure will be determined by using the following
equation:
M = 2 * P * A * eff * u * r * c
with
M: Brake Torque;
P: Brake Pressure;
A: Brake Piston Surface;
eff: Efficiency;
u: Friction Coefficient;
r: effective Friction Radius;
c: Specific Brake factor.
P=M /(2 *A* eff * u * r* c)
It’s simplify the formula for Rear Brake and Front Brake by using
Brake_Factor_Front and Brake_Factor_Rear in order to use the mean value of
both:
Brake_Factor_Front= AF*effF*uF*rF*cF
Brake_Factor_Rear= AR*effR*uR*rR*cR
It’s also to take into account the following ratios: Ratio of Final Drive (iFD) and
Ratio of Transmission (iTR) between the E-Drive and E-Brake.
[P.S. The ratio iTR is equal 1 because there isn't any Gearbox in the model.]
P= M * iFD * iTR/[2 * (Brake_Factor_Front + Brake_Factor_Rear)/2]
The used channel data are then :
a[0] : eDrive_Torque;
a[1] : iFD;
a[2] : iTR;
a[3] : Brake_Factor_Front;
a[4] : Brake_Factor_Rear.
*/
eBrake = a[0]*2*a[1]*a[2]/(2*(a[3]+a[4]));
}
else
{
/* and if the Torque is positive that means E-Motor drives the vehicle the E-Brake
should be equal Zero */
E-Brake = 0.;
}
/*Calculation of reduced mechanical E-Brake Function to determine a maximum
value By applying negative load the E-Drive applies a resistance which generates
electrical power. In this case the mechanical power required for breaking will be
reduced by the power which the E-Motor has recuperated.
That means the vehicle will at first be electrically braked and The rest required
brake energie will be provided by mechanical brakes */
if ((E-Brake+a[5])>0.0)
{
y[0]= (E-Brake+a[5]);
}
else
{
y[0]= 0.0;
}
diagram and data sources in tables as numerical data for the respective component
for analyzing and extracting.
The data sources mode sidebar (Fig. 6.67) has a tree structure where the data
source, data file, and table are shown. In data source view is a list of all curve data
that can be shown.
Energy Flow Graphics feature can be selected in the Result Manager (Fig. 6.68).
An animation shows the flow between the drivetrain components over time of the
following parameters: Current [A], Energy Stored [kJ], Energy [kJ], Power [kW],
Speed [1/min], Velocity [km/h], Torque [Nm], Force [N], Voltage [V].
The animation can be controlled by appropriate switches (play, pause, step-wise,
stop, repeat, velocity). A time-bar, showing the actual simulation time, can be
moved to a certain time-step with the mouse.
Below the flow graphics a monitor is shown which displays the result curves of
selected components. With a zoom-option, the horizontal direction can be scaled so
that long cycles can also be viewed properly.
6.2 Run Simulation 433
maximum power out, and power input (Fig. 6.74), voltage, current, charge, and
State of Charge (Fig. 6.75), and electrical power and power loss (Fig. 6.76), total
input and output energy and charge (Fig. 6.77), speed, braking torque, and brake
6.2 Run Simulation 441
power (Fig. 6.78), brake pressure, braking torque and brake power (Fig. 6.79),
braking torque and rotational acceleration (Fig. 6.80), acceleration velocity and
distance (Fig. 6.81), acceleration pedal travel, brake pedal force and clutch pedal
442 6 Electric Powertrain Configuration Model and Simulation
Fig. 6.74 Electrical power, maximum power out and power inputs
travel (Fig. 6.82), load signal, brake pressure, and clutch release (Fig. 6.83), real
time, current gear and velocity (Fig. 6.84), rotational accelerations (Fig. 6.85),
torques (Fig. 6.86), power loss and efficiency (Fig. 6.87), torque, speed, and
6.2 Run Simulation 443
For all inputs and outputs, a specific graph can be made: switch and resistance
(Fig. 6.92), e-drive torque input (Fig. 6.93), driver brake pressure (Fig. 6.94), input
vehicle velocity (Fig. 6.95), input vehicle acceleration (Fig. 6.96), input load signal
6.2 Run Simulation 445
(Fig. 6.97), input brake pressure (Fig. 6.98), mod load signal (Fig. 6.99), input load
signal cockpit (Fig. 6.100), speed input and output (Fig. 6.101), rotational accel-
eration input and output (Fig. 6.102), efficiency and torque loss (Fig. 6.103),
446 6 Electric Powertrain Configuration Model and Simulation
Fig. 6.82 Acceleration pedal travel, brake pedal force, and clutch pedal travel
air resistance force, resistance force and traction force (Fig. 6.104), wheel loads
front (Fig. 6.105), wheel loads rear (Fig. 6.106), axle loads (Fig. 6.107), friction
coefficient and friction coefficient exploit (Fig. 6.108), longitudinal force and slip
6.2 Run Simulation 447
(Fig. 6.109), desired, current, minimum, and maximum velocity (Fig. 6.110), clutch
pedal travel and acceleration pedal travel (Fig. 6.111), brake pedal travel and
acceleration pedal travel (Fig. 6.112), clutch pedal travel and gear (Fig. 6.113),
brake pedal travel, and power of all brakes (Fig. 6.114).
448 6 Electric Powertrain Configuration Model and Simulation
Fig. 6.104 Air resistance force, resistance force and traction force
References
1. AVL Cruise (2011) Users guide. AVL List GmbH, Graz, Austria, Document no. 04.0104.2011,
Edition 06.2011
2. AVL Cruise (2011) Primer. AVL List GmbH, Graz, Austria, Document no. 04.0103.2011,
Edition 06.2011
Chapter 7
Creating Virtual Road Infrastructure
AVL Road Importer is an AVL Concerto™ option (Fig. 7.1). The purpose of this
option is to convert any GPS data or x, y data of a real road into a format readable
by AVL Concerto™ [1].
This tool is not a GIS tool in the sense that it does not permit to log GPS data nor
does it have any link to any mapping tool. All the data are provided by the user.
Therefore, the generated road is representative of that data. The only given
information about the road are the points, which the user can limit in case they
cannot be done with the given set of parameters.
This subchapter describes the possibilities of AVL Road Importer as well as its
intended use. To navigate, use the command bar available on the left part of the
screen. It is composed of four sections:
Data—in this section, the user is able to either import GPX or KML files to
AVL Concerto™ or simply choose a file in AVL Concerto™ Data Explorer;
Road—in this section, the user is able to choose what data to work with, but also
prepare it before the main algorithm;
Plot—in this section, the user can find supplementary tools to verify the quality
of the road just generated but also you will be able to specify additional elements
such as speed limits;
Export—in this section, the user will find the necessary tools to export to AVL
InMotion or AVL CRUISE format [1].
The Data submenu is accessible by clicking on the button Data in the command bar.
In the submenu there are the following possibilities (Fig. 7.2):
Open Data Explorer—this opens the AVL Concerto™ Data Explorer and the
user is able to access any AVL Concerto™ file (read). This is also where previously
imported KML and GPX files can be accessed.
Import KML File—the user is able to import KML through this. (by choosing
the KML files wish to import and click ok). KML (Keyhole Markup Language) is
an XML language focused on geographic visualization, including annotation of
maps and images. Geographic visualization includes not only the presentation
of graphical data on the globe, but also control of the user’s navigation in the sense
of where to go and where to look. After importing a KML file, it will appear in the
Channel-Browser. The imported files are saved to the directory/Roads/KML/as a
Concerto Transport File (CTF).
Import GPX File—the user will be able to import GPX through this (by
choosing the GPX files you wish to import and click ok). GPX (GPS eXchange
Format) is a lightweight XML data format for interchange of GPS data (waypoints,
routes, and tracks) between applications and web services on the Internet. After
importing a GPX file, it will appear in the Channel-Browser. The imported files are
saved to the directory/Roads/GPX/as a Concerto Transport File (CTF).
The Road submenu is accessible by clicking on the button Road in the command
bar. In the submenu the user will find the following possibilities (Fig. 7.3):
Channel Selection—by clicking on this button, the user can access the channel
selection page (Fig. 7.4) [1]. Each modification of the parameters needs to be saved
by clicking on the Apply button.
• Project—define the project name here. The project name in itself is not used by
the algorithm but it helps predefine the name of the output files;
• Road—define which channel to use for either latitude or longitude or x and y.
Each zone can be activated by checking the corresponding checkbox. If GPS
data is specified, the algorithm will convert it into UTM coordinates;
• Altitude—defines which channel to use for the altitude. A channel can be
chosen directly into (from data) or the channel containing the barometric
pressure (from barometer). If the barometric pressure is chosen, the algorithm
will convert the pressure into an altitude in meters.
• Algorithm—this section contains all options that only concern the algorithm:
i. Road width—this is the only parameter defining the road. This value is used
to determine if the points will be removed or not depending on if the turn is
doable. The larger the value the harder it is to make the turn, hence more
values will be removed from the data.
ii. Loop circuit—this parameter is still experimental. It enables the creation of
a closed track. Activating this option will not only reconnect the endpoint to
the first point but it will also try to correct the altitude to fully connect the
three dimensions of the road.
Execute—this is how to start the algorithm. Depending on the file size this step
may take some time. Once the algorithm finishes a file appears in the data explorer
(Fig. 7.7) [1]. This file should contain:
• The original data in the measurement key O;
• The parameter used for the algorithm in the measurement key P;
• The results of the algorithm in the measurement key R.
To save this file, the user must right click on it and export it as a CTF, ATF, or
ATF-XML (Fig. 7.8).
The Plot submenu is accessible by clicking on the button Plot in the command bar.
In the submenu there are the following possibilities (Fig. 7.9) [1].
• Elevation profile—clicking on this button, will access the elevation profile page
(Fig. 7.10). The purpose of this page is to verify the quality of the altitude.
• Speed limits—by clicking this button, the speed limits page is accessed
(Fig. 7.11). On this page, the specific speed limits for all parts of the road are
defined.
The user can define the speed limits using the band cursor to select a portion of
the road, and then click on a speed to be applied on this portion. By default the user
may choose between 50, 70, 90, or 110 (km/h). The used unit is the one that the
user specified in parameters submenu [1].
If the default values are not enough, the user can define own values in the
textbox on the right and click on Apply.
The graphs on the bottom left of the page allow the user to see the set speed
limits and it is updated every time a change is made.
The speed limits can be set up by automatic definition. In this section several
methods are available to try and automatically define the speed limit based on the
input vehicle speed coming from the data [1]. These methods are:
i. Method percentage—this algorithm will create speed limit zones based on the
maximum vehicle speed witnessed in the data. Four zones will be identified as:
25, 50, 75, and 100 % of the maximum vehicle speed (Fig. 7.12);
ii. Method-specific country—this algorithm will determine the speed limits to
apply based on the vehicle speed and thresholds defined by the user. To define
the thresholds type them in the popup window (Fig. 7.13);
iii. Method max speed break down every X—this algorithm will determine the
cut-off vehicle speed into zones depending on the distance traveled X where X
472 7 Creating Virtual Road Infrastructure
can be either 5, 10, or 50 m. The algorithm will then apply the maximum of
the vehicle speed in the defined zones as the speed limits for the road
(Fig. 7.14).
For each algorithm above, there will be a detection of zero values for the vehicle
speed. Every time a vehicle stop is detected there will be a stop panel added to the
road automatically that will be the duration of the stop event.
The Export submenu is accessible by clicking on the button Export in the command
bar. In the submenu there are the following options (Fig. 7.15) [1].
• InMotion test run—when clicking on InMotion test run, the user will be
prompted for a directory to save the test run. The extension is ROAD, and this
test run only contains information linked to the road. It cannot be used as a
standalone in InMotion. It can, however, be imported directly in the Road menu
of InMotion. (this test run is not a digitalized road, the user has to import the
road segments it contains);
• AVL CRUISE way file—when clicking on AVL CRUISE way file, the user will
be prompted for a directory to save the way file. The extension is WAY, and this
file contains only the altitude profile versus the distance for now.
474 7 Creating Virtual Road Infrastructure
Reference
1. AVL Road Importer (2013) User manual version 1.7, AVL France Technical Documentation
Chapter 8
Loop Powertrain Simulation
steering wheel, stepping on the gas, brake and clutch pedals, shifting gears in a
manual transmission vehicle, etc., are controlled by the virtual driver. There are two
approaches that can be taken into account to perform the driver actions:
• Simple control—with simple control, the actions are performed by specifying in
advance what will happen at a certain time or distance;
• IPGDriver—the course is controlled by a smart computer driver, which tries to
maintain a course on a test track and can also plan ahead. The driver can be
modified to operate within specified limits. For example, a driver may have a
fast or slow reaction time, he may want to maintain a certain speed, and he could
be trying to go as fast as possible [3].
When the virtual vehicle, virtual road, and virtual driver are taken together, it can
refer to it as the Virtual Vehicle Environment (VVE), since the virtual vehicle
“drives” on the virtual road and is “driven” by the virtual driver (Fig. 8.1).
The CIT (CarMaker Interface Toolbox) consists of a number of tools (applica-
tions and utilities) that run on a host computer. The CIT includes the following:
• The CarMaker GUI—this is the main Graphical User Interface which is used to
control the actions of the VVE, select the virtual vehicle parameter data, define
or select the virtual road, set the virtual driver parameters, define or load
maneuvers;
• The Vehicle Data Set Editor—edit any parameter of the vehicle with a conve-
nient user interface. Each individual submodel class is editable in a separate tab;
• The IPGControl—visualization and analysis tools IPGControl can be used to
view selected output quantities in real-time, load post-simulation data files, and
plot and analyze the results;
To start CarMaker under windows press the Start button and select Programs—IPG
—CarMaker.
The CarMaker main GUI pops up and automatically loads the project folders
where the user was working before last shut down of the program. Nonetheless, it
can easily switch the project folder by selecting File—Project Folder and select the
one project for work [2].
The main window of CarMaker is called the main Graphical User Interface GUI
(Fig. 8.2).
flexible body mode, the main body is split into two bodies. The main body/bodies
include all sprung mass of the vehicle [4].
Engine Mount: CarMaker offers the possibility to simulate an elastically
mounted engine. By activating this feature, an additional mass is added to the
model which needs to be connected to the main body by a joint. The required
settings can be specified in the Engine of the Vehicle Data Set Editor tab (Fig. 8.5).
The Engine Body is parameterized similar to the vehicle body. The position of
the center of gravity in the design frame is required, as well as the moments and
products of inertia [2].
Joint Engine section, the connection point from the engine mass to the vehicle
body is specified. Mounting generalized Joint Force Mode, each of the five force
elements Force X, Force Y, Force Z, Torque X, and Torque Y are made of a spring
and damper element. The amplification factor can be used to easily manipulate the
spring and damper characteristics by a multiplication factor given in this field.
The Characteristic Value description: the linear spring and damper characteristic
defined by a single coefficient in (N/m) for the spring stiffness and in (Ns/m) for the
damping.
1D look-up table description: nonlinear spring and damper characteristic defined
by a look-up table. For the spring, the map requires the force in (N) over
482 8 Loop Powertrain Simulation
compression in (m); the damper element is specified by force in (N) over velocity in
(m/s).
Suspensions for a given axle: the spring characteristics are the same for both
sides of the vehicle. The Spring is modeled as a component that generates a force
when it is compressed or stretched (Fig. 8.6).
For each front and rear axle, the stiffness of the Secondary Spring can be defined
with a simple coefficient if the secondary spring force is a linear function of the
wheel travel variation (Fig. 8.7). If the mode not specified is selected, the secondary
spring is deactivated [2].
The stiffness parameter is either a simple coefficient or a table according to the
selected parameter Mode.
Amplification parameter enables to scale the stiffness of the secondary spring
very quickly for test purposes, instead of modifying the whole table. It can also use
this parameter to convert the values to fit to the units required by CarMaker. The
default value is 1.0. If the value is 0.0, the secondary spring is deactivated [2].
For a given axle, the Damper characteristics are the same for both sides of the
vehicle (Fig. 8.8). The damper is modeled as a component that generates a force
when being compressed or deflected (reaction to the change of velocity) [2].
For each front and rear axle, the damper characteristic is defined either by a
coefficient or by a table of values according to the selected mode.
The characteristics of the damper are split into two domains: push and pull. For
each domain a different damping can be defined. The push domain corresponds to a
positive damper speed, which means that the damper is being compressed.
For a given axle, the Buffer characteristics are the same for both sides of the
vehicle (Fig. 8.9). The buffers are used to limit the wheel travel in one direction or
in both up and down. If the wheel traveled far enough to hit one of the buffers, the
buffer acts like an additional spring [2].
That is why the user has to parameterize two buffers per axle (four in total),
stiffness for each buffer, and the wheel travel from which the buffers are activated in
both positive and negative directions.
The buffer stiffness can be parameterized by a coefficient (mode Characteristic
Value) or by a table of values (mode look-up table 1D).
For each front and rear axle, the buffer stiffness are defined either by a coefficient
or by a table of values according to the selected mode.
484 8 Loop Powertrain Simulation
The characteristics of the buffer are split into two domains: push and pull, where
the push domain corresponds to a positive wheel travel.
In the Stabilizer tab, only the properties of the single Anti Roll Bar
(ARB) element are specified. Its angle or length variation is defined in the kine-
matics of the model [2].
Various Kinematics models are available (Fig. 8.10). The kinematics of the
model has to be specified or each of the front and rear axles separately.
The Compliance model takes into account the transition of the wheel when
forces are applied due to elasticity effects of the suspension (Fig. 8.11).
The user can choose between the following compliance descriptions that are
available in the model option on the compliance tab, for others the user has to select
the mode not specified and define them in an external kinematics file [2].
Wheel Bearing friction activates or deactivates the effect of friction at the wheel
bearings. The wheel bearing friction is characterized by a single coefficient μ, which
defines the friction force on the bearing. The bearing radius R enables to evaluate
the friction torque on the bearing based on the friction force. Using the amplifi-
cation factor amp the friction coefficient can be scaled very quickly instead of
modifying the parameter itself [2].
External Forces are calculated by an external model and can be added to the
spring, damper, or stabilizer forces calculated by CarMaker. The model can be
generated by a C code or using MATLAB®/Simulink™.
8.1 IPG CarMaker 485
The hydraulic power steering unit supports the steering wheel motion on the
base of hydraulic pressure generated by an engine-driven pump.
Using the DVA mechanism the user can implement his own assistance model or
a steering without assistance.
Tires are parametrized within the vehicle data set (Fig. 8.15). The purpose of the
tire models is to calculate the slips and the stress applied to the contract point on the
road (patch point) and then to convert them into stresses applied to the wheel center,
since the rest of the vehicle model directly uses those efforts.
Brake model determines the ratio between the forces at the brake pedal to the
braking torque at each wheel (Fig. 8.16). The brake system consists of hydraulic
models and does not contain any controller model.
table 2D engine speed and gas may be used instead if the user wishes to define the
torque values at mid-loads on his own. This is only necessary if you have a
complete and accurate torque map of an engine and if the tests require a very
accurate engine response [2].
For this model two characteristics have to be parameterized: the maximum
engine torque and the drag torque characteristics (at a given speed the gas pedal is
completely released). CarMaker interpolates the mid-load characteristics automat-
ically, depending on the specified exponent.
A Fuel Consumption Map can only be added using the engine model look-up
table (Fig. 8.19). Fuel tank parameters relating to the fuel and fuel tank can be
carried out on this tab [2].
The Clutch model enables the transfer of all or none of the torque produced by
the engine to the driveline part (Fig. 8.20).
speeds [2]. Under about 100 km/h the influence of the aerodynamic forces is so
small that the user can afford to leave the default values (Fig. 8.25).
In the Sensors tab, there is the possibility to define four types of sensors: Sideslip
Angle Sensors, Body Sensors, Driving Assistance Sensors, and Road Property
Sensors (Fig. 8.26).
Sideslip Angle Sensors are located in the vehicle center of gravity by default.
A Body Sensor is an inertial sensor that can be positioned anywhere on the vehicle
to measure the position, velocity, rotational velocity, acceleration and rotational
acceleration. Up to 10 body sensors can be placed on the vehicle [2].
Driving Assistance Sensors (DASensor) can be parameterized independently,
just by selecting the sensor in the general area. The parameters that are specific to
the selected sensor are displayed under the graphical representation of the vehicle.
The Free Space Sensor Module (FSSensor) is an extended DASensor module,
whose sensor beam is subdivided in equiangular horizontal and vertical segments.
The Traffic Sign Sensor (TSSensor) detects all or preselected traffic sign which
falls within its defined range and horizontal/vertical aperture angles. The sensor
cross-checks if the sign faces the sensor and then sorts all detected signs in
ascending order with distance.
The Line Sensor module is used to detect lines and barriers that have been
defined by the user on the road dialog using the road marking and/or traffic barrier
options as a new bump or marker. It can be compared to an idealized camera.
494 8 Loop Powertrain Simulation
The Road Property Sensor collects road information at a user defined point
ahead of the vehicle. Several attributes like the road bend and its deviation, road
marker attributes (speed limit), or the longitudinal and lateral slope are detected.
This sensor can give information about the relative position: deviation distance,
deviation angle, and current driving lane are calculated. This data can be important
for various applications such as lane keeping assistance, lane departure warning,
autonomous driving; sign detection, energy management, pre-scanning, and opti-
mization of fuel consumption or detection of wheel lift [2].
In the Miscellaneous module the user can create his own vehicle graphics.
Supported file types are PNG images (Fig. 8.27). Using PNG pictures, any photo
can be imported.
A picture of the vehicle in top view and another one in side view is required. To
create a PNG using IPGMovie, the picture needs to be scaled to the length of the
vehicle (1 m = 100 pixel) and saved as a PNG image. Using the graphical tool,
everything around the vehicle needs to be deleted to create a transparent
background.
Then, select all and paste it to another file of the size of 560 × 226 pixels (72
dpi). Move the picture to the left side of the window (in case of a side view
additionally align it vertically to the bottom). The top view should be saved in the
same folder and finally the pictures can be selected in the CarMaker GUI [2].
Loads: The following section will present informations about the additional
loads (Fig. 8.28). With the so-called trim loads there is the possibility to add some
extra loads in order to simulate passengers, luggage, and heavy measurement
equipment. Check if the static position of the vehicle using ModelCheck, the effects
of trim, loads will be included in the height of the vehicle’s center of gravity or the
wheel positions at the start of the configuration.
For some specific simulation purposes, the user may want to simulate a TestRun
with an additional load on the vehicle, but still want to check the static state of the
vehicle without it. The picture represents the vehicle body, defined loads (this is
illustrated by a point) [2].
Maneuver is the concept of CarMaker for the driving scenario (Fig. 8.29) and it
is following maneuver definition, which is split into several maneuver steps
(acceleration, braking, etc.).
3 4
Fig. 8.29 Maneuver menu: 1 list of the maneuver steps, 2 specification of a maneuver step, 3
longitudinal dynamics, 4 lateral dynamics
8.1 IPG CarMaker 497
IPGRoad, the road model of CarMaker, enables to calculate the position of any
point on the road surface, including its coefficient of friction (Fig. 8.30). IPGRoad
describes the centerline of the road (x, y, z), which is symbolized by the white
stripes on the picture above.
Once the road geometry is defined, the track width counting from the centerline
to both sides of the road has to be parameterized, too. The same goes for the friction
coefficient of the road surface. In addition to that, it can redefine the width and
friction coefficient only for a certain section of the road. Once the TestRun has
started, the road definition is automatically linked to IPGMovie.
The IPGRoad dialog box can be accessed from the main GUI by selecting
Parameters—Road. The first tab that appears shows the general settings of the road
(Fig. 8.31). The global settings are [5]:
• Start coordinates x, y, z (m) defines the origin of the road;
• Start direction (deg) defines the angle offset from the first segment regarding the
x-axis;
• Car starts at (m) define at which distance the car should start;
• Driving lane should be used: the left or the right one;
• Country is option for the visualization and availability of the traffic signs;
• Track width (m) is defined on both sides of the centerline;
• Margin width (m) the margins are basically longitudinal stripes on both sides of
the entire road length;
• Friction stripe (m) in order to have a variable friction coefficient along the track
width up to two stripes can be defined.
the road is defined either using digitized road data or the segment-based approach
the user can check it with the Bird’s-Eye view (Fig. 8.33).
An IPGMovie interface tab in the road dialog defines a road in terms of seg-
ments, bumps, markers, vehicle start position, etc. (Fig. 8.34).
CarMaker offers the possibility to use real measurements to define the maneu-
vers (Fig. 8.35).
Several data may be read from a file: vehicle speed, steering wheel angle,
steering torque, pedal position, gear number, park brake position, and the pressure
in the master brake cylinder.
This functionality is referred to as Input From File (IFF), and its corresponding
GUI can be reached via the CarMaker GUI—Parameters—Input From File.
By clicking on the GUI using the right mouse button, a context menu opens. The
user can choose between Reset, Load, Save Configuration, Import, and Edit Aliases.
The configurations can be loaded from the project, shared and central directory
folder SimInput/Settings, whereas the saving of the configurations is possible only
within the project directory folder SimInput/Settings.
The file containing the measurement loads is the Input File tab. Only one file can
be chosen, thus if the measurements are split into several files, the user should
gather them into a single one.
This file must be a text file without commas or other punctuation marks and
should consider the following syntax: each column is referred to as a channel. Tabs
are allowed to separate the columns of each channel: the first line specifies the
channel names. In an optional second line, the unit and/or a factor for each channel
can be specified [8].
502 8 Loop Powertrain Simulation
Any input file requires a Time Channel. Without this channel, the IFF func-
tionality cannot work properly. It can be selected from the available channels using
a drop-down menu on the right side of the Channel in File. It is also possible to
define channel aliases for an automatic selection of a channel. This is possible by
clicking the right mouse button on the GUI and choosing Edit Aliases.
In the Starting Conditions tab the conditions for gear step, velocity, and maxi-
mum gas are used as initial values at the start of the maneuver.
In order to have a steady situation at the start of the maneuver, CarMaker starts a
procedure to reach these values. The parameterized maneuver does not start and the
simulation status in the CarMaker GUI remains on Preparation as long as the
Starting Conditions are not reached.
Selecting and formatting the Quantities to overwrite are the list of the quantities
that may be overwritten. In fact, it cannot be overwritten by the IFF functionality.
This channel is used by functionality, the option Speed Profile for the definition of
the longitudinal dynamics in the maneuver GUI.
The four last quantities are referred to as user defined. This means that the user
can link other quantities to the IFF functionality via the C-interface of CarMaker.
IPGDriver calculates a static desired course at the beginning of the simulation.
During this calculation the driver is able to optimize the path with respect to corner
cutting [5]. IPGDriver is a controller for following a course and a speed controller
on a given track (Fig. 8.36).
IPGDriver enables you to add the control actions of a human driver to your
complete vehicle simulation. These actions include the steering, braking, throttle
position, gear shifting, and clutch operation. To sum up, the driver actions are:
• Choice of driving course within the lane borders (corner cutting);
• Steering;
• Choice of the driving speed according to the course and the vehicle behavior;
• Influence on the speed over the accelerator and brake pedal, as well as over the
clutch pedal and gear selection;
• Regard traffic objects and follow.
CarMaker includes two versions of IPGDriver:
• User parameterized driver: this model is part of every version of CarMaker. The
user parameterized driver can perform a driver adaptation to automatically learn
about the vehicle limits by performing a preset sequence of maneuvers;
• Racing driver: this model is available in the pro version of CarMaker and
enables the optimization of the lap time in function of the speed and to perform a
driver adaptation procedure to automatically find these limits.
The option User parameterized Driver allows parameterization of the driver
model in the tab Standard Parameters.
IPGDriver is able to recognize traffic objects driving ahead. If the check box
Traffic is activated, the driver will follow the traffic objects (Fig. 8.37).
Using the Environment GUI (Fig. 8.38), it is possible to define environmental
conditions like the temperature, the time of day, or the wind velocity for the
8.1 IPG CarMaker 503
simulation [6]. If the model takes these parameters into account, they can influence
the results of the simulation.
The generic model uses a simplified approach to calculate the environmental air
temperature, pressure, and density. Other environmental parameters are kept con-
stant. It is possible to define temperature offsets depending on the time of day
and/or distance. Both offsets can be defined via the look-up tables.
To use IPGControl offline, perform a simulation and select save to the store of
results. In IPGControl load the File—Load File—SimOutput—NameOfPC—Date
—model_name.erg. After selecting the result file to be loaded, it appears in the list
of available data sets [2].
Session log window is opened automatically if errors occur. It can be opened
manually over main menu from: Simulation—session log (Fig. 8.41).
Session log is found in the Simulation menu and that is a record of events for
each simulation in a log file. Indeed, it is useful to keep a history of important or
unusual situations and events during the simulation of a TestRun that does not
disappear when the simulation is finished or the user turns off the computer. This
might be of particular interest for test automation, as the user can trace any potential
errors that might have during a given session.
By pressing the Analyze button the user will proceed to the system analysis. Its
results will be displayed and classified according to the following criteria [2]:
506 8 Loop Powertrain Simulation
• Project: a list of all projects (defined with Settings > Result Folder) that displays
the total simulation time and distance for each project.
• User: all users that have run simulations on the system are listed in this section.
Total simulation time and distance for each project are shown as well.
• TestRun: TestRun files that have been loaded and played in the defined time
frame appear in this section with total simulation time and distance.
8.2 CRUISE—CarMaker Co-Simulation 507
If a AVL CRUISE powertrain is chosen in the vehicle data set CarMaker starts the
AVL CRUISE simulation program during the initialization phase of the TestRun
simulation.
A CarMaker working directory is created over main menu File—Project Folder—
Create Project (Fig. 8.42).
Selecting CarMaker working directory over main menu: File—Project Folder—
Select (Fig. 8.43).
set to dynamic, the brake calculation always stays dynamic and never switches to a
kinematic connection.
The most important aspect is that there are exactly the same numbers of Hub
modules as wheels in the model and they are properly connected or positioned like
in the AVL CRUISE stand-alone model. The hub location change is done in the
properties dialog by choosing a proper location from the pull-down menu.
SAM (System Analysis Mode) folder with system analysis task is added in the
calculation task by default (Fig. 8.47).
The system analysis simulation is used to analyze the parts of the driveline or
components in sub-system. It is switched on/off with a check box on the Calculation
Center application form.
There is no need to define standard vehicle specific tasks such as Cycle Run, Full
Load Acceleration, etc. Components Cockpit, Vehicle and task components Driver,
Course and Profile are also not required.
SAM Task will be used in AVL CRUISE powertrain model for calculation of
powertrain quantities which are required and transferred to CarMaker application
[11].
The model with the interface for co-simulation developed in AVL CRUISE will
be saved in the working folder CarMaker/CRUISE, where it will be loaded and
launched as project work in the application CarMaker [11].
Choose the powertrain folder in the vehicle data set from CarMaker. In the
pull-down menu choose AVL_CRUISE (Fig. 8.48).
To include the AVL CRUISE powertrain in the CarMaker simulation environ-
ment, click on the Browse button to select the desired AVL CRUISE powertrain.
Select Open and AVL CRUISE powertrain will be integrated into the CarMaker
simulation environment as shown in the following figure.
In order to really include AVL CRUISE powertrain in CarMaker simulation
environment, press folder button to select desired AVL CRUISE powertrain. In the
Project input field there must be a path to the AVL CRUISE PRJ file. Select Open
and AVL CRUISE powertrain will be integrated into CarMaker Simulation
Environment (Fig. 8.49).
For the purposes of this manual a AVL CRUISE powertrain model will be
constructed from scratches, using only one standard AVL CRUISE model.
Therefore, select that model using the Start AVL CRUISE command.
In order to make the final step toward including AVL CRUISE powertrain in
CarMaker environment is, enough to select already prepared AVL CRUISE pow-
ertrain model through CarMaker File Explorer and press Open.
AVL CRUISE will open and when selecting Vehicle Model, the user can view
and extract it from the AVL CRUISE GUI screen (Fig. 8.50).
This AVL CRUISE model consists of one system and one sub-system, i.e., all
model components are contained in one sub-system of the model. In order to
prepare AVL CRUISE powertrain model for CarMaker co-simulation the user has
to create two systems, each with three sub-system.
514 8 Loop Powertrain Simulation
First step is to add two sub-system in the current system. Right click on
sub-system 001 in sub-system window part and select add new sub-system options
as shown in Fig. 8.51.
Repeat this procedure to create another sub-system in system 001. The result of
this procedure should be one system with three sub-systems (sub-system 001,
sub-system 002, and sub-system 003) [11].
In order to distinguish all further actions in system and sub-system structure
rename the current system into AVL CRUISE, sub-system 001 into Vehicle,
sub-system 002 into Powertrain, and sub-system 003 into CarMaker. After that,
select all components that are not part of the powertrain and right click somewhere in
the Vehicle Model area and then select sub-system—Vehicle. All selected compo-
nents will be added to the Vehicle sub-system of the AVL CRUISE system [11].
Additionally, the powertrain components of Aut FWD model must be added to
the Powertrain sub-system of AVL CRUISE system. To do so, select all remaining
components in Vehicle Model window again, right click somewhere in Vehicle
Model area and then select sub-system—Powertrain (Fig. 8.52).
The selected components for the Vehicle sub-system are: Vehicle, Wheels,
Brakes, Catalyst, Cockpit, Monitor, and ASC.
Hence that CarMaker sub-system of AVL CRUISE system contains no com-
ponents so far. Therefore, we added already mentioned components CM-Car and
hub to the CarMaker sub-system (Fig. 8.53). CM-Car and hub components can be
found in Modules—Special Modules section. Select CarMaker sub-system in
AVL CRUISE system and use drag&drop technique to add one CM-Car component
and four hub components to that sub-system. Since the hub component in fact
represents brake and virtual wheel, mechanical connection must be made so that
front axle hub components are connected to the differential.
Complete CarMaker Sub-System can be loaded with corresponding modules.
With a right mouse button click on sub-system area and from sub-systems popup
516 8 Loop Powertrain Simulation
menu load CM-Car.sub from directory data\sub_system. After, rename and arrange
CarMaker sub-system, select each component and adjust it with move and stretch to
gain position and shape/look of the current model.
Next step is to make a distinction between AVL CRUISE simulation model and
CarMaker specific AVL CRUISE powertrain model. Therefore, we will first create
8.2 CRUISE—CarMaker Co-Simulation 517
a copy of the built system. Right click on AVL CRUISE system, select duplicate
and rename the created system into something more intuitive CarMaker PT.
Deselect activation and visibility of AVL CRUISE sub-system, while it is not
part of the CarMaker PT system and establish missing mechanical connection from
the differential to each of the front hubs.
There are now two systems (AVL CRUISE and CarMaker PT), each with three
sub-systems (Vehicle, Powertrain, and CarMaker). Next step is to group some of
the sub-systems into each system: AVL CRUISE System includes Vehicle and
Powertrain (Fig. 8.54), CarMaker PT system and includes Powertrain and
CarMaker (Fig. 8.55).
The connections in the data bus for CarMaker sub-system are described in
Table 8.1 and in Fig. 8.56.
The connections in the data bus for the drivetrain sub-system are described in
Table 8.2.
After all necessary connections have been established go to Calculation Center
of AVL CRUISE and activate SAM Task and adjust parameters with the following
values:
• End Time: check Infinite option;
• Time Step: 0.001 (s);
• Simulation 5: FSS Heun.
To run the simulation, first the model must be connected with AVL CRUISE,
(connection with the desired model must be checked), and then the road must be
loaded, with the desired maneuvers.
To check the connection with AVL CRUISE, the model must be opened like
shown in Fig. 8.57.
After pressing the start button, like shown in Fig. 8.39, the simulation starts. In
the IPGControl data window, the user can monitor any parameters by selecting
them from the quantities tab.
For the electric model the four main monitored properties are: the Current [A],
the Voltage [V], Energy [kJ], and Energy Stored [kJ]. They are selected from
Quantities—AVL.CRUISE.efg like shown in Fig. 8.58 and monitored in real time,
as seen in Fig. 8.59.
If the user is interested in the powertrain, the parameters that can be monitored
are shown in Fig. 8.60 or if the interest in the Driver, the parameters are shown in
Fig. 8.61.
In order to see the whole variation of the parameters, total fit must be selected
from the display tab.
After the simulation stops, and after pressing total fit, the recorded variations can
be exported either as a. csv file and post-processed further like in Figs. 8.62, 8.63
and 8.64, or they can be saved as they were recorded and printed (Fig. 8.65).
By exporting the Current [A], the Voltage [V], Energy [kJ] and Energy Stored
[kJ], the user can analyze the state of the battery, the recovered energy either versus
elapsed time, or by the track. Basically, the user can conclude which is the optimal
track for recovered energy or underline the optimal range.
Fig. 8.59 Monitoring the electric model in the IPGControl data window
522 8 Loop Powertrain Simulation
322
Voltage [V]
321.9
321.8
321.7
321.6
321.5
0 50 100 150 200 250 300 350
Time [s]
8.00E+03
Energy [kJ]
6.00E+03
4.00E+03
2.00E+03
0.00E+00
-2.00E+03
0 50 100 150 200 250 300 350
Time [s]
300
250
200
150
100
50
0
50-
0 50 100 150 200 250 300 350
Time [s]
524 8 Loop Powertrain Simulation
3.00E+00
2.00E+00
1.50E+00
1.00E+00
5.00E-01
0.00E+00
-5.00E-01
0 50 100 150 200 250 300 350
Time [s]
References
1. IPG CarMaker (2014) Installation guide version 4.5. IPG Automotive, Karlsruhe, Germany
2. IPG CarMaker (2014) User’s guide version 4.5.2. IPG Automotive, Karlsruhe, Germany
3. IPG CarMaker (2014) IPG driver release notes 6.4. IPG Automotive, Karlsruhe, Germany
4. IPG CarMaker (2014) Vehicle data request form data required for parametrization and
validation of a vehicle data set in CarMaker 4.5. IPG Automotive, Karlsruhe, Germany
5. IPG CarMaker (2014) IPG driver user manual 6.4. IPG Automotive, Karlsruhe, Germany
6. IPG CarMaker (2014) Reference manual version 4.5. IPG Automotive, Karlsruhe, Germany
7. IPG CarMaker (2014) Programmer’s guide version 4.5.2. IPG Automotive, Karlsruhe,
Germany
8. IPG CarMaker (2014) Performance report version 4.5. IPG Automotive, Karlsruhe, Germany
9. IPG CarMaker (2014) Quick start guide version 4.5.2. IPG Automotive, Karlsruhe, Germany
10. AVL Cruise (2010) Coupling cruise—drive. AVL list GmbH, Graz, Austria, Document no.
04.0116.2010, Edition 07.2010
11. AVL Cruise (2010) Cruise–CarMaker/TruckMaker co-simulation. AVL List GmbH, Graz,
Austria, Document no. 04.0115.2010, Edition 07.2010
12. AVL Cruise (2011) Interfaces. AVL List GmbH, Graz, Austria, Document no. 04.0110.2011,
Edition 06.2011