核电包壳材料氢脆问题分析最完整讲义

0 vol% hydride,0.075 with stress oriented at 0
Numerical studies of ζ hydride in zirconium: Selection & re-orientation of the trigonal precipitates in a hexagonal matrix under external stress
Lingfei ZHANG, Ludovic THUINET, Alexandre LEGRIS
0 Configurational B(n ) s ik nk ij (n )s 0 jl nl energy
•
Define:
Therefore, minimal elastic strain energy reduces to find minimal configurational energy, or B(n); n vector is normal of precipitate’s habit plane as shown in Figure to the left. NB: Khachaturyan theory has been only able to determine the habit plane on which the n vector is normal to it that hydride will precipitates within it; none of the hydride true morphology or the actual orientation is possible through this approach. E.g., two habit planes could coincide with the same orientation as projected onto the basal plane!
Reference: Mishima, Y., S. Ishino, and H. Kawnishi, Some observations on the dissolution and precipitation of zirconium hydrides in a-zirconium by electron microscopy. Journal Name: pp 489-496 of International Congress on Hydrogen in Metals. Vols. 1 and 2. Paris Editions Science et Industrie, 1972.
1 1 1 d 3k d 3k 0 * 0 0 E Vincl Cijkl ij ij kl kl (V Vincl )Cijkl ij kl s n G [n (q )]s 0 jl nl 3 3 ik k ij 2 2 2 (2 ) (2 )
2D modeling of microstructure dependent strain energy vs. phase field modeling
Strain energy for Hydride volume fraction approach zero Phase field modeling for single hydride Phase field modeling for a population of hydrides Strain energy for Hydride volume fraction >>0
28th of September 2011, Les Renardières EDF R&D
Experimental review
Re-orientation: time-dependent, concentration dependent
Hoop stress definition sh=pr/t, t:thickness, p: internal pressure, r:radius
P M ijkl ,ijkl Elastic tensor of precipitates and matrix respectively

is volume fraction
1 P Gil ijkl n j nk
Therefore the Green function is modified to
Theoretical background
• Khachaturyan developed theory on elastic strain energy due to inhomogeneous solids which gives approach to calculate the energy due to precipitates configuration; The equilibrium state of inhomogeneous solids with coherent lattice should correspond to minimal elastic strain energy
3D modeling of minimum strain energy surface
Conclusion and future works Acknowledgement
Literatures and theoretical background
literature milestones Microelasticity of Khachaturyan’s inhomogeneous solid Limits of microelastic theory
* * Therefore, B( n ) s ik nk Gij ( n )s jl nl
Is modified with the elastic tensor of precipitates when w approach zero this is compared to the B function using elastic tensor of average
* kl
* * Therefore, B( n ) s ik nk Gij ( n )s jl nl
n is unit vector of normal to habit plane and varies from 0-p to calculate B
Limit of microelastic theory
Strain energy under uniform stress
When the system of multicomponent with coherent interaction, the external stress could affect the configurational energy, hence the arrangement of precipitates due to minimal B(n)
Literatures milestones
f: the angle between unit vector and referential axis q: the angle between applied stress and referential axis
The other documents that well studied the hydrides distribution free or under stress is found only in 1970’s by Mishima et al in which the massive hydrides appears to be less affected by the stress and the dissolution instead of re-orientation was found on the basal plane. This means reorientation could occur in prismatic planes that has to be studied theoretically.
* 0 * s ij Sijkl kl
0 Sijkl
is average compliance tensor dependent on elastic stiffness and volume fraction of each materials component is effective strain due to applied uniform stress for instance along axis
One of the key tasks of the B function calculation lying with the determination of the Green function in 2D or 3D. Khachaturyan et al has developed the theoretical approach that inverse Green function is related to the elastic tensor of matrix and precipitates in coherent case by spacial unit vector xi x j therefore the inverse G-1
Higher the hydrides concentration higher stress level is needed to trigger re-orientation
Outline
Literature and theoretical background 2D modeling on basal plane: minimium strain energy and phase field method 2D modeling on prismatic plane: minimum strain energy and phase field method
UniversiBiblioteka Baidué Scientifique et Technologique de Lille
Laboratoire de Métallurgie Physique et Génie des Matériaux USTL Bât C6 – 2e étage 59655 Villeneuve d’Ascq
0 B( n ) s ik nk Gij ( n )s 0 jl nl
The stress term in above equation is the stress free stress; In order to calculate the mimimal B(n) as function of angel respective to the reference axis, one should replace stress-free stress by effective stress.
1 Gil ijkl n j nk
The use of above equation should be proceeded with caution since no previous results have been found as benchmark. Despite the definition of ijkl is also called elastic average by arbitrary P M ijkl ijkl (1 )ijkl