Evaluation of External Magnetic Flux Density in Piezo-Flexomagnetic Nanobeams Using a Hybrid 1D-2D Finite Element Framework

ملخص تقني: تقييم كثافة التدفق المغناطيسي الخارجي في العوارض النانوية ذات الخصائص الكهروضغطية-المرنة المغناطيسية (Piezo-Flexomagnetic)

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

المنهجية
يقترح المؤلفون إطار عمل هجين (1D-2D) يعتمد على العناصر المحدودة لتقييم كثافة التدفق المغناطيسي الخارجي ($B$) وشدة المجال ($H$) الناتج عن انحناء عارضة نانوية كهروضغطية-مرنة مغناطيسية.

• الإطار المقترن: يربط النموذج بين نموذج عارضة "تيموشينكو" أحادي البعد (1D) (الذي يحكم المحور الهيكلي المتعادل) ومسألة مغناطيسية ستاتيكية ثنائية الأبعاد (2D) تشمل جسم العارضة وفضاء الهواء المحيط.
• المعادلات الحاكمة:
  • الميكانيكا: تتضمن حركية العارضة الإزاحة المحورية ($u$)، والانحراف المستعرض ($w$)، والدوران ($\phi$). وتتضمن المعادلات التكوينية الاقتران الكهروضغطي المغناطيسي ($d_{31}, d_{15}$) والاقتران المرن المغناطيسي ($f_{31}$)، مما يربط الانفعال وتدرجات الانفعال بالاستقطاب المغناطيسي والإجهاد.
  • المغناطيسية: يتم حل الجهد المغناطيسي القياسي ($\psi$) في كل من نطاق العارضة ونطاق الهواء. ويفرض الصياغة استمرارية الجهد المغناطيسي القياسي عبر الواجهة بين العارضة والهواء مع تطبيق شروط "ديريكليه" للمجال البعيد.
• التنفيذ العددي: يتم حل النظام باستخدام الحل البرمجي مفتوح المصدر FreeFEM++.
  • التقسيم (Discretization): تُستخدم عناصر لاغرانج التربيعية المستمرة (P2) للمتغيرات الحركية والجهد المغناطيسي. وتُستخدم عناصر خطية غير مستمرة (P1dc) في مرحلة ما بعد المعالجة للتعامل مع قفزات التدرج عند الواجهات.
  • خوارزمية الحل: يتم استخدام مخطط تكرار النقطة الثابتة المتدرج. حيث يتم دمج الجهد المغناطيسي ثنائي الأبعاد عبر سمك العارضة لتوليد أحمال مغناطيسية معممة (عزوم وقوى قص) للمحلل الهيكلي أحادي البعد. وعلى العكس من ذلك، يتم رسم إزاحات الهيكل أحادية البعد في النطاق ثنائي الأبعاد لتعمل كمصادر للجهد المغناطيسي.
• التحقق: تم التحقق من صحة النموذج مقابل الحلول التحليلية للعارضات المعزولة مغناطيسياً عن طريق كبح نفاذية نطاق الهواء، مما أظهر توافقاً جيداً مع الأدبيات الموجودة.

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

• وجود المجال الخارجي: تثبت الدراسة وجود توزيعات تدفق مغناطيسي خارجي ملموسة في الهياض الحرة، حتى في غياب الاقتران الكهروضغطي المغناطيسي (حالة المرونة المغناطيسية الصرفة).
• خصائص المجال:
  • في العوارض المرنة المغناطيسية، يخلق تدرج الانفعال الثابت عبر السمك مصدراً مغناطيسياً مميزاً عند السطح السفلي ومصرفاً عند السطح العلوي.
  • تظهر كثافة التدفق المغناطيسي ($B$) وشدة المجال ($H$) سلوكيات مختلفة داخل العارضة مقارنة بالهواء. فداخل العارضة، ينعكس اتجاه $H$ بالنسبة لـ $B$ بسبب الاقتران المغناطيسي المرن القوي، بينما يميل $B$ إلى التركز داخل منطقة العارضة ذات النفاذية الأعلى، مكوناً أنماطاً حلقية.
  • يكون المجال الخارجي مستمراً في المكون العمودي لـ $B$ ولكن غير مستمر في المكون المماسي عبر الواجهة، وهو ما يتوافق مع الشروط الكهرومغناطيسية الحدودية.
• تحليل الحساسية: يحدد تحليل منهجي المعلمات المادية الأكثر تأثيراً على كثافة التدفق المغناطيسي المستعرض الخارجي ($B_z$):
  • النفاذية ($\mu_0$): تعد نفاذية نطاق الهواء الخارجي العامل الأكثر هيمنة؛ حيث يؤدي الوسط الخارجي الأكثر نفاذية إلى إنتاج $B_z$ أعلى.
  • معاملات الكهروضغطية المغناطيسية: في العوارض الكهروضغطية-المرنة المغناطيسية، يكون معامل اقتران القص ($d_{15}$) أكثر تأثيراً على التدفق الخارجي من اقتران الانفعال العمودي ($d_{31}$). ويعزى ذلك إلى إلغاء تأثيرات الانفعال العمودي عبر السمك أثناء الانحناء، بينما يظل انفعال القص ثابتاً.
  • معامل المرونة المغناطيسية: في العوارض المرنة المغناطيسية الصرفة، يعد المعامل $f_{31}$ (الذي يربط تدرج الانفعال بالتدفق المغناطيسي) هو المعلم المهيمن.

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