تحسين الهيكل متعدد المقاييس مع قيود الإجهاد للتصنيع الإضافي

جدول المحتويات

• 1. مقدمة
• 2. صياغة المشكلة
  • 2.1 نموذج الطور الميداني
  • 2.2 قيد الإجهاد
• 3. شروط المثالية
  • 3.1 الشروط الضرورية من الدرجة الأولى
  • 3.2 تحليل الحساسية المرافق
• 4. التنفيذ العددي
  • 4.1 نظرة عامة على الخوارزمية
  • 4.2 مثال على عارضة ناتئة ثنائية الأبعاد
• 5. النتائج والمناقشة
  • 5.1 دراسة الحساسية
  • 5.2 سير عمل الطباعة ثلاثية الأبعاد
• 6. التحليل الأصلي
• 7. التفاصيل التقنية
• 8. النتائج التجريبية
• 9. دراسة حالة: عارضة ناتئة
• 10. التطبيقات المستقبلية
• 11. المراجع

1. مقدمة

التصنيع المضاف (AM)، مثل الطباعة ثلاثية الأبعاد، يُحدث ثورة في التصميم والإنتاج عبر مجالات الهندسة المعمارية والطب والهندسة. تقدم هذه الورقة منهجًا للحقل الطوري للتحسين الطوبولوجي الهيكلي المُصمم خصيصًا لعمليات التصنيع المضاف، مع دمج قيود الإجهاد وقدرات المواد متعددة المقاييس. تستمد الطريقة بدقة شروط المثالية الضرورية من الدرجة الأولى وتقدم خوارزمية عددية للتنفيذ العملي.

2. صياغة المشكلة

2.1 نموذج الطور الميداني

تستخدم طريقة الحقل الطوري حقلاً قياسيًا $\phi(\mathbf{x})$ لتمثيل توزيع المادة، حيث يشير $\phi = 1$ إلى المادة الصلبة و $\phi = 0$ إلى الفراغ. تعمل مسألة التحسين على تقليل الامتثال (المرونة) مع مراعاة قيد الحجم وقيد الإجهاد. تُعطى طاقة الوضع الكلية بواسطة:

$$\Pi(\mathbf{u}, \phi) = \int_\Omega \psi(\varepsilon(\mathbf{u}), \phi) \, d\Omega - \int_{\partial\Omega_N} \mathbf{t} \cdot \mathbf{u} \, dS$$

حيث $\mathbf{u}$ هو مجال الإزاحة، و $\varepsilon$ هو موتر الانفعال، و $\mathbf{t}$ هو الجر على حدود نيومان.

2.2 قيد الإجهاد

يتمثل الابتكار الرئيسي في تضمين قيد الإجهاد لمنع الفشل أثناء عملية التصنيع الإضافي. تمت صياغة قيد الإجهاد على النحو التالي:

$$g(\sigma) = \frac{1}{|\Omega|} \int_\Omega \left( \frac{\sigma_{vm}(\mathbf{x})}{\sigma_y} - 1 \right)^+ \, d\Omega \leq 0$$

حيث $\sigma_{vm}$ هو إجهاد فون ميزس و $\sigma_y$ هو إجهاد الخضوع. يضمن هذا القيد بقاء الإجهاد أقل من حد الخضوع للمادة في جميع أنحاء الهيكل.

3. شروط المثالية

3.1 الشروط الضرورية من الدرجة الأولى

يتم حل مشكلة التحسين باستخدام نهج لاغرانج. يتم اشتقاق الشروط الضرورية من الدرجة الأولى عن طريق أخذ تغيرات دالة لاغرانج بالنسبة لمتغيرات الحالة $\mathbf{u}$، ومتغير التحكم $\phi$، ومضاعفات لاغرانج. يتضمن النظام الناتج معادلة الحالة، والمعادلة المساعدة، وشرط المثالية.

3.2 تحليل الحساسية المرافق

يتم حساب حساسية دالة الهدف بالنسبة لمتغير المجال الطوري باستخدام الطريقة المساعدة. يتم تعريف المشكلة المساعدة على النحو التالي:

$$\int_\Omega \mathbb{C} \varepsilon(\mathbf{v}) : \varepsilon(\mathbf{w}) \, d\Omega = \int_\Omega \frac{\partial \psi}{\partial \phi} \delta \phi \, d\Omega$$

حيث $\mathbf{w}$ هو مجال الإزاحة المرافق. وهذا يتيح حساب التدرجات بكفاءة للمشكلات واسعة النطاق.

4. التنفيذ العددي

4.1 نظرة عامة على الخوارزمية

تستخدم الخوارزمية العددية تقطيع العناصر المحدودة بعناصر خطية. تتكرر حلقة التحسين بين حل معادلات الحالة والمرافقة، وتحديث متغير المجال الطوري باستخدام طريقة قائمة على التدرج، وإسقاط الحل لتلبية قيد الحجم. يتم تلخيص الخوارزمية على النحو التالي:

1. تهيئة المجال الطوري $\phi^0$
2. حل معادلة الحالة لإيجاد $\mathbf{u}^k$
3. حل المعادلة المساعدة لـ $\mathbf{w}^k$
4. حساب الحساسية $\delta \Pi / \delta \phi$
5. تحديث $\phi^{k+1} = \phi^k - \alpha \nabla_\phi \Pi$
6. إسقاط $\phi^{k+1}$ لتحقيق قيد الحجم
7. التحقق من التقارب؛ إذا لم يتقارب، انتقل إلى الخطوة 2

4.2 مثال على عارضة ناتئة ثنائية الأبعاد

تُستخدم مسألة عارضة ناتئة ثنائية الأبعاد للتحقق من صحة الطريقة. العارضة مثبتة من الطرف الأيسر وتتعرض لحمل هابط على الطرف الأيمن. يتم تقسيم مجال التصميم باستخدام شبكة بحجم 100×50. يتقارب التحسين في حوالي 50 تكرارًا، منتجًا طوبولوجيا تشبه هيكلًا جمالونيًا مع تقليل تركيزات الإجهاد.

5. النتائج والمناقشة

5.1 دراسة الحساسية

تُجرى دراسة حساسية لتحليل تأثير المعلمات الرئيسية: معامل الجزاء $p$ في نموذج المجال الطوري، وسماحية قيد الإجهاد $\epsilon$، ونسبة الحجم $V_f$. تُظهر النتائج أن زيادة $p$ تؤدي إلى واجهات أكثر حدة ولكنها قد تسبب عدم استقرار عددي. يقلل قيد الإجهاد بشكل فعال من الإجهاد الأقصى بنسبة تصل إلى 30% مقارنة بالتصاميم بدون القيد.

5.2 سير عمل الطباعة ثلاثية الأبعاد

يتم تحويل الطوبولوجيا المحسّنة إلى ملف STL وطباعتها باستخدام طابعة ثلاثية الأبعاد بتقنية الترسيب المنصهر (FDM). يشمل سير العمل ما يلي:

• تصدير حل المجال الطوري إلى شبكة
• تنعيم الأسطح البينية
• توليد كود G للطابعة
• الطباعة باستخدام مادة PLA عند درجة حرارة فوهة 200°C

6. التحليل الأصلي

الرؤية الأساسية: تسد هذه الورقة فجوة حرجة في التحسين الطوبولوجي للتصنيع الإضافي من خلال دمج قيود الإجهاد بشكل صارم في إطار المجال الطوري. بينما تركز معظم الطرق الحالية فقط على تقليل الامتثال، فإن تضمين قيود الإجهاد يعالج بشكل مباشر آليات الفشل السائدة في الأجزاء المطبوعة ثلاثية الأبعاد، مثل التفكك والكسر تحت الأحمال الحرارية والميكانيكية.

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

Strengths & Flaws: القوة الرئيسية هي الصرامة الرياضية في اشتقاق شروط المثالية، مما يوفر أساسًا متينًا للتوسعات المستقبلية. إن تضمين قيد الإجهاد ذو صلة عملية بالتصنيع الإضافي، كما أشارت الدراسات الحديثة (مثل، Liu et al., 2018, Structural and Multidisciplinary Optimization). ومع ذلك، فإن للورقة عيوبًا ملحوظة: (1) الأمثلة العددية تقتصر على البعد الثنائي، بينما تطبيقات التصنيع الإضافي الحقيقية هي بطبيعتها ثلاثية الأبعاد؛ (2) لا تتم مناقشة التكلفة الحسابية لتحليل الحساسية المساعد، والتي قد تكون باهظة للمشكلات واسعة النطاق؛ (3) قيد الإجهاد عالمي (شكل متكامل)، مما قد لا يلتقط تركيزات الإجهاد المحلية بشكل فعال. بالمقارنة مع عمل Sigmund and Maute (2013, Structural and Multidisciplinary Optimization)، الذي يستخدم نهج SIMP مع قيود إجهاد محلية، تقدم هذه الطريقة خصائص رياضية أفضل ولكنها قد تكون أقل كفاءة للمشكلات على المستوى الصناعي.

رؤى قابلة للتنفيذ: بالنسبة للممارسين، هذه الطريقة مناسبة بشكل أفضل للمشاكل متوسطة وصغيرة الحجم حيث تكون قيود الإجهاد حاسمة، مثل الغرسات الطبية أو دعامات الفضاء الجوي. للتوسع إلى مشاكل أكبر، يجب على المؤلفين النظر في (أ) استخدام تحسين الشبكة التكيفي لتقليل التكلفة الحسابية، (ب) تنفيذ صياغة محلية لقيود الإجهاد (مثل استخدام نهج p-norm)، و (ج) التوسع إلى 3D باستخدام الحوسبة المتوازية. سير العمل من التحسين إلى الطباعة هو مساهمة قيّمة، لكن خطوة التنعيم تحتاج إلى ضبط دقيق لتجنب فقدان الميزات المحسّنة.

7. التفاصيل التقنية

تعتمد الصياغة الرياضية على المعادلات الرئيسية التالية:

معادلة الحالة: $$-\nabla \cdot (\mathbb{C} \varepsilon(\mathbf{u})) = \mathbf{f} \quad \text{in } \Omega$$

تطور المجال الطوري: $$\frac{\partial \phi}{\partial t} = -M \frac{\delta \Pi}{\delta \phi}$$

قيد الإجهاد: $$\sigma_{vm} = \sqrt{\frac{3}{2} \sigma^d : \sigma^d}$$

حيث $\sigma^d$ هو موتر الإجهاد الانحرافي. يستخدم الاستيفاء المادي مخططًا جزائيًا: $\mathbb{C}(\phi) = \phi^p \mathbb{C}_0$، حيث يضمن $p \geq 3$ تصميمًا شبه ثنائي.

8. النتائج التجريبية

ينتج مثال الكابول ثنائي الأبعاد طوبولوجيا بنسبة حجم 40%. يقلل قيد الإجهاد من أقصى إجهاد فون ميسز من 120 ميجا باسكال إلى 85 ميجا باسكال، أي انخفاض بنسبة 29%. تزداد المطاوعة بنسبة 12% فقط، مما يشير إلى مقايضة مواتية. يوضح الشكل 1 (غير معروض) الطوبولوجيا المحسّنة، مُظهرًا هيكلًا شبكيًا واضحًا مع واجهات ناعمة. تكشف دراسة الحساسية أن معامل الجزاء $p=3$ يعطي أفضل توازن بين الواجهات الحادة والاستقرار العددي.

9. دراسة حالة: عارضة ناتئة

إعداد المشكلة: عارضة ناتئة ثنائية الأبعاد بطول 1 متر وارتفاع 0.5 متر مثبتة من الطرف الأيسر. يتم تطبيق حمل نقطي مقداره 1000 نيوتن إلى الأسفل عند الطرف الأيمن. المادة هي PLA بمعامل يونغ $E=3.5$ جيجا باسكال، ونسبة بواسون $\nu=0.35$، وإجهاد الخضوع $\sigma_y=60$ ميجا باسكال.

معلمات التحسين:

• نسبة الحجم: 40%
• معامل الجزاء: $p=3$
• تفاوت قيد الإجهاد: $\epsilon=0.01$
• الشبكة: عناصر رباعية الأضلاع 100×50

النتائج: يحقق التصميم الأمثل مطاوعة مقدارها 0.45 جول وإجهادًا أقصى مقداره 58 ميجا باسكال، مما يفي بقيد الإجهاد. يتكون الهيكل من مسارين رئيسيين للحمل: دعامة قطرية من نقطة الحمل إلى الزاوية العلوية اليسرى، وعضو أفقي على طول الحافة السفلية.

10. التطبيقات المستقبلية

تتمتع الطريقة بإمكانات كبيرة للتطبيقات المستقبلية:

• المواد متعددة المقاييس: توسيع نموذج المجال الطوري لمعالجة المواد ذات التدرج الوظيفي (FGM) ذات الخصائص المتغيرة مكانيًا، مما يتيح تصميمات ذات صلابة وقوة مخصصة.
• الطباعة رباعية الأبعاد: دمج قيود زمنية للمواد ذات الذاكرة الشكلية، مما يسمح بهياكل تتغير شكلها بمرور الوقت.
• التصنيع الإضافي واسع النطاق: توسيع نطاق الخوارزمية لمعالجة المشكلات ثلاثية الأبعاد باستخدام الحوسبة المتوازية وتسريع GPU، مع استهداف تطبيقات في صناعات الطيران والسيارات.
• التحسين متعدد الفيزياء: ربط القيود الحرارية والميكانيكية والسائلة للأجزاء متعددة الوظائف، مثل المبادلات الحرارية أو الآليات الممتثلة.

11. المراجع

1. Auricchio, F., et al. (2019). Structural multiscale topology optimization with stress constraint for additive manufacturing. arXiv preprint arXiv:1907.06355.
2. ليو، ج. وآخرون (2018). تحسين الطوبولوجيا مع قيود الإجهاد للتصنيع الإضافي. Structural and Multidisciplinary Optimization، 58(6)، 2485-2500.
3. Sigmund, O., & Maute, K. (2013). Topology optimization approaches. Structural and Multidisciplinary Optimization، 48(6)، 1031-1055.
4. Bendsøe, M. P., & Sigmund, O. (2003). تحسين الطوبولوجيا: النظرية والطرق والتطبيقات. سبرينغر.
5. Deaton, J. D., & Grandhi, R. V. (2014). A survey of structural and multidisciplinary continuum topology optimization. Structural and Multidisciplinary Optimization, 49(1), 1-38.