مصنع أشباه الموصلات نموذج مباشر

تدير قسم الطباعة الضوئية في مصنع رقائق 300 mm — ثمانية ماسحات، وغرفة نظيفة، ووصفة يجب أن تبقى ضمن المواصفات. تصطف الرقائق، وتتعطل الأدوات، وتتسلل الجسيمات، ويتحرك المردود مع كل قرار.

ما الذي ستتعلمه

المحاكي

الوقت 0 h
جاهزية الأداة 100% · العمل قيد التنفيذ 20 · فئة ISO للغرفة النظيفة 3.4S1⚙S2▶S3▶S4▶S5▶S6▶S7‖S8‖WIP 20▶ 34 wph
  • منتج
  • انتظار
  • هندسة
  • متوقف (إصلاح)
  • جسيمات محمولة جوًا

الضوابط

الرقائق المطلقة إلى القسم في الساعة.

حصة وقت الماسح الضوئي المخصصة لدفعات رقائق تجريبية تُعلِّم العملية: إنتاج أقل الآن، وعيوب أقل لاحقًا.

يصلح كل فني أداة واحدة في كل مرة.

مزيد من تبديلات الهواء يخفف الجسيمات؛ وتزداد قدرة المروحة مع مكعب السرعة.

كل شخص يرتدي زيّ الغرفة النظيفة ما زال يطلق جسيمات.

رقائق تُقاس في الساعة: مخطط SPC أدق، وسعة أقل قليلًا.

يستغرق 3 h ويضيف شخصين في الغرفة أثناء عملهما.

ساعتان بنصف السعة، ثم تعود العملية إلى الهدف.

المؤشرات

الشرائح الجيدة الخارجة
19547/h
المردود
89.8%
عادي
فئة ISO للغرفة النظيفة
3.4
عادي
زمن دورة القسم
1.6h
عادي
جاهزية الأداة100 %
استخدام الأداة71 %
وقت الهندسة10 %
بدايات الرقائق المُطلقة34 wph
العمل قيد التنفيذ20 wafers
كثافة العيوب0.110 /cm²
قدرة العملية Cpk1.25 · تحذير
خارج المواصفات (مكافئ رقائق)0.0 wph
قدرة القسم1339 kW

الاتجاه

الشرائح الجيدة الخارجة: — /h250000

عينة SPC (انزياح CD)

عينة SPC (انزياح CD): 0 ⚑+3σ−3σ

سيناريوهات الأزمات

المستوى 1 · انحراف الجسيمات

مناوبة روتينية على خط شرائح الجوال. في مكان ما فوق القسم يوشك مرشّح سقفي على الفشل. أبقِ الغرفة النظيفة ضمن فئتها، واحمِ المردود — ولا تترك المراوح بأقصى سرعة أكثر مما تحتاج.

  • عودة فئة ISO إلى ≤ 3.6 في النهاية
  • لا تتجاوز الفئة 4.55 بعد الساعة الأولى
  • متوسط المردود ≥ 87 % في الساعات الأولى من الانحراف
  • متوسط قدرة القسم ≤ 1,380 kW في النصف الثاني

المستوى 2 · انهيار الماسحات

يعمل القسم بحرارة عالية عند 38 بداية رقاقة في الساعة. ثلاث ماسحات على وشك التوقف معًا. امنع انفجار الطابور وأنت تعيدها، وسلّم مع ذلك شرائح جيدة.

  • WIP ≤ 60 رقاقة في النهاية
  • زمن الدورة لا يتجاوز 3.5 h أبدًا بعد الانهيار (نافذة زمن الانتظار)
  • ما لا يقل عن 900,000 شريحة جيدة في 48 h
  • أطلِق ما لا يقل عن 1,600 بداية رقاقة خلال 48 h (خطة البدايات: 34 في الساعة)

المستوى 3 · انحراف الوصفة

دخلت مراجعة جديدة للوصفة الخدمة مساء أمس. لا أحد يعلم بعد أنها تُزيح البعد الحرج ببطء. راقب مخطط SPC: التقط الانحراف وتراجع قبل أن تخرج الرقائق عن المواصفات.

  • ما لا يزيد عن 12 مكافئ رقاقة خارج المواصفات إجمالًا
  • Cpk ≥ 1.2 في النهاية
  • ما لا يقل عن 620,000 شريحة جيدة

المستوى 2 · تصاعد إنتاج منتج جديد

تدخل شريحة مسرّع كبيرة الإنتاج بعملية غير ناضجة: كثافة العيوب عالية وشريحة 600 mm² لا ترحم. يُحمَّل الخط عند 42 بداية رقاقة في الساعة. أمامك 14 يومًا. وازن بين إنتاج اليوم والتعلّم لأجل الغد.

  • كثافة العيوب ≤ 0.08 /cm² في اليوم 14
  • ما لا يقل عن 660,000 شريحة جيدة خلال 14 يومًا

الأساس — النموذج وراء الأرقام

كل علاقة يستخدمها المحاكي مع مصدرها. الثوابت المعلَّمة كافتراضات هي معايرات توضيحية.

تتبع حالات المعدات 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 h يقوم مقام نافذة زمن الانتظار؛ أما النوافذ الحقيقية فتحددها هندسة العمليات خطوة بخطوة، والرقائق التي تتجاوزها تُتلف عادة.
تقريب كينغمان: ينمو زمن الطابور مثل 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 وτ) معايرات توضيحية، وليست بيانات من مصنع حقيقي.
إجمالي الشرائح لكل رقاقة 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

من يعمل في هذا المجال

نموذج تعليمي — ليس لاتخاذ قرارات تشغيلية. تعاير المواقع الحقيقية كل ثابت وفق معداتها وبياناتها.