On the complementary roles of anisotropic crack density and anisotropic crack driving force in phase-field modeling of mixed-mode fracture

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

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

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

وتشمل الفجوات الرئيسية التي تعالجها هذه الدراسة ما يلي:

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

2. المنهجية

قام المؤلفون بتطوير نموذج مجال طوري متسق ديناميكياً حرارياً للكسر متباين الخواص ضمن إطار عمل مرن فائق (hyperelastic) في حدود الانفعال الكبير.

الإطار النظري:

• نوع النموذج: نموذج مجال طوري من نوع AT1 (دالة هندسية خطية $\alpha(d)=d$) مع دالة تدهور تربيعية $g(d)=(1-d)^2$. يضمن هذا وجود عتبة مرونة محدودة، مما يمنع بدء الضرر تحت إجهاد حرج.
• القانون التكويني: طاقة انحراف "نيو-هوكين" المعممة (GNH) مدمجة مع تقسيم حجمي يعتمد على $J$ وطاقة ألياف متباينة الخواص.
• آليات تباين الخواص:
  1. كثافة الشقوق متباينة الخواص ($\gamma_c$): محكومة بموتر هيكلي $A = I + \beta a \otimes a$. وهي تعاقب تدرجات الضرر على طول اتجاه الألياف، مما يزيد فعلياً من مقاومة الكسر عمودياً على الألياف.
  2. قوة الدفع متباينة الخواص ($H$): قوة دفع مختلطة الأنماط ومعيرة معرفة كالتالي:
     H = \left(\frac{\Psi^+_v}{G_{c,I}}\right)^p + \left(\frac{\Psi^+_d}{G_{c,II}}\right)^q + \left(\frac{\Psi^+_{ani}}{G_{c,ani}}\\right)^r
     حيث تمثل $\Psi^+_{ani}$ طاقة الانفعال متباينة الخواص الشدّية المخزنة في الألياف. يلتقط هذا الحد "خلط الأنماط" الناتج عن المادة.
• عدم العكسية: يتم فرضها عبر تقليل مقيد بالحدود ($d_{n} \le d \le 1$) بدلاً من المتغيرات التاريخية، مما يضمن اتساقاً ديناميكياً حرارياً صارماً.
• التنفيذ العددي:
  • المحلل (Solver): مخطط نيوتن/منطقة الثقة المتدرج (التقليل المتناوب) المنفذ في FEniCS.
  • التنظيم (Regularization): تم استخدام تقنيات لضمان التقارب، بما في ذلك صلابة المتبقي، وإزاحات المميز للتحلل الطيفي، وأقواس "ماكولي" الناعمة.

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

1. التحقق من الصحة: لوحة ذات شق حافة ثنائي الأبعاد تحت شد أحادي المحور (مطاط Ecoflex 00-30) للتحقق من البيانات التجريبية (Lu et al., 2021).
2. الدراسة البارامترية 1 (SEN): عينات ذات شق حافة مفرد لعزل آثار تباين كثافة الشقوق ($\beta$) والتباين المرن ($k_1$).
3. الدراسة البارامترية 2 (OHT): عينات شد ذات ثقب مفتوح لاستقصاء نشوء الشقوق وتأثيرات تركيز الإجهاد.

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

1. العزل المنهجي للآليات: أول دراسة تفصل وتقارن بين تأثير كثافة الشقوق متباينة الخواص (مقاومة الكسر) وطاقة الانفعال متباينة الخواص (قوة الدفع) في إطار عمل موحد.
2. اكتشاف التآزر غير الخطي: أظهرت الدراسة أنه عندما تعمل الآليتان معاً، فإن تأثيرهما المشترك على ذروة الحمل وتبديد الطاقة يتجاوز المجموع الخطي لمساهماتهما الفردية.
3. توسع الدور المعتمد على الهندسة: كشفت الدراسة أن دور طاقة الانفعال متباينة الخواص ($\Psi_{ani}$) يعتمد على الهندسة. في العينات ذات الشقوق المسبقة، تؤثر بشكل أساسي على انحراف الشق؛ أما في هندسات تركيز الإجهاد (OHT)، فهي تغير بشكل جوههري الصلابة المرنة العالمية، وذروة الحمل، وإزاحة الكسر.
4. معيار النمط المختلط المعاير: اقترح نموذجاً لقوة دفع شفاف فيزيائياً يفصل بين المساهمات المتناحية (الحجمية/الانحرافية) والمساهمات متباينة الخواص، حيث يتم معايرة كل منها وفقاً لقساوة الكسر الخاصة بها باستخدام أسس قانون القوة القابلة للضبط.

4. النتائج الرئيسية

أ. التحقق من الصحة (لوحة ذات شق حافة)

• نجح النموذج في إعادة إنتاج الترتيب الرتيب لذروة الأحمال ومنحنيات (القوة-الإزاحة) لمختلف زوايا الشق المسبق ($0^\circ, 30^\circ, 60^\circ$) مقارنة بالبيانات التجريبية.
• لوحظ نقص طفيف في التنبؤ بذروة القوة (10–15%) عُزي إلى افتراضات الإجهاد المستوي ومتوسط معاملات قساوة الكسر.

ب. عينات الشق المفرد (SEN)

• تباين كثافة الشقوق ($\beta$):
  • الدور: يتحكم في اتجاه مسار الشق.
  • السلوك: مع زيادة $\beta$، ينحرف الشق نحو زاوية الألياف ($\theta_f$)، ليصل إلى زاوية الشق الكاملة ($\theta_d = \theta_f$) عند قيم $\beta$ العالية.
  • الاستجابة الميكانيكية: لا يغير الصلابة المرنة الأولية (تظل الاستجابة متناحية)، ولكنه يزيد بشكل كبير من ذروة الحمل وإزاحة الكسر (المتانة).
• طاقة الانفعال متباينة الخواص ($k_1$):
  • الدور: يحرف الشق ولكن مع تشبع سريع.
  • السلوك: تؤدي زيادة $k_1$ إلى انحراف فوري (مثلاً من $0^\circ$ إلى $21^\circ$ عند قيم $k_1$ منخفضة جداً)، لكن الزاوية تتشبع بالقرب من $28^\circ$ (فشل في الوصول إلى زاوية الألياف الكاملة $30^\circ$) حتى عند مستويات الصلابة العالية.
  • الاستجابة الميكانيكية: لا يغير ذروة الحمل أو الصلابة بشكل كبير في هندسات SEN بمجرد أن تكون $k_1$ غير صفرية.
• الاستنتاج لـ SEN: تباين كثافة الشقوق هو العامل المهيمن للتحكم في المسار والمتانة؛ بينما للتباين المرن تأثير محدود ومتشبع على انحراف المسار في هذه الهندسة.

ج. عينات شد ذات ثقب مفتوح (OHT)

• تباين كثافة الشقوق فقط:
  • نجح في حرف مسار الشق.
  • نتيجة حاسمة: له تأثير ضئيل على استجابة (القوة-الإزاحة) العالمية (تظل الصلابة وذروة الحمل متطابقة تقريباً عبر زوايا الألياف المختلفة) لأن المجال المرن متناحٍ. التباين الوحيد يأتي من زيادة طول المسار.
• التباين المرن فقط:
  • يحرف الشق بدرجة أقل من آلية الكثافة.
  • نتيجة حاسمة: يغير الاستجابة الميكانيكية العالمية بشكل كبير. فهو يخلق صلابة أولية، وذروة قوة، وإزاحة كسر تعتمد على زاوية الألياف. الاستجابة الأكثر صلابة عند زوايا الألياف العالية تؤدي إلى فشل مبكر.
• التباين المشترك:
  • التآزر: ينتج عن التأثير المشترك مسارات شقوق تتبع اتجاه الألياف بدقة أكبر مما تفعله أي آلية بمفردها.
  • التفاعل غير الخطي: كانت الفجوة في ذروة القوة بين زوايا الأليفة $0^\circ$ و $90^\circ$ في الحالة المشتركة حوالي 1.6 ضعف مجموع الفجوات الملاحظة في حالات الآليات الفردية. وهذا يؤكد وجود تفاعل تآزري غير خطي حيث يعزز كل من المقاومة المعززة وقوة الدفع المعززة بعضهما البعض.

5. الأهمية والآثار المترتبة

• البصيرة الفيزيائية: توضح الدراسة أن تباين كثافة الشقوق يحكم أين يذهب الشق (مقاومة الكسر/المسار)، بينما طاقة الانفعال متباينة الخواص تحكم مدى قوة دفعه، وفي هندسات تركيز الإجهاد، تتحكم في كيفية تشوه المادة مرناً قبل الفشل.
• استراتيجية المعايرة:
  • عينات SEN غير كافية لمعايرة معامل التباين المرن ($k_1$) لأن الاستجابة العالمية غير حساسة له.
  • عينات OHT (تركيز الإجهاد) ضرورية لمعايرة $k_1$ لأنها تكشف عن تأثيره على الصلابة وذروة الحمل.
• أفضل الممارسات في النمذجة: من أجل محاكاة دقيقة للمواد المعززة بالألياف، يجب تضمين كلتا الآليتين. الاعتماد فقط على تباين كثافة الشقوق يغفل خلط الأنماط الناتج عن المادة وتغيرات الصلابة، بينما الاعتماد فقط على طاقة الانفعال يفشل في التقاط المدى الكامل لانحراف مسار الشق.
• القيود والعمل المستقبلي: نموذج AT1 الحالي حساس للمقياس الطولي. سيتعين على العمل المستقبلي استكشاف دوال تدهور عقلانية لضمان الاستقلال عن المقياس الطولي، وفصل التلين (softening) بين المصفوفة والألياف، والتوسع ليشمل عائلات متعددة من الألياف.