Fabrică de semiconductori Model live

Conduci bay-ul de litografie al unei fabrici de wafere de 300 mm — opt scannere, o cameră curată și o rețetă care trebuie să rămână în specificații. Waferele se adună la coadă, utilajele se strică, particulele pătrund și randamentul se schimbă la fiecare decizie.

Ce vei învăța

Simulator

Timp 0 h
Disponibilitatea utilajelor 100% · Lucru în curs 20 · Clasa ISO a camerei curate 3,4S1⚙S2▶S3▶S4▶S5▶S6▶S7‖S8‖WIP 20▶ 34 wph
  • Productiv
  • Standby
  • Inginerie
  • Oprit (reparație)
  • Particule din aer

Comenzi

Waferele lansate în bay pe oră.

Partea din timpul scannerului acordată loturilor experimentale de plachete care învață procesul: mai puțină producție acum, mai puține defecte mai târziu.

Fiecare tehnician repară un utilaj pe rând.

Mai multe schimburi de aer diluează particulele; puterea ventilatoarelor crește cu cubul vitezei.

Fiecare persoană echipată continuă să elibereze particule.

Waferele măsurate pe oră: diagramă SPC mai clară, capacitate puțin mai mică.

Durează 3 h și adaugă 2 persoane în încăpere cât timp lucrează.

2 h la jumătate din capacitate, apoi procesul revine la țintă.

Indicatori

Cipuri bune produse
19547/h
Randament
89,8%
normal
Clasa ISO a camerei curate
3,4
normal
Timp de ciclu al bay-ului
1,6h
normal
Disponibilitatea utilajelor100 %
Utilizarea utilajelor71 %
Timp de inginerie10 %
Wafer starts lansate34 wph
Lucru în curs20 wafers
Densitatea defectelor0,110 /cm²
Capabilitatea procesului Cpk1,25 · avertizare
În afara specificației (echivalent wafere)0,0 wph
Puterea bay-ului1339 kW

Tendință

Cipuri bune produse: — /h250000

Eșantion SPC (abatere CD)

Eșantion SPC (abatere CD): 0 ⚑+3σ−3σ

Scenarii de criză

Nivel 1 · Excursie de particule

O tură de rutină pe linia de cipuri mobile. Undeva deasupra bay-ului un filtru de tavan este pe punctul de a ceda. Menține camera curată în clasă, protejează randamentul — și nu lăsa ventilatoarele la viteză maximă mai mult decât este necesar.

  • Clasa ISO revine la ≤ 3,6 la final
  • Niciodată peste clasa 4,55 după prima oră
  • Randament mediu ≥ 87 % în primele ore ale excursiei
  • Puterea medie a bay-ului ≤ 1.380 kW în a doua jumătate

Nivel 2 · Defectarea scannerelor

Bay-ul funcționează la limită, cu 38 de wafer starts pe oră. Trei scannere sunt pe punctul de a se opri simultan. Nu lăsa coada să explodeze cât timp le readuci și continuă să livrezi cipuri bune.

  • WIP ≤ 60 de wafere la final
  • Timp de ciclu niciodată peste 3,5 h după defectare (fereastra de timp de așteptare)
  • Cel puțin 900.000 de cipuri bune în 48 h
  • Lansează cel puțin 1600 de wafer starts în 48 h (planul de lansări: 34 pe oră)

Nivel 3 · Deriva rețetei

O nouă revizuire a rețetei a intrat în producție aseară. Nimeni nu știe încă că deplasează lent dimensiunea critică. Urmărește diagrama SPC: prinde deriva și revino la rețeta anterioară înainte ca waferele să iasă din specificații.

  • Cel mult 12 echivalente wafer în afara specificației în total
  • Cpk ≥ 1,2 la final
  • Cel puțin 620.000 de cipuri bune

Nivel 2 · Creșterea unui produs nou

Un cip accelerator mare intră în producție cu un proces imatur: densitatea defectelor este mare, iar un cip de 600 mm² este necruțător. Linia este încărcată la 42 de wafer starts pe oră. Ai 14 zile. Echilibrează producția de azi cu învățarea pentru mâine.

  • Densitate de defecte ≤ 0,08 /cm² în ziua 14
  • Cel puțin 660.000 de cipuri bune în 14 zile

Baza — modelul din spatele cifrelor

Fiecare relație folosită de simulator, cu sursa ei. Constantele marcate ca ipoteze sunt calibrări ilustrative.

Stările echipamentelor urmează SEMI E10: productiv, standby, inginerie, oprire programată și neprogramată.
Availability = up tools / 8 · Utilization = productive tool-hours / total[1]
Utilajele se defectează aleatoriu, cu un timp până la defectare exponențial; reparațiile așteaptă un tehnician liber.
P(fail in Δt) = 1 − e^(−Δt/MTBF), MTBF = 250 h; repair ~ Exp(mean 6 h) once a technician is free[17][1]Ipoteză: ordinele de mărime (MTBF, ratele de generare, k, τ) sunt calibrări ilustrative, nu date dintr-o fabrică reală.
Legea lui Little leagă lucrul în curs, debitul și timpul de ciclu.
CT = WIP / TH = T0 + queue / TH, T0 = 1 h[6][5]
Într-o linie stabilă, debitul este egal cu rata de lansare, deci un wafer start nelansat este un wafer care nu iese niciodată. Waferele care așteaptă prea mult între pași depășesc fereastra de timp de așteptare.
TH = release rate while u < 1 ⇒ wafers out ≈ Σ starts; queue wait = CT − T0 = queue / TH ≤ queue-time window (3.5 h CT in semi-tool-crash)[5][16]Ipoteză: limita de 3,5 h a timpului de ciclu ține locul unei ferestre de timp de așteptare (queue-time); ferestrele reale sunt stabilite pas cu pas de ingineria de proces, iar waferele care le depășesc sunt de regulă rebutate.
Aproximarea Kingman: timpul de coadă crește ca u/(1−u) — explodează când utilizarea se apropie de 100 %.
CTq ≈ ((ca² + ce²)/2) · (u/(1−u)) · te[5]Arătat pentru intuiție: coada simulatorului rezultă din sosiri aleatorii și capacitate aleatorie, nu din această formulă.
Limita clasei ISO 14644-1 pentru particule de dimensiune D; la 0,1 µm clasa este logaritmul concentrației.
Cn = 10^N · (0.1/D)^2.08 ⇒ N = log10(C≥0.1µm)[2]
Cameră bine amestecată: concentrație = generare ÷ (schimburi de aer × volum × eficiența filtrului).
C = G / (ACH · V · η)[13][11]Ipoteză: ordinele de mărime (MTBF, ratele de generare, k, τ) sunt calibrări ilustrative, nu date dintr-o fabrică reală.
Legea de afinitate a ventilatoarelor: puterea crește cu cubul vitezei ventilatorului.
P_fan = P_max · (speed)³[10]
Randament binomial negativ: defectele se grupează, deci randamentul scade mai lent decât prezice modelul Poisson.
Y = (1 + A·D/α)^(−α), α = 2 (α→∞: Y = e^(−A·D))[3][4][15][14]
Învățarea randamentului: densitatea defectelor scade spre un nivel matur pe măsură ce se consumă timp de inginerie.
D(t+Δt) = D∞ + (D − D∞)·e^(−Δt/τ), τ = 120 h · (0.10 / engineering share)[15][18]Ipoteză: τ este comprimat în timp de circa 100× față de ratele raportate de învățare din industrie (4–6,5 % pe lună, Leachman), astfel încât un scenariu de 14 zile să arate efectul.
Particulele care aterizează pe wafere adaugă defecte ucigașe proporțional cu concentrația din aer.
D_total = D_learn + k · C, k = 4·10⁻⁶ cm⁻² per particle/m³Ipoteză: o constantă ilustrativă — fabricile reale calibrează ratele defectelor ucigașe din propriile date de inspecție.
Randament parametric: ponderea unui proces normal care se încadrează în limitele specificației; Cpk măsoară marja.
Y_param = Φ((USL−μ)/σ) − Φ((LSL−μ)/σ), Cpk = min(USL−μ, μ−LSL)/(3σ)[7][8][19]
SPC: medii de n eșantioane, cu limite de control ±3σ/√n și reguli de zonă.
x̄ ~ N(μ, σ/√n), control limits ±3σ/√n; zone rules (2 of 3 beyond 2σ, 4 of 5 beyond 1σ, 8 on one side)[8][7][12]
Alte constante de funcționare folosite de model.
8 scanners × 6 wafer-starts/h · metrology −0.5 % capacity per sampled wafer/h · rollback: 2 h at half capacity + 4 wafer-equivalents reworked · drift 0.25 nm/h · crash repair 18 technician-hours · leak crew +2 people for 3 h · tools 150 kW busy / 60 kW idle, fans 150 kW at 100 %Ipoteză: ordinele de mărime (MTBF, ratele de generare, k, τ) sunt calibrări ilustrative, nu date dintr-o fabrică reală.
Cipuri brute (dies) per wafer de 300 mm pentru un cip cu aria A.
DPW = π·d²/(4A) − π·d/√(2A), d = 300 mm[9]

Aleatoriu: un generator mulberry32 cu seed; distribuțiile folosite — uniformă, exponențială (CDF inversă), normală (Box–Muller), Poisson (Knuth). Seed-ul este afișat și poate fi partajat.

Surse

  1. SEMI E10 — Specification for Definition and Measurement of Equipment Reliability, Availability, and Maintainability (RAM) and Utilization — SEMI
  2. ISO 14644-1:2015 Cleanrooms and associated controlled environments — Part 1: Classification of air cleanliness by particle concentration — ISO, 2015
  3. C. H. Stapper, F. M. Armstrong, K. Saji — Integrated circuit yield statistics — Proceedings of the IEEE 71(4), 1983
  4. J. A. Cunningham — The use and evaluation of yield models in integrated circuit manufacturing — IEEE Trans. Semiconductor Manufacturing 3(2), 1990
  5. W. J. Hopp, M. L. Spearman — Factory Physics (3rd ed.), ch. 7–8: Little’s law, Kingman (VUT) equation — Waveland Press, 2008
  6. J. D. C. Little — A Proof for the Queuing Formula L = λW — Operations Research 9(3), 1961
  7. D. C. Montgomery — Introduction to Statistical Quality Control (x̄ charts, process capability Cpk) — Wiley, 2019
  8. NIST/SEMATECH e-Handbook of Statistical Methods — 6.3 Univariate and Multivariate Control Charts; 6.1.6 Process capability — NIST
  9. Dies-per-wafer estimate DPW = πd²/(4S) − πd/√(2S) (de Vries, “Investigation of gross die per wafer formulas”, IEEE TSM 18(1)) — IEEE, 2005
  10. Fan affinity laws: flow ∝ speed, pressure ∝ speed², power ∝ speed³ — U.S. DOE — Improving Fan System Performance: A Sourcebook for Industry
  11. EN 1822-1:2019 High efficiency air filters (EPA, HEPA and ULPA) — classification (U15 ≥ 99.9995 % at MPPS) — CEN, 2019
  12. Western Electric Statistical Quality Control Handbook (1956) — zone rules for control charts — Western Electric / NIST e-Handbook 6.3.2
  13. W. Whyte — Cleanroom Technology: Fundamentals of Design, Testing and Operation (2nd ed.), ch. on dispersion of particles from people; well-mixed room dilution equation — Wiley, 2010
  14. Yu. I. Bogdanov, N. A. Bogdanova, V. L. Dshkhunyan — Statistical Yield Modeling for IC Manufacture: Hierarchical Fault Distributions (§2 Compound Poisson distribution, after Eq. (17), p. 4 of the arXiv PDF — large-area clustering negative binomial model: the cluster parameter’s “typical values approximately range from 0.3 to 7”) — arXiv physics/0303039, 2003
  15. R. C. Leachman — Yield Modeling and Analysis (Poisson, Murphy, Seeds, negative-binomial models; §8: SMLY survey with C. N. Berglund for International SEMATECH, 2002–03 — yield loss fitted as YL(t) = YL(0)·e^(−λt), Table 2 averages 4.4 / 4.0 / 6.5 %/month at 350 / 250 / 180 nm) — UC Berkeley, IEOR 130 course notes (unpublished), 2014
  16. A. Klemmt, L. Mönch — Scheduling jobs with time constraints between consecutive process steps in semiconductor manufacturing (time windows set by process engineering against native oxidation and contamination; jobs that violate them are scrapped) — Proceedings of the 2012 Winter Simulation Conference, 2012
  17. NIST/SEMATECH e-Handbook of Statistical Methods — 8.1.6.1 Exponential distribution (constant failure rate) — NIST
  18. C. Weber — Yield learning and the sources of profitability in semiconductor manufacturing and process development — IEEE Trans. Semiconductor Manufacturing 17(4), 2004
  19. M. Abramowitz, I. Stegun — Handbook of Mathematical Functions, 26.2.17 — NBS, 1964

Cine face asta ca meserie

Model educațional — nu pentru decizii operaționale. Site-urile reale calibrează fiecare constantă la propriile echipamente și date.