A Shifted Cohesive-Zone Method for Non-Interface-Fitted Meshes with Applications to Crystal Plasticity

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

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

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

• محددات الشبكات المتوافقة مع الواجهة: إن توليد شبكات عالية الجودة من نوع رباعي الأضلاع أو سداسي الأوجه تتوافق مع البنى المجهرية ثلاثية الأبعاد المعقدة هو أمر مكلف حاسوبياً، وغالباً ما يؤدي إلى عناصر مشوهة (جاكوبيان سالب) أو يتطلب أدوات تقسيم شبكي معقدة.
• محددات الطرق الحالية غير المتوافقة مع الواجهة: الطرق مثل XFEM أو طرق العناصر المحدودة المغمورة (IFEM) تتطلب غالباً إثراء فضاء التقريب محلياً، مما يؤدي إلى فضاءات دالية غير منتظمة وتعقيد في التنفيذ.
• الفجوة: هناك حاجة إلى طريقة تسمح باستخدام شبكات قياسية غير متوافقة (خاصة رباعية الأضلاع/سداسية الأوجه) لمشكلات الواجهة التماسكية مع الحفاظ على الدقة وتجنب أعباء إنشاء الشبكة.

2. المنهجية: طريقة المنطقة التماسكية المزاحة (SCZM)

قام المؤلفون بتوسيع طريقة الحدود المزاحة (SBM) لتشمل صياغات المنطقة التماسكية. الفلسفة الجوهرية هي فرض شروط الواجهة على واجهة بديلة ($\Gamma_h$) محددة بواسطة أوجه الشبكة، بدلاً من الواجهة الفيزيائية الحقيقية ($\Gamma$).

المكونات النظرية الرئيسية:

• بناء الواجهة البديلة: يتم تقريب الواجهة الحقيقية بواسطة مجموعة من أوجه الشبكة ($\Gamma_h$) التي تتقاطع مع الواجهة الحقيقية.
• الصياغة الضعيفة المزاحة:
  • يتم فرض قانون الشد-الانفصال على $\Gamma_h$ بدلاً من $\Gamma$.
  • التصحيح الهندسي: يتم استخدام توسع نقل السطح من الدرجة الأولى لربط تكامل الواجهة الحقيقية بالواجهة البديلة. يتضمن ذلك عامل تصحيح المساحة الهندسية ($|n \cdot n_h|$) لمراعاة الزاوية بين الناظم الحقيقي ($n$) وناظم الوجه ($n_h$).
  • التصحيح الاتجاهي: إضافة حاسمة لـ SBM القياسية. نظرًا لأن ناظم الوجه يختلف عن الناظم الحقيقي، ينشأ مكون شد مماس؛ لذا اشتق المؤلفون حد تصحيح اتجاهي ($\sigma \tau_h$) لمعالجة هذا عدم المحاذاة، مما يضمن اتساق متجه الشد.
• تقريب قفزة الإزاحة: يتم تقريب قفزة الإزاحة ($\llbracket u \rrbracket$) المطلوبة للقانون التماسكِي عند الموقع البديل باستخدام توسع تايلور من الدرجة الأولى: $\llbracket u \rrbracket|_{\Pi(x)} \approx \llbracket u \rrbracket + \llbracket \nabla u \rrbracket \cdot d$، حيث $d$ هو متجه المسافة إلى الواجهة الحقيقية.
• النماذج التكوينية: يدعم الإطار النماذج المعتمدة على التاريخ، وتحديداً اللدونة البلورية (اللدونة اللزوجية أحادية البلورة ذات الانفعال الكبير مع تصلد Voce) المقترنة بمختلف قوانين الشد-الانفصال (الأسية، وثنائية الطور، والنماذج ثلاثية الأبعاد المزدوجة).

التنفيذ والخوارزميات:

• البرمجيات: تم التنفيذ في إطار العمل مفتوح المصدر MOOSE، مقترناً بـ NEML2 للنمذجة التكوينية.
• تصنيف النقاط: تم تطوير خوارزمية تصنيف نقاط محسنة بـ PCA لتحديد الحبيبة التي ينتمي إليها كل عنصر. بدلاً من إسقاط الأشعة الثابت، يتم إسقاط الأشعة على طول المكون الرئيسي لأقل تباين في هندسة الحدود، مما يقلل بشكل كبير من عمليات التحقق من التقاطع.
• بناء البديل: يتم تعيين معرف حبيبة واحد لكل عنصر بناءً على الكسر الحجمي المهيمن.
• التحقق: يقوم إجراء نقل الحل بإسقاط نتائج SCZM على شبكة مرجعية متوافقة مع الواجهة (IFM) للتصور والمقارنة.

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

1. أول تطبيق على اللدونة البلورية: يوضح SCZM للنماذج التكوينية المعتمدة على التاريخ، وهي خطوة هامة تتجاوز تطبيقات المرونة الخطية السابقة.
2. صياغة ضعيفة مبسطة: اشتق صياغة ضعيفة مختزلة على الواجهة البديلة تبسط تنفيذ البواقي غير الخطية ومصفوفات المماس المتسقة.
3. التصحيح الاتجاهي: قدم حداً تصحيحياً ضرورياً للشبكات غير المتوافقة مع الواجهة للتعامل مع عدم محاذاة الناظم، مما يمنع تذبذبات الإجهاد والأخطاء.
4. التصنيف المدرك للهندسة: طور خوارزمية تعتمد على PCA تسرع بناء النطاقات البديلة، مما يجعل الطريقة قابلة للتوسع للنماذج ثلاثية الأبعاد المعقدة.
5. التوفر مفتوح المصدر: يوفر أول تنفيذ ضمن إطار عمل MOOSE، مما يسمح بالتكامل مع وحدات الفيزياء الموجودة.

4. النتائج العددية والتحقق

تم التحقق من صحة الطريقة من خلال سلسلة تصاعدية من الاختبارات المرجعية:

• الحل المصنع (MMS):
  • أظهر تقارباً من الدرجة الأولى للحالات العامة (المحدودة بالتقريب الهندسي للواجهة).
  • أظهر تقارباً من الدرجة الثانية للحالات المبسطة حيث تتلاشى حدود تصحيح المجال.
  • الطرق القياسية غير المتوافقة مع الواجهة (بدون إزاحة) فشلت في التقارب.
• اختبارات الواجهة الواحدة:
  • المرونة: طابقت نتائج SCZM الحلول المرجعية المتوافقة مع الواجهة (IFM) لقوى التفاعل، وطاقة التحرر، وحقول الإجهاد.
  • اللدونة البلورية: طابقت SCZM مع التصحيح الاتجاهي الحلول المرجعية IFM بشكل وثيق لقوى التفاعل وتوزيعات الإجهاد. أظهر النهج "الساذج" (بدون إزاحة أو تصحيح) أخطاءً كبيرة.
• النماذج المجهرية متعددة البلورات (2D و 3D):
  • ثنائية الأبعاد (2D): تقاربت SCZM على شبكات رباعية الأضلاع غير المتوافقة مع حل IFM مع زيادة دقة الشبكة. لقد رصدت بدقة منحنيات القوة-الإزاحة العالمية، وتطور الضرر، وتركيزات الإجهاد المحلية (مثل نقاط التقاء الثلاثي).
  • ثلاثية الأبعاد (3D): طُبقت على نموذج RVE ثلاثي الأبعاد مكون من 15 حبيبة باستخدام عناصر سداسية الأوجه. أعادت SCZM إنتاج أنماط التشوه، وقوى التفاعل، وتطور الضرر لحل IFM بدقة عالية، رغم التعقيد في هندسة الحبيبات.
• الأداء: أظهرت خوارزمية التصنيف المحسنة بـ PCA زمن تشغيل يقترب من الثبات بغض النظر عن تعقيد الحدود، مما وفر تسريعاً كبيراً مقارنة بطريقة إسقاط الأشعة العشوائية.

5. الأهمية والخاتمة

• الكفاءة: تلغي SCZM الحاجة إلى إنشاء شبكات تستغرق وقتاً طويلاً ومتوافقة مع الواجهة، خاصة في البنى المجهرية ثلاثية الأبعاد المعقدة حيث تُفضل العناصر سداسية الأوجه من أجل الدقة ولكن يصعب إنشاؤها.
• الدقة: من خلال دمج التصحيحات الهندسية والاتجاهية، تحقق SCZM دقة تضاهي الطرق المتوافقة مع الواجهة دون عبء إنشاء الشبكة.
• تعدد الاستخدامات: تتميز الطريقة بالمتانة عبر أنواع مختلفة من العناصر (مثلثية، رباعية الأضلاع، سداسية الأوجه) والنماذج المادية (من المرونة الخطية إلى اللدونة البلورية المعقدة).
• العمل المستقبلي: يقترح المؤلفون توسيع الطريقة لتشمل تصحيحات فجوة هندسية من درجات أعلى (مستوحاة من Gap-SBM) والاقتران مع تحسين الشبكة التكيفي (القائم على octree) لمزيد من الكفاءة.

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