doc - crc · doc - crc numerical modeling of zirconium alloys hot extrusion - alexis gaillac - p.7...
TRANSCRIPT
DOC - CRC
Presentation title – Presenter/ref. - 17 June 2016 - p.1Presentation title – Presenter/ref. - 17 June 2016 - p.1
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion
Alexis Gaillac, Pierre Barberis
18th ASTM Int'l Symp on Zr in the Nuclear Industry
16-19 May 2016 – Hilton Head Island - USA
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.3
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Outline
Introduction
The Hot extrusion process
Numerical modeling of Zr manufacturing processes
Description and validation of the models
Analytical and Finite Element model
Recrystallization model
Validation of the models
Applications
Extrusion loads
Recrystallization
Lubricants and tools optimization
Conclusion
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.4
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Outline
Introduction
The Hot extrusion process
Numerical modeling of Zr manufacturing processes
Description and validation of the models
Analytical and Finite Element model
Recrystallization model
Validation of the models
Applications
Extrusion loads
Recrystallization
Lubricants and tools optimization
Conclusion
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.5
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Tubes hollows
With suitable geometry
for cold pilgering
Without surface defects
Without microstructural
heterogeneities
IntroductionThe Hot Extrusion process
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.6
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Direct hot extrusion through a conical die
Billets machined from forged bars
Extrusion temperature = 600°C
Extrusion ratio = 10
IntroductionThe Hot Extrusion process
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.7
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
IntroductionNumerical modeling of Zr manufacturing processes
Since the 1990’s, AREVA NP uses numerical modeling for
a mastering of the whole Zr manufacturing processes.
Extractive metallurgy
VAR Melting
Forming processes
(Forging, Extrusion, Rolling, … )
Heat treatments, Quenching
US Inspections
VAR Melting
ASTM STP 1529, Chengdu 2010
Cold Pilgering
ASTM STP 1354,
Toronto 1998
J. Mater. Process.
Technol., Vol. 117, 2006
Quenching
ASTM STP 1543,
Hyderabad 2013
Resistance Butt
Welding
ASTM STP 1543, Hyderabad 2013
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.8
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
IntroductionNumerical modeling of Zr Hot Extrusion
For Hot Extrusion, the purpose of numerical modeling is to :
Better understand the effect of extrusion parameters on the
Extrusion load
Microstructural evolutions during extrusion
Stresses and temperatures applied to the tools
Check the capability of extrusion press and tools for a new configuration
Optimize extrusion parameters for a given material and geometry, in order to
ensure a good quality of the extruded products
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.9
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Outline
Introduction
The Hot extrusion process
Numerical modeling of Zr manufacturing processes
Description and validation of the models
Analytical and Finite Element model
Recrystallization model
Validation of the models
Applications
Extrusion loads
Recrystallization
Lubricants and tools optimization
Conclusion
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.10
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Analytical Model (AM)
First simple estimation for extrusion load
Isothermal and homogeneous deformation of the billet is assumed
𝐹 = 𝐹𝑑 + 𝐹𝑓
𝐹𝑑 = 𝜋. 𝑅𝑜𝑏2 − 𝑅𝑖𝑏
2 . 𝜎0. ln 𝜆
𝐹𝑓 = 2. 𝜋. 𝑅𝑜𝑏 + 𝑅𝑖𝑏 . 𝐿𝑏. 𝜎0.ഥ𝑚
3
𝜆 =𝑅𝑜𝑏2 − 𝑅𝑖𝑏
2
𝑅𝑜𝑡2 − 𝑅𝑖𝑡
2
Extrusion load
Deformation load
Friction load
Reduction ratio
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.11
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Analytical Model (AM)
First simple estimation for extrusion load
Isothermal and homogeneous deformation of the billet is assumed
Strain rate
(s-1)
Flow stress (MPa) at Temperature (°C)
500 600 700 800
0.1 224 165 104 74
1 240 180 132 96
10 257 196 169 127
𝑇𝑒𝑥𝑡𝑟𝑢𝑠𝑖𝑜𝑛 =𝑇𝑖𝑛𝑖𝑡 + 𝑇𝑚𝑎𝑥
2
𝑇𝑚𝑎𝑥 = 𝑇𝑖𝑛𝑖𝑡 +ln 𝜆 . 𝜎0𝜌. 𝐶𝑝
𝐹𝑑 = 𝜋. 𝑅𝑜𝑏2 − 𝑅𝑖𝑏
2 . 𝜎0. ln 𝜆
𝐹𝑓 = 2. 𝜋. 𝑅𝑜𝑏 + 𝑅𝑖𝑏 . 𝐿𝑏. 𝜎0.ഥ𝑚
3
Average flow
stress from
torsion and
compression
tests (Zircaloy-4)
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.12
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Finite Element Model (FEM)
2D thermo-mechanical coupled model, Forge NxT software
Initial geometries and meshing
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.13
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Finite Element Model (FEM)
2D thermo-mechanical coupled model, Forge NxT software
Geometries and meshing during extrusion
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.14
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Finite Element Model (FEM)
2D thermo-mechanical coupled model, Forge NxT software
Flow stress as a function of temperature, strain and strain rate
From torsion and compression tests (Zircaloy-4)
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.15
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Recrystallization Model
Model developed by Gaudout [1] and based on the
Gourdet-Montheillet [2] and Avrami [3] approaches
From FEM results (temperature, strain and strain rate)
Dynamic microstructure fragmentation during extrusion (Gourdet-Montheillet)
Post-dynamic grain growth during cooling of the extruded tube (Avrami)
The result is a recrystallized volume fraction ( = surface fraction of equiaxed
grains on the micrographs)
Model parameters were fitted on recrystallized volume fraction after torsion
tests and heat treatments
Schematic representation of the
microstructure during fragmentation [1].
Thick lines are grains boundaries
(crystallographic misorientation > 15°),
fine lines are sub-grain boundaries.
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.16
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Validation of the models
Extrusion loads
Calculated loads are consistent with experimental one
Extrusion load (T)
Experimental Finite Element Model Analytical Model
Min 1600
Max 18001700 1667
Experimental loads
for Zircaloy-4 tubes
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.17
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Validation of the models
Temperature
Calculated temperatures are consistent with experimental ones
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.18
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Validation of the models
Recrystallization
Recrystallized volume fraction ( = surface fraction of equiaxed grains on
the micrographs) after small scaled extrusion tests and heat treatments
Rather good agreement between numerical results and experimental ones
© AREVA NP – Property of AREVA NP
Figures from Gaudout [1]
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.19
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Outline
Introduction
The Hot extrusion process
Numerical modeling of Zr manufacturing processes
Description and validation of the models
Analytical and Finite Element model
Recrystallization model
Validation of the models
Applications
Extrusion loads
Recrystallization
Lubricants and tools optimization
Conclusion
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.20
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Extrusion loads
Extrusion load curve can be decomposed into 4
sequences and explained using FEM results
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.21
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Extrusion loads
Extrusion load curve can be decomposed into 4
sequences and explained using FEM results
Filling of the gaps between billet and tools
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.22
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Extrusion loads
Extrusion load curve can be decomposed into 4
sequences and explained using FEM results
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.23
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Extrusion loads
Extrusion load curve can be decomposed into 4
sequences and explained using FEM results
Peak load (lower temperature and maximum flow stress because of
lower strain), front end of the tube exits the die
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.24
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Extrusion loads
Extrusion load curve can be decomposed into 4
sequences and explained using FEM results
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.25
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Extrusion loads
Extrusion load curve can be decomposed into 4
sequences and explained using FEM results
Maximum temperature and minimal load, then pseudo-steady-state
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.26
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Extrusion loads
Extrusion load curve can be decomposed into 4
sequences and explained using FEM results
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.27
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Extrusion loads
Extrusion load curve can be decomposed into 4
sequences and explained using FEM results
Rising load due to the chilling effect of the tools
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.28
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Recrystallization at front end
RX ratio = surface fraction of equiaxed grains
Outer surface: RX ratio = 95% Mid-thickness: RX ratio = 70%
300mm from front end = steady-state microstructure
© AREVA NP – Property of AREVA NP
R
L
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.29
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Recrystallization at front end
Results from the FEM + recrystallization model
Consistent with experimental ones
RX ratio lower at front end because of lower strain
Following our experience, to avoid cracks or voids during pilgering,
necessary discarded length (RX<50%) is about 100 mm
This model is validated and can be applied to new extrusion configurations
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.30
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Lubricants and tools optimization
Choice of the right tools materials and lubricant
Mandrel: up to 700 MPa at 600°C
Die: up to 680°C at the end of extrusion
Lubricant: temperature range 450 – 650°C
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.31
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Outline
Introduction
The Hot extrusion process
Numerical modeling of Zr manufacturing processes
Description and validation of the models
Analytical and Finite Element model
Recrystallization model
Validation of the models
Applications
Extrusion loads
Recrystallization
Lubricants and tools optimization
Conclusion
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.32
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Conclusion
Conclusion
Models are validated regarding to experimental data
For the standard extrusion configuration, models give better
understanding of the effect of process parameters on the
Extrusion load
Microstructural evolutions during extrusion
Stresses and temperatures applied to the tools
For new configurations (reduction ratio, extruded material, …),
models are great tools to
Check the capability of extrusion press and tools
Optimize extrusion parameters for a given material and geometry, in
order to ensure a good quality of the extruded products
© AREVA NP – Property of AREVA NP
DOC - CRC
Numerical Modeling of Zirconium Alloys Hot Extrusion - Alexis Gaillac - p.33
18th Int'l Symp on Zr in the Nuclear Industry - 16-19 May 2016 – Hilton Head Island - USA
Thanks
Any questions ?
© AREVA NP – Property of AREVA NP