A unified variational framework for phase-field fracture and third-medium contact in finite deformation hyperelasticity

إليك ملخص تقني مفصل للورقة البحثية بعنوان: "إطار عمل متغير موحد لكسر حقل الطور وتلامس الوسط الثالث في المرونة الفائقة ذات التشوه المحدود."

1. بيان المشكلة

تتناول الورقة البحثية التحدي الحسابي المتمثل في محاكاة ظواهر التلامس والكسر المقترنة ضمن إطار مرونة فائقة ذات تشوه محدود. تواجه مناهج ميكانيكا الحوسبة التقليدية عقبات كبيرة عندما تتفاعل هاتان الظاهرتان:

• تعقيد التلامس: تتطلب خوارزميات التلامس الكلاسيكية اكتشافاً صريحاً لأزواج التلامس، وتحديد أسطح "السيد والعبد" (master-slave)، وفرض قيود عدم الاختراق (عبر طرق الجزاء، أو مضاعفات لاغرانج، أو لاغرانج المعززة). وتصبح هذه العمليات مكلفة حسابياً ومعقدة خوارزمياً عندما تتطور الأسطح ديناميكياً (مثل انغلاق وجوه الشقوق أو تفاعل الحطام).
• تعقيد الكسر: يلغي نموذج كسر حقل الطور (PFF) الحاجة إلى تتبع الشقوق الصريح عن طريق تنظيم الشقوق الحادة إلى مجالات ضرر منتشرة. ومع ذلك، فإن نماذج (PFF) القياسية لا تتعامل بطبيعتها مع تلامس وجوه الشقوق المتكونة حديثاً أو التفاعل بين مدخل صلب وركيزة متكسرة دون طبقات خوارجزية إضافية.
• الفجوة: هناك نقص في وجود إطار عمل متغير موحد يربط بشكل طبيعي بين تنظيم الشقوق (PFF) وتنظيم واجهات التلامس (تلامس الوسط الثالث - TMC) لمحاكاة الكسر الناجم عن التلامس دون تتبع صريح للواجهة.

2. المنهجية

يقترح المؤلفون إطار عمل متغير موحد يقوم على مبدأ أن كل من طوبولوجيا الشقوق وواجهات التلامس يمكن تنظيمهما باستخدام أوصاف حقول مستمرة. تدمج المنهجية ثلاثة مكونات متميزة في دالة طاقة كامنة كلية واحدة ($\Pi$):

أ. تقسيم النطاق

يتم تقسيم النطاق الحسابي ($\Omega_0$) إلى ثلاثة نطاقات فرعية غير متداخلة:

1. الركيزة ($\Omega_1$): الجسم القابل للتشوه حيث يحدث الكسر. يتم نمذجته باستخدام كسر حقل الطور (PFF) مع متغير ضرر $d \in [0, 1]$.
2. المدخل ($\Omega_2$): جسم مرن صلب أو شديد الصلابة لا يتعرض للكسر.
3. الوسط الثالث ($\Omega_3$): مادة وهمية مرنة تملأ الفجوة بين $\Omega_1$ و $\Omega_2$. وهي تعمل كواجهة تلامس.

ب. دالة الطاقة الموحدة

يتم تقليل الطاقة الكلية بالنظر إلى الإزاحة ($u$)، والضرر ($d$)، والحقول المساعدة ($p, q$):
$$\Pi(u, d, p, q) = \Pi_{int} + \Pi_{frac} - \Pi_{ext}$$

• طاقة الركيزة ($\Omega_1$): تستخدم تفكيك طاقة الانفعال الطيفي لفصل طاقات الشد ($\Psi^+$) والضغط ($\Psi^-$). يتم إضعاف الجزء الشدي فقط بواسطة دالة الضرر $g(d)$، مما يمنع التئام الشقوق تحت تأثير الضغط.
• طاقة الكسر ($\Omega_1$): تستخدم دالة كثافة الشق AT2 (مشهد طاقة سلس) لتنظيم سطح الشق الحاد إلى حزمة منتشرة بعرض $2\ell$.
• طاقة الوسط الثالث ($\Omega_3$): تُنمذج كمادة مرنة فائقة من نوع نيو-هوكيان (Neo-Hookean) ذات صلابة منخفضة للغاية (معلمة القياس $\gamma \ll 1$). مع اقتراب الأجسام، ينضغط الوسط ($J \to 0$)، مما ينقل قوى التلامس بشكل طبيعي دون اكتشاف صريح.
• تنظيم الحقل المساعد: لمنع تشوه الشبكة وانقلاب العناصر أثناء الضغط الشديد للوسط الثالث، يتم إدخال حقول مساعدة ($p$ للدوران المحلي، و $q للتشوه الحجمي$). تقوم هذه الحقول بتنظيم الكميات الكينماتيكية، مما يضمن الاستقرار العددي.

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

• التقسيم (Discretization): تُستخدم عناصر لاغرانج المستمرة لجميع الحقول.
• استراتيجية الحل: يتم استخدام مخطط متدرج (تقليل بديل) بسبب عدم تقعر الطاقة الكلية.
  1. خطوة الميكانيكا/التلامس: حل الإزاحة ($u$) والحقول المساعدة ($p, q$) بشكل متكامل مع تثبيت الضرر ($d$).
  2. خطوة الضرر: حل الضرر ($d$) مع تثبيت الكينماتيكا، مع مراعاة قيود عدم الارتداد ($\dot{d} \geq 0$).
• البرمجيات: تم التنفيذ باستخدام FEniCS مع مكتبة multiphenics للأشكال المتغيرة ذات البنية الكتلية.

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

1. صياغة متغيرة موحدة: تؤسس الورقة لأول إطار عمل يربط بسلاسة بين PFF و TMC ضمن دالة طاقة واحدة، مما يلغي الحاجة إلى خوارزميات اكتشاف التلامس الصريحة أو تتبع الشقوق.
2. تنظيم الحقل المساعد للتلامس: وسعت المنهجية طريقة TMC باستخدام الحقول المساعدة للتعامل مع التشوه الشديد للشبكة في الوسط الثالث، مما يسمح بمحاكاة قوية للتلامس في التشوهات المحدودة حيث $J \to 0$.
3. التعامل الطبيعي مع الطوبولوجيا المتطورة: يعالج الإطار تلقائياً تلامس وجوه الشقوق المتكونة حديثاً وتطور منطقة التلامس بين المدخل والركيزة، حيث يتكيف الوسط الثالث مع الهندسة المتغيرة دون تدخل خوارزمي.
4. التفكيك الطيفي في التشوه المحدود: نجحت الورقة في دمج تفكيك طاقة الانفعال الطيفي (لمنع التئام الشقوق) مع المرونة الفائقة للتشوه المحدود في منطقة التلامس.

4. النتائج والتحقق

تم التحقق من صحة الإطار من خلال أربع أمثلة عددية:

• اختبار صندوق C ثنائي الأبعاد (التحقق من TMC): أظهر أن الوسط الثالث ينقل قوى التلامس بنجاح ويحافظ على جودة الشبكة تحت الضغط الشديد. وقد تبين أن الحقول المساعدة ضرورية؛ فبدونها، كان الوسط يعمل إما كحاجز صلب أو يتشوه بشكل مفرط.
• اختبار الحافة الواحدة المنحوتة (SEN) ثنائي الأبعاد (التحقق من PFF): أعاد إنتاج مسارات كسر النمط الأول (شق مستقيم) والنمط الثاني (شق ملتوي/منحنٍ)، مما أثبت صحة التفكيك الطيفي وتطور الضرر.
• اختبار الانحناء ثلاثي النقاط: حاكت عملية ضغط مدخل صلب على عارضة بها شق مسبق. التقط النموذج الانتقال من تلامس نقطي إلى منطقة تلامس محدودة مع تكيف المدخل مع الركيزة، مما أدى إلى فتح الشق المسبق ونموه للأعلى.
• اختبار قرص برازيلي ثلاثي الأبعاد:
  • أعاد إنتاج شق الانقسام الرأسي المميز على طول قطر التحميل.
  • النتيجة الرئيسية: التقطت المحاكاة بشكل طبيعي مناطق كسر ثانوية من نوع السحق (crushing-type) بالقرب من مناطق التلامس (منصات التحميل). هذه الظاهرة، الناتجة عن تركيز الإجهاد من النمو التدريجي لمنطقة التلامس، تتوافق مع الملاحظات التجريبية ولكنها مستحيلة الإعادة باستخدام النماذج المبسطة التي تفرض هندسات تلامس ثابتة (مثل أحمال الخطوط).

5. الأهمية والآفاق المستقبلية

• القدرة التنبؤية: يسمح هذا الإطار بالمحاكاة التنبؤية لتفاعلات التلامس والكسر المعقدة (مثل الفشل الناجم عن التغلغل/النقرة، أو التفتت) دون افتراضات مسبقة حول هندسة التلامس أو مسارات الشقوق.
• الدقة التجريبية: من خلال التقاط مناطق السحق الثانوية في اختبار القرص البرازيلي، يسد النموذج الفجوة بين نماذج التحميل النظرية المبسطة والواقع التجريبي المعقد.
• القيود:
  • لا يتم تحقيق تلامس بصفر فجوة تماماً؛ حيث توجد فجوة متبقية صغيرة بسبب السمك المحدود للوسط الثالث (يتحكم فيه $\gamma$).
  • التكلفة الحسابية عالية بسبب درجات الحرية الإضافية (الوسط الثالث + الحقول المساعدة) ومخطط الحل المتدرج.
• العمل المستقبلي: يخطط المؤلفون لتوسيع الإطار ليشمل الاحتكاك، الالتصاق، القصور الذاتي (الاصطدام)، المواد متباينة الخواص، والاقتران الحراري/الكيميائي. أحد التطبيقات المستهدفة هو تحليل قضبان الوقود النووي (التفاعل الميكانيكي بين الحبيبات والغلاف).

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