Étude de l'interaction dislocation - amas de lacunes par simulations numériques

Vacancy clusters have been observed and characterized experimentally in highly pure metals after plastic deformation or after a particular sequence of heat treatments. These clusters hinder the dislocation propagation and can therefore harden the metal.Using numerical simulations we have explored different mecanisms of dislocation propagation through a vacancy-cluster distribution, for several applied shear stress and temperature. At high stresses, the force applied on the dislocation becomes greater than the pinning forces acting on the line. The dislocation gets through the cluster distribution by gliding and shearing the clusters. The dependence of the pinning force with the cluster size is adjusted on our molecular static simulations. In this stress range, the pinning configurations are rare and the thermal activation is sufficient to unpin the line. The probability for the line to pass the pinning configuration depends on the activation enthalpy, a parameter that we have also estimated using an analytical model adjusted on our atomistic results. At lower stresses, when the applied force is below the pinning forces induced by the cluster, the probability that the dislocation unpins by pure glide becomes negligeable. The diffusion of vacancies, emitted preferentially from the vacancy clusters, intervenes and promotes the formation of jogs that contributes to the unpinning of the line. Such a mecanism is the glide assisted by climb. The emission, the absorption and the vacancy migration barriers have been determined by molecular static and are highly dependent on the elastic field and the atomic network distortion induced by the dislocation. This promotes a strong diffusion anisotropy in the vicinity of the dislocations which leads in particular to the pipe diffusion mechanism. The evolution with time of all these mechanisms has been studied using an elastic line model coupled to a kinetic Monte Carlo algorithm in which the parameters come from our atomistic simulations. According to the model assumptions, we obtained an estimation of dislocation velocity as a function of the applied shear stress and the temperature. We used the Orowan's law to estimate the strain rate related to such mechanisms.Des amas de lacunes ont été observés et caractérisés expérimentalement dans les métaux de haute pureté après déformation plastique ou après une suite de traitements thermiques particuliers. Ces amas sont des obstacles à la propagation des dislocations et peuvent par conséquent induire un durcissement du métal.Cette étude par simulations numériques a permis d'explorer différents mécanismes de propagation de dislocations dans une concentration d'amas en fonction de la contrainte de cisaillement appliquée et de la température. À haute contrainte, la force appliquée sur la dislocation devient supérieure aux forces d'ancrage s’exerçant sur la ligne. La dislocation franchit la distribution d'amas en glissant et en cisaillant les amas. La dépendance de la force d'ancrage en fonction de la taille de l'amas est ajustée sur nos simulations de statique moléculaire. Dans ce domaine de contrainte, les configurations d'amas ancrant la dislocation sont rares et l'activation thermique suffit à désancrer la ligne. La probabilité de désancrer la ligne dépend de l'enthalpie d'activation, un paramètre que nous avons également estimé à l'aide d'un modèle analytique ajusté sur nos résultats atomistiques. À plus faible contrainte, lorsque la force appliquée est inférieure aux forces d'ancrage induites par les amas, la probabilité que la dislocation se désancre uniquement par glissement devient faible. La diffusion des lacunes, émises préférentiellement des amas, intervient alors et favorise la formation de crans. Cela contribue au désancrage de la ligne. Ce mécanisme est le glissement assisté par la montée. Les barrières d'émission, d'absorption et de migration de lacunes ont été déterminées par statique moléculaire et sont fortement dépendantes du champ élastique et de la distorsion du réseau atomique générés par la présence de la dislocation. Cela induit une forte anisotropie de diffusion au voisinage des dislocations qui conduit notamment au mécanisme de 'pipe diffusion'. L'évolution au cours du temps de l'ensemble de ces mécanismes a été étudiée à l'aide d'un modèle de ligne élastique couplé à un algorithme de Monte Carlo cinétique dont l'ensemble des barrières d'énergie provient de nos simulations atomistiques. Moyennant les hypothèses du modèle, nous avons alors obtenu une estimation de la vitesse des dislocations en fonction de la contrainte et de la température appliquée. Nous avons ensuite utilisé la loi d'Orowan pour estimer la vitesse de déformation liée à ces mécanismes