1 Uvod Pri poškodbi osrednjega živčnega sistema lahko nastopi popolna ali delna paraliza gornjih ekstremitet, zaradi česar ima bolnik omejene zmožnosti opravljanja vsakdanjih aktivnosti. Obstoječa rehabilitacija ohromljenih gornjih okončin temelji večinoma na fizioterapevtski vadbi, ki z dolgotrajnim ponavljanjem gibov pacientu do neke mere omogoči povrnitev hotene kontrole gibanja. Kljub različnim metodam za oceno funkcionalnosti roke je končna ocena povrnitve funkcionalnosti odvisna od fizioterapevta in njegove izkušenosti, torej je subjektivna. Za oceno senzorimotorične kontrole med prijemom, kjer nastopajo sile manjše od maksimalne sile prijema, so se za primerne izkazale naloge sledenja sile [1, 2]. Pri izvajanju naloge sledenja sile oseba s silo prijema sledi želeni vrednosti referenčnega signala, ki je prikazan na zaslonu. Poleg ocene prijema lahko naloge sledenja uporabimo tudi za vadbo. V študijah je bilo pokazano, da uporaba sistema za sledenje sile izboljša točnost kontrole sile prijema in poveča maksimalno silo prijema [2, 3]. Kot dodatek k navadni terapiji je mogoče koristno uporabiti tudi funkcionalno električno stimulacijo (FES) [4-6]. FES je metoda, pri kateri z električnimi impulzi vzdražimo periferni živec, ki oživčuje določeno mišično skupino, s čimer umetno vzbudimo mišično kontrakcijo in dosežemo funkcionalen gib. Večina sistemov FES, ki so danes uporabljeni v klinični praksi, je vodenih odprtozančno, kar pomeni, da so vzorci stimulacije prednastavljeni oz. jih spreminja terapevt ali uporabnik. Če želimo imeti avtomatsko vodenje FES, je treba uporabiti zaprtozančno voden sistem FES. Pričujoča študija predstavlja sistem za sledenje sile prijema v kombinaciji z zaprtozančno vodeno FES. Sistem je bil razvit za vadbo prijema bolnikov s poškodbami osrednjega živčnega sistema, ki imajo delno ohranjeno hoteno kontrolo roke. Zaprtozančni sistem FES je pri tem uporabljen kot pomoč bolniku, da lažje in bolje sledi referenčni vrednosti sile prijema. V delu so v drugem poglavju opisane glavne komponente sistema. Tretje poglavje podaja korake s katerimi prilagodimo sistem uporabniku, v četrtem pa je podana metodologija evalvacije sistema v praksi. V petem poglavju so predstavljeni rezultati, v zadnjem poglavju pa podani sklepi. 2 Opis sistema Na sliki 1 je prikazana blokovna shema funkcionalnih sklopov, ki sestavljajo sistem za krepitev senzomotoričnih sposobnosti gornjih okončin. Ti so naprava za merjenje sile prijema, električni stimulator ter osebni računalnik s prikazovalnikom. - stimulator naprava za merjenje sile prijema elektrode povratna vizualna informacija regulator referenca sile PC PI Slika 1: Blokovna shema sistema Figure 1. Conceptual scheme of the system Za zajemanje sil prijema je bila razvita naprava, ki je prikazana na sliki 2. Njeno ogrodje je sestavljeno iz aluminijastih profilov, na katere sta pritrjena dva senzorja za merjenje sil in navorov (JR3, 50M31A-I25, Inc., Woodland, CA, USA) ter podpora za podlaket. Kot je razvidno s slike 2, je levi senzor namenjen merjenju sil palca. Palec je pritrjen na senzor prek podpore in Velcro® traku. Desni senzor je namenjen merjenju sile preostalih štirih prstov roke. Pritrditev za prste je izvedena z dvema vzporednima aluminijastima profiloma, s katerima popolnoma omejimo gibanje prstov v smeri fleksije in ekstenzije. Vse podpore za prste so obložene z neoprensko gumo. Za preprečitev premikanja roke med vadbo podlaket s trakovi pritrdimo v oporo za podlaket. Pri napravi za merjenje sile prijema je mogoče nastaviti pozicijo obeh senzorjev in opore za podlaket, kar nam omogoča, da napravo priredimo za merjenje desne ali leve roke. Podatki z obeh senzorjev se zajemajo s vmesniško kartico PCI s frekvenco 100 Hz ter so filtrirani v realnem času z na kartici intgriranim filtrom z mejno frekvenco 31,25 Hz. Za stimulacijo je uporabljen štirikanalni električni stimulator s površinskimi elektrodami, ki omogoča računalniško vodeno spreminjanje jakosti stimulacije. Izhodi stimulatorja so monofazni tokovni impulzi z maksimalnim trajanjem 1,29 ms in maksimalno amplitudo 100 mA. Trajanje posameznega impulza lahko določimo na 10 µs, amplitudo pa na 0,4 mA natančno. Jakost stimulacije je vodena s proporcionalno-integralnim (PI) regulatorjem, ki spreminja pulzno širino stimulacijskih impulzov na podlagi razlike med želeno in dejansko vrednostjo sile prijema, pri čemer je amplituda stimulacijskih impulzov konstantna. Vrednosti o amplitudi in trajanju vsakega stimulacijskega impulza so stimulatorju posredovane s serijsko povezavo RS232 s frekvenco 33 Hz, ki je hkrati tudi frekvenca stimulacijskih impulzov. Vodenje sistema FES za krepitev senzomotoričnih sposobnosti gornjih ekstremitet 189 Slika 2: Naprava za merjenje sile prijema a) JR3 senzorja sil in navorov, b) opora za podlakt, c) in d) pritrditvi za prste Figure 2. Hand-force measuring device a) two JR3 force/torque sensors, b) forearm support c) and d) finger supports for fixation of fingers Razvita programska oprema teče na osebnem računalniku pod operacijskim sistemom Windows XP. S programsko opremo zajemamo dejansko silo prijema, generiramo referenčno silo, prikazujemo dejansko in referenčno silo na prikazovalniku ter vodimo stimulator. Poleg izvajanja naloge sledenja sile nam program omogoča tudi merjenje maksimalne sile prijema. Pri sili prijema upoštevamo le z komponento vektorja sile, ki deluje pravokotno na senzor sile v smeri fleksije/ekstenzije prstov. 3 Prilagoditev sistema uporabniku Ker je odziv mišice na električno stimulacijo odvisen od postavitve elektrod, je treba pri vsaki uporabi sistem najprej prilagoditi uporabniku. Prilagoditev sistema je izvedena v več korakih. 3.1 Namestitev elektrod Prvi korak je namestitev površinskih elektrod na podlaket. Eno od elektrod namestimo v bližini zapestja, drugo pa više na podlakti, in sicer tako, da stimulacija povzroči fleksijo/ekstenzijo prstov brez hkratne fleksije/ekstenzije zapestja. Med nameščanjem elektrod je pomembno, da oseba drži roko v enakem položaju, kot se bo le ta nahajala med vadbo. S spremembo položaja roke se namreč spremeni tudi položaj elektrod na koži glede na mišice s tem pa tudi odziv na električno stimulacijo. Po namestitvi elektrod prste in podlaket namestimo in pritrdimo v napravo. 3.2 Identifikacija modela mišice V naslednjem koraku sledi identifikacija modela električno stimulirane mišice v izometričnih razmerah. Za model mišice smo uporabili Hammersteinov model mišice [7], ki je prikazan na sliki 3. Model tvorita statična nelinearnost mišice fr (angl. recruitment curve), ki predstavlja nelinearen odziv mišice na električno stimulacijo, in linearna diskretna prenosna funkcija, ki opisuje dinamični odziv mišice. Signali v predstavljenem modelu so: pw(t) - pulzna širina stimulacijskih pulzov ã(t) - aktivacijski nivo mišice f(t) - sila, ki je posledica kontrakcije mišice q B (q ) - -1k f A (q )f -1 fr pw(t) ã(t) f(t) mišica Slika 3. Hammersteinov model mišice Figure 3. Hammerstein model of the muscle V linerani diskretni prenosni funkciji s parametrom k ( 1≥k ) določimo diskretno zakasnitev, Af in Bf pa sta polinoma operatorja zakasnitve q-1 in imata naslednjo obliko: na nf qaqaqA −−− +⋅⋅⋅++= 11 1 1)( (1) nb nf qbqbbqB −−− +⋅⋅⋅++= 110 1)( (2) pri čemer je stopnja polinoma Af ponavadi za ena manjša kot stopnja Bf ( 1−= nanb ). Iz literature [7] sledi, da dinamiko mišice najbolje opisujeta strukturi prenosne funkcije (na, nb, k) = (3, 2, 1) in (2, 1, 1). Pri modeliranju smo, po preizkusu obeh modelov, izbrali prvo možnost: 3 3 2 2 1 1 2 2 1 10 1 1 11 1 )( )( )( −−− −−− − −− +++ ++ = qaqaqa qbqbbq qA qBq f f (3) Za določitev statične nelinearnosti mišice smo uporabili stimulacijski signal pri kateremu smo povečevali trajanje stimulacijskih impulzov od 0 do 500 µs, po korakih 10 µs. Frekvenca stimulacijskih impulzov je bila 33 Hz, amplituda pa konstantna. Meritev izvedemo petkrat z vmesnimi premori, ki trajajo 15 s. Približek statične nelinearnosti mišice izračunamo s povprečenjem rezultatov vseh petih meritev. Po določitvi statične nelinearnosti mišice je treba določiti parametre prenosne funkcije (3). Te določimo z identifikacijskim postopkom odziva na PRBS (angl. pseudo-random binary sequence) signal [7]. Ustrezno 190 Perdan, Kamnik, Obreza, Bajd, Munih skaliran signal PRBS določa v identifikacijskem postopku širino stimulacijskih impulzov. Vhodni signal PRBS vodimo v stimulator prek nasprotne statične nelinearne karakteristike mišice, s čimer izničimo vpliv nelinearnega dela modela mišice. Odziv mišice izmerimo z napravo za merjenje sile prijema. Vhodni signal PRBS in izmerjeni odziv nato uporabimo kot vhodna parametra funkcije arx programskega paketa Matlab. Rešitev funkcije arx so koeficienti ai in bi polinomov Af in Bf . 3.3 Nastavitev ojačenj KI in KP PI regulatorja Za nastavitev ojačenja KP proporcionalnega dela in ojačenja KI integralnega dela PI regulatorja uporabimo orodje za optimizacijo v programskem okolju Matlab- Simulink. Optimizacijsko orodje v simulacijskih tekih na podlagi razlik med referenčnim signalom in odzivom sistema spreminja ojačenja regulatorja za zagotovitev najmanjšega odstopanja. Shematično je postopek optimizacije regulatorja prikazan na sliki 4. Za metodo iskanja je bila izbrana metoda Latin hypercube, začetni vrednosti ojačenj pa sta bili KI = 0,1 in KP = 1. - Model šicemi Referenca sile PI regulator Signal Constraint blok - K , KP I Optimizacija sledenja Slika 4: Shematičen prikaz nastavitve regulatorja Figure 4. Conceptual scheme of controller tuning 4 Metodologija evalvacije sistema Pri evalvaciji sistema je sodelovalo pet zdravih (spol moški, starost 30,6 ± 4,9 let) in dve tetraparetični osebi (AD: spol moški, starost 28 let, poškodba C5/C6; AS: spol moški, starost 15 let, poškodba C3/C4), ki so sistem uporabljali štiri dni zaporedoma. Vsaka oseba je vsak dan opravila tri vaje sledenja za fleksorje in tri za ekstenzorje prstov desne roke. Prvi trije dnevi so bili namenjeni privajanju, sistem vodenja FES pa je bil preizkušen četrti dan. Preizkus je bil izveden v dveh načinih delovanja. Pri zdravih osebah je bilo preizkušeno avtomatsko vodenje sile prijema s FES. Pri tem zdrave osebe niso smele s hoteno kontrakcijo mišic vplivati na stimulacijo FES, kar je bilo doseženo tako, da jim ni bila posredovana vizualna povratna informacija o referenčni in dejanski sili. Sila prijema je bila dosežena izključno s pomočjo delovanja zaprtozančnega sistema FES. Pri obeh tetraparetičnih osebah pa je bilo preizkušeno delovanje zaprtozančnega sistema FES v interakciji s hoteno aktivnostjo osebe. V tem načinu sta k sili prijema pomagala delna hotena aktivnost osebe kot tudi umetno vzbujena kontrakcija mišic. Referenčni signal želene sile je bil sestavljen iz štirih posameznih period sinusnega signala frekvenc 0,025, 0,033, 0,05 in 0,1 Hz, ki so si sledile ena za drugo z vmesnimi premori po 15 s. Amplituda signala je bila nastavljena na 50 odstotkov maksimalne sile prijema, ki jo je oseba sposobna doseči s hoteno aktivnostjo. Za evalvacijo sledenja je bila izračunana relativna napaka sledenja rrmse (angl. relative root mean square error), ki podaja relativno napako med referenčno Fref in izmerjeno Fizm silo: ∑ = = − = Tt t ref izmref F tFtF T rrmse 0 2 2 )max( ))()((1 (4) 5 Rezultati Na sliki 5 je prikazan primer izmerjenih statičnih nelinearnosti mišic fleksorjev in ekstenzorjev prstov za zdravo osebo ter obe tetraparetični osebi. Iz slike je razvidno, da se karakteristike med seboj razlikujejo tako po doseženi sili v odvisnosti od impulzne širine (pw), kot tudi po širini mrtvega območja, območja naraščanja in območja nasičenja. Tetraparetična oseba AD je imela zaradi spastičnih fleksorjev vidno slabši odziv ekstenzorjev prstov (slika 5, spodnji graf). 50 100 150 200 250 300 350 400 450 500 0 20 40 60 80 Stati na nelinearnost mišice: fleksorjič pw [µs] F [ N ] tetrapareti na oseba ADč tetrapareti na oseba ASč zdrava oseba 0 100 200 300 400 500 0 20 40 60 Stati na nelinearnost mišice: ekstenzorjič pw [µs] F [ N ] tetrapareti na oseba ADč tetrapareti na oseba ADč zdrava oseba Slika 5: Primer izmerjenih statičnih nelinearnosti mišic fleksorjev (zgoraj) in ekstenzorjev prstov (spodaj) Figure 5. Example of static recruitment nonlinearity for finger flexors (top graph) and extensors (bottom graph) Na sliki 6 je prikazan primer izmerjenega odziva na signal PRBS in izračunan odziv modela na isti signal. Slika 7 prikazuje primer rezultatov naloge sledenja sile pri zdravi osebi. Zgornji graf prikazuje rezultat sledenja v primeru, ko je zdrava oseba sledila želeni sili le s hoteno kontrolo mišic in pri tem dosegla napako sledenja rrmse = 0,19. Na spodnjem grafu slike 7 je prikazano sledenje referenčni sili v primeru, ko je ta Vodenje sistema FES za krepitev senzomotoričnih sposobnosti gornjih ekstremitet 191 dosežena samo z avtomatsko vodeno FES. Pri tem je bila dosežena napaka sledenja rrmse = 0,29. Amplituda referenčnega signala je bila pri avtomatsko vodeni FES znižana na 20 % maksimalne sile prijema zaradi lažjega prenašanja bolečine, ki lahko nastopi pri električni stimulaciji. 0 2 4 6 8 10 12 0 10 20 30 40 50 60 t [s] F [ N ] model mišica Slika 6: Izmerjeni odziv mišic fleksorjev prstov in odziv modela mišice na signal PRBS Figure 6. Measured response of the muscle to PRBS and model response to PRBS Slika 7: Primer sledenja referenčne sile zdrave osebe za fleksorje prstov (a) hotena kontrola brez pomoči FES, (b) samo FES Figure 7. Example of force tracking in a healthy subject for finger flexors (a) only voluntary activity (b) only FES Primer sledenja tetraparetične osebe AD je prikazan na sliki 8. Zgornji graf prikazuje primer, ko oseba ni mogla dobro slediti referenčnemu signalu. Večja odstopanja so opazna predvsem kadar je oseba morala dalj časa vzdrževati večjo silo prijema. Manjša odstopanja so vidna tudi pri kratkotrajnem vzdrževanju sile prijema. Sredinski graf na sliki 8 prikazuje sledenje referenčni sili s pomočjo FES in kot je razvidno, je sledenje znatno boljše. V spodnjem grafu slike 8 je prikazano spreminjanje impulzne širine stimulacijskih impulzov pri sledenju referenčni sili s pomočjo FES. Slika 8: Primer sledenja referenčne sile pri tetraparetični osebi Figure 8. Force tracking task of an incomplete tetraplegic person V tabeli 1 so zbrani podatki povprečnih napak sledenja rrmse za četrti dan merjenj zdravih oseb. Iz rezultatov je razvidno, da regulator pri nalogi sledenja za fleksorje slabše sledi referenci kot zdrave osebe, pri ekstenzorjih prstov pa so rezultati primerljivi. rrmse hotena kontrola FES fleksorji 0,22 ± 0,06 0,47 ± 0,15 ekstenzorji 0,41 ± 0,17 0,46 ± 0,33 Tabela 1: Povprečna vrednost rrmse za zdrave osebe pri hoteni in FES vodeni nalogi sledenja Tabel 1. Average rrmse value for healthy subjects: voluntary control and a FES-driven tracking task rrmse hotena kontrola FES + hotena fleksorji 0,64 ± 0,30 0,56 ± 0,35 ekstenzorji 0,97 ± 1,03 0,90 ± 0,47 Tabela 2: Povprečna vrednost rrmse za tetraparetični osebi pri hoteno ter FES + hoteno vodeni nalogi sledenja Tabel 2. Average rrmse value for incomplete tetraplegic subjects: voluntary and FES + voluntary driven tracking task V tabeli 2 so zbrani podatki povprečnih napak sledenja rrmse za četrti dan merjenj obeh tetraparetičnih oseb. Iz rezultatov je razvidno, da pri tetraparetičnih osebah FES rahlo pripomore k boljšemu rezultatu sile prijema. 192 Perdan, Kamnik, Obreza, Bajd, Munih 6 Sklep V pričujoči študiji je predstavljen sistem za krepitev senzomotoričnih sposobnosti gornjih okončin oseb z delno kontrolo prijema. Sistem je namenjen vadbi mišic fleksorjev in ekstenzorjev prstov roke v izometričnih razmerah. Vadba temelji na izvajanju naloge sledenja sile s pomočjo zaprtozančno vodene FES. Članek predstavlja sistem, postopek prilagoditve uporabniku in rezultate evalvacije zaprtozančnega sistema vodenja FES. V postopku prilagoditve sistema uporabniku je zaprtozančni sistem vodenja FES po postavitvi elektrod uglašen na odziv mišic prek identifikacije modela mišice in optimizacije parametrov regulatorja glede na kakovost sledenja referenčnega signala. V evalvaciji zaprtozančnega sistema vodenja FES je sodelovalo pet zdravih in dve tetraparetični osebi. Testne osebe so po treh dneh privajalne vadbe izvedle testno sledenje sile s hoteno aktivnostjo in s pomočjo FES. Rezultati pri zdravih osebah, ki so nalogo sledenja izvajale samo s hoteno kontrolo prijema ali samo s pomočjo FES, so pokazali, da so napake sledenja sile pri prijemu, vodenim s FES, pri ekstenzorjih primerljive s tistimi, ki jih dosegajo zdrave osebe s hoteno kontrolo prijema. Nekoliko slabši rezultati sledenja s pomočjo FES so bili doseženi pri fleksiji prstov. Rezultati testiranja pri tetraparetičnih osebah, kjer je bil zaprtozančni sistem vodenja FES preizkušen v kombinaciji s hoteno aktivacijo mišic, kažejo, da je sistem FES tudi v tem primeru kljub slabšim rezultatom v primerjavi z zdravimi osebami rahlo izboljšal kakovost sledenja sile. Predstavljeni rezultati kažejo, da je postopek uglasitve zaprtozančnega sistema FES ustrezen in da zaprtozančni sistem deluje robustno tudi pri motnjah, ki nastanejo pri dodatni hoteni aktivnosti trenirane osebe. Ob delni hoteni kontroli zaprtozančni FES sistem pripomore k boljši kakovosti sledenja referenčne sile, kar poleg večje mišične aktivnosti dodatno motivira trenirano osebo. Na tej podlagi je mogoče povzeti, da je sistem za krepitev senzomotoričnih sposobnosti gornjih okončin primeren za uporabo v klinični praksi. 7 Literatura [1] G. Kurillo, A. Zupan, T. Bajd, Force tracking system for the assessment of grip force control in patients with neuromuscular diseases. Clinical Biomechanics, 19, 1014–1021, 2004. [2] G. Kriz, J. Hermsdörfer, C. Marquardt, N. Mai, Feedback-based training of grip force control in patients with brain damage, Archives of Physical Medicine and Rehabilitation, 76, 653–659, 1995. [3] G. Kurillo, M. Gregorič, N. Goljar, T. Bajd, Grip force tracking system for assessment and rehabilitation of hand function. Technology and Health Care, 13, 137–149, 2005. [4] B. R. Bowman, L. L. Baker, R. L. Waters, Positional feedback and electrical stimulation: an automated treatment for the hemiplegic wrist, Archives of Physical Medicine and Rehabilitation, 60, 497–501, 1979. [5] M. R. Popović, A. Curt, T. Keller, V. Dietz, Functional electrical stimulation for grasping and walking: indications and limitations, Spinal Cord, 39, 403–412, 2001. [6] K. L. Kilgore, P. H. Peckham. G. B. Thorpe, M. W. Keith, K. Gallaher-Stone, Synthesis of hand grasp using functional neuromuscular stimulation, IEEE Transactions on Biomedical Engineering, 36, 7, 761– 770, 1989. [7] K. J. Hunt, M. Munih, N. N. Donaldson, F. M. D. Barr, Investigation of the Hammerstein hypothesis in the modeling of electrically stimulated muscle, IEEE Transactions on Biomedical Engineering, 45, 8, 998– 1008, 1998. Jernej Perdan je na Fakulteti za elektrotehniko Univerze v Ljubljani mladi raziskovalec podjetja Iskra Medical. Njegovo raziskovalno področje so rehabilitacijski sistemi gornjih okončin. Roman Kamnik je docent na Fakulteti za elektrotehniko Univerze v Ljubljani. Njegovo raziskovalno področje so biomehanika, rehabilitacijska in mobilna robotika. Pavla Obreza je fizioterapevtka na Inštitutu Republike Slovenije za rehabilitacijo. Tadej Bajd je redni profesor na Fakulteti za elektrotehniko Univerze v Ljubljani in izredni član SAZU. Njegovo raziskovalno področje dela je analiza in sinteza gibanja pri človeku in stroju. Marko Munih je redni profesor na Fakulteti za elektrotehniko Univerze v Ljubljani, kjer je predstojnik Laboratorija za robotiko in biomedicinsko tehniko. Njegova področja dela so rehabilitacijska in industrijska robotika ter mikroračunalniški in biomedicinski tehnični sistemi.