Portfolio

What is coupled thermomechanical analysis in disc brakes?

Coupled thermomechanical analysis is an advanced finite element workflow that simultaneously calculates structural deformation and frictional heat generation during a braking event. At BanuMusa R&D, our Finite Element Analysis (FEA) Consultancy team specializes in simulating these extreme transient events. This case study details our advanced simulation methodology for the coupled thermomechanical analysis of disc and brake pads using Dassault Systèmes SIMULIA Abaqus, aiming to predict thermal fatigue, optimize heat dissipation, and prevent catastrophic structural failure.

Client & Industry Background

In the automotive and heavy machinery sectors, the braking system operates as the most critical safety mechanism. During emergency or harsh braking events, the massive kinetic energy of a vehicle is converted almost entirely into thermal energy through dry friction between the brake pads and the rotor disc. Our client, a heavy-machinery component manufacturer, required a highly accurate predictive model to evaluate rotor survivability under extreme, repeated braking loads.

The Engineering Challenge: Focal Hot Spots & Phase Transformation

The primary failure mechanism in heavy-duty disc brakes is the generation of localized thermal gradients. When non-uniform contact occurs between the brake pad and the disc—often caused by high local pressure and a high rate of heat flux generation—Critical Focal Hot Spots begin to form.

When these hot spots become macroscopic, the localized surface temperature can rise exponentially, often exceeding 1100°C. At this extreme thermal threshold, the cast iron surface of the disc undergoes a detrimental phase change, transforming into brittle martensite upon rapid cooling (Also, we built an in-house Abaqus UMAT subroutine for phase transformation that you can check it here). This microstructural phase change has a severely negative effect, acting as an initiation site for thermal cracks that propagate rapidly through the cast iron matrix.

To avoid this, engineers must optimize the structural design to reduce the maximum surface temperature and ensure a uniform contact pressure distribution.

FEA Simulation Methodology in Abaqus

To accurately capture the interaction between structural deformation and heat generation, our engineers utilize the fully coupled temperature-displacement solvers in Abaqus/Standard.

Frictional Heat Flux Generation

The core of the simulation relies on accurately converting frictional work into heat. The rate of frictional heat generation (q) at the sliding interface is calculated using the standard formulation:

q = η * μ * P * V

Where:

• η is the fraction of dissipated energy converted into heat (typically 0.9 or 90%).

• μ is the temperature-dependent coefficient of friction.

• P is the local contact pressure.

• V is the relative sliding velocity between the pad and the disc.

Material Calibration & Thermal Properties

Various industrial studies have shown that the coefficient of friction, specific heat, and thermal conductivity have the greatest effect on contact temperature prediction. If the material’s Young’s modulus is too low, the pad may deform excessively; if it is too high, stress concentrations will form.

To eliminate trial-and-error in defining these temperature-dependent properties, our engineers rely on the MatLib for Abaqus Material Library. MatLib provides pre-calibrated, unit-consistent thermal-mechanical data for various grades of cast irons, composite friction materials, and steel alloys, ensuring high-fidelity simulation results.

Validation & Transient Thermal Response

A crucial aspect of our Engineering as a Service (EaaS) workflow is rigorous Verification and Validation (V&V). The graph above demonstrates the transient thermal response of the brake disc at a specific radial distance (r = 110 mm) during an emergency braking cycle.

1. High-Fidelity Validation: The Abaqus simulation results (blue line) show an exceptional correlation with the established analytical/experimental benchmark by Gao (2002) (orange line). This validates our methodology regarding boundary conditions and convection coefficients.

2. Sawtooth Heat Flux Pattern: The distinct zigzag fluctuations are not numerical errors; they represent the true physics of rotary frictional heating. Every time the brake pad passes over the tracking node (r = 110 mm), a sudden heat flux shock causes a rapid temperature spike. As the pad moves away, the exposed surface cools briefly via air convection, creating the sawtooth pattern until the vehicle comes to a complete stop.

3. Martensite Prevention: Our analysis confirms that the peak focal hot spot temperature reaches approximately 280°C. This provides a massive safety margin against the 1100°C martensitic phase transformation threshold, ensuring the disc remains ductile and free from thermal cracking.

Conclusion & Fatigue Life Prediction

The thermomechanical coupling of disc and brake pads is a highly non-linear problem that cannot be solved using traditional hand calculations. By implementing this advanced FEA methodology, BanuMusa R&D successfully determined the low-cycle thermal fatigue life of the disc, identified the exact nodes prone to crack growth, and optimized the brake system design for high-performance racing applications.

Simulation vs. Reality: Disc Brake Validation

Watch our side-by-side comparison of physical testing and Abaqus FEA results for thermal coning and frictional heat flux.