Technische Samenvatting: Over het Ontstaan van een Pendelstructuur bij Meervoudige Contactlocomotie

Probleemstelling

Het Lineair Omgekeerde Pendelmodel (LIPM) is een alomtegenwoordige aanname in de besturing van benen locomotie, vaak opgelegd als een modelbeperking in plaats van dat het natuurlijk voortvloeit uit het doel van de regelaar. Hoewel er extensies bestaan voor variabele hoogte of hoekimpuls, wordt de pendelvorm doorgaans "met de hand" toegevoegd. Dit artikel onderzoekt een smallere vraag: Welke specifieke doelstellingsfunctie zorgt ervoor dat het pendelkrachtenpatroon ontstaat als de optimale oplossing, en waar breekt dit beeld af?

De auteurs analyseren een klein centroidaal Optimaal Besturingsprobleem (OCP) dat de veranderingssnelheid van hoekimpuls ($\|\dot{H}_G\|$) bestraft, naast standaardtermen voor contactkrachten en normale versnelling. Het doel is om te bepalen of het minimaliseren van de hoekimpuls-snelheid het systeem van nature naar de pendelvariëteit stuurt (waar $\dot{H}_G = 0$) en om de grenzen van dit gedrag te karakteriseren onder verschillende standconfiguraties (volledige rang vs. $N=2$) en taakeisen.

Methodologie

Het artikel formuleert een Net-Wrench OCP (Definitie 1) dat een kostfunctie minimaliseert bestaande uit:

1. Straf voor Hoekimpuls: $\alpha \|\dot{H}_G\|^2$
2. Taakvolging: $\lambda \|\dot{H}_G - \dot{H}_{task}\|^2$ (voor niet-balans taken)
3. Straf voor Normale Versnelling: $\beta (n^\top \ddot{c})^2$
4. Krachtregularisatie: $\gamma \sum \|f_i\|^2$

De analyse verloopt in drie theoretische fasen, gevalideerd door numerieke experimenten op een puntmassa-vierpoot en de Unitree Go1-robot in MuJoCo (zowel open-lus QP als gesloten-lus koppel-niveau MPC):

1. Analyse van Volledige Rang Stand: Met behulp van Singuliere Waarde Decompositie (SVD) van de moment-Jacobaan leiden de auteurs het asymptotische gedrag af naarmate het gewicht voor hoekimpuls $\alpha \to \infty$.
2. Analyse van $N=2$ Stand: Onderzoek van het rang-deficiënte geval (bijv. diagonale galop), waarbij de momentoperator rang 1 heeft, wat leidt tot onvermijdelijke residuen en wrijvingskegelbeperkingen.
3. Taak-gedreven Analyse: Onderzoek naar hoe het toevoegen van een taakvolgingsterm ($\lambda > 0$) het optimum verschuift weg van de pendelvariëteit.

Belangrijkste Bijdragen en Theoretische Resultaten

1. Ontstaan van Pendelstructuur in Volledige Rang Stand

• Stelling 1: Voor standen met volledige rang stort de optimale oplossing asymptotisch in naar het pendelkrachtenpatroon ($\dot{H}_G \to 0$) met een snelheid van $O(1/\alpha)$.
• Geometrische Constante: De evenredigheidsconstante wordt bepaald door de SVD van de moment-Jacobaan, wat direct overeenkomt met de voet-span geometrie.
• Verbinding met Klassieke Modellen: Op deze pendelvariëteit is het Nul Moment Punt (ZMP) natuurlijk gelijk aan het virtuele draaipunt, en komt de divergente component van beweging (DCM) orbitale vergelijking voort als een gevolg, in plaats van een aanname.

2. Fundamentele Grenzen in $N=2$ Stand (Diagonale Galop)

• Stelling 5 (Geometrische Vloer): In $N=2$ standen is de moment-Jacobaan rang-1. Een residu van hoekimpuls-snelheid blijft bestaan langs de richting van de lijn die de twee voeten verbindt ($\hat{D}$), ongeacht krachthervorming. Dit creëert een niet-nul ondergrens (vloer) voor $\|\dot{H}_G\|$.
• Stelling 6 (Knik in Wrijvingskegel): Een tweede belemmering ontstaat door de wrijvingskegel. Als de kracht die nodig is om de loodrechte momentcomponent te annuleren buiten de wrijvingskegel valt, treedt er een "niet-gladde knik" op bij een kritieke horizontale versnelling ($a_x^*$). Voorbij dit punt overschrijdt de bereikbare $\|\dot{H}_G\|$ strikt de geometrische vloer.
• Implicatie: In diagonale galop kantelt de pendelrichting vaak buiten de wrijvingskegel (bijv. $40.6^\circ$ versus $31^\circ$ voor $\mu=0.6$), waardoor de pendeloplossing onhaalbaar wordt en een $\alpha$-onafhankelijke vloer ontstaat.

3. Taak-gedreven Afwijking

• Stelling 7: Wanneer een taakterm ($\lambda > 0$) wordt toegevoegd, wordt de optimale hoekimpuls-snelheid een gladde interpolatie: $\dot{H}_G^* = \frac{\lambda}{\alpha + \lambda} \dot{H}_{task}^*$.
• Fasediagram: Het systeem gaat soepel over van een "balans-gedomineerd" regime (pendel, $\lambda/\alpha \ll 1$) naar een "taak-gedomineerd" regime (niet-pendel, $\lambda/\alpha \gg 1$). Dit verklaart waarom dezelfde robot LIPM-achtig kan lijken in statische balans maar niet-pendel tijdens wendbare manoeuvres zonder dat de regelaar wordt gewisseld.

Experimentele Resultaten

De auteurs hebben deze claims gevalideerd met de Unitree Go1 en een puntmassamodel:

• Snelheidsverificatie (Puntmassa): Naarmate $\alpha$ toenam, daalde de afwijking van het pendelkrachtenpatroon met een factor 227, met een log-log helling van $-1$, wat de $O(1/\alpha)$ snelheid bevestigt (Stelling 1).
• Geometrische Constante (Go1 N=4): In een vier-voet stand kwam de empirische schaalconstante ($K_e = 8.4$) overeen met de analytische voorspelling afgeleid van voetspans ($K_a = 9.8$) binnen $\sim 16\%$. De singuliere waarden van de moment-Jacobaan kwamen overeen met de fysieke voetspans tot vier decimalen.
• Vloerverificatie (Go1 N=2): In diagonale galop zakte $\|\dot{H}_G\|$ uit op een voorspelde geometrische vloer van $0.2835 \, \text{m}^2/\text{s}^2$, wat overeenkwam met de analytische ondergrens tot vier cijfers.
• Gesloten-lus Realiteit: In een gesloten-lus koppel-niveau MPC bleef de vloer bestaan maar was deze hoger ($\approx 1.71 \, \text{m}^2/\text{s}^2$) dan de statische voorspelling. De auteurs schrijven dit gat toe aan dynamische effecten (zwaai-impacten, contactovergangen, servo-fouten) die niet worden vastgelegd in de quasi-statische analyse.
• ZMP Identiteit: De numerieke experimenten bevestigden dat op de pendelvariëteit de op kracht gebaseerde ZMP convergeert naar het virtuele draaipunt, wat de theoretische link tussen LIPM en het OCP-optimum valideert.

Betekenis en Aanspraken

Het artikel claimt de link tussen LIPM en optimale besturing expliciet en controleerbaar te maken. Het betoogt dat:

1. LIPM een gevolg is, niet alleen een keuze: Het pendelpatroon ontstaat natuurlijk bij het minimaliseren van de hoekimpuls-snelheid in standen met volledige rang, mits het gewicht $\alpha$ voldoende hoog is.
2. Grenzen voorspelbaar zijn: Het afbreken van de pendelaanname in $N=2$ standen is geen falen van de regelaar maar een fundamentele geometrische en wrijvingsbeperking. Het artikel biedt gesloten-vorm expressies voor de ondergrenzen en de kritieke versnelling waarbij wrijvingsbeperkingen binden.
3. Gefuseerd perspectief op Balans en Taak: De verhouding $\lambda/\alpha$ biedt een continue parameter die de overgang tussen balans-gecentreerd (LIPM-achtig) en taak-gecentreerd (wendbaar) gedrag verklaart.

De auteurs blijven bescheiden wat betreft de gesloten-lus resultaten, en merken op dat hoewel het kwalitatieve gedrag (het bestaan van een vloer) geldt in dynamische simulaties, het kwantitatieve gat tussen statische theorie en gesloten-lus realiteit significant is en waarschijnlijk te wijten is aan niet-gemodelleerde dynamica (compliantie, slip, eindige besturingsfrequenties). Zij claimen niet deze dynamische gaten op te lossen, maar eerder de statische grenzen te definiëren waarbinnen die dynamica opereert.