Matlab Model
Matlab Model
Eduardo Perondi
Universidade Federal do Rio Grande do Sul
Filiada à
AV. PAULISTA, 2073 - HORSA II - CJ. 1003 - CEP 01311-940 - SÃO PAULO – SP
Downloaded from SAE International by Univ of California Berkeley, Wednesday, August 01, 2018
2005-01-3969
that could be used to demonstrate it. An analysis in this The approach adopted in this work is alternative
sense should involve requirements of vehicle potency as a because we propose the use of a flow-mechanic system.
function of the transport cycle to determine the energy ratio This study is based on the analogy with the electric
that can be absorbed by the brakes. regenerative systems, aiming costs decreasing and faster
Wicks and Donnelly (2002) determine the work cycle and easier application.
for an urban bus and perform a regenerative brakes
application potentiality analysis. They show that the portion 3. REGENERATIVE BRAKE DESIGN
of 91% of total energy is wasted by the brakes, being the
remaining energy dissipated by the air and tire-soil This work reports the design of a pneumatic
frictions. It is also verified that 59% of the work produced regenerative brake system to be applied in a commercial
by the motor is consumed by the brakes, being, in an ideal vehicle. The proposed pneumatic system must be installed
situation, the energy saving potential. These results indicate in parallel to the traditional brake, seeking global energy
a clear justification for the research in this area. efficiency improvement. Therefore, it is necessary to
Today, the most advanced studies are in the application accomplish the determination of a system with potential
of electric regenerative systems. This technology is mainly application viability and achieve the energetic efficiency
based in the use of the speed difference between a rotaries study with the regenerative system integration. The
axis and a fixed magnetic field that induces the efficiency study is done through simulations using a vehicle
manifestation of electric current. Previously, in the past dynamics model integrated with a regenerative system
systems, the current induced during the activation of the model. A private transport vehicle (urban zone bus) which
brake in electric motors was dissipated by Joule effect in works with compressed air to execute some actions as
benches of resistors, wasting the energy. Due to the automatic door driving and in the passive brake potency
technological progress of electronics and control system (air brake) was used.
techniques, it is possible at present to absorb and store this The solution of the equations set seeking the system
electrical energy. dynamic behavior determination is accomplished with the
In Flinders, Mathew and Oghanna (1995) it is presented use of the numeric integration software Matlab® (2002).
the implantation in trains of regenerative systems based on Through the Simulink module, the blocks diagrams are
modern tiristors application witch allows high potency built and the set of associated differential equations is
rectification. The regenerative system, used in north of solved. The simulation results allow the observation the
New Zealand railroad, provided 13% of reduction in energy influence of the regenerative system in the dynamic
waste. In the same paper, it is also shown the measurement behavior and the energetic efficiency.
of energy losses during the application of dynamic brakes The parameters of the dissipative brake were defined so
in the source line of Central Queensland. In this railroad, that the system behavior performs accordingly to the
the coal is transported by a distance of 400km. The total Brazilian National Standards (ABNT, 2004). The
energy wasted during the activation of the brake is of regenerative brake system was incorporated in a way that
8,7MWh. Taking into account the whole process, the modified system maintains the same operational and
approximately 100GWh are wasted. functional characteristics of the previous system. The
Besides its application in trains, regenerative systems simulations were accomplished using a private transport
are used in hybrid vehicles (HEV, hybrid electric vehicles), cycle model. The results were analyzed considering the
in which there is an auxiliary electric motor coupled to the energy savings in the braking process and in the global
combustion motor. This vehicle conception presents context, aiming at the evaluation of the economical
advantages, as smaller dimensions and weights and the potential of the regenerative system application.
appropriate rotation control of the combustion motor, Automotive vehicles can be used in several
causing the economy of fuel and emissions reduction as environmental conditions, so their requirements can vary
said by He and Hodgson (2002). In Cikanek and Bailey significantly. The energy consumed in a vehicle is mainly
(2002), a Regenerative Braking System (RBS) is proposed due to the resistance of the air, to the friction of the tires
to be assembled to a Parallel Hybrid Electric Vehicle with the soil and the internal dumping of the movable parts,
(PHEV). Simulation results from a dynamic model of the besides the amount consumed by the brakes. Therefore, it is
PHEV exhibiting the regenerative braking performance not appropriate the use of a regenerative system for
improvement are presented. vehicles that travel broad distances because the braking
Panagiotidis, Delagrammatikas and Assanis (2000) action is not frequent in this operational situation. Smaller
report fuel savings of 15% in some HEV vehicles vehicles have space restriction, hindering the application of
configurations. These are values much more conservatives the proposed pneumatic regenerative system. So, the most
than the ideal 59% presented previously, but are still very appropriate vehicles for the system implantation are buses
significant concerning energy efficiency performance. The and urban trucks that have intense accelerations cycles and,
advanced technology involving hybrid vehicles limits its consequently, considerable efficiency losses caused by the
application to few commercial models. frequent activation of brakes. This is the main situation
usually addressed in the electric regenerative brakes
Downloaded from SAE International by Univ of California Berkeley, Wednesday, August 01, 2018
d 2T dT (4)
Jt 2
T f Cang Ftx .rt 0
dt dt
Equations (1) to (4) allows the composition of the calculated using the tangent speed of the vehicle and
block diagram elaborated for Simulink/Matlab equivalent tangent speeds of the front and back wheels and
implementation. Figure 3 presents the system main block the horizontal force of the vehicle. The vertical reactions
diagram. are calculated in each wheel using the model of suspension
of Maclay (1997), where the stiffness and dumping
coefficients of the suspension and the gravity center
position are considered. The reactions are multiplied by the
friction coefficients and the result is used to calculate the
front and rear friction forces. The zero value is imposed to
the sliding movement when the speed arrives to null value
to avoid singularity in Eq. (6).
The transmission system is elaborated based on the
generalization of the Eq. (4). In this system, the torque due
to the action of the brakes, Tf, is calculated in a subsystem
block (the “brake subsystem block”). In the transmission
block the back and front horizontal forces are used to
calculate the equivalent angular speeds in each wheel
which depend on the action of the dissipative brakes, on the
Figure 3 – Main block diagram of the model. loss for the internal friction in the two wheels and on the
action of the regenerative brake in the back wheels.
The friction coefficient µ is a function of many The brake subsystem block diagram is a subsystem of
parameters like environment conditions, vertical loads, the transmission system. In this subsystem, the hydraulic
surface types, translation speeds and internal pressures of pressure of the brake action is multiplied by the actuator
the tires. For analysis simplification purposes, in this study, cylinder piston area, resulting in the equivalent force. This
only the initial speed is tested with different values. The force, through the equation for drum brake from Gillespie
other parameters are considered constants. (1992), results in the equivalent torque. The drum brake
action is then added to the regenerative torque which
The braking action is simulated with dry asphalt. The calculus is based on the analysis of the free body diagram
determination of the coefficient is done through the presented in Fig. 5.
correlation graph between relative sliding and friction
coefficient, supplied by Gillespie (1992) and reproduced in
Fig. 4.
1
0,9
0,8
Coeficiente de Atrito
friction coeficient
0,7
0,6
0,5
Figure 5 – Free body diagram of the air compression
0,4
0,3
system.
0,2
0,1 The relationship between the angular speed wa of the
0 axis and the translational speed of the piston vp given by
0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100%
Mabie (1980) is presented in Eq. (7):
sliding ratio
Deslizamento (%) (%)
Relativo
R p Lp sen2T
vp R p wa ( senT ) (7)
Figure 4 – Relative sliding versus Friction coefficient 2
The sliding is defined by Gillespie (1992) as where Rp is the radius of the "crank", Lp is the length of the
"connecting rod" and ș is the angle formed between the
v wr radius Rp and the normal vector of the traverse area of the
desl (6)
v piston. The relationship between the mass flow m dm / dt
and the translational speed is presented in Eq. (8):
where v is the tangent velocity, w is the angular velocity of
the wheel and r is the radius of the wheel. In the sliding m UA p v p (8)
subsystem the horizontal front and back forces are
Downloaded from SAE International by Univ of California Berkeley, Wednesday, August 01, 2018
Where ȡ is the atmospheric air density and Ap is the Table 1 – List of Parameters
traverse area of the piston.
The variation of the pressure in the tank (Eq. (9))
caused by the injection of this mass flow, being considered M = 12000 [kg] n = 0,150 [m] Kf = 480000 [Nm]
an adiabatic isentropic process with small temperature
variation, is given by Perondi and Guenther (2003).
Car = 40 [Ns/m] m = 0,150 [m] Kt = 720000 [Nm]
dp RrT (9)
m
dt Vt J1 = 25 [kgm²] rr = 0,304 [m] Cvf =40000 [Ns/m]
This torque is then added along the time to the torque rt = 0,49 [m] d1 = 2,95 [m] R = 1,4
given by the traditional dissipative brakes. The pressure
generated by the regenerative system is accumulated in the
vessel. Af = 0,005 [m²] d2 = 2,5 [m] Rp = 0,07 [m]
The elastic potential energy accumulated in the tank
through the pressure increasing the saved is expressed by
Eq. (11):
mif = 0,41 d3 = 1,5 [m] Lp = 0,20 [m]
E p.Vt (11)
e = 0,300 [m] Iyy p0 = 101325 [Pa]
=88500[kgm²]
5. RESULTS AND ANALYSIS
Figure 8 – Acceleration of the vehicle. Figure 12 shows the vessel internal pressure for the
ramp braking curve while Fig. 13 presents the vessel
. Note that there is the presence of oscillations in the internal pressure in another test case in which a pulsed
blue curve caused by the alternative loading of the braking action is simulated.
compressor action. The translational speed is presented in
Fig. 9.
Figure 9 – Translational speed. Figure 12 shows that the pressure in the vessel at the
end of compression process is 0,40MPa (equivalent to 4,0
The back wheel angular speed of the regenerative bar). This pressure is accumulated in the 0,1m³ vessel. The
system decreases faster than the system with traditional stored energy is calculated through Eq. (11), and it is equal
brakes, as can be observed in Fig. 11. This effect is caused to 40kJ. The total initial kinetic energy present in the
Downloaded from SAE International by Univ of California Berkeley, Wednesday, August 01, 2018
vehicle is of 1653kJ, that represents a conservation of in vehicle accessories, as the windshield wiper, the
2,45% of energy. compressor of the air conditioner, in the doors drives, etc.
There is also the possibility combustion motor efficiency
improvement through compressed air injection, increasing
the thermal efficiency by a major compression ratio.
Further studies can accomplish the design of a
compressor actuation smoothing system to be applied to
avoid oscillations in the acceleration curve, as shown in
Fig. 8. Also, anti-lock braking systems (ABS) can be
adopted, in parallel to the regenerative system, to avoid the
wheels blocking. Alternatives, as the use of small and
multiple dispersed pressure vases along the structure of the
vehicle should be considered for the warranty of the
equality of the operation of the brake. Other necessary
Figure 13 – Tank pressure (“pulse” action case). study can accomplish the global transport system with
emphasis in the use of the pneumatic energy saving. This
In the simulation presented in Fig. 13, the brake is study must involve the quantification of the amount of air
actuated in intervals of 2s with period of 1s. As it can be consumed by the accessories in most of the urban transport
verified in Fig. 3, the final pressure is approximately 0,98 vehicles to investigate the economical viability of the
MPa, corresponding to 98kJ in accumulated energy. This proposed system. The design and construction of a
final pressure represents the conservation of 5,92% of the prototype for experimental tests is also necessary to verify
available kinetic initial energy. the system viability.
The appearance of the ISSN code at the bottom of this page indicates SAE’s
consent that copies of the paper may be made for personal or internal use of
specific clients. This consent is given on the condition however, that the copier pay a
$ 7.00 per article copy fee through the Copyright Clearance Center, Inc. Operations
Center, 222 Rosewood Drive, Danvers, MA 01923 for copying beyond that permitted
by Sections 107 or 108 of U.S. Copyright Law. This consent does not extend to
other kinds of copying such as copying for general distribution, for advertising or
promotional purposes, for creating new collective works, or for resale.
SAE routinely stocks printed papers for a period of three years following date of
publication. Direct your orders to SAE Customer Sales and Satisfaction Department.
Quantity reprint rates can be obtained from the Customer Sales and Satisfaction
Department.
ISSN 0148-7191
© Copyright 2005 Society of Automotive Engineers, Inc.
Positions and opinions advanced in this paper are those of the author(s) and not
necessarily those of SAE. The author is solely responsible for the content of the
paper. A process is available by which discussions will be printed with the paper if it
is published in SAE Transactions. For permission to publish this paper in full or in
part, contact the SAE Publications Group.
Persons wishing to submit papers to be considered for presentation or publication through SAE
should send the manuscript or a 300 word abstract of a proposed manuscript to: Secretary,
Engineering Meetings Board, SAE.