Finite compressibility and strain hardening in elasto-plastic models of amorphous matter

إليك ملخص تقني مفصل للورقة البحثية بعنوان "القدرة على الانضغاط المحدودة وتصلد الانفعال في النماذج المرنة-اللدنة للمادة غير المتبلورة" لـ إلغايلياني، وفانديمبروك، ومالوني.

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

تُحكم عملية الخضوع اللدن والمرونة في المواد غير المتبلورة (مثل الزجاج المعدني، والزجاج الغرواني، والمستحلبات) بواسطة تحولات قص موضعية. وهناك عامل حاسم، غالباً ما يتم تجاهله في النماذج النظرية، وهو القدرة على الانضغاط المحدودة، والتي تتميز بالنسبة بين معامل الحجم ($K$) ومعامل القص ($\mu$).

• الفجوة: بينما ترتبط نسبة بواسون ($\nu$) أو النسبة $K/\mu$ تجريبياً بليونة الزجاج المعدني (حيث ترتبط النسبة العالية لـ $K/\mu$ بليونة أعلى)، فإن معظم النماذج الميزوسكوبية (متوسطة المقياس) المرنة-اللدنة الحديثة تعامل المواد كأنها غير قابلة للانضغاط أو تتجاهل الاعتماد الكمي لإعادة توزيع الإجهاد على القدرة على الانضغاط.
• السؤال: كيف يؤثر تغيير القدرة على الانضغافة اللامتناهية ($K/\mu$) على سلوك الخضوع العياني، وتحديداً الانتقال من الأنظمة المرنة إلى الأنظمة اللدنة في القص الدوري وفي تدفق القص الأمامي؟

2. المنهجية

يستخدم المؤلفون نموذجاً آلياً مرناً-لدناً ميزوسكوبياً يعتمد على رصف (tiling) لشبكة ممتلئة من مواقع تحول القص المتطابقة.

• بناء النموذج:
  • النظام عبارة عن شبكة مربعة ثنائية الأبعاد ($L \times L$) حيث يمتلك كل موقع دالة طاقة انفعال تربيعية مجزأة.
  • يحدث الخضوع الموضعي عندما يتجاوز الانفعال المرن الموضعي عتبة معينة (المحددة بالواحد الصحيح). وعند الخضوع، يخضع الموقع لزيادة منفصلة في الانفعال اللدن (الانفعال الذاتي).
  • تجزئة الشبكة المتدرجة (Staggered Grid Discretization): لحل معادلات المرونة الخطية، استخدم المؤلفون مخطط شبكة متدرجة (شائع في علم الزلازل) لضمان التوافق الكينماتيكي. وهذا يتجنب السلوكيات المرضية التي تظهر في مخططات الفرق المركزي القياسية فيما يتعلق بمكونات الإجهاد غير القطرية.
• مشكلة تضمين إسليفيه (Eshelby Inclusion Problem):
  • يدمج النموذج إجهاد الخلف لإسليفيه ($\sigma_0$). فعندما يخضع موقع ما للخضوع، فإنه يحفز مجال إجهاد في المصفوفة المحيطة.
  • يعتمد مقدار إعادة توزيع الإجهاد هذا على المعامل غير البعدي $\alpha = \frac{2K}{K+\mu}$.
  • يتراوح $\alpha$ من 0 (أقصى قدرة على الانضغاط) إلى 2 (غير قابل للانضغاط).
• بروتوكولات المحاكاة:
  1. القص الدوري: يتعرض النظام لقص تذبذبي بسعة $\Gamma$. وتركز الدراسة على دورات الحد النهائي.
  2. تدفق القص الأمامي المستقر: يتم إخضاع النظام للقص بشكل أحادي الاتجاه حتى يصل إلى تدفق مستقر.

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

1. الربط الكمي بين القدرة على الانضغاط والخضوع: تثبت الورقة أن إجهاد الخلف لإسليفيه ($\sigma_0$) هو الآلية الأساسية التي تربط القدرة على الانضغاط ($K/\mu$) بالخضوع العياني. ومع زيادة $K/\mu$ (أي تصبح المادة أقل قدرة على الانضغاط)، يصبح $\sigma_0$ أكثر سلبية.
2. التصلد الناشئ دون قواعد مضافة (Ad-hoc): يظهر النموذج سلوكيات تصلد معقدة (متناحية وحركية) ناتجة بحت من التفاعلات المرنة وإجهاد الخلف لإسليفيه، دون إدخال معاملات أو قواعد تصلد مجهرية.
3. قوانين القياس الموحدة: يستنتج المؤلفون علاقات قياس لـ السعة الدورية الحرجة ($\Gamma_0$) وإجهاد التدفق في الحالة المستقرة ($\sigma_y$) بدلالة $\sigma_0$ و $\alpha$، مما يوضح كيفية اعتماد هذه الكميات العيانية على المعاملات المرنة المجهرية.

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

أ. سلوك القص الدوري

يظهر النظام ثلاثة أنظمة مع زيادة سعة الانفعال $\Gamma$:

1. النظام المرن (ER): $\Gamma < \Gamma_0$. لا توجد إعادة ترتيب لدنة؛ ويعود النظام إلى حالته الأصلية.
2. نظام اللدونة العكوس (RPR): $\Gamma_0 < \Gamma < \Gamma_y$. يدخل النظام في دورات حدية تتضمن إعادة ترتيب لدنة ولكنها عكوسة مجهرياً (تعود إلى التكوينات السابقة).
3. النظام الانتشاري (DR): $\Gamma > \Gamma_y$. حركة انتشارية غير عكوسة.

النتائج المتعلقة بـ $\Gamma_0$ (انتقال التلاكؤ):

• الاعتماد على القدرة على الانضغاط: تزدزد سعة الانتقال $\Gamma_0$ مع انخفاض $K/\mu$ (أي أن الأنظمة الأكثر قدرة على الانضغاط تخضع في وقت مبكر). وعلى العكس من ذلك، فإن الأنظمة الأقل قدرة على الانضغاط (التي تمتلك $K/\mu$ أعلى) هي "أكثر صلابة" فعلياً، حيث تظهر نطاقاً مرناً أكبر.
• العلاقة: يتم تحديد $\Gamma_0$ بواسطة إجهاد الخلف:
  $$\Gamma_0 = \frac{1 + \sigma_0}{2}$$
  بما أن $\sigma_0$ يصبح أكثر سلبية مع زيادة $K/\mu$، فإن $\Gamma_0$ يزداد.
• بداية اللدونة: عند $\Gamma > \Gamma_0$، تبدأ اللدونة عند انفعال $\gamma^* = \Gamma_0 - (\Gamma - \Gamma_0)$. وتقفز معدلات الانفعال اللدي بشكل حاد إلى $\approx 0.2$ عند $\gamma^*$.
• توزيع الإجهاد: يطور توزيع احتمالية الإجهادات الموضعية، $P(\sigma)$، "كتفاً" مميزاً عند $\sigma_0$. هذه الميزة هي التي تحدد بدء اللدونة.

ب. سلوك تدفق القص الأمامي المستقر

• إجهاد التدفق ($\sigma_y$): ينخفض إجهاد التدفق في الحالة المستقرة مع زيادة $K/\mu$.
• القياس (Scaling): على غرار القص الدوري، يتم التحكم في إجهاد التدفق بواسطة إجهاد الخلف:
  $$\sigma_y = \sigma'_y + \sigma_0(1 - \sigma'_y)$$
  حيث $\sigma'_y \approx 0.54$ هو ثابت عالمي مشتق من توزيع الإجهاد المعاد قياسه.
• انهيار البيانات (Data Collapse): من خلال إعادة قياس الإجهاد والانفعال باستخدام $\sigma_0$ و $\alpha$، تنهار منحنيات التحميل لمختلف درجات القدرة على الانضغاط لتندمج في منحنيات رئيسية، مما يؤكد عالمية الآلية.

ج. التصلد الناشئ عن الانفعال

يعيد النموذج طبيعياً إنتاج نظامين متميزين من أنظمة التصلد الملاحظة في اللدونة المتصلة:

1. التصلد المتناحي (Isotropic Hardening) ($\Gamma < \Gamma_0$): يتوسع النطاق المرن بشكل متماثل أثناء تدوير النظام.
2. التصلد الحركي (Kinematic Hardening) ($\Gamma > \Gamma_0$): يترجم النطاق المرن في فضاء الإجهاد دون تغيير حجمه. ويحدث الانتقال بين هذين النظامين تماماً عند السعة الحرجة $\Gamma_0$.

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

5. الأهمية والتداعيات

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

باختصار، توضح الورقة أن نسبة بواسون (عبر $K/\mu$) تتحكم في سعة إجهاد الخلف لإسليفيه، والذي بدوره يحدد حدود الخضوع العيانية، وإجهاد التدفق، وطبيعة تصلد الانفعال في المواد غير المتبلورة.