Bonded-particle model for magneto-elastic rods

ملخص تقني لورقة بحثية بعنوان: "نموذج الجسيمات المترابطة للقضبان المغناطيسية المرنة"

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

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

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

تواجه النهج الحالية قيوداً:

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

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

2. المنهجية

قام المؤلفون بتطوير نموذج جسيمات مترابطة (BPM) مخصص للقضبان المغناطيسية المرنة، تم تنفيذه ضمن إطار عمل LAMMPS (المحاكي الذري الجزيئي الضخم واسع النطاق).

أ. صياغة القضيب المرن (BPM الدوراني المشترك)

• التقطيع (Discretization): يتم نمذجة القضيب كسلسلة من الجسيمات الكروية القابلة للتوجيه المتصلة بروابط مرنة.
• الكينماتيكا: يمتلك كل جسيم موضعاً وتوجهاً (يُمثل بـ "كواتيرنيون" - quaternion - وحدة واحدة).
• التفكيك المتماثل: على عكس صياغات BPM السابقة التي تستخدم جسيماً واحداً كإطار مرجعي، يقدم هذا النموذج إطاراً مركزياً ($C$) يتم الحصول عليه عبر متوسط توجهات الجسيمين المترابطين عبر الاستيفاء الخطي الكروي للكواتيرنيون.
• حساب القوة/عزم الدوران: يتم تفكيك الإزاحة النسبية والدوران إلى أربعة أنماط، يرتبط كل منها بنبض (spring) خطي:
  1. التمدد ($K_r$): قوة محورية.
  2. القص ($K_s$): قوة عرضية (مشتقة من الإزاحة النسبية العمودية على الرابط).
  3. الانحناء ($K_b$): عزم دوران مشتق من مكون "الأرجحة" (swing) للدوران.
  4. الالتواء ($K_t$): عزم دوران مشتق من مكون "الالتواء" (twist) للدوران.
• ثوابت النابض: تم اشتقاق الثوابت لتطابق نظرية عارضة أويلر-برنولي ($K_r = EA/L$، $K_s = 12EI/L^3$، $K_t = GJ/L$، $K_b = EI/L$).
• التماثل: يضمن الإطار المركزي صياغة متماثلة تحت تبادل الجسيمات، مما يلغي عيوب القص الدوراني الشائعة في الصياغات غير المتماثلة.

ب. الاقتران المغناطيسي المرن

• تعيين ثنائي القطب: يحمل كل جسيم عزم ثنائي قطب مغناطيسي ($\mathbf{m}$) يدور مع توجه الجسيم ($\mathbf{m} = \mathbf{R}\tilde{\mathbf{m}}$).
• التفاعلات:
  • المجالات الخارجية: تُحسب القوى وعزوم الدوران من طاقة زيمان ($U = -\mathbf{m} \cdot \mathbf{B}$).
  • ثنائي القطب-ثنائي القطب: تُحسب التفاعلات الزوجية باستخدام جهود ثنائي القطب-ثنائي القطب القياسية، مما يسمح بالتفاعلات المغناطيسية بعيدة المدى دون تعديل صياغة الرابط الهيكلي.
• التنفيذ: يستخدم حزمة DIPOLE في LAMMPS وأنماط روابط مخصصة (bpm/corotational).

ج. الديناميكا والتلامس

• التكامل: يتم دمج معادلات الحركة باستخدام خوارزمية "Velocity-Verlet" للترجمة وتحديث "Richardson" للكواتيرنيونات.
• التلامس: يمنع نموذج تلامس "Hertz-Mindlin" الحبيبي التقاطع الذاتي. ولضمان تلامس سلس، تم ضبط أقطار الجسيمات لتكون أكبر قليلاً من طول الرابط.
• التخميد: يتعامل التخميد على مستوى الرابط مع التذبذبات عالية التردد؛ بينما يساعد تخميد "Langevin" العالمي على التقارب نحو التوازن شبه الساكن.

د. التفاعل بين السائل والمنشأ (FSI)

• الاقتران: يتم ربط النموذج بمحلل سوائل يعتمد على طريقة لاتييس بولتزمان (LBM) باستخدام طريقة الحدود المغمورة (IBM).
• الآلية: تتبادل الجسيمات "لاغرانجية" الزخم مع شبكة السائل "الأويلرية" عبر دالة دلتا منتظمة. وتعمل قوى الجزاء (penalty forces) على دفع سرعات الجسيمات نحو سرعات السائل المحلية.
• الدورات الفرعية (Sub-cycling): للحفاظ على الاستقرار، يتم إجراء خطوات متعددة من طريقة العناصر المنفصلة (DEM) لكل خطوة زمنية واحدة من LBM.

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

1. إطار عمل موحد: أول تنفيذ لنموذج جسيمات مترابطة للقضبان المغناطيسية المرنة يوحد التشوهات المرنة الكبيرة، والتلامس، والتفاعلات المغناطيسية بعيدة المدى في كود واحد للعناصر المنفصلة (LAMMPS).
2. صياغة دورانية مشتركة متماثلة: تقديم إطار مرجعي مركزي لتشوه الرابط يضمن الثبات تحت التبادل والاستقرار تحت الالتواء والانحناء الشديدين.
3. القدرة على تعدد الفيزياء: تكامل سلس للتحريك المغناطيسي وديناميكا السوائل (LBM-IBM) دون تغيير الصياغة الهيكلية الأساسية، مما يسمح بمحاكة الأنظمة المعقدة المستوحاة حيوياً.
4. القابلية للتوسع: يستفيد من كفاءة التوازي والبنية التحتية الموجودة في LAMMPS للكهرباء الساكنة طويلة المدى (المكيفة للمغناطيسية الساكنة) وحلول السوائل.

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

تم التحقق من صحة النموذج مقابل ثلاث مشكلات مرجعية وتطبيق واحد متعدد الفيزياء:

• الالتواء الشديد (القضبان المرنة):
  • محاكاة الانبعاج وتشكل الـ plectoneme للقضبان المستقيمة والمنحنية طبيعياً تحت الالتواء الشديد.
  • النتيجة: توافق ممتاز مع البيانات التجريبية (Lazarus et al.) والحلول شبه التحليلية فيما يتعلق بزوايا الالتواء الحرجة وملفات الانحراف.
• التحريك المغناطيسي (المجالات غير المنتظمة):
  • محاكاة العوارض الممغنطة في مجالات منتظمة ومجالات ذات تدرج ثابت.
  • النتيجة: طابقت انحرافات الطرف نتائج طريقة العناصر المحدودة (FEM) (Yan et al.) بدقة عالية، مما أظهر قدرة النموذج على التعامل مع تدرجات المجال المغناطيسي المعقدة.
• تفاعلات ثنائي القطب-ثنائي القطب (القضبان الحلزونية):
  • محاكاة قضبان حلزونية ذات تجاذب مغناطيسي داخلي لدراسة التلاكؤ الميكانيكي (hysteresis).
  • النتيجة: رصد السلوك التلاكتي (عتبات الانكماش مقابل الانفصال) الملاحظ في التجارب، رغم ملاحظة بعض التفاوتات الكمية بسبب حدود التقطيع في تقريب ثنائي القطب النقطي.
• التفاعل بين السائل والمنشأ:
  • التدفق التذبذبي: محاكاة خيوط مرنة في تدفق تذبذبي، مع مطابقة بيانات السعة والتردد التجريبية.
  • الأهداب المغناطيسية: محاكة ضخ السوائل بواسطة مصفوفة من الأهداب المحفزة مغناطيسياً ذات الموجات الميتوكرونية (metachronal waves).
  • النتيجة: إعادة إنتاج أنماط التدفق الدوامي بنجاح وإثبات أن كفاءة الضخ تنخفض مع زيادة الطول الموجي، وهو ما يتوافق مع الملاحظات التجريبية.

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

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

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