Adaptive Learning of Periodic Solutions for Nonlinear Systems by Physics-informed Gaussian Process
تقترح هذه الورقة إطار عمل جديد للعملية الغاوسية المستنيرة بالفيزياء (PIGP) يستخدم نوى دورية مخصصة واستراتيجية نقاط تجميع تكيفية لحساب الحلول الدورية للأنظمة الديناميكية غير الخطية بكفاءة ودقة، مما يظهر أداءً حوسبيًا متفوقًا مقارنة بطريقة التوازن التوافقي.
المؤلفون الأصليون:Yixin Li, Zhong-Rong Lu, Dahao Yang, Zechang Zheng, Jike Liu, Li Wang
في عالم الهندسة، لا تسلك العديد من الهياكل مسارات بسيطة ومستقيمة. فعندما يتأرجح جسر بفعل الرياح، أو يدور شفرة توربين بسرعة عالية، أو يهتز مبنى أثناء زلزال، فإن القوى المعنية غالبًا ما تتغير بطرق معقدة اعتمادًا على مقدار حركة الجسم أو سرعته. تُعرف هذه بـ "الأنظمة غير الخطية". وخلافًا للنابض البسيط الذي يتمدد بنفس المقدار مهما سحبته بقوة، يمكن لهذه الأنظمة أن تلتوي، أو تنكسر، أو تستقر في أنماط غريبة ومتكررة يصعب التنبؤ بها. ويعد العثور على الحركة المحددة والمتكررة — التي تسمى الحل الدوري — التي سيستقر فيها النظام أمرًا بالغ الأهمية للمهندسين؛ فإذا لم يتمكنوا من التنبؤ بهذه الحركات، فلن يتمكنوا من تصميم هياكل آمنة أو فعالة. ولعقود من الزمن، اعتمد العلماء على أدوات رياضية راسخة لإيجاد هذه الأنماط المتكررة، لكن هذه الأدوات غالبًا ما تعاني عندما يصبح النظام شديد التعقيد أو "متعرجًا" للغاية في سلوكه، مما يتطلب قدرة حوسبية هائلة للوصول إلى إجابة دقيقة.
لقد طور فريق من الباحثين من جامعة سون يات سين وجامعة لييج نهجًا جديدًا لحل هذه المشكلة، وهو نهج يمزج بين قوانين الفيزياء ونوع من تعلم الآلة يسمى "العملية الغاوسية". وبدلاً من حشر النظام في صندوق رياضي جامد، يعامل منهجهم، الذي يسمونه "العملية الغاوسية المستنيرة بالفيزياء"، القوانين الفيزيائية التي تحكم النظام كمجموعة من القرائن. تخيل محاولة تخمين شكل جسم مخفي عن طريق تحسس سطحه في نقاط محددة قلي تعداد؛ فكلما زادت النقاط التي تتحسسها، أصبح الشكل أكثر وضوحًا. في هذا المنهج الجديد، يستخدم الباحثون معادلات الحركة الفيزيائية كـ "أحاسيس" في لحظات زمنية محددة. يبدأون بتخمين تقريبي لكيفية حركة النظام، ثم يستخدمون عملية ذكية وتكرارية لتحسين هذا التخمين. يتعلم النظام من أخطائه، ويعدل فهمه للحركة حتى يتطابق تمامًا مع القواعد الفيزيائية.
ما يجعل هذا النهج قويًا بشكل خاص هو قدرته على التعلم أين يحتاج إلى البحث. فقد اكتشف الباحثون أن ليست كل اللحظات الزمنية متساوية في الأهمية لفهم الحركة. ففي الأنظمة التي تكون فيها الحركة سلسة، تتوزع المعلومات بشكل متساوٍ. أما في الأنظمة التي تكون فيها الحركة خشنة أو مفاجئة — مثل جزء يصطدم بحاجز أو مادة تتغير صلابتها فجأة — فإن المعلومات الأكثر أهمية تتركز في لحظات زمنية صغيرة ومحددة. يكتشف المنهج الجديد تلقائيًا هذه اللحظات الصعبة ويركز جهده الحوسبي هناك، حيث يضيف المزيد من "نقاط الاستشعار" بالضبط حيث تكون درجة عدم اليقين في أعلى مستوياتها. تسمح هذه الاستراتيجية التكيفية للحاسوب بإيجاد الحل باستخدام عدد أقل بكثير من نقاط البيانات مقارنة بالطرق التقليدية، مما يجعل العملية أسرع وأكثر كفاءة بشكل ملحوظ.
اختبر الباحثون منهجهم الجديد على أنواع مختلفة من الأنظمة، بدءًا من المذبذبات البسيطة ذات الكتلة الواحدة وصولًا إلى نماذج المباني المعقدة متعددة الطوابق. وقارنوا نتائجهم بتقنية قياسية مستخدمة على نطاق واسع تُعرف باسم "طريقة التوازن التوافقي"، والتي كانت الأداة الأساسية لعقود. وفي كل اختبار، أثبت المنهج الجديد دقة عالية، وقدرة على التعامل مع الحركات السلسة وتلك التغيرات المفاجئة والمتعرجة. ولعل الأمر الأكثر إثارة للإعجاب هو أنه عندما طلب الباحثون من كلا الطريقتين تحقيق نفس المستوى من الدقة، أتم المنهج المستنير بالفيزياء المهمة باستمرار في وقت أقل. وقد جاء هذا التفوق في السرعة من قدرة المنهج على التكيف؛ فهو لا يضيع الوقت في حساب التفاصيل حيث لا داعي لها، بل يوجه موارده إلى الأجزاء الأكثر صعوبة في الحل.
وبعيدًا عن السرعة والدقة، يقدم المنهج الجديد ميزة فريدة تفتقر إليها الأدوات التقليدية: وهي توفير مقياس داخلي للثقة. فعندما يحسب المنهج حلاً، فإنه ينتج أيضًا خريطة لعدم اليقين، توضح بدقة مدى تأكده من النتيجة في كل لحظة زمنية. هذا يشبه وجود توقعات جوية لا تخبرك بدرجة الحرارة فحسب، بل تعطيك أيضًا نطاقًا واضحًا لمدى تباين تلك الدرجة. وبالنسبة للمهندسين الذين يصممون بنية تحتية حيوية، فإن معرفة حدود حساباتهم لا تقل أهمية عن الحسابات نفسها. ووجد الباحثون أنه مع إضافة المزيد من نقاط الاستشعار إلى النظام، تتقلص خريطة عدم اليقين هذه، مما يؤكد أن الحل أصبح أكثر موثوقية.
بينما يظهر المنهج إمكانات كبيرة، فقد كان الباحثون حذرين في تحديد حدوده. فقد وجدوا أنه بالنسبة للأنظمة ذات التغيرات الحادة والمفاجئة للغاية، يمكن للحل أن يتذبذب أحيانًا قليلاً بالقرب من تلك النقاط الحادة، وهو تحدٍ يواجهه حتى المنهج التقليدي. ومع ذلك، من خلال إثبات أن نهجهم يمكنه التعامل مع مجموعة متنوعة من المشكلات الصعبة بكفاءة أكبر وموثوقية داخلية، فتح الفريق مسارًا جديدًا لتحليل الأنظمة الميكانيكية المعقدة. ويشير عملهم إلى أنه من خلال السماح لتعلم الآلة بالتعلم مباشرة من قوانين الفيزياء، وليس فقط من البيانات، يمكن للمهندسين حل بعض أكثر المشكلات استعصاءً في الديناميكا بوضوح وسرعة كانا بعيدي المنال في السابق.
ملخص تقني: التعلم التكيفي للحلول الدورية للأنظمة غير الخطية بواسطة العمليات الغاوسية الموجهة بالفيزياء
بيان المشكلة تنتشر الأنظمة الديناميكية غير الخطية في الهندسة، حيث تُظهر سلوكيات معقدة مثل عدم الاستقرار، والتفرع (bifurcations)، والحلول المتعددة المتزامنة. وتتمثل المهمة المركزية في تحليل هذه الأنظمة في الحصول على حلولها الدورية. تواجه الطرق الحالية قيوداً كبيرة:
الطرق التحليلية (مثل طرق الاضطراب، والمقاييس المتعددة) تقتصر عموماً على الأنظمة ضعيفة اللاخطية.
الطرق العددية مثل طريقة "رونج-كوتا" (RK) تعاني مع الحلول غير المستقرة وتتطلب أوقات تكامل طويلة للغاية للوصول إلى استجابات الحالة المستقرة. أما "طريقة الرماية" (shooting method)، ورغم قدرتها على التقاط الحلول غير المستقرة، إلا أنها تتطلب غالباً عمليات تكامل زمني متكررة ومكثفة.
طريقة التوازن التوافقي (HBM)، وهي نهج شبه تحليلي واسع الاستخدام، تقرب الحلول عبر متسلسلات فوريه. وبينما تعد فعالة للأنظمة السلسة، فإن معدل تقاربها يتدهور بشكل ملحو-ظور في الأنظمة غير السلسة، مما يتطلب عدداً كبيراً من الحدود التوافقية ويؤدي إلى تكاليف حوسبة عالية. علاوة على ذلك، تفتقر طريقة HBM إلى معايير صريحة لاختيار الرتبة التوافقية المثلى ولا توفر بطبيعتها تقديراً لعدم اليقين.
المنهجية يقترح البحث إطار عمل للعمليات الغاوسية الموجهة بالفيزياء (PIGP) للحساب التكيفي للحلول الدورية. تدمج المنهجية بين تعلم الآلة البايزي والمعادلات الفيزيائية الحاكمة من خلال ثلاثة مكونات رئيسية:
الخطيّة والحل التكراري: بما أن العمليات الغاوسية (GP) مهيأة طبيعياً للمعادلات الخطية، يتم خطيّة المعادلة الحاكمة غير الخطية عند كل تكرار باستخدام توسيع تايلور من الدرجة الأولى (تكرار نيوتن). يحول هذا العملية غير الخطية إلى سلسلة من المعادلات الديناميكية الخطية المتغيرة زمنياً.
صياغة العملية الغاوسية الموجهة بالفيزياء (PIGP):
النواة والدوال المتوسطة (Kernel and Mean Functions): يتم بناء دالة نواة دورية مخصصة لفرض دورية الحل. بالإضافة إلى ذلك، يتم إدخال دالة متوسطة تتضمن حدود فوريه للتردد الأساسي لتسهيل الاستخراج المباشر لمعاملات فوريه.
البيانات الرصدية: تُعامل المعادلات الحاكمة الخطية، عند تقييمها في مجموعة من نقاط التجميع (collocation points)، كبيانات رصدية.
الاستدلال البايزي: يتم إنشاء توزيع غاوسي مشترك بين الحل والبيانات الرصدية. يوفر المتوسط البعدي (posterior mean) الحل المحدث، بينما يقيس التباين البعدي (posterior covariance) عدم اليقين في الحل. ويتم تحسين المعلمات الفائقة (hyperparameters) ومعلمات الدالة المتوسطة عبر تقدير الاحتمال الأقصى (MLE) في كل تكرار.
استراتيجية التجميع التكيفية: لتعزيز الكفاءة والدقة، يطور البحث استراتيجية تكيفية لاختيار نقاط التجميع. يتم إضافة نقاط جديدة تدريجياً إلى المجموعة بناءً على مؤشرات الخطأ:
المعتمد على التباين: اختيار النقاط التي يتم فيها تعظيم أثر مصفوفة التباين البعدي.
المعتمد على البواقي (Residual-based): اختيار النقاط التي يتم فيها تعظيم مربع الباقي المتوقع للمعادلة الحاكمة. وقد وجدت الدراسة أن المعيار المعتمد على البواقي يحقق باستمرار دقة أعلى وتقارباً أكثر استقراراً.
المساهمات الرئيسية
إطار عمل مبتكر: تطوير PIGP مصمم خصيصاً للحلول الدورية للأنظمة الديناميكية غير الخطية العامة، متجاوزاً بذلك القصور الحالي في تطبيقات PIGP المقتصرة على المسائل الخطية.
تقدير عدم اليقين: على عكس الطرق العددية التقليدية، توفر PIGP بطبيعتها تباينًا بعديًا، مما يقدم مقياساً موثوقاً وصريحاً لدقة الحل وعدم اليقين.
الكفاءة التكيفية: تقديم استراتيجية تجميع تكيفية تركز الموارد الحوسبية على مناطق الخطأ العالي (مثل النقاط غير السلسة)، مما يحسن الكفاءة بشكل كبير مقارنة بأخذ العينات الموحد.
الأداء المقارن: يوضح البحث أن PIGP تحقق كفاءة حوسبية متفوقة مقارنة بطريقة HBM لمستوى معين من الدقة، لا سيما في الأنظمة غير السلسة.
النتائج أُجريت تجارب عددية على أربعة أنظمة: متذبذب "دوفينج" (Duffing oscillator) سلس ذو درجة حرية واحدة (SDOF)، ونظام SDOF غير سلس يحتوي على حلول غير مستقرة ومضاعفة للدورة (period-doubling)، ونظام مقترن سلس ذو درجتي حرية، ونموذج هيكل قص ثلاثي الطوابق مع قوى استعادة غير سلسة (صلابة ثنائية الخطية وتخميد لزوجي غير خطي).
الدقة والتقارب: أعادت PIGP إنتاج الحلول المرجعية (التي تم الحصول عليها عبر طريقة الرماية) بدقة لكل من الأنظمة السلسة وغير السلسة، والمستقرة وغير المستقرة. انخفض خطأ L2 النسبي والتباين البعدي باستمرار مع زيادة عدد نقاط التجميع.
السلوك التكيفي: بالنسبة للأنظمة غير السلسة، نجحت الاستراتيجية التكيفية في تركيز نقاط التجميع حول السمات غير السلسة (مثل حيث تكون الإزاحة أو السرعة صفراً)، بينما ظلت النقاط موزعة بانتظام في الأنظمة السلسة.
المقارنة مع HBM:
التقارب: بينما أظهرت طريقة HBM خطأً أقل قليلاً لنفس عدد المجاهيل (بسبب مثالية قطع فوريه في معيار L2-norm)، أظهر كلا الطريقتين معدلات تقارب نوعية متشابهة (أسية للأنظمة السلسة، ومتعددة الحدود للأنظمة غير السلسة).
الكفاءة: كانت PIGP أكثر كفاءة حوسبية باستمرار من HBM. وعبر مستويات دقة مختلفة، تطلبت PIG_P وقتاً حوسبياً أقل بكثير (غالباً أقل من 60% من وقت HBM) لتحقيق نفس مستوى الدقة.
الأهمية والادعاءات يدعي البحث أن الـ PIGP المقترحة تشكل أداة قوية وفعالة قائمة على البيانات لحل الأنظمة الدورية غير الخطية. وتكمن أهميتها الأساسية في:
التعامل مع اللاخطية: معالجة كل من اللاخطية السلسة وغير السلسة بنجاح دون الحاجة إلى مجموعة محددة مسبقاً من دوال الأساس.
المتانة: توفير الحلول للمدارات الدورية غير المستقرة والأنظمة ذات التفرعات.
الكفاءة: التفوق على طريقة HBM الراسخة في التكلفة الحوسبية مقابل الدقة المتكافئة، ويرجع ذلك أساساً إلى الاختيار التكيفي لنقاط التجميع.
الموثوقية: تقديم تقدير متأصل لعدم اليقين، والذي يعمل كمؤشر خطأ مدمج، مما يعالج قصوراً رئيسياً في الطرق التقليدية.
يشير المؤلفون بتواضع إلى وجود قيد: بالنسبة للمسائل شديدة عدم السلاسة، قد يظهر الحل تذبذبات بالقرب من النقاط غير السلسة، وهو تحدٍ تشترك فيه طريقة HBM. ويقترحون أن العمل المستقبلي يمكن أن يتضمن تصميم دوال متوسط ونواة مخصصة مع احتمالات (priors) غير سلسة للتخفيف من هذه المشكلة بشكل أكبر.