Maximizing Task Endurance through Muscle Fatigue Minimization with Consideration of Muscle Fatigability and Task Contribution: an in-Silico FES Study of Handcycling

Technische Samenvatting: Maximaliseren van Taakduur door Spiervermoeidheid te Minimaliseren bij FES-gestuurde Handfietsen

Probleemstelling
Functionele Elektrische Stimulatie (FES) is een revalidatietechniek die wordt gebruikt om de motorische functie te herstellen na neuromusculoskeletale aandoeningen. De klinische bruikbaarheid hiervan wordt echter aanzienlijk gehinderd door het voortijdige optreden van perifere spiervermoeidheid. In tegenstelling tot natuurlijke rekrutering, die bij uithoudingstaken bij voorkeur langzame vezels activeert, activeert FES synchroon motorische eenheden onder de elektroden met beperkte vezeltype-selectiviteit, wat leidt tot snelle vermoeidheid. Hoewel er diverse strategieën bestaan om vermoeidheid op spierniveau te mitigeren (bijv. het moduleren van frequentie of pulswijdte), houden deze vaak geen rekening met de mechanische redundantie van het bewegingsapparaat. Huidige optimale controlemethoden voor FES behandelen vermoeidheid vaak als een verstoring in plaats van een toestand die geminimaliseerd moet worden, en bestaande kostenfuncties (meestal gericht op het minimaliseren van elektrische lading of pulswijdte) zijn niet rigoureus vergeleken op hun vermogen om de taakduur te verlengen. Deze studie adresseert de kloof in het identificeren van kostenfuncties en weegstrategieën die de taakduur maximaliseren in een redundant systeem van de bovenste ledematen.

Methodologie
De studie maakte gebruik van een in-silico benadering met een musculoskeletaal model van 2 vrijheidsgraden en 4 spieren van de rechter bovenarm (anterieure/posterieure deltoïdeus, biceps, triceps), aangepast voor een handfiets-taak bij 60 RPM. De systeemdynamica combineerde een musculoskeletaal model met de fenomenologische FES-pulswijdte- en vermoeidheidsmodellen van Ding et al. (2003, 2007).

De kern van de methodologie bestond uit een Receding Horizon Optimal Control (RHOC) schema. Vanwege de computationele onhandelbaarheid van het oplossen van een enkele optimalisatie van 20 minuten, werd het probleem opgedeeld in vensters van 2 tot 5 krukcycli. De studie verliep in drie afzonderlijke fasen:

1. Vergelijking van Kostenfuncties: Zeventien kostenfuncties werden geëvalueerd, waarbij vier genormeerde archetypen (gemiddelde, kwadratisch, cubisch, oneindig) werden gecombineerd met vijf fysiologische/stimulatievariabelen (pulswijdte, spierkracht, stress, vermogen en de vermoeidheidstoestand $a'$). Deze werden vergeleken met een standaard klinisch-achtige stimulatiepatroon.
2. Bayesiaanse Optimalisatie: Een Bayesiaanse optimalisatiebenadering werd gebruikt als een externe loop black-box methode om spierspecifieke gewichten binnen een kwadratische vermoeidheidskostenfunctie af te stemmen. Dit had als doel het aantal voltooide cycli voordat taalfalen optrad te maximaliseren, waarbij de zoekruimte van gewichten voor de vier spieren werd verkend.
3. Afleiding van Fysiologische Gewichten: Om de hoge computationele kosten van Bayesiaanse optimalisatie (135 uur) aan te pakken, werd een originele, interpreteerbare methode voorgesteld. Deze methode berekende spiergewichten op basis van twee metrieken: spiervermoeibaarheid (snelheid van krachtverlies onder hoge belasting) en mechanische kritikaliteit (unieke koppelbijdrage en ondersteuning vóór negatieve koppelzones).

Belangrijkste Bijdragen

• Systematische Vergelijking van Kostenfuncties: De studie biedt een uitgebreide evaluatie van 17 variaties in kostenfuncties, waarbij wordt aangetoond dat de keuze van de doelvariabele een significant grotere impact heeft op de duur dan het wiskundige archetype van de norm.
• Gewogen Optimalisatiestrategie: Het introduceert een kader voor het integreren van spierspecifieke gewichten in FES-optimale controle, waarbij verder wordt gegaan dan uniforme kostenfuncties om rekening te houden met fysiologische verschillen (vermoeibaarheid) en mechanische rollen.
• Interpreteerbare Wegingsheuristiek: Het papier stelt een snelle (<1 seconde) berekeningsmethode voor voor spiergewichten op basis van vermoeibaarheid en mechanische kritikaliteit, wat een praktisch alternatief biedt voor de computationeel dure Bayesiaanse optimalisatie terwijl het vergelijkbare prestaties levert.
• Identificatie van Faalmechanismen: De studie identificeert dat taalfalen bij handfietsen niet louter toe te schrijven is aan globale vermoeidheid, maar wordt gedreven door specifieke spieren (met name de anterieure deltoïdeus en de biceps) die verzadiging bereiken, wat de noodzaak benadrukt voor een adaptieve herverdeling van inspanning.

Resultaten

• Prestaties van de Kostenfunctie: De kostenfunctie voor het minimaliseren van vermoeidheid ($\phi_{a'}^2$) leverde het beste resultaat op wat betreft duur, waarbij 74,7% meer cycli werden behaald dan de standaard pulswijdte-minimalisatie ($\phi_{pw}^1$) en 42,9% meer dan de minimalisatie van spierstress. De geoptimaliseerde variabele (vermoeidheid) verklaarde 49,8% van de variantie in duur, terwijl het normarchetype 20,4% verklaarde.
• Bayesiaanse Optimalisatie: Het optimaliseren van spierspecifieke gewichten via Bayesiaanse methoden verhoogde de duur met nog eens 40,6% (bereikend van 1.729 cycli) vergeleken met de ongegewogen vermoeidheidsminimalisatie. Deze strategie slaagde erin de inspanning te herverdelen naar de meer vermoeidheidsbestendige triceps en weg van de snel vermoeiende anterieure deltoïdeus. Deze benadering vereiste echter 135 uur aan berekeningen, wat de directe klinische toepasbaarheid beperkt.
• Fysiologische Weging: De voorgestelde heuristische methode, die zowel spiervermoeibaarheid als mechanische kritikaliteit incorporeerde, behaalde een toename van 37,6% in duur (1.692 cycli) vergeleken met de ongegewogen baseline. Dit resultaat was bijna identiek aan het Bayesiaans-geoptimaliseerde resultaat, maar werd in minder dan één seconde berekend.
• Stimulatiepatronen: Geoptimaliseerde strategieën verschoven van een uniforme verdeling naar adaptieve patronen. Zo bevorderde vermoeidheidsminimalisatie korte stimulatiebursts met grotere pulswijdtes, terwijl pulswijdte-minimalisatie de voorkeur gaf aan lage, verdeelde controle. De gewogen strategieën brachten de progressie van vermoeidheid in balans en vertraagden het punt waarop kritieke spieren (anterieure deltoïdeus, biceps) 50% van hun krachtcapaciteit bereikten.

Significantie
Het artikel betoogt dat het maximaliseren van de taakduur in FES-revalidatie vereist dat men verder kijkt dan eenvoudige lading-minimalisatie. De primaire significantie ligt in het aantonen dat het rekening houden met de veelzijdige aard van uithoudingsvermogen — specifiek de interactie tussen spierspecifieke vermoeibaarheid en mechanische kritikaliteit — essentieel is voor het verlengen van FES-gestuurde activiteiten.

De studie concludeert dat hoewel het minimaliseren van de vermoeidheidsvariabele het meest effectieve kostenfunctie-archetype is, de prestaties ervan sterk afhankelijk zijn van hoe de controller de inspanning verdeelt over redundante spieren. De voorgestelde weegmethode biedt een praktische, interpreteerbare en computationeel efficiënte oplossing om stimulatiestrategieën te personaliseren zonder de prohibitieve kosten van iteratieve optimalisatie. Door het voortijdige uitputting van mechanisch kritieke of zeer vermoeibare spieren te voorkomen, ondersteunen deze strategieën het potentieel voor langere, effectievere revalidatiesessies, waarmee een belangrijke barrière voor de klinische adoptie van FES wordt aangepakt.