💠 โรงงานเซมิคอนดักเตอร์ โมเดลสด

พื้นฐาน — โมเดลเบื้องหลังตัวเลข

ทุกความสัมพันธ์ที่เครื่องจำลองใช้ พร้อมแหล่งที่มา ค่าคงที่ที่ระบุว่าเป็นสมมติฐานคือการปรับเทียบตัวอย่าง

สถานะอุปกรณ์เป็นไปตาม SEMI E10: ผลิต สแตนด์บาย วิศวกรรม หยุดตามแผน และหยุดไม่ตามแผน

Availability = up tools / 8 · Utilization = productive tool-hours / total[1]

เครื่องมือเสียแบบสุ่มด้วยเวลาก่อนเสียแบบเอ็กซ์โพเนนเชียล งานซ่อมรอช่างว่าง

P(fail in Δt) = 1 − e^(−Δt/MTBF), MTBF = 250 h; repair ~ Exp(mean 6 h) once a technician is free[17][1]สมมติฐาน: ขนาดของค่า (MTBF อัตราการเกิด k τ) เป็นการปรับเทียบตัวอย่าง ไม่ใช่ข้อมูลจากโรงงานจริง

กฎของลิตเติลเชื่อมงานระหว่างทำ อัตราการผลิต และรอบเวลาเข้าด้วยกัน

CT = WIP / TH = T0 + queue / TH, T0 = 1 h[6][5]

ในสายการผลิตที่เสถียร ปริมาณงานออกเท่ากับอัตราการปล่อย ดังนั้นเวเฟอร์ที่ไม่ถูกปล่อยเข้าสายการผลิตคือเวเฟอร์ที่ไม่มีวันออกมา เวเฟอร์ที่รอนานเกินไประหว่างขั้นตอนจะละเมิดหน้าต่างเวลาในคิว

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]สมมติฐาน: ขีดจำกัดรอบเวลา 3.5 ชม. ใช้แทนหน้าต่างเวลาในคิว (queue-time window) หน้าต่างจริงกำหนดโดยวิศวกรรมกระบวนการทีละขั้นตอน และเวเฟอร์ที่ละเมิดมักถูกทิ้งเป็นของเสีย

การประมาณของคิงแมน: เวลาในคิวเพิ่มเหมือน u/(1−u) — พุ่งทะยานเมื่ออัตราใช้งานใกล้ 100 %

CTq ≈ ((ca² + ce²)/2) · (u/(1−u)) · te[5]แสดงเพื่อสร้างความเข้าใจ: คิวของเครื่องจำลองเกิดจากการมาแบบสุ่มและกำลังผลิตแบบสุ่ม ไม่ใช่จากสูตรนี้

ขีดจำกัดชั้น ISO 14644-1 สำหรับอนุภาคขนาด D ที่ 0.1 µm ชั้นคือลอการิทึมของความเข้มข้น

Cn = 10^N · (0.1/D)^2.08 ⇒ N = log10(C≥0.1µm)[2]

ห้องผสมสม่ำเสมอ: ความเข้มข้น = อัตราการเกิด ÷ (การเปลี่ยนอากาศ × ปริมาตร × ประสิทธิภาพฟิลเตอร์)

C = G / (ACH · V · η)[13][11]สมมติฐาน: ขนาดของค่า (MTBF อัตราการเกิด k τ) เป็นการปรับเทียบตัวอย่าง ไม่ใช่ข้อมูลจากโรงงานจริง

กฎความคล้ายของพัดลม: กำลังเพิ่มตามกำลังสามของความเร็วพัดลม

P_fan = P_max · (speed)³[10]

ผลผลิตแบบทวินามลบ: ข้อบกพร่องรวมกลุ่ม ผลผลิตจึงลดช้ากว่าที่โมเดลปัวซงทำนาย

Y = (1 + A·D/α)^(−α), α = 2 (α→∞: Y = e^(−A·D))[3][4][15][14]

การเรียนรู้ผลผลิต: ความหนาแน่นข้อบกพร่องลดลงเข้าสู่ระดับที่โตเต็มที่เมื่อใช้เวลาวิศวกรรม

D(t+Δt) = D∞ + (D − D∞)·e^(−Δt/τ), τ = 120 h · (0.10 / engineering share)[15][18]สมมติฐาน: τ ถูกบีบเวลาราว 100 เท่าเทียบกับอัตราการเรียนรู้ของอุตสาหกรรมที่รายงาน (4–6.5 % ต่อเดือน, Leachman) เพื่อให้สถานการณ์ 14 วันแสดงผลได้

อนุภาคที่ตกบนเวเฟอร์เพิ่มข้อบกพร่องร้ายแรงตามสัดส่วนความเข้มข้นในอากาศ

D_total = D_learn + k · C, k = 4·10⁻⁶ cm⁻² per particle/m³สมมติฐาน: ค่าคงที่ตัวอย่าง — โรงงานจริงปรับเทียบอัตราข้อบกพร่องร้ายแรงจากข้อมูลการตรวจสอบของตนเอง

ผลผลิตเชิงพารามิเตอร์: สัดส่วนของกระบวนการปกติที่อยู่ในขีดจำกัดข้อกำหนด Cpk วัดส่วนเผื่อ

Y_param = Φ((USL−μ)/σ) − Φ((LSL−μ)/σ), Cpk = min(USL−μ, μ−LSL)/(3σ)[7][8][19]

SPC: ค่าเฉลี่ยของตัวอย่าง n ตัว พร้อมขีดจำกัดควบคุม ±3σ/√n และกฎโซน

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]

ค่าคงที่การดำเนินงานอื่น ๆ ที่โมเดลใช้

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 %สมมติฐาน: ขนาดของค่า (MTBF อัตราการเกิด k τ) เป็นการปรับเทียบตัวอย่าง ไม่ใช่ข้อมูลจากโรงงานจริง

จำนวนได (die) รวมต่อเวเฟอร์ 300 mm สำหรับไดที่มีพื้นที่ A

DPW = π·d²/(4A) − π·d/√(2A), d = 300 mm[9]

ความสุ่ม: ตัวสร้าง mulberry32 แบบมีซีด การแจกแจงที่ใช้ — สม่ำเสมอ เอ็กซ์โพเนนเชียล (CDF ผกผัน) ปกติ (Box–Muller) ปัวซง (Knuth) ซีดแสดงให้เห็นและแชร์ได้

แหล่งที่มา

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