Comportement de déformation et mécanisme d'endommagement du 3C-SiC sous nanoindentation +

PAM-XIAMEN fournit des plaquettes 3C-SiC à des fins de recherche. Pour plus de détails, veuillez contacter :[email protected]

Le carbure de silicium cubique (3C-SiC) présente un grand potentiel dans des domaines de pointe tels que les systèmes microélectromécaniques, l'informatique quantique et les dispositifs optoélectroniques, en raison de sa large bande interdite, de sa conductivité thermique élevée, de son intensité de champ de claquage critique élevée et de son excellente inertie chimique. Cependant, la dureté ultra-élevée et la fragilité intrinsèque du 3C-SiC rendent son usinage ultra-précis extrêmement difficile, et les dommages superficiels/subsurfaces induits par l'usinage sont devenus un goulot d'étranglement majeur limitant ses applications haut de gamme. Par conséquent, une compréhension approfondie du comportement de déformation et des mécanismes d’endommagement du 3C-SiC sous charge de contact revêt une grande importance théorique et technique pour optimiser ses techniques de traitement.

Ces dernières années, les chercheurs ont exploré en profondeur la réponse mécanique du 3C-SiC en combinant des expériences de nanoindentation avec des simulations de dynamique moléculaire. Zhu et coll. a étudié systématiquement le comportement de nanoindentation du 3C-SiC sur les plans cristallins (010), (110) et (111) à l'aide de simulations de dynamique moléculaire (MD), révélant les caractéristiques anisotropes de la déformation élastique et de l'évolution des dislocations, et a introduit le tétraèdre de Thompson pour expliquer le mécanisme intrinsèque de la déformation plastique. Parallèlement, la température, en tant que facteur clé dans les environnements de service réels, affecte de manière significative les propriétés mécaniques du 3C-SiC. Kourbani et coll. découvert grâce à des simulations MD que l'augmentation de la température réduit la dureté et le module de Young du matériau et favorise l'amorphisation. D'un autre côté, les méthodes traditionnelles de mécanique des milieux continus ont des limites dans le traitement de problèmes discontinus tels que l'initiation et la propagation de fissures. Pour résoudre ce problème, Yang et Zhu ont utilisé la péridynamique pour simuler l'évolution des dommages du 3C-SiC sous indentation et rayure, reproduisant ainsi le phénomène d'apparition et la formation de fissures radiales. Ces études fournissent des références multi-échelles pour comprendre le comportement mécanique et la théorie de l'usinage du 3C-SiC.

1. Déformation élastique anisotrope et comportement plastique précoce du 3C-SiC

Le 3C-SiC a une structure de mélange de zinc et sa réponse mécanique dépend fortement de l'orientation des cristaux. Grâce à des simulations MD, Zhu et al. ont constaté que lors de la nanoindentation sur les plans (010), (110) et (111), les régions de déformation élastique présentent des symétries nettement différentes : le plan (010) présente une symétrie quadruple, avec une zone de déformation pyramidale s'étendant le long de la direction <110> ; le plan (110) présente une symétrie centrale, la zone de déformation se scindant en deux parties ; le plan (111) présente une triple symétrie, avec trois zones de distorsion élastique s'étendant également le long de la direction <110>. Ce résultat indique que la direction <110> est la trajectoire optimale pour le transfert de contrainte élastique.

Fig. 1 Morphologie de la zone de distorsion élastique du 3C-SiC à une profondeur d'indentation de 1,1 nm : plan (a, d) (010) ; (b,e) (110) avion ; (c,f) (111) plan ; les flèches rouges indiquent la direction de propagation de la déformation élastique

D’après les courbes charge-profondeur, la longueur de l’étage élastique varie considérablement selon les différents plans cristallins. La courbe Ph pour le plan (010) correspond très bien au modèle de contact Hertz pour h <2,0 nm, alors que les étages élastiques pour les plans (110) et (111) sont significativement plus courts, indiquant une élasticité plastique plus précoce. Notamment, les événements pop-in apparaissent dans les courbes pour les plans (010) et (110), mais ne sont pas observés pour le plan (111). Grâce à l'analyse du déplacement atomique, Zhu et al. a confirmé que l'événement pop-in correspond à un mouvement directionnel d'un grand nombre d'atomes le long de la direction <110>, conduisant à une libération soudaine de l'énergie élastique stockée. Par exemple, lorsque la profondeur d'indentation sur le plan (110) augmente de 2,2 nm à 2,3 nm, le nombre d'atomes dont le déplacement est supérieur à 0,1 nm augmente de 9 558 à 13 767 ; cette migration atomique collective est le mécanisme microscopique direct du pop-in.

Figure 2 (ac) Courbes de profondeur de charge de nanoindentation sur des surfaces 3C-SiC (010), (110) et (111) ; (d) Comparaison de trois courbes de pH et du modèle de contact élastique Hertz

Pour décrire quantitativement l'influence de la répartition des contraintes sur l'initiation plastique, Zhu et al. calculé la distribution des contraintes de cisaillement de von Mises. À une profondeur d'indentation de 2,0 nm, le gradient de contrainte diminue dans la zone de distorsion élastique ; lorsque la profondeur atteint 2,1 nm et que les luxations commencent à se nucléer, la contrainte dans la région de dislocation atteint 8 à 10 GPa, tandis que la contrainte entre les régions déformées et non déformées chute à 5 à 8 GPa. Cela indique que l'expansion des dislocations dissipe efficacement les contraintes locales élevées, permettant à une partie de la structure déformée de revenir à une structure de diamant parfaite.

2. Dislocation Nucléation, propagation et analyse du tétraèdre de Thompson

La déformation plastique du 3C-SiC est principalement dominée par le glissement de dislocation. En utilisant un algorithme d'extraction de dislocations, Zhu et al. identifié que la dislocation initiale est une boucle de dislocation parfaite avec le vecteur de Burgers b = 1/2[1-10], nucléée à partir de la région de contrainte de cisaillement élevée située sous le pénétrateur. Le glissement de dislocation transforme la séquence d'empilement atomique de l'ABCABC d'origine en ABABC, où les atomes de silicium s'écartent de leurs positions d'origine et forment de nouvelles liaisons avec les atomes de carbone adjacents dans la direction <110>. Ce processus marque une transition irréversible de l’élasticité à la plasticité.

Fig. 3 Vues en coupe de la zone de dislocation lors de la nanoindentation : (a) vue locale ; (b) section transversale ; les flèches en pointillés rouges indiquent le glissement des atomes de silicium de leurs positions de liaison d'origine vers de nouvelles positions

À mesure que la profondeur d'indentation augmente, les dislocations embryonnaires évoluent progressivement vers des rubans de failles empilées et des boucles de dislocations prismatiques. En prenant le plan (010) comme exemple : à une profondeur de 2,1 nm, deux défauts d'empilement (SF1 et SF2) situés sur les plans (111) se croisent, formant une luxation en tige d'escalier (SRD) de 1/2 [1-10], qui se connecte ensuite dans une boucle de dislocation en ruban ; lorsque la profondeur augmente jusqu'à 2,6 nm, la dislocation « se détache », formant des boucles de dislocation prismatiques isolées. De même, sur les plans (110) et (111), la séquence d'évolution de la luxation embryonnaire → dislocation du ruban → luxation prismatique est également observée, mais la vitesse d'évolution et la morphologie géométrique diffèrent en raison de l'anisotropie.

Fig. 4 Evolution des dislocations avec profondeur d'indentation sur le plan (010) : (a) h = 2,1 nm ; (b) h = 2,2 nm ; (c) h = 2,5 nm ; (d) h = 2,6 nm

Pour expliquer systématiquement cette anisotropie, Zhu et al. a introduit le modèle du tétraèdre de Thompson. Pour l'indentation sur le plan {100}, les défauts d'empilement sur les quatre plans de glissement {111} (ACE, AEG, ACG, CEG) se croisent, chacun produisant des dislocations en tige d'escalier le long de la direction <110>, la réaction étant représentée par δβ = δE + Eβ, etc. Pour le plan {110}, certains plans de glissement sont perpendiculaires à la surface d'indentation, provoquant une expansion des boucles de dislocation prismatique perpendiculairement à la superficie ; pour le plan {111}, trois dislocations en ruban présentent une triple symétrie, formant finalement trois boucles de dislocation prismatiques. Le tétraèdre de Thompson fournit un cadre géométrique concis mais puissant pour comprendre l'anisotropie plastique des céramiques à structure cubique en diamant.

3. Effet de la température sur les propriétés mécaniques et le comportement de déformation

Dans les applications pratiques, les dispositifs 3C-SiC peuvent fonctionner sur une large plage de températures, de sorte que les effets de la température ne peuvent être ignorés. Kourbani et coll. ont étudié la réponse de nanoindentation du 3C-SiC sur le plan (010) sur une plage de températures de 50K à 1000K à l'aide de simulations MD. Ils ont constaté qu'à mesure que la température augmente, la charge d'indentation maximale diminue de 2 555,4 nN (50 K) à 2 197,4 nN (1 000 K). En utilisant la méthode Oliver – Pharr, le module d'Young diminuait linéairement de 328 GPa (50 K) à 293,5 GPa (1 000 K), avec une valeur à température ambiante (300 K) d'environ 318 GPa, proche de la valeur expérimentale de 314 GPa ; la dureté diminue de la même manière de 45,3 GPa à 37,9 GPa. Cet effet adoucissant résulte de l'augmentation de l'espacement interatomique et de la réduction de l'énergie de la liaison covalente Si – C avec l'augmentation de la température.

La température affecte également de manière significative le comportement des luxations. Kourbani et coll. analysé les morphologies de luxations à différentes températures à l'aide de DXA. Des luxations parfaites avec b = 1/2<110> et des luxations partielles de Shockley avec b = 1/6<112> ont été observées à toutes les températures, les luxations parfaites étant dominantes. À mesure que la température augmente, la longueur totale de la dislocation augmente considérablement et la formation de boucles de dislocation prismatiques de type « lasso » devient plus active. Cela indique que l'activation thermique favorise la nucléation et le glissement des dislocations, améliorant ainsi la déformabilité plastique.

Fig. 5 Analyse DXA et IDS : structures de dislocation du 3C-SiC après indentation jusqu'à la profondeur maximale (h = 2,5 nm) à différentes températures (50K, 100K, 300K, 500K, 800K, 1000K)

Le rôle de la température dans la régulation de la transformation de phase et de l'amorphisation est particulièrement intéressant. En utilisant la méthode d'identification de la structure du diamant, Kourbani et al. ont découvert que des atomes avec une structure hexagonale en diamant (c'est-à-dire une structure wurtzite B4) apparaissent dans la zone d'indentation, comme le montre la figure 6. À mesure que la température augmente, le nombre « d'autres » atomes (représentant l'état amorphe) diminue, tandis que le nombre d'atomes HD (structure wurtzite) diminue également à 1 000 K. L'analyse de la fonction de distribution radiale montre qu'après l'indentation, les pics d'ordre à longue distance disparaissent, confirmant qu'une amorphisation se produit effectivement et devient plus prononcée à des températures plus élevées. L'analyse du nombre de coordination révèle en outre qu'à mesure que la température augmente, la proportion d'atomes avec CN = 4 diminue, tandis que le nombre d'atomes à faible coordination (CN = 1, 2, 3) augmente, indiquant que la température favorise la transition des atomes du réseau parfait à une structure désordonnée. Ces résultats révèlent le double effet de la température sur la stabilité structurelle du 3C-SiC : d’une part, l’activation thermique favorise le glissement des dislocations ; d’autre part, cela exacerbe les dégâts d’amorphisation.

Fig. 6 Répartition des atomes déformés dans le 3C-SiC après indentation à différentes températures (atomes de structure diamant cubique exclus) : HD, HD (1er voisin), HD (2ème voisin) représentent les atomes liés à la structure wurtzite ; « autre » représente des atomes amorphes

4. Simulation des dommages et des fractures basée sur la péridynamique

Bien que les simulations MD puissent résoudre la déformation plastique à l’échelle atomique, leur échelle de calcul est généralement limitée à la plage nanométrique (des dizaines de nanomètres), ce qui rend difficile la comparaison directe avec des expériences à l’échelle micrométrique. De plus, MD est inefficace pour gérer la propagation des fissures lors d’une rupture fragile. Pour résoudre ces problèmes, Yang et Zhu ont simulé les processus d'indentation et de grattage du 3C-SiC en utilisant une approche péridynamique (PD). La péridynamique est une théorie non locale qui remplace les équations aux dérivées partielles par des équations intégrales, ce qui la rend naturellement adaptée au traitement de problèmes discontinus tels que l'initiation et la propagation de fissures.

Dans les simulations d'indentation, Yang et Zhu ont observé un phénomène d'apparition clair dans la courbe charge-profondeur, ce qui est très cohérent avec le comportement d'éclatement plastique observé dans les simulations MD. Les courbes de chargement à différentes vitesses d'indentation (1 à 20 m/s) se chevauchent largement, mais lors du déchargement, des vitesses plus élevées entraînent des chutes de charge plus rapides. Plus important encore, les simulations PD ont reproduit avec succès la morphologie des fissures induites par l'indentation : lorsque la profondeur d'indentation atteint 15 μm à une vitesse de 5 m/s, des fissures radiales, des fissures latérales et des fissures médianes peuvent être clairement observées, ce qui est cohérent avec la morphologie expérimentale des fissures d'indentation. Cela démontre que la PD peut capturer les caractéristiques de rupture des matériaux fragiles à l’échelle du continuum.

Fig. 7 Morphologie des dommages par indentation du 3C-SiC à une vitesse d'indentation de 5 m/s : (a – c) vues de dessus ; (d – i) vues en coupe

Dans des simulations de grattage, Yang et Zhu ont examiné les variations de la force de friction et de la force normale à différentes profondeurs de grattage (1 à 4 μm). À des profondeurs de 1 μm et 2 μm, aucune fissure n’apparaît sur les côtés de la rainure de rayure et l’enlèvement de matière est principalement plastique ; lorsque la profondeur augmente jusqu'à 3 μm et 4 μm, des fissures radiales apparaissent sur les bords des rainures et les forces de friction et normales présentent des fluctuations périodiques. L'énergie de coupe spécifique (énergie consommée par unité de volume de matériau enlevé) augmente de manière non linéaire à mesure que la profondeur de grattage diminue, atteignant jusqu'à 22 GPa à 1 μm, démontrant un effet de taille significatif.

Yang et Zhu ont ensuite simulé le processus de grattage à double pénétrateur pour imiter l'interaction de plusieurs abrasifs lors du meulage réel. Lorsque l'espacement entre les deux pénétrateurs est de 50 μm, la largeur de la rayure est d'environ 92 μm, soit plus de deux fois la largeur d'une rayure à un seul pénétrateur, et le volume d'enlèvement de matière est le plus grand ; à un espacement de 60 µm, les deux rayures sont indépendantes et le volume d'extraction n'est que le double de celui d'un seul pénétrateur ; à un espacement de 40 µm, les zones endommagées se chevauchent, mais le volume total éliminé diminue. Ce phénomène est attribué à l'effet de superposition du champ de déplacement : à un espacement approprié, des fissures se connectent entre les deux pénétrateurs, conduisant à un enlèvement de matière coopératif ; un espacement trop petit entraîne une suppression mutuelle des champs de déplacement, tandis qu'un espacement trop grand ne produit aucune interaction. Cette découverte fournit une base théorique pour optimiser les paramètres du processus de broyage.

Fig. 8 Morphologie des dommages causés par les rayures à double pénétrateur (profondeur d'indentation 3 μm) : (a, d) espacement de 40 μm, rayures irrégulières ; (b, e) espacement de 50 µm, largeur de rayure considérablement augmentée, volume d'enlèvement de matière maximal ; (c,f) espacement 60μm, deux rayures indépendantes

Fig. 9 Vues latérales en coupe du grattage à double pénétrateur : (a) espacement de 40 μm ; (b) espacement de 50 μm ; (c) espacement de 60 μm

Grâce à une combinaison de simulations MD et PD, les chercheurs ont révélé la déformation anisotrope, l'évolution des dislocations, les effets de la température et les mécanismes de dommages/fractures du 3C-SiC sous nanoindentation. L'étude révèle que les comportements élastiques et plastiques dépendent fortement de l'orientation des cristaux, avec des dislocations se nucléant le long du système de glissement {111}<110> et évoluant en boucles prismatiques ; l'augmentation de la température réduit la dureté et le module tout en favorisant l'amorphisation. La méthode PD reproduit avec succès le phénomène de pop-in lors de l'indentation et de la propagation des fissures lors des rayures, et identifie un espacement optimal entre les doubles pénétrateurs pour améliorer l'efficacité de l'enlèvement de matière. Les travaux futurs pourraient développer davantage des modèles de couplage multi-échelles et des fonctions potentielles apprises par machine pour réaliser des prédictions quantitatives depuis les mécanismes atomiques jusqu'aux processus d'usinage.

Que vous ayez besoin d'une plaquette 3C-SiC pour la recherche ou pour des applications industrielles, veuillez nous contacter par e-mail à[email protected]et[email protected].

 

Les références:

1. Zhu, B., Zhao, D., Zhang, Z., Niu, Y., Zhang, Z., Zhao, J.,… et Zhao, H. (2024). Etude atomique du comportement en déformation et de l'effet d'anisotropie du 3C-SiC sous nanoindentation. Journal de recherche et de technologie des matériaux, 28, 2636-2647.
2. Kourbani, H., Bogtob, S., Samiri, A., Hassani, A. et Hasnaoui, A. (2025). Nanoindentation thermique du carbure de silicium (3C) 3C-SiC par simulation de dynamique moléculaire. Communications à semi-conducteurs, 403, 116013.
3. Xu, Y. et Zhu, P. (2022). Simulations péridynamiques des dommages en indentation et grattage du 3C-SiC. Journal de recherche sur les matériaux, 37(23), 4381-4391.