Bridging perturbation and variational approaches in brittle fracture

ملخص تقني: الجسر بين نهج الاضطراب والنهج التبايني في الكسر الهش

بيان المشكلة
تعد نمذجة انتشار الشقوق ثلاثية الأبعاد (3D) في المواد الصلبة غير المتجانسة والهشة تماماً عملية مكلفة حاسوبياً بسبب تركيزات الإجهاد المنفردة عند جبهة الشق والطبيعة متعددة المقاييس للكسر في الأوساط المضطربة. تواجه الطرق الحالية مقايضات كبيرة:

• طرق المجال الطوري (Phase-field methods): (بناءً على فرانسفورت وماريغو، 1998) تتميز بالمتانة في الانتشار ثلاثي الأبعاد في الأوساط المتقلبة، لكنها تتطلب شبكات حجمية مكلفة وتدخل طول مقياس ضرر منتشر ($\ell$) يمكن أن يتفاعل مع أحجام عدم التجانس ويغير فيزياءها.
• نهج الاضطراب (Perturbation approaches): (بناءً على رايس، 1985) تحافظ على شقوق حادة وتوفر كفاءة حاسوبية عالية من خلال إنشاء شبكة لجبهة الشق فقط. ومع ذلك، فهي تعتمد تقليدياً على تنظيم لزوجي "ad hoc" لمعيار غريفيث، مما قد يؤدي إلى عيوب عددية تؤثر على ديناميكيات الكسر. علاوة على ذلك، تركز معظم دراسات الاضطراب على الشقوق شبه اللانهائية في النمط الأول (Mode I)، مما يترك تأثيرات الحجم المحدود، وكثافة الاضطراب، والتحميل المختلط (I+II+III) على الشقوق ثلاثية الأبعاد دون فهم كافٍ.

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

1. الصيغة الاضطرابية التباينية:

   • يتم تقليل الطاقة الكلية $\Pi_{tot}$ مع مراعاة قيد عدم ال reversibility (الشق لا يمكن أن يلتئم).
   • الطاقة الكامنة ($\Pi_{pot}$): بدلاً من توسيع معدل إطلاق الطاقة ($G$) كما في طرق الاضطراب التقليدية، يشتق المؤلفون توسعاً تقاربياً لـ الطاقة الكامنة نفسها لجبهة شق شبه دائرية مضطربة من نصف قطر مرجعي $a_0$. يستخدم هذا الاشتقاق نظرية دالة الوزن لـ "بويكنر-رايس" ونظرية "غاو" (1988) لعوامل شدة الإجهاد (SIFs) في الأنماط المختلطة.
   • الطاقة المبددة ($\Pi_{dis}$): تُحسب عن طريق تكامل مجال طاقة الكسر غير المتجانس $G_c$ عبر السطح المشقوق.
   • التقييم الفعال: يتم تقييم الطاقة الكامنة ومشتقاتها باستخدام تحويلات فوريه السريعة (FFT)، مستفيدين من الطبيعة غير الموضعية للتفاعلات المرنة (الممثلة بمؤثر لابلاس الكسري وتحويل هيلبرت).

2. التنفيذ العددي:

   • تمت صياغة المشكلة كمسألة تقليل غير محدبة ذات قيود صندوقية (box-constrained minimization).
   • المحلل (Solver): تم استخدام خوارزمية نيوتن-كونيغ (Newton-conjugate gradient) خالية من المصفوفات مع منطقة ثقة (trust region).
   • القيود:
     • القيود الصندوقية تفرض عدم عكسية الشق ($a_i \ge a_i^{prev}$).
     • قيود منطقة الثقة (حيث يتم تحديد نصف القطر بواسطة طول الارتباط لعدم التجانس $d$) تمنع المحلل من القفز فوق حواجز الطاقة، مما يضمن اختيار حالات التوازن شبه المستقرة ذات المعنى الفيزيائي بدلاً من القيم الدنيا العالمية التعسفية.
   • التحسين المسبق (Preconditioning): يتم تطبيق مُحسن مسبق قائم على الفيزياء، مشتق من مصفوفة "هيسيان" العكسية لمسألة متجانسة مكافئة، لتسريع التقارب، خاصة في حالات تباين عدم التجانس المنخفض.
   • البرمجيات: تم التنفيذ باستخدام لغة بايثون مع مكتبات petsc4py و JAX للتمايز التلقائي والترجمة الفورية (JIT).

المساهمات الرئيسية

• الاشتقاق النظري: اشتقاق صارم للتوسع التقاربي للطاقة الكامنة لشق شبه دائري تحت تحميل مختلط (I+II+III)، مما يوفر صيغة مغلقة تربط المبادئ التباينية بنظرية الاضطراب.
• الإطار الخوارزمي: محلل قوي خالٍ من المصفوفات يفرض القيود الفيزيائية (عدم العكسية، حواجز الطاقة، التفاعلات بعيدة المدى) مع الحفاظ على الكفاءة الحسابية (تعقيد $O(N \log N)$ لكل تكرار).
• التحقق: تم التحقق من التنفيذ مقابل حلول تحليلية مشتقة حديثاً للشقوق القصية في الأوساط المتجانسة وغير المتجانسة، مما أظهر دقة في التنبؤ بتشوهات الجبهة وتطور أنماط فوريه.

النتائج
أجرى المؤلفون 116,000 محاكاة واسعة النطاق لانتشار الشقوق الشدّية (النمط I) والقصية (النمط II+III) في أوساط مضطربة لقياس تأثير آثار الحجم المحدود، وكثافة الاضطراب، واختلاط الأنماط.

1. الانتقال إلى التقطع (Intermittency): تعيد عمليات المحاكاة إنتاج الانتقال من النمو المستمر السلس إلى النمو المتقطع (تتبادل مراحل التثبيت والانهيارات الجليدية/avalanches). يتم التحكم في هذا الانتقال من خلال النسبة بين محيط الشق وطول لاركن، والذي يعتمد على كثافة الاضطراب ومقياس عدم التجانس.
2. تأثيرات اختلاط الأنماط (Mode Mixity):
   • لخلط الأنماط تأثير محدود على بداية التقطع.
   • ومع ذلك، في التحميل المختلط (النمط II+III) مع نسبة بواسون غير صفرية ($\nu \neq 0$)، تطور جبهة الشق شكلاً شبه بيضاوي بسبب اختلاف الصلابة بين الأنماط، حيث تتقدم قطاعات النمط III بشكل أبطأ من قطاعات النمط II.
3. تبديد الطاقة والطاقة المشعة: تكشف عمليات المحاكاة أن الديناميكيات المتقطعة تؤدي إلى تقلبات متعددة المقاييس في تبديد الطاقة والطاقة الكامنة. يعمل الفرق بين انخفاضات الطاقة الكامنة وارتفاعات الطاقة المبددة كوكيل للطاقة المشعة، مما يظهر أحداثاً تشبه انفجارات إطلاق العزم الزلزالي.
4. التقوية/الإضعاف المعتمد على الحجم:
   • الشقوق الصغيرة: تسبب عدم التجانس "إضعافاً" ظاهرياً ($\Delta S > 0$) حيث تنحني جبهة الشق نحو المناطق ذات المتانة المنخفضة.
   • الشقوق الكبيرة: يحدث انتقال إلى "التقوية" ($\Delta S < 0$) مع دخول الشق في نظام تثبيت قوي، حيث تسبب عدم الاستقرارات جعل الجبهة تختار تفضيلياً المناطق ذات المتانة العالية.
   • يزداد حجم هذا التأثير مع مربع كثافة الاضطراب ($\sigma/G_c^0$).

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

يؤكد المؤلفون أن نهجهم يسمح بالتقييم المباشر لاستقرار الشقوق وقياس آثار الحجم المحدود واختلاط الأنماط، وهي أمور لم تكن مفهومة تماماً سابقاً في الكسور ثلاثية الأبعاد غير المتجانسة. تسلط النت