A variational multi-phase model for elastoplastic materials with microstructure evolution

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

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

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

غالباً ما تعجز النماذج الحالية عن وصف:

• الطبيعة الفورية للتحولات الطورية في الزمن.
• التطور المستمر للكسور الحجمية للأطوار الناشئة.
• التفاعل بين التشوه اللدن وخلط الأطوار ضمن إطار متسق ثيرموديناميكياً.

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

2. المنهجية

تستند المنهجية المقترحة إلى ميكانيكا الأوساط المتصلة والمبادئ المتغيرة للمواد غير المرنة. ويتمثل جوهر النهج في تقليل دالة لاغرانج (Lagrangian functional) المكونة من معدل الطاقة الحرة وإمكانات التبديد.

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

• الإطار الثيرموديناميكي: يستخدم النموذج كثافة طاقة هلمهولتز الحرة ($\Psi$) وإمكانية التبديد ($\Delta$). وهو يلتزم بالقانون الثاني للديناميكا الحرارية، وتحديداً مبدأ الحد الأدنى لإمكانية التبديد للمواد غير المرنة المستقلة عن المعدل.
• صياغة متعددة الأطوار:
  • يتم نمذجة المادة كمزيج من $k$ من الأطوار المحتملة بكسور حجمية $\lambda_i$.
  • التحول العشوائي: تُحكم التحولات الطورية بواسطة مصفوفة معدل انتقال $g_{ij}$، والتي تصف التقدم العشوائي من الطور $i$ إلى $j$.
  • مسافة التبديد: لنمذجة التحولات الطورية الفورية، استخدم المؤلفون مفهوم "مسافة التبديد" ($D$)، المشتق من إمكانية التبديد. وهذا يسمح بقفزة فورية في متغيرات الحالة.
  • مقاييس يونج (Young Measures): لضمان التطور المستمر للبنية المجهرية عبر الزمن، يتم تمثيل الكسور الحجمية للأطوار الجديدة باستخدام مقاييس يونج المعممة. ويتعامل هذا مع توزيع المتغيرات الداخلية (مثل الانفعال اللدن) عبر الأطوار.
• الطاقة المسترخية: يتم اشتقاق جهد طاقة حرة مسترخٍ (فعال) ($\Psi_{rel}$) عن طريق تقليل إجمالي الطاقة مع مراعاة قيود الكسر الحجمي. وينتج عن ذلك موترات صلابة فعالة ($C_{eff}$) وانفعالات لدنة فعالة ($\varepsilon^p_{eff}$)، باستخدام صيغة الحد الأدنى لـ Reuß.
• معادلات التطور: من خلال تقليل لاغرانج بالنسبة للمتغيرات الحركية (معدلات الانفعال اللدن $\dot{\varepsilon}^p_i$ ومعدلات الانتقال $g_{ij}$)، يشتق المؤلفون:
  • قواعد التدفق: التي تحكم التشوه اللدن داخل كل طور.
  • قواعد الانتقال: التي تحكم تطور الكسور الحجمية ($\dot{\lambda}_i$).
  • شروط الخضوع: حيث يتم تعريف دالات خضوع متميزة للتدفق اللدن ($\phi_i$) وللتحولات الطورية ($\phi_{ij}$).
• التنظيم (Regularization): تُضاف حدود تنظيم لزوجية إلى إمكانية التبديد لضمان السلامة الرياضية والمتانة العددية.

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

• نموذج متغير موحد: تقدم الورقة إطاراً جديداً يجمع بين مسافة التبديد (للتحولات الفورية) ومقاييس يونج (للتطور المستمر للكسر الحجمي) ضمن إطار متغير واحد.
• مصفوفة الانتقال العشوائي: يسمح إدخال مصفوفة الاحتمالات ($g_{ij}$) بنمذجة التفاعلات الطورية المعقدة ومتعددة الاتجاهات وبدء أطوار جديدة بناءً على القوى الدافعة الثيرموديناميكية.
• التنفيذ الخوارزمي: طور المؤلفون روتيناً للتحديث التدريجي (الخوارزمية 1) يعالج بدء الأطوار الجديدة، وحساب الخصائص الفعالة، وتحديث حالات الإجهاد.
• التحقق العددي: تم تنفيذ النموذج باستخدام طريقة العناصر المحدودة (FEM) مع أدوات مفتوحة المصحة (Gmsh و Julia Ferrite)، مما يثبت قابليته للتطبيق على مسائل القيم الحدية.

4. النتائج

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

أ. سلوك النقطة المادية (اختبار دوري 1D)

• الإعداد: نقطة مادية واحدة مكونة من ثلاثة أطوار تعرضت لإجهاد شد دوري ($\pm 6.5\%$ انفعال).
• الملاحظات:
  • بدأت الأطوار الجديدة ($\lambda_1, \lambda_2$) فور الوصول إلى إجهاد الخضوع.
  • أظهر النموذج سلوك الاستقرار (shakedown behavior): فبعد التشوه اللدن الأولي وتطور الطور، استقر النظام في استجابة مرنة بحتة مع كسور حجمية ثابتة.
  • تقاربت الكسور الحجمية نحو حالة مستقرة، مما أظهر قدرة النموذج على التقاط تشبع التغيرات في البنية المجهرية.

ب. اختبار مرجعي ثنائي الأبعاد

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

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

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

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