Finite Element Analysis of Disc Brake Considering Heat

Flow

Noman Bashir Muhammad Azhar Saeed

Department of Mechatronics Engineering Department of Mechatronics Engineering
National University of Sciences & Technology National University of Sciences & Technology
Peshawar Road, Rawalpindi Pakistan Peshawar Road, Rawalpindi Pakistan
nomanbashir28@gmail.com mazharsaeed166@gmail.com
Abstract—This research aims to provide a better performance discs are less susceptible to the "brake fade”. Disc brakes are
to the most critical part of an automotive system namely brake. now getting popular among manufacturers as it provides better
Brakes are responsible for slowing down of vehicle by dissipating stopping power and are more self-adjusting than drum brakes.
the kinetic energy. When the brake system is brought into action The basic mechanism of disk brakes goes like when the
a significant amount of heat runs and get dissipated through the
brake paddle is pressed the hydraulic fluid fills the caliper with
system. For covering all these aspects designers use number of
pressure causing the piston to move in outward direction
procedures/experiments to choose suitable material that provides
a data about temperature distribution, stresses, contact pressure,
making brake pads press on rotor that is in motion as it has also
strain distribution, deformation, shear stress and strain. We aim been attached to wheel via Hub.
to use simulations on ANSYS for the above mentioned aspects Joe Y.G. cha explained disc brake instability analysis
with different materials and predict the most preferable from because of the oscillation produced in course of braking action
these results. and how they flow causing certain stresses in system.
Thermal management contributes to a large extent to the
Keywords—ANSYS; SolidWorks; Finite Element Analysis; wear of the rotor and the absolute braking performance. Every
Disc Brakes; time the brakes are used, the friction between the brake pads
that rub against the rotor creates incredible amounts of heat.
I. INTRODUCTION If the temperature rises too high, there is a chance that a
Brakes are somewhat the most delicate and sensitive part of thin layer of brake pads will form on the surface of the rotor,
any automobile. As the technology progresses the automobile causing the brakes to fade. The braking power of the vehicle is
manufacturers are getting in a race to get faster and faster reduced once the brake is affected by high temperatures.
vehicles on roads. These advancements in line require a similar
progressive braking technology and mechanism for vehicle.
Every braking system tends to meet the following
requirements:
 The braking force should be high enough to stop
vehicle in minimum time and distance.
 Provide maximum control on vehicle during
braking action.
 The system should be made of anti-fade materials
whose friction characteristics do not drop/fade.
 Brakes should not get really heated up that may
cause system to malfunction
Currently the most efficient and most widely used braking
technology is disc brakes system. Following are components
of a Disc Braking System: Figure 1-How a Disc Brake Works (7)
 Brake pads
 Caliper
 Rotor
In 1902, Frederick William Lanchester patented the first
caliper type automobile disc brake. He used it successfully on
Lanchester cars. In comparison with earlier used drum brakes,
disc brakes provides better braking action, due to the fact that
the disc cools down at more rapid rate. As a consequence,
 Numerical solution for a variety of engineering problems material of the pads as well as disc which generates heat.
can be found by using FEM. Any complex geometry made of Since this generated heat is very small, so we are not taking it
any material can be handled by using this method. under into account. The friction heat flux generated in the disc-pad
different conditions. The general finite element method interface is as follows [10].
satisfies the requirement of analysis of the complex design and
design systems of today, in which the equilibrium equation q (x,y,t) = µP(x,y) * v(t)= µP(x,y)ɷ(t).r (1)
solution is normally unusable. In finite element method,
“nodes” are formed. These nodes form the design. Connected Where is the co-efficient of friction, v is the velocity of
with the nodes are finite elements that form the network of sliding that depends on angular velocity ɷ and radius r of the
finite elements that are related to structural and material disc. P is the pressure of contact. Sum of the heat generated in
features of the geometry. These elements determine the disc qD and that generated in the pads q P is equal to total heat
behavior of the model under certain conditions. Finite element produced q. γ is the relative braking energy which is absorbed
density which is known as mesh, is different at different points by the disc.
in the material. It depends on stress level in the given area.
Areas of higher load variations generally needs higher mesh γ = qD/q = qD/(qD+qP) =1/1+(ρPcPkP/ρDcDkD) (2)
density than with little or least load Variation.
Where q is the density, c is the specific heat and k is the
A. Problem
thermal conductivity. P and D in suffix is for pads and disc
The main problem with the disc brake is that when brake respectively.
is applied, a frictional force arises between disc and pads. This
frictional force causes temperature variation from 20oC to
more than 700oC within a few seconds. This paper presents a
three-dimensional Finite Element Model (FEM) of the disc
and two pads.
B. Finite Element Model
In FE Model, the first step was the preparation of a
structural model of disc and pads. ANSYS, a finite element
software, was used to perform this task. Next step was to
initialize the calculations in ANSYS. For this purpose, certain
boundary conditions were defined, and meshing was done on
Figure 2 -Disc Pad Isolated from Wheel Assembly
finite element model. Our finite element model comprising of
two pads and a single disc is shown in the figure 7. Area of
contact of disc and pads is illustrated I n figure 8. Properties D. Material Properties
of mesh are given the table 1. Three different materials were The materials mentioned in table 2 will be used in the FEM
selected for the disc which include Carbon-Ceramic, Grey analysis and prediction for the best possible material.
Cast Iron and Stainless-Steel Grade 316. Properties of each of
these materials are given in the table 2.
E. Cast Iron
C. Thermal Analysis
Cast iron usually refers to grey cast iron, but identifies a large
Disc-brake system of a car separated from wheel assembly group of ferrous alloys, which solidify with a eutectic. Iron
is shown in the figure 1. When brake is applied, temperature accounts for more than 95%, while the main alloying elements
rises at the contact of disc and pads and heat thus produced, is are carbon and silicon. The amount of carbon in cast iron is
dissipated into the disc and pads. Heat that gets wasted into the range 2.1-4%, as ferrous alloys with less are denoted
the pads can be calculated correctly by estimating that carbon steel by definition. Cast irons contain appreciable
produced at disc-pad interface. Heat generated at the interface amounts of silicon, normally 1-3%, and consequently these
is equal to heat dissipated into disc and pads. Because of very alloys should be considered ternary Fe-C-Si alloys. Here
high angular speed, so we can ignore temperature gradient in graphite is present in the form of flakes. Disc brake discs are
theta-direction. Based on this, we can assume that heat flux is commonly manufactured out of a material called grey cast
applied on whole circumferential domain. This method, we iron. [15]. for cast iron,
used, is approved by experimental investigations and widely  Young’s modulus (E) = 125 GPa
used in analytical studies.  Poisson’s ratio (v) =0.25

Vehicle in motion basically, posses two kinds of energies that

is kinetic energy and potential energy. On application of
brakes, these energies are converted into thermal energies at
disc-pad interface. There is also slight wear and tear of the
 Table 1- Material Properties
Model Generations: Material Yield Tensile Elastic Shear
Strength Strength Modulus Modulus
Two different approaches can be used to generate a the N/m2 N/m2 N/m2 N/m2
geometry for subject analysis: Al Alloy 1.65 3 7e+011 3.1e+010
 Direct generation (Modeling in ANSYS [10]
 Solid modeling Stainless 290e+6 500e+6 205e+9 152e+9
We first designed the model in Solidworks and saved it with Steel
.IGS extension and later imported the geometry in ANSYS Grade
F. Mesh Model 316
The model is first made in SolidWorks ® and then a mesh [9]
network was generated in ANSYS 18.1 Workbench The finite Grey 98e+6 250e+6 105e+9 69e+9
element model of the rotor is carried out with a mesh of 20351 Cast
elements for a total of 39208 nodes, where the mesh of the Iron
disc and pad resulting from ANSYS software is presented in [8]
Figure 13. Rotors are made of cast iron for following three reasons [13]:
& It is relatively hard and resists wear. & It is cheaper than
steel or aluminum. & It absorbs and dissipates heat well to
cool the brakes

Mainly three types of mechanical stress is subjected on disc

brake.  Traction force, caused by centrifugal effect and it
occurs when wheel is rotating and no brake force is applied to
the disc.  Compressive force, when the brake is applied due
to action of the force, exerted by pressing the pad
perpendicular onto the surface of the disc.  Due to braking
Figure 3-Mesh Network action caused by rubbing on the brake pad against the surface
of the disc. It acts in opposite direction of the disc rotation.[5]

.[5]
ANALYSIS:
Ali. Belhocine, Abd. Rahim Abu Bakar , Mostefa Bouchetara,
Ansys workbench provides with a large range of analysis that
Structural and contact analysis of disc brake assembly during
can be performed on a certain geometry. Following are the
single stop braking event, American Journal of Mechanics and
steps needed to find a solution or result of any type of design
Applications.
on ANSYS.
1. Open ANSYS Workbench.
2. Select the type of analysis that you want to perform.
3. Choose the material from Engineering Data sources P = Fdisc÷ (As*µ) = 10000 Pa (3)
that you need your geometry to be made of.
4. Then you need to import the required geometry file In equation (3) Fdisc is the force that is being applied on the
on which the analysis is going to be performed. rotor disc. As is the surface area in contact with the brake pads
5. After Importing the file . you need to click on Setup, and µ is the coefficient of the friction. So the total
This will open Mechanical [ANSYS Mechanical Pressure due to the brake pad pressing on the rotor is
Enterprise]. considered to be 1 KPa.
6. After this generate suitable mesh and apply suitable
conditions and constraints. 6.1.
7. Last step is to solve for obtaining results.

Large deformations are occurred at the outer radius of the disc.

The smoothness of the mesh increases the precision of the
solution.
Ventilation provoides a better help for the rotor to let it keep
itself cooler and temperature affects get decreased .
But too many ventilation holes will make the rotor weak and
certain chance/ probability to get fractured.
 Engineering Research and Development (IJMPERD) ISSN(P): 2249-
6890; ISSN(E): 2249-8001 Vol. 4, Issue 1, Feb 2014, 95-112
[7] HowStuffWorks.com<https://auto.howstuffworks.com/auto-
parts/brakes/brake-type/disc-brake.htm> 14 April 2018
[8] http://www.makeitfrom.com/material-properties/Grey_Cast_Iron
asm.matweb.com/search/SpecificMaterial.asp?bassnum=mq316a10
REFERENCES
Ffd
[9] Newcomb, Transient temp attain on disk brakes, Br. J. Applied Physics.
[1] Rakesh Jaiswal, Anupam Raj Jha, Anush Karki, Debayan Das, Pawan 10 (1959) 339e340.
Jaiswal, Saurav Rajgadia and Ankit “Structural And Thermal Analysis [10] D. Majcherczak, P. Dufrenoy, Y. Berthier, Tribologically, thermally and
Of Disc Brake Using Solidworks And Anysys” International Journal of mechanical coupling aspects of the dry sliding contact, Tribool.
Mechanical Engineering and Technology (IJMET) Volume 7, Issue 1, International. 40 (2007) 834e843001.
Jan-Feb 2016, pp. 67-77, Article ID: IJMET_07_01_008 [11] http://scholar.google.com/scholar?
hl=en&as_sdt=0%2c5&q=strain+in+the+material+of+disc+brake+on+b
[2] Abhishek Kumar Tiwari, Akhilesh Kumar Tiwari, Pramod Yadav,
raking&btnG=#d=gs_qabs&p=&u=%23p%3DkIOX-11varQJ
Harigovind Singh Yadav, Shyam Bihari Lal” Finite Element Analysis of
Disc Brake for Aluminum Alloys” International Journal of Scientific &
Engineering Research, Volume 5, Issue 4, April-2014, ISSN 2229-5518 [15] V. M. M. Thilak, R. Krishnaraj, Dr. M. Sakthivel, K. Kanthavel, Deepan
[3] A. Belhocine , O. I. Abdullah, “ Finite Element Analysis of Automotive Marudachalam and M.G, R.Palani, (2011),” Transient Thermal and Structural
Disc Brake and Pad in Frictional Model Contact”, Int J Advanced Analysis of the Rotor Disc of Disc Brake”, International Journal of Scientific
Design and Manufacturing Technology, Vol. 7/ No. 4/ December – 2014 & Engineering Research Volume 2, Issue 8 .
[4] Vineet Kumar Dwivedi, Arvind Singh,” FEA Analysis of Brake using
Ansys - A Review” , 2017 IJSRST | Issue 8 | Print ISSN:2395-6011 |
Online ISSN: 2395-602X
[5] Ali Belhocine ,”FE prediction of thermal performance and stresses in an
automotive disc braking”, INTL J Adv Manuf Technol (2017) :3563-
3578 DOI 10.1007/s001 70-016-9357-y, Published Online : 23 August
2016 © Springer- Verlag London 2016
[6] Mahmood hasan dakhil, a. K. Rai, p. Ravinder Reddy & Ahmed
Abdulhussein Jabbar, “Design And Structural Analysis Of Disc Brake In
Automobiles”,International Journal of Mechanical and Production