Spotlight #67: Inżynieria Materiałowa — Dyslokacje, Diagramy Fazowe i Wzrost Ziaren

I. Poślizg dyslokacji krawędziowej — dlaczego metale są słabe (i jak to naprawiamy)

Oto zagadka, która przez dekady zastanawiała fizyków: teoretyczna wytrzymałość na ścinanie idealnego kryształu metalu — naprężenie potrzebne, by przesunąć jedną połowę kryształu względem drugiej jednocześnie — wynosi około G / 6 do G / 30, gdzie G to moduł ścinania. Dla miedzi daje to w przybliżeniu 1–10 GPa. A jednak prawdziwa miedź ulega odkształceniu plastycznemu już przy 50–100 MPa. To dwa do trzech rzędów wielkości słabiej, niż przewiduje teoria. Powodem są dyslokacje.

Dyslokacja krawędziowa to defekt liniowy w krysztale: dodatkowa półpłaszczyzna atomów, która kończy się wewnątrz sieci. Zamiast przesuwać całą płaszczyznę atomów jednocześnie — proces wymagający GPa naprężenia — dyslokacja porusza się, zrywając i odtwarzając jedno wiązanie na raz wzdłuż swojej linii. Dyslokacja przesuwa się o jedną stałą sieci (jeden wektor Burgersa) na każde zerwane wiązanie, przemiatając płaszczyznę poślizgu jak fałda na dywanie. Naprężenie potrzebne do przesunięcia dyslokacji to naprężenie Peierlsa–Nabarro:

tau_PN = (2G / (1 - nu)) * exp(-2*pi*d / (b * (1 - nu)))

  G   : modul scinania (GPa)
  nu  : wspolczynnik Poissona (~0.33 dla wiekszosci metali)
  d   : odstep miedzy sasiednimi plaszczyznami poslizgu (m)
  b   : wielkosc wektora Burgersa (stala sieci, ~0.25 nm dla Cu)

Dla Cu (FCC): tau_PN ~ 10^-4 * G  ~ 4 MPa  (poslizg na plaszczyznach {111})
Dla Si (DC):  tau_PN ~ 10^-2 * G  ~ 2 GPa  (tlumaczy kruchosc krzemu)

Zależność wykładnicza od d/b jest kluczowym spostrzeżeniem. Metale FCC, takie jak miedź, aluminium i złoto, mają gęsto upakowane płaszczyzny {111} o dużym odstępie i małych wektorach Burgersa, co daje bardzo niskie naprężenie Peierlsa i doskonałą ciągliwość. Diamentowy krzem ma mały odstęp płaszczyzn i duże wektory Burgersa, więc dyslokacje są w temperaturze pokojowej niemal nieruchome — dlatego krzem się rozpada, zamiast się zginać.

Mechanizmy umacniania

Jeśli to dyslokacje sprawiają, że metale są słabe, to blokowanie ruchu dyslokacji jest sposobem na ich ponowne wzmocnienie. Każdy mechanizm umacniania w metalurgii strukturalnej działa poprzez utrudnianie poślizgu dyslokacji:

• Umacnianie przez zgniot — dyslokacje spiętrzają się i splątują, blokując się nawzajem. Walcowanie na zimno podwaja lub potraja granicę plastyczności metalu.
• Umacnianie roztworowe — atomy domieszki lokalnie zniekształcają sieć, tworząc pola naprężeń, które przypinają dyslokacje. Dlatego stal nierdzewna (Fe + Cr + Ni) jest znacznie mocniejsza niż czyste żelazo.
• Umacnianie wydzieleniowe — drobne cząstki drugiej fazy (np. Al&sub3;Cu w duraluminium) działają jak przeszkody, które dyslokacje muszą przeciąć lub obejść (mechanizm Orowana).
• Umacnianie granicami ziaren — granice ziaren całkowicie blokują ruch dyslokacji, co opisuje zależność Halla–Petcha: sigma_y = sigma_0 + k / sqrt(d), gdzie d to średnica ziarna.

II. Diagramy fazowe stopów dwuskładnikowych — odczytywanie przepisu na mikrostrukturę

Diagram fazowy to mapa. Dla stopu dwuskładnikowego (dwa składniki, powiedzmy cyna i ołów) mówi, które fazy — ciecz, roztwór stały alfa, roztwór stały beta, faza międzymetaliczna — są stabilne przy każdej kombinacji składu i temperatury. To pojedynczy najważniejszy diagram w metalurgii: każdy proces odlewania, obróbki cieplnej i specyfikacja lutowania na nim polega.

Najbardziej charakterystyczną cechą jest punkt eutektyczny: jedyny skład, przy którym stop ma najniższą możliwą temperaturę topnienia i gdzie ciecz zamarza bezpośrednio w dwa ściśle wymieszane ciała stałe, bez przechodzenia przez dwufazową strefę papkowatą. W układzie Sn–Pb eutektyk znajduje się przy 61,9% wag. Sn i 183 °C — podstawa tradycyjnego lutu elektronicznego 60/40, wybranego właśnie dlatego, że krzepnie ostro (bez zakresu papkowatego) i w najniższej możliwej temperaturze.

Reakcja eutektyczna:
  L  -->  alfa + beta     (przy T_eutektycznej, jeden sklad)

Dla Sn-Pb:
  T_eutektyczna = 183 °C,  C_eutektyczny = 61.9 % wag. Sn
  faza alfa: roztwor staly bogaty w Pb (maks. 19.2% wag. Sn przy T_eutektycznej)
  faza beta: roztwor staly bogaty w Sn (maks. 2.5% wag. Pb przy T_eutektycznej)

Regula dzwigni (w obszarze dwufazowym alfa + L):
  Przy calkowitym skladzie C_0, temperaturze T:
    x_alfa = (C_L - C_0) / (C_L - C_alfa)    [udzial fazy stalej]
    x_L    = (C_0 - C_alfa) / (C_L - C_alfa) [udzial fazy ciekej]

  C_alfa : sklad fazy alfa (odczytany z lewej granicy)
  C_L    : sklad fazy ciekej (odczytany z prawej granicy)
  "Dzwignia": C_0 to punkt podparcia; dlugosci daja udzialy faz.

Reguła dźwigni nazywa się tak, ponieważ dwufazowa linia wiążąca działa jak dźwignia: udział każdej fazy jest odwrotnie proporcjonalny do jej odległości od całkowitego składu, tak jak cięższa masa musi znajdować się bliżej punktu podparcia, by zrównoważyć lżejszą. Jest to bezpośrednia konsekwencja zachowania masy i dotyczy każdego obszaru dwufazowego na dowolnym równowagowym diagramie fazowym.

Odczytywanie ścieżki chłodzenia

Rozważmy podeutektyczny stop Sn–Pb o składzie 40% wag. Sn, powoli chłodzony od 300 °C. Ścieżka chłodzenia przecina likwidus przy około 240 °C, gdzie zarodkuje pierwsza faza stała (bogata w Pb faza alfa). Gdy temperatura spada przez dwufazowy obszar L + alfa, udział fazy stałej rośnie, a obie fazy zmieniają skład, podążając za granicami likwidusu i solidusu. Przy 183 °C pozostała ciecz (teraz dokładnie o składzie 61,9% wag. Sn) przechodzi przemianę eutektyczną: zamarza w drobną warstwową mieszaninę alfa i beta. Ostateczna mikrostruktura to pierwotne dendryty alfa otoczone matrycą eutektyczną — dokładnie to, co widać w przekroju starego lutu hydraulicznego.

Symulacja pozwala umieścić kursor w dowolnym punkcie diagramu fazowego i zobaczyć na żywo obliczoną regułę dźwigni, z kolorowo oznaczonymi udziałami i składami faz aktualizowanymi w trakcie przeciągania. Animowany tryb chłodzenia pokazuje zdarzenia zarodkowania, ewolucję reguły dźwigni i ostateczną zamrożoną mikrostrukturę w panelu obok siebie.

III. Wzrost ziaren — model Pottsa i dlaczego wyżarzane metale mięknieją

Gdy zgniatany na zimno metal jest podgrzewany — wyżarzany — zachodzą kolejno trzy procesy: zdrowienie (defekty punktowe anihilują), rekrystalizacja (nowe, wolne od naprężeń ziarna zarodkują i rosną) i wreszcie normalny wzrost ziaren (duże ziarna rosną kosztem małych, napędzane wyłącznie redukcją powierzchni granic ziaren). To właśnie ten ostatni proces modeluje nasza symulacja.

Siłą napędową jest termodynamika. Granice ziaren to obszary o wysokiej energii — atomy na granicy znajdują się w mniej korzystnym środowisku niż atomy wewnątrz ziarna, z nadmiarem energii rzędu 0,5–1 J m⁻² dla granic o dużym kącie w metalach. Całkowita energia granic ziaren polikryształu jest proporcjonalna do całkowitej powierzchni granic, która skaluje się jak 1/d (d = średnia średnica ziarna). Rozrost ziaren zmniejsza tę energię, eliminując granice, powiększając średni rozmiar ziarna w czasie.

Model Pottsa reprezentuje każde miejsce sieci jako zmienną spinową Q ∈ {1, 2, …, Q_max}, gdzie każda wartość reprezentuje odrębną orientację ziarna. Miejsca o tym samym Q należą do tego samego ziarna; granice istnieją tam, gdzie sąsiadujące miejsca się różnią. W każdym kroku Monte Carlo losowo wybierane jest miejsce sieci, przypisywana jest mu próbna orientacja od sąsiadów, a próba jest akceptowana lub odrzucana zgodnie z kryterium Metropolisa:

Hamiltonian:
  H = J * sum_{<i,j>} (1 - delta(Q_i, Q_j))
  J  : energia granicy na pare miejsc (J > 0 karze roznych sasiadow)
  delta : delta Kroneckera (1 jesli to samo ziarno, 0 jesli rozne)

Regula akceptacji Metropolisa:
  delta_H = H_nowe - H_stare
  P_akceptacji = 1            jesli delta_H <= 0  (ruch korzystny)
  P_akceptacji = exp(-delta_H / kT)  w przeciwnym razie

Prawo wzrostu sredniej powierzchni ziarna (prawo paraboliczne):
  <A>(t) = <A>(0) + 2*M*gamma*t
  lub rownowaznie:  d(t)^2 - d(0)^2 = K*t
  K : stala szybkosci (proporcjonalna do ruchliwosci granicy M i energii gamma)
  d(t) : srednia srednica ziarna w czasie t
  wykladnik n = 0.5 w idealnym (normalnym) wzroscie ziaren

Paraboliczne prawo wzrostu d ∝ t^0,5 jest wynikiem pola średniego: zakłada, że wszystkie ziarna są identycznymi kulami (lub okręgami w 2D), a wszystkie granice poruszają się z tą samą prędkością. Prawdziwe metale często pokazują d ∝ t^n z n < 0,5, ponieważ ciągnięcie roztworowe, przypinanie wydzieleniami (przypinanie Zenera) lub efekty teksturowe selektywnie spowalniają ruch granic. Symulacja mierzy n bezpośrednio z nachylenia wykresu log–log średniej powierzchni względem czasu Monte Carlo.

Przypinanie Zenera

Dodanie do symulacji dyspersji drobnych, nieruchomych cząstek (obojętne miejsca, które nigdy nie zmieniają orientacji) zatrzymuje wzrost ziaren: granice wyginają się wokół cząstek i w końcu zostają przypięte. Limit Zenera dla maksymalnej średnicy ziarna wynosi d_Z = (4r) / (3f), gdzie r to promień cząstki, a f to udział objętościowy. Dlatego drobne cząstki tlenkowe są celowo wprowadzane do stali umacnianych dyspersyjnie tlenkami (ODS), stosowanych w reaktorach jądrowych: zamrażają one strukturę ziaren w temperaturze eksploatacji, zapobiegając utracie wytrzymałości, która w przeciwnym razie wynikałaby z rozrostu ziaren.

Wypróbuj sam

Wszystkie trzy symulacje działają w przeglądarce bez instalacji:

Myśl na zakończenie

Inżynieria materiałowa jest, w swej istocie, dyscypliną łączącą zachowanie w skali atomowej z wydajnością inżynierską. Trzy tematy tego spotlightu obejmują siedem rzędów wielkości w skali długości — od wektora Burgersa dyslokacji miedzi o wielkości 0,25 nm, przez mikrometrowe lamele eutektycznego lutu, po milimetrowe ziarna wyżarzonej blachy stalowej — a jednak wszystkie trzy są rządzone tą samą podstawową zasadą: systemy minimalizują energię swobodną, a tempo, w jakim to robią, determinuje wszystko, od granicy plastyczności po trwałość eksploatacyjną. Symulacje w tym artykule czynią tę podróż przez skale namacalną, pozwalając w kilka minut zbadać to, czego zrozumienie zajęło metalurgom dekady.