An extended deep energy method for thermo-mechanical crack propagation
تقدم هذه الورقة طريقة طاقة عميقة ممتدة تستخدم شبكتين عصبيتين وتمثيلاً خطياً متعدد الأضلاع حاداً لحل مشكلات انتشار الشقوق الحرارية الميكانيكية المقترنة، محققةً دقة عالية في استخراج عامل شدة الإجهاد وتنبؤ مسار الشق عبر ظروف تحميل متنوعة دون الحاجة إلى طول تنظيم.
المؤلفون الأصليون:Han Zhang, Mehrisadat Makki Alamdari, Babak Shahbodagh, Mohammad Vahab, Cosmin Anitescu, Timon Rabczuk, Elena Atroshchenko
في عالم الهندسة، غالبًا ما تفشل المواد ليس فقط بسبب شدها بقوة، ولكن بسبب تسخينها أو تبريدها بشكل غير متساوٍ. فعندما يُغمر طبق من السيراميك في ماء بارد، أو يواجه مكون معدني في محرك نفاث حرارة شديدة، تؤدي فروق درجات الحرارة إلى نشوء إجهادات داخلية يمكن أن تكون قوية بما يكفي لكسر المادة. ونادرًا ما تسير هذه الشقوق في خطوط مستقيمة؛ بل تنحني وتتفرع وتلتوي في أنماط يصعب التنبؤ بها. وبالنسبة للمهندسين الذين يصممون كل شيء، بدءًا من المفاعلات النووية وصولًا إلى المركبات الفضائية، فإن فهم كيفية نمو الشق بدقة تحت تأثير هذه القوى الحرارية والميكانيكية مجتمعة هو مسألة تتعلق بالسلامة والبقاء. وغالبًا ما تواجه الطرق الحاسوبية التقليدية للتنبؤ بهذا السلوك صعوبة، لأنها يجب أن تعيد تشكيل شبكتها الرقمية باستمرار لتتبع المسار المتعرج لشق متحرك، وهي عملية مكلفة حاسوبيًا وعرضة للخطأ.
لقد طور فريق من الباحثين طريقة جديدة لمحاكاة هذا السلوك المعقد باستخدام الذكاء الاصطناعي، وتحديدًا تقنية تعامل المشكلة كمسألة تقليل للطاقة. وبدلاً من إجبار الكمبيوتر على حل سلسلة من المعادلات الصعبة خطوة بخطوة، تطلب الطريقة الجديدة من شبكة عصبية إيجاد حالة الطاقة الأدنى للنظام بأكمله دفعة واحدة. وقد بنى الباحثون نظامًا تعمل فيه شبكتان عصبيتان منفصلتان بالتوازي: إحداهما تتنبأ بكيفية تدفق الحرارة عبر المادة، والأخرى تتنبأ بكيفية تشوه المادة تحت الإجهاد. ومن الأهمية بمكان أن كلتا الشبكتين تُعطيان "خريطة" خاصة لموقع الشق. وتسمح هذه الخريطة للشبكات بفهم أن المادة منقسمة إلى جزأين، مما يمكنها من تمثيل القفزة المفاجئة في درجة الحرارة وانفصال أوجه المادة دون الحاجة إلى إعادة رسم الشبكة الحاسوبية في كل مرة يتحرك فيها الشق.
اختبر الباحثون هذا النهج في عدة سيناريوهات صعبة، بما في ذلك لوح معدني به ثلمة تم تسخينه ثم سحبه بعيدًا، وعينة على شكل صليب تعرضت لكل من الحرارة والقوة. وفي كل حالة، نجحت الطريقة في التنبؤ بالمسار الذي سيسلكه الشق. وأظهرت عمليات المحاكاة أن الشق سينمو في الاتجاه الذي يحرر أكبر قدر من الطاقة، وهي قاعدة تسري في العالم الحقيقي. كما تحقق الفريق من أن طريقتهم يمكنها حساب شدة الإجهاد عند طرف الشق بدقة، وهي قيمة حرجة تحدد ما إذا كان الشق سيتوقف أم سيستمر في النمو. ومن خلال مقارنة نتائجهم بالحلول المعروفة والنماذج الحاسوبية المتقدمة الأخرى، وجدوا أن تنبؤاتهم تطابقت بدقة عالية، وغالبًا في غضون جزء من المائة في المئة.
أحد الجوانب الأكثر أهمية لهذا العمل هو كيفية تعامله مع الحافة الحادة للشق. فالعديد من طرق المحاكاة الحديثة تقوم بتنعيم الشق وتحويله إلى منطقة ضبابية لتسهيل الرياضيات، ولكن هذا يمكن أن يخفي تفاصيل مهمة حول كيفية سلوك طرف الشق. أما الطريقة الجديدة فتحافظ على الشق كخط حاد ومتميز طوال المحاكاة بأكملها. وهذا يسمح للباحثين بقراءة قيم الإجهاد مباشرة من الطرف دون الحاجة إلى التخمين أو التقدير من منطقة ضبابية. كما تأخذ الطة في الاعتبار حقيقة أن مقاومة المادة للتشقق يمكن أن تتغير اعتمادًا على مدى سخونة طرف الشق. ومن خلال حل تدفق الحرارة والإجهاد الميكانيكي معًا، يلتقط النظام حلقة التغذية الراجعة حيث تسبب الحرارة الإجهاد، ويغير الشق الناتج طريقة تدفق الحرارة.
تؤكد الدراسة أن نهج الشبكة العصبية القائم على الطاقة هو أداة قوية لفهم الكسر الحراري الميكانيكي. فهو يتجنب الحاجة إلى شبكات حاسوبية معقدة ومتغيرة باستمرار، ويوفر طريقة مباشرة ودقيقة لرؤية كيفية سلوك الشقوق في ظل الظروف القاسية الموجودة في الهندسة الواقعية. وقد أظهر الباحثون أن الطريقة تعمل لكل من المواد البسيطة وتلك ذات الخصائص المتغيرة، مثل المواد المتدرجة وظيفيًا المصممة لتحمل البيئات الحرارية القاسية. ومن خلال النجاح في التنبؤ بمسارات الشقوق في سيناريوهات الشد، والقص، والتحميل المشترك، يقدم هذا العمل طريقة موثوقة جديدة للمحاكاة، مما قد يساعد المهندسين في تصميم هياكل أكثر أمانًا يمكنها تحمل التهديدات المزدوجة للحرارة والقوة.
ملخص تقني: طريقة الطاقة العميقة الممتدة للانتشار الحراري-الميكانيكي للشقوق
بيان المشكلة يتضمن الكسر الحراري-الميكانيكي اقتران التوصيل الحراري العابر على نطاق مشقوق مع انتشار الشق المدفوع بالإجهادات الحرارية والميكانيكية المتطورة. وبينما تم توسيع نطاق حلول الطاقة العصبية مؤخرًا لتمثيل الشقوق الحادة في الظروف المتساوية الحرارة، إلا أنها لم تعالج بعد المعالجة المتزامنة لتوصيل الحرارة على نطاق مشقوق وانتشار الشق تحت الإجهادات الحرارية الناتجة. تعتمد النهج الحالية غالبًا على طرق "حقل الطور" (phase-field)، التي تنظم الشق إلى حزمة ضرر منتشرة، أو الطرق التقليدية القائمة على الشبكات (مثل XFEM أو مناطق التماسك) التي تتطلب تتبعًا هندسيًا صريحًا وإعادة بناء للشبكة. توجد فجوة كبيرة في الحلول القائمة على الشبكات العصبية التي يمكنها التعامل مع شق حاد، وحل حقل درجة الحرارة على النطاق المشقوق، ونشر الشق بناءً على الحقول المدربة دون الحاجة إلى تنظيم طول أو تحديثات صريحة للشبكة.
المنهجية يقترح المؤلفون طريقة الطاقة العميقة الممتدة (EDM) التي تحافظ على الشق كخط متعدد الأضلاع حاد طوال المحاكاة. تدمج المنهجية ثلاثة مكونات أساسية:
تضمين عدم الاستمرارية: يتم ترميز هندسة الشق عبر دالة تضمين قياسية، ρ(x)، تكون غير مستمرة عبر أوجه الشق ومستمرة في مكان آخر. تُضاف هذه الدالة كمدخل لشبكتين عصبيتين منفصلتين: واحدة لدرجة الحرارة (T) والأخرى للإزاحة (u). يسمح هذا لشبكة عالمية واحدة بتمثيل قفزات الإزاحة وقيم درجات الحرارة المتميزة على أوجه الشق المتقابلة دون تقسيم النطاق.
تقليل الطاقة المتدرج: تسير العملية وفق تسلسل متدرج عند كل خطوة تحميل:
الحل الحراري: يتم الحصول على حقل درجة الحرارة عن طريق تقليل دالة تراكمية لتوصيل الحرارة العابر (باستخدام مخطط أويلر الخلفي). تُعامل أوجه الشق كحدود أديباتية (كظومة)، وهو شرط يتم استيفاؤه طبيعيًا بواسطة الصيغة الضعيفة لدالة الطاقة، مضافًا إليها حد اتساق.
الحل الميكانيكي: يتم الحصول على حقل الإزاحة عن طريق تقليل طاقة الوضع المرنة-الحرارية. يتم إثراء شبكة الإزاحة بالقرب من طرف الشق باستخدام توسيع ويليامز (Williams expansion) مع سعات قابلة للتدريب (aI,aII)، على غرار طريقة XFEM، لالتقاط تفرد الإجهاد من رتبة r−1/2.
التكامل: تُقدر الطاقات عبر تكامل مونت كارلو على نقاط موزعة فوق عناصر خلفية مربعة. ومن الأهمية بمكان إعادة رسم هذه النقاط دوريًا أثناء التدريب لمنع الحلول الزائفة (مثل طبقات الانتقال المختبئة بين نقاط ثابتة) وتكثيفها بالقرب من طرف الشق.
انتشار الشق والقوة الدافعة: يُحكم نمو الشق بمعايير ميكانيكا الكسر المرنة الخطية (LEFM) التي يتم تقييمها على الحقول المدربة:
استخراج عوامل شدة الإجهاد (SIFs): يتم استخراج عوامل شدة الإجهاد باستخدام تكامل التفاعل. وتحت التحميل الحراري، يتضمن ذلك حد مساحة ويلسون وييو (Wilson and Yu area term)، وهو أمر ضروري لاستقلال المسار في حقول درجات الحرارة غير المنتظمة.
فحوصات القابلية: لضمان الموثوقية، يتم التحقق من استخراج SIF عن طريق مسح نصف قطر المحيط. لا يتم قبول القراءة إلا إذا كان الانحراف عبر أنصاف الأقطار ضمن قدر معين من التفاوت، مما يؤكد وجود التفرد الصحيح.
معيار النمو: يتقدم الشق عندما يصل معدل تحرر الطاقة في اتجاه أقصى إجهاد حلقي إلى قيمة حرجة تعتمد على درجة الحرارة (Gc(T)). يتم تحديد الاتجاه بواسطة معيار أقصى إجهاد حلقي، ويتم تقييم المقاومة عند درجة حرارة طرف الشق.
المساهمات الرئيسية
حل موحد حراري-ميكانيكي: يقدم هذا العمل حلاً يعتمد على الشبكات العصبية يحل حقل درجة الحرارة على نطاق مشقوق وينشر شقًا حادًا بناءً على معايير الكسر التي يتم تقييمها على تلك الحقول، وهو مزيج لم يسبق تسجيله في الأدبيات.
تمثيل الشق الحاد بدون تنظيم: من خلال الجمع بين تضمين عدم الاستمرارية والإثراء بالقرب من الطرف، تمثل هذه الطريقة شقًا حادًا وتصل إلى قيم SIF مباشرة، متجنبة طول التنظيم وتكلفة "حقول الطور" المنتشرة.
استخراج متين للكسر الحراري: تقدم الورقة بروتوكولًا صارمًا لاستخراج SIF تحت الحمل الحراري، وتحديدًا من خلال دمج حد مساحة ويلسون وييو وفرض فحص قابلية المسح الشعاعي لتصفية القراءات غير الموثوقة.
التعامل مع الهندسات المعقدة: يتم بناء دالة التضمين للتعامل مع الانثناءات والزوايا في الخط متعدد الأضلاع للشق، مما يضمن بقاء عدم الاستمرارية محددة جيدًا حتى مع تطور هندسة الشق.
النتائج تم تقييم الطريقة عبر أربعة اختبارات مرجعية:
شق حافة حراري ثابت: وافق SIF المستخرج قيمة منشورة بنسبة 0.11% على أدق شبكة، مما أثبت صحة إجراء الاستخراج بما في ذلك حد المساحة.
شق مشقوق أحادي الجانب تحت الشد (متجانس): تحت ثلاثة أحمال حرارية (ΔT=−50,0,+50 كلفن)، طابقت مسارات الشق المتوقعة والتحميل الذروي مراجع حقل الطور.
شق مشقوق أحادي الجانب تحت القص (متدرج وظيفيًا): أعادت الطريقة إنتاج مسارات بدء الانتشار بما يتوافق مع مراجع حقل الطور ومع حل مستقل للشق الحاد باستخدام العناصر المحدودة.
عينة صليبية مشقوقة: تحت التحميل الميكانيكي والحراري والمشترك، اتبعت مسارات الشق المتوقعة الحلول المنشورة من ست دراسات مستقلة.
في جميع الحالات، نجحت الطريقة في إعادة إنتاج مسارات الشق وتأثير التحميل الحراري على البدء والحمل الذروي، مع الحفاظ على الشق ككيان هندسي حاد.
الأهمية تزعم الورقة أن هذا النهج يسد الفجوة بين كفاءة الحلول العصبية القائمة على الطاقة والدقة الفيزيائية المطلوبة للكسر الحراري-الميكانيكي للشقوق الحادة. ومن خلال الاحتفاظ بالشق الحاد والوصال المباشر لكميات ميكانيكا الكسر (SIFs، ودرجة حرارة طرف الشق)، تتجنب الطريقة الاسترداد غير المباشر لهذه الكميات المتأصل في مناهج حقل الطور. علاوة على ذلك، فإن الإجراءات العددية المقترحة — وتحديدًا إعادة أخذ عينات نقاط التكامل والتحقق من قابلية مسح نصف القطر لـ SIF — تعالج مشكلات الاستقرار والدقة التي أعاقت تطبيق طرق الطاقة العميقة في مشا المسائل المعقدة والمقترنة للكسر. يوضح العمل أن الشبكات العصبية يمكنها بفعالية حل التوصيل الحراري العابر والالاستقرار المرن شبه الساكن على نطاقات متقطعة متطورة دون الحاجة إلى تتبع صريح للشبكة أو تنظيم.