ملخص تقني: حول ظهور البنية البندولية في الحركة متعددة الاتصالات

بيان المشكلة

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

يحلل المؤلفون مسألة تحسين مركزية (OCP) صغيرة لمركز الكتلة، تعاقب معدل تغير الزخم الزاوي ($\|\dot{H}_G\|$) إلى جانب حدود قياسية لقوى الاتصال والتسارع العمودي. الهدف هو تحديد ما إذا كان تقليل معدل الزخم الزاوي يدفع النظام طبيعياً نحو المنطوي البندولي (حيث $\dot{H}_G = 0$) وتوصيف حدود هذا السلوك تحت تكوينات الوقوف المختلفة (كامل الرتبة مقابل $N=2$) ومتطلبات المهام.

المنهجية

يصيغ البحث مسألة تحسين صافي العزم (Net-Wrench OCP) (التعريف 1) التي تقلل دالة تكلفة تتكون من:

1. عقوبة الزخم الزاوي: $\alpha \|\dot{H}_G\|^2$
2. تتبع المهمة: $\lambda \|\dot{H}_G - \dot{H}_{task}\|^2$ (للمهام غير المتعلقة بالتوازن)
3. عقوبة التسارع العمودي: $\beta (n^\top \ddot{c})^2$
4. تنظيم القوة: $\gamma \sum \|f_i\|^2$

يتقدم التحليل عبر ثلاث مراحل نظرية، تم التحقق منها تجريبياً على نموذج كتلة نقطية وروبوت Unitree Go1 في بيئة MuJoCo (كلاهما في تحسين تربيعي مفتوح الحلقة وMPC للتحكم في عزم الدوران في الحلقة المغلقة):

1. تحليل الوقوف كامل الرتبة: باستخدام تفكيك القيم المفردة (SVD) لمصفوفة جاكوبي العزم، يستنتج المؤلفون السلوك التقاربي مع زيادة وزن الزخم الزاوي $\alpha \to \infty$.
2. تحليل الوقوف $N=2$: استقصاء الحالة ناقصة الرتبة (مثل مشية التروت القطري)، حيث تكون مؤثرات العزم برتبة 1، مما يؤدي إلى بقايا لا يمكن تجنبها وقيود مخروط الاحتكاك.
3. التحليل المدفوع بالمهمة: فحص كيفية تأثير إضافة حد تتبع المهمة ($\lambda > 0$) في إزاحة الحل بعيداً عن المنطوي البندولي.

المساهمات الرئيسية والنتائج النظرية

1. ظهور البنية البندولية في الوقوف كامل الرتبة

• النظرية 1: بالنسبة لحالات الوقوف كاملة الرتبة، ينهار الحل الأمثل تقاربياً نحو نمط القوة البندولية ($\dot{H}_G \to 0$) بمعدل $O(1/\alpha)$.
• الثابت الهندسي: يتم تحديد ثابت التناسب من خلال SVD لمصفوفة جاكوبي العزم، والذي يتوافق مباشرة مع هندسة اتساع القدمين.
• الارتباط بالنماذج الكلاسيكية: على هذا المنطوي البندولي، يتساوى نقطة الصفر للعزم (ZMP) طبيعياً مع نقطة الارتكاز الافتراضية، وتنبثق معادلة مدار مكون التباعد للحركة (DCM) كنتيجة وليس كافتراض.

2. الحدود الأساسية في وقفة $N=2$ (التروت القطري)

• النظرية 5 (الأرضية الهندسية): في حالات الوقوف $N=2$، تكون مصفوفة جاكوبي العزم برتبة 1. يظل هناك معدل زخم زاوي متبقٍ على طول اتجاه الخط الواصل بين القدمين ($\hat{D}$)، بغض النظر عن إعادة توزيع القوى. هذا يخلق حداً أدنى غير صفري لـ $\|\dot{H}_G\|$.
• النظرية 6 (انثناء مخروط الاحتكاك): تظهر عقبة ثانية من مخروط الاحتكاك. إذا كانت القوة المطلوبة لإلغاء مكون العزم العمودي تقع خارج مخروط الاحتكاك، يحدث "انثناء غير سلس" عند تسارع أفقي حرج ($a_x^*$). بعد هذه النقطة، يتجاوز $\|\dot{H}_G)$ المحقق بشكل صارم الأرضية الهندسية.
• الآثار المترتبة: في مشية التروت القطري، يميل الاتجاه البندولي غالباً خارج مخروط الاحتكاك (مثلاً $40.6^\circ$ مقابل $31^\circ$ لمعامل احتكاك $\mu=0.6$)، مما يجعل الحل البندولي غير ممكن ويخلق أرضية مستقلة عن $\alpha$.

3. الانحراف المدفوع بالمهمة

• النظرية 7: عند إضافة حد للمهمة ($\lambda > 0$)، يصبح معدل الزخم الزاوي الأمثل عبارة عن استكمال سلس: $\dot{H}_G^* = \frac{\lambda}{\alpha + \lambda} \dot{H}_{task}^*$.
• مخطط الطور: ينتقل النظام بسلاسة من نظام "يهيمن عليه التوازن" (بندولي، $\lambda/\alpha \ll 1$) إلى نظام "تهيمن عليه المهمة" (غير بندولي، $\lambda/\alpha \gg 1$). وهذا يفسر لماذا يمكن لنفس الروبوت أن يبدو شبيهاً بـ LIPM في التوازن الساكن ولكنه غير بندولي أثناء المناورات الرشيقة دون تغيير المتحكم.

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

تحقق المؤلفون من هذه الادعاءات باستخدام Unitree Go1 ونموذج كتلة نقطية:

• التحقق من المعدل (كتلة نقطية): مع زيادة $\alpha$، انخفض الانحراف عن نمط القوة البندولية بمقدار 227 ضعفاً، مع ميل لوغاريتمي-لوغاريتمي قدره $-1$، مما يؤكد معدل$O(1/\alpha)$ (النظرية 1).
• الثابت الهندسي (Go1 N=4): في حالة الوقوف بأربع قدمين، طابق الثابت التجريبي ($K_e = 8.4$) التنبؤ التحليلي المستمد من اتساع القدمين ($K_a = 9.8$) بنسبة خطأ $\sim 16\%$. كما طابقت القيم المفردة لمصفوفة جاكوبي العزم اتساع القدمين الفيزيائي حتى أربعة أرقام عشرية.
• التحقق من الأرضية (Go1 N=2): في مشية التروت القطري، تشبعت $\|\dot{H}_G\|$ عند أرضية هندسية متوقعة قدرها $0.2835 \, \text{m}^2/\text{s}^2$، مطابقةً الحد التحليلي حتى أربعة أرقام.
• واقع الحلقة المغلقة: في نظام MPC للتحكم في عزم الدوران في الحلقة المغلقة، استمرت الأرضية ولكن بقيمة أعلى ($\approx 1.71 \, \text{m}^2/\text{s}^2$) من التنبؤ الساكن. يعزو المؤلفون هذه الفجوة إلى التأثيرات الديناميكية (اصطدامات الأرجحة، انتقالات الاتصال، أخطاء السيرفو) التي لا تشملها التحليلات شبه الساكنة.
• تطابق ZMP: أكدت التجارب العددية أنه في المنطوي البندولي، تتقارب ZMP القائمة على القوة مع نقطة الارتكاز الافتراضية، مما يثبت الارتباط النظري بين LIPM والحل الأمثل لـ OCP.

الأهمية والادعاءات

يزعم البحث جعل الربط بين LIPM والتحكم الأمثل صريحاً وقابلاً للتحقق. ويجادل بأن:

1. LIPM هو نتيجة وليس مجرد خيار: يظهر النمط البندولي طبيعياً عند تقليل معدل الزخم الزاوي في حالات الوقوف كاملة الرتبة، بشرط أن يكون الوزن $\alpha$ مرتفعاً بما يكفي.
2. الحدود قابلة للتنبؤ: إن انهيار افتراض البندول في حالات الوقوف $N=2$ ليس فشلاً للمتحكم بل هو قيد هندسي واحتكاكي أساسي. يوفر البحث تعبيرات مغلقة للأرضيات والحد الحرج للتسارع حيث تتقيد قيود الاحتكاك.
3. رؤية موحدة للتوازن والمهمة: يوفر النسبة $\lambda/\alpha$ معاملاً مستمراً يشرح الانتقال بين السلوك المتمحور حول التوازن (الشبه بـ LIPM) والسلوك المتمحور حول المهمة (الرشيق).

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