Spotlight #26 — Inżynieria materiałowa

Inżynieria materiałowa leży na przecięciu fizyki, chemii i inżynierii. Każdą właściwość materiału — twardość, przewodnictwo, ciągliwość, temperaturę topnienia — można wyjaśnić strukturą atomową i defektami. Ten spotlight przeprowadza przez osiem symulacji, które czynią te powiązania namacalnymi: możesz obracać sieci krystaliczne, obserwować dyfuzję atomów wzdłuż granic ziaren, zginać wirtualny pręt stalowy poza granicę plastyczności i odkrywać, dlaczego krzem ma przerwę energetyczną, a miedź nie.

1. Struktury krystaliczne — BCC, FCC, HCP i diament

Każdy metal w stanie stałym organizuje swoje atomy w jedną z czternastu sieci Bravais'go. Trzy najczęstsze to sieć regularna przestrzennie centrowana (BCC: żelazo w temperaturze pokojowej, wolfram, chrom), sieć regularna ściennie centrowana (FCC: aluminium, miedź, złoto, nikiel) i sieć heksagonalna gęsto upakowana (HCP: tytan, magnez, cynk). Struktura diamentu (krzem, german, diament) leży u podstaw całej nowoczesnej elektroniki.

Struktura   Atomy/komórkę   Liczba koord.   APF     Przykład
BCC         2               8               0,680   Fe-α, W, Cr
FCC         4               12              0,740   Al, Cu, Au
HCP         6 (efekt. 2)    12              0,740   Ti, Mg, Zn
Diament     8               4               0,340   Si, Ge, C(diament)

APF = współczynnik upakowania = objętość atomów / objętość komórki elementarnej
FCC i HCP są oba gęsto upakowane (APF 0,740, różnią się tylko ułożeniem ABCABC vs ABABAB)

Parametr sieci a (FCC z promienia atomowego r): a = 2√2 · r
Odległość międzypłaszczyznowa (wskaźniki Millera hkl, sieć regularna): d_hkl = a / √(h²+k²+l²)
          

Przeglądarka 3D struktur krystalicznych pozwala obracać wszystkie sześć popularnych struktur — w tym NaCl (jonową) i diament — z renderowaniem sortowanym według głębokości (algorytm malarski), więc widoczne są też atomy wewnętrzne. Przełącz superkomórkę 2×2×2, by policzyć odległości sfer koordynacyjnych, i obserwuj, jak panel liczby koordynacyjnej aktualizuje się przy przełączaniu struktur.

2. Dyfuzja atomowa i prawa Ficka

Atomy w ciele stałym nie są nieruchome — przeskakują między miejscami sieci przez wakanse, pozycje międzywęzłowe lub granice ziaren. Opis makroskopowy to drugie prawo Ficka: ∂C/∂t = D·∇²C. Sam współczynnik dyfuzji D zależy jednak od temperatury zgodnie z zależnością Arrheniusa — wykładniczą barierą aktywacji wynikającą z energii potrzebnej, by przecisnąć się obok sąsiednich atomów.

Pierwsze prawo Ficka (strumień w stanie ustalonym):
  J = −D · ∂C/∂x        [atomy m⁻² s⁻¹]

Drugie prawo Ficka (zależne od czasu):
  ∂C/∂t = D · ∂²C/∂x²   (1D, stałe D)

Współczynnik dyfuzji Arrheniusa:
  D(T) = D₀ · exp(−Qd / kT)    lub równoważnie exp(−Qd / RT)
  D₀ = czynnik przedwykładniczy [m² s⁻¹]
  Qd = energia aktywacji dyfuzji [eV lub J/mol]

Rozwiązanie nawęglania (ciało półnieskończone, stałe stężenie powierzchniowe Cs):
  (Cs − C(x,t)) / (Cs − C₀) = erf( x / (2√(Dt)) )
  gdzie C₀ = początkowe stężenie w masie materiału

Typowe energie aktywacji (samodyfuzja):
  Fe w Fe-α (BCC):  0,84 eV
  Cu w Cu (FCC):    2,04 eV
  B w Si:           3,46 eV  (powolna — istotna przy produkcji tranzystorów)
          

Symulacja Dyfuzji w krysztale implementuje 2D model mechanizmu wakansowego na sieci kwadratowej, ze współczynnikami przeskoku ustalonymi przez prawo Arrheniusa. Podnieś suwak temperatury i obserwuj, jak fronty dyfuzji rozchodzą się przez początkowo rozdzielony materiał — profil Ficka pojawia się samoistnie z milionów pojedynczych skoków atomowych.

3. Zachowanie mechaniczne — naprężenie, odkształcenie i pęknięcie

Ciągnij pręt metalowy ze wzrastającą siłą, a najpierw rozciągnie się sprężyście (prawo Hooke'a, odwracalnie), potem przejdzie w odkształcenie plastyczne (trwałe, napędzane ruchem dyslokacji), umocni się przez zgniot, gdy dyslokacje się zaplątują, osiągnie wytrzymałość na rozciąganie, wytworzy przewężenie i wreszcie pęknie. Cała krzywa koduje mikrostrukturę stopu: drobniejsze ziarno podnosi granicę plastyczności (prawo Halla-Petcha), podczas gdy cząstki wydzieleń blokują dyslokacje (umacnianie wydzieleniowe, podstawa stopów aluminiowych w lotnictwie).

Zakres sprężysty:     σ = E · ε            (moduł Younga E, zwykle 70–210 GPa dla metali)
Współczynnik Poissona:    ν = −ε_boczne / ε_osiowe   (ν ≈ 0,3 dla większości metali)
Kryterium plastyczności (von Misesa): σ_ef = √½·[(σ₁−σ₂)²+(σ₂−σ₃)²+(σ₃−σ₁)²] = σ_y

Umocnienie potęgowe (plastyczne): σ = σ_y + K · ε_p^n
  K = współczynnik umocnienia, n = wykładnik umocnienia odkształceniowego (0,1–0,5)

Halla-Petcha (umocnienie wielkością ziarna): σ_y = σ₀ + k / √d
  d = średnia średnica ziarna, k ≈ 0,7 MPa·m^½ (stal)

Kryterium pęknięcia Griffitha:
  σ_f = √(2Eγ / πa)    dla pęknięcia przelotowego o połowie długości a
  Odporność na pękanie:   K_Ic = σ_f · Y · √(πa)   [MPa·m^½]
          

Symulacja Naprężenie-odkształcenie animuje próbkę rozciąganą z przewężeniem i pęknięciem w czasie rzeczywistym. Wybierz spośród stali (E = 200 GPa), aluminium (E = 69 GPa), gumy (hipersprężystej), kości (anizotropowej) lub polimeru (pełzającego lepkosprężyście). Pole odporności — całka σdε — aktualizuje się na żywo, więc możesz porównać pochłanianie energii przed pęknięciem. Symulacja Mechaniki pękania dodaje koncentrację naprężeń na czubku pęknięcia: obserwuj, jak pęknięcie rozprzestrzenia się w miarę kumulowania się zmęczenia materiału pod cyklicznym obciążeniem.

4. Dyslokacje — dlaczego metale są dużo słabsze, niż przewiduje teoria

Teoretyczne naprężenie ścinające potrzebne, by przesunąć jedną płaszczyznę krystaliczną względem drugiej, wynosi τ_teor = G/(2π) ≈ 10–30 GPa. Ale rzeczywiste metale ulegają odkształceniu przy 50–500 MPa — 100 razy słabiej. Wyjaśnienie, znalezione niezależnie przez Taylora, Orowana i Polanyi'ego w 1934 roku, to dyslokacja: defekt liniowy, w którym poślizg rozpoczął się, ale jeszcze się nie zakończył. Tylko atomy w pobliżu rdzenia dyslokacji poruszają się w danym momencie; dyslokacja przesuwa się przez kryształ jak fałda wędrująca przez dywan.

Wektor Burgersa b: kierunek i wielkość zniekształcenia sieci
  |b| = a/2 · √(h²+k²+l²) dla systemu poślizgu FCC ⟨110⟩{111}

Pole naprężeń dyslokacji krawędziowej:
  σ_xx = − Gb/(2π(1−ν)) · y(3x²+y²)/(x²+y²)²
  σ_yy =   Gb/(2π(1−ν)) · y(x²−y²)/(x²+y²)²
  τ_xy =   Gb/(2π(1−ν)) · x(x²−y²)/(x²+y²)²

Naprężenie Peierlsa-Nabarro (tarcie sieciowe):
  τ_PN ≈ 2G · exp(−2πw/b)  gdzie w = szerokość dyslokacji

Umocnienie dyslokacyjne (umocnienie Taylora):
  Δσ = M · α · G · b · √ρ
  ρ = gęstość dyslokacji [m⁻²], α ≈ 0,3, M = współczynnik Taylora ≈ 3,06

Namnażanie źródła Franka-Reada: dyslokacje wyginają się między dwoma punktami zakotwiczenia
  i zataczają pełne pętle, gdy τ ≥ Gb/(2L),  L = odstęp punktów zakotwiczenia
          

Symulacja Poślizgu dyslokacji pokazuje dyslokacje krawędziowe i śrubowe przesuwające się przez kryształ 2D pod przyłożonym naprężeniem ścinającym. Podnieś naprężenie i obserwuj, jak dyslokacje się zarodkują, mnożą przez źródła Franka-Reada, spiętrzają się na granicach ziaren i ostatecznie powodują odkształcenie plastyczne. Nakładka kolorystyczna pokazuje lokalną koncentrację naprężeń ścinających wokół każdego rdzenia dyslokacji.

5. Diagramy fazowe — stopy dwuskładnikowe

Dwuskładnikowy diagram fazowy mapuje, które fazy (stała, ciekła, mieszana) są termodynamicznie stabilne w zależności od temperatury i składu. Punkt eutektyczny to para skład-temperatura, przy której ciecz zamarza bezpośrednio w dwie fazy stałe jednocześnie — najniższa temperatura topnienia w układzie. Reguła dźwigni podaje udział każdej fazy w obszarze dwufazowym bez rozwiązywania żadnych równań.

Reguła faz Gibbsa:   F = C − P + 2
  F = liczba stopni swobody (zmienne intensywne swobodnie regulowalne)
  C = liczba składników, P = liczba faz

Przy stałym ciśnieniu (izobarycznie): F = C − P + 1

Dwuskładnikowy układ eutektyczny (C=2):
  Obszar jednofazowy (ciecz lub ciało stałe): F = 2 (temperatura + skład)
  Obszar dwufazowy (likwidus + solidus): F = 1 (ustal T → oba składy ustalone)
  Punkt eutektyczny (3 fazy, ciecz + 2 ciała stałe): F = 0 (niezmienny)

Reguła dźwigni (udział fazy α przy składzie X_0 między Xα i Xβ):
  f_α = (X_β − X_0) / (X_β − X_α)
  f_β = (X_0 − X_α) / (X_β − X_α)

Przemiana martenzytyczna: bezdyfuzyjne ścinanie (austenit Fe-C → martenzyt)
  Temperatura początkowa: Ms ≈ 539 − 423(%C) − 30,4(%Mn) − 17,7(%Ni) − 12,1(%Cr) [°C]
          

Symulacja Diagramu fazowego stopu rysuje pełny dwuskładnikowy diagram fazowy typu eutektycznego Pb-Sn. Kliknij dowolne miejsce na diagramie, a obliczenie reguły dźwigni pojawi się automatycznie — względne długości odcinków linii wiążącej skalują udziały faz. Możesz też symulować ścieżkę chłodzenia: podążaj pionowo w dół wzdłuż danego składu i obserwuj, jak mikrostruktura formuje się na animowanej wstawce.

6. Struktura pasmowa półprzewodników

Dlaczego krzem przewodzi prąd tylko po domieszkowaniu, miedź przewodzi swobodnie, a diament nie przewodzi wcale? Odpowiedź tkwi w strukturze pasmowej elektronów — dozwolonych poziomach energetycznych dla elektronów w sieci periodycznej. W metalach poziom Fermiego znajduje się wewnątrz częściowo wypełnionego pasma. W półprzewodnikach przerwa energetyczna oddziela wypełnione pasmo walencyjne od pustego pasma przewodnictwa; energia cieplna lub fotony (albo atomy domieszki) mogą przekroczyć tę przerwę. W izolatorach przerwa jest zbyt szeroka.

Przerwa energetyczna prawie swobodnego elektronu na granicy strefy k = π/a:
  E_przerwa = 2|V_G|   (V_G = składowa Fouriera potencjału periodycznego)

Rozkład Fermiego-Diraca (prawdopodobieństwo obsadzenia stanu E):
  f(E) = 1 / (exp((E − E_F) / kT) + 1)

Stężenie nośników samoistnych:
  n_i = √(N_c · N_v) · exp(−E_g / 2kT)
  N_c = 2(2πm_e*kT/h²)^(3/2),  N_v analogicznie z m_h*

Domieszkowanie typu n (atomy donorowe z poziomem E_d blisko pasma przewodnictwa):
  Wymrożenie → nasycenie → reżim samoistny wraz ze wzrostem T
  n ≈ N_d (przy umiarkowanym T, pełna jonizacja)

Współczynnik Halla (określenie znaku nośników):
  R_H = −1/(ne)  (typ n)    R_H = +1/(pe) (typ p)

Przerwa prosta a skośna:
  Prosta (GaAs, GaN): przejścia optyczne z zachowaniem pędu → LED, laser
  Skośna (Si, Ge): przejście wymaga fononu → słaby emiter
          

Symulacja Pasm półprzewodnika rysuje dyspersję E–k dla jednowymiarowego modelu Kroniga-Penneya i pozwala regulować głębokość studni potencjału, odstęp i temperaturę. Ogon Fermiego-Diraca wypełnia się w miarę zwiększania temperatury — możesz zaobserwować przejście od reżimu domieszkowego (gęstość nośników ustalona przez domieszki) do samoistnego (gęstość nośników ustalona przez przerwę energetyczną i kT). Osobno panel domieszkowania pokazuje poziomy donorowe i akceptorowe typu n i p leżące tuż wewnątrz przerwy energetycznej.

Dlaczego inżynieria materiałowa zasługuje na więcej symulacji

Inżynieria materiałowa to jedna z dziedzin nauk fizycznych najbardziej przyjaznych symulacji. Dynamika molekularna, analiza metodą elementów skończonych, teoria funkcjonału gęstości, modele pola fazowego — wszystkie te metody są rutynowo stosowane w badaniach materiałowych, a ich wizualne wyniki są od razu intuicyjne. Dyslokacja poruszająca się przez kryształ wygląda tak, jak powinna. Granica fazowa migrująca przez ziarno wygląda tak, jak powinna.

Mimo to inżynieria materiałowa często bywa pomijana w edukacji fizycznej online. Większość kursów wprowadzających przeskakuje bezpośrednio od orbitali atomowych do termodynamiki, pomijając po drodze ciało krystaliczne. Ten spotlight próbuje wypełnić tę lukę symulacjami zakotwiczonymi w prawdziwych parametrach materiałowych i rzeczywistych równaniach — a nie zabawkach.

W przyszłych falach: ewolucja gęstości dyslokacji podczas obróbki na zimno, potencjał pary Lennarda-Jonesa łączący siły międzyatomowe ze stałymi sprężystości, wzrost ziaren podczas wyżarzania (normalny i anormalny) oraz pełna symulacja złącza p-n pokazująca obszar zubożenia i krzywą prądowo-napięciową.