Axisymmetric Analysis of Bolted Disc Brake Assembly To Evaluate Thermal Stresses
Axisymmetric Analysis of Bolted Disc Brake Assembly To Evaluate Thermal Stresses
Axisymmetric Analysis of Bolted Disc Brake Assembly To Evaluate Thermal Stresses
Abstract:
Typically thermo-mechanical analysis including complexities such as contacts and bolt preloads are
carried out using three dimensional models. These analyses require significant time and effort in FE model
building, analysis setup, solution, and results processing. It also requires special effort to ensure it is error
free.
In order to get stable and accurate results element size and time step selection is very important in
transient analysis. These aspects are discussed in this paper.
This paper also talks about simplified yet almost equally accurate modeling and analysis method for
thermo-mechanical analysis using brake fade test simulation as an example. This methodology is based on
use of ABAQUS Axisymmetric analysis technique modified to represent effect of discrete bolting, bolt
preloads, and contacts within various components of the assembly.
Analysis results as well as analysis turnaround times are compared between this new method and the
conventional method. Up to 80% time can be saved with significant improvement in the accuracy of the
results.
Repetitive braking of the vehicle leads to heat generation during each braking event. The resulting rise in
temperatures has very significant role in the performance of the braking system. Problems such as
premature wear of brake pads and thermal cracking of brake discs are attributed to high temperatures.
Consequently controlling the temperature profiles and thermo-mechanical stresses are critical to proper
functioning of the braking system. CAE simulations are often used to for evaluating the brake disc design
using thermo-mechanical analysis techniques. Conventional approach is to use three dimensional FE
models of the brake discs. This approach has major drawbacks of higher pre and post processing as well as
solution times. Need is felt to develop a quick and reliable method to evaluate the thermal stresses in brake
discs. This paper describes one such approach based on modified axisymmetry analysis.
Figure 1: Brake thermal analysis process for a vehicle under a given braking schedule
First, a transient heat transfer analysis is carried out in order to predict the temperature distribution over the
disc during the brake fade test (repeated acceleration and braking events). Subsequently, a mechanical
analysis is performed introducing the computed temperature as thermal loads and considering the bolt
pretension.
A need was felt to reduce the analysis cycle times through development of alternate simulation
methodology. This has been achieved through using axisymmetry analysis method modified to include
representation of discretely bolted assembly. It also incorporates bolt preloading and contacts.
Figure 2: FE models of a disc brake system used in the modified axisymmetry analysis.
The same FE mesh has been used for the both thermal and mechanical analysis. For the transient thermal
analysis, second order CAX3T elements are used for discretization of the disc and hub whereas CAX6M
(second order element with modified integration) are used for the structural analysis. The bolts are modeled
as CPS6 (second order plane stress elements). Matching mesh pattern is required between the mating
components at the interface.
In order to get stable and accurate results element size and time step has been selected using Time
Integration method which is relation between thermal diffusivity and minimum time.
∆x = [4 * (a * tmin) 0.5] / N, where ‘a’ is thermal diffusivity and tmin is minimum time required to model
physical event accurately.
To account for all eight bolts used in the joint, the combined cross-sectional areas of the shank and the head
of the bolts is calculated and redistributed to the bolt mesh appropriately using the area attributes for the
solid elements.
Following procedure is used to calculate effective bolt element layer thicknesses, contact surfaces involved
in interaction between plain stress bolt elements and axisymmetric disc elements and material properties for
hole region in axisymmetric model.
Figure 3 below shows the schematic representation of discretization of the bolt using plane stress elements.
Each plane stress element represents a volume that extends out of the – plane. Element A represents a
volume calculated as (HA) x (Area A). Likewise, element B represents a volume calculated as (HB) x
(Area B).The sectional area in the – plane pertaining to a given element can be calculated as
Where R is the bolt head radius or the shank radius, depending on the element location, and X1 and X2 are
-coordinates of the left and right side of the given element, respectively. This sectional area is divided by
the respective element width, and then we obtain representative element thickness. Multiplying each
element thickness by the number of bolts in the model gives the effective element thicknesses.
Sectional areas that are associated with bolt head elements located on the model's contact surfaces are used
to calculate the surface areas of the nodes used in defining the node-based surfaces of the model.
A the presence of the bolt holes in Axisymmetric disc brake system is represented by assigning smear
material properties in the bolt hole area of the mesh and using inhomogeneous material properties that
correspond to a weaker material in this region.
Thermal Analysis:
In this thermal analysis a time-dependent heat source is activated for 15 braking events. The heat flux, heat
transfer coefficients, test conditions such as initial operating temperature of the disc and duty cycle define
the loading conditions for this analysis.
A transient heat transfer solution is used to predict the temperature distribution under these loading
conditions. The comparison of temperature results between “Conventional approach” and “Modified
axisymmetry approach “in the disc is shown in Figure 4.
Conclusion
• The modified axisymmetry analysis used is significantly faster. The time to process a typical brake
disc analysis has been reduces substantially from several hours to minutes. Upto 80-85% savings
in time are seen.
• It is shown that even the bolt stresses predicted in axisymmetric method are accurate.
2007 Abaqus India Regional Users’ Meet Page 6 of 7
Acknowledgments
We would like to thank Mr. Ashok Joshi (Practice Head) and Mr. D.P.Salunkhe (Project Manager) for
giving us the opportunity to work on this challenging project and to present this paper at the ABAQUS
INDIA user conference.
We would also like to thank co-workers involved in the vehicle testing group and vehicle design group for
giving us the desired technical inputs at the appropriate time of the project.
References
2. SAE - 2005- 2006-01-3558: FE Prediction of Thermal Performance and Stresses in a Disc Brake
System by Amol Apte and H Ravi.