Stability analysis of admittance control using asymmetric stiffness matrix

Titel: Stabiliteitsanalyse van admittantiecontrole met een asymmetrische stijfheidsmatrix

1. Probleemstelling

Bij robotica-taken die veel contact met de omgeving vereisen (zoals assemblage of "peg-in-hole"), is het instellen van de stijfheid van het regelsysteem cruciaal voor de prestaties. Traditionele methoden voor compliance-controle (zoals impedantie- en admittantiecontrole) gebruiken doorgaans een symmetrische stijfheidsmatrix.

• Beperking: Hoewel het gebruik van een asymmetrische stijfheidsmatrix de bereikbare bewegingsruimte kan vergroten en de robot meer flexibiliteit kan bieden bij complexe taken, is de stabiliteit van dergelijke systemen niet bewezen.
• Uitdaging: Asymmetrische matrices voldoen vaak niet aan de passiviteitseis (een veelgebruikte methode om stabiliteit te garanderen). Dit kan leiden tot ongewenste oscillaties, specifiek "spiraalvormige oscillaties" (spiral oscillations), waarbij energie in het systeem kan toenemen in plaats van te dissiperen. Er ontbreekt tot nu toe een theoretische onderbouwing voor het veilig instellen van parameters in dergelijke asymmetrische systemen.

2. Methodologie

De auteurs analyseren de convergentiestabiliteit van een planair admittantie-model voor een robotarm, waarbij wordt aangenomen dat de interne positie-regellus stabiel is.

• Modelopbouw:
  • Het admittantie-model wordt beschreven door: $F_{ext} = M\ddot{x} + D\dot{x} + Kx$.
  • De stijfheidsmatrix $K$ wordt opgesplitst in een symmetrisch deel ($K_s$) en een asymmetrisch deel ($K_a$).
  • Het asymmetrische deel creëert een "curl-veld" (krulveld) dat een rotatiekracht rond het evenwichtspunt genereert, wat niet-conservatief is en oscillaties kan opwekken.
• Stabiliteitsanalyse (Root Locus):
  • In plaats van puur te vertrouwen op passiviteit (wat moeilijk is bij asymmetrie), gebruiken de auteurs de root locus-methode (plaats van de nulpunten) om de eigenwaarden van het systeem te analyseren.
  • Ze leiden voorwaarden af voor de dempingscoëfficiënt ($d$) en de stijfheidsparameters zodat alle eigenwaarden van het systeemmatrix negatief reële delen hebben (stabiliteit).
• Kernvoorwaarde:
  • De auteurs tonen aan dat spiraalvormige oscillaties optreden als de eigenwaarden van de stijfheidsmatrix complexe getallen zijn (d.w.z. een imaginaire component hebben).
  • Ze leiden een ongelijkheid af die bepaalt wanneer de eigenwaarden reëel blijven. Als de eigenwaarden reëel zijn, kan kritische demping worden bereikt en worden spiraaloscillaties onderdrukt.
  • De afgeleide voorwaarde voor de asymmetrische component $k_a$ is:
    $$|k_a| < \sqrt{k_s^2 + \frac{1}{4}(k_x - k_y)^2}$$
    Waarbij $k_s$ de symmetrische kruiskoppeling is en $k_x, k_y$ de diagonale stijfheden.

3. Belangrijkste Bijdragen

1. Identificatie van Spiraaloscillaties: De auteurs hebben de exacte conditie geklaard waaronder asymmetrische stijfheidsmatrijzen spiraaloscillaties opwekken. Ze tonen aan dat het ontbreken van deze oscillaties equivalent is aan het hebben van uitsluitend reële eigenwaarden voor de stijfheidsmatrix.
2. Stabiliteitsgrenzen: Ze hebben een methode ontwikkeld om de stabiliteitsgrenzen voor de dempingsparameters af te leiden op basis van de eigenwaarden van het systeem, zelfs wanneer passiviteit niet kan worden gegarandeerd.
3. Parametervoorwaarden: Er zijn specifieke ongelijkheden afgeleid die aangeven hoe de demping ($d$) moet worden ingesteld in relatie tot de stijfheidsparameters om convergentie te garanderen.

4. Resultaten

De theorie werd gevalideerd via simulaties en experimenten:

• Simulaties:
  • Vergelijking tijdrespons: Simulaties toonden aan dat bij een asymmetrische matrix met complexe eigenwaarden (imaginaire componenten), het systeem oscilleert en spiraalvormige banen volgt, zelfs bij hoge demping.
  • Kritische demping: Wanneer de parameters voldeden aan de afgeleide voorwaarde (zodat eigenwaarden reëel waren), kon het systeem worden gedempt tot kritische demping zonder spiraaloscillaties.
  • Passiviteit: Hoewel het systeem lokaal niet passief was (energie kon tijdelijk toenemen bij lage snelheden door het curl-veld), bleek het systeem globaal stabiel en convergent zolang de afgeleide dempingsvoorwaarde werd gehaald.
• Experimenten:
  • Een robotarm werd gebruikt om contact te maken met een metalen plaat.
  • Er werd een vergelijking gemaakt tussen een symmetrische en een asymmetrische (driehoekige) stijfheidsmatrix.
  • Uitkomst: Bij gebruik van de asymmetrische matrix was de afwijking tussen de gewenste en daadwerkelijke positie kleiner (6 mm vs 8 mm bij symmetrisch). Bovendien werd de contactkrak hoger gehandhaafd omdat de beweging meer geconcentreerd was in de gewenste as, wat de efficiëntie van de assemblage-taak verbeterde.

5. Betekenis en Conclusie

Dit onderzoek is significant omdat het de theoretische barrière voor het gebruik van asymmetrische stijfheidsmatrijzen in admittantiecontrole wegneemt.

• Praktische toepassing: Robotontwerpers kunnen nu asymmetrische stijfheidsmatrices veilig gebruiken om de bewegingsvrijheid en aanpassingsvermogen van robots te vergroten, zolang ze de afgeleide dempingsvoorwaarden in acht nemen.
• Stabiliteit zonder passiviteit: Het paper toont aan dat stabiliteit kan worden gegarandeerd via eigenwaarde-analyse (root locus), zelfs als het systeem niet strikt passief is.
• Toekomstperspectief: De theorie is momenteel beperkt tot 2D-systemen, maar de auteurs suggereren dat deze kan worden uitgebreid naar 3D-systemen en variabele stijfheidscontrole, mogelijk met behulp van het concept van "energy tanks".

Kortom, dit paper biedt een wiskundig onderbouwde methode om de stabiliteit van robuuste, contactrijke robottaken te waarborgen door het slim en veilig gebruik van asymmetrische stijfheidsregeling.