Advanced Manufacturing

Hot stamping simulation of martensite, cycle time and tool cooling before try-out

A press-hardened part gets its strength from cooling in the die, not from the coil. How a coupled thermo-mechanical-metallurgical model, a measured heat transfer coefficient and a conjugate tool-cooling study turn that into strength, cycle time and tool cost you can plan.

Stacks of pressed sheet-metal parts with flanged edges and punched holes, in black and white
Library photograph of stacked pressed sheet-metal parts. Illustrative only: it is not a press-hardened part and not one of our results.

Hot stamping simulation predicts, in one coupled model, how an austenitised boron steel blank cools through transfer, forming and the closed-die quench, how much of it becomes martensite and where, and how long the press must hold. We use it to set tool cooling and cycle time and to give the crash team a believable strength map before the tool is cut.

A press-hardened pillar or rail is bought for its strength, and that strength is not in the coil. It is made in the die, by cooling. When try-out shows a soft corner, a hold time that holds the line below its parts-per-hour target, or a part that will not sit in the gauge, the remedy is in tool steel: channels re-drilled, an insert replaced, a launch date moved. The design-review questions are plain. Will it reach full hardness everywhere? How fast can we run it? Will the cooling channels we drew do their job?

Why is hot stamping simulation harder than cold forming simulation?

Three problems are solved in the same elements at the same time. Thermal: the thin blank loses heat to air during transfer and to the tool once the die closes, ρcp∂T/∂t = ∇·(k∇T) + ΣLi∂fi/∂t, where the last term is the latent heat of the transformations, and heat crosses the interface as q″ = h(Tblank − Ttool). Mechanical: flow stress depends on temperature, strain rate and strain, and the blank is formed while it cools. Metallurgical: the phases that form decide hardness, and they also add strain and heat. Total strain carries thermal and transformation strain and transformation-induced plasticity on top of the elastic and plastic parts, so the cooling path sets distortion as well as strength.

The model has to follow the whole sequence: furnace austenitisation, transfer, forming, quench under closing force, hold and removal. Karbasian and Tekkaya (2010), A review on hot stamping, Journal of Materials Processing Technology 210 reviews the process, and Merklein et al. (2016), Hot stamping of boron steel sheets with tailored properties: a review, Journal of Materials Processing Technology 228 reviews how tailored properties are produced.

How is the martensite fraction predicted in hot stamping?

Whether a point ends up hard depends on whether its own cooling path misses the nose of the diffusional curve on the continuous-cooling transformation (CCT) diagram.

Cooling paths on a continuous-cooling transformation diagram for press hardening A sketch of a continuous-cooling transformation diagram for a boron steel, temperature upward and time on a logarithmic axis to the right. A horizontal line marks where austenite becomes stable at the top, and another marks the martensite start temperature near the bottom. A C-shaped curve bounds the region where ferrite, pearlite and bainite can form by diffusion. Four cooling paths start from the austenitised blank. After a short, slow air-cooling segment during transfer, the fast path in a closed, well-contacted die runs steeply to the left of the C-curve and crosses the martensite start line, giving martensite. The slow path, from poor contact or a wide gap, crosses the C-curve and ends in soft diffusional products. The tailored path, from a heated tool segment, enters the lower flank of the curve and holds at a bainite temperature. The drawing teaches that part strength is decided by whether the cooling path misses the nose, and that the hot-tool route is how a soft zone is made on purpose. Two textbook relations sit below. CONTINUOUS-COOLING TRANSFORMATION (CCT) SKETCH Ac3 Ms austenite ferrite · pearlite · bainite diffusional products, soft martensite, hard TEMPERATURE TIME, LOG SCALE 1 2 3 4 5 6 fm = 1 − exp[ −α (Ms − T) ] diffusional: f = 1 − exp(−k tⁿ), summed over the path by the additivity rule
  1. Austenitised blank and transfer. The blank leaves the furnace fully austenitic and cools slowly in air on its way to the press. The longer the transfer, the lower the temperature at which the die closes and the closer the starting point sits to the nose.
  2. The nose of the diffusional curve. Ferrite, pearlite and bainite form by diffusion, so they need time. The curve bounds where they can start, and its nose is the shortest time at which they can. It moves with chemistry, austenite grain size and any strain held in the austenite.
  3. Fast path: contact in the die. Pressed hard against a cold tool, the blank cools steeply, passes to the left of the nose and reaches the martensite start temperature with austenite intact. Below it, martensite forms by shear, following the Koistinen-Marburger relation.
  4. Slow path: gap or poor contact. Where the blank springs off the tool, meets a wide clearance or sits on scale, the interface coefficient falls, the path bends into the curve and soft diffusional products form. The part is soft exactly where the contact was worst.
  5. Tailored path: a heated tool segment. A tool section held warm, or without cooling, slows the quench on purpose. The path enters the lower flank of the curve and holds at a bainite temperature, giving a softer, more ductile zone in a chosen place.
  6. Ms line: latent heat and strain. Below the martensite start temperature the austenite shears and expands, releasing latent heat and adding a transformation strain that, with transformation-induced plasticity, is what distorts the part and warms it locally.
Three cooling paths from the same austenitised blank: the die decides which one each region takes. Illustrative schematic.

There is no single critical cooling rate to take from a handbook. It depends on chemistry, austenite grain size and austenitising practice, and strain held in the austenite can shift the diffusional transformations, so forming and quench cannot be treated as separate steps. A calibrated model carries the CCT behaviour of your coil instead. Ferrite, pearlite and bainite form by diffusion, with Avrami-type kinetics, f = 1 − exp(−k·tn), where k and n depend on temperature; the additivity rule advances them along a changing temperature by summing the fraction of the incubation time used up until it reaches one. Martensite is displacive: it depends on undercooling, not time, and follows fm = fγ·{1 − exp[−α(Ms − T)]}, the relation of Koistinen and Marburger (1959) in Acta Metallurgica 7, with fγ the austenite left at Ms. Hardness is commonly taken as a rule of mixtures over the phase fractions, calibrated against hardness traverses or tensile tests on a formed part.

Three contour plots of a bent structural part at quench times of 1, 2 and 7 seconds, with the colours changing as a phase transformation spreads through it
Library image of a quench phase-transformation contour on a part at three times. Illustrative only, not one of our results.

Why does the heat transfer coefficient decide the result?

Every prediction above inherits the interface. The coefficient h rises with contact pressure and falls with any gap; it changes with coating or oxide scale and with the surface temperature of the tool. And the blank does not sit flat on the tool everywhere: side walls, radii and thin sections see lower pressure and more clearance. A constant h gives a smooth, optimistic answer, and the slow region is rarely where you expect it. The remedy is to measure it. Instrumented-tool trials, with thermocouples in tool and blank, give h as a function of pressure by inverse analysis, and the model is then checked on a different part or condition than the one it was calibrated on.

Which hot stamping question needs which model
QuestionModelUse it when
Is the blank too cold when the die closes?Thermal model of transfer, with emissivity and oxideTransfer time or furnace-to-press layout is in question
Does every region reach full martensite?Coupled forming and quench, calibrated kinetics, h(p, gap)Part strength or soft spots are the worry
Where do channels go, and how long is the hold?Conjugate tool and coolant model over several cyclesTool design, cycle time and tool cost
Where will it warp, and why does it not fit the gauge?Thermo-mechanical model with transformation strains, through to room temperatureFit, joining and assembly tolerances
Where should a soft zone sit?Tailored tool segments with temperature or hold as parametersThe crash design asks for a ductile region
How will the formed part behave in a crash?Hardness map passed to the crash meshRelease of the crash model for sign-off

How do tool cooling channels set hot stamping cycle time?

After forming, the cycle is mostly waiting. The press holds until the slowest region of the part has cooled through its transformation and is stiff enough to leave the die, and transfer, closing and opening add the rest. Shortening the hold is therefore a tool-cooling problem, and tool cooling is conjugate heat transfer: conduction in the tool steel, convection in the channel, coupled through the interface. Distance from channel to surface and spacing between channels set how evenly the surface cools. Flow rate sets the coolant-side film coefficient, but once that film is not the bottleneck, more flow buys pumping power and little else; the conductivity of the tool steel and the interface decide the rest. Coolant also warms along a series circuit, and the tool itself heats over repeated strokes, so the study must run several cycles. It is the same physics as in our heat exchanger article, applied to a die.

Hot stamping die in cross-section: where the heat goes A cross-section of a closed hot stamping die with a hat-shaped part. A hot blank sits between an upper die and a lower die. Cooling channels are drilled in both dies at a distance below the working surface, with a spacing between neighbours, and coolant flows through them. Heat arrows run from the blank across the contact interface into the tool steel and on to the channels. On the left, a tailored zone replaces the channels with heater cartridges and an insulating slit, so the blank there cools slowly and stays soft. On the right a magnified circle shows the contact interface: two rough surfaces touch at a few points under pressure, so heat passes mainly through the contacts and poorly across the gaps. Below, the interface heat-flux relation and the cycle-time sum. CLOSED DIE, SECTION THROUGH A HAT-SHAPED PART (QUALITATIVE) UPPER DIE LOWER DIE 5 d s 3 1 4 in out 6 INTERFACE, MAGNIFIED p 2 q″ = h(p, gap, coating) · (Tblank − Tdie) h is the dominant unknown: it changes over the part and over the cycle cycle time ≈ transfer + closing + hold until cooled + opening
  1. The hot blank. An austenitic blank is formed and then held under pressure in the closed die. Everything after closing is heat extraction: the blank has to drop through the transformation range fast enough, everywhere, before the die can open.
  2. The contact interface. Two rough surfaces touch at a few points, so heat flows mainly through the contacts. The coefficient h rises with contact pressure and falls with any gap, and a coating or oxide layer adds resistance. It is the least certain input in the whole model.
  3. Cooling channels. Distance from channel to working surface and spacing between channels set how evenly and how fast the tool surface cools: closer channels extract faster but print stripes of temperature into the surface and take away tool material.
  4. Heat path through the tool steel. Heat crosses the interface, conducts through the tool steel and is picked up by the coolant. The tool steel conductivity and the path length are as important as the coolant flow rate.
  5. Tailored zone. Heater cartridges and an insulating slit replace channels in one section, so the blank there cools slowly, transforms to softer phases, and the part gets a ductile zone by design.
  6. Coolant circuit. Coolant warms along its length, so a long series of channels cools the last channel less than the first. The tool and the coolant are one conjugate heat transfer problem, and the cycle time follows from it.
Heat leaves the blank through a few contact points, crosses the tool steel and is carried away by the coolant. Illustrative schematic.

Do not tune channels on a uniform heat transfer coefficient

A layout tuned with one constant h is tuned for a contact that does not exist. Run the cooling study with h as a function of pressure and gap, and judge the layout on the region with the worst contact, not on the average.

Soft zones, coatings and distortion: what else must the model get right?

  • Tailored properties. A heated or uncooled tool segment, or a slower local quench, gives a softer zone on purpose. The model must place the transition between hard and soft where the crash design wants it, and keep it there as the tool warms.
  • Coating and oxidation. A coated blank and a bare blank differ in emissivity during transfer and in interface resistance during the quench, so each needs its own card. Scale adds a resistive layer at the interface.
  • Distortion. Thermal gradients and the expansion at Ms leave residual stress and shape error in the part. Our sheet metal forming and stress-relief and heat-treatment work cover the same chain for springback and residual stress, and the metal forming simulation section places this process among the others.

How does the strength map reach the crash model?

The output the crash team needs is not a colour plot but a property map on their mesh. We map hardness or phase fraction onto the crash model and convert it to yield and tensile strength with a correlation measured on your coil. Rate dependence and failure still need their own tests, as described in our crash simulation credibility checklist, and soft-zone transitions are the places to check first.

Where does AI help hot stamping simulation, and where does it not?

  • Calibration. Fitting h(p, gap) and the kinetics parameters to instrumented-tool and dilatometry data is an inverse problem. Bayesian calibration returns the parameters with their uncertainty, which then carries through to the prediction. See hybrid physics-AI modelling.
  • A surrogate of the martensite and hardness map against transfer time, austenitising, hold time, tool temperature and interface scale, trained on a planned set of coupled runs. See surrogate models for CFD and FEA.
  • Channel-layout optimisation with a validated surrogate in the loop, so hundreds of layouts are screened and the finalists are confirmed in the full model.
  • An anomaly gate that compares line data, such as die temperature, hold time and coolant outlet temperature, with the prediction, and flags a drifting tool before parts go soft. See AI surrogate models you can sign off.

Where AI does not help: replacing the transformation physics when the data covers one steel and one tool. A surrogate cannot invent the CCT diagram of a coil it has not seen, and it extrapolates badly to a new coating or a new tool. Physics supplies the structure, data supplies the parameters, and the surrogate is stated with the range it was trained on.

From one hot stamped part to a press-hardening capability you keep

One part studied once answers one question. The larger return comes when the set-up becomes a workflow your own team runs, so the next part, coil or tool revision costs days of calculation, not weeks of try-out, and every answer is checked the same way.

From one simulation study to a permanent workflow A loop of six steps: physics model, validation, automated runs, AI surrogate, decision tool and field data, which feeds back into the physics model. YOUR WORKFLOW Better with every run physics + data + AI, kept in-house 1Physics model 2Validation 3Automated runs 4AI surrogate 5Decision tool 6Field data CFD · FEA TESTS · TRIALS DOE · HPC CLOUD OPTIMISE DECIDE MONITOR RECALIBRATE
  1. Physics model. CFD and FEA set up once, as templates: geometry clean-up, mesh, solver settings and post-processing that every later run reuses.
  2. Validation. Mesh and time-step studies, and comparison with tests or full-scale data, so the model's uncertainty is known before it is trusted.
  3. Automated runs. Scripted runs over the parameters and operating conditions that matter, on your hardware or in the cloud.
  4. AI surrogate. A machine-learning model trained on those runs answers in seconds and says how sure it is; low-confidence questions go back to the solver.
  5. Decision tool. Optimisation, trade-off charts or operating guidance, in a form your engineers and operators use every day.
  6. Field data. Test, trial and monitoring data flow back to recalibrate the model, so it improves with every project instead of being rebuilt.
Physics, validation, automation and AI in one loop that improves with every run. Illustrative schematic.
  • A validated material and transformation card library. Flow stress, CCT-based kinetics, latent heat and transformation strains for each coil and coating you run, with the dilatometry or trial data behind them and the date they were last checked.
  • A templated coupled model for transfer, forming and quench in your forming solver, with the interface coefficient as a function of pressure and gap, so a new part is set up from the template rather than from scratch.
  • Automated cooling-channel studies over distance, spacing, flow and tool temperature, run over several press cycles, reporting hold time and the worst-cooled region.
  • A surrogate with a stated range for the martensite and hardness map against cycle parameters, with a solver gate that sends out-of-range requests back to the full model.
  • A hand-over to the crash team: the hardness map mapped onto the crash mesh with your property correlation, scripted and versioned.
  • Training and handover, so your engineers can extend the workflow themselves. See our training and software development services.

Questions engineers ask about hot stamping simulation

What is hot stamping, or press hardening?
A blank of boron steel such as 22MnB5 is heated until it is austenitic, transferred to a press, formed and quenched in the closed, cooled die. The fast cooling turns the austenite to martensite, which is where the strength of the part comes from.
Can simulation predict hardness and martensite fraction?
Yes, if the model carries calibrated transformation kinetics for your coil and an interface heat transfer coefficient measured for your tool. Without those two, the plot is a guess in colour.
How do you shorten hot stamping cycle time?
Most of it is the hold until the slowest region has transformed. Better contact, more even channel distance and spacing, and lower tool temperature shorten it, and a conjugate model over several cycles shows which of them pays.
Why is the heat transfer coefficient so uncertain?
It depends on contact pressure, gap, coating or scale and tool surface temperature, and all of them vary over the part and over the cycle. We measure it in instrumented trials and calibrate it as a function of pressure, with an uncertainty.
What do you need from us to start?
The part and tool geometry, the coil data or certificate, any CCT or dilatometry data, the cooling layout and any hardness traverses or cycle logs from trials.

Our sheet metal forming work, the metal forming simulation section and the wider thermal analysis practice cover this ground. If a press-hardened part is on your desk, send us one part number, its coil data and the tool cooling layout, together with the hardness traverse or cycle log from your last trial. We will reply with how we would build the model, which measurements would calibrate it and what the first results would let you decide. Start with a conversation.

OLTRION Engineering

CAE, CFD and simulation team

The OLTRION engineering team runs CFD, FEA and multiphysics simulation programmes for aerospace, automotive, energy and industrial clients, and builds the machine-learning tools that speed them up. These articles come from that project work.

Where this applies

Related capabilities

Have a model that will not converge?

Talk to the engineers who wrote this

Send the geometry, the physics and the decision it has to support. We answer with an approach, not a sales call.

Talk to an engineer