A unified rate-dependent anisotropic hysteresis model for grain-oriented, non-oriented, and directional electrical steels

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

يقدم هذا البحث طريقة موحدة جديدة لنمذجة هذه الرقصة المغناطيسية. قام المؤلفان، ماريو بندرا وفلوريان بولتشاك، ببناء محرك رياضي متطور يعتمد على مفهوم يسمى معادلة فوكر-بلانك (Fokker–Planck equation). فكر في هذا كأنه توقعات جوية عالية الدقة للأسهم المغناطيسية. بدلاً من التخمين أين تذهب، يتتبع هذا النموذج احتمالية توجيه كل سهم في كل اتجاه ممكن على سطح كرة. وهو يجمع بين ثلاثة "قوى" توجه الأسهم: قوة جذب المغناطيس الخارجي (Zeeman)، وتفضيل اتجاه بلوري معين (uniaxial)، والبنية المعقدة التي تشبه الشبكة للمعدن نفسه (cubic). والنتيجة هي مجموعة واحدة مشتركة من القواعد التي تعمل لكل من الاتجاهين "الأمامي" و"الجانبي" للصلب عالي الأداء، دون الحاجة إلى تعديل الخصائص الأساسية للمادة لكل اختبار.

المشكلة: الصلب "ذو الوجهين"

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

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

الحل الجديد: خريطة ثلاثية الأبعاد للأسهم المغناطيسية

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

إليك سحر نهجهما:

1. مجموعة واحدة من القواعد: استخدما مجموعة واحدة مشتركة من المعلمات الفيزيائية (مثل مغناطيسية التشبع) لكل من الاتجاهين الأمامي والجانبي. لم يحتاجا لتغيير "ارتفاع" مغناطيسية الصلب اعتماداً على الاتجاه الذي يدفعانه فيه.
2. رؤية "عمق القشرة": لاحظا أنه مع زيادة تردد المجال المغناطيسي (من 5ز هرتز إلى 200 هرتز)، يتغير السلوك. فأدخلا تعديلاً ذكياً يعتمد على التردد في الجزء "المكعب" (cubic) من مشهد الطاقة. هذا التعديل يتوافق تماماً مع فيزياء كيفية اختراق الكهرباء للمعدن ("عمق القشرة" - skin depth). يبدو الأمر كما لو أن المعدن يصبح "أكثر صلابة" قليلاً تجاه التغيرات المغناطيسية عند السرعات العالية، ونموذجهما يلتقط هذا بشكل طبيعي.
3. لغز "الركبتين": عند دفع الصلب جانبياً، يكون المنحنى المغناطيسي ذا شكل غريب يشبه "الركبتين" (ينحني مرتين قبل أن يستقر). لم يستطع نموذج الدوران الأساسي رؤية ذلك. لذا، أضاف المؤلفان "تراكب هيستيرون" (hysteron overlay) — وهي طبقة رياضية إضافية صغيرة تعمل مثل مجموعة من المفاتيح الصغيرة المنفصلة. هذه الطبوة تندمج تماماً في المكان الذي تظهر فيه القياسات الواقعية تبدل الجدران المغنية الداخلية للصلب.

ما وجداه

عندما اختبرا نموذجهما الجديد مقابل بيانات واقعية من صلب عالي الجودة يسمى Hi-B 27M-OH، كانت النتائج مبهرة.

• الدقة: حققا متوسط خطأ قدره 0.069 تسلا فقط عبر ثمانية اختبارات مختلفة (أمامي وجانبي عند أربعة ترددات مختلفة).
• الإصلاح: نجحا في إزالة أثر "الاعتماد على الاتجاه بنسبة 35%". وجد نموذجهما قيمة واحدة لمغناطيسية التشبع تبلغ 1.41×10⁶ أمبير متر/متر عملت لكلا الاتجاهين.
• المقايضة: للحصول على هذا التوحيد المثالي، اضطرا لقبول قيمة أقل قليلاً لـ "التباين المكعب" (وهو مقياس لمدى مقاومة البنية البلورية للدوران) عما هو موجود في بلورة واحدة مثالية. يوضح البحث أن هذا ليس خطأً؛ بل لأن الصلب يتكون من حبيبات صغيرة كثيرة، والنسيج المحدد للصلب "يخفي" بعض الاتجاهات البلورية، مما يقلل فعلياً من المقاومة الإجمالية التي يراها النموذج.

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

ما ليس عليه (وما هو التالي)

المؤلفان واضحان جداً بشأن ما لا يفعله نموذجهما بعد.

• ليس حلاً سحرياً لكل الزوايا: يعمل النموذج بشكل مثالي للاتجاهات الأمامية والجانبية الرئيسية، لكن التنبؤ بالسلوك عند زاوية 45 درجة لا يزال تحدياً مفتوحاً.
• ليس محاكاة فيزيائية كاملة للوحة بأكملة: عند الترددات العالية جداً (فوق 150 هرتز)، لا يخترق المجال المغناطيسي المعدن بشكل متساوٍ (تأثير القشرة). يتعامل النموذج مع هذا باستخدام اختصار رياضي ذكي، لكن المؤلفين يعترفان أنه للحصول على أعلى دقة، ستحتاج في النهاية إلى ربطه بمحلل كهرومغناطيسي كامل لرؤية المجال داخل شريحة المعدن طبقة تلو الأخرى.
• ليس مثالياً لكل التفاصيل الدقيقة: بينما يتعامل بشكل رائع مع "الحلقات الكبرى" الرئيسية للمغناطيسية، فإنه لا يستطيع بعد التقاط التأثيرات المعقدة للذاكرة الخاصة بـ "الحلقات الصغرى" (الحركات الصغيرة ذهاباً وإياباً) أو تأثيرات الانحياز المستمر (DC bias).

الخلاصة

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