منهجية تصميم التجارب لتطوير عمليات الطباعة ثلاثية الأبعاد: من الفرز المتعامد ونمذجة سطح الاستجابة إلى معايرة نافذة العملية

مقدمة: لماذا يدور تحسين المعلمات في حلقة مفرغة

في الضبط الميداني لعمليات صهر سرير المسحوق بالليزر للمعادن (LPBF) وتلبيد البوليمرات، يتكرر السيناريو التالي: يرفع المهندس قدرة الليزر من 250 W إلى 280 W، فتزداد الكثافة بنسبة 0.3%، ثم يواصل الرفع؛ وعند 300 W تظهر ظاهرة ثقب المفتاح، فيعود إلى خفض القدرة. وبعد جولة كاملة تضم أكثر من عشرة عينات، تكون النتيجة منحنى خبرة مجزأ يجب البدء به من جديد عند تغيير الماكينة أو دفعة المسحوق. لا تكمن المشكلة في نقص جهد المهندس، بل في أن منهج تغيير عامل واحد بالتناوب لا يستطيع التعامل مع التأثيرات التفاعلية؛ إذ تحدد القدرة وسرعة المسح معاً كثافة الطاقة الحجمية، كما يغير سمك طبقة المسحوق ظروف تبديد الحرارة في بركة الانصهار. هذه العوامل الثلاثة مترابطة، ولن تكشف معالجة متغير واحد في كل مرة الصورة الكاملة.

يوفر تصميم التجارب (Design of Experiments، DoE) المسار الهندسي لمعالجة هذا النوع من الاقتران متعدد العوامل: تقدير التأثيرات الرئيسية، وتأثيرات التفاعل، وسطح الاستجابة دفعة واحدة بأقل عدد من العينات، مع تقديم نافذة عملية قابلة للاعتماد مصحوبة بفواصل ثقة. تقدم هذه المقالة منهجية تنفيذية متكاملة، موجهة لتطوير عمليات الطباعة ثلاثية الأبعاد، بدءاً من الفرز المتعامد ونمذجة سطح الاستجابة وصولاً إلى التحقق من النافذة.

أولاً: التمييز بين الأنواع الثلاثة للتجارب: الفرز والنمذجة والتحقق

لا يكتمل تطوير العملية بتجربة واحدة، بل يتقدم عادة عبر ثلاث طبقات، تختلف كل منها في الهدف وعدد العوامل وعدد العينات.

الطبقة الأولى هي تجارب الفرز، وهدفها انتقاء 2 إلى 4 عوامل مؤثرة فعلاً من بين 6 إلى 10 عوامل مرشحة. تشمل العوامل النموذجية قدرة الليزر، وسرعة المسح، وتباعد مسارات المسح، وسمك الطبقة، ودرجة حرارة التسخين المسبق، واستراتيجية المسح (عرض الشريط)، وزاوية الدوران بين الطبقات، ومعدل تدفق الغاز الواقي. تستخدم مرحلة الفرز تصميماً عاملياً جزئياً أو تصميم بلاكيت-بورمان، ويمكن إنجازها في 8 إلى 12 تجربة، مع الاكتفاء بالكثافة وخشونة السطح كمؤشرين عامين للاستجابة.

الطبقة الثانية هي تجارب النمذجة، وهدفها بناء نموذج سطح استجابة للعاملين أو العوامل الثلاثة الرئيسية التي تم فرزها. ويُستخدم عادة تصميم مركب مركزي (CCD) أو تصميم بوكس-بهنكن. يتضمن تصميم CCD ثلاثي العوامل 6 نقاط محورية و8 نقاط مكعبة و6 مكررات للنقطة المركزية، أي 20 تجربة إجمالاً. وتتوسع الاستجابات لتشمل الكثافة، ومقاومة الشد، والانحراف الأبعادي، وخشونة السطح Ra.

الطبقة الثالثة هي تجارب التحقق، حيث تُطبع 3 إلى 5 نقاط تأكيدية مكررة ضمن المجال الأمثل الذي يتنبأ به النموذج، للتحقق من أن الانحراف بين القيم المتوقعة والمقاسة يقع ضمن ±5%. وبعد نجاح التحقق، تُضيّق حدود النافذة بنسبة 10% لتصبح نطاق المعلمات في بطاقة عملية الإنتاج الكمي.

ثانياً: التجارب المتعامدة: تغطية أكبر فضاء بأقل عدد من العينات

عندما يكون عدد العوامل محدوداً وعدد المستويات متساوياً، يكون الجدول المتعامد الخيار الأكثر اقتصاداً. فمثلاً، في LPBF للفولاذ المقاوم للصدأ 316L، عند دراسة ثلاثة عوامل بثلاثة مستويات هي قدرة الليزر P، وسرعة المسح v، وتباعد المسارات h، يتطلب التصميم العاملي الكامل 27 تجربة، بينما لا يحتاج جدول L9(3^4) المتعامد إلا إلى 9 تجارب.

التصميم المحدد هو: P عند 200/240/280 W، وv عند 800/1000/1200 mm/s، وh عند 0.09/0.11/0.13 mm. بعد ترتيب التجارب وفق L9، تُطبع مجموعة من عينات مكعبة بحجم 10×10×10 mm لكل توليفة مستويات، وتُقاس الكثافة بطريقة أرخميدس مع تسجيل خشونة السطح Ra في الوقت نفسه.

يمكن استخدام تحليل المدى لترتيب العوامل بسرعة: يُحسب متوسط الاستجابة لكل عامل عند كل مستوى، ويكون المدى R = max(متوسط الاستجابة) − min(متوسط الاستجابة). وكلما كبر R كان تأثير العامل أكثر أهمية. فإذا كانت النتائج R(P)=2.1% وR(v)=1.4% وR(h)=0.6%، يكون ترتيب الأهمية: القدرة > السرعة > التباعد. ثم يُجرى تحليل التباين (ANOVA)، ويُستخدم اختبار F للحكم على الدلالة؛ إذ تنتقل العوامل التي تزيد قيمة F فيها على القيمة الحرجة وتقل قيمة p فيها عن 0.05 إلى جولة النمذجة التالية.

يتمثل قيد التجارب المتعامدة في افتراضها عدم وجود تفاعل بين العوامل. وإذا وُجد اشتباه باقتران قوي بين القدرة والسرعة، لأن كثافة الطاقة الحجمية E = P/(v·h·t) هي بحد ذاتها علاقة ضرب، فيجب الانتقال إلى تصميم سطح الاستجابة في مرحلة النمذجة بدلاً من الاستمرار في تطبيق الجدول المتعامد قسراً.

ثالثاً: منهجية سطح الاستجابة: من النقاط المنفصلة إلى سطح مستمر

جوهر منهجية سطح الاستجابة (RSM) هو استخدام متعدد حدود تربيعي لملاءمة العلاقة بين الاستجابة والعوامل:

Y = β0 + Σβi·xi + Σβii·xi² + Σβij·xi·xj + ε

ويمثل الحد المتقاطع βij التعبير الكمي عن تأثير التفاعل. في تصميم CCD مثلاً، تكون المعلمة الرئيسية في مخطط العوامل الثلاثة هي المسافة المحورية α. ولتصميم قابل للدوران، تكون α = 2^(k/4) = 2^(3/4) ≈ 1.682، حيث k هو عدد العوامل. وتقع النقاط المحورية عند ±1.682 في الفضاء المرمّز، والنقاط المكعبة عند ±1، بينما تتكرر النقطة المركزية 6 مرات لتقدير الخطأ الصرف.

بعد إتمام الملاءمة، يجب إجراء ثلاثة اختبارات: أولاً اختبار عدم الملاءمة (lack-of-fit)، حيث تشير قيمة p الأكبر من 0.05 إلى عدم وجود انحراف معنوي بين صيغة النموذج والبيانات؛ ثانياً معامل التحديد R² ومعامل R² المعدل، وينبغي أن يكون الفرق بينهما أقل من 0.2، وإلا فقد يكون هناك فرط ملاءمة؛ وثالثاً فحص طبيعية البواقي وتجانس التباين. ومن الحالات الشائعة عملياً بلوغ R² قيمة 0.95 بينما لا يتجاوز R² المعدل 0.86؛ وعندها ينبغي حذف الحدود غير المعنوية وإعادة الملاءمة.

يشكل تحسين الاستجابات المتعددة تحدياً في الطباعة ثلاثية الأبعاد: يجب تعظيم الكثافة، وتقليل الخشونة، وجعل الانحراف الأبعادي قريباً من الصفر، وهي أهداف غالباً ما تتعارض فيما بينها. هندسياً، تُستخدم دالة الرغوبية (desirability function) لتحويل كل استجابة إلى المجال 0—1، ثم يؤخذ المتوسط الهندسي للحصول على الرغوبية الكلية D، وتُحدد توليفة المعلمات عند تعظيم D. وتوزع الأوزان وفق أولوية قبول العميل؛ فمثلاً يمكن للأجزاء الجوية تخصيص وزن 0.5 للكثافة، و0.3 للخشونة، و0.2 للانحراف الأبعادي.

رابعاً: معايرة نافذة العملية: تحويل النقطة المثلى إلى نطاق قابل للاعتماد

العثور على النقطة المثلى ليس سوى الخطوة الأولى، إذ يحتاج الإنتاج الكمي إلى نطاق. ومن الأساليب العملية استخدام مخطط الخطوط الكنتورية المتراكبة (overlay contour plot): تُرسم حدود القبول لكل استجابة على شكل خطوط كنتورية، مثل الكثافة ≥ 99.5%، وRa ≤ 12 μm، والانحراف الأبعادي ≤ ±0.1 mm، ثم تُؤخذ منطقة تداخل جميع الخطوط باعتبارها نافذة العملية.

بعد تحديد النافذة، يجب ترك هامش أمان. والممارسة الموصى بها هي تقليص حدود النافذة إلى الداخل بنسبة 10%. فعلى سبيل المثال، إذا أعطى النموذج نطاق قدرة مقبولاً من 235 إلى 265 W، فيُكتب في بطاقة العملية 238 إلى 262 W. ويرجع ذلك إلى أن تقلب دفعات المسحوق، وانجراف نقاء الغاز الواقي، وتناقص قدرة النظام البصري قد تحرك نقطة العملية الفعلية تدريجياً؛ وترك الهامش يمكن أن يرفع مؤشر قدرة العملية Cpk من 1.0 إلى أكثر من 1.33.

يُنفذ التحقق من النافذة عبر ثلاث دفعات تأكيدية، تضم كل دفعة 5 قطع: تُطبع مجموعة عند مركز النافذة، وأخرى قرب الجانب الداخلي من الحد العلوي، وثالثة قرب الجانب الداخلي من الحد السفلي. وإذا كانت الدفعات الثلاث المتتالية كلها مطابقة، وكان انجراف متوسطات الدفعات أقل من ثلث الانحراف المعياري للاستجابة، تُعد النافذة مستقرة. كما تُسجل بيانات التحقق في سجل العملية لتكون خط أساس للتتبع اللاحق.

خامساً: الأخطاء الشائعة وقائمة التحقق

الخطأ الأول: اعتبار عدد العينات كامل التكلفة. في الواقع، تكون تكلفة 9 تجارب متعامدة إضافة إلى 20 تجربة CCD أقل بكثير من 30 جولة عشوائية لتغيير عامل واحد، كما أن الناتج يكون نموذجاً قابلاً لإعادة الاستخدام.

الخطأ الثاني: إهمال تأثيرات الكتل. إذا أُنجزت تجارب CCD العشرون خلال ثلاثة أيام، فقد تختلط تغيرات حالة المعدة ودفعة المسحوق ودرجة الحرارة والرطوبة المحيطتين مع البواقي. والطريقة الصحيحة هي إدراج التاريخ كعامل كتلة في التصميم، أو ترتيب مكررات للنقطة المركزية داخل كل كتلة لمراقبة الانجراف.

الخطأ الثالث: النظر إلى الكثافة فقط. فقد لا تكون نقطة العملية التي تحقق كثافة 99.5% مطابقة من حيث الخواص الميكانيكية، وقد تُقلل اختبارات الكثافة من تقدير عيوب عدم الانصهار. ويوصى بقياس مقاومة الشد والاستطالة بالتزامن على الأقل، وإضافة الفحص المعدني المجهري والفحص العيني بالتصوير المقطعي الصناعي CT عند الحاجة.

قائمة تحقق التنفيذ: هل تغطي مستويات العوامل نطاق التشغيل المستقر للمعدة؛ هل رُتبت مكررات للنقطة المركزية لتقدير الخطأ الصرف؛ هل أُجري اختبار عدم الملاءمة؛ هل تتوافق أوزان الاستجابات المتعددة مع معايير قبول العميل؛ هل تحتوي النافذة على هامش 10%؛ وهل تغطي دفعات التحقق حدود النافذة.

الخلاصة

يحوّل تصميم التجارب تطوير عمليات الطباعة ثلاثية الأبعاد من محاولات تجريبية قائمة على الخبرة إلى نشاط هندسي قابل للقياس والتكرار: تنجز التجارب المتعامدة فرز العوامل في 9 إلى 12 تجربة، بينما تبني منهجية سطح الاستجابة نموذجاً تربيعياً يتضمن تأثيرات التفاعل في نحو 20 تجربة، وتعالج دالة الرغوبية تعارض الأهداف المتعددة، وتقدم الخطوط الكنتورية المتراكبة نافذة عملية قابلة للاعتماد بهامش أمان. وبالنسبة إلى منصات إدارة دورة الحياة الكاملة مثل lantu3D، التي تمتد من التصميم إلى تسليم المنتج الفعلي، لا تقتصر قيمة DoE على ضبط المعلمات؛ بل إنها ترسخ معرفة العملية في نماذج معلمات قابلة للبحث وإعادة الاستخدام والتحقق، بحيث يصبح كل تطوير للعملية نقطة انطلاق للتطوير التالي، لا بداية جديدة من الصفر.