🔬 condensed matter

Unjamming in a 3D Granular System: The Micromechanical Role of Friction in Force Distributions and Rheological Properties

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

المؤلفون الأصليون: Vicente Salinas, Héctor Alarcón, Eduardo Rojas, Pablo Gutiérrez, Gustavo Castillo

نُشر 2026-04-23

📖 5 دقيقة قراءة🧠 قراءة متعمّقة

CC BY 4.0

المؤلفون الأصليون: Vicente Salinas, Héctor Alarcón, Eduardo Rojas, Pablo Gutiérrez, Gustavo Castillo

البحث الأصلي مرخَّص بموجب CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ هذا شرح مولَّده بالذكاء الاصطناعي للبحث أدناه. لم يكتبه المؤلفون ولم يصادقوا عليه. وللتحقق من الدقة التقنية، يرجى الرجوع إلى البحث الأصلي. اقرأ إخلاء المسؤولية الكامل

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

عادةً، عندما نفكر في الرمل أو الحبيبات، نعتبرها إما صلبة (مثل قلعة رملية) أو سائلة (مثل سكب الرمل من دلو). لكن هناك منطقة وسط غريبة تسمى "الانسداد" (Jamming). فكر في الأمر كأنه ازدحام مروري؛ عندما تكون السيارات متراصة بإحكام، لا يمكنها الحركة (صلب). وإذا أزلت بعض السيارات، فقد تظل عالقة. ولكن في النهاية، إذا أزلت ما يكفي منها، سيبدأ الخط بأكمله في التدفق (سائل).

أراد الباحثون معرفة كيف "يفك" كوم من الرمل انسداده (يتحول من صلب إلى سائل) عندما تسحب قطعاً منه من الأسفل، بدلاً من هزّه أو سكبه.

إليك قصة تجربتهم، مشروحة ببساطة:

1. الإعداد: "قلعة الرمل الرقمية"

بدلاً من استخدام رمل حقيقي (الذي قد يكون فوضوياً ويصعب تتبعه)، بنوا عالماً رقمياً باستخدام الكمبيوتر.

أنشأوا صندوقاً يحتوي على 30,000 كرة افتراضية (مثل الكرات الزجاجية).
تركوا الكرات تسقط وتستقر تحت تأثير الجاذبية، تماماً مثل الرمل الحقيقي.
ثم بدأوا عملية "هدم محكومة". كل بضع ثوانٍ، كانوا يختارون عشوائياً 30 كرة من النصف السفلي من الصندوق ويحذفونها.
الهدف: مراقبة رد فعل الكرات المتبقية. هل ستظل مكانها؟ هل ستنزلق للأسفل؟ هل ستنهار؟

2. فصول المسرحية (الفصلان)

أظهرت المحاكاة الحاسوبية أن الرمل يتصرف في "فصلين" متميزين:

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

3. المكون السري: الاحتكاك (اللزوجة/التلاصق)

لعب العلماء بمتغير يسمى الاحتكاك.

احتكاك منخفض (كرات ناعمة): تخيل كرات مغطاة بالتفلون. تنزلق بجانب بعضها بسهولة.
احتكاك عالٍ (كرات خشنة): تخيل كرات مغطاة بورق الصنفرة. هي تتشبث ببعضها بقوة.

الاكتشاف:

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

وجد الباحثون قاعدة رياضية: كلما كانت الكرات أكثر "لزوجة" (احتكاكاً)، انخفضت الكثافة التي يمكن أن يصل إليها النظام قبل أن ينهار بالكامل.

4. "سلاسل القوة": الهيكل العظمي غير المرئي

هذا هو الجزء الأكثر إثارة. في كومة الرمل، لا تحمل جميع الكرات نفس الوزن.

تخيل مجموعة من الناس يقفون في دائرة. إذا استند شخص على آخر، فإن ذلك الشخص يستند على التالي، مما يخلق سلسلة من الدعم.
في الرمل، تسمى هذه "سلاسل القوة" (Force Chains). وهي خطوط غير مرئية من الكرات تحمل الحمل الثقيل، بينما تكتفي الكرات الأخرى بالتواجد دون فعل شيء.

ماذا حدث عندما سحبوا الكرات؟

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

5. الصورة الكبيرة: لماذا يهم هذا؟

قد تتساءل، "من يهتم بسحب الكرات من صندوق كمبيوتر؟"

هذا يساعدنا على فهم الكوارث في العالم الحقيقي:

الانهيارات الأرضية: عندما يفقد التربة تماسكه وينزلق من التل.
البالوعات الأرضية: عندما تنهار الأرض فجأة.
الكويكبات: كيف تتصرف أكوام الصخور السائبة في الفضاء عند تعرضها لضربة.

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

الخلاصة

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

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

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

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

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

يعالج المؤلفون التحديات التالية:

كيف يفقد النظام استقراره الميكانيكي عند إزالة الجسيمات بدلاً من تعريضه للقص؟
كيف يؤثر معامل الاحتكاك بين الجسيمات ( $\mu_p$ ) على الكثافة الحجمية الحرجة ( $\phi_c$ ) وعدد التنسيق ( $Z_c$ ) عند نقطة فك الارتباط؟
ما هي التغيرات الميكانيكية الدقيقة في شبكات القوة (سلاسل القوة) وعدم تجانس الإجهاد مع اقتراب النظام من الاستقرار الحدي؟

2. المنهجية

تستخدم الدراسة طريقة العناصر المنفصلة (DEM) باستخدام برنامج MercuryDPM مفتوح المصدر.

إعداد النظام:
- الجسيمات: $N_0 = 30,000$ جسيم كروي (القطر $d=1$ ملم، الكثافة $\rho=2500$ كجم/م³).
- الهندسة: محصور في صندوق مستطيل ( $L_x = L_y = 24d$ ) مع جدران صلبة.
- الفيزياء: تتفاعل الجسيمات عبر نموذج "نابض-مخمد" خطي (طولي وعرضي) مع احتكاك كولوم. تم تضمين احتكاك التدحرج لمراعاة خشونة السطح.
- الجاذبية: تعمل عموديًا، مما يخلق حالة إجهاد متباينة الخواص.
بروتوكول القيادة (استخراج الجسيمات):
- بدلاً من القص، يتم دفع النظام خارج حالة التوازن عن طريق إزالة الجسيمات عشوائيًا من النصف السفلي للحاوية على فترات زمنية منتظمة ( $\Delta t_e = 0.02$ ثانية).
- يحاكي بروتوكول "الاستخراج المدفوع" هذا فقدان الدعم بشكل محكوم، مما يسمح للنظام بالتطور من حالة مرتبطة (تشبه الصلب) نحو حالة مفككة الارتباط (تشبه السائل).
- يتم تغيير معامل الاحتكاك الساكن بين الجسيمات $\mu_p$ كمعامل تحكم رئيسي ( $\mu_p \in [0.2, 0.9]$ ).
المقاييس الرئيسية:
- الكلية (Macroscopic): عدد الجسيمات النشطة ( $N_K$ )، ارتفاع مركز الكتلة ( $z_{cm}$ )، الكثافة الحجمية ( $\phi$ )، وعدد التنسيق ( $Z$ ).
- الدقيقة (Microscopic): توزيعات القوة العمودية ( $P(f_n)$ )، مؤشر اللدونة ( $\zeta = f_t / \mu_p f_n$ )، وإحصائيات سلاسل القوة.
- الأداة الإحصائية: يتم تطبيق معامل جيني (Gini coefficient) على مؤشر اللدونة وكثافة سلاسل القوة لقياس عدم تجانس توزيع الإجهاد.

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

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

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

أ. الأنظمة الديناميكية

يُظهر النظام نظامين متميزين أثناء استخراج الجسيمات:

النظام الأولي (المرتبط): حركة الجسيمات محدودة. يرتفع مركز الكتلة ( $z_{cm}$ ) خطيًا (بسبب الإزالة من الأسفل)، وتنخفض الكثافة الحجمية ( $\phi$ ) خطيًا. تكون إعادة الترتيب موضعية.
النظام المتأخر (فك الارتباط/عدم الاستقرار): يدخل النظام في حالة من إعادة الترتيب الجماعي المستدام. يبدأ $z_{cm}$ في الانخفاض خطيًا (انهيار)، ويصل عدد الجسيمات النشطة ( $N_K$ ) إلى حالة التشبع. يحافظ النظام على كثافة حجمية حرجة ثابتة ( $\phi_c$ ) رغم استمرار فقدان الجسيمات.

ب. المعلمات الهيكلية الحرجة

الكثافة الحجمية الحرجة ( $\phi_c$ ): مع زيادة $\mu_p$ ، تنخفض $\phi_c$ . تتبع العلاقة دالة كسرية:
$\phi_c = \frac{p_1 \mu_p + p_2}{\mu_p + q_1}$
يشير هذا إلى أن الاحتكاك الأعلى يسمح للنظام بالبقاء مستقرًا عند كثافات أقل (قريبة من التعبئة العشوائية الرخوة).
عدد التنسيق الحرج ( $Z_c$ ): تنخفض $Z_c$ مع زيادة $\mu_p$ . وتتماشى العلاقة بين $Z_c$ و $\phi_c$ مع توقعات التعبئة العشوائية الرخوة (RLP) لـ Song et al، مما يؤكد أن الاحتكاك يقلل من عدد التلامسات المطلوبة للاستقرار.

ج. توزيعات القوة وعدم التجانس

القوى العمودية: ينتقل توزيع القوى العمودية ( $P(f_n)$ ) من التحلل الأسي (الذي يميز الحالات الكثيفة) إلى "توزيع أسي ممدد" ذي ذيل أوسع ( $c \approx 0.8$ ) مع اقتراب فك الارتباط. يشير هذا إلى إعادة تنظيم حيث تصبح القوى المتبقية أكثر عدم تجانسًا.
مؤشر اللدونة ( $\zeta$ ): ينخفض متوسط مؤشر اللدونة $\langle \zeta \rangle$ مع زيادة $\mu_p$ . يتناقص معامل جيني لمؤشر اللدونة ( $G_\zeta$ ) بمرور الوقت ولكنه يستقر عند قيمة حرجة تعتمد على $\mu_p$ ، مما يشير إلى أن الاحتكاك الأعلى يؤدي إلى توزيع أحمال أكثر عدم تجانسًا (مزيد من السلاسل "القوية").
سلاسل القوة:
- يتناقص عدد سلاسل القوة ومتوسط شدتها مع اقتراب النظام من فك الارتباط.
- التشبع العالمي: يتقارب معامل جيني لكثافة سلاسل القوة ( $G_F$ ) نحو قيمة عالمية قدرها $\approx 0.284$ لجميع معاملات الاحتكاك عند نقطة فك الارتباط. يشير هذا إلى أن الحالة النهائية لأقصى عدم تجانس في الإجهاد محكومة بالقيود الهندسية وفقدان الاستقرار بدلاً من التبديد الاحتكاكي.

5. الأهمية

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

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