🔢 mathematics

A Variational Formulation for Deformable Particle Simulations and its Level Set Discrete Element Method Implementation

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

المؤلفون الأصليون: Thomas Henzel, Konstantinos Karapiperis

نُشر 2026-02-16

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

CC BY 4.0

المؤلفون الأصليون: Thomas Henzel, Konstantinos Karapiperis

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

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

فكر في طريقة (DEM) التقليدية كأنك تلعب بـ كرات الزجاج (الماربل). يمكنك دحرجتها، وجعلها ترتد، وتكديسها. الكمبيوتر يعرف بالضبط مكان كل كرة وكيف تصطدم ببعضها البعض. لكن هناك مشكلة: في هذه المحاكاة، الكرات صلبة تماماً. إذا ضغطت كرتين معاً، فإنهما لا تنبعجان؛ بل تتوقفان ببساطة. وإذا أردت رؤية كرة وهي تتسطح أو تنثني، فإن البرامج القديمة لا يمكنها فعل ذلك دون تعطل المحاكاة أو استغراق ملايين السنين في الحسابات.

تقدم هذه الورقة البحثية الجديدة ترقية ذكية: طريقة العناصر المنفصلة القابلة للتشوه (Deformable DEM). الأمر يشبه تحويل كرات الزجاج الصلبة تلك إلى كرات ضغط (Stress balls) أو قطع حلوى "جيمي بير" (Gummy bears). الآن، عندما تصطدم، فإنها تنبعج، وتتسطح، وترتد، تماماً مثل المادة الحقيقية.

إليك كيف فعلوا ذلك، مقسماً إلى مفاهيم بسيطة:

1. "المغير السحري للشكل" (الصيغة التباينية - The Variational Formulation)

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

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

لقد أدركوا أن معظم الأجسام لا تلتوي في أشكال غريبة وعشوائية، بل تنثني أو تتمدد أو تنضغط بطرق يمكن التنبؤ بها. لذا، بدلاً من تتبع ملايين النقاط، هم يتتبعون مجرد عدد قليل من "الأنماط" (Modes) (مثل النوتات الموسيقية).

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

2. "المخطط الشبح" (طريقة المستوى - Level Set Method)

عندما تنبعج كرة، يتغير شكلها. كيف يعرف الكمبيوتر أين أصبح الحرف الآن؟

يستخدمون ما يسمى بـ "المستوى" (Level Set).

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

3. لماذا هذا مهم: حل "غولدي لوكس" (الحل المثالي)

قبل هذا، كان لدى العلماء خياران سيئان:

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

هذه الطريقة الجديدة هي حل "غولدي لوكس" (الحل المثالي الذي يجمع بين السرعة والواقعية):

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

أمثلة من الواقع من الورقة البحثية

اختبر المؤلفون هذا في سيناريوهين:

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

الخلاصة

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

باختصار: لقد علموا الكمبيوتر كيف يلعب بـ قطع الحلوى (Gummy bears) بدلاً من كرات الزجاج، باستخدام اختصار ذكي يجعل المحاكاة تعمل بسرعة كافية لنمذجة أكوام كاملة من الرمل، أو الصخور، أو حتى الأنسجة البيولوجية دون الحاجة إلى كمبيوتر خارق.

إليك ملخص تقني مفصل للورقة البحثية بعنوان: "صياغة تباينية لمحاكاة الجسيمات القابلة للتشكل وتطبيقها باستخدام طريقة العناصر المنفصلة لمستوى الطقم (Level Set Discrete Element Method)."

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

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

تواجه النهج الحالية لنمذجة الجسيمات القابلة للتشكل مقايضات كبيرة:

الاقتران بين طريقة العناصر المحدودة وطريقة العناصر المنفصلة (FEM-DEM): تحل مسائل القيم الحدية الداخلية لكل جسيم باستخدام طريقة العناصر المحدودة (FEM). ورغم دقتها، إلا أنها مكلفة حاسوبياً للغاية للأنظمة الكبيرة بسبب الحاجة إلى شبكات (meshes) فردية وإعادة بناء الشبكة.
نماذج العناقيد/الزنبركات (Cluster/Spring Models): تمثل الجسيمات كمجموعات من الوحدات الصلبة الأصغر المتصلة بنوابض. تتطلب هذه النماذج العديد من الجسيمات الفرعية لتحقيق تشوه واقعي، مما يؤدي إلى تكاليف حاسوبية عالية.
النماذج ذات الرتبة المختزلة (Reduced-Order Models): تستخدم كينماتيكا مبسطة أو مصفوفات امتثال، لكنها غالباً ما تفتقر إلى اتصال صارم بالميكانيكا الأساسية، وتكون مقيدة بأشكال هندسية محددة، أو تفشل في التعامل مع اللاخطية بشفافية.

التحدي يكمن في تطوير إطار عمل يمكنه التقاط التشوه الحجمي على مستوى الحبيبة مع الحفاظ على الكفاءة الحسابية، والقابلية للتوسع، والوضوح الهندسي لـ DEM الكلاسيكية للأجسام الصلبة.

2. المنهجية

يقترح المؤلفون إطار عمل الـ DEM التبايني القابل للتشكل ضمن صياغة طريقة مستوى الطقم (LS-DEM). تتكون المنهجية من ثلاث ركائز أساسية:

أ. الصياغة التباينية (مبدأ لاغرانج-دالمبير)

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

متغيرات الحالة: يتم وصف كل حبيبة عبر الموقع ( $X$ )، والدوران ( $\Theta$ )، ومجموعة من الإحداثيات النمطية المعممة ( $\varepsilon_\alpha$ ) التي تمثل التشوهات المرنة الداخلية.
الكينماتيكا المختزلة: يتم تقريب تشوه الجسيم كتركيبة خطية من دوال شكل مكانية ثابتة (أنماط) $\Phi_\alpha(x)$ مع سعات تعتمد على الزمن $\varepsilon_\alpha(t)$ :
$u_\varepsilon(x, t) = \varepsilon_\alpha(t) \Phi_\alpha(x)$
الطاقات: يشتق الإطار مصطلحات الطاقة الحركية والجهدية لهذه الأنماط التشوهية.
- الطاقة الحركية: تُعرف عبر مصفوفة كتلة معممة $M_{\alpha\beta}$ .
- طاقة الجهد: تُعرف عبر مصفوفة صلابة معممة $K_{\alpha\beta}$ (للمرونة الخطية) أو دالة طاقة اختيارية $U_\varepsilon(\varepsilon_\alpha)$ (لاستيعاب اللاخطية).
المعادلات الحاكمة: بتطبيق مبدأ لاغرانج-دالمبير، نحصل على نظام من المعادلات التفاضلية العادية (ODEs) التي تربط حركة الجسم الصلب بالتطور النمطي:
$M_\varepsilon \ddot{\varepsilon} + K_\varepsilon \varepsilon = F_\varepsilon$
حيث تمثل $F_\varepsilon$ القوى المعممة المشتقة من تفاعلات التلامس.

ب. تنفيذ مستوى الطقم (LS-DEM)

يستخدم إطار العمل صياغة LS-DEM، حيث يتم تمثيل هندسة الجسيم بواسطة دالة مستوى الطقم $\phi(x)$ على شبكة وعقد سطحية.

تحديث الهندسة: مع تشوه الجسيمات، يجب أن تتطور كل من عقد السطح ومجال مستوى الطقم.
استراتيجيات تطور مستوى الطकम: يقيم البحث عدة طرق لتحديث مستوى الطقم تحت تأثير التشوه:
1. إعادة الحساب (Recomputation): دقيقة ولكنها مكلفة.
2. الادراف (Advection): يحل معادلة نقل؛ يتطلب استقراء السرعة.
3. الادراف شبه لاغرانجي (المختار): يقوم برسم خرائط لنقاط الشبكة الحالية إلى التكوين المرجعي باستخدام مجال الإزاحة العكسي ( $x_{ref} = x - u(x,t)$ ) ويقوم باستكمال مستوى الطقم الأولي. هذا فعال حاسوبياً ويتجنب حل المعادلات التفاضلية الجزئية في كل خطوة.
4. الاستكمال المسبق (Precomputed Interpolation): يستخدم مستويات طقم محسوبة مسبقاً لحالات تشوه محددة (تكلفة ذاكرة عالية).

ج. استخراج الأنماط

يدعم الإطار كلاً من الأنماط التحليلية (تعبيرات مغلقة، مثل انحناء العوارض) والأنماط العددية المستخرجة من محاكاة FEM عالية الدقة أو التجارب.

للأنظمة المعقدة، يتم استخراج الأنماط عبر تحليل التفكيك المتعامد للمكونات الرئيسية (POD) لمجالات الإزاحة من محاكاة الخلية الوحدة أو عنصر الحجم التمثيلي (RVE).
يتم التعامل مع تأثيرات التلامس غير الخطية عبر بناء دالة طاقة غير خطية أو جدول بحث (lookup table) للعلاقة بين القوة والإزاحة بناءً على الأنماط المستخرجة.

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

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

4. النتائج والتحقق

تحقق المؤلفون من إطار العمل من خلال عدة اختبارات مرجعية:

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

5. الأهمية

يجسر هذا العمل الفجوة بين ميكانيكا المتصل عالية الدقة (FEM) والطرق المنفصلة الفعالة (DEM).

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