🔬 condensed matter

multisphere: a Python implementation of the Multi Sphere Shape generator (MSS) for DEM simulations

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

المؤلفون الأصليون: Felix Buchele, Patric Müller, Thorsten Pöschel

نُشر 2026-04-23

📖 4 دقيقة قراءة☕ قراءة في استراحة قهوة

CC BY 4.0

المؤلفون الأصليون: Felix Buchele, Patric Müller, Thorsten Pöschel

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

تخيل أنك تحاول بناء نسخة طبق الأصل من جسم معقد — مثل بقرة، أو حبة رمل، أو أرنب — ولكن يُسمح لك فقط باستخدام كرات رخامية مستديرة تماماً.

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

إليك قصة multisphere، الأداة الموصوفة في الورقة، مشروحة بكلمات بسيطة:

1. المشكلة: محدودية "الكرات الرخامية المستديرة"

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

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

2. الحل: "المترجم من البكسل إلى الرخام"

ابتكر المؤلفون أداة بلغة بايثون تسمى multisphere. فكر فيها كأنها طابعة ثلاثية الأبعاد لا تستخدم سوى الكرات الرخامية.

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

3. لماذا هذه الأداة مميزة؟ (الجزء "الذكي")

هناك أدوات أخرى تحاول القيام بذلك، لكن multisphere مميزة لسببين:

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

4. كيف تعمل؟ (الخطوات الثلاث)

تتبع الأداة وصفة بسيطة مكونة من ثلاث خطوات:

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

5. لماذا يجب أن نهتم؟

تخيل أنك تدرس كيفية انهيار جسر تحت حركة مرور كثيفة، أو كيفية تحرك انزلاق أرضي أسفل جبل.

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

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

إليك ملخص تقني مفصل للورقة البحثية المتعلقة بحزمة البرمجيات multisphere بلغة بايثون.

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

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

التحدي: محاكاة الجسيمات غير الكروية مكلفة حاسوبياً لأن اكتشاف التلامس وتقييم القوى المعقدة يصعب حسابهما مباشرة.
الحل الحالي: يقوم "نهج تعدد الكرات" (multi-sphere approach) بتقريب الأشكال المعقدة ( $S$ ) كunion (اتحاد) لـ $n$ من الكرات المتداخلة ( $\tilde{S}$ ). وهذا يقلل التفاعلات المعقدة إلى تفاعلات بسيطة بين كرة وأخرى.
الفجوة: إن إيجاد المجموعة المثلى من مواضع الكرات وأنصاف أقطارها لتقليل عدم التطابق الهندسي بين الشكل المستهدف والتقريب هو مشكلة تحسين صعبة. غالباً ما تتطلب المولدات الموجودة عدداً كبيراً من الكرات لتحقيق دقة عالية، أو تفتقر إلى الحفاظ على التماثل، أو لا يمكن دمجها بسهء في سير عمل DEM الحديث القائم على بايثون.

2. المنهجية

تقدم الورقة البحثية multisphere، وهي حزمة بايثون مفتوحة المصدر تنفذ خوارزمية مولد "شكل تعدد الكرات" (MSS).

الخوارزمية الأساسية (MSS)

تحل خوارزمية MSS مشكلة التحسين المتمثلة في تقليل الفرق المتماثل (عدم التطابق) بين الشكل المستهدف $S$ وتقريب تعدد الكرات $\tilde{S}$ :
$|S \triangle \tilde{S}| \equiv |(S \cup \tilde{S}) \setminus (S \cap \tilde{S})|$
تشمل الميزات الرئيسية لخوارزمية MSS ما يلي:

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

بنية البرمجيات

تتبع الحزمة سير عمل يتكون من ثلاث مراحل:

المعالجة المسبقة:
- تقبل المدخلات كـ أسطح مثلثية (STL, OBJ, PLY) أو أحجام فوكسلية (مصفوفات NumPy).
- يتم تحويل الأسطح داخلياً إلى VoxelGrid (شبكة فوكسل ثلاثية الأبعاد) مع حجم فوكسل ( $\Delta$ ) مشتق من صندوق الإحاطة (bounding box).
المعالجة:
- تقوم الدالة الأساسية multisphere_from_voxels (أو multispher_from_mesh) بتوليد مجموعة الكرات.
- معايير الإنهاء: يتوقف إعادة البناء بناءً على حدود يحددها المستخدم: الدقة المستهدفة (معامل Dice-Sørensen)، أو الحد الأدنى لنصف قطر الكرة المقبول، أو الحد الأقصى لعدد الكرات.
- حصر الشبكة (Mesh Confinement): يمنع العلم (flag) الاختياري (confine_mesh) التقريب من الامتداد خارج سطح الشبكة الأصلية عن طريق تقليل أنصاف أقطار الكرات، وهو أمر بالغ الأهمية لمحاكاة تكسر الجسيمات لتجنب التداخلات غير الفيزيائية.
المعالجة اللاحقة:
- تُخرج كائن SpherePack يحتوي على مواضع الكرات وأنصاف أقطارها بالوحدات الفيزيائية.
- توجد أدوات لتحويل مجموعة الكرات مرة أخرى إلى شبكة سطحية (create_multisphere_mesh) لأغراض التصور وتقييم الدقة.

التبعيات التقنية

اللغة: Python $\ge$ 3.10.
المكتبات الأساسية: numpy, scipy, scikit-image, trimesh.
التصور: pyvista.
عمليات الشبكة: manifold3D.

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

تنفيذ مفتوح المصدر: توفر أول تنفيذ متاح بلغة بايثون لخوارزمية MSS، مما يجعلها متاحة عبر PyPI (pip install multisphere).
واجهة موحدة: تقدم سير عمل سلساً لتحويل كل من شبكات CAD/التصوير المقطعي والبيانات الفوكسلية إلى تمثيلات تعدد كرات جاهزة لمحاكاة DEM.
التحسين: تقدم أشكال جسيمات عالية الدقة باستخدام عدد أقل من الكرات مقارنة بالطرق الأخرى، مما يقلل مباشرة من التكلفة الحسابية لمحاكاة DEM (حيث يمثل عدد تفاعلات كرة-كرة عنق زجاجة).
أدوات تقييم قوية: تتضمن وظائف مدمجة لحساب معامل Dice-Sørensen ( $D$ ) للتقييم الكمي للدقة وأدوات للتصور ثلاثي الأبعاد التفاعلي.

4. النتائج والتحقق من الصحة

تتحقق الورقة من الحزمة باستخدام دراستي حالة رئيسيتين:

النماذج الهندسية القياسية: تمت إعادة بناء "Stanford Bunny" بأعداد كرات متغيرة ( $n=10, 50, 629$ ). تُظهر عمليات التصور التحسن التدريجي في دقة الشكل مع زيادة $n$ .
المواد الحبيبية من العالم الحقيقي: تم اختبار الحزمة على 402 جسيماً من "Sand Atlas" (بما في ذلك رمال هامبورغ، والكسكسي، وأنواع مختلفة من الرمال).
- الدقة: تم رسم معامل Dice ( $D$ ) مقابل عدد الكرات. تُظهر النتائج اتجاهاً واضحاً حيث تزداد الدقة مع زيادة عدد الكرات، مع حساب القيمة المتوسطة والانحراف المعياري عبر جميع أنواع الجسيمات.
- الأداء: يشير تحليل وقت التشغيل إلى أن الخوارزمية فعالة حاسوبياً، وهي مناسبة للتطبيقات التي تتطلب إعادة بناء متكررة للجسيمات (مثل محاكاة التفتت).

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

تعالج حزمة multisphere عنق زجاجة حرج في محاكاة DEM: وهي المقايضة بين دقة شكل الجسيم والتكلفة الحسابية.

الكفاءة: من خلال توليد تمثيلات دقيقة بعدد أقل من الكرات، فإنها تقلل بشكل كبير من عدد حسابات التلامس المطلوبة، مما يسرع المحاكاة التي تتضمن أعداداً كبيرة من الجسيمات غير المنتظمة.
تعدد الاستخدامات: تدعم المحاكاة الديناميكية حيث تتغير أشكال الجسيمات (مثل التشوه اللدن أو التكسر)، حيث يمكن لخوارزمية MSS إعادة بناء الأشكال الجديدة بكفاءة أثناء التشغيل.
سهولة الوصول: من خلال تغليف خوارزمية هندسة حاسوبية معقدة في حزمة بايثون سهلة الاستخدام مع دعم لتنسيقات الملفات القياسية (STL, VTK, CSV)، فإنها تخفض حاجز الدخول أمام الباحثين في علوم المواد والميكانيكا الجيولوجية لإجراء محاكاة DEM عالية الدقة ومرتبطة بالشكل.