Design and Simulate An Aerodynamic Urban Concept Car Body For The Shell Eco-Marathon
Design and Simulate An Aerodynamic Urban Concept Car Body For The Shell Eco-Marathon
Design and Simulate An Aerodynamic Urban Concept Car Body For The Shell Eco-Marathon
URABNCONCEPT CAR
BODY FOR THE SHELL ECO-MARATHON WITH
LESS COEFFICIENT OF DRAG
A Thesis Proposal
Presented to the
School of Mechanical and Manufacturing Engineering
Mapua Institute of Technology
In Partial Fulfillment
of the Requirements for the Degree of
Bachelor of Science in Mechanical Engineering
By:
Da Silva, Elisario M.
2009108764
Diwa II, Jose S.
2009151238
Pimentel, Marc Desie D.
2010100349
December 2014
TABLE OF CONTENTS
Bibliography .....................................................................................................................70
LIST OF FIGURES
Figure 2-1 Effect of the wind angle to the drag and lift coefficient ..........................6
Figure 2-2 Center Pressure (Anderson 2001) .............................................................9
Figure 2-3 Example of Flow Similarity (Anderson 2001) ........................................10
Figure 2-4 Division of Flow (Anderson 2001) ...........................................................12
Figure 2-5 Subsonic Flow (Anderson 2001) ..............................................................13
Figure 3-1 Bernoullis Principle .................................................................................16
Figure 3-2 Bernoullis Principle on wing with Air as fluid......................................13
Figure 3-3 Drag Force generated by motion and air flow .......................................18
Figure 3-4 Various shapes with coefficient of drag ..................................................20
Figure 3-5 Down force on an automobile ..................................................................21
Figure 3-6 Rolling resistance on automobile tire ......................................................23
Figure 3-7 Wind Speed graph in Manila for the month of February .....................25
Figure 4-1 Process Flowchart of the study ................................................................26
Figure 4-2 HABAGAT ................................................................................................28
Figure 4-3 TVR TUSCAN ..........................................................................................28
Figure 4-4 Haribon CAD using Autodesk Inventor .................................................29
Figure 4-5 Flow Simulation of Haribon using Autodesk Flow Design ...................30
Figure 4-6 Habagat CAD using Autodesk Inventor .................................................31
Figure 4-7 Flow simulation of Habagat using Autodesk Flow Design....................32
LIST OF TABLES
ii
ACKNOWLEDGEMENT
This thesis is the end of our journey in obtaining our BS degree in Mechanical
Engineering. We have not travelled alone in this journey. The contributions of many different
people, in their different ways, have made this study possible.
Foremost, we would like to extend our appreciation and thank God for the wisdom and
perseverance that He has bestowed upon us throughout this study.
Moreover, we would like to grab this opportunity to express our gratitude to everyone
who supported us throughout this journey. We are thankful for their aspiring guidance,
invaluably constructive criticism and indeed friendly advice during the study. We are sincerely
thankful to them for sharing their truthful and illuminating views on a number of issues related to
this thesis project.
We would like to express our sincere gratitude to our advisor Prof. Sherwin S. Magon for
the continues support in this research, for his patience, motivations and immense knowledge.
This work would have not been possible without his guidance, support and encouragement.
Under his guidance we overcame many difficulties and learned a lot.
Besides our advisor, we would also like to thank our panels; Engr. Igmedio Isla, Engr. Jaime
Honra and Engr. Hans Felix Bosshard for their helpful criticisms, insightful comments and
valuable suggestions.
Last but not the least, we would like to thank our families especially our parents for
supporting us spiritually throughout our lifes.
iii
ABSTRACT
iv
CHAPTER 1
INTRODUCTION
This is a study to design and simulate an aerodynamic urbanconcept exterior body for Shell
Eco- Marathon Asia 2013 to attain low coefficient of drag and help the vehicle attain less fuel
consumption. This chapter discusses the background, statement of the problem, objectives,
significance and scope and limitations of the study.
1.1 Background of the study
The Royal Dutch Shell sponsors an annual competition about the future of
transportation and mobility titled Shell Eco-marathon. It is a worldwide competition
wherein students are challenged to design, build and drive the most energy-efficient car.
The competition is focused on maximizing a vehicles mileage with a given amount of
fuel rather than achieving high speeds. The winner of the competition would be the
vehicle that could travel the most distance given one liter of fuel.
This is the third time that the Mapa Institute of Technology School of Mechanical
and Manufacturing Engineering will be competing for the challenge under the
urbanconcept category. Urbanconcept is a prototype car that contains all the features of
todays commercially available cars. The institutes first urbanconcept entry in 2011 was
named Habagat. Due to some constraints and difficulties, it was unable to finish the
race.
CHAPTER 2
REVIEW OF RELATED LITERATURE
The presentation of both thesis and journal articles in this part shows to expound further how
aerodynamic concepts can significantly affect the output of a vehicle.
2.1 Review of Related Studies
The presentation of both thesis and journal articles in this part shows to expound
further how aerodynamic concepts can significantly affect the output of a vehicle.
2.1.1 Computational Study of Flow Around a Car body ( Z. Zheng 2009).
Z.Zheng made a study about the flow around a car body with different relative
angle wind to car axis, varying it from 0 to 30 as shown in figure 2-1. The scale of
the automobile he uses is 435 mm x 168 mm x 148 mm (length x width x height) with
360,000 mesh size. In this study, the drag and lift coefficients are used in order to
compare the different effects of the front window angle on the performance of the
vehicle.
F
Figure 2-1. Effect of the wind angle to the drag and lift coefficient (Zheng 2009)
6
2.1.2 Design and Construction of the Urban-concept Car Exterior for Shell Ecomarathon Asia 2011 (M. Bernabe et. al, 2011)
Previous vehicle entry of Mapa for Shell Eco-Marathon Asia urbanconcept car
made a research about its exterior body. The researcher design and fabricate the
exterior body of an urbanconcept that compete in Malaysia on 2011. The design of
the exterior body has average coefficient of drag of 9.7. Tear drop shape of an
exterior body was highly recommended by the researchers since it only has 0.05
coefficient of drag.
2.1.2 Simulation and analysis of drag and lift coefficient between Sedan and
Hatcback car (Salleh 2009)
The thesis showed difference coefficients barriers and lift coefficients for two
basic design types of sedan and hatchback cars. Process simulation and analysis for
14 both the model design was conducted with computer-aided drawing software and
analyzed using COSMOS Floworks software. Restriction coefficient and lift
coefficient for the hatchback design is lower than the sedan design. With this
observation, the hatchback design becomes more efficient and its aerodynamic value
increases.
2.2.2 Aerodynamics
The dynamic of gases, especially atmospheric interactions with moving objects.
(Anderson 2001)
2.2.2.1 Dynamic Pressure
Dynamic pressure is a property of moving gas. It is also used in the lift
coefficient and drag coefficient theory. Dynamic pressure is directly
proportional to density and velocity which have a unit of pressure. (Anderson
2001).
2.2.2.2 Center of Pressure
Center of pressure is located where the resultant of a distributed load
effectively acts on the body. Center of pressure is just like center of gravity
where center of gravity is located on a body due to weight while center of
pressure is located on a body due to pressure. When moment is applied to the
center of pressure the acting distributed load will be equal to zero. (Anderson
2001)
11
12
13
CHAPTER 3
THEORETICAL CONSIDERATION
This chapter focuses on the concepts involved in an aerodynamic design for a vehicle
exterior body. The vehicle body plays an important role in increasing the efficiency of a vehicle
since it is in primary contact with the air flow. The parameters that are set by Shell EcoMarathon Challenge in the design of the vehicle exterior body are also discussed in this chapter.
3.1 Shell Eco-marathon Official Rules
The primary consideration in the design of the vehicle exterior body should comply
with the rules and regulations of the Shell Eco-Marathon (SEM). There are set of
specifications for the vehicle dimensions and safety issues that should be followed.
3.2 Aerodynamics
The aerodynamics is one of the important parameters in the design and development
of automobile. Under road conditions, an automobiles performance is affected by air
resistance. When a car is in motion, air flow in opposite direction and in some cases it
flows across the vehicle. Therefore, to be able to overcome this resisting flow of air, a car
should be driven with greater amount of power to move forward.
3.2.1 Bernoullis Principle
The basic principle in consideration for aerodynamic design is the Bernoullis
principle. This theory states that as the fluid velocity increase, the pressure exerted by
that fluid decreases, various flows can be examined using this theory. The Bernoullis
14
equation will be used in the calculation of fluid flow. Also, the Bernoullis principle
came from Thermodynamics.
3.2.1.1 Bernoullis Equation
Energy Balance, Assuming zero friction loss
From Energy Balance:
Energy In = Energy Out
(Eq. 3.1)
1 +1 +1 +1 =2 +2 +2 +2
(Eq. 3.2)
(2 1 2 +1 + +1 = 2 2 2 +2 + +2 )
1 2
2
1
2
+1 + +
= 2 2 2 +2 + +
1 2
+1 + 1 =
2
1 2
2
2 2
2
+2 + 2
+1 +1 =
2 2
2
+2+2
(Eq. 3.3)
(Eq. 3.4)
(Eq. 3.5)
(Eq. 3.6)
Note: Force (F) is equal to the product of Pressure (P) and Area (A). There
is also not change in temperature (T=0)
Equation 3.5 and 3.6 are called Bernoullis equation. Either of the two
equations can be used to calculate problems in fluid dynamics.
The second law of thermodynamics states that when energy is transferred,
that energy cannot be conserved and some energy must be reduced to some lower
value. From the second law of thermodynamics, the Bernoullis equation (Eq.
3.5) is modified in consideration of friction losses ( ).
1 2
+1 + 1 =
2
2 2
2
+2 + 2 +
(Eq. 3.7)
15
2 2 1 2
2
(Eq. 3.8)
16
From the law of mass conservation and continuity equation, mass does not
change with respect to position, nor is it affected by pressure, temperature or
motion.
1 = 2
(Eq. 3.9)
1 1 1 = 2 2 2
(Eq. 3.10)
1 = 2
(Eq. 3.11)
Q=
(Eq. 3.12)
(Eq. 3.13)
(Eq. 3.14)
17
Where:
= Frontal Pressure
= Drag Force
= Frontal Area
3.2.3 Drag Force
Drag Force is a type of force that is relevant to motion. This force is generated by
interaction between a solid body and a fluid. There are two requirements to generate
drag. First, the body should be in contact with the fluid to generate drag. Since drag
force is acting between a body and fluid, no fluid means zero drag. Second, there
should be a difference in velocity between the body and the fluid.
18
(Eq. 3.15)
Where:
=
12
2
= 2
(Eq. 3.16)
Where:
= Drag Force
= Mass density of the fluid
19
20
depends on the function of the vehicle under design. The equation below shows the
formula for down force:
F = 0.5* * * * 2
(Eq. 3.17)
Where:
= Total Area of a vehicle
= Coefficient of Drag
= Density of air depending on temperature and altitude
2 = Speed of the car
21
= KE
1
= 2 2
1
= 2 2
1
= 2 [( )( )()( 3 )]
(Eq. 3.18)
(Eq. 3.19)
(Eq. 3.20)
(Eq. 3.21)
Where:
= Energy absorbed that cause energy loss and drag
KE = Kinetic Energy in the air transferred in the vehicle
= Mass of air
v = Velocity of the air relative to the vehicle
= Frontal Area of the vehicle
D = Distance travelled by the vehicle
= Coefficient of Drag
22
(Eq. 3.22)
Where:
= Power needed to overcome rolling resistance
= Rolling coefficient of the tires
m = Mass of the body
g = Gravitational Acceleration
v = Velocity of the car
(Eq. 3.23)
Where:
= Power save due to reduction of weight
m = Mass of the Body
g = Gravitational Acceleration
v = Velocity of the Car
3.6 Resistance due to Inertia
Inertia resistance is defined by Mass and Inertia. Inertia is difficulty to change the
momentum of a body. The relationship of momentum, mass and velocity is given by:
Momentum (p) = Mass (m) x Velocity (v)
(Eq. 3.24)
Also, From Newtons Law of motion mass (m) can be derived from:
Force (F) = Mass (m) x Acceleration (a)
(Eq. 3.25)
24
Figure 3-7 Wind Speed graph in Manila for the month of February
(Source: https://weatherspark.com/averages/33313/2/Metro-Manila-Philippines)
25
CHAPTER 4
METHODOLOGY
This chapter discussed the methods and procedures of designing and fabricating the exterior
body of urbantype vehicle entry of Mapa Institute of Technology for Shell Eco-Marathon Asia
2013. This chapter also discussed the different software used to help the researcher to conduct
the study. The process flowchart of the study is in figure 4-1:
Simulation of Aerodynamic
Characteristic through Computer
Software
28
29
31
32
Chapter 5
DISCUSSION AND ANALYSIS OF RESULTS
This chapter shows the results of the simulation and discusses the effect of the results to the
vehicle body.
5.1 Results
5.1.1 Data of Drag Force and Coefficient of Drag were calculated using the
Autodesk Flow Design Simulation Software.
Habagat
Trial
Wind Speed
Coefficient of
Drag Force
Highest Front
(m/s)
Drag
(N)
pressure
(Pa)
7.22
0.47
14.130
26.512
10.56
0.47
30.359
60.388
13.89
0.47
51.602
104.778
33
Haribon
Trial
Wind Speed
Coefficient of
Drag Force
Highest Front
(m/s)
Drag
(N)
Pressure
(Pa)
7.22
0.31
0.074
11.374
10.56
0.30
0.155
22.357
13.89
0.30
0.272
39.983
CAD software is used in determining the vehicle total frontal area and Flow Design is
utilized to calculate the drag forced given the wind speed. Highest front pressure is
calculated as well. The coefficient of drag is determined by the shape of the design. Drag
force is affected mainly by velocity. Increasing velocity yields higher value for drag force
and higher value for frontal pressure.
35
Chapter 6
CONCLUSION AND RECOMMENDATION
This chapter presents the conclusion of the study and recommendations future improvement
of studies that will be conducted.
6.1 Conclusion
To be able to compete, the Shell Eco-marathon 2013 rules and regulations should be
considered. The parameters in the design and the final dimensions are shown on appendix
K.
Frontal area should be reduced if one is to reduce the value of drag coefficient. The
new urbanconcept body for Haribon, has less frontal area as compared to Habagat. The
reduction in front area is 10.59%. Computation for reduction in area is shown in appendix
L.
The shape of an object affects the numerical value of the coefficient of drag. The
frontal area of the vehicle body is the reference in the simulation of the coefficient of
drag. Minimizing frontal area is the key in reducing drag coefficient. The computed
coefficient of drag via simulation is less compared to commercial cars since its size is
minimized and the body is streamlined.
Drag force, meanwhile is also reduced with minimized frontal area of the design.
Although drag coefficient is determined by an objects shape, drag force can still change
36
depending on the numerical value of the velocity. From the data results, it is clear that
increasing velocity would yield increasing value for the drag force.
6.2 Recommendation
Detailed designing and accurate molding of the vehicle body is a must since the shape
of an object greatly affects the coefficient of drag. The design should be conceptualized
based on existing aerodynamic models such as current trend on car shapes streamlining.
A software simulation is a must in designing the vehicle body so that the team can cut in
determining the optimum shape of the vehicle. Simulation gives the advantage of
skipping construction of prototypes and can give the team more time on fabrication.
Time is limited for the team to produce a vehicle for the Shell Eco-marathon
competition; design phase should be immediately started as soon as a new team is formed
so that design details can thoroughly checked and more time can be provided on
fabrication phase.
The fabrication of the vehicle body should also consider light-weight materials since
it also affects vehicle performance. Rules and Regulations of the Shell Eco-marathon
competition should always be taken into account in the design and fabrication of the
vehicle body.
Besides the above given, teamwork and sharing of ideas gives best result in the
design and construction of a good vehicle.
37
APPENDICES
Appendix A. Gannt Chart
Appendix B. Expenses
Description
Fabrication of Mold
Fabrication of Exterior Body
Total Expenses
Price
P 90,000
P 90,000
P 180,000
38
39
40
FRONT VIEW
TOP VIEW
41
SIDE VIEW
REAR VIEW
42
43
44
FRONT VIEW
TOP VIEW
45
SIDE VIEW
REAR VIEW
46
SIDE VIEW
FRONT VIEW
47
TOP VIEW
48
FRONT VIEW
SIDE VIEW
49
TOP VIEW
50
Stop time = 2 sec (at least); say 5 sec stop time per lap.
Total stop time, tstop-total = 5 sec x 10 laps = 50 sec; say 1 min per trial
Total time to finish, ttotal = 29 min (includes stopping time)
Max time for car to finish at track, tmax = 29min 1 min = 28 min
51
1000
12 1
=
=
= 25.71; 26
1
max
28 60
= 50
50 + 26
= 38
2
= 26
1000 1
= 7.22 /
1 3600
= 38
1000 1
= 10.56 /
1 3600
= 50
1000 1
= 13.89 /
1 3600
52
53
54
55
56
Import CAD File, in this study the researcher import locally a CAD file
57
The file imported is saved as a STL; STL is the neutral file format of CAD
On the toggle bar you can edit the XYZ orientation of the model
58
Since the Air Flow is from Left side to right, fix the desired position of the model.
On the toggle bar you can edit the wind speed from 0-100 m/s
59
Different characteristics of the flow lines can be edit on the flow lines settings located at the
toggle bar.
Click the drag plot on the toggle bar to view the drag graph
60
61
62
% Difference =
% Difference =
..
.+.
= 10.59%
63
65
Appendix N. Sponsors
Autodesk
BOSCH
Decal Republic
ELEWELD
66
SHELL
SolidWorks
TIARA
Malabon Diesel
SMART
SUN Cellular
TOYOTA
67
BIBLIOGRAPHY
68